GDAL Python 过滤Shape Polygon中的面积小于某个阈值的小图斑

原创

于 2023-09-06 11:10:31 发布 · 367 阅读

2 ·

CC 4.0 BY-SA版权

文章标签：

#python #开发语言

该篇文章详细描述了如何使用Python和GDAL库，根据面积大小筛选出地理空间矢量文件中的大图斑，适合数据预处理需求。

# -*- coding: utf-8 -*-
# !/usr/bin/mgdal_env
# @Time : 2023/9/6 9:36
# @Author : Hexk
# 过滤矢量文件中的面积小于某个阈值的小图斑

from osgeo import ogr, gdal, osr
import os


def ShapeFiltratePitch(_input_path, _output_path, _area_threshold):
    """
    过滤POLYGON Shape中的细小图斑，根据面积来过滤。
    :param _input_path: 输入文件路径
    :param _output_path: 输出文件的路径，不包括shp名称
    :param _area_threshold: 设定过滤面积阈值，单位平方米
    :return: None
    """
    driver = ogr.GetDriverByName('ESRI SHAPEFILE')
    gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES")
    gdal.SetConfigOption("SHAPE_ENCODING", "UTF8")

    src_ds