iwt2d.m

来自「该代码为图像融合的vc源程序,对学习图像融合有很重要的价值」· M 代码 · 共 47 行

function y=iwt2d(w,lpfy,hpfy,is2y,lpfx,hpfx,is2x,F_zero,F_one)

%$Id: iwt2d.m,v 1.1.1.1 2000/10/16 09:30:40 mattb Exp $
%===========================================================================
%                   Matt Bowers - University of Bath
%
%                        Copyright (c) 1996
%===========================================================================
%
%  Two dimensional inverse wavelet transform of w
%
%  edge affects handling of edges:
%
%       = 1 - use reflection appropriate to odd or even length
%       = 0 - use wraparound
%
%===========================================================================


%edge affects handling of edges:
%		= 1 - use reflection appropriate to odd or even length
%		= 0 - use wraparound

[lpfx_rows lpfx_cols]=size(lpfx);
[lpfy_rows lpfy_cols]=size(lpfy);
[hpfx_rows hpfx_cols]=size(hpfx);
[hpfy_rows hpfy_cols]=size(hpfy);

[i_rows i_cols]=size(w);

%Column transforms
%disp ('Processing columns...')
for c=1:i_cols
  col=[w(:,c) ; 0]; %Column + zero token
  w(:,c)=(reshape(col(lpfy),lpfy_rows,lpfy_cols)*F_zero) + ((reshape(-col(hpfy),hpfy_rows,hpfy_cols).*is2y)*F_one);
end


%Row transforms
%disp('Processing rows...')
for r=1:i_rows
  row=[w(r,:) 0];	%Row + zero token
  w(r,:)=(reshape(row(lpfx),lpfx_rows,lpfx_cols)*F_zero)' + ((reshape(-row(hpfx),hpfx_rows,hpfx_cols).*is2x)*F_one)';
end

y=w;