続 二点間距離 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 偉い!ってことです