Extracting projection from grib files

Specific Domains Geo
3

Another way is with GMT.jl (that also uses GDAL for this)

using GMT

getproj("era5-2mt-2019-03-uk.grib")
"GEOGCS[\"Coordinate System imported from GRIB file\",DATUM[\"unnamed\",SPHEROID[\"Sphere\",6367470,0]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AXIS[\"Latitude\",NORTH],AXIS[\"Longitude\",EAST]]"

prj = wkt2proj(getproj("era5-2mt-2019-03-uk.grib"))
"+proj=longlat +R=6367470 +no_defs"
2 Likes