I’m trying to spatially subset a GFS 0.5 degree forecast read form a GRIB file downloaded from NOAA website (https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-004-0.5-degree/analysis/)
I can load the file using rasterio with cfgrib backend, but when it comes to clipping the dataset to a selected region, rasterio crashes with “No valid geometry objects found for rasterize” error message. Identical clipping operation works on MODIS files, so so I don’t expect problems here.