使用wkbMultiPolygon()函数在Python中生成多边形图层

要使用wkbMultiPolygon()函数在Python中生成多边形图层，首先需要导入必要的模块和函数。

首先，导入ogr模块和wkbMultiPolygon函数：

from osgeo import ogr
from osgeo import ogr

接下来，定义一个函数来创建多边形图层并保存到矢量文件中：

def create_multi_polygon_layer(file_path):
    # 创建数据源
    driver = ogr.GetDriverByName('ESRI Shapefile')
    data_source = driver.CreateDataSource(file_path)

    # 创建图层
    layer = data_source.CreateLayer('multi_polygon', geom_type=ogr.wkbMultiPolygon)

    # 添加字段
    field_defn = ogr.FieldDefn('id', ogr.OFTInteger)
    layer.CreateField(field_defn)

    # 创建多边形
    multi_polygon = ogr.Geometry(ogr.wkbMultiPolygon)

    # 添加多边形
    polygon1 = ogr.Geometry(ogr.wkbPolygon)
    ring1 = ogr.Geometry(ogr.wkbLinearRing)
    ring1.AddPoint(0, 0)
    ring1.AddPoint(1, 0)
    ring1.AddPoint(1, 1)
    ring1.AddPoint(0, 1)
    polygon1.AddGeometry(ring1)
    multi_polygon.AddGeometry(polygon1)

    polygon2 = ogr.Geometry(ogr.wkbPolygon)
    ring2 = ogr.Geometry(ogr.wkbLinearRing)
    ring2.AddPoint(2, 2)
    ring2.AddPoint(3, 2)
    ring2.AddPoint(3, 3)
    ring2.AddPoint(2, 3)
    polygon2.AddGeometry(ring2)
    multi_polygon.AddGeometry(polygon2)

    # 创建要素
    feature_defn = layer.GetLayerDefn()
    feature = ogr.Feature(feature_defn)
    feature.SetGeometry(multi_polygon)
    feature.SetField('id', 1)

    # 添加要素
    layer.CreateFeature(feature)

    # 保存数据源和释放资源
    data_source = None

最后，调用这个函数来创建多边形图层并保存到矢量文件中：

create_multi_polygon_layer('multi_polygon.shp')

这个例子创建了一个名为'multi_polygon'的图层，包含两个多边形。每个多边形都由一个线性环组成，每个线性环由一系列点构成。图层还包含一个整数类型的'id'字段。

多边形1的顶点坐标为(0, 0)，(1, 0)，(1, 1)，(0, 1)。

多边形2的顶点坐标为(2, 2)，(3, 2)，(3, 3)，(2, 3)。

最后，将多边形图层保存为'multi_polygon.shp'矢量文件。