diff options
| author | Glynn Clements <glynn@gclements.plus.com> | 2010-05-14 18:55:26 +0100 |
|---|---|---|
| committer | Elliott Sales de Andrade <quantum.analyst@gmail.com> | 2015-05-29 15:32:39 -0400 |
| commit | 35c0972b0c958d45ad1e3a375be2d3a8b951586a (patch) | |
| tree | 3696cace74be68a85b7c1ac948597173b233afce /src | |
| parent | bc68930723d48f2669d6bd9ebe176614e9dd893d (diff) | |
| download | PROJ-35c0972b0c958d45ad1e3a375be2d3a8b951586a.tar.gz PROJ-35c0972b0c958d45ad1e3a375be2d3a8b951586a.zip | |
Fix inverse stereo projection on 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.
Diffstat (limited to 'src')
| -rw-r--r-- | src/PJ_stere.c | 2 |
1 files changed, 0 insertions, 2 deletions
diff --git a/src/PJ_stere.c b/src/PJ_stere.c index 52ec1df1..a0524627 100644 --- a/src/PJ_stere.c +++ b/src/PJ_stere.c @@ -192,8 +192,6 @@ setup(PJ *P) { /* general initialization */ } 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; |