reg_dir = "/Users/isa/Dropbox/MultiModalRegPaper" results_dir = reg_dir + "/DenseResults/Aerial2Aerial" gt_dir = reg_dir + "/ground_truth" result_summary_file = results_dir + "/error_summary.txt" Tfos = open(result_summary_file, 'w') bh_fname = results_dir + "/BH_vsi_06_finer.txt" gt_bh_fname = gt_dir + "/aerial-aerial/bh_vsi_06_Hs.txt" geo_bh_fname = gt_dir + "/aerial-lidar-demeaned/BH_2006/Hs_geo.txt" error = reg3d.transformation_error_general(fname = bh_fname, gt_fname = gt_bh_fname, geo_fname = geo_bh_fname ) Tfos.write("-----------------------------------------------------\n") Tfos.write("*********** Error BH_VSI_2006 ***********************\n") Tfos.write("-----------------------------------------------------\n") np.savetxt( Tfos , error.reshape((1,3)) , fmt='%.7g', delimiter=' ') capitol_fname = results_dir + "/capitol_11_06_finer.txt" gt_capitol_fname = gt_dir + "/aerial-aerial/capitol_11_06_Hs.txt" geo_capitol_fname = gt_dir + "/aerial-lidar-demeaned/capitol_2006/Hs_geo.txt" error = reg3d.transformation_error_general(fname = capitol_fname, gt_fname = gt_capitol_fname, geo_fname = geo_capitol_fname )
inital_error = np.array([[ 0.0, 3.64969371, 3.881038 ], [ 0.0, 3.71314362, 6.67001667], [ 0.0, 3.71530045, 7.37906836], [ 0.0, 4.11227278, 6.20087232], [ 0.0, 3.84084958, 2.51475658]]) for site_idx in range(0, len(sites_vj)): for trial in range(0,5): print 'Trial: ' + str(trial) results_dir = dense_dir + "/" + sites_vj[site_idx ] t_fname = results_dir + "/" + file_base[site_idx] + "_H" + str(trial) + ".txt" gt_fname = reg_dir + "/ground_truth/aerial-lidar-demeaned/" + sites_isa[site_idx] + "/Hs_geo_matrix.txt" error[trial,:, site_idx ] = reg3d.transformation_error_general(fname = t_fname, gt_fname = gt_fname) Tfos = open(result_summary_file, 'w') for site_idx in range(0, len(sites_vj)): print "-----------------------------------------------------" print "*********** Errors " + sites_vj[site_idx] + " *************" print "-----------------------------------------------------" print error[:,:,site_idx] print 'Mean:' mean_error = np.mean(abs(error[:,:,site_idx]), axis=0) print mean_error
#!/usr/bin/env python # encoding: utf-8 """ Author: Isabel Restrepo Script to find the average alignment error of manual registrationπ """ import os import sys import numpy as np from scipy import stats import reg3d ntrials = 31 error = np.zeros((ntrials-1, 3)); for trial in range(1,ntrials): print 'Trial: ' + str(trial) fname = "/data/lidar_providence/manual_error/Hs" + str(trial) + ".txt" gt_fname = "/data/lidar_providence/manual_error/Hs-identity.txt" error[trial-1,:] = reg3d.transformation_error_general(fname = fname, gt_fname = gt_fname) print 'Errors:' print error print 'Mean:' print np.mean(abs(error), axis=0) print 'Std:' print np.std(abs(error), axis=0)