QGIS根据行政区边界,批量裁剪矢量要素

要使用QGIS编程实现根据行政区边界批量裁剪矢量要素,可以借助PyQGIS来完成。简单的示例代码:


from qgis.core import QgsVectorLayer, QgsProject, QgsGeometry, QgsFeature

# 设置行政区边界文件路径
boundary_file = '/path/to/admin_boundary.shp'
# 设置待裁剪的矢量要素文件夹路径
input_folder = '/path/to/input_features/'
# 设置输出文件夹路径
output_folder = '/path/to/output_features/'

# 加载行政区边界图层
boundary_layer = QgsVectorLayer(boundary_file, 'Admin Boundary', 'ogr')

# 遍历待裁剪的矢量要素文件夹中的所有文件
for filename in os.listdir(input_folder):
    if filename.endswith('.shp'):  # 只处理Shapefile格式的文件
        input_file = os.path.join(input_folder, filename)
        output_file = os.path.join(output_folder, filename)
        
        # 加载待裁剪的矢量要素图层
        input_layer = QgsVectorLayer(input_file, 'Input Features', 'ogr')
        
        # 创建输出图层,与输入图层具有相同的属性表结构
        output_layer = QgsVectorLayer(input_layer.source(), 'Output Features', 'memory')
        output_provider = output_layer.dataProvider()
        
        # 将行政区边界限制为当前待裁剪的矢量要素的空间范围
        boundary_layer.setSubsetString(input_layer.extent().asWktPolygon())
        
        # 对待裁剪的矢量要素进行裁剪
        for feature in input_layer.getFeatures():
            # 获取矢量要素的空间几何形状
            geom = feature.geometry()
            
            # 判断矢量要素是否与行政区边界相交
            if not geom.intersects(boundary_layer.extent()):
                continue
            
            # 裁剪矢量要素
            clipped_geom = geom.intersection(boundary_layer.extent())
            
            # 创建裁剪后的要素
            clipped_feature = QgsFeature(feature)
            clipped_feature.setGeometry(clipped_geom)
            
            # 将裁剪后的要素添加到输出图层中
            output_provider.addFeature(clipped_feature)
        
        # 保存输出图层为Shapefile文件
        QgsVectorFileWriter.writeAsVectorFormat(output_layer, output_file, 'utf-8', output_layer.crs(), 'ESRI Shapefile')
 

请根据实际情况修改代码中的文件路径。注意,该代码针对Shapefile格式的数据进行裁剪,如果输入数据是其他格式,则需要相应调整相关代码。此外,还需要确保已经安装并正确配置了QGIS Python环境。

你可能感兴趣的:(QGIS,QGIS)