add script to extract districts from units geojson