前言

本程序主要应用于边界线上点的稀释，本文提供随机数和交互式输入两种方式实现数据载入。

一、算法实现

# -*- coding = utf-8 -*-
# @time:2021/6/5 21:39
# Author:scw
# @File:DP.py
# @Software:PyCharm
# Environment:Python3.6
from math import sqrt
import numpy as np
import random as r_seed = 10def PLDistance(Point, Point_Beg, Point_End):p, begP, endP = np.array(Point), np.array(Point_Beg), np.array(Point_End)area = abs(np.cross(endP - begP, p - begP))bottom = sqrt((endP[0] - begP[0]) ** 2 + (endP[1] - begP[1]) ** 2)return area / bottomclass DP_S(object):def __init__(self):self.dslimit = dslimitself.qualify_list = list()self.disqualify_list = list()def dp(self, point_list):if len(point_list) < 3:self.qualify_list.extend(point_list[::-1])else:maxds_index, maxds = 0, 0for index, point in enumerate(point_list):if index in [0, len(point_list) - 1]:continuedistance = PLDistance(point, point_list[0], point_list[-1])if distance > maxds:maxds_index = indexmaxds = distanceif maxds < self.dslimit:self.qualify_list.append(point_list[-1])self.qualify_list.append(point_list[0])else:sequence_Beg = point_list[:maxds_index]sequence_End = point_list[maxds_index:]for sequence in [sequence_Beg, sequence_End]:if len(sequence) < 3 and sequence == sequence_End:self.qualify_list.extend(sequence[::-1])else:self.disqualify_list.append(sequence)def main(self, point_list):self.dp(point_list)while len(self.disqualify_list) > 0:self.dp(self.disqualify_list.pop())print(self.qualify_list)print(len(self.qualify_list))if __name__ == '__main__':dslimit = eval(input("请输入容差:"))ptcnt1 = eval(input("请输入点的个数:"))point_list = [[0] * ptcnt1] * ptcnt1ls = [[0] * ptcnt1] * ptcnt1flag = input("是否使用随机数产生坐标值：[Y/N]")if flag in ['y', 'Y']:print("随机数种子为:{}".format(_seed))r.seed(_seed)for i in range(ptcnt1):ls[i] = round(r.uniform(1, 100), 2)ls = sorted(ls, reverse=True)for i in range(ptcnt1):point_list[i] = [ls[i], round(r.uniform(1, 100), 2)]# point_list[i] = [float(j) for j in point_list[i]]else:print("请输入坐标：[同行数字用空格分隔，不同点则用回车换行]")for i in range(ptcnt1):point_list[i] = list(map(float, input().split(" ")))point_list.sort(key=lambda x: x[0], reverse=True)d = DP_S()d.main(point_list)

二、输出结果

（1）随机数方式

请输入容差:10
请输入点的个数:20
是否使用随机数产生坐标值：[Y/N]y
随机数种子为:10
[[5.41, 84.56], [21.4, 7.43], [25.75, 44.17], [33.45, 32.12], [38.78, 89.92], [43.46, 37.25], [46.23, 30.98], [52.55, 94.16], [57.57, 14.26], [58.23, 1.4], [60.72, 5.38], [67.82, 96.97], [81.52, 98.26], [86.16, 14.16], [95.33, 66.52], [99.66, 68.9]]
16

（2）交互式输入

请输入容差:2
请输入点的个数:8
是否使用随机数产生坐标值：[Y/N]n
请输入坐标：[同行数字用空格分隔，不同点则用回车换行]
0 0
1 1
2 0
4 1
6 5
7 1
10 0
20 1
[[0.0, 0.0], [6.0, 5.0], [7.0, 1.0], [20.0, 1.0]]
4

道格拉斯普克(Douglas-Peuker)算法python实现相关推荐

道格拉斯-普克 Douglas-Peuker(DP算法) python java实现
1.道格拉斯-普克抽稀算法说明道格拉斯-普克抽稀算法是用来对大量冗余的图形数据点进行压缩以提取必要的数据点. 该算法实现抽稀的过程是: 1)对曲线的首末点虚连一条直线,求曲线上所有点与直线的距离,并 ...
【Java代码】道格拉斯-普克 Douglas-Peucker 抽稀算法（算法流程图解+使用JDK8方法实现+详细注解源码）
1.算法说明道格拉斯-普克算法 Douglas-Peucker Algorithm 简称 D-P 算法,亦称为拉默-道格拉斯-普克算法.迭代适应点算法.分裂与合并算法,是将曲线近似表示为一系列点 ...
2021-03-20 GPS抽稀之道格拉斯-普克（Douglas-Peuker)算法
GPS抽稀之道格拉斯-普克(Douglas-Peuker)算法道格拉斯-普克算法是我们常用的一种轨迹点的抽稀算法,抽稀出来的点可以尽可能的维持原先轨迹点的大体轮廓,剔除一些非必要的点. 道格拉斯-普 ...
GPS定位轨迹抽稀之道格拉斯-普克（Douglas-Peuker)算法详解
目录 1. 抽稀 2. 原理 3. 具体思路 4. 代码示例道格拉斯-普克算法是我们常用的一种轨迹点的抽稀算法,抽稀出来的点可以尽可能的维持原先轨迹点的大体轮廓,剔除一些非必要的点 2. 原理假设 ...
道格拉斯-普克 Douglas-Peuker(DP算法)
道格拉斯-普克抽稀算法,是用来对大量冗余的图形数据点进行压缩以提取必要的数据点. 该算法实现抽稀的过程是: 1)对曲线的首末点虚连一条直线,求曲线上所有点与直线的距离,并找出最大距离值dmax,用dm ...
道格拉斯—普克(Douglas一Peukcer)算法
Douglas一Peukcer算法由D.Douglas和T.Peueker于1973年提出,简称D一P算法,是目前公认的线状要素化简经典算法.现有的线化简算法中,有相当一部分都是在该算法基础上进行改进 ...
道格拉斯-普克算法(Douglas–Peucker algorithm)
道格拉斯-普克算法(Douglas–Peucker algorithm,亦称为拉默-道格拉斯-普克算法.迭代适应点算法.分裂与合并算法)是将曲线近似表示为一系列点,并减少点的数量的一种算法.该算法的原 ...
【Python】道格拉斯-普克抽稀算法
常用的地图点压缩 1.算法应用道格拉斯-普克抽稀算法,是用来对大量冗余的图形数据点进行压缩以提取必要的数据点. 2.算法步骤对每一条曲线的首末点虚连一条直线,求所有点与直线的距离, 并找出最大距离 ...
opencv 凸包convexHull、道格拉斯-普克算法Douglas-Peucker algorithm、approxPloyDP 函数
凸包convexHull.道格拉斯-普克算法Douglas-Peucker algorithm.approxPloyDP 函数道格拉斯-普克算法(Douglas–Peucker algorithm) ...
OpenCV 学习笔记03 凸包convexHull、道格拉斯-普克算法Douglas-Peucker algorithm、approxPloyDP 函数...
凸形状内部的任意两点的连线都应该在形状里面. 1 道格拉斯-普克算法 Douglas-Peucker algorithm 这个算法在其他文章中讲述的非常详细,此处就详细撰述. 下图是引用维基百科的.ε ...

道格拉斯普克(Douglas-Peuker)算法python实现

前言

一、算法实现

二、输出结果

道格拉斯普克(Douglas-Peuker)算法python实现相关推荐

最新文章

热门文章