Mastering Geospatial Analysis with Python

Writing polygon data to a newly created shapefile

Our current polygon only exists in memory. We can create a new shapefile and write the polygon geometry we created earlier to this shapefile. The script consists of the following steps:

1. Import the modules and set the spatial reference (in this case, World Geodetic System 1984 (WGS1984)).
2. Create the shapefile, then the layer using polygon geometry. Next, the geometry is put inside a feature and the feature in a layer. Notice that the script directly references the polygon from the earlier example.
3. The catch is to use the right geometry type in the first line of code, which in this case should be wkbPolygon.
4. The polygon geometry from our earlier example is referenced in this step and put into ...