alt = line[25:30]
#print dt
#print "%s, %s, %s" % (lat, lon, alt)
pos = Position(lat, lon, alt)
flight.add_gps_position(dt, pos)
return flight
wereld = World()
for filename in sys.argv[1:]:
current_flight = read_igc_file(filename)
wereld.add_flight(current_flight)
begin_time = datetime.datetime(2009, 05, 30, 6)
end_time = datetime.datetime(2009, 05, 31, 22)
for time in daterange(begin_time, end_time, datetime.timedelta(seconds=1)):
distances = []
for [i,flight] in enumerate(wereld.get_flights()):
flight_cur_pos = flight.get_position(time)
if flight_cur_pos is not None:
current_distances = wereld.get_distance_to_others(flight_cur_pos, time, 40.0)
distances += current_distances
if len(distances) > 0:
[mean, stdev] = mean_stdev(distances)
print str(time.strftime("%s")) + "," + str(mean) + "," + str(stdev)
sys.stdout.flush()
