Advertisement
Guest User

Untitled

a guest
Feb 18th, 2020
112
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.21 KB | None | 0 0
  1. import math
  2. import ast
  3.  
  4. def haversine2(l1, l2, l3, l4):
  5. R = 6372800 # Earth radius in meters
  6. lat1 = l1
  7. lon1 = l2
  8. lat2 = l3
  9. lon2 = l4
  10.  
  11. phi1, phi2 = math.radians(lat1), math.radians(lat2)
  12. dphi = math.radians(lat2 - lat1)
  13. dlambda = math.radians(lon2 - lon1)
  14.  
  15. a = math.sin(dphi/2)**2 + \
  16. math.cos(phi1)*math.cos(phi2)*math.sin(dlambda/2)**2
  17.  
  18. return 2*R*math.atan2(math.sqrt(a), math.sqrt(1 - a))
  19.  
  20. suma = 0
  21. i = 0
  22. f = open("data2.txt", "r")
  23. for line in f:
  24. c = line.split(',')
  25. l1 = float(c[0])
  26. l2 = float(c[1])
  27. a = f.readline()
  28. if len(a.strip()) != 0 :
  29. a = ast.literal_eval(a)
  30. #print(ast.literal_eval(a))
  31. #a = f.readline().split(',')
  32. if i == 1:
  33. # print the distance between our signals
  34. distance2 = haversine2(l1,l2,l3,l4)
  35. print("distance: ",distance2)
  36. l3 = l1
  37. l4 = l2
  38. print("array length: ",tam)
  39. d = float(distance2/tam)
  40. print("each call: ",d)
  41. for i in a:
  42. print(i)
  43. #m = float(i)
  44. #print("mmmmm", m)
  45. #if float(i) < 900:
  46. # sum = sum + d
  47. #else:
  48. # sum = 0
  49. print("suma: ",suma)
  50. tam = len(a)
  51. else:
  52. l3 = l1
  53. l4 = l2
  54. i = 1
  55. tam = len(a)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement