diff options
| -rw-r--r-- | src/projections/bonne.cpp | 11 | ||||
| -rw-r--r-- | test/gie/builtins.gie | 25 |
2 files changed, 31 insertions, 5 deletions
diff --git a/src/projections/bonne.cpp b/src/projections/bonne.cpp index 0e9bae79..289eb23d 100644 --- a/src/projections/bonne.cpp +++ b/src/projections/bonne.cpp @@ -26,9 +26,14 @@ static PJ_XY e_forward (PJ_LP lp, PJ *P) { /* Ellipsoidal, forward */ double rh, E, c; rh = Q->am1 + Q->m1 - pj_mlfn(lp.phi, E = sin(lp.phi), c = cos(lp.phi), Q->en); - E = c * lp.lam / (rh * sqrt(1. - P->es * E * E)); - xy.x = rh * sin(E); - xy.y = Q->am1 - rh * cos(E); + if (fabs(rh) > EPS10) { + E = c * lp.lam / (rh * sqrt(1. - P->es * E * E)); + xy.x = rh * sin(E); + xy.y = Q->am1 - rh * cos(E); + } else { + xy.x = 0.; + xy.y = 0.; + } return xy; } diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index 91f0cb22..5f78c979 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -618,7 +618,7 @@ Bonne (Werner lat_1=90) =============================================================================== ------------------------------------------------------------------------------- -operation +proj=bonne +ellps=GRS80 +lat_1=0.5 +lat_2=2 +operation +proj=bonne +ellps=GRS80 +lat_1=0.5 ------------------------------------------------------------------------------- tolerance 0.1 mm accept 2 1 @@ -641,7 +641,18 @@ accept -200 -100 expect -0.001796698 0.499095631 ------------------------------------------------------------------------------- -operation +proj=bonne +R=6400000 +lat_1=0.5 +lat_2=2 +operation +proj=bonne +ellps=GRS80 +lat_1=90 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 90 +expect 0 0 + +direction inverse +accept 0 0 +expect 0 90 + +------------------------------------------------------------------------------- +operation +proj=bonne +R=6400000 +lat_1=0.5 ------------------------------------------------------------------------------- tolerance 0.1 mm accept 2 1 @@ -663,6 +674,16 @@ expect -0.001790562 0.500895246 accept -200 -100 expect -0.001790561 0.499104753 +------------------------------------------------------------------------------- +operation +proj=bonne +R=6400000 +lat_1=90 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 90 +expect 0 0 + +direction inverse +accept 0 0 +expect 0 90 =============================================================================== Cal Coop Ocean Fish Invest Lines/Stations |