#'C:\\Users\\pkhomchuk\\My Projects\\UAS - FAA Software\\Tracker\\Radar_3D_4' #'C:\\Users\\pkhomchuk\\My Projects\\UAS - FAA Delivery\\UAS - FAA Software\\Tracker\\Maneuver', #'C:\\Users\\pkhomchuk\\My Projects\\UAS - FAA Delivery\\UAS - FAA Software\\Tracker\\Straight', #'C:\\Users\\pkhomchuk\\My Projects\\UAS - FAA Delivery\\UAS - FAA Software\\Tracker\s\55', '' ) TGT_BEACON = 2332 NMI2FT = 6076.11548 def listFilesInDir( dir, ext ) : os.chdir( dir ) return glob.glob( '*.' + ext ) errors = [] for ind, folder in enumerate( output_folders ) : if folder != '' : residual_rng_rate = [] output_file_names = listFilesInDir( folder, 'txt' ) truth_ind = 0 for filename in output_file_names : output = OutputFile.read( folder + "\\" + filename ) trk = OutputFile.Track( output, beacon = TGT_BEACON, sensor = 'radar' ); if trk.dataIsGood : residuals = np.ndarray.tolist( NMI2FT * np.transpose( trk.residuals )) residual_rng_rate += residuals[ 0 ] plt.figure( 1 ) plt.boxplot( residual_rng_rate ) plt.show()
truth_file_names = listFilesInDir(folder, 'txt') for filename in truth_file_names: truth = TruthFile.read(folder + '\\' + filename) tgt_true[ind].append(TruthFile.Aircraft(truth, beacon=TGT_BEACON)) own_true[ind].append(TruthFile.Aircraft(truth, beacon=OWN_BEACON)) errors = [] for ind, folder in enumerate(output_folders): if folder != '': hist_errors = TrkErrors.HistErrors() output_file_names = listFilesInDir(folder, 'txt') truth_ind = 0 for filename in output_file_names: output = OutputFile.read(folder + "\\" + filename) trk = OutputFile.Track(output, rack_id=1, track_type=2) if trk.dataIsGood: truth_ind = int( math.floor( (int(filename.split('.')[0][3:]) - 1) / N_MONTE_CARLO)) err = TrkErrors.Errors(own_true[ind][truth_ind], tgt_true[ind][truth_ind], trk, True) hist_errors = hist_errors + err.hist_errors #if ind % N_MONTE_CARLO == 0 and ind > 0 : # truth_ind += 1 errors.append(hist_errors) n = len(errors) # Plot horizontal position error gs = gridspec.GridSpec(1, n)