[ Web Proxy ]
URL:
Viewing: https://raw.githubusercontent.com/GeospatialPython/geospatialpython/master/subset.py [Back]  [Original]

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") 

Web Proxy Viewer  |  New URL  |  Original Page