Fixed issue #194 (Bug with the stereographic projection at equator) t…

…hat multiplied equatorial coordinates by 2
OSGeo · Oct 18, 2016 · 93024d9 · 93024d9
1 parent 26607f1
commit 93024d9
Showing 1 changed file with 5 additions and 5 deletions.
diff --git 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[] = {