tidenoisehard2.m

来自「denoising image with DWT」· M 代码 · 共 41 行

function out = TIDenoiseHard2(Noisy,L,qmf,thresh)
%  TIDenoiseHard2 -- Translation-Invariant Hard Wavelet Thresholding of
%		Images.
%  Usage 
%    out=TIDenoiseHard2(Noisy,L,qmf,sigma)
%  Inputs
%    Noisy      Noisy image.
%    L  	Deepness of decomposition.
%    qmf    	Quadrature Mirror Filter for Wavelet Transform
%    thresh     Threshold to apply on Wavelet Coefficients
%  Outputs 
%    out    	Restored image.
%

[n,J]=quadlength(Noisy);
out=zeros(n,n);
lim=8;
for i=1:lim,
	for j=1:lim,
		Noistrans = cyclespin2(Noisy,i,j);
		dout = HT2(Noistrans,L,qmf,thresh);
		dout = cyclespin2(dout,n-i,n-j);
		out = out+dout;
	end
end
out=out/(lim*lim);

% Written by Maureen Clerc and Jerome Kalifa, 1997
% [e-mail], [e-mail]

    
    
 
 
%
%  Part of Wavelab Version 850
%  Built Tue Jan  3 13:20:40 EST 2006
%  This is Copyrighted Material
%  For Copying permissions see COPYING.m
%  Comments? e-mail [e-mail]