TIME_PRIOR_SOLVE inverse solver to account for time differences img= time_prior_solve( 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 both data1 and data2 may be matrices (MxT) each of M measurements at T times if either data1 or data2 is a vector, then it is expanded to be the same size matrix
0001 function img= time_prior_solve(inv_model, varargin ) 0002 % TIME_PRIOR_SOLVE inverse solver to account for time differences 0003 % img= time_prior_solve( 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 % both data1 and data2 may be matrices (MxT) each of 0010 % M measurements at T times 0011 % if either data1 or data2 is a vector, then it is expanded 0012 % to be the same size matrix 0013 0014 % TODO: This function really should be calling the proper 0015 % prior calculator functions, and not reimplementing 0016 % them internally 0017 0018 % (C) 2005 Andy Adler. License: GPL version 2 or version 3 0019 % $Id: time_prior_solve.m 5700 2018-03-19 15:37:06Z aadler $ 0020 0021 warning('EIDORS:deprecated','TIME_PRIOR_SOLVE is deprecated as of 08-Jun-2012. Use INV_SOLVE_TIME_PRIOR instead.'); 0022 0023 if isfield(inv_model,'time_prior_solve'); 0024 inv_model.inv_solve_time_prior = inv_model.time_prior_solve; 0025 end 0026 if isfield(inv_model,'time_smooth_prior'); 0027 inv_model.prior_time_smooth = inv_model.time_smooth_prior; 0028 end 0029 0030 img = inv_solve_time_prior(inv_model, varargin{:});