main@rasterio.groups.io | Output to grib2?

Home Messages Hashtags Subgroups

Date Date 1 - 7 of 7

Output to grib2?

Shane Mill - NOAA Affiliate #154 Is there anyway to output to grib2 after performing map algebra? For example, If you are performing a calculation between two grib2 bands, which rasterio does quite well, are you able to transform the numpy array into a grib2 file? I have been able to do this with geotiff but not grib2.
More All Messages By This Member
Sean Gillies #155 Hi Shane, Rasterio's rio-convert command can convert from GeoTIFF to GRIB, so it should be possible $ rio convert -f GRIB ~/code/rasterio/tests/data/RGB.byte.tif /tmp/rgb.grb2 $ rio info /tmp/rgb.grb2 {"bounds": [101985.0, 2611486.28, 339316.64, 2826915.0], "colorinterp": ["undefined", "undefined", "undefined"], "count": 3, "crs": "PROJCS[\"unnamed\",GEOGCS[\"Coordinate System imported from GRIB file\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",0],PARAMETER[\"central_meridian\",-2072.483648],PARAMETER[\"scale_factor\",0.9996],PARAMETER[\"false_easting\",500000],PARAMETER[\"false_northing\",0],UNIT[\"Metre\",1]]", "descriptions": ["0[-] RESERVED(0) (Reserved)", "0[-] RESERVED(0) (Reserved)", "0[-] RESERVED(0) (Reserved)"], "driver": "GRIB", "dtype": "float64", "height": 718, "indexes": [1, 2, 3], "lnglat": [84.75845370919801, 24.561589206961006], "mask_flags": [["nodata"], ["nodata"], ["nodata"]], "nodata": 0.0, "res": [300.04, 300.04], "shape": [718, 791], "tiled": false, "transform": [300.04, 0.0, 101985.0, 0.0, -300.04, 2826915.0, 0.0, 0.0, 1.0], "units": [null, null, null], "width": 791} I'm not a regular GRIB user, but it looks to me that GDAL produces a GRIB2 file when you write using the GRIB format driver, and that the format driver reads either version without a need to specify. GRIB2 writing requires GDAL 2.3.0 or newer. The latest rasterio wheels on PyPI include GDAL 2.4.0, and should be sufficient. On Fri, Mar 15, 2019 at 11:53 AM Shane Mill - NOAA Affiliate via Groups.Io <shane.mill=noaa.gov@groups.io> wrote: toggle quoted messageShow quoted text Is there anyway to output to grib2 after performing map algebra? For example, If you are performing a calculation between two grib2 bands, which rasterio does quite well, are you able to transform the numpy array into a grib2 file? I have been able to do this with geotiff but not grib2. -- Sean Gillies
More All Messages By This Member
Shane Mill - NOAA Affiliate #156 On Fri, Mar 15, 2019 at 11:02 AM, Sean Gillies wrote: rio convert -f GRIB ~/code/rasterio/tests/data/RGB.byte.tif /tmp/rgb.grb2 Sean, thanks for the quick response! You are correct, that does work. I think my issue may have been that I was using gdalinfo to check the resulting grib2 file, and my gdal version was older. When I check with rio info it seems to be correct. Thanks again! -Shane
More All Messages By This Member
Shane Mill - NOAA Affiliate #159 Hey Sean, I know that you were saying that you don't use grib2 often, but I wanted to follow up with the previous question. I am able to output the numpy array to GRIB. Say that you have one band, you can read the tags of that one band with src.tags(1). These tags will be the defaults that come from: https://gdal.org/frmt_grib.html (if you scroll down to the section "GRIB2 Write Support"). The problem with this is that for the grib2 to be readable, you need to be able to update these tags. I'm wondering how to change the tags either when originally creating the GRIB or after the fact when using "update_tags" explained here: https://github.com/mapbox/rasterio/blob/master/docs/topics/tags.rst It appears that update_tags works fine for GTiff but not for GRIB. Let me know if you have any thoughts on this. Thanks as always! -Shane On Fri, Mar 15, 2019 at 2:16 PM Shane Mill - NOAA Affiliate via Groups.Io <shane.mill=noaa.gov@groups.io> wrote: toggle quoted messageShow quoted text On Fri, Mar 15, 2019 at 11:02 AM, Sean Gillies wrote: rio convert -f GRIB ~/code/rasterio/tests/data/RGB.byte.tif /tmp/rgb.grb2 Sean, thanks for the quick response! You are correct, that does work. I think my issue may have been that I was using gdalinfo to check the resulting grib2 file, and my gdal version was older. When I check with rio info it seems to be correct. Thanks again! -Shane -- Shane Mill Meteorological Application Developer, AceInfo Solutions Meteorological Development Laboratory (NWS) Phone: 301-427-9452
More All Messages By This Member
Shane Mill - NOAA Affiliate #162 Okay so I think that I can update the GRIB specific options using rio warp ie: rio warp 1000_500_thick.grb 1000_500_thick_new2.grb --format GRIB --profile BAND_1_PDS_PDTN= with BAND_1_PDS_PDTN being a GRIB specific option. I'm going to explore this more today. Thanks, Shane On Tue, Mar 19, 2019 at 3:50 PM Shane Mill - NOAA Affiliate via Groups.Io <shane.mill=noaa.gov@groups.io> wrote: toggle quoted messageShow quoted text Hey Sean, I know that you were saying that you don't use grib2 often, but I wanted to follow up with the previous question. I am able to output the numpy array to GRIB. Say that you have one band, you can read the tags of that one band with src.tags(1). These tags will be the defaults that come from: https://gdal.org/frmt_grib.html (if you scroll down to the section "GRIB2 Write Support"). The problem with this is that for the grib2 to be readable, you need to be able to update these tags. I'm wondering how to change the tags either when originally creating the GRIB or after the fact when using "update_tags" explained here: https://github.com/mapbox/rasterio/blob/master/docs/topics/tags.rst It appears that update_tags works fine for GTiff but not for GRIB. Let me know if you have any thoughts on this. Thanks as always! -Shane On Fri, Mar 15, 2019 at 2:16 PM Shane Mill - NOAA Affiliate via Groups.Io <shane.mill=noaa.gov@groups.io> wrote: On Fri, Mar 15, 2019 at 11:02 AM, Sean Gillies wrote: rio convert -f GRIB ~/code/rasterio/tests/data/RGB.byte.tif /tmp/rgb.grb2 Sean, thanks for the quick response! You are correct, that does work. I think my issue may have been that I was using gdalinfo to check the resulting grib2 file, and my gdal version was older. When I check with rio info it seems to be correct. Thanks again! -Shane -- Shane Mill Meteorological Application Developer, AceInfo Solutions Meteorological Development Laboratory (NWS) Phone: 301-427-9452 -- Shane Mill Meteorological Application Developer, AceInfo Solutions Meteorological Development Laboratory (NWS) Phone: 301-427-9452
More All Messages By This Member
Sean Gillies #164 Hi Shane, I think I understand now. BAND_1_PDS_PDTN isn't a tag, it's a "creation option". You can apply these at the time you open a dataset for writing like with rasterio.open("file.grb2", "w", driver="GRIB", ...., BAND_1_PDS_PDTN=something) as dataset: dataset.write(data) On Fri, Mar 22, 2019 at 9:06 AM Shane Mill - NOAA Affiliate via Groups.Io <shane.mill=noaa.gov@groups.io> wrote: toggle quoted messageShow quoted text Okay so I think that I can update the GRIB specific options using rio warp ie: rio warp 1000_500_thick.grb 1000_500_thick_new2.grb --format GRIB --profile BAND_1_PDS_PDTN= with BAND_1_PDS_PDTN being a GRIB specific option. I'm going to explore this more today. Thanks, Shane On Tue, Mar 19, 2019 at 3:50 PM Shane Mill - NOAA Affiliate via Groups.Io <shane.mill=noaa.gov@groups.io> wrote: Hey Sean, I know that you were saying that you don't use grib2 often, but I wanted to follow up with the previous question. I am able to output the numpy array to GRIB. Say that you have one band, you can read the tags of that one band with src.tags(1). These tags will be the defaults that come from: https://gdal.org/frmt_grib.html (if you scroll down to the section "GRIB2 Write Support"). The problem with this is that for the grib2 to be readable, you need to be able to update these tags. I'm wondering how to change the tags either when originally creating the GRIB or after the fact when using "update_tags" explained here: https://github.com/mapbox/rasterio/blob/master/docs/topics/tags.rst It appears that update_tags works fine for GTiff but not for GRIB. Let me know if you have any thoughts on this. Thanks as always! -Shane On Fri, Mar 15, 2019 at 2:16 PM Shane Mill - NOAA Affiliate via Groups.Io <shane.mill=noaa.gov@groups.io> wrote: On Fri, Mar 15, 2019 at 11:02 AM, Sean Gillies wrote: rio convert -f GRIB ~/code/rasterio/tests/data/RGB.byte.tif /tmp/rgb.grb2 Sean, thanks for the quick response! You are correct, that does work. I think my issue may have been that I was using gdalinfo to check the resulting grib2 file, and my gdal version was older. When I check with rio info it seems to be correct. Thanks again! -Shane -- Shane Mill Meteorological Application Developer, AceInfo Solutions Meteorological Development Laboratory (NWS) Phone: 301-427-9452 -- Shane Mill Meteorological Application Developer, AceInfo Solutions Meteorological Development Laboratory (NWS) Phone: 301-427-9452 -- Sean Gillies
More All Messages By This Member
Shane Mill - NOAA Affiliate #165 Ah I see, I was confused then just because these come up when you "rio insp" and then "src.tags(1)", these creation options show up under the tags for band1. That's very helpful so thanks for the response! Shane toggle quoted messageShow quoted text On Fri, Mar 22, 2019 at 12:09 PM Sean Gillies <sean.gillies@...> wrote: Hi Shane, I think I understand now. BAND_1_PDS_PDTN isn't a tag, it's a "creation option". You can apply these at the time you open a dataset for writing like with rasterio.open("file.grb2", "w", driver="GRIB", ...., BAND_1_PDS_PDTN=something) as dataset: dataset.write(data) On Fri, Mar 22, 2019 at 9:06 AM Shane Mill - NOAA Affiliate via Groups.Io <shane.mill=noaa.gov@groups.io> wrote: Okay so I think that I can update the GRIB specific options using rio warp ie: rio warp 1000_500_thick.grb 1000_500_thick_new2.grb --format GRIB --profile BAND_1_PDS_PDTN= with BAND_1_PDS_PDTN being a GRIB specific option. I'm going to explore this more today. Thanks, Shane On Tue, Mar 19, 2019 at 3:50 PM Shane Mill - NOAA Affiliate via Groups.Io <shane.mill=noaa.gov@groups.io> wrote: Hey Sean, I know that you were saying that you don't use grib2 often, but I wanted to follow up with the previous question. I am able to output the numpy array to GRIB. Say that you have one band, you can read the tags of that one band with src.tags(1). These tags will be the defaults that come from: https://gdal.org/frmt_grib.html (if you scroll down to the section "GRIB2 Write Support"). The problem with this is that for the grib2 to be readable, you need to be able to update these tags. I'm wondering how to change the tags either when originally creating the GRIB or after the fact when using "update_tags" explained here: https://github.com/mapbox/rasterio/blob/master/docs/topics/tags.rst It appears that update_tags works fine for GTiff but not for GRIB. Let me know if you have any thoughts on this. Thanks as always! -Shane On Fri, Mar 15, 2019 at 2:16 PM Shane Mill - NOAA Affiliate via Groups.Io <shane.mill=noaa.gov@groups.io> wrote: On Fri, Mar 15, 2019 at 11:02 AM, Sean Gillies wrote: rio convert -f GRIB ~/code/rasterio/tests/data/RGB.byte.tif /tmp/rgb.grb2 Sean, thanks for the quick response! You are correct, that does work. I think my issue may have been that I was using gdalinfo to check the resulting grib2 file, and my gdal version was older. When I check with rio info it seems to be correct. Thanks again! -Shane -- Shane Mill Meteorological Application Developer, AceInfo Solutions Meteorological Development Laboratory (NWS) Phone: 301-427-9452 -- Shane Mill Meteorological Application Developer, AceInfo Solutions Meteorological Development Laboratory (NWS) Phone: 301-427-9452 -- Sean Gillies -- Shane Mill Meteorological Application Developer, AceInfo Solutions Meteorological Development Laboratory (NWS) Phone: 301-427-9452
More All Messages By This Member

1 - 7 of 7

1

Previous Topic Next Topic

Home Hashtags Subgroups About Features Pricing Updates Terms Help Twitter