diff options
| author | Kristian Evers <kristianevers@gmail.com> | 2016-10-18 14:19:44 +0200 |
|---|---|---|
| committer | Kristian Evers <kristianevers@gmail.com> | 2016-10-18 14:19:44 +0200 |
| commit | 93024d9c167ef503a62992a3f266ba377536398f (patch) | |
| tree | 3fdd10c6976ed04823e27ee370ecc6ddcdd71034 /src | |
| parent | 26607f19d9a8c1442bf2c30744fd1a11a4791c25 (diff) | |
| download | PROJ-93024d9c167ef503a62992a3f266ba377536398f.tar.gz PROJ-93024d9c167ef503a62992a3f266ba377536398f.zip | |
Fixed issue #194 (Bug with the stereographic projection at equator) that multiplied equatorial coordinates by 2
Diffstat (limited to 'src')
| -rw-r--r-- | src/PJ_stere.c | 10 |
1 files changed, 5 insertions, 5 deletions
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[] = { |