1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
|
def countdistance(long1,la1,long2,la2): long1, la1, long2, la2 = map(radians, [long1, la1, long2, la2]) dlong = long2 - long1 dla = la2 - la1 a = sin(dla/2)**2 + cos(la1) * cos(la2) * sin(dlong/2)**2 c = 2 * asin(sqrt(a)) r = 6371.0088 lineardistance=c * r * 1000 return lineardistance def getLD(arrLike): long1 = arrLike['s_long'] la1 = arrLike['s_la'] long2 = arrLike['e_long'] la2 = arrLike['e_la'] LD = countdistance(long1,la1,long2,la2) return LD Thurdata['LD'] = Thurdata.apply(getLD,axis=1)
|