如何使用JTS或NTS查找曲线中超过指定角度的转向点？

嘿，这个需求其实挺常见的，咱们核心思路就是盯着曲线上连续的三个顶点——计算中间那个点和前后点形成的转向角度，只要这个角度超过你设定的阈值，那它就是你要找的转向点。下面分别给你Java（JTS）和C#（NTS）的实现方案：

JTS（Java）实现方法 #

步骤拆解 #

1. 从你的LineString（也就是那条曲线）里提取所有顶点坐标，用getCoordinates()就能拿到坐标数组。
2. 遍历数组里的内顶点——也就是从第2个点到倒数第2个点（因为得有前、中、后三个点才能算转向）。
3. 对每组连续的三个点p0（前）、p1（中）、p2（后），计算p0→p1和p1→p2这两个方向的夹角，取实际转向的角度大小（比如左转90度和右转90度都算90度）。
4. 把计算出的角度和你的阈值对比，超过的话就把p1加入转向点列表。

代码示例 #

import org.locationtech.jts.geom.Coordinate;
import org.locationtech.jts.geom.LineString;
import org.locationtech.jts.algorithm.Angle;
import java.util.ArrayList;
import java.util.List;

public class TurnPointDetector {
    // 阈值参数是弧度，如果你习惯用角度，记得用Math.toRadians(阈值角度)转换后传入
    public static List<Coordinate> findTurnPoints(LineString lineString, double thresholdRadians) {
        List<Coordinate> turnPoints = new ArrayList<>();
        Coordinate[] coords = lineString.getCoordinates();
        
        // 少于3个点的话根本不存在转向点，直接返回空列表
        if (coords.length < 3) {
            return turnPoints;
        }
        
        for (int i = 1; i < coords.length - 1; i++) {
            Coordinate p0 = coords[i-1];
            Coordinate p1 = coords[i];
            Coordinate p2 = coords[i+1];
            
            // 计算两个向量的方向角
            double dirFromP0ToP1 = Angle.angle(p0, p1);
            double dirFromP1ToP2 = Angle.angle(p1, p2);
            
            // 计算转向角度的绝对值，再取最小的那个（确保范围在0到π之间）
            double turnAngle = Math.abs(Angle.diff(dirFromP0ToP1, dirFromP1ToP2));
            turnAngle = Math.min(turnAngle, 2 * Math.PI - turnAngle);
            
            // 超过阈值就记录这个点
            if (turnAngle > thresholdRadians) {
                turnPoints.add(p1.copy()); // 复制坐标避免原对象被意外修改
            }
        }
        return turnPoints;
    }
}

NTS（C#）实现方法 #

NTS是JTS的C#移植版，逻辑和Java完全一致，只是API命名稍有不同：

代码示例 #

using NetTopologySuite.Geometries;
using NetTopologySuite.Algorithm;
using System.Collections.Generic;

public static class TurnPointDetector
{
    // 阈值参数是弧度，角度转弧度用 Math.PI * 角度值 / 180 即可
    public static List<Coordinate> FindTurnPoints(LineString lineString, double thresholdRadians)
    {
        List<Coordinate> turnPoints = new List<Coordinate>();
        Coordinate[] coords = lineString.Coordinates;
        
        if (coords.Length < 3)
            return turnPoints;
        
        for (int i = 1; i < coords.Length - 1; i++)
        {
            Coordinate p0 = coords[i-1];
            Coordinate p1 = coords[i];
            Coordinate p2 = coords[i+1];
            
            double dirFromP0ToP1 = AngleUtility.Angle(p0, p1);
            double dirFromP1ToP2 = AngleUtility.Angle(p1, p2);
            
            double turnAngle = System.Math.Abs(AngleUtility.Difference(dirFromP0ToP1, dirFromP1ToP2));
            turnAngle = System.Math.Min(turnAngle, 2 * System.Math.PI - turnAngle);
            
            if (turnAngle > thresholdRadians)
            {
                turnPoints.Add(new Coordinate(p1)); // 复制坐标
            }
        }
        return turnPoints;
    }
}

额外小贴士 #

• 如果你的曲线有很多密集的细碎点，建议先用抽稀算法（比如Douglas-Peucker）简化曲线，再查找转向点，避免得到一堆无意义的结果。JTS/NTS里都有现成的DouglasPeuckerSimplifier类可以用。
• 要是你需要区分左转还是右转，可以通过计算叉积的正负来判断——叉积为正就是左转，为负就是右转，这个可以在代码里额外加判断逻辑。

内容的提问来源于stack exchange，提问作者DuckQueen