<div dir="ltr">Using GDAL 1.11.2 (Anaconda Python osgeo binstar install).<div><br></div><div>I have a WKT projection: <font face="monospace, monospace">'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]]'</font><br><div><br></div><div>I import this into a srs object and attempt to export to Proj4.  This is throwing:</div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">RuntimeError: No translation for Mercator to PROJ.4 format is known.</font><br></div><div><br></div><div>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:</div><div><br></div><div><font face="monospace, monospace">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 '.</font><br></div><div><br></div><div>Which leads me to try:</div><div><br></div><div><div><font face="monospace, monospace">gdalsrsinfo -V '+proj=merc +lon_0=0 +k=0 +x_0=0 +y_0=0 +a=3396190 +b=3376200 +units=m +no_defs'</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">Validate Succeeds</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">PROJ.4 : '+proj=merc +lon_0=0 +k=0 +x_0=0 +y_0=0 +a=3396190 +b=3376200 +units=m +no_defs '</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">OGC WKT :</font></div><div><font face="monospace, monospace">PROJCS["unnamed",</font></div><div><font face="monospace, monospace">    GEOGCS["unnamed ellipse",</font></div><div><font face="monospace, monospace">        DATUM["unknown",</font></div><div><font face="monospace, monospace">            SPHEROID["unnamed",3396190,169.894447223611]],</font></div><div><font face="monospace, monospace">        PRIMEM["Greenwich",0],</font></div><div><font face="monospace, monospace">        UNIT["degree",0.0174532925199433]],</font></div><div><font face="monospace, monospace">    PROJECTION["Mercator_1SP"],</font></div><div><font face="monospace, monospace">    PARAMETER["central_meridian",0],</font></div><div><font face="monospace, monospace">    PARAMETER["scale_factor",0],</font></div><div><font face="monospace, monospace">    PARAMETER["false_easting",0],</font></div><div><font face="monospace, monospace">    PARAMETER["false_northing",0],</font></div></div><div><font face="monospace, monospace">    UNIT["Meter",1]]</font><br></div><div><br></div><div>So gdalsrsinfo is able to validate the proj4 string.  </div><div><br></div><div>Any ideas?</div><div><br></div></div></div>