Workaround for newer GEOS versions.

See also #484
matplotlib · Aug 4, 2020 · a9c274e · a9c274e
1 parent 0216074
commit a9c274e
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/lib/mpl_toolkits/basemap/__init__.py b/lib/mpl_toolkits/basemap/__init__.py
@@ -1459,6 +1459,10 @@ def _readboundarydata(self,name,as_polygons=False):
                         # convert polygons to line segments
                         poly = _geoslib.LineString(poly.boundary)
                     else:
+                        # this is a workaround to avoid
+                        # GEOS_ERROR: CGAlgorithmsDD::orientationIndex encountered NaN/Inf numbers
+                        b[np.isposinf(b)] = 1e20
+                        b[np.isneginf(b)] = -1e20
                         poly = Shape(b)
                     # this is a workaround to avoid
                     # "GEOS_ERROR: TopologyException: