Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import math
- import ast
- def haversine2(l1, l2, l3, l4):
- R = 6372800 # Earth radius in meters
- lat1 = l1
- lon1 = l2
- lat2 = l3
- lon2 = l4
- phi1, phi2 = math.radians(lat1), math.radians(lat2)
- dphi = math.radians(lat2 - lat1)
- dlambda = math.radians(lon2 - lon1)
- a = math.sin(dphi/2)**2 + \
- math.cos(phi1)*math.cos(phi2)*math.sin(dlambda/2)**2
- return 2*R*math.atan2(math.sqrt(a), math.sqrt(1 - a))
- suma = 0
- i = 0
- f = open("data2.txt", "r")
- for line in f:
- c = line.split(',')
- l1 = float(c[0])
- l2 = float(c[1])
- a = f.readline()
- if len(a.strip()) != 0 :
- a = ast.literal_eval(a)
- #print(ast.literal_eval(a))
- #a = f.readline().split(',')
- if i == 1:
- # print the distance between our signals
- distance2 = haversine2(l1,l2,l3,l4)
- print("distance: ",distance2)
- l3 = l1
- l4 = l2
- print("array length: ",tam)
- d = float(distance2/tam)
- print("each call: ",d)
- for i in a:
- print(i)
- #m = float(i)
- #print("mmmmm", m)
- #if float(i) < 900:
- # sum = sum + d
- #else:
- # sum = 0
- print("suma: ",suma)
- tam = len(a)
- else:
- l3 = l1
- l4 = l2
- i = 1
- tam = len(a)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement