multipolygon fixes for issue #3

mander/districts.py

@@ -1,15 +1,14 @@
 import geopandas as gpd
+import utils
 
+class District(object):
 
-class district(object):
-
-    def __init__(self, pathToGeojsonOrShapefile, epsg='2163'):
+    def __init__(self, pathToGeojsonOrShapefile, epsg='2163', **kwargs):
         self.epsg = epsg
         self.gdf = gpd.read_file(pathToGeojsonOrShapefile)
         if not self.gdf.crs:
             self.gdf.crs = {'init': 'epsg:4326'}
         self.gdf = self.gdf.to_crs(epsg=self.epsg)
         self.area = self.gdf.area.values[0]
         self.perimeter = self.gdf.length.values[0]
-        self.x, self.y = self.gdf.geometry[0].exterior.coords.xy
-        self.coordPairs = zip(self.x, self.y)
+        self.coordPairs = utils.getCoordPairs(self.gdf)

mander/utils.py

@@ -9,3 +9,17 @@ def getMinBoundingCircleArea(district):
     x, y, radius = mbc.make_circle(coords)
     area = math.pi * radius**2
     return area
+
+
+def getCoordPairs(gdf):
+
+    coordPairs = []
+    if gdf.geom_type[0] == 'Polygon':
+        x, y = gdf.geometry[0].exterior.coords.xy
+        coordPairs = zip(x, y)
+    elif gdf.geom_type[0] == 'MultiPolygon':
+        for polygon in gdf.geometry[0]:
+            x, y = polygon.exterior.coords.xy
+            coordPairs += zip(x, y)
+    return coordPairs
+