[gdal-dev] Proj4 Transformation Error 6

Jay L. jlaura at asu.edu
Fri Jun 12 09:59:21 PDT 2015


Using GDAL 1.11.2 (Anaconda Python osgeo binstar install).

I have a WKT projection:
'PROJCS["Mercator",GEOGCS["GCS_Moon_2000",DATUM["D_Moon_2000",SPHEROID["Moon_2000_IAU_IAG",1737400.0,0.0]],PRIMEM["Reference_Meridian",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",180.0],PARAMETER["Standard_Parallel_1",0.0],UNIT["Meter",1.0]]'

I import this into a srs object and attempt to export to Proj4.  This is
throwing:

RuntimeError: No translation for Mercator to PROJ.4 format is known.

Along the same lines, I am trying to compute pixel to latlon and the
inverse.  So I need a transformation object.  No problem setting that up
until:

RuntimeError: Failed to initialize PROJ.4 with `+proj=merc +lon_0=0 +k=0
+x_0=0 +y_0=0 +a=3396190 +b=3376200 +units=m +no_defs '.

Which leads me to try:

gdalsrsinfo -V '+proj=merc +lon_0=0 +k=0 +x_0=0 +y_0=0 +a=3396190
+b=3376200 +units=m +no_defs'

Validate Succeeds

PROJ.4 : '+proj=merc +lon_0=0 +k=0 +x_0=0 +y_0=0 +a=3396190 +b=3376200
+units=m +no_defs '

OGC WKT :
PROJCS["unnamed",
    GEOGCS["unnamed ellipse",
        DATUM["unknown",
            SPHEROID["unnamed",3396190,169.894447223611]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Mercator_1SP"],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",0],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]

So gdalsrsinfo is able to validate the proj4 string.

Any ideas?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20150612/f015b3bb/attachment.html>


More information about the gdal-dev mailing list