PDIPM_DIFF inverse solver for difference data using Primal/Dual interior point method img= ab_pdipm( inv_model, data1, data2) img => output image (or vector of images) inv_model => inverse model struct data1 => differential data at earlier time data2 => differential data at later time inv_model.pdipm_diff.norm_data 1 or 2 (DEFAULT 2) inv_model.pdipm_diff.norm_image 1 or 2 (DEFAULT 2) inv_model.pdipm_diff.beta (default 1e-6) Parameters: max_iters = inv_model.parameters.max_iterations (default 10) Max number of iterations before stopping min change = inv_model.parameters.min_change (default 0) Min Change in objective fcn (norm(y-Jx)^2 + hp*TV(x)) before stopping beta is the parameter that smooths the TV functional
0001 function img=pdipm_diff( varargin ) 0002 % PDIPM_DIFF inverse solver for difference data using Primal/Dual interior point method 0003 % img= ab_pdipm( inv_model, data1, data2) 0004 % img => output image (or vector of images) 0005 % inv_model => inverse model struct 0006 % data1 => differential data at earlier time 0007 % data2 => differential data at later time 0008 % 0009 % inv_model.pdipm_diff.norm_data 1 or 2 (DEFAULT 2) 0010 % inv_model.pdipm_diff.norm_image 1 or 2 (DEFAULT 2) 0011 % inv_model.pdipm_diff.beta (default 1e-6) 0012 % 0013 % Parameters: 0014 % max_iters = inv_model.parameters.max_iterations (default 10) 0015 % Max number of iterations before stopping 0016 % min change = inv_model.parameters.min_change (default 0) 0017 % Min Change in objective fcn (norm(y-Jx)^2 + hp*TV(x)) before stopping 0018 % beta is the parameter that smooths the TV functional 0019 0020 % (C) 2008 Andrea Borsic. License: GPL version 2 or version 3 0021 % $Id: pdipm_diff.m 3544 2012-07-25 12:51:47Z aadler $ 0022 0023 0024 warning('EIDORS:deprecated','PDIPM_DIFF is deprecated as of 08-Jun-2012. Use INV_SOLVE_DIFF_PDIPM instead.'); 0025 0026 if isfield(varargin{1},'pdipm_diff'); 0027 varargin{1}.inv_solve_diff_pdipm = varargin{1}.pdipm_diff; 0028 end 0029 0030 img = inv_solve_diff_pdipm(varargin{:});