基于OSGEO的两个面状矢量文件空间操作：相交、交集取反、相减和合并

John H.

已于 2022-09-29 10:21:27 修改

阅读量1.3k

点赞数

分类专栏： OSGEO 文章标签： python

于 2022-03-24 09:16:06 首次发布

本文链接：https://blog.youkuaiyun.com/weixin_42778485/article/details/123702359

版权

本文介绍了如何使用OSGeo库在Python环境中进行面状矢量文件的空间操作，包括相交、交集取反、相减和合并等核心功能，详细阐述了相关代码实现和步骤。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

import os, shutil, sys
from osgeo import ogr, gdal, gdal_array, ogr, osr
import math, numpy


def shapefile_geometric_relations_operations(_input_shapefile_1, _input_shapefile_2, result_shapfile, action=1):
    '''
    _input_shapefile_1第一个矢量文件的绝对路径
    _input_shapefile_2第二个矢量文件的绝对路径
    result_shapfile输出矢量文件绝对路径
    action=1表示取两个矢量的交集intersection
          =2表示取两个矢量交集的反集sysdifference
          =3表示两个矢量相减difference
          =4表示取两个矢量的并集union
    这个函数对输入的矢量文件没有坐标系的要求，输出的矢量文件和输入的矢量文件坐标系相同
    '''
    _ds_1, _ds_2 = ogr.Open(_input_shapefile_1), ogr.Open(_input_shapefile_2)
    _lyr_1, _lyr_2 = _ds_1.GetLayer(), _ds_2.GetLayer()
    shpdriver = ogr.GetDriverByName('ESRI Shapefile')
    if os.path.exists(result_shapfile):
        shpdriver.DeleteDataSource(result_shapfile)
    
    outDifference = shpdriver.CreateDataSource(result_shapfile)
    outDifferenceLyr = outDifference