关于空间分析JTS的总结——筛点

筛点:往往库里返回的数据的点非常密集,筛点后会甩掉大量的较小偏移的临近点,并根据点与点之前距离的设定(例如:距离50米以内的点甩掉)和最大甩点数(例如:最多甩到10000个点就不甩了),减小精度来增加空间运算速度。

实例代码:

distanceToleranceMeter:50 --最小甩点距离
maxCoordinateCount:10000  --最大甩点数

 

复制代码
@Override
    public <T> List<Coordinate> simplify(List<T> lngLatList, Function<T, Coordinate> transformer) {
        // 地图匹配最小距离 100m
        double distanceTolerance = GeometryUtil.getEarthArcDistance(distanceToleranceMeter);
        List<Coordinate> coordinates = Lists.newArrayList();

        Coordinate lastCoordinate = null;
        for (T lonLat : lngLatList) {
            // Longitude;Latitude;z,
            Coordinate coordinate = transformer.apply(lonLat);
            if (coordinate == null) {
                continue;
            }
            if (lastCoordinate == null || !coordinate.equals2D(lastCoordinate, 0.000001)) {
                coordinates.add(coordinate);
                lastCoordinate = coordinate;
            }
        }
        if (coordinates.isEmpty()) {
            return coordinates;
        }
        Coordinate[] pts = coordinates.toArray(new Coordinate[coordinates.size()]);
        coordinates.clear();
        double currentDistanceTolerance = distanceTolerance;
        do {
            // 精简点算法在参数【distanceTolerance】定值情况下,点不会再减少,故增加距离,避免死循环
            pts = DouglasPeuckerLineSimplifier.simplify(pts, currentDistanceTolerance);
            currentDistanceTolerance += GeometryUtil.getEarthArcDistance(10D);
        } while (pts.length > maxCoordinateCount);
        return Lists.newArrayList(pts);
    }
复制代码
复制代码
@Override
    public List<Coordinate> buildCoordinates(List<String> lngLatList, boolean inverseGeo) {
        List<Coordinate> coordinates  = simplify(lngLatList, new Function<String, Coordinate>() {
            public Coordinate apply(String coord) {
                try {
                    String[] arr = StringUtils.split(coord,";");
                    String longitude = arr[0];
                    String latitude = arr[1];
                    if (longitude != null && latitude != null) {
                        return new Coordinate(Double.parseDouble(longitude), Double.parseDouble(latitude));
                    }
                } catch (NumberFormatException e) {
                    logger.error(e.getMessage(), e);
                }
                return null;
            }
        });
        if(inverseGeo){//是否批量转偏移
            List<String[]> stringList = coordinates.stream().map(t -> new String[]{t.getX() + "", t.getY() + ""}).collect(Collectors.toList());
            List<String> result = mapService.inverseGeoConvert(stringList);
            coordinates = result.stream().map(t -> new Coordinate(Double.parseDouble(t.split(",")[0]), Double.parseDouble(t.split(",")[1]))).collect(Collectors.toList());
        }
        return coordinates;
    }
复制代码

 

posted @   蔡徐坤1987  阅读(709)  评论(0编辑  收藏  举报
编辑推荐:
· AI与.NET技术实操系列(二):开始使用ML.NET
· 记一次.NET内存居高不下排查解决与启示
· 探究高空视频全景AR技术的实现原理
· 理解Rust引用及其生命周期标识(上)
· 浏览器原生「磁吸」效果!Anchor Positioning 锚点定位神器解析
阅读排行:
· DeepSeek 开源周回顾「GitHub 热点速览」
· 物流快递公司核心技术能力-地址解析分单基础技术分享
· .NET 10首个预览版发布:重大改进与新特性概览!
· AI与.NET技术实操系列(二):开始使用ML.NET
· 单线程的Redis速度为什么快?
点击右上角即可分享
微信分享提示