Reproject BedMachine Antarctica to lat/lon

Use this Forum to find information on, or ask a question about, NASA Earth Science data.
Post Reply
EarthdataForumContributor
Posts: 284
Joined: Thu Jul 11, 2019 4:32 pm America/New_York
Answers: 3

Reproject BedMachine Antarctica to lat/lon

by EarthdataForumContributor » Thu Apr 28, 2022 4:27 pm America/New_York

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!

Tags:

NSIDC - danicalinda.cantarero
User Services
User Services
Posts: 4
Joined: Wed Jan 19, 2022 11:45 am America/New_York
Answers: 0

Re: Reproject BedMachine Antarctica to lat/lon

by NSIDC - danicalinda.cantarero » Thu Apr 28, 2022 4:31 pm America/New_York

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).

Post Reply