本文主要是介绍用python 求两地的直线距离，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！
通过百度，搜索到了球面直线距离计算公式
比如你要求point1(lng1,lat1) 和 point2(lng2,lat2)的距离Distance，那么
Distance = R*arccos(cos(lat1*pi/180 )*cos(lat2*pi/180)*cos(lng1*pi/180 -lng2*pi/180)+sin(lat1*pi/180 )*sin(lat2*pi/180))
其中，R=6370996.81 是地球半径，pi为圆周率π，(lng1,lat1),(lng2,lat2)分别是point1和point2的经纬度
R = 6371004 #The radius of the earth
PI = 3.1415926
MAX_NUM = 10000
def get_distance(points):length = len(points)if length <= 1:raise PathException('The length of points must be greater than or equal to 2')distance = [[0 for col in range(length)] for row in range(length)]for i in range(0,length-1):for j in range(i+1,length): f1 = points[i].find('(') + 1f2 = points[i].find(',')f3 = points[i].find(')')lnga = float(points[i][f1:f2]) * PI /180lngb = float(points[j][f1:f2]) * PI /180lata = float(points[i][f2+1:f3]) * PI /180latb = float(points[j][f2+1:f3]) * PI /180 C1 = math.cos(lata)*math.cos(latb)*math.cos(lnga-lngb)C2 = C1 + math.sin(lata)*math.sin(latb)D = R * math.acos(C2)distance[i][j] = round(D,2)distance[j][i] = distance[i][j] distance[i][i] = MAX_NUM distance[length-1][length-1] = MAX_NUM return distance
其中，输入的参数points，格式如下：
points = ['Point(123.4201360000,41.7713180000)',#致用路 1
'Point(123.4264960000,41.7752190000)',#北门 2
'Point(123.4262980000,41.7736180000)',#信息科学与工程学院 3
'Point(123.4262620000,41.7722200000)',#恩承图书馆 4
'Point(123.4279150000,41.7706590000)',#材料与冶金学院 5
'Point(123.4250770000,41.7693010000)',#求实路 6
'Point(123.4289570000,41.7684940000)',#知行楼 7
'Point(123.4234060000,41.7722200000)'#汉卿北路 左8
'Point(123.4292630000,41.7719100000)' #汉卿北路 右8
]
这篇关于用python 求两地的直线距离的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！
