如何沿两点最短距离方向切割环形LineString？

我有一个农田的Polygon，在其边界上标记两个Point以测量两点间距离。多数情况该方法运行正常，但有时切割出的LineString会走“远路”。

核心问题：如何沿两点最短距离方向切割LineString？

背景信息 #

• 先将Polygon转换为LineString
• 使用Turf函数在起止点切割线段

import { polygonToLine } from '@turf/polygon-to-line'
import lineSlice from '@turf/line-slice'

const slicedLineString = lineSlice(
    startPoint,
    endPoint,
    polygonToLine(polygonFeature)
)

示意图 #

• 农田空白示意图：
• 正确示例（黑线为最短距离）：
• 错误示例（黑线未沿最短距离）：

问题根源是lineSlice会严格遵循LineString的原始坐标顺序切割路径。当两点在Polygon边界上的位置导致沿长路径切割时，就会出现“远路”情况。解决思路是生成两条可能的路径（正向和反向），选择长度更短的那条。

实现代码 #

import { polygonToLine } from '@turf/polygon-to-line'
import lineSlice from '@turf/line-slice'
import lineSliceReverse from '@turf/line-slice-reverse'
import length from '@turf/length'

// 将多边形转为边界线
const polygonBoundary = polygonToLine(polygonFeature);
// 生成两条候选路径
const pathForward = lineSlice(startPoint, endPoint, polygonBoundary);
const pathReverse = lineSliceReverse(startPoint, endPoint, polygonBoundary);

// 选择更短的路径作为结果
const shortestPath = length(pathForward) < length(pathReverse) ? pathForward : pathReverse;

兼容处理 #

如果你的Turf版本没有lineSliceReverse，可以通过反转LineString坐标实现同样效果：

import { polygonToLine } from '@turf/polygon-to-line'
import lineSlice from '@turf/line-slice'
import length from '@turf/length'
import reverse from '@turf/reverse'

const polygonBoundary = polygonToLine(polygonFeature);
const pathForward = lineSlice(startPoint, endPoint, polygonBoundary);
// 反转边界线后切割得到反向路径
const pathReverse = lineSlice(startPoint, endPoint, reverse(polygonBoundary));

const shortestPath = length(pathForward) < length(pathReverse) ? pathForward : pathReverse;

注意事项 #

• 该方法适用于无孔洞的简单Polygon，如果是带孔洞的复杂农田边界，需要先筛选出外层边界再进行切割。
• 若两点距离极近（几乎重合），可以直接返回两点连线作为最短路径，避免不必要的计算。

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