www.pudn.com > Tikhonov_Regularization_for_super_resolution.rar > WTWZ.m, change:2010-07-06,size:729b


function G1 = WTWZ(LR, Zn, M, B, S) 
% computer the term W'*W*Zn 
% input£º 
% LR: the sequence of low resolution images 
% Zn: the current high resolution image 
% M£º registration parameter 
% B£ºPSF 
% S: sampling factor 
 
% bluring 
Xn = imfilter(Zn, B, 'symmetric'); 
 
% shift and sampling 
Htemp = zeros(size(LR)); 
for k=1:size(LR,3) 
  Htemp(:,:,k)=Xn(M(k, 2):S:(size(LR, 1) - 1) * S + M(k, 2), M(k, 1):S:(size(LR, 2) -1) * S + M(k, 1)); 
end 
 
% shift and sampling after transformation 
HRsd = zeros([size(Zn) size(LR, 3)]); 
for k=1:size(LR,3) 
  HRsd(M(k, 2):S:(size(LR, 1) - 1) * S + M(k,2), M(k, 1):S:(size(LR, 2) - 1) * S + M(k, 1), k) = Htemp(:,:,k); 
end 
 
% debluring 
G1 = imfilter(HRsd, B', 'symmetric');