楕円体上の距離

id:yellow_73:20080502#p2のつづきというか自前で計算できないかと思って、楕円体上の距離についてもう少し考えてみる。
日本国内の任意の点について、その点から1キロ四方以内のポイントデータを抜き出せているとして、その点に最も近いポイントデータを得る、というもの。
前掲の地質調査所月報の付録2によると、Sodanoの式というのがあるそう。数式はそこに出ていました。
ただ、ITRF三次元直交座標系に変換して、それの直線距離を出しても、1キロ未満なら問題が無いのではないかと思い出した今日この頃。
まあ極端なこと言ったら、対象の中心について dX/dλ, dY/dφ を出して、これを係数として、線形化してしまってもとか。