JACOBIAN_FILTERED: J= jacobian_filtered( fwd_model, img) Filter a jacobian matrix by a specified filter function INPUT: fwd_model = forward model fwd_model.jacobian_filtered.jacobian = actual jacobian function (J0) fwd_model.jacobian_filtered.filter = Filter Matrix (F) img = image background for jacobian calc OUTPUT: J = Jacobian matrix = F*J0
0001 function J= jacobian_filtered( fwd_model, img) 0002 % JACOBIAN_FILTERED: J= jacobian_filtered( fwd_model, img) 0003 % 0004 % Filter a jacobian matrix by a specified filter function 0005 % INPUT: 0006 % fwd_model = forward model 0007 % fwd_model.jacobian_filtered.jacobian = actual jacobian function (J0) 0008 % fwd_model.jacobian_filtered.filter = Filter Matrix (F) 0009 % img = image background for jacobian calc 0010 % OUTPUT: 0011 % J = Jacobian matrix = F*J0 0012 0013 % (C) 2009 Andy Adler. License: GPL version 2 or version 3 0014 % $Id: jacobian_filtered.m 3101 2012-06-08 14:34:08Z bgrychtol $ 0015 0016 J0 = feval(fwd_model.jacobian_filtered.jacobian, ... 0017 fwd_model, img); 0018 F = fwd_model.jacobian_filtered.filter; 0019 0020 J= F*J0; 0021