下载app送18 88彩金

書名： QGIS Python Programming Cookbook
作者名： Joel Lawhead
本章字數： 354字
更新時間： 2021-07-23 19:48:54

Buffering a feature intermediate

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.

Getting ready

Once again, we'll use the same New York City Museums layer. You can download the layer from https://geospatialpython.googlecode.com/svn/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.

How to do it...

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 feature:
```
ft = fts.next()
```
Select this feature:
```
lyr.setSelectedFeatures([ft.id()])
```
Create the buffer:
```
buff = ft.geometry().buffer(.2,8)
```

Set up a memory layer for the buffer's geometry:

buffLyr = QgsVectorLayer('Polygon?crs=EPSG:4326', 'Buffer' , 'memory')

Access the layer's data provider:
```
pr = buffLyr.dataProvider()
```
Create a new feature:
```
b = QgsFeature()
```
Set the feature's geometry with the buffer geometry:
```
b.setGeometry(buff)
```
Add the feature to the data provider:
```
pr.addFeatures([b])
```
Update the buffer layer's extents:
```
buffLyr.updateExtents()
```
Set the buffer layer's transparency so that you can see other features as well:
```
buffLyr.setLayerTransparency(70)
```

Add the buffer layer to the map:

QgsMapLayerRegistry.instance().addMapLayers([buffLyr])

Verify that your map looks similar to this screenshot:

How it works...

The interesting portion of this recipe starts with Step 6, which creates the buffer geometry. The parameters for the buffer() method are the distance in map units for the buffer followed by the number of straight line segments used to approximate curves. The more segments you specify, the more the buffer appears like a circle. However, more segments equals greater geometric complexity and therefore slower rendering, as well as slower geometry calculations. The other interesting feature of this recipe is Step 13, in which we set the transparency of the layer to 70 percent. We also introduce the way to create a new layer, which is done in memory. Later chapters will go more in depth on creating data.

官术网_书友最值得收藏!

QGIS Python Programming Cookbook

Buffering a feature intermediate

Getting ready

How to do it...

How it works...