<div dir="ltr">Hi,<br>I am trying to wrap some EASE-grid data from NSIDC for global datasets. NSIDC says it is based on the Cylindrical Equal Area, and they do provide both an EPSG code for it (3410), as well as a Proj string:<br>
"+proj=cea +lat_0=0 +lon_0=0 +lat_ts=30 +a=6371228.0 +units=m" (see <<a href="http://nsidc.org/data/atlas/epsg_3410.html">http://nsidc.org/data/atlas/epsg_3410.html</a>>. It appears that GDAL supports it, and the relevant line for it is in $GDAL_DATA/pcs.csv:<br>
3410,NSIDC EASE-Grid Global,9001,4053,19869,9834,1,0,8807,0,9001,8823,30,9102,8802,0,9102,8806,0,9001,,,,,,,,,<br><br>So far so good, but warping an image into something else doesn't work:<br>$ gdalwarp --debug on -of GTiff -s_srs "EPSG:3410" -t_srs 'EPSG:4326' INPUT OUTPUT.tif<br>
EPSG: No WKT support for projection method 9834.<br>ERROR 6: EPSG PCS/GCS code 3410 not found in EPSG support files. Is this a valid<br>EPSG coordinate system?<br>ERROR 1: Translating source or target SRS failed:<br>EPSG:3410<br>
<br>$ gdalwarp --debug on -of GTiff -s_srs "+proj=cea +lat_0=0 +lon_0=0 +lat_ts=30 +a=6371228.0 +units=m" -t_srs 'EPSG:4326' INPUT OUTPUT.tif<br>OGR_PROJ4: Can't find ellipse definition in:<br>+proj=cea +lat_0=0 +lon_0=0 +lat_ts=30 +a=6371228.0 +units=m<br>
ERROR 1: Translating source or target SRS failed:<br>+proj=cea +lat_0=0 +lon_0=0 +lat_ts=30 +a=6371228.0 +units=m<br><br>Hmmm.... Does anyone have any ideas around this?<br><br>Thanks!<br>Jose<br><br clear="all"><br>-- <br>
Centre for Terrestrial Carbon Dynamics<br>Department of Geography, University College London<br>Gower Street, London WC1E 6BT, UK<br>
</div>