緯度経度のメートル換算

Question

緯度経度の変化を計算し、メートルに換算しそれをUnity内座標に反映
させたいです。

例　35.681302, 139.766965(開始地点の緯度経度) 35.681343, 139.766920(移動後の緯度経度)

Updateで現在地は更新し続けるのであくまでも一例です。

この2点の差を緯度経度別に計算します。
開始地点p1　移動後の点p2とする。

緯度　p2 - p1 ＝ 0.000041
経度 p2 - p1 = -0.000045

負の値も絶対値をとって正の値とする。移動方向はのちに考えることとする。
この得られた差をメートル換算したいわけです。

地球の演習を約40,000,000mとみなし360度で割ると1度あたりの変化量が求まる。
1度:111km
0.1度:11km
0.01:1111m
0.001:111m
0.0001:11m
0.00001:1m
0.000001:11cm
0.0000001:1cm

プログラミング言語はC#で実現したい換算方法を書きます。
緯度差0.000041の場合
整数部
0 × 111000(m)
小数部
0 × 11000(m)
0 × 1111(m)
・・・・

度数の変化量に応じた値を掛け算して最終的に整数部と小数部を一つの変数に格納したいのですが、プログラムがわからないので質問しました。
各桁に違う値を掛け算するので、値を取り出さなければなりません。

整数部と小数部にわけて計算したらよいのはわかりますが、小数部を切り出そうとして%1するとfloatやdoubleの性質上余計な値が混入します。

Answer 1

緯度1度あたりの距離というのは質問文にある通りですが、経度1度あたりというのは決められません。緯度によります。

ではどう計算するかですが、 https://qiita.com/yuba/items/4372944ce0f6a0bf6cb5 こちらにまとめました。

Comments

[KSwordOfHaste]

あ、そうですね！

Answer 2

整数部と小数部に分けて計算しても正確な値がでるわけではないと思います。例えば地球が完全な球と仮定していますし、円周も約40,000,000メートルと仮定しておられますので。計算の元とするデータに精度がないなら途中計算や計算結果に精度を求めても意味はないです。

ゆえに平易に全てdoubleで演算してよいと思います。

例えば移動前の緯度がla0, 異動後がla1だとすれば

(la1-la0)*40E6/360

でよろしいのではないでしょうか。

Comments

[KSwordOfHaste]

経度についてはyubaさんコメントのとおりです。赤道直下なら単純に上記でよいですが緯度が変化するに従い距離は短くなるので三角関数で計算しないといけませんね。自分のコメントは「整数部・小数部にわけて計算する」点についてのみの指摘と受け取ってくださいませ。