For Mars use its ellipsoid axis:  +a=3396200 +b=3376189<br>or the mean sphere: +R=3389508<br><br>Data from <a href="http://www.hnsky.org/iau-iag.htm">http://www.hnsky.org/iau-iag.htm</a><br><br>Regards<br><br>José Gonçalves<br>
<br><br><br><br><div class="gmail_quote">Em 10 de abril de 2012 16:16, Fischer, Robert P. (GISS-6110)[COLUMBIA UNIVERSITY] <span dir="ltr">&lt;<a href="mailto:robert.p.fischer-1@nasa.gov">robert.p.fischer-1@nasa.gov</a>&gt;</span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">How do I change the radius of the (spherical) Earth in proj.4?  I&#39;m using the proj.4 string:<br>
<br>
    &quot; +proj=laea +lon_0=-40.000000 +lat_0=74.000000 +x_0=1000000 +y_0=1700000 +ellps=sphere&quot;<br>
<br>
..and have discovered that proj.4 assumes a radius of exactly 6370997m (from pj_ellps.c):<br>
<br>
    &quot;sphere&quot;,   &quot;a=6370997.0&quot;,  &quot;b=6370997.0&quot;, &quot;Normal Sphere (r=6370997)&quot;,<br>
<br>
We are running code on a spherical Earth, with a radius (slightly) different from that above.  For consistency sake, proj.4 will need to use the same Earth radius as the code that calls it.  How do I set up a projection on a spherical Earth with a custom radius?<br>

<br>
Also... if someone wanted to use proj.4 for mapping on Mars, how would they do that?<br>
<br>
Thanks,<br>
-- Bob<br>
_______________________________________________<br>
Proj mailing list<br>
<a href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a><br>
<a href="http://lists.maptools.org/mailman/listinfo/proj" target="_blank">http://lists.maptools.org/mailman/listinfo/proj</a><br>
</blockquote></div><br>