Saturday, 11 May 2019

gdal - How to process GRIB2 Data for Visualization?


My goal is to extract wind information from data provided to me in the GRIB2 format as it has been done in the attached image.


I was curious if anybody had any experience with this format. I'm not sure if I am going in the right direction or what format I should be exporting it to. My end goal is to display this information through openlayers. So far, I have used gdal_translate to extract specific information ex. Temperature from band 25 (see command below).


gdal_translate -of Gtiff -b 25 weather.grib2 band25.tiff

I also have access to FME Desktop professional edition.



Answer



I was able to figure out my issue



I used the WGRIB2 utility in order to convert my GRIB2 data to a CSV and subsequently into my POSTGIS database. I used the following command to take the VGRD, UGRD, and TMP at 100mb:


wgrib2 -match ":(VGRD|UGRD|TMP): 100mb" weather.grib2 -lon lon1 lat1 -csv csvfile

I then loaded my CSV to my POSTGIS database and overlayed it on top of my baselayer.


Here is my final result for my temperature layer at 100mb http://i.imgur.com/KEnz0.jpg


No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...