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

Adding a polygon feature to a vector layer

In this recipe, we'll add a polygon to a layer. A polygon is the most complex kind of geometry. However, in QGIS, the API is very similar to a line.

Getting ready

For this recipe, we'll use a simple polygon shapefile, which you can download as a ZIP file from https://geospatialpython.googlecode.com/files/polygon.zip.

Extract this shapefile to a folder called polygon in your /qgis_data directory.

How to do it...

This recipe will follow the standard PyQGIS process of loading a layer, building a feature, and adding it to the layer's data provider, as follows:

  1. Start QGIS.
  2. From the Plugins menu, select Python Console.
  3. First, load the layer and validate it:
    vectorLyr = QgsVectorLayer('/qgis_data/polygon/polygon.shp', 'Polygon' , "ogr")
    vectorLyr.isValid()
    
  4. Next, access the layer's data provider:
    vpr = vectorLyr.dataProvider()
    
  5. Now, build a list of points for the polygon:
    points = []
    points.append(QgsPoint(-123.26,49.06))
    points.append(QgsPoint(-127.19,43.07))
    points.append(QgsPoint(-120.70,35.21))
    points.append(QgsPoint(-115.89,40.02))
    points.append(QgsPoint(-113.04,48.47))
    points.append(QgsPoint(-123.26,49.06))
    
  6. Next, create a geometry object and ingest the points as a polygon. We nest our list of points in another list because a polygon can have inner rings, which will consist of additional lists of points being added to this list:
    poly = QgsGeometry.fromPolygon([points])
    
  7. Next, build the feature object and add the points:
    f = QgsFeature()
    f.setGeometry(poly)
    
  8. Finally, add the feature to the layer's data provider and update the extents:
    vpr.addFeatures([f])
    

How it works...

Adding a polygon is very similar to adding a line, with one key difference that is a common pitfall. The last point must be identical to the first point in order to close the polygon. If you don't repeat the first point, you won't receive any errors, but the polygon will not be displayed in QGIS, which can be difficult to troubleshoot.

主站蜘蛛池模板: 察隅县| 高密市| 马尔康县| 资阳市| 皮山县| 望都县| 武城县| 昌宁县| 抚松县| 巴彦县| 丰镇市| 浙江省| 上犹县| 高密市| 陵川县| 故城县| 张家港市| 娄烦县| 岗巴县| 肇庆市| 威海市| 临安市| 石阡县| 石首市| 丹阳市| 古田县| 连南| 麦盖提县| 潮安县| 额尔古纳市| 噶尔县| 秦安县| 齐河县| 綦江县| 宁阳县| 璧山县| 仁怀市| 白银市| 京山县| 宁明县| 泰宁县|