[gdal-dev] EASE-grid projection

[Jose Gomez-Dans]

Hi,
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:
"+proj=cea +lat_0=0 +lon_0=0 +lat_ts=30 +a=6371228.0 +units=m" (see <
http://nsidc.org/data/atlas/epsg_3410.html>. It appears that GDAL supports
it, and the relevant line for it is in $GDAL_DATA/pcs.csv:
3410,NSIDC EASE-Grid
Global,9001,4053,19869,9834,1,0,8807,0,9001,8823,30,9102,8802,0,9102,8806,0,9001,,,,,,,,,

So far so good, but warping an image into something else doesn't work:
$ gdalwarp --debug on -of GTiff -s_srs "EPSG:3410" -t_srs 'EPSG:4326' INPUT
OUTPUT.tif
EPSG: No WKT support for projection method 9834.
ERROR 6: EPSG PCS/GCS code 3410 not found in EPSG support files.  Is this a
valid
EPSG coordinate system?
ERROR 1: Translating source or target SRS failed:
EPSG:3410

$ 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
OGR_PROJ4: Can't find ellipse definition in:
+proj=cea +lat_0=0 +lon_0=0 +lat_ts=30 +a=6371228.0 +units=m
ERROR 1: Translating source or target SRS failed:
+proj=cea +lat_0=0 +lon_0=0 +lat_ts=30 +a=6371228.0 +units=m

Hmmm.... Does anyone have any ideas around this?

Thanks!
Jose


-- 
Centre for Terrestrial Carbon Dynamics
Department of Geography, University College London
Gower Street, London WC1E 6BT, UK
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20081010/f7a962fd/attachment.html