Mercurial Hosting > traffic-intelligence
view scripts/safety-analysis.py @ 930:7db0f2853bfd
corrected projection back to image space
author | Nicolas Saunier <nicolas.saunier@polymtl.ca> |
---|---|
date | Thu, 13 Jul 2017 00:31:54 -0400 |
parents | 063d1267585d |
children | a2f3f3ca241e |
line wrap: on
line source
#! /usr/bin/env python import storage, prediction, events, moving import sys, argparse, random import matplotlib.pyplot as plt import numpy as np # todo: very slow if too many predicted trajectories # add computation of probality of unsucessful evasive action parser = argparse.ArgumentParser(description='The program processes indicators for all pairs of road users in the scene') parser.add_argument('--cfg', dest = 'configFilename', help = 'name of the configuration file', required = True) parser.add_argument('-n', dest = 'nObjects', help = 'number of objects to analyse', type = int) # TODO analyze only parser.add_argument('--prediction-method', dest = 'predictionMethod', help = 'prediction method (constant velocity (cvd: vector computation (approximate); cve: equation solving; cv: discrete time (approximate)), normal adaptation, point set prediction)', choices = ['cvd', 'cve', 'cv', 'na', 'ps']) parser.add_argument('--pet', dest = 'computePET', help = 'computes PET', action = 'store_true') parser.add_argument('--display-cp', dest = 'displayCollisionPoints', help = 'display collision points', action = 'store_true') parser.add_argument('--nthreads', dest = 'nProcesses', help = 'number of processes to run in parallel', type = int, default = 1) args = parser.parse_args() params = storage.ProcessParameters(args.configFilename) # parameters for prediction methods if args.predictionMethod is not None: predictionMethod = args.predictionMethod else: predictionMethod = params.predictionMethod def accelerationDistribution(): return random.triangular(-params.maxNormalAcceleration, params.maxNormalAcceleration, 0.) def steeringDistribution(): return random.triangular(-params.maxNormalSteering, params.maxNormalSteering, 0.) if predictionMethod == 'cvd': predictionParameters = prediction.CVDirectPredictionParameters() if predictionMethod == 'cve': predictionParameters = prediction.CVExactPredictionParameters() elif predictionMethod == 'cv': predictionParameters = prediction.ConstantPredictionParameters(params.maxPredictedSpeed) elif predictionMethod == 'na': predictionParameters = prediction.NormalAdaptationPredictionParameters(params.maxPredictedSpeed, params.nPredictedTrajectories, accelerationDistribution, steeringDistribution, params.useFeaturesForPrediction) elif predictionMethod == 'ps': predictionParameters = prediction.PointSetPredictionParameters(params.maxPredictedSpeed) # no else required, since parameters is required as argument # evasiveActionPredictionParameters = prediction.EvasiveActionPredictionParameters(params.maxPredictedSpeed, # params.nPredictedTrajectories, # params.minExtremeAcceleration, # params.maxExtremeAcceleration, # params.maxExtremeSteering, # params.useFeaturesForPrediction) objects = storage.loadTrajectoriesFromSqlite(params.databaseFilename, 'object', args.nObjects, withFeatures = (params.useFeaturesForPrediction or (predictionMethod == 'ps'))) # if params.useFeaturesForPrediction: # features = storage.loadTrajectoriesFromSqlite(params.databaseFilename,'feature') # needed if normal adaptation # for obj in objects: # obj.setFeatures(features) interactions = events.createInteractions(objects) for inter in interactions: inter.computeIndicators() inter.computeCrossingsCollisions(predictionParameters, params.collisionDistance, params.predictionTimeHorizon, params.crossingZones, nProcesses = args.nProcesses) if args.computePET: for inter in interactions: inter.computePET(params.collisionDistance) storage.saveIndicatorsToSqlite(params.databaseFilename, interactions) if args.displayCollisionPoints: plt.figure() allCollisionPoints = [] for inter in interactions: for collisionPoints in inter.collisionPoints.values(): allCollisionPoints += collisionPoints moving.Point.plotAll(allCollisionPoints) plt.axis('equal')