Home > Matching_SYM_LSM > src > 2D_LSM_62 > 2D_LSM_62_sym_Functions > image_pair_analyse.m

image_pair_analyse

PURPOSE ^

% analyse result of matching image pair

SYNOPSIS ^

function [outputArg1,outputArg2] = image_pair_analyse(A_a,R_a,Ai_est,Red,est_sigma_0,est_x,plot_type,test_symmetry)

DESCRIPTION ^

% analyse result of matching image pair

 wf 7/2020

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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 %% analyse result of matching image pair
0010 %
0011 % wf 7/2020
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     %% provide result
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

Generated on Sun 19-Jul-2020 23:00:25 by m2html © 2005