地理空间-Java实现航迹稀释

Java实现航迹点稀释算法(Douglas - Peucker算法)的示例代码,该算法可在保证航迹整体形状变化不大的情况下减少航迹点数量:


import java.util.ArrayList;

import java.util.List;

class Point {

double x;

double y;

public Point(double x, double y) {

this.x = x;

this.y = y;

}

}

public class TrackThinning {

public static List<Point> douglasPeucker(List<Point> points, double epsilon) {

if (points.size() < 3) {

return points;

}

int index = -1;

double dmax = 0;

int end = points.size();

for (int i = 1; i < end - 1; i++) {

double d = perpendicularDistance(points.get(i), points.get(0), points.get(end - 1));

if (d > dmax) {

index = i;

dmax = d;

}

}

List<Point> result = new ArrayList<>();

if (dmax > epsilon) {

List<Point> recursiveResult1 = douglasPeucker(points.subList(0, index + 1), epsilon);

List<Point> recursiveResult2 = douglasPeucker(points.subList(index, points.size()), epsilon);

result.addAll(recursiveResult1.subList(0, recursiveResult1.size() - 1));

result.addAll(recursiveResult2);

} else {

result.add(points.get(0));

result.add(points.get(points.size() - 1));

}

return result;

}

private static double perpendicularDistance(Point point, Point start, Point end) {

double dx = end.x - start.x;

double dy = end.y - start.y;

double numerator = Math.abs((dy * point.x - dx * point.y) + (end.x * start.y - end.y * start.x));

double denominator = Math.sqrt(dy * dy + dx * dx);

return numerator / denominator;

}

public static void main(String\[\] args) {

List<Point> points = new ArrayList<>();

points.add(new Point(0, 0));

points.add(new Point(1, 1));

points.add(new Point(2, 2));

points.add(new Point(3, 2));

points.add(new Point(4, 3));

points.add(new Point(5, 4));

double epsilon = 0.5;

List<Point> thinnedPoints = douglasPeucker(points, epsilon);

for (Point p : thinnedPoints) {

System.out.println("(" + p.x + ", " + p.y + ")");

}

}

}

相关推荐
晓晓_za89866819 小时前
开源 GEO 优化源码二次开发:贴牌改造与业务模块扩展实践
java·运维·服务器·开发语言·性能优化·开源
别慌,让我先缓缓19 小时前
Java入门:从编写到运行第一步
java
liangshanbo121519 小时前
Nginx 进程管理命令深度面试题解析
java·nginx·面试
Mr. zhihao19 小时前
死锁排查实战:JVM 唯一会“自动报案“的问题(场景 B5)
java·jvm·gc
LINgZone219 小时前
后端Java编码规范
java·开发语言
weixin_3993806920 小时前
TongWeb7.0.4.9M11一键安装脚本(by sy)
java·spring
CallFay云起未来20 小时前
AI客服能不能减少人工回复?从重复咨询到人机协同的落地分析
java·大数据·人工智能·文心一言
我命由我1234520 小时前
Android 开发 - 获取当前设备屏幕的旋转角度
android·java·java-ee·android studio·android jetpack·android-studio·android runtime
Elastic 中国社区官方博客20 小时前
OpenTelemetry Java 扩展:无需分叉 agent 即可自定义追踪
java·大数据·运维·开发语言·数据库·人工智能·elasticsearch
京东云开发者21 小时前
上游给空、下游拿到 -1,我把故障一直追到了 commons-beanutils 的构造函数
java·ai编程