投稿

続 二点間距離 Python

 今回調べるにあたってPythonで書きましたが どこやらで、日本語が使えると見て。。漢字を多用しているので。。やった事だけ。。 なにせ。   def ラジアン(緯度2,経度2):  def 印刷(内容1,内容2,旗): 長半径 * acos(sin(緯度1)*sin(緯度2)+cos(緯度1)*cos(緯度2)*cos(二点の経度差))/1000  sqrt((子午線曲率半径 * 二点の緯度差)**2 + (卯酉線曲半径 * 二点の経度差 * cos(二点の緯度平均))**2)/1000 こんな感じになってるもんで、邪道極めた-;)  i=5 while i < 50000000: #(m)     roop(i)     i *= 10 で、5000kmまで の処理 def roop(i):     azi=0     while azi < 360:   y2 = geod.Direct(la1,lo1,azi,i) #geographiclibを使用して距離と角度で緯度経度を取る  処理  azi += 10 #角度を10度足す です。geographiclib 偉い!ってことです

続 二点間距離

イメージ
 前日の物を、5000kmまで範囲を広げmy mapで表示してみた 4点の角度のスクリーンショット  Hubeny式では、確かに180度をまたぐとですね 結論は、自分に合った式を使えばいいですね^^; 1km未満ならコストの安い 球面三角、GPSで国内使用なら Hubeny式、技術者の方は、Karney式でしょうか?  Vincenty式も精度高いですが、対蹠地やその付近では処理回数が増えます

1000km地点まで

 どれを使うべきか?悩んだので調べてみた  34.68749750521538 135.52589615498832     角度 0 距離 0.001 km wikipe       0.0010013258907237497 Hubeny       0.0009999999997330844 Vincenty     0.0010000000012002728 Karney       0.0010000000008480198 34.687497368267906 135.52589805001     角度 10 距離 0.001 km wikipe       0.001005826254481053 Hubeny       0.0009999999985602706 Vincenty     0.0009999998439811167 Karney       0.0010000000001367118 34.68749696158658 135.5258998874524     角度 20 距離 0.001 km wikipe       0.0010013258907237497 Hubeny       0.0009999999988614258 Vincenty     0.0009999993962663063 Karney       0.0009999999990664122 34.68749629752823 135.5259016114857     角度 30 距離 0.001 km wikipe       0...