Mastering Geospatial Analysis with Python

Creating polygon geometry with OGR

OGR lets you write vector geometries such as points, lines, mulitipoints, multilinestrings, multipolygons and geometry collections. You can give these geometry values in coordinates or meters if you plan to project them later. All geometries you create follow the same procedure, separate points are defined and then strung together as lines or polygons. You define separate entities in numbers, encode them in well-known binary (WKB), and the final polygon is translated to well-known text (WKT). A Jupyter Notebook will return the coordinates of the polygon but won't plot it automatically, for this, we'll use Shapely later in this chapter:

In: from osgeo import ogr    r = ogr.Geometry(ogr.wkbLinearRing) r.AddPoint(1,1) ...