Hello, A follow up with some more information. According to Zachary Flamig, there are no issues in GDAL 2.2.0. I encountered these issues in 3.0.4. In 3.1.0, the warp works correctly, but gdal_translate still fails with an out of memory issue. It creates a TIF of 1.8gb before it errors out.
Thomas On Fri, Jun 5, 2020 at 11:42 AM Thomas Horner <[email protected]> wrote: > Hello list, > > The US National Weather Service recently released version 3.2 of their NBM > weather model blend. I have not been able to successfully to use gdalwarp > or gdal_translate on these GRIB2 files. For gdal_translate, I get an > "Error: Out of memory" error, even though it's only 80mb and I have 16gb of > RAM. > > This is what happens when using gdalwarp to warp to EPSG:4326: > > ------- > > thomas@eolus:~/model_processing/tempWxData$ gdalwarp > nbm_20200605_16Z_t001.grib2 nbm.tif -t_srs EPSG:4326 -multi > > Processing nbm_20200605_16Z_t001.grib2 [1/1] : 0Using internal nodata > values (e.g. 9999) for image nbm_20200605_16Z_t001.grib2. > > ERROR 1: Too many points (529 out of 529) failed to transform, unable to > compute output bounds. > Warning 1: Unable to compute source region for output window > 0,0,2345,1597, skipping. > -------- > > gdalinfo reveals the following information about the CRS: > > ------- > Driver: GRIB/GRIdded Binary (.grb, .grb2) > Files: nbm_20200605_16Z_t001.grib2 > Size is 2345, 1597 > Coordinate System is: > PROJCRS["unnamed", > BASEGEOGCRS["Coordinate System imported from GRIB file", > DATUM["unnamed", > ELLIPSOID["Sphere",6371200,0, > LENGTHUNIT["metre",1, > ID["EPSG",9001]]]], > PRIMEM["Greenwich",0, > ANGLEUNIT["degree",0.0174532925199433, > ID["EPSG",9122]]]], > CONVERSION["unnamed", > METHOD["Lambert Conic Conformal (2SP)", > ID["EPSG",9802]], > PARAMETER["Latitude of false origin",25, > ANGLEUNIT["degree",0.0174532925199433], > ID["EPSG",8821]], > PARAMETER["Longitude of false origin",265, > ANGLEUNIT["degree",0.0174532925199433], > ID["EPSG",8822]], > PARAMETER["Latitude of 1st standard parallel",25, > ANGLEUNIT["degree",0.0174532925199433], > ID["EPSG",8823]], > PARAMETER["Latitude of 2nd standard parallel",25, > ANGLEUNIT["degree",0.0174532925199433], > ID["EPSG",8824]], > PARAMETER["Easting at false origin",0, > LENGTHUNIT["Metre",1], > ID["EPSG",8826]], > PARAMETER["Northing at false origin",0, > LENGTHUNIT["Metre",1], > ID["EPSG",8827]]], > CS[Cartesian,2], > AXIS["easting",east, > ORDER[1], > LENGTHUNIT["Metre",1]], > AXIS["northing",north, > ORDER[2], > LENGTHUNIT["Metre",1]]] > Data axis to CRS axis mapping: 1,2 > Origin = (-3272421.457337170839310,3790842.106035435572267) > Pixel Size = (2539.702999999999975,-2539.702999999999975) > Corner Coordinates: > Upper Left (-3272421.457, 3790842.106) (138d23'39.07"W, 53d 3'24.76"N) > Lower Left (-3272421.457, -265063.585) (126d17' 8.03"W, 19d12'55.22"N) > Upper Right ( 2683182.078, 3790842.106) ( 59d 1'17.54"W, 54d22'46.84"N) > Lower Right ( 2683182.078, -265063.585) ( 69d11'54.98"W, 20d19' 6.30"N) > Center ( -294619.690, 1762889.261) ( 98d21'20.81"W, 40d37' 0.08"N) > ------- > > If you want a file to test with, download from here: > > https://nomads.ncep.noaa.gov/pub/data/nccf/com/blend/prod/blend.20200605/16/grib2/blend.t16z.master.f001.co.grib2 > > That file will be removed after a couple days though. > > I don't know the details of coordinate transformation, so I have no idea > how to solve this problem. Can anyone help? > > Thanks, > Thomas > > >
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
