def trajectory_from_dictionary(dictionary): """Returns a trajectory constructed from the given dictionary. Args: dictionary: the dictionary to convert into a trajectory """ #verify domain is valid and import appropriate domain try: domain = importlib.import_module("tracktable.domain."+ dictionary['domain'].lower()) except ImportError: raise ValueError("Error: invalid domain name: "+dictionary['domain'].lower()) dimension = domain.DIMENSION #verify each coordinate matches dimension points = [] numPoints = len(dictionary['coordinates']) for point in dictionary['coordinates']: if len(point) != dimension: raise ValueError("Error: point " + str(point) + " has "+ str(len(point)) + " coordinate(s), expected " + str(dimension)+".") #verify point properties values lists are of equal length for (name, attributes) in dictionary['point_properties'].items(): if len(attributes['values']) != numPoints: raise ValueError("Error: "+name+" property has only " + str(len(attributes['values'])) + " values, but there are " + str(numPoints) + " points in the trajectory.") #verify there are the right number of timestamps if len(dictionary['timestamps']) != numPoints: raise ValueError("Error: Only " + str(len(dictionary['timestamps'])) + " timestamp values, but there are " + str(numPoints) + " points in the trajectory.") #verify object_id is a string if not isinstance(dictionary['object_id'], str): raise ValueError("Error: object_id must be a string, but got a value of type "+type(dictionary['object_id']).__name__) #generate points / position list for i in range(numPoints): point = domain.TrajectoryPoint(dictionary['coordinates'][i]) point.object_id = dictionary['object_id'] point.timestamp = Timestamp.from_string(dictionary['timestamps'][i]) for (name, attributes) in dictionary['point_properties'].items(): if attributes['type'] == "datetime" or attributes['type'] == "timestamp": #okay to support both? todo point.set_property(name, Timestamp.from_string(attributes['values'][i], format_string='%Y-%m-%d %H:%M:%S')) else: point.set_property(name, attributes['values'][i]) points.append(point) #make trajectory trajectory = domain.Trajectory.from_position_list(points) #add trajectory properties for (name, attributes) in dictionary['trajectory_properties'].items(): if attributes['type'] == "datetime" or attributes['type'] == "timestamp": #okay to support both? todo trajectory.set_property(name, Timestamp.from_string(attributes['value'], format_string='%Y-%m-%d %H:%M:%S')) else: trajectory.set_property(name, attributes['value']) return trajectory
def test_no_time_zone_iso(self): timestring = '2014-12-25T13:12:45' parsed_timestamp = Timestamp.from_string(timestring) constructed_timestamp = datetime(2014, 12, 25, hour=13, minute=12, second=45, tzinfo=pytz.utc) self.assertTrue(parsed_timestamp.tzinfo is not None) self.assertTrue(parsed_timestamp == constructed_timestamp)
def test_with_time_zone_iso(self): timestring = '2014-12-25T10:12:45-0300' parsed_timestamp = Timestamp.from_string( timestring, format_string='%Y-%m-%dT%H:%M:%S%z') constructed_timestamp = datetime(2014, 12, 25, hour=13, minute=12, second=45, tzinfo=pytz.utc) self.assertTrue(parsed_timestamp.tzinfo is not None) self.assertTrue(parsed_timestamp == constructed_timestamp)
def gen_json_and_trajectory(self): domain = tracktable.domain.terrestrial json = "{\"coordinates\": [[26.995, -81.9731], [27.0447, -81.9844], [27.1136, -82.0458]], \"domain\": \"terrestrial\", \"object_id\": \"AAA001\", \"point_properties\": {\"altitude\": {\"type\": \"int\", \"values\": [2700, 4200, 6700]}, \"heading\": {\"type\": \"float\", \"values\": [108.1, 108.2, 225.3]}, \"note\": {\"type\": \"str\", \"values\": [\"hello\", \"world\", \"!\"]}, \"time2\": {\"type\": \"datetime\", \"values\": [\"2004-01-01 00:00:01\", \"2004-01-01 00:00:02\", \"2004-01-01 00:00:03\"]}}, \"timestamps\": [\"2004-12-07 11:36:18\", \"2004-12-07 11:37:56\", \"2004-12-07 11:39:18\"], \"trajectory_properties\": {\"percent\": {\"type\": \"float\", \"value\": 33.333}, \"platform\": {\"type\": \"str\", \"value\": \"Boeing 747\"}, \"start\": {\"type\": \"datetime\", \"value\": \"2004-12-07 11:36:00\"}, \"tailNum\": {\"type\": \"int\", \"value\": 3878}}}" # Manually set up a matching Trajectory object p1 = domain.TrajectoryPoint() p1[0] = 26.995 if domain.DIMENSION > 1: p1[1] = -81.9731 if domain.DIMENSION > 2: p1[2] = -100.0 p1.set_property('altitude', 2700) p1.set_property('heading', 108.1) p1.set_property('note', "hello") p1.set_property('time2', Timestamp.from_string('2004-01-01 00:00:01')) p1.object_id = 'AAA001' p1.timestamp = Timestamp.from_string('2004-12-07 11:36:18', format_string='%Y-%m-%d %H:%M:%S') p2 = domain.TrajectoryPoint() p2[0] = 27.0447 if domain.DIMENSION > 1: p2[1] = -81.9844 if domain.DIMENSION > 2: p2[2] = -101.0 p2.set_property('altitude', 4200) p2.set_property('heading', 108.2) p2.set_property('note', "world") p2.set_property('time2', Timestamp.from_string('2004-01-01 00:00:02')) p2.object_id = 'AAA001' p2.timestamp = Timestamp.from_string('2004-12-07 11:37:56', format_string='%Y-%m-%d %H:%M:%S') p3 = domain.TrajectoryPoint() p3[0] = 27.1136 if domain.DIMENSION > 1: p3[1] = -82.0458 if domain.DIMENSION > 2: p3[2] = -102.0 p3.set_property('altitude', 6700) p3.set_property('heading', 225.3) p3.set_property('note', "!") p3.set_property('time2', Timestamp.from_string('2004-01-01 00:00:03')) p3.object_id = 'AAA001' p3.timestamp = Timestamp.from_string('2004-12-07 11:39:18', format_string='%Y-%m-%d %H:%M:%S') trajectory = domain.Trajectory.from_position_list([p1, p2, p3]) trajectory.set_property('percent', 33.333) trajectory.set_property('platform', "Boeing 747") trajectory.set_property('start', Timestamp.from_string('2004-12-07 11:36:00')) trajectory.set_property('tailNum', 3878) return json, trajectory
def test_no_time_zone_iso(self): timestring = '2014-12-25T13:12:45' parsed_timestamp = Timestamp.from_string(timestring) constructed_timestamp = datetime(2014, 12, 25, hour=13, minute=12, second=45, tzinfo=pytz.utc) self.assertTrue(parsed_timestamp.tzinfo is not None) self.assertTrue(parsed_timestamp == constructed_timestamp)
def test_with_time_zone_iso(self): timestring = '2014-12-25T10:12:45-0300' parsed_timestamp = Timestamp.from_string( timestring, format_string='%Y-%m-%dT%H:%M:%S%z') constructed_timestamp = datetime(2014, 12, 25, hour=13, minute=12, second=45, tzinfo=pytz.utc) self.assertTrue(parsed_timestamp.tzinfo is not None) self.assertTrue(parsed_timestamp == constructed_timestamp)
def gen_json_and_trajectory(self): domain = tracktable.domain.terrestrial json = "{\"coordinates\": [[26.995, -81.9731], [27.0447, -81.9844], [27.1136, -82.0458]], \"domain\": \"terrestrial\", \"object_id\": \"AAA001\", \"point_properties\": {\"altitude\": {\"type\": \"int\", \"values\": [2700, 4200, 6700]}, \"heading\": {\"type\": \"float\", \"values\": [108.1, 108.2, 225.3]}, \"note\": {\"type\": \"str\", \"values\": [\"hello\", \"world\", \"!\"]}, \"time2\": {\"type\": \"datetime\", \"values\": [\"2004-01-01 00:00:01\", \"2004-01-01 00:00:02\", \"2004-01-01 00:00:03\"]}}, \"timestamps\": [\"2004-12-07 11:36:18\", \"2004-12-07 11:37:56\", \"2004-12-07 11:39:18\"], \"trajectory_properties\": {\"percent\": {\"type\": \"float\", \"value\": 33.333}, \"platform\": {\"type\": \"str\", \"value\": \"Boeing 747\"}, \"start\": {\"type\": \"datetime\", \"value\": \"2004-12-07 11:36:00\"}, \"tailNum\": {\"type\": \"int\", \"value\": 3878}}}" # Manually set up a matching Trajectory object p1 = domain.TrajectoryPoint() p1[0] = 26.995; if domain.DIMENSION > 1: p1[1] = -81.9731 if domain.DIMENSION > 2: p1[2] = -100.0 p1.set_property('altitude', 2700) p1.set_property('heading', 108.1) p1.set_property('note', "hello") p1.set_property('time2', Timestamp.from_string('2004-01-01 00:00:01')) p1.object_id = 'AAA001' p1.timestamp = Timestamp.from_string('2004-12-07 11:36:18', format_string='%Y-%m-%d %H:%M:%S') p2 = domain.TrajectoryPoint() p2[0] = 27.0447; if domain.DIMENSION > 1: p2[1] = -81.9844 if domain.DIMENSION > 2: p2[2] = -101.0 p2.set_property('altitude', 4200) p2.set_property('heading', 108.2) p2.set_property('note', "world") p2.set_property('time2', Timestamp.from_string('2004-01-01 00:00:02')) p2.object_id = 'AAA001' p2.timestamp = Timestamp.from_string('2004-12-07 11:37:56', format_string='%Y-%m-%d %H:%M:%S') p3 = domain.TrajectoryPoint() p3[0] = 27.1136; if domain.DIMENSION > 1: p3[1] = -82.0458 if domain.DIMENSION > 2: p3[2] = -102.0 p3.set_property('altitude', 6700) p3.set_property('heading', 225.3) p3.set_property('note', "!") p3.set_property('time2', Timestamp.from_string('2004-01-01 00:00:03')) p3.object_id = 'AAA001' p3.timestamp = Timestamp.from_string('2004-12-07 11:39:18', format_string='%Y-%m-%d %H:%M:%S') trajectory = domain.Trajectory.from_position_list([p1,p2,p3]) trajectory.set_property('percent', 33.333) trajectory.set_property('platform', "Boeing 747") trajectory.set_property('start', Timestamp.from_string('2004-12-07 11:36:00')) trajectory.set_property('tailNum', 3878) return json, trajectory
def test_with_time_zone_iso(self): if sys.version_info[0] == 3: timestring = '2014-12-25T10:12:45-0300' parsed_timestamp = Timestamp.from_string( timestring, format_string='%Y-%m-%dT%H:%M:%S%z') constructed_timestamp = datetime(2014, 12, 25, hour=13, minute=12, second=45, tzinfo=pytz.utc) self.assertTrue(parsed_timestamp.tzinfo is not None) self.assertTrue(parsed_timestamp == constructed_timestamp) else: print("Time zone parsing is only supported in Python 3") self.assertTrue(True)
def gen_dictionary_and_trajectory(self, domainString): domain = importlib.import_module("tracktable.domain." + domainString.lower()) # Manually set up a dictionary which contains trajectory info dictionary = { 'domain': domainString.lower(), 'object_id': 'AAA001', 'trajectory_properties': { 'percent': { 'type': 'float', 'value': 33.333 }, 'platform': { 'type': 'str', 'value': "Boeing 747" }, 'start': { 'type': 'datetime', 'value': "2004-12-07 11:36:00" }, 'tailNum': { 'type': 'str', 'value': '3878' } }, 'timestamps': [ '2004-12-07 11:36:18', '2004-12-07 11:37:56', '2004-12-07 11:39:18' ], 'point_properties': { 'altitude': { 'type': 'float', 'values': [2700, 4200, 6700] }, 'heading': { 'type': 'float', 'values': [108.1, 108.2, 225.3] }, 'note': { 'type': 'str', 'values': ["hello", "world", "!"] }, 'time2': { 'type': 'datetime', 'values': [ "2004-01-01 00:00:01", "2004-01-01 00:00:02", "2004-01-01 00:00:03" ] } } } if domain.DIMENSION == 1: dictionary.update( {'coordinates': [(26.995), (27.0447), (27.1136)]}) if domain.DIMENSION == 2: dictionary.update({ 'coordinates': [(26.995, -81.9731), (27.0447, -81.9844), (27.1136, -82.0458)] }) if domain.DIMENSION == 3: dictionary.update({ 'coordinates': [(26.995, -81.9731, -100.0), (27.0447, -81.9844, -101.0), (27.1136, -82.0458, -102.0)] }) # Manually set up a matching Trajectory object p1 = domain.TrajectoryPoint() p1[0] = 26.995 if domain.DIMENSION > 1: p1[1] = -81.9731 if domain.DIMENSION > 2: p1[2] = -100.0 p1.set_property('altitude', 2700) p1.set_property('heading', 108.1) p1.set_property('note', "hello") p1.set_property('time2', Timestamp.from_string('2004-01-01 00:00:01')) p1.object_id = 'AAA001' p1.timestamp = Timestamp.from_string('2004-12-07 11:36:18', format_string='%Y-%m-%d %H:%M:%S') p2 = domain.TrajectoryPoint() p2[0] = 27.0447 if domain.DIMENSION > 1: p2[1] = -81.9844 if domain.DIMENSION > 2: p2[2] = -101.0 p2.set_property('altitude', 4200) p2.set_property('heading', 108.2) p2.set_property('note', "world") p2.set_property('time2', Timestamp.from_string('2004-01-01 00:00:02')) p2.object_id = 'AAA001' p2.timestamp = Timestamp.from_string('2004-12-07 11:37:56', format_string='%Y-%m-%d %H:%M:%S') p3 = domain.TrajectoryPoint() p3[0] = 27.1136 if domain.DIMENSION > 1: p3[1] = -82.0458 if domain.DIMENSION > 2: p3[2] = -102.0 p3.set_property('altitude', 6700) p3.set_property('heading', 225.3) p3.set_property('note', "!") p3.set_property('time2', Timestamp.from_string('2004-01-01 00:00:03')) p3.object_id = 'AAA001' p3.timestamp = Timestamp.from_string('2004-12-07 11:39:18', format_string='%Y-%m-%d %H:%M:%S') trajectory = domain.Trajectory.from_position_list([p1, p2, p3]) trajectory.set_property('percent', 33.333) trajectory.set_property('platform', "Boeing 747") trajectory.set_property('start', Timestamp.from_string('2004-12-07 11:36:00')) trajectory.set_property('tailNum', '3878') return dictionary, trajectory
def gen_dictionary_and_trajectory(self, domainString): domain = importlib.import_module("tracktable.domain."+domainString.lower()) # Manually set up a dictionary which contains trajectory info dictionary = {'domain' : domainString.lower(), 'object_id' : 'AAA001', 'trajectory_properties' : {'percent' : {'type': 'float', 'value': 33.333}, 'platform': {'type': 'str', 'value': "Boeing 747"}, 'start' : {'type': 'datetime', 'value': "2004-12-07 11:36:00"}, 'tailNum' : {'type': 'int', 'value': 3878}}, 'timestamps' : ['2004-12-07 11:36:18',#-0000', '2004-12-07 11:37:56',#-0000', '2004-12-07 11:39:18'],#-0000'], 'point_properties' : {'altitude': {'type': 'int', 'values': [2700, 4200, 6700]}, 'heading' : {'type': 'float', 'values': [108.1, 108.2, 225.3]}, 'note' : {'type': 'str', 'values': ["hello", "world", "!"]}, 'time2' : {'type': 'datetime', 'values': ["2004-01-01 00:00:01", "2004-01-01 00:00:02", "2004-01-01 00:00:03"]} } } if domain.DIMENSION == 1: dictionary.update({'coordinates' : [(26.995), (27.0447), (27.1136)]}) if domain.DIMENSION == 2: dictionary.update({'coordinates' : [(26.995, -81.9731), (27.0447, -81.9844), (27.1136, -82.0458)]}) if domain.DIMENSION == 3: dictionary.update({'coordinates' : [(26.995, -81.9731, -100.0), (27.0447, -81.9844, -101.0), (27.1136, -82.0458, -102.0)]}) # Manually set up a matching Trajectory object p1 = domain.TrajectoryPoint() p1[0] = 26.995; if domain.DIMENSION > 1: p1[1] = -81.9731 if domain.DIMENSION > 2: p1[2] = -100.0 p1.set_property('altitude', 2700) p1.set_property('heading', 108.1) p1.set_property('note', "hello") p1.set_property('time2', Timestamp.from_string('2004-01-01 00:00:01')) p1.object_id = 'AAA001' p1.timestamp = Timestamp.from_string('2004-12-07 11:36:18', format_string='%Y-%m-%d %H:%M:%S') p2 = domain.TrajectoryPoint() p2[0] = 27.0447; if domain.DIMENSION > 1: p2[1] = -81.9844 if domain.DIMENSION > 2: p2[2] = -101.0 p2.set_property('altitude', 4200) p2.set_property('heading', 108.2) p2.set_property('note', "world") p2.set_property('time2', Timestamp.from_string('2004-01-01 00:00:02')) p2.object_id = 'AAA001' p2.timestamp = Timestamp.from_string('2004-12-07 11:37:56', format_string='%Y-%m-%d %H:%M:%S') p3 = domain.TrajectoryPoint() p3[0] = 27.1136; if domain.DIMENSION > 1: p3[1] = -82.0458 if domain.DIMENSION > 2: p3[2] = -102.0 p3.set_property('altitude', 6700) p3.set_property('heading', 225.3) p3.set_property('note', "!") p3.set_property('time2', Timestamp.from_string('2004-01-01 00:00:03')) p3.object_id = 'AAA001' p3.timestamp = Timestamp.from_string('2004-12-07 11:39:18', format_string='%Y-%m-%d %H:%M:%S') trajectory = domain.Trajectory.from_position_list([p1,p2,p3]) trajectory.set_property('percent', 33.333) trajectory.set_property('platform', "Boeing 747") trajectory.set_property('start', Timestamp.from_string('2004-12-07 11:36:00')) trajectory.set_property('tailNum', 3878) return dictionary, trajectory
def run_test(): # Define three sets of points: ABQ to San Diego, San Diego to # Seattle, Denver to NYC. ABQ->SAN->SEA should break into two # trajectories because of a timestamp break in San Diego. The # flight to Denver will begin right when the flight to Seattle # ends so we expect to break that one based on the distance # threshold. print("Beginning run_test()") from tracktable.domain.terrestrial import TrajectoryPoint albuquerque = TrajectoryPoint( -106.5, 35.25 ) albuquerque.timestamp = Timestamp.from_string('2010-01-01 12:00:00') albuquerque.object_id = 'flight1' san_diego1 = TrajectoryPoint( -117.16, 32.67 ) san_diego1.timestamp = Timestamp.from_string('2010-01-01 15:00:00') san_diego1.object_id = 'flight1' san_diego2 = TrajectoryPoint( -117.16, 32.67 ) san_diego2.timestamp = Timestamp.from_string('2010-01-01 16:00:00') san_diego2.object_id = 'flight1' seattle = TrajectoryPoint( -122.31, 47.60 ) seattle.timestamp = Timestamp.from_string('2010-01-01 19:00:00') seattle.object_id = 'flight1' denver = TrajectoryPoint( -104.98, 39.79 ) denver.timestamp = Timestamp.from_string('2010-01-01 19:01:00') denver.object_id = 'flight1' new_york = TrajectoryPoint( -74.02, 40.71 ) new_york.timestamp = Timestamp.from_string('2010-01-02 00:00:00') new_york.object_id = 'flight1' # Now we want sequences of points for each flight. abq_to_sd = TrajectoryPointSource() abq_to_sd.start_point = albuquerque abq_to_sd.end_point = san_diego1 abq_to_sd.num_points = 180 sd_to_sea = TrajectoryPointSource() sd_to_sea.start_point = san_diego2 sd_to_sea.end_point = seattle sd_to_sea.num_points = 360 # flying very slowly denver_to_nyc = TrajectoryPointSource() denver_to_nyc.start_point = denver denver_to_nyc.end_point = new_york denver_to_nyc.num_points = 600 # wow, very densely sampled print("Done creating point sources") all_points = list(itertools.chain( abq_to_sd.points(), sd_to_sea.points(), denver_to_nyc.points() )) trajectory_assembler = AssembleTrajectoryFromPoints() trajectory_assembler.input = all_points trajectory_assembler.separation_time = timedelta(minutes=30) trajectory_assembler.separation_distance = 100 trajectory_assembler_minimum_length = 10 print("Done instantiating assembler") all_trajectories = list(trajectory_assembler.trajectories()) print("Assembler statistics: {} points, {} valid trajectories, {} invalid trajectories".format( trajectory_assembler.valid_trajectory_count, trajectory_assembler.invalid_trajectory_count, trajectory_assembler.points_processed_count )) print("Done creating trajectories. Found {}.".format(len(all_trajectories))) test_point_proximity = geomath.sanity_check_distance_less_than(1) def test_timestamp_proximity(time1, time2): return ( (time2 - time1).total_seconds() < 1 ) error_count = 0 if len(all_trajectories) != 3: sys.stdout.write('ERROR: test_trajectory_assembly: Expected 3 trajectories but got {}\n'.format(len(all_trajectories))) error_count += 1 if not test_point_proximity(all_trajectories[0][0], albuquerque): sys.stdout.write('ERROR: test_trajectory_assembly: Expected point 0 of first trajectory to be Albuquerque ({}) but it is instead {}\n'.format(albuquerque, str(all_trajectories[0][0]))) error_count += 1 if not test_point_proximity(all_trajectories[0][-1], san_diego1): sys.stdout.write('ERROR: test_trajectory_assembly: Expected last point of first trajectory to be San Diego ({}) but it is instead {}\n'.format(san_diego, str(all_trajectories[0][-1]))) error_count += 1 if not test_point_proximity(all_trajectories[1][0], san_diego2): sys.stdout.write('ERROR: test_trajectory_assembly: Expected point 0 of second trajectory to be San Diego ({}) but it is instead {}\n'.format(san_diego, str(all_trajectories[1][0]))) error_count += 1 if not test_point_proximity(all_trajectories[1][-1], seattle): sys.stdout.write('ERROR: test_trajectory_assembly: Expected last point of second trajectory to be Seattle ({}) but it is instead {}\n'.format(seattle, str(all_trajectories[1][-1]))) error_count += 1 if not test_point_proximity(all_trajectories[2][0], denver): sys.stdout.write('ERROR: test_trajectory_assembly: Expected first point of third trajectory to be Denver ({}) but it is instead {}\n'.format(denver, str(all_trajectories[2][0]))) error_count += 1 if not test_point_proximity(all_trajectories[2][-1], new_york): sys.stdout.write('ERROR: test_trajectory_assembly: Expected last point of third trajectory to be New York ({}) but it is instead {}\n'.format(new_york, str(all_trajectories[2][-1]))) error_count += 1 if not test_timestamp_proximity(all_trajectories[0][0].timestamp, albuquerque.timestamp): sys.stdout.write('ERROR: test_trajectory_assembly: Expected timestamp at beginning of trajectory 0 to be {} but it is instead {}\n'.format(albuquerque.timestamp, all_trajectories[0][0].timestamp)) error_count += 1 if not test_timestamp_proximity(all_trajectories[0][-1].timestamp, san_diego1.timestamp): sys.stdout.write('ERROR: test_trajectory_assembly: Expected timestamp at end of trajectory 0 to be {} but it is instead {}\n'.format(san_diego1.timestamp, all_trajectories[0][-1].timestamp)) error_count += 1 if not test_timestamp_proximity(all_trajectories[1][0].timestamp, san_diego2.timestamp): sys.stdout.write('ERROR: test_trajectory_assembly: Expected timestamp at beginning of trajectory 1 to be {} but it is instead {}\n'.format(san_diego2.timestamp, all_trajectories[1][0].timestamp)) error_count += 1 if not test_timestamp_proximity(all_trajectories[1][-1].timestamp, seattle.timestamp): sys.stdout.write('ERROR: test_trajectory_assembly: Expected end at beginning of trajectory 1 to be {} but it is instead {}\n'.format(seattle.timestamp, all_trajectories[1][-1].timestamp)) error_count += 1 if not test_timestamp_proximity(all_trajectories[2][0].timestamp, denver.timestamp): sys.stdout.write('ERROR: test_trajectory_assembly: Expected timestamp at beginning of trajectory 2 to be {} but it is instead {}\n'.format(denver.timestamp, all_trajectories[2][0].timestamp)) error_count += 1 if not test_timestamp_proximity(all_trajectories[2][-1].timestamp, new_york.timestamp): sys.stdout.write('ERROR: test_trajectory_assembly: Expected timestamp at end of trajectory 2 to be {} but it is instead {}\n'.format(new_york.timestamp, all_trajectories[2][-1].timestamp)) error_count += 1 print("Done checking proximities") return error_count