B
y
0
C
∣
A
2
B
2
d= \frac{|Ax_{0}+By_{0}+C|}{ \sqrt{A{2}+B{2}}}
d=A2+B2
∣Ax0+By0+C∣
2.2 源码
一下是优化前的代码,优化后的代码请查看《道格拉斯-普克 Douglas-Peucker 抽稀算法近1000倍的速度提升》两个对象封装类【为了简洁注解未使用DOC规范】:
// 类1 坐标数据封装
@Data
@NoArgsConstructor
@AllArgsConstructor
public class PointData {
// 标记是否删除(0保留 1删除)
private int isDelete;
// 数据点的经度
private double longitudeEx;
// 数据点的纬度
private double latitudeEx;
}
// 类2 点INDEX及到线距离封装
@Data
@NoArgsConstructor
@AllArgsConstructor
@Builder
public class DistanceData {
// 当前点的 INDEX 值
private int index;
// 当前点到直线的距离值
private double distance;
}
算法核心代码及测试代码【为了简洁部分注解未使用DOC规范】:
public class Douglas {
// 矢量曲线的点数据
private static List<PointData> points = new ArrayList<>();
// 距离阈值(值越大效果越好但丢失细节就越多)
private static final double DISTANCE_THRESHOLD = 0.00001;
/\*\*
\* 标记要删除的点
\*
\* @param from 矢量曲线的起始点
\* @param to 矢量曲线的终止点
\*/
private static void markDeletePoints(PointData from, PointData to) {
// 计算由起始点和终止点构成直线方程一般式的系数
double fromLat = from.getLatitudeEx();
double fromLng = from.getLongitudeEx();
double toLat = to.getLatitudeEx();
double toLng = to.getLongitudeEx();
// 求解点到直线距离方程的参数
double sqrtPowAddVal = Math.sqrt(Math.pow((fromLat - toLat), 2) + Math.pow((fromLng - toLng), 2));
double parameterA = (fromLat - toLat) / sqrtPowAddVal;
double parameterB = (toLng - fromLng) / sqrtPowAddVal;
double parameterC = (fromLng \* toLat - toLng \* fromLat) / sqrtPowAddVal;
int fromIndex = points.indexOf(from);
int toIndex = points.indexOf(to);
if (toIndex == fromIndex + 1) {
return;
}
// 起止点之间的点到起止点连线的距离数据集合
List<DistanceData> distanceDataList = new ArrayList<>();
double powAddRes = Math.pow(parameterA, 2) + Math.pow(parameterB, 2);
for (int i = fromIndex + 1; i < toIndex; i++) {
double lng = points.get(i).getLongitudeEx();
double lat = points.get(i).getLatitudeEx();
// 点到直线的距离
double distance = Math.abs(parameterA \* (lng) + para