diff options
| author | Even Rouault <even.rouault@mines-paris.org> | 2019-04-02 11:53:52 +0200 |
|---|---|---|
| committer | GitHub <noreply@github.com> | 2019-04-02 11:53:52 +0200 |
| commit | 3272718cdc48af4a97e5c6b41c50ad3aeffcb81e (patch) | |
| tree | 41062cd523b5e699d5097999f0cc02ff4a150bd3 /src/projections/bonne.cpp | |
| parent | b8a62522b0eca6b4b806471cdff984a28ea1c8fe (diff) | |
| parent | d3eada4ec2a1b825aaacc1e3ad9fc1cd4da09ee5 (diff) | |
| download | PROJ-3272718cdc48af4a97e5c6b41c50ad3aeffcb81e.tar.gz PROJ-3272718cdc48af4a97e5c6b41c50ad3aeffcb81e.zip | |
Merge pull request #1401 from rouault/ossfuzz_14003_14010
Ossfuzz 14003 14010
Diffstat (limited to 'src/projections/bonne.cpp')
| -rw-r--r-- | src/projections/bonne.cpp | 11 |
1 files changed, 8 insertions, 3 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; } |