Leaflet中测地线与Turf多边形投影不一致致相交判断异常

我使用Leaflet和turf.js的turf.booleanIntersects()方法判断线段是否穿过多边形，代码可正常运行，但出现不一致情况：方法判定附件图中的线段穿过两个多边形，但地图上绘制的线段看起来并未穿过上方的多边形。想确认是地图显示有误（若有误该如何修复）、线段实际并未穿过多边形，还是测地线与GeoJSON多边形在Leaflet地图上的投影不一致导致的问题？

核心代码片段 #

LayerGroup = L.layerGroup().addTo(mymap);
var line = new L.Geodesic([[point1, point2]], geodesicOptions).addTo(LayerGroup);
LinesList.push(line);


var line = turf.lineString([
        [LinesList[i].points[0][0].lng, LinesList[i].points[0][0].lat],
        [LinesList[i].points[0][1].lng, LinesList[i].points[0][1].lat]
        ]);

// 多个点组成polygonlines数组
var point= [end.lng, end.lat];
polygonlines.push(point);

var geojsonPolygon =
    {
    "type": "Feature",
    "properties": {},
    "geometry": {
    "type": "Polygon",
    "coordinates": [polygonlines]
    }
}

var turfpolygon = turf.polygon(geojsonPolygon.geometry.coordinates);
if (turf.booleanIntersects(line, turfpolygon) == true) {
   const p2 = L.geoJSON(turfpolygon).addTo(mymap);
}                               

原因分析 #

1. 投影不一致是核心问题
   Turf.js的所有空间计算都是基于WGS84经纬度（EPSG:4326）的平面计算，完全不考虑地球曲率；但你用L.Geodesic绘制的是测地线（大圆航线），它是地球曲面上的最短路径，在Web墨卡托投影（EPSG:3857）的Leaflet地图上会显示为曲线。

   你传给Turf.js的线段是原始两点组成的平面直线（WGS84坐标系下的直线，对应地球曲面上的等角航线），和地图上显示的测地线并非同一条路径。Turf判断的是这条等角航线是否穿过多边形，而你视觉上看到的是测地线，两者路径差异导致了判断结果和显示不一致。

2. 地图显示本身无错误
   Leaflet默认采用Web墨卡托投影，L.Geodesic正确绘制了地球曲面上的最短路径，多边形则是按WGS84坐标直接投影到Web墨卡托上显示，显示逻辑本身没有问题。

修复方案 #

要让Turf的判断结果和地图显示的测地线一致，需要将测地线的离散坐标点传给Turf，而非原始的两点直线：

修改步骤 #

1. 获取L.Geodesic生成的测地线所有离散点
   L.Geodesic会将测地线拆分成多个线段模拟曲线，可通过getLatLngs()方法获取这些点：

   // 替换原有的turf.lineString创建代码
   const geodesicLatLngs = line.getLatLngs()[0]; // Geodesic的点为二维数组结构
   const turfLineCoords = geodesicLatLngs.map(latlng => [latlng.lng, latlng.lat]);
   var line = turf.lineString(turfLineCoords);
   

2. 使用该离散点组成的Turf线串执行相交判断
   这样Turf判断的就是地图上显示的测地线是否穿过多边形，结果会和视觉呈现一致。

额外注意事项 #

• 若多边形坐标基于Web墨卡托，需先转换为WGS84经纬度再传给Turf（Leaflet的L.geoJSON默认会处理坐标转换，经纬度创建的多边形无需额外操作）。
• 若测地线离散点数量过少可能存在误差，可调整L.Geodesic的steps参数增加点数，提升判断精度。

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