<div dir="ltr">Super cool.  Thanks!</div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jun 12, 2015 at 10:29 AM, Kyle Shannon <span dir="ltr"><<a href="mailto:kyle@pobox.com" target="_blank">kyle@pobox.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Jay,<br>
<div><div class="h5"><br>
<br>
On Fri, Jun 12, 2015 at 11:21 AM, Kyle Shannon <<a href="mailto:kyle@pobox.com">kyle@pobox.com</a>> wrote:<br>
> Jay,<br>
><br>
><br>
> On Fri, Jun 12, 2015 at 10:59 AM, Jay L. <<a href="mailto:jlaura@asu.edu">jlaura@asu.edu</a>> wrote:<br>
>> Using GDAL 1.11.2 (Anaconda Python osgeo binstar install).<br>
>><br>
>> I have a WKT projection:<br>
>> '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]]'<br>
>><br>
>> I import this into a srs object and attempt to export to Proj4.  This is<br>
>> throwing:<br>
>><br>
>> RuntimeError: No translation for Mercator to PROJ.4 format is known.<br>
><br>
> OGRSpatialReference needs a specific Mercator projection such as:<br>
> Mercator_1SP, Mercator_2SP, etc.  Use Mercator_1SP in your case for<br>
> the PROJECTION parameter.<br>
><br>
>><br>
>> Along the same lines, I am trying to compute pixel to latlon and the<br>
>> inverse.  So I need a transformation object.  No problem setting that up<br>
>> until:<br>
>><br>
>> RuntimeError: Failed to initialize PROJ.4 with `+proj=merc +lon_0=0 +k=0<br>
>> +x_0=0 +y_0=0 +a=3396190 +b=3376200 +units=m +no_defs '.<br>
>><br>
>> Which leads me to try:<br>
>><br>
>> gdalsrsinfo -V '+proj=merc +lon_0=0 +k=0 +x_0=0 +y_0=0 +a=3396190 +b=3376200<br>
>> +units=m +no_defs'<br>
>><br>
>> Validate Succeeds<br>
>><br>
>> PROJ.4 : '+proj=merc +lon_0=0 +k=0 +x_0=0 +y_0=0 +a=3396190 +b=3376200<br>
>> +units=m +no_defs '<br>
>><br>
>> OGC WKT :<br>
>> PROJCS["unnamed",<br>
>>     GEOGCS["unnamed ellipse",<br>
>>         DATUM["unknown",<br>
>>             SPHEROID["unnamed",3396190,169.894447223611]],<br>
>>         PRIMEM["Greenwich",0],<br>
>>         UNIT["degree",0.0174532925199433]],<br>
>>     PROJECTION["Mercator_1SP"],<br>
>>     PARAMETER["central_meridian",0],<br>
>>     PARAMETER["scale_factor",0],<br>
>>     PARAMETER["false_easting",0],<br>
>>     PARAMETER["false_northing",0],<br>
>>     UNIT["Meter",1]]<br>
>><br>
>> So gdalsrsinfo is able to validate the proj4 string.<br>
><br>
> gdalsrsinfo can't validate the WKT though, without the Mercator_1SP<br>
><br>
>><br>
>> Any ideas?<br>
>><br>
>><br>
>> _______________________________________________<br>
>> gdal-dev mailing list<br>
>> <a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
>> <a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
><br>
><br>
><br>
> --<br>
> Kyle<br>
<br>
</div></div>Mercator is a valid key in the ESRI WKT dialect:<br>
<br>
Type "copyright", "credits" or "license()" for more information.<br>
>>> from osgeo import osr<br>
>>> wkt = '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]]'<br>
>>> srs1 = osr.SpatialReference()<br>
>>> srs1.ImportFromWkt(wkt)<br>
0<br>
>>> srs1.ExportToProj4()<br>
''<br>
>>> srs1.MorphFromESRI()<br>
0<br>
>>> srs1.ExportToProj4()<br>
'+proj=merc +lon_0=180 +k=1 +x_0=0 +y_0=0 +a=1737400 +b=1737400<br>
+units=m +no_defs '<br>
<br>
The osr.SpatialReference.MorphFromESRI() should do the job too.<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
Kyle<br>
</font></span></blockquote></div><br></div>