From 93024d9c167ef503a62992a3f266ba377536398f Mon Sep 17 00:00:00 2001 From: Kristian Evers Date: Tue, 18 Oct 2016 14:19:44 +0200 Subject: Fixed issue #194 (Bug with the stereographic projection at equator) that multiplied equatorial coordinates by 2 --- src/PJ_stere.c | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) (limited to 'src') diff --git a/src/PJ_stere.c b/src/PJ_stere.c index cc1773da..cfea73d8 100644 --- a/src/PJ_stere.c +++ b/src/PJ_stere.c @@ -52,7 +52,7 @@ static XY e_forward (LP lp, PJ *P) { /* Ellipsoidal, forward */ goto xmul; /* but why not just xy.x = A * cosX; break; ? */ case EQUIT: - A = 2. * Q->akm1 / (1. + cosX * coslam); + A = Q->akm1 / (1. + cosX * coslam); xy.y = A * sinX; xmul: xy.x = A * cosX; @@ -320,10 +320,10 @@ int pj_stere_selftest (void) { }; XY e_fwd_expect[] = { - { 445289.70910023432, 221221.76694834774}, - { 445289.70910023432, -221221.76694835056}, - {-445289.70910023432, 221221.76694834774}, - {-445289.70910023432, -221221.76694835056}, + { 222644.8545501172, 110610.8834741739}, + { 222644.8545501172, -110610.8834741739}, + {-222644.8545501172, 110610.8834741739}, + {-222644.8545501172, -110610.8834741739}, }; XY s_fwd_expect[] = { -- cgit v1.2.3