[Proj] Proj Digest, Vol 72, Issue 4

[massimo costantini]

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.
Infact:
#>./invproj +proj=stere +lat_0=0.0 +lon_0=0.0
110000 110000
0d59'17.505''E 0d59'40.934''N

#>./proj +proj=stere +lat_0=0.0 +lon_0=0.0
1 1
222627.90 221171.23
0.5 0.5
111318.11 110577.11



On Mon, May 17, 2010 at 11:23 AM, massimo costantini <
massimo.costantini at gmail.com> wrote:

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