repo/traffic-intelligence: python/storage.py annotate

annotate python/storage.py @ 935:0e63a918a1ca

updated classify-objects

author	Nicolas Saunier <nicolas.saunier@polymtl.ca>
date	Fri, 14 Jul 2017 16:30:57 -0400
parents	c030f735c594
children	fbf12382f3f8

rev	line source
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1 #! /usr/bin/env python
208 d9855499fc88 Added functions to write trajectories through sqlite Francois Belisle <belisle.francois@gmail.com> parents: 204 diff changeset	2 # -- coding: utf-8 --
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	3 '''Various utilities to save and load data'''
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	4
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	5 import utils, moving, events, indicators, shutil
665 15e244d2a1b5 corrected bug with circular import for VideoFilenameAddable, moved to base module Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 664 diff changeset	6 from base import VideoFilenameAddable
343 74e437ab5f11 first version of indicator loading code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 342 diff changeset	7
813 ef8795dba5ed avoid crash for missing configuration filename Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 812 diff changeset	8 from os import path
919 7b3f2e0a2652 saving and loading prototype trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 918 diff changeset	9 from copy import copy
417 a2ff03a52b73 added basic logging capability for debugging Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 394 diff changeset	10 import sqlite3, logging
901 753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	11 from numpy import log, min as npmin, max as npmax, round as npround, array, sum as npsum, loadtxt, floor as npfloor, ceil as npceil, linalg
756 a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	12 from pandas import read_csv, merge
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	13
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	14
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	15 commentChar = '#'
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	16
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	17 delimiterChar = '%';
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	18
50 7b06d649122b re-arrangement Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 48 diff changeset	19 ngsimUserTypes = {'twowheels':1,
204 966c2cd2bd9f added code to load object trajectories (average of features) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 203 diff changeset	20 'car':2,
966c2cd2bd9f added code to load object trajectories (average of features) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 203 diff changeset	21 'truck':3}
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	22
908 b297525b2cbf added options to the prototype cluster algorithm, work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 901 diff changeset	23 tableNames = {'feature':'positions',
b297525b2cbf added options to the prototype cluster algorithm, work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 901 diff changeset	24 'object': 'objects',
b297525b2cbf added options to the prototype cluster algorithm, work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 901 diff changeset	25 'objectfeatures': 'positions'}
b297525b2cbf added options to the prototype cluster algorithm, work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 901 diff changeset	26
259 8ab76b95ee72 added code to save collision points and crossing zones in txt files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 246 diff changeset	27 #########################
8ab76b95ee72 added code to save collision points and crossing zones in txt files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 246 diff changeset	28 # Sqlite
8ab76b95ee72 added code to save collision points and crossing zones in txt files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 246 diff changeset	29 #########################
8ab76b95ee72 added code to save collision points and crossing zones in txt files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 246 diff changeset	30
491 343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	31 # utils
343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	32 def printDBError(error):
343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	33 print('DB Error: {}'.format(error))
343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	34
343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	35 def dropTables(connection, tableNames):
343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	36 'deletes the table with names in tableNames'
343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	37 try:
343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	38 cursor = connection.cursor()
343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	39 for tableName in tableNames:
343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	40 cursor.execute('DROP TABLE IF EXISTS '+tableName)
343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	41 except sqlite3.OperationalError as error:
343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	42 printDBError(error)
343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	43
917 89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	44 def deleteFromSqlite(filename, dataType):
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	45 'Deletes (drops) some tables in the filename depending on type of data'
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	46 if path.isfile(filename):
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	47 connection = sqlite3.connect(filename)
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	48 if dataType == 'object':
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	49 dropTables(connection, ['objects', 'objects_features'])
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	50 elif dataType == 'interaction':
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	51 dropTables(connection, ['interactions', 'indicators'])
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	52 elif dataType == 'bb':
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	53 dropTables(connection, ['bounding_boxes'])
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	54 elif dataType == 'pois':
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	55 dropTables(connection, ['gaussians2d', 'objects_pois'])
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	56 elif dataType == 'prototype':
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	57 dropTables(connection, ['prototypes', 'objects_prototypes'])
917 89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	58 else:
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	59 print('Unknown data type {} to delete from database'.format(dataType))
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	60 connection.close()
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	61 else:
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	62 print('{} does not exist'.format(filename))
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	63
919 7b3f2e0a2652 saving and loading prototype trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 918 diff changeset	64 def tableExists(connection, tableName):
740 867bab9f317a function to check the existence of tables Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 736 diff changeset	65 'indicates if the table exists in the database'
867bab9f317a function to check the existence of tables Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 736 diff changeset	66 try:
919 7b3f2e0a2652 saving and loading prototype trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 918 diff changeset	67 #connection = sqlite3.connect(filename)
740 867bab9f317a function to check the existence of tables Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 736 diff changeset	68 cursor = connection.cursor()
867bab9f317a function to check the existence of tables Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 736 diff changeset	69 cursor.execute('SELECT COUNT(*) FROM SQLITE_MASTER WHERE type = \'table\' AND name = \''+tableName+'\'')
867bab9f317a function to check the existence of tables Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 736 diff changeset	70 return cursor.fetchone()[0] == 1
867bab9f317a function to check the existence of tables Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 736 diff changeset	71 except sqlite3.OperationalError as error:
867bab9f317a function to check the existence of tables Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 736 diff changeset	72 printDBError(error)
867bab9f317a function to check the existence of tables Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 736 diff changeset	73
830 2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	74 def createTrajectoryTable(cursor, tableName):
919 7b3f2e0a2652 saving and loading prototype trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 918 diff changeset	75 if tableName.endswith('positions') or tableName.endswith('velocities'):
830 2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	76 cursor.execute("CREATE TABLE IF NOT EXISTS "+tableName+" (trajectory_id INTEGER, frame_number INTEGER, x_coordinate REAL, y_coordinate REAL, PRIMARY KEY(trajectory_id, frame_number))")
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	77 else:
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	78 print('Unallowed name {} for trajectory table'.format(tableName))
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	79
879 f9ea5083588e Initial commit of object DB storage update pstaub parents: 872 diff changeset	80 def createObjectsTable(cursor):
f9ea5083588e Initial commit of object DB storage update pstaub parents: 872 diff changeset	81 cursor.execute("CREATE TABLE IF NOT EXISTS objects (object_id INTEGER, road_user_type INTEGER, n_objects INTEGER, PRIMARY KEY(object_id))")
f9ea5083588e Initial commit of object DB storage update pstaub parents: 872 diff changeset	82
913 1cd878812529 work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 910 diff changeset	83 def createAssignmentTable(cursor, objectType1, objectType2, columnName1, columnName2):
1cd878812529 work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 910 diff changeset	84 cursor.execute("CREATE TABLE IF NOT EXISTS "+objectType1+"s_"+objectType2+"s ("+columnName1+" INTEGER, "+columnName1+" INTEGER, PRIMARY KEY("+columnName1+","+columnName2+"))")
1cd878812529 work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 910 diff changeset	85
1cd878812529 work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 910 diff changeset	86 def createObjectsFeaturesTable(cursor): # same as
1cd878812529 work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 910 diff changeset	87 cursor.execute("CREATE TABLE IF NOT EXISTS objects_features (object_id INTEGER, trajectory_id INTEGER, PRIMARY KEY(object_id, trajectory_id))")
879 f9ea5083588e Initial commit of object DB storage update pstaub parents: 872 diff changeset	88
f9ea5083588e Initial commit of object DB storage update pstaub parents: 872 diff changeset	89
830 2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	90 def createCurvilinearTrajectoryTable(cursor):
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	91 cursor.execute("CREATE TABLE IF NOT EXISTS curvilinear_positions (trajectory_id INTEGER, frame_number INTEGER, s_coordinate REAL, y_coordinate REAL, lane TEXT, PRIMARY KEY(trajectory_id, frame_number))")
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	92
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	93 def createFeatureCorrespondenceTable(cursor):
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	94 cursor.execute('CREATE TABLE IF NOT EXISTS feature_correspondences (trajectory_id INTEGER, source_dbname VARCHAR, db_trajectory_id INTEGER, PRIMARY KEY(trajectory_id))')
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	95
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	96 def createInteractionTable(cursor):
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	97 cursor.execute('CREATE TABLE IF NOT EXISTS interactions (id INTEGER PRIMARY KEY, object_id1 INTEGER, object_id2 INTEGER, first_frame_number INTEGER, last_frame_number INTEGER, FOREIGN KEY(object_id1) REFERENCES objects(id), FOREIGN KEY(object_id2) REFERENCES objects(id))')
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	98
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	99 def createIndicatorTable(cursor):
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	100 cursor.execute('CREATE TABLE IF NOT EXISTS indicators (interaction_id INTEGER, indicator_type INTEGER, frame_number INTEGER, value REAL, FOREIGN KEY(interaction_id) REFERENCES interactions(id), PRIMARY KEY(interaction_id, indicator_type, frame_number))')
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	101
2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	102 def insertTrajectoryQuery(tableName):
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	103 return "INSERT INTO "+tableName+" VALUES (?,?,?,?)"
831 a8ff35e6fb43 forgotten change Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 830 diff changeset	104
879 f9ea5083588e Initial commit of object DB storage update pstaub parents: 872 diff changeset	105 def insertObjectQuery():
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	106 return "INSERT INTO objects VALUES (?,?,?)"
879 f9ea5083588e Initial commit of object DB storage update pstaub parents: 872 diff changeset	107
f9ea5083588e Initial commit of object DB storage update pstaub parents: 872 diff changeset	108 def insertObjectFeatureQuery():
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	109 return "INSERT INTO objects_features VALUES (?,?)"
879 f9ea5083588e Initial commit of object DB storage update pstaub parents: 872 diff changeset	110
714 d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	111 def createIndex(connection, tableName, columnName, unique = False):
711 523eda2fafd4 added function to create index Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 704 diff changeset	112 '''Creates an index for the column in the table
523eda2fafd4 added function to create index Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 704 diff changeset	113 I will make querying with a condition on this column faster'''
523eda2fafd4 added function to create index Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 704 diff changeset	114 try:
714 d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	115 #connection = sqlite3.connect(filename)
711 523eda2fafd4 added function to create index Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 704 diff changeset	116 cursor = connection.cursor()
714 d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	117 s = "CREATE "
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	118 if unique:
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	119 s += "UNIQUE "
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	120 cursor.execute(s+"INDEX IF NOT EXISTS "+tableName+"_"+columnName+"_index ON "+tableName+"("+columnName+")")
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	121 connection.commit()
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	122 #connection.close()
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	123 except sqlite3.OperationalError as error:
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	124 printDBError(error)
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	125
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	126 def getNumberRowsTable(connection, tableName, columnName = None):
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	127 '''Returns the number of rows for the table
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	128 If columnName is not None, means we want the number of distinct values for that column
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	129 (otherwise, we can just count(*))'''
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	130 try:
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	131 cursor = connection.cursor()
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	132 if columnName is None:
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	133 cursor.execute("SELECT COUNT(*) from "+tableName)
711 523eda2fafd4 added function to create index Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 704 diff changeset	134 else:
714 d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	135 cursor.execute("SELECT COUNT(DISTINCT "+columnName+") from "+tableName)
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	136 return cursor.fetchone()[0]
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	137 except sqlite3.OperationalError as error:
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	138 printDBError(error)
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	139
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	140 def getMinMax(connection, tableName, columnName, minmax):
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	141 '''Returns max/min or both for given column in table
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	142 minmax must be string max, min or minmax'''
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	143 try:
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	144 cursor = connection.cursor()
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	145 if minmax == 'min' or minmax == 'max':
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	146 cursor.execute("SELECT "+minmax+"("+columnName+") from "+tableName)
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	147 elif minmax == 'minmax':
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	148 cursor.execute("SELECT MIN("+columnName+"), MAX("+columnName+") from "+tableName)
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	149 else:
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	150 print("Argument minmax unknown: {}".format(minmax))
d6c69d3d09e5 developing python performance tests for SQLite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 712 diff changeset	151 return cursor.fetchone()[0]
711 523eda2fafd4 added function to create index Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 704 diff changeset	152 except sqlite3.OperationalError as error:
523eda2fafd4 added function to create index Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 704 diff changeset	153 printDBError(error)
523eda2fafd4 added function to create index Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 704 diff changeset	154
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	155 def getObjectCriteria(objectNumbers):
585 aab2242ea88c minor modification of objectNumbers to None Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 584 diff changeset	156 if objectNumbers is None:
aab2242ea88c minor modification of objectNumbers to None Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 584 diff changeset	157 query = ''
aab2242ea88c minor modification of objectNumbers to None Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 584 diff changeset	158 elif type(objectNumbers) == int:
852 45a53542e046 updated unnecessary complicated query Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 850 diff changeset	159 query = '<= {0}'.format(objectNumbers-1)
263 c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	160 elif type(objectNumbers) == list:
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	161 query = 'in ('+', '.join([str(n) for n in objectNumbers])+')'
5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	162 else:
5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	163 print('objectNumbers {} are not a known type ({})'.format(objectNumbers, type(objectNumbers)))
5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	164 query = ''
263 c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	165 return query
c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	166
834 119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	167 def loadTrajectoriesFromTable(connection, tableName, trajectoryType, objectNumbers = None, timeStep = None):
239 93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	168 '''Loads trajectories (in the general sense) from the given table
93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	169 can be positions or velocities
93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	170
93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	171 returns a moving object'''
143 436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	172 cursor = connection.cursor()
204 966c2cd2bd9f added code to load object trajectories (average of features) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 203 diff changeset	173
215 5e2983b05d4e created first doctest tests for storage Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 212 diff changeset	174 try:
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	175 objectCriteria = getObjectCriteria(objectNumbers)
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	176 queryStatement = None
215 5e2983b05d4e created first doctest tests for storage Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 212 diff changeset	177 if trajectoryType == 'feature':
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	178 queryStatement = 'SELECT * from '+tableName
834 119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	179 if objectNumbers is not None and timeStep is not None:
119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	180 queryStatement += ' WHERE trajectory_id '+objectCriteria+' AND frame_number%{} = 0'.format(timeStep)
119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	181 elif objectNumbers is not None:
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	182 queryStatement += ' WHERE trajectory_id '+objectCriteria
834 119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	183 elif timeStep is not None:
119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	184 queryStatement += ' WHERE frame_number%{} = 0'.format(timeStep)
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	185 queryStatement += ' ORDER BY trajectory_id, frame_number'
215 5e2983b05d4e created first doctest tests for storage Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 212 diff changeset	186 elif trajectoryType == 'object':
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	187 queryStatement = 'SELECT OF.object_id, P.frame_number, avg(P.x_coordinate), avg(P.y_coordinate) from '+tableName+' P, objects_features OF WHERE P.trajectory_id = OF.trajectory_id'
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	188 if objectNumbers is not None:
834 119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	189 queryStatement += ' AND OF.object_id '+objectCriteria
119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	190 if timeStep is not None:
119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	191 queryStatement += ' AND P.frame_number%{} = 0'.format(timeStep)
712 21aeadcfbabb added script to test SQLite performance and impact of indices Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 711 diff changeset	192 queryStatement += ' GROUP BY OF.object_id, P.frame_number ORDER BY OF.object_id, P.frame_number'
588 c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	193 elif trajectoryType in ['bbtop', 'bbbottom']:
587 cf578ba866da added code to load bounding box corners as trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 586 diff changeset	194 if trajectoryType == 'bbtop':
cf578ba866da added code to load bounding box corners as trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 586 diff changeset	195 corner = 'top_left'
cf578ba866da added code to load bounding box corners as trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 586 diff changeset	196 elif trajectoryType == 'bbbottom':
cf578ba866da added code to load bounding box corners as trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 586 diff changeset	197 corner = 'bottom_right'
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	198 queryStatement = 'SELECT object_id, frame_number, x_'+corner+', y_'+corner+' FROM '+tableName
834 119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	199 if objectNumbers is not None and timeStep is not None:
119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	200 queryStatement += ' WHERE object_id '+objectCriteria+' AND frame_number%{} = 0'.format(timeStep)
119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	201 elif objectNumbers is not None:
712 21aeadcfbabb added script to test SQLite performance and impact of indices Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 711 diff changeset	202 queryStatement += ' WHERE object_id '+objectCriteria
834 119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	203 elif timeStep is not None:
119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	204 queryStatement += ' WHERE frame_number%{} = 0'.format(timeStep)
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	205 queryStatement += ' ORDER BY object_id, frame_number'
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	206 else:
777 ef6dd60be2e1 added function to save feature trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 773 diff changeset	207 print('Unknown trajectory type {}'.format(trajectoryType))
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	208 if queryStatement is not None:
417 a2ff03a52b73 added basic logging capability for debugging Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 394 diff changeset	209 cursor.execute(queryStatement)
a2ff03a52b73 added basic logging capability for debugging Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 394 diff changeset	210 logging.debug(queryStatement)
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	211 except sqlite3.OperationalError as error:
491 343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	212 printDBError(error)
215 5e2983b05d4e created first doctest tests for storage Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 212 diff changeset	213 return []
239 93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	214
143 436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	215 objId = -1
436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	216 obj = None
436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	217 objects = []
436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	218 for row in cursor:
436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	219 if row[0] != objId:
436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	220 objId = row[0]
834 119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	221 if obj is not None and (obj.length() == obj.positions.length() or (timeStep is not None and npceil(obj.length()/timeStep) == obj.positions.length())):
143 436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	222 objects.append(obj)
636 3058e00887bc removed all issues because of tests with None, using is instead of == or != Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 628 diff changeset	223 elif obj is not None:
588 c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	224 print('Object {} is missing {} positions'.format(obj.getNum(), int(obj.length())-obj.positions.length()))
143 436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	225 obj = moving.MovingObject(row[0], timeInterval = moving.TimeInterval(row[1], row[1]), positions = moving.Trajectory([[row[2]],[row[3]]]))
436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	226 else:
436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	227 obj.timeInterval.last = row[1]
436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	228 obj.positions.addPositionXY(row[2],row[3])
436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	229
834 119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	230 if obj is not None and (obj.length() == obj.positions.length() or (timeStep is not None and npceil(obj.length()/timeStep) == obj.positions.length())):
143 436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	231 objects.append(obj)
636 3058e00887bc removed all issues because of tests with None, using is instead of == or != Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 628 diff changeset	232 elif obj is not None:
588 c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	233 print('Object {} is missing {} positions'.format(obj.getNum(), int(obj.length())-obj.positions.length()))
143 436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	234
239 93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	235 return objects
93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	236
850 c724a51d4f5f modified function signature Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 849 diff changeset	237 def loadUserTypesFromTable(cursor, objectNumbers):
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	238 objectCriteria = getObjectCriteria(objectNumbers)
910 b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	239 queryStatement = 'SELECT object_id, road_user_type FROM objects'
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	240 if objectNumbers is not None:
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	241 queryStatement += ' WHERE object_id '+objectCriteria
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	242 cursor.execute(queryStatement)
588 c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	243 userTypes = {}
c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	244 for row in cursor:
c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	245 userTypes[row[0]] = row[1]
c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	246 return userTypes
c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	247
919 7b3f2e0a2652 saving and loading prototype trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 918 diff changeset	248 def loadTrajectoriesFromSqlite(filename, trajectoryType, objectNumbers = None, withFeatures = False, timeStep = None, tablePrefix = None):
777 ef6dd60be2e1 added function to save feature trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 773 diff changeset	249 '''Loads the trajectories (in the general sense,
ef6dd60be2e1 added function to save feature trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 773 diff changeset	250 either features, objects (feature groups) or bounding box series)
ef6dd60be2e1 added function to save feature trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 773 diff changeset	251 The number loaded is either the first objectNumbers objects,
ef6dd60be2e1 added function to save feature trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 773 diff changeset	252 or the indices in objectNumbers from the database'''
588 c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	253 connection = sqlite3.connect(filename)
239 93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	254
919 7b3f2e0a2652 saving and loading prototype trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 918 diff changeset	255 if tablePrefix is None:
7b3f2e0a2652 saving and loading prototype trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 918 diff changeset	256 prefix = ''
7b3f2e0a2652 saving and loading prototype trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 918 diff changeset	257 else:
7b3f2e0a2652 saving and loading prototype trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 918 diff changeset	258 prefix = tablePrefix + '_'
7b3f2e0a2652 saving and loading prototype trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 918 diff changeset	259 objects = loadTrajectoriesFromTable(connection, prefix+'positions', trajectoryType, objectNumbers, timeStep)
7b3f2e0a2652 saving and loading prototype trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 918 diff changeset	260 objectVelocities = loadTrajectoriesFromTable(connection, prefix+'velocities', trajectoryType, objectNumbers, timeStep)
239 93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	261
93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	262 if len(objectVelocities) > 0:
93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	263 for o,v in zip(objects, objectVelocities):
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	264 if o.getNum() == v.getNum():
239 93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	265 o.velocities = v.positions
343 74e437ab5f11 first version of indicator loading code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 342 diff changeset	266 o.velocities.duplicateLastPosition() # avoid having velocity shorter by one position than positions
239 93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	267 else:
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	268 print('Could not match positions {0} with velocities {1}'.format(o.getNum(), v.getNum()))
239 93c26e45efd8 modified functions to read velocities from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 235 diff changeset	269
263 c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	270 if trajectoryType == 'object':
c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	271 cursor = connection.cursor()
c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	272 try:
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	273 # attribute feature numbers to objects
852 45a53542e046 updated unnecessary complicated query Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 850 diff changeset	274 queryStatement = 'SELECT trajectory_id, object_id FROM objects_features'
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	275 if objectNumbers is not None:
910 b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	276 queryStatement += ' WHERE object_id '+getObjectCriteria(objectNumbers)
852 45a53542e046 updated unnecessary complicated query Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 850 diff changeset	277 queryStatement += ' ORDER BY object_id' # order is important to group all features per object
910 b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	278 logging.debug(queryStatement)
417 a2ff03a52b73 added basic logging capability for debugging Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 394 diff changeset	279 cursor.execute(queryStatement)
263 c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	280
c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	281 featureNumbers = {}
c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	282 for row in cursor:
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	283 objId = row[1]
a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	284 if objId not in featureNumbers:
263 c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	285 featureNumbers[objId] = [row[0]]
c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	286 else:
c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	287 featureNumbers[objId].append(row[0])
c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	288
c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	289 for obj in objects:
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	290 obj.featureNumbers = featureNumbers[obj.getNum()]
a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	291
a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	292 # load userType
850 c724a51d4f5f modified function signature Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 849 diff changeset	293 userTypes = loadUserTypesFromTable(cursor, objectNumbers)
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	294 for obj in objects:
a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	295 obj.userType = userTypes[obj.getNum()]
685 94b291a5f933 several updates for display Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 680 diff changeset	296
94b291a5f933 several updates for display Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 680 diff changeset	297 if withFeatures:
94b291a5f933 several updates for display Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 680 diff changeset	298 nFeatures = 0
94b291a5f933 several updates for display Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 680 diff changeset	299 for obj in objects:
94b291a5f933 several updates for display Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 680 diff changeset	300 nFeatures = max(nFeatures, max(obj.featureNumbers))
834 119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	301 features = loadTrajectoriesFromSqlite(filename, 'feature', nFeatures+1, timeStep = timeStep)
685 94b291a5f933 several updates for display Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 680 diff changeset	302 for obj in objects:
94b291a5f933 several updates for display Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 680 diff changeset	303 obj.setFeatures(features)
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	304
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	305 except sqlite3.OperationalError as error:
491 343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	306 printDBError(error)
588 c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	307 objects = []
263 c71540470057 reorganized loading trajectories and objects, added loading feature numbers for objects (the set of features grouped as one moving object Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 259 diff changeset	308
143 436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	309 connection.close()
436b87d4b992 wrote the code to load positions from sqlite database Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 78 diff changeset	310 return objects
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	311
910 b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	312 def loadObjectFeatureFrameNumbers(filename, objectNumbers = None):
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	313 'Loads the feature frame numbers for each object'
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	314 connection = sqlite3.connect(filename)
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	315 cursor = connection.cursor()
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	316 try:
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	317 queryStatement = 'SELECT OF.object_id, TL.trajectory_id, TL.length FROM (SELECT trajectory_id, max(frame_number)-min(frame_number) AS length FROM positions GROUP BY trajectory_id) TL, objects_features OF WHERE TL.trajectory_id = OF.trajectory_id'
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	318 if objectNumbers is not None:
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	319 queryStatement += ' AND object_id '+getObjectCriteria(objectNumbers)
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	320 queryStatement += ' ORDER BY OF.object_id, TL.length DESC'
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	321 logging.debug(queryStatement)
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	322 cursor.execute(queryStatement)
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	323 objectFeatureNumbers = {}
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	324 for row in cursor:
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	325 objId = row[0]
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	326 if objId in objectFeatureNumbers:
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	327 objectFeatureNumbers[objId].append(row[1])
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	328 else:
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	329 objectFeatureNumbers[objId] = [row[1]]
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	330 return objectFeatureNumbers
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	331 except sqlite3.OperationalError as error:
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	332 printDBError(error)
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	333 return None
b58a1061a717 loading is faster for longest object features Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 909 diff changeset	334
780 1b22d81ef5ff cleaned and checked storage with functions for curvilinear trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 778 diff changeset	335 def addCurvilinearTrajectoriesFromSqlite(filename, objects):
1b22d81ef5ff cleaned and checked storage with functions for curvilinear trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 778 diff changeset	336 '''Adds curvilinear positions (s_coordinate, y_coordinate, lane)
1b22d81ef5ff cleaned and checked storage with functions for curvilinear trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 778 diff changeset	337 from a database to an existing MovingObject dict (indexed by each objects's num)'''
1b22d81ef5ff cleaned and checked storage with functions for curvilinear trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 778 diff changeset	338 connection = sqlite3.connect(filename)
778 bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	339 cursor = connection.cursor()
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	340
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	341 try:
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	342 cursor.execute('SELECT * from curvilinear_positions order by trajectory_id, frame_number')
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	343 except sqlite3.OperationalError as error:
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	344 printDBError(error)
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	345 return []
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	346
784 30bd0f2223b7 tweaking curvilinear trajectory loading Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 783 diff changeset	347 missingObjectNumbers = []
778 bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	348 objNum = None
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	349 for row in cursor:
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	350 if objNum != row[0]:
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	351 objNum = row[0]
784 30bd0f2223b7 tweaking curvilinear trajectory loading Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 783 diff changeset	352 if objNum in objects:
783 f1040d739bbf added check to avoid issue if there are more cuvilinear trajectories than objects Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 780 diff changeset	353 objects[objNum].curvilinearPositions = moving.CurvilinearTrajectory()
f1040d739bbf added check to avoid issue if there are more cuvilinear trajectories than objects Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 780 diff changeset	354 else:
784 30bd0f2223b7 tweaking curvilinear trajectory loading Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 783 diff changeset	355 missingObjectNumbers.append(objNum)
30bd0f2223b7 tweaking curvilinear trajectory loading Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 783 diff changeset	356 if objNum in objects:
30bd0f2223b7 tweaking curvilinear trajectory loading Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 783 diff changeset	357 objects[objNum].curvilinearPositions.addPositionSYL(row[2],row[3],row[4])
30bd0f2223b7 tweaking curvilinear trajectory loading Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 783 diff changeset	358 if len(missingObjectNumbers) > 0:
30bd0f2223b7 tweaking curvilinear trajectory loading Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 783 diff changeset	359 print('List of missing objects to attach corresponding curvilinear trajectories: {}'.format(missingObjectNumbers))
778 bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	360
918 3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	361 def saveTrajectoriesToTable(connection, objects, trajectoryType, tablePrefix = None):
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	362 'Saves trajectories in table tableName'
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	363 cursor = connection.cursor()
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	364 # Parse feature and/or object structure and commit to DB
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	365 if(trajectoryType == 'feature' or trajectoryType == 'object'):
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	366 # Extract features from objects
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	367 if trajectoryType == 'object':
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	368 features = []
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	369 for obj in objects:
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	370 if obj.hasFeatures():
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	371 features += obj.getFeatures()
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	372 if len(features) == 0:
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	373 print('Warning, objects have no features') # todo save centroid trajectories?
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	374 elif trajectoryType == 'feature':
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	375 features = objects
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	376 # Setup feature queries
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	377 if tablePrefix is None:
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	378 prefix = ''
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	379 else:
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	380 prefix = tablePrefix+'_'
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	381 createTrajectoryTable(cursor, prefix+"positions")
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	382 createTrajectoryTable(cursor, prefix+"velocities")
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	383 positionQuery = insertTrajectoryQuery(prefix+"positions")
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	384 velocityQuery = insertTrajectoryQuery(prefix+"velocities")
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	385 # Setup object queries
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	386 if trajectoryType == 'object':
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	387 createObjectsTable(cursor)
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	388 createObjectsFeaturesTable(cursor)
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	389 objectQuery = insertObjectQuery()
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	390 objectFeatureQuery = insertObjectFeatureQuery()
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	391 for feature in features:
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	392 num = feature.getNum()
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	393 frameNum = feature.getFirstInstant()
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	394 for p in feature.getPositions():
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	395 cursor.execute(positionQuery, (num, frameNum, p.x, p.y))
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	396 frameNum += 1
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	397 velocities = feature.getVelocities()
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	398 if velocities is not None:
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	399 frameNum = feature.getFirstInstant()
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	400 for v in velocities[:-1]:
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	401 cursor.execute(velocityQuery, (num, frameNum, v.x, v.y))
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	402 frameNum += 1
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	403 if trajectoryType == 'object':
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	404 for obj in objects:
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	405 if obj.hasFeatures():
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	406 for feature in obj.getFeatures():
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	407 featureNum = feature.getNum()
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	408 cursor.execute(objectFeatureQuery, (obj.getNum(), featureNum))
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	409 cursor.execute(objectQuery, (obj.getNum(), obj.getUserType(), 1))
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	410 # Parse curvilinear position structure
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	411 elif(trajectoryType == 'curvilinear'):
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	412 createCurvilinearTrajectoryTable(cursor)
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	413 curvilinearQuery = "INSERT INTO curvilinear_positions VALUES (?,?,?,?,?)"
918 3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	414 for obj in objects:
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	415 num = obj.getNum()
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	416 frameNum = obj.getFirstInstant()
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	417 for p in obj.getCurvilinearPositions():
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	418 cursor.execute(curvilinearQuery, (num, frameNum, p[0], p[1], p[2]))
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	419 frameNum += 1
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	420 else:
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	421 print('Unknown trajectory type {}'.format(trajectoryType))
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	422 connection.commit()
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	423
3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	424 def saveTrajectoriesToSqlite(outputFilename, objects, trajectoryType):
829 0ddcc41663f5 renaming Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 813 diff changeset	425 '''Writes features, ie the trajectory positions (and velocities if exist)
778 bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	426 with their instants to a specified sqlite file
780 1b22d81ef5ff cleaned and checked storage with functions for curvilinear trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 778 diff changeset	427 Either feature positions (and velocities if they exist)
879 f9ea5083588e Initial commit of object DB storage update pstaub parents: 872 diff changeset	428 or curvilinear positions will be saved at a time'''
778 bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	429
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	430 connection = sqlite3.connect(outputFilename)
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	431 try:
918 3a06007a4bb7 modularized save trajectories, added slice to Trajectory, etc Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 917 diff changeset	432 saveTrajectoriesToTable(connection, objects, trajectoryType, None)
778 bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	433 except sqlite3.OperationalError as error:
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	434 printDBError(error)
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	435 connection.close()
bd684e57c431 integrated code from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 777 diff changeset	436
834 119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	437 def loadBBMovingObjectsFromSqlite(filename, objectType = 'bb', objectNumbers = None, timeStep = None):
768 f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	438 '''Loads bounding box moving object from an SQLite
f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	439 (format of SQLite output by the ground truth annotation tool
f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	440 or Urban Tracker
588 c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	441
768 f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	442 Load descriptions?'''
f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	443 connection = sqlite3.connect(filename)
f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	444 objects = []
f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	445
f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	446 if objectType == 'bb':
834 119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	447 topCorners = loadTrajectoriesFromTable(connection, 'bounding_boxes', 'bbtop', objectNumbers, timeStep)
119c4efe6398 added option to load subsampled trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 831 diff changeset	448 bottomCorners = loadTrajectoriesFromTable(connection, 'bounding_boxes', 'bbbottom', objectNumbers, timeStep)
850 c724a51d4f5f modified function signature Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 849 diff changeset	449 userTypes = loadUserTypesFromTable(connection.cursor(), objectNumbers) # string format is same as object
588 c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	450
c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	451 for t, b in zip(topCorners, bottomCorners):
c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	452 num = t.getNum()
c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	453 if t.getNum() == b.getNum():
768 f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	454 annotation = moving.BBMovingObject(num, t.getTimeInterval(), t, b, userTypes[num])
f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	455 objects.append(annotation)
588 c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	456 else:
768 f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	457 print ('Unknown type of bounding box {}'.format(objectType))
588 c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	458
c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	459 connection.close()
768 f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	460 return objects
588 c5406edbcf12 added loading ground truth annotations (ground truth) from polytrack format Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 587 diff changeset	461
341 2f39c4ed0b62 first version of indicator saving to sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 340 diff changeset	462 def saveInteraction(cursor, interaction):
2f39c4ed0b62 first version of indicator saving to sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 340 diff changeset	463 roadUserNumbers = list(interaction.getRoadUserNumbers())
343 74e437ab5f11 first version of indicator loading code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 342 diff changeset	464 cursor.execute('INSERT INTO interactions VALUES({}, {}, {}, {}, {})'.format(interaction.getNum(), roadUserNumbers[0], roadUserNumbers[1], interaction.getFirstInstant(), interaction.getLastInstant()))
341 2f39c4ed0b62 first version of indicator saving to sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 340 diff changeset	465
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	466 def saveInteractionsToSqlite(filename, interactions):
340 1046b7346886 work in progress on storing indicator values Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 335 diff changeset	467 'Saves the interactions in the table'
1046b7346886 work in progress on storing indicator values Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 335 diff changeset	468 connection = sqlite3.connect(filename)
1046b7346886 work in progress on storing indicator values Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 335 diff changeset	469 cursor = connection.cursor()
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	470 try:
4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	471 createInteractionTable(cursor)
4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	472 for inter in interactions:
4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	473 saveInteraction(cursor, inter)
4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	474 except sqlite3.OperationalError as error:
491 343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	475 printDBError(error)
341 2f39c4ed0b62 first version of indicator saving to sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 340 diff changeset	476 connection.commit()
340 1046b7346886 work in progress on storing indicator values Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 335 diff changeset	477 connection.close()
1046b7346886 work in progress on storing indicator values Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 335 diff changeset	478
341 2f39c4ed0b62 first version of indicator saving to sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 340 diff changeset	479 def saveIndicator(cursor, interactionNum, indicator):
2f39c4ed0b62 first version of indicator saving to sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 340 diff changeset	480 for instant in indicator.getTimeInterval():
2f39c4ed0b62 first version of indicator saving to sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 340 diff changeset	481 if indicator[instant]:
2f39c4ed0b62 first version of indicator saving to sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 340 diff changeset	482 cursor.execute('INSERT INTO indicators VALUES({}, {}, {}, {})'.format(interactionNum, events.Interaction.indicatorNameToIndices[indicator.getName()], instant, indicator[instant]))
2f39c4ed0b62 first version of indicator saving to sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 340 diff changeset	483
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	484 def saveIndicatorsToSqlite(filename, interactions, indicatorNames = events.Interaction.indicatorNames):
340 1046b7346886 work in progress on storing indicator values Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 335 diff changeset	485 'Saves the indicator values in the table'
1046b7346886 work in progress on storing indicator values Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 335 diff changeset	486 connection = sqlite3.connect(filename)
341 2f39c4ed0b62 first version of indicator saving to sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 340 diff changeset	487 cursor = connection.cursor()
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	488 try:
4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	489 createInteractionTable(cursor)
830 2a5856961933 first working version of feature merging (works with feature grouping) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 829 diff changeset	490 createIndicatorTable(cursor)
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	491 for inter in interactions:
4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	492 saveInteraction(cursor, inter)
4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	493 for indicatorName in indicatorNames:
4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	494 indicator = inter.getIndicator(indicatorName)
636 3058e00887bc removed all issues because of tests with None, using is instead of == or != Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 628 diff changeset	495 if indicator is not None:
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	496 saveIndicator(cursor, inter.getNum(), indicator)
4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	497 except sqlite3.OperationalError as error:
491 343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	498 printDBError(error)
341 2f39c4ed0b62 first version of indicator saving to sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 340 diff changeset	499 connection.commit()
2f39c4ed0b62 first version of indicator saving to sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 340 diff changeset	500 connection.close()
340 1046b7346886 work in progress on storing indicator values Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 335 diff changeset	501
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	502 def loadInteractionsFromSqlite(filename):
482 f6415f012640 adding functionalities (save images directly to display trajectories to create movies Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 417 diff changeset	503 '''Loads interaction and their indicators
343 74e437ab5f11 first version of indicator loading code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 342 diff changeset	504
74e437ab5f11 first version of indicator loading code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 342 diff changeset	505 TODO choose the interactions to load'''
74e437ab5f11 first version of indicator loading code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 342 diff changeset	506 interactions = []
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	507 connection = sqlite3.connect(filename)
4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	508 cursor = connection.cursor()
4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	509 try:
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	510 cursor.execute('SELECT INT.id, INT.object_id1, INT.object_id2, INT.first_frame_number, INT.last_frame_number, IND.indicator_type, IND.frame_number, IND.value from interactions INT, indicators IND WHERE INT.id = IND.interaction_id ORDER BY INT.id, IND.indicator_type, IND.frame_number')
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	511 interactionNum = -1
4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	512 indicatorTypeNum = -1
343 74e437ab5f11 first version of indicator loading code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 342 diff changeset	513 tmpIndicators = {}
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	514 for row in cursor:
628 977407c9f815 corrected bugs in loading interactions (index shifted) and added functionalities to play/plot interactions Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 619 diff changeset	515 if row[0] != interactionNum:
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	516 interactionNum = row[0]
628 977407c9f815 corrected bugs in loading interactions (index shifted) and added functionalities to play/plot interactions Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 619 diff changeset	517 interactions.append(events.Interaction(interactionNum, moving.TimeInterval(row[3],row[4]), row[1], row[2]))
977407c9f815 corrected bugs in loading interactions (index shifted) and added functionalities to play/plot interactions Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 619 diff changeset	518 interactions[-1].indicators = {}
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	519 if indicatorTypeNum != row[5] or row[0] != interactionNum:
695 957126bfb456 corrected bug with indicator loading(also now correctly loading mostsevereismax) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 693 diff changeset	520 indicatorTypeNum = row[5]
628 977407c9f815 corrected bugs in loading interactions (index shifted) and added functionalities to play/plot interactions Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 619 diff changeset	521 indicatorName = events.Interaction.indicatorNames[indicatorTypeNum]
977407c9f815 corrected bugs in loading interactions (index shifted) and added functionalities to play/plot interactions Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 619 diff changeset	522 indicatorValues = {row[6]:row[7]}
695 957126bfb456 corrected bug with indicator loading(also now correctly loading mostsevereismax) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 693 diff changeset	523 interactions[-1].indicators[indicatorName] = indicators.SeverityIndicator(indicatorName, indicatorValues, mostSevereIsMax = not indicatorName in events.Interaction.timeIndicators)
343 74e437ab5f11 first version of indicator loading code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 342 diff changeset	524 else:
74e437ab5f11 first version of indicator loading code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 342 diff changeset	525 indicatorValues[row[6]] = row[7]
693 5ee22bf7e4d5 corrected bug when loading indicator time intervals and updated how queries are created for better legibility Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 690 diff changeset	526 interactions[-1].indicators[indicatorName].timeInterval.last = row[6]
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	527 except sqlite3.OperationalError as error:
491 343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	528 printDBError(error)
343 74e437ab5f11 first version of indicator loading code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 342 diff changeset	529 return []
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	530 connection.close()
343 74e437ab5f11 first version of indicator loading code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 342 diff changeset	531 return interactions
74e437ab5f11 first version of indicator loading code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 342 diff changeset	532 # load first and last object instants
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	533 # CREATE TEMP TABLE IF NOT EXISTS object_instants AS SELECT OF.object_id, min(frame_number) as first_instant, max(frame_number) as last_instant from positions P, objects_features OF WHERE P.trajectory_id = OF.trajectory_id group by OF.object_id order by OF.object_id
342 4d69486869a5 work on loading indicators Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 341 diff changeset	534
390 12be4a0cb9aa sql code to create bounding boxes in image space Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 382 diff changeset	535 def createBoundingBoxTable(filename, invHomography = None):
482 f6415f012640 adding functionalities (save images directly to display trajectories to create movies Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 417 diff changeset	536 '''Create the table to store the object bounding boxes in image space
390 12be4a0cb9aa sql code to create bounding boxes in image space Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 382 diff changeset	537 '''
12be4a0cb9aa sql code to create bounding boxes in image space Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 382 diff changeset	538 connection = sqlite3.connect(filename)
12be4a0cb9aa sql code to create bounding boxes in image space Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 382 diff changeset	539 cursor = connection.cursor()
12be4a0cb9aa sql code to create bounding boxes in image space Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 382 diff changeset	540 try:
394 6567fee37c16 renamed table for bounding boxes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 393 diff changeset	541 cursor.execute('CREATE TABLE IF NOT EXISTS bounding_boxes (object_id INTEGER, frame_number INTEGER, x_top_left REAL, y_top_left REAL, x_bottom_right REAL, y_bottom_right REAL, PRIMARY KEY(object_id, frame_number))')
6567fee37c16 renamed table for bounding boxes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 393 diff changeset	542 cursor.execute('INSERT INTO bounding_boxes SELECT object_id, frame_number, min(x), min(y), max(x), max(y) from '
390 12be4a0cb9aa sql code to create bounding boxes in image space Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 382 diff changeset	543 '(SELECT object_id, frame_number, (x{}+y{}+{})/w as x, (x{}+y{}+{})/w as y from '
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	544 '(SELECT OF.object_id, P.frame_number, P.x_coordinate as x, P.y_coordinate as y, P.x_coordinate{}+P.y_coordinate{}+{} as w from positions P, objects_features OF WHERE P.trajectory_id = OF.trajectory_id)) '.format(invHomography[0,0], invHomography[0,1], invHomography[0,2], invHomography[1,0], invHomography[1,1], invHomography[1,2], invHomography[2,0], invHomography[2,1], invHomography[2,2])+
390 12be4a0cb9aa sql code to create bounding boxes in image space Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 382 diff changeset	545 'GROUP BY object_id, frame_number')
12be4a0cb9aa sql code to create bounding boxes in image space Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 382 diff changeset	546 except sqlite3.OperationalError as error:
491 343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	547 printDBError(error)
390 12be4a0cb9aa sql code to create bounding boxes in image space Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 382 diff changeset	548 connection.commit()
12be4a0cb9aa sql code to create bounding boxes in image space Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 382 diff changeset	549 connection.close()
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	550
586 ff4f0ce46ca6 modified name for loading bounding boxes (only for display) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 585 diff changeset	551 def loadBoundingBoxTableForDisplay(filename):
647 458890c0437c cleaned functions to lead bounding boxes (for display) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 646 diff changeset	552 '''Loads bounding boxes from bounding_boxes table for display over trajectories'''
393 eaf7765221d9 added code to create bounding boxes and display them (inc scripts) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 390 diff changeset	553 connection = sqlite3.connect(filename)
eaf7765221d9 added code to create bounding boxes and display them (inc scripts) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 390 diff changeset	554 cursor = connection.cursor()
eaf7765221d9 added code to create bounding boxes and display them (inc scripts) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 390 diff changeset	555 boundingBoxes = {} # list of bounding boxes for each instant
eaf7765221d9 added code to create bounding boxes and display them (inc scripts) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 390 diff changeset	556 try:
394 6567fee37c16 renamed table for bounding boxes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 393 diff changeset	557 cursor.execute('SELECT name FROM sqlite_master WHERE type=\'table\' AND name=\'bounding_boxes\'')
768 f8e0a8ea8402 updated the bounding box code (the name so that it is general for ground truth and UT) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 759 diff changeset	558 result = cursor.fetchall()
393 eaf7765221d9 added code to create bounding boxes and display them (inc scripts) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 390 diff changeset	559 if len(result) > 0:
394 6567fee37c16 renamed table for bounding boxes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 393 diff changeset	560 cursor.execute('SELECT * FROM bounding_boxes')
393 eaf7765221d9 added code to create bounding boxes and display them (inc scripts) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 390 diff changeset	561 for row in cursor:
eaf7765221d9 added code to create bounding boxes and display them (inc scripts) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 390 diff changeset	562 boundingBoxes.setdefault(row[1], []).append([moving.Point(row[2], row[3]), moving.Point(row[4], row[5])])
eaf7765221d9 added code to create bounding boxes and display them (inc scripts) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 390 diff changeset	563 except sqlite3.OperationalError as error:
491 343cfd185ca6 minor changes and reaarrangements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 482 diff changeset	564 printDBError(error)
393 eaf7765221d9 added code to create bounding boxes and display them (inc scripts) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 390 diff changeset	565 return boundingBoxes
eaf7765221d9 added code to create bounding boxes and display them (inc scripts) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 390 diff changeset	566 connection.close()
eaf7765221d9 added code to create bounding boxes and display them (inc scripts) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 390 diff changeset	567 return boundingBoxes
eaf7765221d9 added code to create bounding boxes and display them (inc scripts) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 390 diff changeset	568
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	569 #########################
917 89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	570 # saving and loading for scene interpretation: POIs and Prototypes
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	571 #########################
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	572
921 630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	573 def savePrototypesToSqlite(filename, prototypes):
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	574 '''save the prototypes (a prototype is defined by a filename, a number and type'''
917 89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	575 connection = sqlite3.connect(filename)
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	576 cursor = connection.cursor()
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	577 try:
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	578 cursor.execute('CREATE TABLE IF NOT EXISTS prototypes (prototype_filename VARCHAR, prototype_id INTEGER, trajectory_type VARCHAR CHECK (trajectory_type IN (\"feature\", \"object\")), nmatchings INTEGER, PRIMARY KEY (prototype_filename, prototype_id, trajectory_type))')
921 630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	579 for p in prototypes:
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	580 cursor.execute('INSERT INTO prototypes VALUES(?,?,?,?)', (p.getFilename(), p.getNum(), p.getTrajectoryType(), p.getNMatchings()))
917 89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	581 except sqlite3.OperationalError as error:
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	582 printDBError(error)
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	583 connection.commit()
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	584 connection.close()
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	585
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	586 def savePrototypeAssignmentsToSqlite(filename, objects, labels, prototypes):
c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	587 connection = sqlite3.connect(filename)
c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	588 cursor = connection.cursor()
c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	589 try:
c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	590 cursor.execute('CREATE TABLE IF NOT EXISTS objects_prototypes (object_id INTEGER, prototype_filename VARCHAR, prototype_id INTEGER, trajectory_type VARCHAR CHECK (trajectory_type IN (\"feature\", \"object\")), PRIMARY KEY(object_id, prototype_filename, prototype_id, trajectory_type))')
c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	591 for obj, label in zip(objects, labels):
c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	592 proto = prototypes[label]
c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	593 cursor.execute('INSERT INTO objects_prototypes VALUES(?,?,?,?)', (obj.getNum(), proto.getFilename(), proto.getNum(), proto.getTrajectoryType()))
c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	594 except sqlite3.OperationalError as error:
c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	595 printDBError(error)
c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	596 connection.commit()
c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	597 connection.close()
917 89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	598
920 499154254f37 improved prototype loading Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 919 diff changeset	599 def loadPrototypesFromSqlite(filename, withTrajectories = True):
917 89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	600 'Loads prototype ids and matchings (if stored)'
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	601 connection = sqlite3.connect(filename)
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	602 cursor = connection.cursor()
921 630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	603 prototypes = []
920 499154254f37 improved prototype loading Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 919 diff changeset	604 objects = []
917 89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	605 try:
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	606 cursor.execute('SELECT * FROM prototypes')
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	607 for row in cursor:
921 630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	608 prototypes.append(moving.Prototype(row[0], row[1], row[2], row[3]))
920 499154254f37 improved prototype loading Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 919 diff changeset	609 if withTrajectories:
921 630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	610 for p in prototypes:
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	611 p.setMovingObject(loadTrajectoriesFromSqlite(p.getFilename(), p.getTrajectoryType(), [p.getNum()])[0])
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	612 # loadingInformation = {} # complicated slightly optimized
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	613 # for p in prototypes:
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	614 # dbfn = p.getFilename()
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	615 # trajType = p.getTrajectoryType()
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	616 # if (dbfn, trajType) in loadingInformation:
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	617 # loadingInformation[(dbfn, trajType)].append(p)
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	618 # else:
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	619 # loadingInformation[(dbfn, trajType)] = [p]
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	620 # for k, v in loadingInformation.iteritems():
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	621 # objects += loadTrajectoriesFromSqlite(k[0], k[1], [p.getNum() for p in v])
917 89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	622 except sqlite3.OperationalError as error:
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	623 printDBError(error)
89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	624 connection.close()
921 630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	625 if len(set([p.getTrajectoryType() for p in prototypes])) > 1:
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	626 print('Different types of prototypes in database ({}).'.format(set([p.getTrajectoryType() for p in prototypes])))
630934595871 work in progress with prototype class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 920 diff changeset	627 return prototypes
917 89cc05867c4c reorg and work in progress Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 915 diff changeset	628
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	629 def savePOIsToSqlite(filename, gmm, gmmType, gmmId):
871 6db83beb5350 work in progress to update gaussian mixtures Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 862 diff changeset	630 '''Saves a Gaussian mixture model (of class sklearn.mixture.GaussianMixture)
6db83beb5350 work in progress to update gaussian mixtures Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 862 diff changeset	631 gmmType is a type of GaussianMixture, learnt either from beginnings or ends of trajectories'''
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	632 connection = sqlite3.connect(filename)
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	633 cursor = connection.cursor()
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	634 if gmmType not in ['beginning', 'end']:
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	635 print('Unknown POI type {}. Exiting'.format(gmmType))
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	636 import sys
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	637 sys.exit()
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	638 try:
872 c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	639 cursor.execute('CREATE TABLE IF NOT EXISTS gaussians2d (poi_id INTEGER, id INTEGER, type VARCHAR, x_center REAL, y_center REAL, covariance VARCHAR, covariance_type VARCHAR, weight, precisions_cholesky VARCHAR, PRIMARY KEY(poi_id, id))')
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	640 for i in xrange(gmm.n_components):
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	641 cursor.execute('INSERT INTO gaussians2d VALUES(?,?,?,?,?,?,?,?,?)', (gmmId, i, gmmType, gmm.means_[i][0], gmm.means_[i][1], str(gmm.covariances_[i].tolist()), gmm.covariance_type, gmm.weights_[i], str(gmm.precisions_cholesky_[i].tolist())))
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	642 connection.commit()
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	643 except sqlite3.OperationalError as error:
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	644 printDBError(error)
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	645 connection.close()
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	646
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	647 def savePOIAssignmentsToSqlite(filename, objects):
915 13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	648 'save the od fields of objects'
13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	649 connection = sqlite3.connect(filename)
13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	650 cursor = connection.cursor()
13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	651 try:
13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	652 cursor.execute('CREATE TABLE IF NOT EXISTS objects_pois (object_id INTEGER, origin_poi_id INTEGER, destination_poi_id INTEGER, PRIMARY KEY(object_id))')
13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	653 for o in objects:
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	654 cursor.execute('INSERT INTO objects_pois VALUES(?,?,?)', (o.getNum(), o.od[0], o.od[1]))
915 13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	655 connection.commit()
13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	656 except sqlite3.OperationalError as error:
13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	657 printDBError(error)
13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	658 connection.close()
13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	659
927 c030f735c594 added assignment of trajectories to prototypes and cleanup of insert queries Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 922 diff changeset	660 def loadPOIsFromSqlite(filename):
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	661 'Loads all 2D Gaussians in the database'
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	662 from sklearn import mixture # todo if not avalaible, load data in duck-typed class with same fields
872 c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	663 from ast import literal_eval
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	664 connection = sqlite3.connect(filename)
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	665 cursor = connection.cursor()
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	666 pois = []
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	667 try:
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	668 cursor.execute('SELECT * from gaussians2d')
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	669 gmmId = None
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	670 gmm = []
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	671 for row in cursor:
872 c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	672 if gmmId is None or row[0] != gmmId:
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	673 if len(gmm) > 0:
871 6db83beb5350 work in progress to update gaussian mixtures Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 862 diff changeset	674 tmp = mixture.GaussianMixture(len(gmm), covarianceType)
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	675 tmp.means_ = array([gaussian['mean'] for gaussian in gmm])
871 6db83beb5350 work in progress to update gaussian mixtures Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 862 diff changeset	676 tmp.covariances_ = array([gaussian['covar'] for gaussian in gmm])
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	677 tmp.weights_ = array([gaussian['weight'] for gaussian in gmm])
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	678 tmp.gmmTypes = [gaussian['type'] for gaussian in gmm]
872 c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	679 tmp.precisions_cholesky_ = array([gaussian['precisions'] for gaussian in gmm])
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	680 pois.append(tmp)
872 c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	681 gaussian = {'type': row[2],
c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	682 'mean': row[3:5],
c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	683 'covar': array(literal_eval(row[5])),
c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	684 'weight': row[7],
c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	685 'precisions': array(literal_eval(row[8]))}
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	686 gmm = [gaussian]
872 c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	687 covarianceType = row[6]
c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	688 gmmId = row[0]
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	689 else:
872 c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	690 gmm.append({'type': row[2],
c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	691 'mean': row[3:5],
c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	692 'covar': array(literal_eval(row[5])),
c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	693 'weight': row[7],
c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	694 'precisions': array(literal_eval(row[8]))})
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	695 if len(gmm) > 0:
871 6db83beb5350 work in progress to update gaussian mixtures Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 862 diff changeset	696 tmp = mixture.GaussianMixture(len(gmm), covarianceType)
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	697 tmp.means_ = array([gaussian['mean'] for gaussian in gmm])
871 6db83beb5350 work in progress to update gaussian mixtures Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 862 diff changeset	698 tmp.covariances_ = array([gaussian['covar'] for gaussian in gmm])
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	699 tmp.weights_ = array([gaussian['weight'] for gaussian in gmm])
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	700 tmp.gmmTypes = [gaussian['type'] for gaussian in gmm]
872 c70adaeeddf5 solved issue with latest version of scikit-learn Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 871 diff changeset	701 tmp.precisions_cholesky_ = array([gaussian['precisions'] for gaussian in gmm])
805 180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	702 pois.append(tmp)
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	703 except sqlite3.OperationalError as error:
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	704 printDBError(error)
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	705 connection.close()
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	706 return pois
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	707
180b6b0231c0 added saving/loading points of interests Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 784 diff changeset	708 #########################
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	709 # saving and loading for scene interpretation (Mohamed Gomaa Mohamed's PhD)
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	710 #########################
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	711
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	712 def writePrototypesToSqlite(prototypes,nMatching, outputFilename):
915 13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	713 ''' prototype dataset is a dictionary with keys== routes, values== prototypes Ids '''
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	714 connection = sqlite3.connect(outputFilename)
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	715 cursor = connection.cursor()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	716
915 13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	717 cursor.execute('CREATE TABLE IF NOT EXISTS prototypes (prototype_id INTEGER,routeIDstart INTEGER,routeIDend INTEGER, nMatching INTEGER, PRIMARY KEY(prototype_id))')
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	718
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	719 for route in prototypes.keys():
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	720 if prototypes[route]!=[]:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	721 for i in prototypes[route]:
915 13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	722 cursor.execute('insert into prototypes (prototype_id, routeIDstart,routeIDend, nMatching) values (?,?,?,?)',(i,route[0],route[1],nMatching[route][i]))
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	723
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	724 connection.commit()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	725 connection.close()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	726
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	727 def readPrototypesFromSqlite(filename):
915 13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	728 '''
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	729 This function loads the prototype file in the database
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	730 It returns a dictionary for prototypes for each route and nMatching
915 13434f5017dd work to save trajectory assignment to origin and destinations Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 913 diff changeset	731 '''
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	732 prototypes = {}
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	733 nMatching={}
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	734
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	735 connection = sqlite3.connect(filename)
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	736 cursor = connection.cursor()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	737
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	738 try:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	739 cursor.execute('SELECT * from prototypes order by prototype_id, routeIDstart,routeIDend, nMatching')
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	740 except sqlite3.OperationalError as error:
773 bf4a1790cfac minor bug and improvements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 768 diff changeset	741 printDBError(error)
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	742 return []
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	743
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	744 for row in cursor:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	745 route=(row[1],row[2])
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	746 if route not in prototypes.keys():
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	747 prototypes[route]=[]
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	748 prototypes[route].append(row[0])
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	749 nMatching[row[0]]=row[3]
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	750
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	751 connection.close()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	752 return prototypes,nMatching
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	753
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	754 def writeLabelsToSqlite(labels, outputFilename):
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	755 """ labels is a dictionary with keys: routes, values: prototypes Ids
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	756 """
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	757 connection = sqlite3.connect(outputFilename)
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	758 cursor = connection.cursor()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	759
777 ef6dd60be2e1 added function to save feature trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 773 diff changeset	760 cursor.execute("CREATE TABLE IF NOT EXISTS labels (object_id INTEGER,routeIDstart INTEGER,routeIDend INTEGER, prototype_id INTEGER, PRIMARY KEY(object_id))")
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	761
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	762 for route in labels.keys():
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	763 if labels[route]!=[]:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	764 for i in labels[route]:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	765 for j in labels[route][i]:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	766 cursor.execute("insert into labels (object_id, routeIDstart,routeIDend, prototype_id) values (?,?,?,?)",(j,route[0],route[1],i))
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	767
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	768 connection.commit()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	769 connection.close()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	770
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	771 def loadLabelsFromSqlite(filename):
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	772 labels = {}
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	773
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	774 connection = sqlite3.connect(filename)
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	775 cursor = connection.cursor()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	776
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	777 try:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	778 cursor.execute('SELECT * from labels order by object_id, routeIDstart,routeIDend, prototype_id')
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	779 except sqlite3.OperationalError as error:
773 bf4a1790cfac minor bug and improvements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 768 diff changeset	780 printDBError(error)
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	781 return []
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	782
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	783 for row in cursor:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	784 route=(row[1],row[2])
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	785 p=row[3]
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	786 if route not in labels.keys():
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	787 labels[route]={}
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	788 if p not in labels[route].keys():
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	789 labels[route][p]=[]
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	790 labels[route][p].append(row[0])
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	791
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	792 connection.close()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	793 return labels
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	794
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	795 def writeSpeedPrototypeToSqlite(prototypes,nmatching, outFilename):
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	796 """ to match the format of second layer prototypes"""
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	797 connection = sqlite3.connect(outFilename)
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	798 cursor = connection.cursor()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	799
777 ef6dd60be2e1 added function to save feature trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 773 diff changeset	800 cursor.execute("CREATE TABLE IF NOT EXISTS speedprototypes (spdprototype_id INTEGER,prototype_id INTEGER,routeID_start INTEGER, routeID_end INTEGER, nMatching INTEGER, PRIMARY KEY(spdprototype_id))")
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	801
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	802 for route in prototypes.keys():
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	803 if prototypes[route]!={}:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	804 for i in prototypes[route]:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	805 if prototypes[route][i]!= []:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	806 for j in prototypes[route][i]:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	807 cursor.execute("insert into speedprototypes (spdprototype_id,prototype_id, routeID_start, routeID_end, nMatching) values (?,?,?,?,?)",(j,i,route[0],route[1],nmatching[j]))
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	808
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	809 connection.commit()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	810 connection.close()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	811
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	812 def loadSpeedPrototypeFromSqlite(filename):
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	813 """
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	814 This function loads the prototypes table in the database of name <filename>.
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	815 """
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	816 prototypes = {}
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	817 nMatching={}
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	818 connection = sqlite3.connect(filename)
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	819 cursor = connection.cursor()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	820
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	821 try:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	822 cursor.execute('SELECT * from speedprototypes order by spdprototype_id,prototype_id, routeID_start, routeID_end, nMatching')
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	823 except sqlite3.OperationalError as error:
773 bf4a1790cfac minor bug and improvements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 768 diff changeset	824 printDBError(error)
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	825 return []
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	826
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	827 for row in cursor:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	828 route=(row[2],row[3])
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	829 if route not in prototypes.keys():
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	830 prototypes[route]={}
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	831 if row[1] not in prototypes[route].keys():
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	832 prototypes[route][row[1]]=[]
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	833 prototypes[route][row[1]].append(row[0])
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	834 nMatching[row[0]]=row[4]
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	835
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	836 connection.close()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	837 return prototypes,nMatching
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	838
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	839
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	840 def writeRoutesToSqlite(Routes, outputFilename):
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	841 """ This function writes the activity path define by start and end IDs"""
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	842 connection = sqlite3.connect(outputFilename)
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	843 cursor = connection.cursor()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	844
777 ef6dd60be2e1 added function to save feature trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 773 diff changeset	845 cursor.execute("CREATE TABLE IF NOT EXISTS routes (object_id INTEGER,routeIDstart INTEGER,routeIDend INTEGER, PRIMARY KEY(object_id))")
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	846
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	847 for route in Routes.keys():
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	848 if Routes[route]!=[]:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	849 for i in Routes[route]:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	850 cursor.execute("insert into routes (object_id, routeIDstart,routeIDend) values (?,?,?)",(i,route[0],route[1]))
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	851
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	852 connection.commit()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	853 connection.close()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	854
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	855 def loadRoutesFromSqlite(filename):
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	856 Routes = {}
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	857
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	858 connection = sqlite3.connect(filename)
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	859 cursor = connection.cursor()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	860
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	861 try:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	862 cursor.execute('SELECT * from routes order by object_id, routeIDstart,routeIDend')
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	863 except sqlite3.OperationalError as error:
773 bf4a1790cfac minor bug and improvements Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 768 diff changeset	864 printDBError(error)
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	865 return []
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	866
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	867 for row in cursor:
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	868 route=(row[1],row[2])
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	869 if route not in Routes.keys():
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	870 Routes[route]=[]
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	871 Routes[route].append(row[0])
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	872
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	873 connection.close()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	874 return Routes
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	875
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	876 def setRoutes(filename, objects):
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	877 connection = sqlite3.connect(filename)
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	878 cursor = connection.cursor()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	879 for obj in objects:
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	880 cursor.execute('update objects set startRouteID = {} WHERE object_id = {}'.format(obj.startRouteID, obj.getNum()))
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	881 cursor.execute('update objects set endRouteID = {} WHERE object_id = {}'.format(obj.endRouteID, obj.getNum()))
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	882 connection.commit()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	883 connection.close()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	884
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	885 def setRoadUserTypes(filename, objects):
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	886 '''Saves the user types of the objects in the sqlite database stored in filename
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	887 The objects should exist in the objects table'''
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	888 connection = sqlite3.connect(filename)
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	889 cursor = connection.cursor()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	890 for obj in objects:
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	891 cursor.execute('update objects set road_user_type = {} WHERE object_id = {}'.format(obj.getUserType(), obj.getNum()))
736 967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	892 connection.commit()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	893 connection.close()
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	894
967d244968a4 work in progress on saving/loading prototypes Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 722 diff changeset	895 #########################
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	896 # txt files
a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	897 #########################
a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	898
525 7124c7d2a663 first draft of loading from VISSIM file Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 524 diff changeset	899 def openCheck(filename, option = 'r', quitting = False):
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	900 '''Open file filename in read mode by default
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	901 and checks it is open'''
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	902 try:
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	903 return open(filename, option)
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	904 except IOError:
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	905 print 'File %s could not be opened.' % filename
525 7124c7d2a663 first draft of loading from VISSIM file Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 524 diff changeset	906 if quitting:
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	907 from sys import exit
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	908 exit()
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	909 return None
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	910
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	911 def readline(f, commentCharacters = commentChar):
36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	912 '''Modified readline function to skip comments
36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	913 Can take a list of characters or a string (in will work in both)'''
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	914 s = f.readline()
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	915 while (len(s) > 0) and s[0] in commentCharacters:
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	916 s = f.readline()
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	917 return s.strip()
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	918
722 49e99ca34a7d corrected bugs in old ubc code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 718 diff changeset	919 def getLines(f, delimiterChar = delimiterChar, commentCharacters = commentChar):
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	920 '''Gets a complete entry (all the lines) in between delimiterChar.'''
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	921 dataStrings = []
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	922 s = readline(f, commentCharacters)
722 49e99ca34a7d corrected bugs in old ubc code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 718 diff changeset	923 while len(s) > 0 and s[0] != delimiterChar:
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	924 dataStrings += [s.strip()]
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	925 s = readline(f, commentCharacters)
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	926 return dataStrings
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	927
829 0ddcc41663f5 renaming Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 813 diff changeset	928 def saveList(filename, l):
514 1ba618fb0f70 corrected bug from merging and argument issue in display-trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 509 diff changeset	929 f = openCheck(filename, 'w')
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	930 for x in l:
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	931 f.write('{}\n'.format(x))
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	932 f.close()
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	933
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	934 def loadListStrings(filename, commentCharacters = commentChar):
514 1ba618fb0f70 corrected bug from merging and argument issue in display-trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 509 diff changeset	935 f = openCheck(filename, 'r')
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	936 result = getLines(f, commentCharacters)
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	937 f.close()
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	938 return result
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	939
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	940 def getValuesFromINIFile(filename, option, delimiterChar = '=', commentCharacters = commentChar):
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	941 values = []
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	942 for l in loadListStrings(filename, commentCharacters):
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	943 if l.startswith(option):
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	944 values.append(l.split(delimiterChar)[1].strip())
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	945 return values
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	946
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	947 class FakeSecHead(object):
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	948 '''Add fake section header [asection]
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	949
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	950 from http://stackoverflow.com/questions/2819696/parsing-properties-file-in-python/2819788#2819788
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	951 use read_file in Python 3.2+
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	952 '''
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	953 def __init__(self, fp):
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	954 self.fp = fp
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	955 self.sechead = '[main]\n'
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	956
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	957 def readline(self):
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	958 if self.sechead:
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	959 try: return self.sechead
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	960 finally: self.sechead = None
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	961 else: return self.fp.readline()
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	962
892 f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	963 def loadPemsTraffic(filename):
f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	964 '''Loads traffic data downloaded from the http://pems.dot.ca.gov clearinghouse
f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	965 into pandas dataframe'''
f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	966 f = openCheck(filename)
f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	967 l = f.readline().strip()
f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	968 items = l.split(',')
f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	969 headers = ['time', 'station', 'district', 'route', 'direction', 'lanetype', 'length', 'nsamples', 'pctobserved', 'flow', 'occupancy', 'speed', 'delay35', 'delay40', 'delay45', 'delay50', 'delay55', 'delay60']
f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	970 nLanes = (len(items)-len(headers))/3
f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	971 for i in xrange(nLanes):
f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	972 headers += ['flow{}'.format(i+1), 'occupancy{}'.format(i+1), 'speed{}'.format(i+1)]
f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	973 f.close()
f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	974 return read_csv(filename, delimiter = ',', names = headers)
f766fe0995f4 added function for PeMS data Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 882 diff changeset	975
649 df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	976 def generatePDLaneColumn(data):
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	977 data['LANE'] = data['LANE\LINK\NO'].astype(str)+'_'+data['LANE\INDEX'].astype(str)
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	978
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	979 def convertTrajectoriesVissimToSqlite(filename):
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	980 '''Relies on a system call to sqlite3
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	981 sqlite3 [file.sqlite] < import_fzp.sql'''
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	982 sqlScriptFilename = "import_fzp.sql"
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	983 # create sql file
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	984 out = openCheck(sqlScriptFilename, "w")
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	985 out.write(".separator \";\"\n"+
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	986 "CREATE TABLE IF NOT EXISTS curvilinear_positions (t REAL, trajectory_id INTEGER, link_id INTEGER, lane_id INTEGER, s_coordinate REAL, y_coordinate REAL, speed REAL, PRIMARY KEY (t, trajectory_id));\n"+
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	987 ".import "+filename+" curvilinear_positions\n"+
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	988 "DELETE FROM curvilinear_positions WHERE trajectory_id IS NULL OR trajectory_id = \"NO\";\n")
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	989 out.close()
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	990 # system call
753 3d48e34db846 switched to subprocess.check_call Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 752 diff changeset	991 from subprocess import check_call
3d48e34db846 switched to subprocess.check_call Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 752 diff changeset	992 out = openCheck("err.log", "w")
3d48e34db846 switched to subprocess.check_call Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 752 diff changeset	993 check_call("sqlite3 "+utils.removeExtension(filename)+".sqlite < "+sqlScriptFilename, stderr = out, shell = True)
3d48e34db846 switched to subprocess.check_call Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 752 diff changeset	994 out.close()
3d48e34db846 switched to subprocess.check_call Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 752 diff changeset	995 shutil.os.remove(sqlScriptFilename)
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	996
754 782e8fd3672c added function to find object ids going through some vissim links Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 753 diff changeset	997 def loadObjectNumbersInLinkFromVissimFile(filename, linkIds):
782e8fd3672c added function to find object ids going through some vissim links Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 753 diff changeset	998 '''Finds the ids of the objects that go through any of the link in the list linkIds'''
782e8fd3672c added function to find object ids going through some vissim links Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 753 diff changeset	999 connection = sqlite3.connect(filename)
782e8fd3672c added function to find object ids going through some vissim links Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 753 diff changeset	1000 cursor = connection.cursor()
782e8fd3672c added function to find object ids going through some vissim links Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 753 diff changeset	1001 queryStatement = 'SELECT DISTINCT trajectory_id FROM curvilinear_positions where link_id IN ('+','.join([str(id) for id in linkIds])+')'
782e8fd3672c added function to find object ids going through some vissim links Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 753 diff changeset	1002 try:
782e8fd3672c added function to find object ids going through some vissim links Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 753 diff changeset	1003 cursor.execute(queryStatement)
782e8fd3672c added function to find object ids going through some vissim links Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 753 diff changeset	1004 return [row[0] for row in cursor]
782e8fd3672c added function to find object ids going through some vissim links Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 753 diff changeset	1005 except sqlite3.OperationalError as error:
782e8fd3672c added function to find object ids going through some vissim links Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 753 diff changeset	1006 printDBError(error)
782e8fd3672c added function to find object ids going through some vissim links Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 753 diff changeset	1007
759 a05b70f307dd added function to count vehicles per VISSIM link Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 756 diff changeset	1008 def getNObjectsInLinkFromVissimFile(filename, linkIds):
a05b70f307dd added function to count vehicles per VISSIM link Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 756 diff changeset	1009 '''Returns the number of objects that traveled through the link ids'''
a05b70f307dd added function to count vehicles per VISSIM link Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 756 diff changeset	1010 connection = sqlite3.connect(filename)
a05b70f307dd added function to count vehicles per VISSIM link Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 756 diff changeset	1011 cursor = connection.cursor()
a05b70f307dd added function to count vehicles per VISSIM link Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 756 diff changeset	1012 queryStatement = 'SELECT link_id, COUNT(DISTINCT trajectory_id) FROM curvilinear_positions where link_id IN ('+','.join([str(id) for id in linkIds])+') GROUP BY link_id'
a05b70f307dd added function to count vehicles per VISSIM link Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 756 diff changeset	1013 try:
a05b70f307dd added function to count vehicles per VISSIM link Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 756 diff changeset	1014 cursor.execute(queryStatement)
a05b70f307dd added function to count vehicles per VISSIM link Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 756 diff changeset	1015 return {row[0]:row[1] for row in cursor}
a05b70f307dd added function to count vehicles per VISSIM link Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 756 diff changeset	1016 except sqlite3.OperationalError as error:
a05b70f307dd added function to count vehicles per VISSIM link Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 756 diff changeset	1017 printDBError(error)
754 782e8fd3672c added function to find object ids going through some vissim links Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 753 diff changeset	1018
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1019 def loadTrajectoriesFromVissimFile(filename, simulationStepsPerTimeUnit, objectNumbers = None, warmUpLastInstant = None, usePandas = False, nDecimals = 2, lowMemory = True):
525 7124c7d2a663 first draft of loading from VISSIM file Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 524 diff changeset	1020 '''Reads data from VISSIM .fzp trajectory file
755 f3aeb0b47eff comment improvement: vissim time is multiplied to get integers similar to frame numbers Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 754 diff changeset	1021 simulationStepsPerTimeUnit is the number of simulation steps per unit of time used by VISSIM (second)
f3aeb0b47eff comment improvement: vissim time is multiplied to get integers similar to frame numbers Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 754 diff changeset	1022 for example, there seems to be 10 simulation steps per simulated second in VISSIM,
f3aeb0b47eff comment improvement: vissim time is multiplied to get integers similar to frame numbers Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 754 diff changeset	1023 so simulationStepsPerTimeUnit should be 10,
527 37830a831818 loading from VISSIM trajectory data works Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 526 diff changeset	1024 so that all times correspond to the number of the simulation step (and can be stored as integers)
660 994dd644f6ab corrected impact of warmUpLastInstant Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 649 diff changeset	1025
994dd644f6ab corrected impact of warmUpLastInstant Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 649 diff changeset	1026 Objects positions will be considered only after warmUpLastInstant
994dd644f6ab corrected impact of warmUpLastInstant Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 649 diff changeset	1027 (if the object has no such position, it won't be loaded)
525 7124c7d2a663 first draft of loading from VISSIM file Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 524 diff changeset	1028
752 14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1029 Assumed to be sorted over time
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1030 Warning: if reading from SQLite a limited number of objects, objectNumbers will be the maximum object id'''
525 7124c7d2a663 first draft of loading from VISSIM file Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 524 diff changeset	1031 objects = {} # dictionary of objects index by their id
7124c7d2a663 first draft of loading from VISSIM file Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 524 diff changeset	1032
636 3058e00887bc removed all issues because of tests with None, using is instead of == or != Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 628 diff changeset	1033 if usePandas:
718 2cd245cb780d added option to set low_memory = False for pandas.read_csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 714 diff changeset	1034 data = read_csv(filename, delimiter=';', comment='*', header=0, skiprows = 1, low_memory = lowMemory)
649 df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1035 generatePDLaneColumn(data)
642 932f96c89212 added pandas to read vissim fzp (more robust with respect to column names Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 636 diff changeset	1036 data['TIME'] = data['$VEHICLE:SIMSEC']*simulationStepsPerTimeUnit
660 994dd644f6ab corrected impact of warmUpLastInstant Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 649 diff changeset	1037 if warmUpLastInstant is not None:
994dd644f6ab corrected impact of warmUpLastInstant Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 649 diff changeset	1038 data = data[data['TIME']>=warmUpLastInstant]
642 932f96c89212 added pandas to read vissim fzp (more robust with respect to column names Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 636 diff changeset	1039 grouped = data.loc[:,['NO','TIME']].groupby(['NO'], as_index = False)
704 f83d125d0c55 cleaning of storage.py and addition of type conversion for VISSIM files (from Laurent Gauthier) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 702 diff changeset	1040 instants = grouped['TIME'].agg({'first': npmin, 'last': npmax})
642 932f96c89212 added pandas to read vissim fzp (more robust with respect to column names Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 636 diff changeset	1041 for row_index, row in instants.iterrows():
932f96c89212 added pandas to read vissim fzp (more robust with respect to column names Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 636 diff changeset	1042 objNum = int(row['NO'])
932f96c89212 added pandas to read vissim fzp (more robust with respect to column names Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 636 diff changeset	1043 tmp = data[data['NO'] == objNum]
932f96c89212 added pandas to read vissim fzp (more robust with respect to column names Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 636 diff changeset	1044 objects[objNum] = moving.MovingObject(num = objNum, timeInterval = moving.TimeInterval(row['first'], row['last']))
932f96c89212 added pandas to read vissim fzp (more robust with respect to column names Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 636 diff changeset	1045 # positions should be rounded to nDecimals decimals only
704 f83d125d0c55 cleaning of storage.py and addition of type conversion for VISSIM files (from Laurent Gauthier) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 702 diff changeset	1046 objects[objNum].curvilinearPositions = moving.CurvilinearTrajectory(S = npround(tmp['POS'].tolist(), nDecimals), Y = npround(tmp['POSLAT'].tolist(), nDecimals), lanes = tmp['LANE'].tolist())
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1047 if objectNumbers is not None and objectNumbers > 0 and len(objects) >= objectNumbers:
752 14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1048 objects.values()
636 3058e00887bc removed all issues because of tests with None, using is instead of == or != Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 628 diff changeset	1049 else:
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1050 if filename.endswith(".fzp"):
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1051 inputfile = openCheck(filename, quitting = True)
636 3058e00887bc removed all issues because of tests with None, using is instead of == or != Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 628 diff changeset	1052 line = readline(inputfile, '*$')
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1053 while len(line) > 0:#for line in inputfile:
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1054 data = line.strip().split(';')
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1055 objNum = int(data[1])
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1056 instant = float(data[0])*simulationStepsPerTimeUnit
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1057 s = float(data[4])
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1058 y = float(data[5])
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1059 lane = data[2]+'_'+data[3]
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1060 if objNum not in objects:
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1061 if warmUpLastInstant is None or instant >= warmUpLastInstant:
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1062 if objectNumbers is None or len(objects) < objectNumbers:
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1063 objects[objNum] = moving.MovingObject(num = objNum, timeInterval = moving.TimeInterval(instant, instant))
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1064 objects[objNum].curvilinearPositions = moving.CurvilinearTrajectory()
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1065 if (warmUpLastInstant is None or instant >= warmUpLastInstant) and objNum in objects:
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1066 objects[objNum].timeInterval.last = instant
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1067 objects[objNum].curvilinearPositions.addPositionSYL(s, y, lane)
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1068 line = readline(inputfile, '*$')
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1069 elif filename.endswith(".sqlite"):
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1070 connection = sqlite3.connect(filename)
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1071 cursor = connection.cursor()
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1072 queryStatement = 'SELECT t, trajectory_id, link_id, lane_id, s_coordinate, y_coordinate FROM curvilinear_positions'
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1073 if objectNumbers is not None:
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1074 queryStatement += ' WHERE trajectory_id '+getObjectCriteria(objectNumbers)
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1075 queryStatement += ' ORDER BY trajectory_id, t'
752 14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1076 try:
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1077 cursor.execute(queryStatement)
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1078 for row in cursor:
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1079 objNum = row[1]
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1080 instant = row[0]*simulationStepsPerTimeUnit
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1081 s = row[4]
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1082 y = row[5]
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1083 lane = '{}_{}'.format(row[2], row[3])
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1084 if objNum not in objects:
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1085 if warmUpLastInstant is None or instant >= warmUpLastInstant:
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1086 if objectNumbers is None or len(objects) < objectNumbers:
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1087 objects[objNum] = moving.MovingObject(num = objNum, timeInterval = moving.TimeInterval(instant, instant))
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1088 objects[objNum].curvilinearPositions = moving.CurvilinearTrajectory()
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1089 if (warmUpLastInstant is None or instant >= warmUpLastInstant) and objNum in objects:
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1090 objects[objNum].timeInterval.last = instant
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1091 objects[objNum].curvilinearPositions.addPositionSYL(s, y, lane)
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1092 except sqlite3.OperationalError as error:
14963a9c3b09 debug finished Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 751 diff changeset	1093 printDBError(error)
750 6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1094 else:
6049e9b6902c work in progress storage vissim sqlite Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 744 diff changeset	1095 print("File type of "+filename+" not supported (only .sqlite and .fzp files)")
751 79405a938407 corrected bug Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 750 diff changeset	1096 return objects.values()
645 5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1097
649 df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1098 def selectPDLanes(data, lanes = None):
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1099 '''Selects the subset of data for the right lanes
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1100
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1101 Lane format is a string 'x_y' where x is link index and y is lane index'''
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1102 if lanes is not None:
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1103 if 'LANE' not in data.columns:
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1104 generatePDLaneColumn(data)
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1105 indices = (data['LANE'] == lanes[0])
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1106 for l in lanes[1:]:
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1107 indices = indices \| (data['LANE'] == l)
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1108 return data[indices]
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1109 else:
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1110 return data
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1111
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1112 def countStoppedVehiclesVissim(filename, lanes = None, proportionStationaryTime = 0.7):
646 6680a69d5bf3 added fast function to detect VISSIM simulations with errors Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 645 diff changeset	1113 '''Counts the number of vehicles stopped for a long time in a VISSIM trajectory file
6680a69d5bf3 added fast function to detect VISSIM simulations with errors Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 645 diff changeset	1114 and the total number of vehicles
6680a69d5bf3 added fast function to detect VISSIM simulations with errors Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 645 diff changeset	1115
6680a69d5bf3 added fast function to detect VISSIM simulations with errors Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 645 diff changeset	1116 Vehicles are considered finally stationary
649 df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1117 if more than proportionStationaryTime of their total time
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1118 If lanes is not None, only the data for the selected lanes will be provided
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1119 (format as string x_y where x is link index and y is lane index)'''
756 a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1120 if filename.endswith(".fzp"):
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1121 columns = ['NO', '$VEHICLE:SIMSEC', 'POS']
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1122 if lanes is not None:
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1123 columns += ['LANE\LINK\NO', 'LANE\INDEX']
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1124 data = read_csv(filename, delimiter=';', comment='*', header=0, skiprows = 1, usecols = columns, low_memory = lowMemory)
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1125 data = selectPDLanes(data, lanes)
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1126 data.sort(['$VEHICLE:SIMSEC'], inplace = True)
646 6680a69d5bf3 added fast function to detect VISSIM simulations with errors Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 645 diff changeset	1127
756 a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1128 nStationary = 0
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1129 nVehicles = 0
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1130 for name, group in data.groupby(['NO'], sort = False):
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1131 nVehicles += 1
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1132 positions = array(group['POS'])
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1133 diff = positions[1:]-positions[:-1]
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1134 if npsum(diff == 0.) >= proportionStationaryTime*(len(positions)-1):
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1135 nStationary += 1
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1136 elif filename.endswith(".sqlite"):
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1137 # select trajectory_id, t, s_coordinate, speed from curvilinear_positions where trajectory_id between 1860 and 1870 and speed < 0.1
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1138 # pb of the meaning of proportionStationaryTime in arterial network? Why proportion of existence time?
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1139 pass
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1140 else:
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1141 print("File type of "+filename+" not supported (only .sqlite and .fzp files)")
646 6680a69d5bf3 added fast function to detect VISSIM simulations with errors Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 645 diff changeset	1142
6680a69d5bf3 added fast function to detect VISSIM simulations with errors Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 645 diff changeset	1143 return nStationary, nVehicles
6680a69d5bf3 added fast function to detect VISSIM simulations with errors Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 645 diff changeset	1144
744 ed6ff2ec0aeb bug correction from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 740 diff changeset	1145 def countCollisionsVissim(filename, lanes = None, collisionTimeDifference = 0.2, lowMemory = True):
645 5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1146 '''Counts the number of collisions per lane in a VISSIM trajectory file
5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1147
5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1148 To distinguish between cars passing and collision,
5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1149 one checks when the sign of the position difference inverts
649 df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1150 (if the time are closer than collisionTimeDifference)
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1151 If lanes is not None, only the data for the selected lanes will be provided
df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1152 (format as string x_y where x is link index and y is lane index)'''
718 2cd245cb780d added option to set low_memory = False for pandas.read_csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 714 diff changeset	1153 data = read_csv(filename, delimiter=';', comment='*', header=0, skiprows = 1, usecols = ['LANE\LINK\NO', 'LANE\INDEX', '$VEHICLE:SIMSEC', 'NO', 'POS'], low_memory = lowMemory)
649 df9ddeaee4a6 added ability to select lanes to count collisions or stationary vehicles in VISSIM files Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 647 diff changeset	1154 data = selectPDLanes(data, lanes)
744 ed6ff2ec0aeb bug correction from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 740 diff changeset	1155 data = data.convert_objects(convert_numeric=True)
ed6ff2ec0aeb bug correction from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 740 diff changeset	1156
645 5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1157 merged = merge(data, data, how='inner', left_on=['LANE\LINK\NO', 'LANE\INDEX', '$VEHICLE:SIMSEC'], right_on=['LANE\LINK\NO', 'LANE\INDEX', '$VEHICLE:SIMSEC'], sort = False)
5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1158 merged = merged[merged['NO_x']>merged['NO_y']]
5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1159
5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1160 nCollisions = 0
5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1161 for name, group in merged.groupby(['LANE\LINK\NO', 'LANE\INDEX', 'NO_x', 'NO_y']):
744 ed6ff2ec0aeb bug correction from Laurent Gauthier Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 740 diff changeset	1162 diff = group['POS_x']-group['POS_y']
704 f83d125d0c55 cleaning of storage.py and addition of type conversion for VISSIM files (from Laurent Gauthier) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 702 diff changeset	1163 # diff = group['POS_x']-group['POS_y'] # to check the impact of convert_objects and the possibility of using type conversion in read_csv or function to convert strings if any
f83d125d0c55 cleaning of storage.py and addition of type conversion for VISSIM files (from Laurent Gauthier) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 702 diff changeset	1164 if len(diff) >= 2 and npmin(diff) < 0 and npmax(diff) > 0:
645 5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1165 xidx = diff[diff < 0].argmax()
5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1166 yidx = diff[diff > 0].argmin()
5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1167 if abs(group.loc[xidx, '$VEHICLE:SIMSEC'] - group.loc[yidx, '$VEHICLE:SIMSEC']) <= collisionTimeDifference:
5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1168 nCollisions += 1
756 a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1169
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1170 # select TD1.link_id, TD1.lane_id from temp.diff_positions as TD1, temp.diff_positions as TD2 where TD1.link_id = TD2.link_id and TD1.lane_id = TD2.lane_id and TD1.id1 = TD2.id1 and TD1.id2 = TD2.id2 and TD1.t = TD2.t+0.1 and TD1.diff*TD2.diff < 0; # besoin de faire un group by??
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1171 # create temp table diff_positions as select CP1.t as t, CP1.link_id as link_id, CP1.lane_id as lane_id, CP1.trajectory_id as id1, CP2.trajectory_id as id2, CP1.s_coordinate - CP2.s_coordinate as diff from curvilinear_positions CP1, curvilinear_positions CP2 where CP1.link_id = CP2.link_id and CP1.lane_id = CP2.lane_id and CP1.t = CP2.t and CP1.trajectory_id > CP2.trajectory_id;
a73f43aac00e moved pandas imports, sql in comments Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 755 diff changeset	1172 # SQL select link_id, lane_id, id1, id2, min(diff), max(diff) from (select CP1.t as t, CP1.link_id as link_id, CP1.lane_id as lane_id, CP1.trajectory_id as id1, CP2.trajectory_id as id2, CP1.s_coordinate - CP2.s_coordinate as diff from curvilinear_positions CP1, curvilinear_positions CP2 where CP1.link_id = CP2.link_id and CP1.lane_id = CP2.lane_id and CP1.t = CP2.t and CP1.trajectory_id > CP2.trajectory_id) group by link_id, lane_id, id1, id2 having min(diff)*max(diff) < 0
645 5ed2118c959d added function to count collisions in vissim Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 643 diff changeset	1173 return nCollisions
524 1dced8932b08 corrected bugs Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 514 diff changeset	1174
173 319a04ba9033 function name change Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 143 diff changeset	1175 def loadTrajectoriesFromNgsimFile(filename, nObjects = -1, sequenceNum = -1):
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1176 '''Reads data from the trajectory data provided by NGSIM project
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1177 and returns the list of Feature objects'''
7 ffddccfab7f9 loading shell objects from NGSIM works Nicolas Saunier <nico@confins.net> parents: 0 diff changeset	1178 objects = []
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1179
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	1180 inputfile = openCheck(filename, quitting = True)
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1181
7 ffddccfab7f9 loading shell objects from NGSIM works Nicolas Saunier <nico@confins.net> parents: 0 diff changeset	1182 def createObject(numbers):
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1183 firstFrameNum = int(numbers[1])
7 ffddccfab7f9 loading shell objects from NGSIM works Nicolas Saunier <nico@confins.net> parents: 0 diff changeset	1184 # do the geometry and usertype
ffddccfab7f9 loading shell objects from NGSIM works Nicolas Saunier <nico@confins.net> parents: 0 diff changeset	1185
72 575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1186 firstFrameNum = int(numbers[1])
575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1187 lastFrameNum = firstFrameNum+int(numbers[2])-1
575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1188 #time = moving.TimeInterval(firstFrameNum, firstFrameNum+int(numbers[2])-1)
575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1189 obj = moving.MovingObject(num = int(numbers[0]),
575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1190 timeInterval = moving.TimeInterval(firstFrameNum, lastFrameNum),
575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1191 positions = moving.Trajectory([[float(numbers[6])],[float(numbers[7])]]),
575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1192 userType = int(numbers[10]))
78 99e807c29753 added loading other information from NGSIM Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 72 diff changeset	1193 obj.userType = int(numbers[10])
99e807c29753 added loading other information from NGSIM Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 72 diff changeset	1194 obj.laneNums = [int(numbers[13])]
99e807c29753 added loading other information from NGSIM Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 72 diff changeset	1195 obj.precedingVehicles = [int(numbers[14])] # lead vehicle (before)
99e807c29753 added loading other information from NGSIM Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 72 diff changeset	1196 obj.followingVehicles = [int(numbers[15])] # following vehicle (after)
99e807c29753 added loading other information from NGSIM Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 72 diff changeset	1197 obj.spaceHeadways = [float(numbers[16])] # feet
99e807c29753 added loading other information from NGSIM Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 72 diff changeset	1198 obj.timeHeadways = [float(numbers[17])] # seconds
327 42f2b46ec210 added class for trajectories in curvilinear coordinates Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 326 diff changeset	1199 obj.curvilinearPositions = moving.CurvilinearTrajectory([float(numbers[5])],[float(numbers[4])], obj.laneNums) # X is the longitudinal coordinate
72 575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1200 obj.speeds = [float(numbers[11])]
575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1201 obj.size = [float(numbers[8]), float(numbers[9])] # 8 lengh, 9 width # TODO: temporary, should use a geometry object
7 ffddccfab7f9 loading shell objects from NGSIM works Nicolas Saunier <nico@confins.net> parents: 0 diff changeset	1202 return obj
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1203
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	1204 numbers = readline(inputfile).strip().split()
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1205 if (len(numbers) > 0):
7 ffddccfab7f9 loading shell objects from NGSIM works Nicolas Saunier <nico@confins.net> parents: 0 diff changeset	1206 obj = createObject(numbers)
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1207
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	1208 for line in inputfile:
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1209 numbers = line.strip().split()
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	1210 if obj.getNum() != int(numbers[0]):
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1211 # check and adapt the length to deal with issues in NGSIM data
72 575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1212 if (obj.length() != obj.positions.length()):
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	1213 print 'length pb with object %s (%d,%d)' % (obj.getNum(),obj.length(),obj.positions.length())
72 575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1214 obj.last = obj.getFirstInstant()+obj.positions.length()-1
575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1215 #obj.velocities = utils.computeVelocities(f.positions) # compare norm to speeds ?
7 ffddccfab7f9 loading shell objects from NGSIM works Nicolas Saunier <nico@confins.net> parents: 0 diff changeset	1216 objects.append(obj)
ffddccfab7f9 loading shell objects from NGSIM works Nicolas Saunier <nico@confins.net> parents: 0 diff changeset	1217 if (nObjects>0) and (len(objects)>=nObjects):
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1218 break
7 ffddccfab7f9 loading shell objects from NGSIM works Nicolas Saunier <nico@confins.net> parents: 0 diff changeset	1219 obj = createObject(numbers)
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1220 else:
327 42f2b46ec210 added class for trajectories in curvilinear coordinates Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 326 diff changeset	1221 obj.laneNums.append(int(numbers[13]))
72 575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1222 obj.positions.addPositionXY(float(numbers[6]), float(numbers[7]))
542 a3add9f751ef added differentiate function for curvilinear trajectories and modified the addPosition functions Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 541 diff changeset	1223 obj.curvilinearPositions.addPositionSYL(float(numbers[5]), float(numbers[4]), obj.laneNums[-1])
72 575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1224 obj.speeds.append(float(numbers[11]))
78 99e807c29753 added loading other information from NGSIM Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 72 diff changeset	1225 obj.precedingVehicles.append(int(numbers[14]))
99e807c29753 added loading other information from NGSIM Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 72 diff changeset	1226 obj.followingVehicles.append(int(numbers[15]))
99e807c29753 added loading other information from NGSIM Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 72 diff changeset	1227 obj.spaceHeadways.append(float(numbers[16]))
99e807c29753 added loading other information from NGSIM Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 72 diff changeset	1228 obj.timeHeadways.append(float(numbers[17]))
7 ffddccfab7f9 loading shell objects from NGSIM works Nicolas Saunier <nico@confins.net> parents: 0 diff changeset	1229
72 575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1230 if (obj.size[0] != float(numbers[8])):
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	1231 print 'changed length obj %d' % (obj.getNum())
72 575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1232 if (obj.size[1] != float(numbers[9])):
329 a70c205ebdd9 added sqlite code, in particular to load and save road user type Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 327 diff changeset	1233 print 'changed width obj %d' % (obj.getNum())
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1234
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	1235 inputfile.close()
7 ffddccfab7f9 loading shell objects from NGSIM works Nicolas Saunier <nico@confins.net> parents: 0 diff changeset	1236 return objects
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1237
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	1238 def convertNgsimFile(inputfile, outputfile, append = False, nObjects = -1, sequenceNum = 0):
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1239 '''Reads data from the trajectory data provided by NGSIM project
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1240 and converts to our current format.'''
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1241 if append:
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	1242 out = openCheck(outputfile,'a')
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1243 else:
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	1244 out = openCheck(outputfile,'w')
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1245 nObjectsPerType = [0,0,0]
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1246
564 36605d843be5 modified bug for reading vissim files, cleaned use of readline with multiple type of characters for comments (to ignore) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 552 diff changeset	1247 features = loadNgsimFile(inputfile, sequenceNum)
0 aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1248 for f in features:
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1249 nObjectsPerType[f.userType-1] += 1
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1250 f.write(out)
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1251
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1252 print nObjectsPerType
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1253
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1254 out.close()
aed8eb63cdde initial commit with non-functional python code for NGSIM Nicolas Saunier <nico@confins.net> parents: diff changeset	1255
895 739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1256 def loadPinholeCameraModel(filename, tanalystFormat = True):
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1257 '''Loads the data from a file containing the camera parameters
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1258 (pinhole camera model, http://docs.opencv.org/2.4/modules/calib3d/doc/camera_calibration_and_3d_reconstruction.html)
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1259 and returns a dictionary'''
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1260 if tanalystFormat:
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1261 f = openCheck(filename, quitting = True)
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1262 content = getLines(f)
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1263 cameraData = {}
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1264 for l in content:
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1265 tmp = l.split(':')
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1266 cameraData[tmp[0]] = float(tmp[1].strip().replace(',','.'))
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1267 return cameraData
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1268 else:
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1269 print('Unknown camera model (not tanalyst format')
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1270 return None
739acd338cc0 added script to extract camera info from tacal file by Lund University (T analyst) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 893 diff changeset	1271
829 0ddcc41663f5 renaming Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 813 diff changeset	1272 def savePositionsToCsv(f, obj):
335 3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1273 timeInterval = obj.getTimeInterval()
3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1274 positions = obj.getPositions()
3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1275 curvilinearPositions = obj.getCurvilinearPositions()
3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1276 for i in xrange(int(obj.length())):
3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1277 p1 = positions[i]
3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1278 s = '{},{},{},{}'.format(obj.num,timeInterval[i],p1.x,p1.y)
636 3058e00887bc removed all issues because of tests with None, using is instead of == or != Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 628 diff changeset	1279 if curvilinearPositions is not None:
335 3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1280 p2 = curvilinearPositions[i]
3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1281 s += ',{},{}'.format(p2[0],p2[1])
3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1282 f.write(s+'\n')
3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1283
829 0ddcc41663f5 renaming Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 813 diff changeset	1284 def saveTrajectoriesToCsv(filename, objects):
514 1ba618fb0f70 corrected bug from merging and argument issue in display-trajectories Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 509 diff changeset	1285 f = openCheck(filename, 'w')
335 3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1286 for i,obj in enumerate(objects):
829 0ddcc41663f5 renaming Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 813 diff changeset	1287 savePositionsToCsv(f, obj)
335 3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1288 f.close()
3950bfe22768 added functions to export trajectories to csv Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 330 diff changeset	1289
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1290
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1291 #########################
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1292 # Utils to read .ini type text files for configuration, meta data...
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1293 #########################
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1294
812 21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1295 class ClassifierParameters(VideoFilenameAddable):
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1296 'Class for the parameters of object classifiers'
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1297 def loadConfigFile(self, filename):
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1298 from ConfigParser import ConfigParser
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1299
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1300 config = ConfigParser()
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1301 config.readfp(FakeSecHead(openCheck(filename)))
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1302 self.sectionHeader = config.sections()[0]
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1303
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1304 self.pedBikeCarSVMFilename = config.get(self.sectionHeader, 'pbv-svm-filename')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1305 self.bikeCarSVMFilename = config.get(self.sectionHeader, 'bv-svm-filename')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1306 self.percentIncreaseCrop = config.getfloat(self.sectionHeader, 'percent-increase-crop')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1307 self.minNPixels = config.getint(self.sectionHeader, 'min-npixels-crop')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1308 x = config.getint(self.sectionHeader, 'hog-rescale-size')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1309 self.hogRescaleSize = (x, x)
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1310 self.hogNOrientations = config.getint(self.sectionHeader, 'hog-norientations')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1311 x = config.getint(self.sectionHeader, 'hog-npixels-cell')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1312 self.hogNPixelsPerCell = (x, x)
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1313 x = config.getint(self.sectionHeader, 'hog-ncells-block')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1314 self.hogNCellsPerBlock = (x, x)
893 ff92801e5c54 updated hog to scikit-image 0.13 (needed to add a block_norm attribute in classifier.cfg) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 892 diff changeset	1315 self.hogBlockNorm = config.get(self.sectionHeader, 'hog-block-norm')
812 21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1316
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1317 self.speedAggregationMethod = config.get(self.sectionHeader, 'speed-aggregation-method')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1318 self.nFramesIgnoreAtEnds = config.getint(self.sectionHeader, 'nframes-ignore-at-ends')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1319 self.speedAggregationQuantile = config.getint(self.sectionHeader, 'speed-aggregation-quantile')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1320 self.minSpeedEquiprobable = config.getfloat(self.sectionHeader, 'min-speed-equiprobable')
899 1466a63dd1cf added a new classification parameter Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 895 diff changeset	1321 self.maxPercentUnknown = config.getfloat(self.sectionHeader, 'max-prop-unknown-appearance')
812 21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1322 self.maxPedestrianSpeed = config.getfloat(self.sectionHeader, 'max-ped-speed')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1323 self.maxCyclistSpeed = config.getfloat(self.sectionHeader, 'max-cyc-speed')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1324 self.meanPedestrianSpeed = config.getfloat(self.sectionHeader, 'mean-ped-speed')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1325 self.stdPedestrianSpeed = config.getfloat(self.sectionHeader, 'std-ped-speed')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1326 self.locationCyclistSpeed = config.getfloat(self.sectionHeader, 'cyc-speed-loc')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1327 self.scaleCyclistSpeed = config.getfloat(self.sectionHeader, 'cyc-speed-scale')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1328 self.meanVehicleSpeed = config.getfloat(self.sectionHeader, 'mean-veh-speed')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1329 self.stdVehicleSpeed = config.getfloat(self.sectionHeader, 'std-veh-speed')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1330
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1331 def __init__(self, filename = None):
813 ef8795dba5ed avoid crash for missing configuration filename Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 812 diff changeset	1332 if filename is not None and path.exists(filename):
812 21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1333 self.loadConfigFile(filename)
813 ef8795dba5ed avoid crash for missing configuration filename Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 812 diff changeset	1334 else:
ef8795dba5ed avoid crash for missing configuration filename Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 812 diff changeset	1335 print('Configuration filename {} could not be loaded.'.format(filename))
812 21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1336
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1337 def convertToFrames(self, frameRate, speedRatio = 3.6):
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1338 '''Converts parameters with a relationship to time in 'native' frame time
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1339 speedRatio is the conversion from the speed unit in the config file
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1340 to the distance per second
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1341
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1342 ie param(config file) = speedRatio x fps x param(used in program)
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1343 eg km/h = 3.6 (m/s to km/h) x frame/s x m/frame'''
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1344 denominator = frameRate*speedRatio
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1345 #denominator2 = denominator**2
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1346 self.minSpeedEquiprobable = self.minSpeedEquiprobable/denominator
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1347 self.maxPedestrianSpeed = self.maxPedestrianSpeed/denominator
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1348 self.maxCyclistSpeed = self.maxCyclistSpeed/denominator
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1349 self.meanPedestrianSpeed = self.meanPedestrianSpeed/denominator
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1350 self.stdPedestrianSpeed = self.stdPedestrianSpeed/denominator
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1351 self.meanVehicleSpeed = self.meanVehicleSpeed/denominator
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1352 self.stdVehicleSpeed = self.stdVehicleSpeed/denominator
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1353 # special case for the lognormal distribution
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1354 self.locationCyclistSpeed = self.locationCyclistSpeed-log(denominator)
854 33d296984dd8 rework and more info on speed probabilities for classification Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 852 diff changeset	1355 #self.scaleCyclistSpeed = self.scaleCyclistSpeed # no modification of scale
812 21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1356
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1357
664 455f9b93819c added capability to set a videofilename to movingobject and interaction, renames interactiontype to collision in interaction Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 660 diff changeset	1358 class ProcessParameters(VideoFilenameAddable):
537 6c264b914846 work on classification parameters Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 536 diff changeset	1359 '''Class for all parameters controlling data processing: input,
812 21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1360 method parameters, etc. for tracking and safety
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1361
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1362 Note: framerate is already taken into account'''
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1363
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1364 def loadConfigFile(self, filename):
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1365 from ConfigParser import ConfigParser
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1366
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1367 config = ConfigParser()
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1368 config.readfp(FakeSecHead(openCheck(filename)))
861 f9c9457b60c2 modification of storage of intrinsic camera and distortion parameters Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 854 diff changeset	1369
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1370 self.sectionHeader = config.sections()[0]
537 6c264b914846 work on classification parameters Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 536 diff changeset	1371 # Tracking/display parameters
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1372 self.videoFilename = config.get(self.sectionHeader, 'video-filename')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1373 self.databaseFilename = config.get(self.sectionHeader, 'database-filename')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1374 self.homographyFilename = config.get(self.sectionHeader, 'homography-filename')
882 4749b71aa7fb corrected bugs in storage.py when having configuration files and other files in different directories: everything should be relative to directory of command line when running the program (see FAQ) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 880 diff changeset	1375 if path.exists(self.homographyFilename):
4749b71aa7fb corrected bugs in storage.py when having configuration files and other files in different directories: everything should be relative to directory of command line when running the program (see FAQ) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 880 diff changeset	1376 self.homography = loadtxt(self.homographyFilename)
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1377 else:
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1378 self.homography = None
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1379 self.intrinsicCameraFilename = config.get(self.sectionHeader, 'intrinsic-camera-filename')
882 4749b71aa7fb corrected bugs in storage.py when having configuration files and other files in different directories: everything should be relative to directory of command line when running the program (see FAQ) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 880 diff changeset	1380 if path.exists(self.intrinsicCameraFilename):
4749b71aa7fb corrected bugs in storage.py when having configuration files and other files in different directories: everything should be relative to directory of command line when running the program (see FAQ) Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 880 diff changeset	1381 self.intrinsicCameraMatrix = loadtxt(self.intrinsicCameraFilename)
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1382 else:
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1383 self.intrinsicCameraMatrix = None
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1384 distortionCoefficients = getValuesFromINIFile(filename, 'distortion-coefficients', '=')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1385 self.distortionCoefficients = [float(x) for x in distortionCoefficients]
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1386 self.undistortedImageMultiplication = config.getfloat(self.sectionHeader, 'undistorted-size-multiplication')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1387 self.undistort = config.getboolean(self.sectionHeader, 'undistort')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1388 self.firstFrameNum = config.getint(self.sectionHeader, 'frame1')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1389 self.videoFrameRate = config.getfloat(self.sectionHeader, 'video-fps')
377 2aed569f39e7 added utils Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 344 diff changeset	1390
812 21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1391 self.classifierFilename = config.get(self.sectionHeader, 'classifier-filename')
21f10332c72b moved the classification parameters from tracking.cfg to a new classifier.cfg and made all classification parameters apparent Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 805 diff changeset	1392
537 6c264b914846 work on classification parameters Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 536 diff changeset	1393 # Safety parameters
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1394 self.maxPredictedSpeed = config.getfloat(self.sectionHeader, 'max-predicted-speed')/3.6/self.videoFrameRate
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1395 self.predictionTimeHorizon = config.getfloat(self.sectionHeader, 'prediction-time-horizon')*self.videoFrameRate
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1396 self.collisionDistance = config.getfloat(self.sectionHeader, 'collision-distance')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1397 self.crossingZones = config.getboolean(self.sectionHeader, 'crossing-zones')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1398 self.predictionMethod = config.get(self.sectionHeader, 'prediction-method')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1399 self.nPredictedTrajectories = config.getint(self.sectionHeader, 'npredicted-trajectories')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1400 self.maxNormalAcceleration = config.getfloat(self.sectionHeader, 'max-normal-acceleration')/self.videoFrameRate**2
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1401 self.maxNormalSteering = config.getfloat(self.sectionHeader, 'max-normal-steering')/self.videoFrameRate
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1402 self.minExtremeAcceleration = config.getfloat(self.sectionHeader, 'min-extreme-acceleration')/self.videoFrameRate**2
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1403 self.maxExtremeAcceleration = config.getfloat(self.sectionHeader, 'max-extreme-acceleration')/self.videoFrameRate**2
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1404 self.maxExtremeSteering = config.getfloat(self.sectionHeader, 'max-extreme-steering')/self.videoFrameRate
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1405 self.useFeaturesForPrediction = config.getboolean(self.sectionHeader, 'use-features-prediction')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1406
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1407 def __init__(self, filename = None):
813 ef8795dba5ed avoid crash for missing configuration filename Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 812 diff changeset	1408 if filename is not None and path.exists(filename):
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1409 self.loadConfigFile(filename)
813 ef8795dba5ed avoid crash for missing configuration filename Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 812 diff changeset	1410 else:
ef8795dba5ed avoid crash for missing configuration filename Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 812 diff changeset	1411 print('Configuration filename {} could not be loaded.'.format(filename))
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1412
901 753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1413 def processVideoArguments(args):
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1414 '''Loads information from configuration file
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1415 then checks what was passed on the command line
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1416 for override (eg video filename and database filename'''
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1417 if args.configFilename is not None: # consider there is a configuration file
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1418 params = ProcessParameters(args.configFilename)
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1419 videoFilename = params.videoFilename
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1420 databaseFilename = params.databaseFilename
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1421 if params.homography is not None:
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1422 invHomography = linalg.inv(params.homography)
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1423 else:
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1424 invHomography = None
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1425 intrinsicCameraMatrix = params.intrinsicCameraMatrix
935 0e63a918a1ca updated classify-objects Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 927 diff changeset	1426 distortionCoefficients = array(params.distortionCoefficients)
901 753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1427 undistortedImageMultiplication = params.undistortedImageMultiplication
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1428 undistort = params.undistort
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1429 firstFrameNum = params.firstFrameNum
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1430 else:
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1431 invHomography = None
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1432 undistort = False
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1433 intrinsicCameraMatrix = None
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1434 distortionCoefficients = []
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1435 undistortedImageMultiplication = None
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1436 undistort = False
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1437 firstFrameNum = 0
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1438
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1439 # override video and database filenames if present on command line
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1440 if args.videoFilename is not None:
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1441 videoFilename = args.videoFilename
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1442 else:
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1443 videoFilename = params.videoFilename
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1444 if args.databaseFilename is not None:
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1445 databaseFilename = args.databaseFilename
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1446 else:
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1447 databaseFilename = params.databaseFilename
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1448
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1449 return params, videoFilename, databaseFilename, invHomography, intrinsicCameraMatrix, distortionCoefficients, undistortedImageMultiplication, undistort, firstFrameNum
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1450
753a081989e2 factorized some argument handling code Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 899 diff changeset	1451 # deprecated
665 15e244d2a1b5 corrected bug with circular import for VideoFilenameAddable, moved to base module Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 664 diff changeset	1452 class SceneParameters(object):
509 935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1453 def __init__(self, config, sectionName):
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1454 from ConfigParser import NoOptionError
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1455 from ast import literal_eval
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1456 try:
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1457 self.sitename = config.get(sectionName, 'sitename')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1458 self.databaseFilename = config.get(sectionName, 'data-filename')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1459 self.homographyFilename = config.get(sectionName, 'homography-filename')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1460 self.calibrationFilename = config.get(sectionName, 'calibration-filename')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1461 self.videoFilename = config.get(sectionName, 'video-filename')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1462 self.frameRate = config.getfloat(sectionName, 'framerate')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1463 self.date = datetime.strptime(config.get(sectionName, 'date'), datetimeFormat) # 2011-06-22 11:00:39
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1464 self.translation = literal_eval(config.get(sectionName, 'translation')) # = [0.0, 0.0]
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1465 self.rotation = config.getfloat(sectionName, 'rotation')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1466 self.duration = config.getint(sectionName, 'duration')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1467 except NoOptionError as e:
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1468 print(e)
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1469 print('Not a section for scene meta-data')
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1470
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1471 @staticmethod
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1472 def loadConfigFile(filename):
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1473 from ConfigParser import ConfigParser
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1474 config = ConfigParser()
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1475 config.readfp(openCheck(filename))
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1476 configDict = dict()
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1477 for sectionName in config.sections():
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1478 configDict[sectionName] = SceneParameters(config, sectionName)
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1479 return configDict
935430b1d408 corrected mask bug in feature tracking, updated display-trajectories to display on undistorted image Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 491 diff changeset	1480
382 ba813f148ade development for clustering Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 377 diff changeset	1481
215 5e2983b05d4e created first doctest tests for storage Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 212 diff changeset	1482 if __name__ == "__main__":
5e2983b05d4e created first doctest tests for storage Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 212 diff changeset	1483 import doctest
5e2983b05d4e created first doctest tests for storage Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 212 diff changeset	1484 import unittest
5e2983b05d4e created first doctest tests for storage Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 212 diff changeset	1485 suite = doctest.DocFileSuite('tests/storage.txt')
5e2983b05d4e created first doctest tests for storage Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 212 diff changeset	1486 unittest.TextTestRunner().run(suite)
72 575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1487 # #doctest.testmod()
575340e6fce3 corrected most of the method to load NGSIM data (adapted to the current MovingObject class Nicolas Saunier <nicolas.saunier@polymtl.ca> parents: 50 diff changeset	1488 # #doctest.testfile("example.txt")

Mercurial Hosting > traffic-intelligence

annotate python/storage.py @ 935:0e63a918a1ca