[gdal-dev] GCTP projection -s_srs=???

[Nikolaos Hatzopoulos]

thanks for the respond,

Yes you are right it is +proj=cea
I got confused with the GCTP codes
about the datum=WGS84 I am not quite sure the specs that I find from
HDFview are:

Projection=GCTP_CEA
ProjParams=(6371228,0,0,0,0,30000000,0,0,0,0,0,0,0)
SphereCode=-1


--Nikos

On Wed, Jan 23, 2013 at 10:57 AM, Frank Warmerdam <warmerdam at pobox.com>wrote:

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