Luis Torres lptorres
lptorres
/ multi2single.py
Created
December 4, 2013 01:55
A simple script that converts shapefiles with multipart polygons into singlepart polygons, and copies the fields of the source shapefile. This script is based on a script by Paolo Corti, found here: https://gist.github.com/anonymous/735330 The original script only converts multipolygons by creating a new feature for each part of a multipolygon. …
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| import os | |
| from osgeo import gdal, ogr, osr | |
| import sys | |
| def initFields(in_lyr, out_lyr): | |
| #Arbitrarily get the first feature | |
| feat = in_lyr[0] | |
| #loop over each field |
lptorres
/ multi2single_inplace.py
Last active
September 14, 2021 06:16
This script is similar to multi2single.py, but it converts the shapefile in place, instead of creating a new shapefile.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| from osgeo import gdal, ogr | |
| import sys | |
| def multipoly2poly(layer): | |
| for feature in layer: | |
| fid = feature.GetFID() | |
| geom = feature.GetGeometryRef() | |
| if geom.GetGeometryName() == 'MULTIPOLYGON': |