Thank you,<br>now it works right.<br><br>Massimo Costantini<br><div class="gmail_quote">On Sat, May 15, 2010 at 7:00 PM, <span dir="ltr"><<a href="mailto:proj-request@lists.maptools.org">proj-request@lists.maptools.org</a>></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">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 'help' to<br>
<a href="mailto:proj-request@lists.maptools.org">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">proj-owner@lists.maptools.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of Proj digest..."<br>
<br>
<br>
Today'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 <<a href="mailto:glynn@gclements.plus.com">glynn@gclements.plus.com</a>><br>
Subject: Re: [Proj] Problem with inverse stereo projection<br>
To: "PROJ.4 and general Projections Discussions"<br>
<<a href="mailto:proj@lists.maptools.org">proj@lists.maptools.org</a>><br>
Message-ID: <<a href="mailto:19437.36550.337500.172870@cerise.gclements.plus.com">19437.36550.337500.172870@cerise.gclements.plus.com</a>><br>
Content-Type: text/plain; charset=us-ascii<br>
<br>
<br>
massimo costantini wrote:<br>
<br>
> someone can tell me why if I set a center o stereo projection on lat 0.0 and<br>
> 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'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| < 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'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->akm1 = 2. * P->k0;<br>
- break;<br>
case OBLIQ:<br>
t = sin(P->phi0);<br>
X = 2. * atan(ssfn_(P->phi0, t, P->e)) - HALFPI;<br>
<br>
<br>
--<br>
Glynn Clements <<a href="mailto:glynn@gclements.plus.com">glynn@gclements.plus.com</a>><br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Sat, 15 May 2010 13:45:03 +0200<br>
From: Maciej Sieczka <<a href="mailto:msieczka@sieczka.org">msieczka@sieczka.org</a>><br>
Subject: [Proj] PROJ.4 4.7.1 any time soon?<br>
To: "PROJ.4 and general Projections Discussions"<br>
<<a href="mailto:proj@lists.maptools.org">proj@lists.maptools.org</a>><br>
Message-ID: <<a href="mailto:4BEE893F.2060203@sieczka.org">4BEE893F.2060203@sieczka.org</a>><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 "ETRS89-based<br>
CRSs missing datum definition?" [1] and "lacking towgs84 for Pulkovo<br>
1942(58), Poland"?<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">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>