[Proj] Proj Digest, Vol 72, Issue 4

[massimo costantini]

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.maptools.org/pipermail/proj/attachments/20100517/a36f76fb/attachment.htm