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

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.

主站蜘蛛池模板: 临城县| 扬州市| 黄陵县| 南阳市| 马公市| 开平市| 泰州市| 惠来县| 集安市| 台江县| 六安市| 崇礼县| 平凉市| 上高县| 花莲县| 揭阳市| 遵义县| 丁青县| 泗阳县| 栖霞市| 朔州市| 蚌埠市| 会东县| 洛南县| 平顶山市| 宁武县| 淮南市| 沛县| 临高县| 小金县| 祁东县| 肥西县| 若尔盖县| 得荣县| 镇平县| 晋城| 米林县| 德江县| 胶州市| 西城区| 新河县|