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

simulated_prepare_true_values

PURPOSE ^

% prepare true and approximate values

SYNOPSIS ^

function par = simulated_prepare_true_values(par)

DESCRIPTION ^

% prepare true and approximate values

 wf 7/2020

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function par = simulated_prepare_true_values(par)
0002 %% prepare true and approximate values
0003 %
0004 % wf 7/2020
0005 
0006 % --- generate variance function (LUT) ----------------------------------
0007 range_g         = 0:255;
0008 par.vg          = var_g(range_g(:),par.stg);
0009 par.vh          = var_h(range_g(:),par.sth);
0010 % include rounding
0011 par.vg = (par.vg+1/12*ones(256,1));
0012 par.vh = (par.vh+1/12*ones(256,1));
0013 % maximum sigma (for determining range of true image)
0014 par.sigma_max   = sqrt(max([par.vg(:);par.vh(:)]));
0015 
0016 % --- true affinity referring to centre of windows ----------------------
0017 par.A_true   = [par.A_t; 0 0 1];
0018 % half transformation A=B^2
0019 par.B_true   = sqrtm(par.A_true);
0020 par.Bi_true  = inv(par.B_true);
0021 
0022 % --- true radiometric parameters R=S^2 ---------------------------------
0023 par.a7_true = sqrt(par.R_true(1));
0024 par.a8_true = par.R_true(2)/(1+par.a7_true);
0025 par.S_true  = [par.a7_true,par.a8_true];
0026 
0027 % --- true parameters ---------------------------------------------------
0028 par.x_true = [reshape(par.A_true(1:2,:),6,1);par.R_true'];
0029 
0030 % --- approximate values ------------------------------------------------
0031 par.A_a  = par.A_true;
0032 par.Ai_a = inv(par.A_a);
0033 par.R_a  = par.R_true;
0034 par.Ri_a =  1/(par.R_a(1,1))*[1,-par.R_a(2)];
0035 
0036 end
0037

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