How to create lats and lons from projected data (Lambert Conic Conformal)
Added by Ruth Mottram about 5 years ago
Hi all,
I have some data from a high resolution NWP model that is run on a Lambert Conic Conformal grid. I'd like ot convert these to lats and lons.
More specifically, I'd like to extract data from some specific points (I have geographic coordinates) where I have weather station observations.
Anyone have any ideas?
cdo griddes gives me this:
#- gridID 1
#
gridtype = projection
gridsize = 1356561
xsize = 1069
ysize = 1269
xname = x
xunits = "m"
yname = y
yunits = "m"
xfirst = 0
xinc = 2500
yfirst = 0
yinc = 2500
grid_mapping = Lambert_Conformal
grid_mapping_name = lambert_conformal_conic
standard_parallel = 72.
longitude_of_central_meridian = -36.
latitude_of_projection_origin = 72.
earth_radius = 6367470.
longitudeOfFirstGridPointInDegrees = -57.097
latitudeOfFirstGridPointInDegrees = 55.81
Best wishes
Ruth
Replies (2)
RE: How to create lats and lons from projected data (Lambert Conic Conformal) - Added by Uwe Schulzweida about 5 years ago
If your CDO installation was compiled with libproj (https://proj.org/) then you can use remapnn to select a grid point by lon/lat coordinates. Here is an example:
cdo remapnn,lon=10_lat=60 infile outfile
RE: How to create lats and lons from projected data (Lambert Conic Conformal) - Added by Uwe Schulzweida about 5 years ago
Use setgridtype,curvilinear to convert the Lambert coordinates to geographical lon/lat coordinates:
cdo setgridtype,curvilinear infile outfile