Page 1 of 1
Reproject BedMachine Antarctica to lat/lon
Posted: Thu Apr 28, 2022 4:27 pm America/New_York
by 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!
Re: Reproject BedMachine Antarctica to lat/lon
Posted: Thu Apr 28, 2022 4:31 pm America/New_York
by NSIDC - 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).