def cross(A,B): return A[0] * B[1] - A[1] * B[0] def vectorMinus( a , b): return ( (a[0] - b[0] )*1000,(a[1] - b[1] )*1000) def getLTDis( A, B ): lon1, lat1, lon2, lat2 = map(radians, [A[0], A[1], B[0], B[1]]) dlon = lon2 - lon1 dlat = lat2 - lat1 a…