def main(): #Load all files and initilize Simple Tracks os.chdir("/home/moyano/Projects/CreateTracks/trips/") all_points = [] for trip in os.listdir("."): trip_data = load_file(trip) all_points += trip_data #tracks.append(Track(trip_data)) boundries = max_bounding_rect(all_points) # depth = depth(boundries, 0.00035) #print "Nesting Level: %i" % depth qtree = QuadTree(10, boundries) #Make the QTree for coord in all_points: qtree.add_point(coord) qtree.traverse() nodes = qtree.leaves #Load Trips trips = [] for trip in os.listdir("."): if not trip.startswith('.'): gps_data = load_file(trip) trips.append(Trip(gps_data,trip)) routes(trips, qtree) #Weighted Points os.chdir("/home/moyano/Projects/CreateTracks/edges") test_file = open("edges.txt", "w") test_file.write("latitude, longitude, ocurrences, color") # print children for node in nodes: p = node._center_of_mass() count = len(node.items) if count > 2: test_file.write("\n") test_file.write(str(p.latitude) + "," + str(p.longitude) + "," + str(count) + "," + get_color(count)) #All points os.chdir("/home/moyano/Projects/CreateTracks/edges") test_file = open("edges2.csv", "w") test_file.write("latitude, longitude") test_file.write("#") # print children xpoints = [] for node in nodes: p = node._center_of_mass() count = len(node.items) if count > 100: for _ in xrange(count): xpoints.append((p.latitude, p.longitude)) test_file.write(str(p.latitude) + ", " + str(p.longitude)) test_file.write("#")