基于OSGEO对单个矢量文件的多个feature（图斑）分别导出（不使用arcpy）

John H.

已于 2022-09-29 10:26:06 修改

阅读量229

点赞数

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

于 2022-09-23 11:39:34 首次发布

本文链接：https://blog.csdn.net/weixin_42778485/article/details/127007812

版权

本文介绍如何利用OSGEO库在Python中操作单个矢量文件，无需使用arcpy，对文件中的多个feature进行独立导出。详细步骤和代码示例帮助理解这一非arcpy解决方案。

摘要由CSDN通过智能技术生成

def dispatchShapefileSplit(dispatchShapefile):
    """
    把一个扰动图斑矢量文件包含多个图斑的进行拆分
    比如原图斑名字为name，包括3个图斑
    拆分后的文件名字分别为name_1, name_2, name_3
    """
    outputPath = os.path.dirname(dispatchShapefile)
    name = os.path.basename(dispatchShapefile).split('.')[0]

    driver = ogr.GetDriverByName('ESRI Shapefile')

    inds = ogr.Open(dispatchShapefile, 0)
    if inds is None:
        print('Could not open this shpfile, named as ', os.path.basename(dispatchShapefile), '.')
        return

    inlayer = inds.GetLayer()
    spatial = inlayer.GetSpatialRef()

    num_features = inlayer.GetFeatureCount()
    if num_features <= 1:
        return
    else:
        feature = inlayer.GetNextFeature()
        i = 0
        while feature:
            feat = inlayer.GetFeature(i)
            geom = feat.geometry()
            output_shpfile = os