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

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:

  1. First, load the layer:
    lyr = QgsVectorLayer("/qgis_data/nyc/NYC_MUSEUMS_GEO.shp", "Museums", "ogr")
    
  2. Next, visualize the layer on the map:
    QgsMapLayerRegistry.instance().addMapLayers([lyr])
    
  3. Access the layer's features:
    fts = lyr.getFeatures()
    
  4. Grab the first feature:
    ft = fts.next()
    
  5. Select this feature:
    lyr.setSelectedFeatures([ft.id()])
    
  6. Create the buffer:
    buff = ft.geometry().buffer(.2,8)
    
  7. Set up a memory layer for the buffer's geometry:
    buffLyr = QgsVectorLayer('Polygon?crs=EPSG:4326', 'Buffer' , 'memory')
    
  8. Access the layer's data provider:
    pr = buffLyr.dataProvider()
    
  9. Create a new feature:
    b = QgsFeature()
    
  10. Set the feature's geometry with the buffer geometry:
    b.setGeometry(buff)
    
  11. Add the feature to the data provider:
    pr.addFeatures([b])
    
  12. Update the buffer layer's extents:
    buffLyr.updateExtents()
    
  13. Set the buffer layer's transparency so that you can see other features as well:
    buffLyr.setLayerTransparency(70)
    
  14. Add the buffer layer to the map:
    QgsMapLayerRegistry.instance().addMapLayers([buffLyr])
    

Verify that your map looks similar to this screenshot:

How to do it...

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.

主站蜘蛛池模板: 邯郸县| 东城区| 伊宁市| 永寿县| 双柏县| 周口市| 永吉县| 易门县| 大丰市| 广汉市| 彭水| 镇平县| 南昌县| 广西| 长沙市| 怀柔区| 巩义市| 藁城市| 佛山市| 靖远县| 泰宁县| 确山县| 扶绥县| 湄潭县| 古浪县| 新兴县| 织金县| 二手房| 盈江县| 礼泉县| 卢氏县| 北票市| 通辽市| 乌兰浩特市| 衢州市| 望奎县| 东乡族自治县| 青冈县| 比如县| 永康市| 博白县|