Remap LonLat Grid to Lambert Conformal Projection with constant 5.5km spacing, but retain LatLon Information?
Added by Simon Hentschel 10 months ago
Hello,
I'm currently trying to remap ERA5 data to a Lambert Conformal Projection with 5.5km spacing, similar to that of CERRA, but moved to cover North America instead.
For this I'm using the "cdo remapbic" command with the grid-description from below.
While the reprojection itself works fine, for further processing I need a mapping from each grid cell to the coordinates it represents.
The converted file itself sadly only contains the x/y variables with the position encoded in 'meters' from the first grid point, instead of the lat/lon values.
Is there any way to retain the lat,lon coordinate map (similar to the original CERRA files)?
Thank you very much for your help!
Simon
#- gridID similar to CERRA (moved to cover North America)
#
gridtype = projection
gridsize = 1142761
xsize = 1069
ysize = 1069
xname = x
xunits = "m"
yname = y
yunits = "m"
xfirst = 0
xinc = 5500
yfirst = 0
yinc = 5500
grid_mapping = Lambert_Conformal
grid_mapping_name = lambert_conformal_conic
standard_parallel = 50.
earth_radius = 6371229.
longitude_of_central_meridian = -94.514057
latitude_of_projection_origin = 50.
longitudeOfFirstGridPointInDegrees = -120
latitudeOfFirstGridPointInDegrees = 20.292281
longitudeOfSouthernPoleInDegrees = 0.
latitudeOfSouthernPoleInDegrees = -90.
Replies (1)
RE: Remap LonLat Grid to Lambert Conformal Projection with constant 5.5km spacing, but retain LatLon Information? - Added by Uwe Schulzweida 10 months ago
Hello Simon,
You get lon/lat values if you add setgridtype,curvilinear after the interpolation:
cdo -setgridtype,curvilinear -remapbil,<targetgrid> infile outfileCheers,
Uwe