1 files changed, 9 insertions, 1 deletions
diff --git a/src/geocent.cpp b/src/geocent.cpp
index c023bdd3..dd14af35 100644
--- a/src/geocent.cpp
+++ b/src/geocent.cpp
@@ -399,7 +399,15 @@ void pj_Convert_Geocentric_To_Geodetic (GeocentricInfo *gi,
  */
     CT = Z/RR;
     ST = P/RR;
-    RX = 1.0/sqrt(1.0-gi->Geocent_e2*(2.0-gi->Geocent_e2)*ST*ST);
+    const double denominator = 1.0-gi->Geocent_e2*(2.0-gi->Geocent_e2)*ST*ST;
+    if( denominator == 0 )
+    {
+        *Latitude = HUGE_VAL;
+        *Longitude = HUGE_VAL;
+        *Height = HUGE_VAL;
+        return;
+    }
+    RX = 1.0/sqrt(denominator);
     CPHI0 = ST*(1.0-gi->Geocent_e2)*RX;
     SPHI0 = CT*RX;
     iter = 0;