I have shapely multipolygon, how can I convert that to list of polygons?
polygons = [Polygon(cham_geom), Polygon(neighbor_geom)]
boundary = cascaded_union(polygons) # results in multipolygon sometimes
if boundary.geom_type == 'MultiPolygon':
# extract polygons out of multipolygon