[Gdal-dev] gdal_translate adding ESRI projection bug?

[Frank Warmerdam]

joosthoek wrote:
> Hi,
> 
> I use gdal_translate with -a_srs "ESRI::projection.prj", with the following
> projection.prj:
> PROJCS["Mars_Sinusoidal_clon93",GEOGCS["MOLA_MEGDR_Sphere",DATUM["MOLA_MEGDR_Sphere",SPHEROID["MOLA_MEGDR_Sphere",3396000.0,0.0]],PRIMEM["Reference_Meridian",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Sinusoidal"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",93.0],UNIT["Meter",1.0]]
> 
> When I open the resulting file in ArcGIS the Central_Meridian is somehow set
> to zero. I tried multiple formats but all put the image in the zero clon
> location. Somehow the 93 isn't read correctly during the gdal_translate (is
> my conclusion). Could this be a bug?

Jelmer,

What format are you translating to?

I tried "testepsg ESRI::test.prj" with your definition and the central
meridian was properly preserved. So I *suspect* we are either writing
out the wrong sinsoidal definition for ArcGIS or there is some particular
issue with the fileformat through which things are being written.

I would suggest you file a specific ticket on the issue with a small
file demonstrating the problem.

Best regards,
-- 
---------------------------------------+--------------------------------------
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    | President OSGeo, http://osgeo.org