I'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>#>./invproj +proj=stere +lat_0=0.0 +lon_0=0.0<br>110000 110000<br>0d59'17.505''E 0d59'40.934''N<br>
<br>#>./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"><<a href="mailto:massimo.costantini@gmail.com">massimo.costantini@gmail.com</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;">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"><<a href="mailto:proj-request@lists.maptools.org" target="_blank">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" 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 'help' 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 "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" target="_blank">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" target="_blank">proj@lists.maptools.org</a>><br>
Message-ID: <<a href="mailto:19437.36550.337500.172870@cerise.gclements.plus.com" target="_blank">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" target="_blank">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" target="_blank">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" target="_blank">proj@lists.maptools.org</a>><br>
Message-ID: <<a href="mailto:4BEE893F.2060203@sieczka.org" target="_blank">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" 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>