I got it:<br><br><a href="http://nsidc.org/data/atlas/epsg_3410.html">http://nsidc.org/data/atlas/epsg_3410.html</a><br><br><br><br><div class="gmail_quote">On Wed, Jan 23, 2013 at 11:46 AM, Nikolaos Hatzopoulos <span dir="ltr"><<a href="mailto:nhatzop@gmail.com" target="_blank">nhatzop@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">thanks for the respond,<br><br>Yes you are right it is +proj=cea <br>I got confused with the GCTP codes<br>about the datum=WGS84 I am not quite sure the specs that I find from HDFview are:<br>
<br>Projection=GCTP_CEA<br>ProjParams=(6371228,0,0,0,0,30000000,0,0,0,0,0,0,0)<br>
SphereCode=-1<br><br><br>--Nikos<div class="HOEnZb"><div class="h5"><br><br><div class="gmail_quote">On Wed, Jan 23, 2013 at 10:57 AM, Frank Warmerdam <span dir="ltr"><<a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Nikos,<br>
<br>
I see that GCTP Projection 97 is Cylindrical Equal Area.  This can<br>
have a standard parallel and a central meridian (from<br>
<a href="http://www.remotesensing.org/geotiff/proj_list/cylindrical_equal_area.html" target="_blank">http://www.remotesensing.org/geotiff/proj_list/cylindrical_equal_area.html</a>).<br>
 Without knowing what these are it is hard to do much, but you might<br>
try it assuming they are zero and zero in which case you would use<br>
something like:<br>
<br>
-s_srs "+proj=cea +datum=WGS84"<br>
<br>
If you examine the HDF metadata more carefully you might find more<br>
information.  Initially see what metadata makes it through GDAL with<br>
gdalinfo.  If that doesn't show anything, perhaps try low level HDF<br>
tools.<br>
<br>
If you make the datasource available, and filed a bug it might be<br>
possible for us to fix ogr_srs_usgs.cpp to handle CEA.  In the<br>
comments I see:<br>
<br>
        // FIXME: CEA --- Cylindrical Equal Area skipped (Grid corners<br>
set in meters for EASE grid)<br>
<br>
Best regards,<br>
Frank<br>
<div><div><br>
<br>
On Wed, Jan 23, 2013 at 10:20 AM, Nikolaos Hatzopoulos<br>
<<a href="mailto:nhatzop@gmail.com" target="_blank">nhatzop@gmail.com</a>> wrote:<br>
> I have this projection for my nasa hdf file:<br>
><br>
> Coordinate System is:<br>
> LOCAL_CS["GCTP projection number 97",<br>
>     UNIT["Meter",1]]<br>
> Origin = (-17334194.000000000000000,7344784.825000000186265)<br>
> Pixel Size = (25067.525334417932754,-25067.525298634809587)<br>
><br>
><br>
> I tried to reporject my image using gdalwrap and I got:<br>
><br>
> ERROR 1: No PROJ.4 translation for source SRS, coordinate<br>
> transformation initialization has failed.<br>
><br>
> I am reading this:<br>
> <a href="http://www.gdal.org/ogr/classOGRSpatialReference.html#a04ab13a2ec4ceffd6f802422f95e8f15" target="_blank">http://www.gdal.org/ogr/classOGRSpatialReference.html#a04ab13a2ec4ceffd6f802422f95e8f15</a><br>
><br>
> from the array I understand:<br>
><br>
> 9 Transverse Mercator<br>
> 7 FN<br>
><br>
><br>
><br>
> how I have to define my -s_srs=?????<br>
><br>
> --Nikos Hatzopoulos<br>
><br>
</div></div>> _______________________________________________<br>
> gdal-dev mailing list<br>
> <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
> <a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
<span><font color="#888888"><br>
<br>
<br>
--<br>
---------------------------------------+--------------------------------------<br>
I set the clouds in motion - turn up   | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a><br>
light and sound - activate the windows | <a href="http://pobox.com/%7Ewarmerdam" target="_blank">http://pobox.com/~warmerdam</a><br>
and watch the world go round - Rush    | Geospatial Software Developer<br>
</font></span></blockquote></div><br>
</div></div></blockquote></div><br>