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环境。