Reproject BedMachine Antarctica to lat/lon

[EarthdataForumContributor]

I have downloaded MEaSUREs BedMachine Antarctica data, and I would like to reproject it from Antarctic Polar Stereographic projection to geodetic latitude/longitude. I would appreciate any guidance on how to do this. Thanks!

[NSIDCx - danicalinda.cantarero]

You can take a look at the Python (xy2ll.py) and Matlab (xy2ll.m) scripts in this GitHub repository: https://github.com/nsidc/nsidc0756-scripts to convert between polar stereographic coordinates to geodetic lat/lon.

Here's an example Python script to use for the reprojection:

Code: Select all

from xy2ll import xy2ll
import numpy as np

xs = np.linspace(-3333000, -500, 6666)
ys = np.linspace(500, -3333000, 6666)

to_lat = []
to_lon = []

for x, y in zip(xs,ys):
	lat, lon = xy2ll(x, y, -1, 0, 71)
	to_lon.append(lon)
	to_lat.append(lat)

The GitHub repository mentioned above also includes scripts to convert from geodetic lat/lon to polar stereographic coordinates (ll2xy.m and ll2xy.py). Matlab and Python scripts for interpolating BedMachine Antarctica data onto user-defined coordinates are also available (interpBedmachineAntarctica.m and interp_bedmachine_antarctica.py).