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

[Nikolaos Hatzopoulos]

I got it:

http://nsidc.org/data/atlas/epsg_3410.html



On Wed, Jan 23, 2013 at 11:46 AM, Nikolaos Hatzopoulos <nhatzop at gmail.com>wrote:

> 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/bef1218a/attachment.html>