0001 function [outputArg1,outputArg2] = image_pair_analyse(...
0002 A_a, ...
0003 R_a, ...
0004 Ai_est, ...
0005 Red, ...
0006 est_sigma_0, ...
0007 est_x, ...
0008 plot_type,test_symmetry)
0009
0010
0011
0012
0013 if Red > 0
0014 if plot_type > 1
0015 A_est= [reshape(est_x.x(1:6),2,3);0 0 1]
0016 A_estimated_corrections = A_est - A_a
0017 R_est = est_x.x(7:8)'
0018 R_estimated_corrections = R_est -R_a
0019 end
0020 if plot_type > 0
0021 if test_symmetry
0022 check_symmetry_BBi_I = A_est*Ai_est-eye(3)
0023 check_symmetry_SSi_I = Rm_est*Rmi_est-eye(2)
0024 end
0025 end
0026
0027
0028
0029
0030 Redundancy = Red
0031 est_sigma_0 = est_sigma_0
0032 disp('Geometric transformation [1, 1, [pixel]]')
0033 sigmas = diag(sqrt(est_x.C))';
0034 A_est = reshape(est_x.x(1:6),2,3)
0035 A_corrections_wrt_approximate_values = A_est - A_a(1:2,:)
0036 Sigmas_A = reshape(sigmas(1:6),2,3)
0037 est_Sigmas_A = Sigmas_A*est_sigma_0
0038 disp('Intensity transformation [1,[gr]]')
0039 R_est = [est_x.x(7),255*est_x.x(8)]
0040 R_corrections_wrt_approximate_values = R_est - R_a
0041 Sigmas_R = [sigmas(7),255*sigmas(8)]
0042 est_Sigmas_R = est_sigma_0*[sigmas(7),255*sigmas(8)]
0043
0044
0045 else
0046 disp('No solution found')
0047 end
0048 end
0049
0050