TIKHONOV_IMAGE_PRIOR calculate image prior Reg= tikhonov_image_prior( inv_model ) Reg => output regularization term inv_model => inverse model struct
0001 function Reg= tikhonov_image_prior( inv_model ); 0002 % TIKHONOV_IMAGE_PRIOR calculate image prior 0003 % Reg= tikhonov_image_prior( inv_model ) 0004 % Reg => output regularization term 0005 % inv_model => inverse model struct 0006 0007 % (C) 2005 Andy Adler. License: GPL version 2 or version 3 0008 % $Id: tikhonov_image_prior.html 2819 2011-09-07 16:43:11Z aadler $ 0009 0010 %pp= aa_fwd_parameters( inv_model.fwd_model ); 0011 0012 if isfield( inv_model.fwd_model, 'coarse2fine' ) 0013 no_dof = size(inv_model.fwd_model.coarse2fine,2); 0014 else 0015 no_dof = size(inv_model.fwd_model.elems,1); 0016 end 0017 0018 Reg = speye( no_dof ); 0019