
QGIS Python Programming Cookbook, Second Edition
By :

Buffering a feature creates a polygon around a feature as a selection geometry or just a simple visualization. In this recipe, we'll buffer a point in a point feature and add the returned polygon geometry to the map.
Once again, we'll use the same New York City Museums layer. You can download the layer from https://github.com/GeospatialPython/Learn/raw/master/NYC_MUSEUMS_GEO.zip.
Unzip that file and place the shapefile's contents in a directory named nyc
within your qgis_data
directory, within your root or home directory.
This recipe involves both a spatial operation and multiple visualizations. To do this, perform the following steps:
First, load the layer:
lyr = QgsVectorLayer("/qgis_data/nyc/NYC_MUSEUMS_GEO.shp", "Museums", "ogr")
Next, visualize the layer on the map:
QgsMapLayerRegistry.instance().addMapLayers([lyr])
Access the layer's features:
fts = lyr.getFeatures()
Grab the first...