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()