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

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.

主站蜘蛛池模板: 怀安县| 沅江市| 河西区| 来安县| 松阳县| 大同市| 两当县| 綦江县| 海门市| 谢通门县| 泾阳县| 建德市| 武胜县| 敖汉旗| 南充市| 武安市| 通许县| 同仁县| 呼玛县| 方城县| 涞源县| 涞水县| 东海县| 沾化县| 达州市| 四川省| 兴文县| 雷波县| 全州县| 绵竹市| 鄂州市| 衢州市| 南皮县| 宣武区| 东乡| 顺义区| 县级市| 紫阳县| 荣昌县| 潜山县| 礼泉县|