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)
