From 683d3097ff2cabd573a82757e8bef6d8f0447d37 Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Sun, 23 Aug 2020 18:11:43 +0200 Subject: lcc: avoid harmless division by zero in untypical case. Fixes https://bugs.chromium.org/p/oss-fuzz/issues/detail?id=24442 --- src/projections/lcc.cpp | 5 +++++ 1 file changed, 5 insertions(+) (limited to 'src') diff --git a/src/projections/lcc.cpp b/src/projections/lcc.cpp index b78b439c..91ffc511 100644 --- a/src/projections/lcc.cpp +++ b/src/projections/lcc.cpp @@ -140,6 +140,11 @@ PJ *PROJECTION(lcc) { Q->n = log(cosphi / cos(Q->phi2)) / log(tan(M_FORTPI + .5 * Q->phi2) / tan(M_FORTPI + .5 * Q->phi1)); + if( Q->n == 0 ) { + // Likely reason is that phi1 / phi2 are too close to zero. + // Can be reproduced with +proj=lcc +a=1 +lat_2=.0000001 + return pj_default_destructor(P, PJD_ERR_CONIC_LAT_EQUAL); + } Q->c = cosphi * pow(tan(M_FORTPI + .5 * Q->phi1), Q->n) / Q->n; Q->rho0 = (fabs(fabs(P->phi0) - M_HALFPI) < EPS10) ? 0. : Q->c * pow(tan(M_FORTPI + .5 * P->phi0), -Q->n); -- cgit v1.2.3