Split a GeoJSON MultiPolygon FeatureCollection into GeoJSON Polygons

split-multi.py

# split-multi.py
# open source, MIT license

import json

js = open('multipolygon.geojson', 'r').read()
gj = json.loads(js)

output = { "type": "FeatureCollection", "features": [] }

for feature in gj['features']:
    if (feature['geometry'] is not None) and (feature['geometry']['type'] == 'MultiPolygon'):
        for poly in feature['geometry']['coordinates']:
            xfeature = { "type": "Feature", "properties": {}, "geometry": { "type": "Polygon" } }
            xfeature['geometry']['coordinates'] = poly
            output['features'].append(xfeature)

open('polygons.geojson', 'w').write(json.dumps(output))