【677】道格拉斯-普克 抽稀算法 simplify

参考：wikipedia 道格拉斯-普克算法

参考：https://www.osgeo.cn/pygis/shapely-sa.html

　　该方法可以在 shapely 里面通过 simplify 来实现

 

　　算法的基本思路是：对每一条曲线的首末点虚连一条直线，求所有点与直线的距离，并找出最大距离值dmax，用dmax与限差D相比：若dmax <D，这条曲线上的中间点全部舍去；若dmax ≥D，保留dmax 对应的坐标点，并以该点为界，把曲线分为两部分，对这两部分重复使用该方法。

 

 

tolerance 表示的距离，因此根据不同的图形会有所变化

object.simplify(tolerance, preserve_topology=True)

p = Point(0.0, 0.0)
>>> x = p.buffer(1.0)
>>> x.area
3.1365484905459398
>>> len(x.exterior.coords)
65
>>> s = x.simplify(0.05, preserve_topology=False)
>>> s.area
3.061467458920719
>>> len(s.exterior.coords)
17