こんにちは。
回転楕円体面上の二点間距離を求める Bowring (1996) の計算方法は全球で高精度ですが1、
- Total Inverse Solutions for the Geodesic and Great Elliptic, B. R. Bowring, Survey Review, 33, 261 (July 1996), 461–476.
- 経緯度を用いた2地点間の測地線長、方位角を求める計算(国土地理院)
- 330k/geodistance-js/gsidistance.js (GitHub) ← 下記をマージしてもらいました
しかしながら、率直に実装し浮動小数点計算すると、対蹠点近傍で j1 > 1 が起き NaN が生じる場合があり、最小限改善しました。