[gdal-dev] NetCDF and EASE-2

Marsh, Chris chris.marsh at usask.ca
Fri Jan 31 20:07:16 PST 2020


Hi,
I'm trying to pull some data out of a netcdf and convert it to a geotiff. The data are the snowdepths from this paper

https://www.nature.com/articles/s41467-019-12566-y

with one of the data files here for testing:
https://drive.google.com/open?id=1qZBYJ4UEbYYHXpSV4gyFBkZYdYKtZGTR


The data ends up weirdly projected when using gdal_transform. I've tried the method outlined here
https://gis.stackexchange.com/questions/249841/parameters-of-projection-for-smap-data#comment392952_249841

but that ended up with an incorrect messy output. My attempt was
gdal_translate -a_ullr -17367530.45 7314540.11 17367530.45 -7314540.11 -a_srs "+proj=cea +lon_0=0 +lat_ts=30 +ellps=WGS84 +units=m" -a_nodata -9999 NETCDF:"SD_20160901.nc":snd_upd  out.tif

However I'm not even sure if that ullr is right, or how I would confirm

I was wondering if anyone had a suggestion on how to proceed.

Meta data from the nc is given here:

Spatial sampling:       1 km global cylindrical 1-km Equal-Area Scalable Earth version 2 (EASE-2)
Dimensions:             4500 (rows) x 34703 (columns)
Domain:                 Mountain ranges (Global Mountain Biodiversity Assessment (GMBA) inventory) from 20°N to 90°N in latitude
No data value:          -9999

Many thanks,

Cheers
Chris

Chris Marsh, PhD
University of Saskatchewan
chrismarsh.ca<http://chrismarsh.ca/>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20200201/0cf9b9cc/attachment.html>


More information about the gdal-dev mailing list