經緯度計算距離

根據經緯度計算距離是有固定公式的,按照公式寫個函數很簡單。

import mathdef cal_dis(lat1, lon1,lat2, lon2):    latitude1 = (math.pi/180)*lat1    latitude2 = (math.pi/180)*lat2    longitude1 = (math.pi/180)*lon1    longitude2= (math.pi/180)*lon2    #因此AB兩點的球面距離為:{arccos[sinb*siny+cosb*cosy*cos(a-x)]}*R    #地球半徑    R = 6378.137    d = math.acos(math.sin(latitude1)*math.sin(latitude2)+ math.cos(latitude1)*math.cos(latitude2)*math.cos(longitude2-longitude1))*R    return dif __name__ == '__main__':    print cal_dis(23.0,101.1,23.06,113.34)

注意,經緯度是角度,而三角函數的輸入是弧度,角度與弧度的關系 180=pi

有一次我在一個項目中發現 matlab 的計算距離和上面這個函數有些出入,雖然相差不大,但是確實不一樣,當時也折騰了好久,后來發現 python 有個自帶的函數計算結果和 matlab 相同。

from geopy.distance import geodesicgeodesic((30.28708,120.12802999999997), (28.7427,115.86572000000001)).m