如何用MapLibre GL判断两点间路径是否经过指定坐标？

核心思路 #

MapLibre GL 本身没有内置路线途经点检测的函数，我们可以通过几何计算实现需求：将生成的路线拆解为连续的线段，计算每个收费站坐标到这些线段的最短距离，当距离小于设定阈值时，判定路线经过该收费站。

实现步骤 #

1. 获取路线的坐标序列

从 MapLibre 路线服务返回的 GeoJSON 结果中，提取 geometry.coordinates（LineString 类型的坐标数组，格式为 [经度, 纬度]）。

2. 实现点到线段的距离计算

由于地理坐标是球面坐标，推荐使用球面距离公式计算点到线段的最短距离，也可以用轻量地理库简化开发。

方法1：纯 JS 实现球面距离计算

// 计算两点间的球面距离（单位：米）
function haversineDistance(coords1, coords2) {
  const R = 6371000; // 地球半径（米）
  const [lon1, lat1] = coords1;
  const [lon2, lat2] = coords2;
  
  const dLat = (lat2 - lat1) * Math.PI / 180;
  const dLon = (lon2 - lon1) * Math.PI / 180;
  
  const a = Math.sin(dLat/2) * Math.sin(dLat/2) +
            Math.cos(lat1 * Math.PI / 180) * Math.cos(lat2 * Math.PI / 180) *
            Math.sin(dLon/2) * Math.sin(dLon/2);
  const c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));
  
  return R * c;
}

// 计算点到线段的最短球面距离（单位：米）
function pointToSegmentDistance(point, segStart, segEnd) {
  const distToStart = haversineDistance(point, segStart);
  const distToEnd = haversineDistance(point, segEnd);
  
  // 线段长度
  const segLength = haversineDistance(segStart, segEnd);
  
  // 若点到两端点的距离均大于线段长度，取最小端点距离
  if (distToStart >= segLength && distToEnd >= segLength) {
    return Math.min(distToStart, distToEnd);
  }
  
  // 计算投影点到起点的比例（近似平面投影，小范围足够准确）
  const dot = ((point[0] - segStart[0]) * (segEnd[0] - segStart[0])) + 
              ((point[1] - segStart[1]) * (segEnd[1] - segStart[1]));
  const projectionRatio = Math.max(0, Math.min(1, dot / (segLength * segLength)));
  
  // 计算投影点坐标
  const projLon = segStart[0] + projectionRatio * (segEnd[0] - segStart[0]);
  const projLat = segStart[1] + projectionRatio * (segEnd[1] - segStart[1]);
  
  return haversineDistance(point, [projLon, projLat]);
}

方法2：使用 Turf.js 简化开发

如果允许引入轻量库，Turf.js 提供了直接的点到线距离计算函数：

import * as turf from '@turf/turf';

// 将路线转为 Turf LineString 对象
const routeLine = turf.lineString(routeCoordinates);

// 计算单个收费站到路线的距离（单位：米）
const distance = turf.pointToLineDistance(turf.point(tollCoords), routeLine, {units: 'meters'});

3. 批量检测收费站

遍历所有收费站坐标，逐个检查是否在路线的距离阈值内：

// 从MapLibre路由结果获取的路线坐标数组
const routeCoordinates = [...];
// 已获取的墨西哥收费站坐标数组（格式：[[lon1, lat1], [lon2, lat2], ...]）
const tollGates = [...];
// 距离阈值（可根据道路宽度调整，比如高速路设15米）
const distanceThreshold = 15;

// 筛选出路线经过的收费站
const passedTolls = tollGates.filter(toll => {
  for (let i = 0; i < routeCoordinates.length - 1; i++) {
    const segStart = routeCoordinates[i];
    const segEnd = routeCoordinates[i+1];
    const dist = pointToSegmentDistance(toll, segStart, segEnd);
    if (dist <= distanceThreshold) {
      return true;
    }
  }
  return false;
});

console.log('路线经过的收费站：', passedTolls);

注意事项 #

• 坐标格式统一：确保路线坐标和收费站坐标都是 [经度, 纬度]（MapLibre 标准格式）。
• 阈值灵活调整：根据道路类型设定阈值，普通道路可设5-10米，高速路设15-20米。
• 性能优化：如果路线极长或收费站数量多，可先通过边界框过滤（只检测路线 bounding box 内的收费站），减少计算量。

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