Remap data from Polar Stereographic projection to Lat-Lon

Question

I have a NetCDF file with storm track density data over south polar stereographic projection. Both X and Y coordinates range from -1.924274 to 1.88504 (I don't know the units here). The link for the .nc file is here.

When I open the file with Xarray and plot the data, I have the following:

Is there any way to remap this to lat-lon grid in Python or CDO? Thanks!

You could try rioxarray: https://corteva.github.io/rioxarray/stable/examples/reproject.html — snowman2, Dec 24 '20 at 23:12

score 3 · Accepted Answer · answered Dec 26 '20 at 08:46

3

You can probably do this using my package nctoolkit (https://nctoolkit.readthedocs.io/en/latest/), which uses CDO as a back end.

import nctoolkit as nc
data = nc.open_data("infile.nc")
data.to_latlon(lon = [lon_min, lon_max], lat = [lat_min, lat_max])
data.plot()

You will need to define lon_min etc.

answered Dec 26 '20 at 08:46

Robert Wilson

3,192
11
19

1

I had to make some adjustments to my data, but it worked, thanks! – Maartaud Dec 30 '20 at 03:40

Remap data from Polar Stereographic projection to Lat-Lon

1 Answers1

Linked