GIS编程：利用Arcpy实现道格拉斯-普克算法（核心代码已用类包装，复制粘贴即可用）

火卫二德莫斯

已于 2022-05-23 08:58:34 修改

阅读量1.7k

点赞数 5

分类专栏： GIS 道格拉斯算法文章标签：算法 python

于 2022-03-12 10:38:06 首次发布

本文链接：https://blog.csdn.net/weixin_47967519/article/details/123438987

版权

本文介绍了一种通过递归实现的 Douglas-Peucker 算法，用于简化多边形，提供了Python代码示例。用户可以自定义输入shapefile，并设置容差阈值，最终输出抽稀后的polyline。

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

道格拉斯-普克算法是一种通过递归实现折线&多边形抽稀简化的算法，原理很简单，大多数文章都有讲解，这里不多赘述，直接上代码。

# -*- coding:utf-8 -*-
import arcpy
import sys
from arcpy import env
import numpy as np
import math

reload(sys)
sys.setdefaultencoding('utf-8')
env.overwriteOutput = True

# 读取一个shapefile, 用户自定义
# input a shapefile customized by user
shp = r'C:\Users\Lenovo\Desktop\11111\polyline.shp'


# 从.shp文件中获取多段线顶点坐标列表
# access the nodes' coordinate of each single feature
def getPolyXY(shape):
    with arcpy.da.SearchCursor(shape, ["SHAPE@"]) as cursor:
        Plines = []
        for row in cursor:
            coordXY = row[0].getPart()
            record = []
            for v in range(row[0].pointCount - 1):
                pnt = coordXY.getObject(0).getObject(v)
                record.append([pnt.X, pnt.Y])
            Plines.append(record)
        return Plines


class Douglas(object):
    def