geospatialpython/subset.py at master · GeospatialPython/geospatialpython · GitHub

20 lines (19 loc) · 667 Bytes

import shapefile

# Create a reader instance

r = shapefile.Reader("buildingfootprints/Building_Footprint")

# Create a writer instance

w = shapefile.Writer(shapeType=shapefile.POLYGON)

# Copy the fields to the writer

w.fields = list(r.fields)

# Grab the geometry and records from all features

# with the correct county name

selection = []

for rec in enumerate(r.records()):

if rec[1][1].startswith("Hancock"):

selection.append(rec)

# Add the geometry and records to the writer

for rec in selection:

w._shapes.append(r.shape(rec[0]))

w.records.append(rec[1])

# Save the new shapefile

w.save("buildingfootprints/HancockFootprints")

Expand file tree