I&#39;ve tryed to use proj and invproj with stereo projection with the code that you send to me, but there is always something wrong. <br>Infact:<br>#&gt;./invproj +proj=stere +lat_0=0.0 +lon_0=0.0<br>110000 110000<br>0d59&#39;17.505&#39;&#39;E 0d59&#39;40.934&#39;&#39;N<br>
<br>#&gt;./proj +proj=stere +lat_0=0.0 +lon_0=0.0<br>1 1<br>222627.90 221171.23<br>0.5 0.5<br>111318.11 110577.11<br><br> <br><br><div class="gmail_quote">On Mon, May 17, 2010 at 11:23 AM, massimo costantini <span dir="ltr">&lt;<a href="mailto:massimo.costantini@gmail.com">massimo.costantini@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Thank you,<br>now it works right.<br><font color="#888888"><br>Massimo Costantini</font><div>
<div></div><div class="h5"><br><div class="gmail_quote">On Sat, May 15, 2010 at 7:00 PM,  <span dir="ltr">&lt;<a href="mailto:proj-request@lists.maptools.org" target="_blank">proj-request@lists.maptools.org</a>&gt;</span> wrote:<br>

<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Send Proj mailing list submissions to<br>
        <a href="mailto:proj@lists.maptools.org" target="_blank">proj@lists.maptools.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://lists.maptools.org/mailman/listinfo/proj" target="_blank">http://lists.maptools.org/mailman/listinfo/proj</a><br>
or, via email, send a message with subject or body &#39;help&#39; to<br>
        <a href="mailto:proj-request@lists.maptools.org" target="_blank">proj-request@lists.maptools.org</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:proj-owner@lists.maptools.org" target="_blank">proj-owner@lists.maptools.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than &quot;Re: Contents of Proj digest...&quot;<br>
<br>
<br>
Today&#39;s Topics:<br>
<br>
   1. Re: Problem with inverse stereo projection (Glynn Clements)<br>
   2. PROJ.4 4.7.1 any time soon? (Maciej Sieczka)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Fri, 14 May 2010 18:56:22 +0100<br>
From: Glynn Clements &lt;<a href="mailto:glynn@gclements.plus.com" target="_blank">glynn@gclements.plus.com</a>&gt;<br>
Subject: Re: [Proj] Problem with inverse stereo projection<br>
To: &quot;PROJ.4 and general Projections Discussions&quot;<br>
        &lt;<a href="mailto:proj@lists.maptools.org" target="_blank">proj@lists.maptools.org</a>&gt;<br>
Message-ID: &lt;<a href="mailto:19437.36550.337500.172870@cerise.gclements.plus.com" target="_blank">19437.36550.337500.172870@cerise.gclements.plus.com</a>&gt;<br>
Content-Type: text/plain; charset=us-ascii<br>
<br>
<br>
massimo costantini wrote:<br>
<br>
&gt; someone can tell me why if I set a center o stereo projection on lat 0.0 and<br>
&gt; lon 0.0 the inverse projection (stereo to geo), return always 0E 0N.<br>
<br>
It works for a sphere, but not an ellipsoid.<br>
<br>
The problem is that the setup function doesn&#39;t initialise all of the<br>
values that it needs for the ellipsoid, inverse, equatorial case.<br>
<br>
In the forward and spherical inverse cases, the equatorial case<br>
(|phi0| &lt; 1e-10) is a simplification of the oblique case, and requires<br>
fewer parameters. However, for the ellipsoid inverse case (e_inverse),<br>
the oblique calculation is used, but the setup function doesn&#39;t<br>
initialise all of the parameters.<br>
<br>
Modifying the setup() function in PJ_stere.c so that EQUIT and OBLIQ<br>
use the same code should fix the problem.<br>
<br>
--- PJ_stere.c  2001-02-07 17:55:05.000000000 +0000<br>
+++ PJ_stere.c.tmp      2010-05-14 18:55:26.000000000 +0100<br>
@@ -195,8 +195,6 @@<br>
                        }<br>
                        break;<br>
                case EQUIT:<br>
-                       P-&gt;akm1 = 2. * P-&gt;k0;<br>
-                       break;<br>
                case OBLIQ:<br>
                        t = sin(P-&gt;phi0);<br>
                        X = 2. * atan(ssfn_(P-&gt;phi0, t, P-&gt;e)) - HALFPI;<br>
<br>
<br>
--<br>
Glynn Clements &lt;<a href="mailto:glynn@gclements.plus.com" target="_blank">glynn@gclements.plus.com</a>&gt;<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Sat, 15 May 2010 13:45:03 +0200<br>
From: Maciej Sieczka &lt;<a href="mailto:msieczka@sieczka.org" target="_blank">msieczka@sieczka.org</a>&gt;<br>
Subject: [Proj] PROJ.4 4.7.1 any time soon?<br>
To: &quot;PROJ.4 and general Projections Discussions&quot;<br>
        &lt;<a href="mailto:proj@lists.maptools.org" target="_blank">proj@lists.maptools.org</a>&gt;<br>
Message-ID: &lt;<a href="mailto:4BEE893F.2060203@sieczka.org" target="_blank">4BEE893F.2060203@sieczka.org</a>&gt;<br>
Content-Type: text/plain; charset=UTF-8; format=flowed<br>
<br>
Hi,<br>
<br>
When can we expect a new PROJ.4 release including fixes to &quot;ETRS89-based<br>
CRSs missing datum definition?&quot; [1] and &quot;lacking towgs84 for Pulkovo<br>
1942(58), Poland&quot;?<br>
<br>
[1]<a href="http://trac.osgeo.org/proj/ticket/11" target="_blank">http://trac.osgeo.org/proj/ticket/11</a><br>
[2]<a href="http://trac.osgeo.org/gdal/ticket/3579" target="_blank">http://trac.osgeo.org/gdal/ticket/3579</a><br>
<br>
Regards,<br>
Maciek<br>
<br>
--<br>
Maciej Sieczka<br>
<a href="http://www.sieczka.org" target="_blank">http://www.sieczka.org</a><br>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
Proj mailing list<br>
<a href="mailto:Proj@lists.maptools.org" target="_blank">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>
<br>
End of Proj Digest, Vol 72, Issue 4<br>
***********************************<br>
</blockquote></div><br>
</div></div></blockquote></div><br>