[gdal-dev] Standard Parallel 1 - being changed but not displayed?

Jonathan Moules jonathan-lists at lightpear.com
Tue Nov 17 08:55:15 PST 2015


Hi List,
 I have a Geotiff which includes this projection:
 
PROJ.4 : '+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
'


OGC WKT :
PROJCS["Mercator",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Mercator_1SP"],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]




If I load this raster into ArcGIS, it displays in the wrong place (a few thousand kilometres North).


I then run it through gdal_translate (GDAL 1.11.1), with no flags:
 * gdal_trainslate input.tif output.tif
 
For output.tif, GDALSRSInfo shows that the projection is identical, but now the file loads correctly in ArcGIS.
The same file works fine in QGIS both before and after the "translation".


Looking at the projection info in ArcGIS, it displays one difference:
Before (not working):
Standard_parallel_1 = 60


After (working):
Standard_parallel_1 = 0


But I don't see anything about those in either of the GDALSRSInfo outputs.


So my questions:
 - What is gdal_translate doing to the file to "fix" it?
 - If it is something to do with Standard Parallel 1 - why isn't this component of the projection exposed by GDAL?
 
Thoughts welcome.
Thanks,
Jonathan

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20151117/d5cc1844/attachment.html>


More information about the gdal-dev mailing list