Pythonで緯度経度をインプットしたときに、複数の2地点間距離を求めるコードを作成しています。
公式をもとに距離を求める関数は作成できたのですが、今のコードだと各地点の緯度経度を都度インプットしてそれぞれ距離を出さないといけないので、地点数が増えたときに膨大な手間がかかります。
ですので添付画像のようなリストをインポートして、経度緯度情報をもとに総当たりで2地点間距離を出せるようにしたいのですが、どのような書き方をすれば良いかご教示お願いできませんでしょうか。
コードは以下
#モジュールimport
from math import sin, cos, acos, radians
#地球の半径を6371kmとする
earth_rad = 6371
##関数定義
def latlng_to_xyz(lat, lng):
rlat, rlng = radians(lat), radians(lng)
coslat = cos(rlat)
return coslat*cos(rlng), coslat*sin(rlng), sin(rlat)
def dist_on_sphere(pos0, pos1, radius=earth_rad):
xyz0, xyz1 = latlng_to_xyz(*pos0), latlng_to_xyz(*pos1)
return acos(sum(x * y for x, y in zip(xyz0, xyz1)))*radius
Osaka = 34.702113, 135.494807
Tokyo = 35.681541, 139.767103
London = 51.476853, 0.0
print(dist_on_sphere(Osaka, Tokyo)) # 403.63km
print(dist_on_sphere(London, Tokyo)) # 9571.22km