Turf.js中rhumbBearing方法在特定坐标下返回NaN的问题分析

问题背景

在使用Turf.js地理空间分析库时，开发者发现其rhumbBearing方法在某些特定坐标输入下会返回NaN值。这个问题最初出现在使用transformRotate方法旋转多边形时被发现的，因为transformRotate内部调用了rhumbBearing方法。

问题复现

通过以下简单代码可以复现该问题：

const pivot = [0, 91];  // 经度0，纬度91
const pointCoords = [0, 0];  // 经度0，纬度0
console.log(turf.rhumbBearing(pivot, pointCoords)); // 输出NaN

技术分析

rhumbBearing方法用于计算两点之间的恒向线方位角（也称为罗盘方位角）。该方法在内部实现时，会进行一些三角函数和对数运算来计算方位角。

当输入坐标超出常规地理坐标范围时（经度-180到180，纬度-90到90），会导致计算过程中出现无效的数学运算。例如：

1. 纬度值超过90度会导致某些三角函数计算失效
2. 在对数运算中，当参数为负数时也会导致NaN结果

根本原因

Turf.js的设计初衷是处理标准地理坐标系统（WGS84）下的空间数据，其所有方法都假设输入坐标在以下范围内：

• 经度：-180到180度
• 纬度：-90到90度

当输入坐标超出这些范围时，虽然Turf.js没有直接报错，但会导致内部数学计算失效，从而返回NaN。

解决方案建议

1. 数据预处理：在使用Turf.js方法前，确保所有坐标点都在有效范围内
2. 坐标规范化：对超出范围的坐标进行规范化处理（如经度超过180时减去360）
3. 替代方案选择：如果是纯二维平面几何计算，建议使用专门的2D几何库而非地理空间库

最佳实践

对于需要在二维平面上进行几何变换（如旋转）的应用场景，建议：

1. 评估是否真的需要使用地理空间库
2. 考虑使用专门的2D几何处理库，如flatten-js等
3. 如果必须使用Turf.js，确保所有坐标数据都经过规范化处理

总结

Turf.js作为地理空间分析库，其方法对输入坐标有明确的取值范围要求。开发者在使用时需要注意数据预处理，确保坐标值在有效范围内。对于纯二维几何计算场景，选择更适合的工具库可能是更好的解决方案。