dtcwtmse.zip_wavelet threshold_双数复_双数复小波

% function [dtcwtpsnr]=dtcwtmse(im, cth); clear % im=imread ('cameraman.tif'); % im=imread ('lena.png'); % im=imread ('rice.png'); % im=imread ('peppers.png'); im=imread ('barbara.png'); % im=imread ('boat.jpg'); % im=imread ('baboon.jpg'); % im=imread ('peppers.png'); sigma=10; im=double(im); nim=im+sigma*randn(size(im)); [Faf, Fsf] = FSfarras; [af, sf] = dualfilt1; % J=4; J=3; cth=3*sigma; y = dualtree2D(nim, J, Faf, af); % y1={y{1,9}{1,1},y{1,8}{1,1},y{1,7}{1,1},y{1,6}{1,1},y{1,5}{1,1},y{1,4}{1,1},y{1,3}{1,1},y{1,2}{1,1},y{1,1}{1,1}}; % y2={y{1,9}{1,2},y{1,8}{1,2},y{1,7}{1,2},y{1,6}{1,2},y{1,5}{1,2},y{1,4}{1,2},y{1,3}{1,2},y{1,2}{1,2},y{1,1}{1,2}}; % y1={y{1,5}{1,1},y{1,4}{1,1},y{1,3}{1,1},y{1,2}{1,1},y{1,1}{1,1}}; % y2={y{1,5}{1,2},y{1,4}{1,2},y{1,3}{1,2},y{1,2}{1,2},y{1,1}{1,2}}; y1={y{1,4}{1,1},y{1,3}{1,1},y{1,2}{1,1},y{1,1}{1,1}}; y2={y{1,4}{1,2},y{1,3}{1,2},y{1,2}{1,2},y{1,1}{1,2}}; [c1, s1] = pdfb2vec(y1); [c2, s2] = pdfb2vec(y2); c1(s1(1, 3) * s1(1, 4):end)=c1(s1(1, 3) * s1(1, 4):end).*((abs(c1(s1(1, 3) * s1(1, 4):end)) > cth)); c2(s2(1, 3) * s2(1, 4):end)=c2(s2(1, 3) * s2(1, 4):end).*((abs(c2(s2(1, 3) * s2(1, 4):end)) > cth)); y1 = vec2pdfb(c1, s1); y2 = vec2pdfb(c2, s2); % a={y1{1,9},y2{1,9}}; % b={y1{1,8},y2{1,8}}; % c={y1{1,7},y2{1,7}}; % d={y1{1,6},y2{1,6}}; % e={y1{1,5},y2{1,5}}; f={y1{1,4},y2{1,4}}; g={y1{1,3},y2{1,3}}; h={y1{1,2},y2{1,2}}; i={y1{1,1},y2{1,1}}; % y0={a,b,c,d,e,f,g,h,i}; % y0={e,f,g,h,i}; y0={f,g,h,i}; x = idualtree2D(y0, J, Fsf, sf); % mse =sqrt( sum( sum( (x - double(x0)).^2 ) )/ prod(size(x))); e=im-x; MSE = sqrt(sum(e(:).*e(:))/prod(size(e))); % consnr=SNR(double(x0), x); % dtcwtsnr=SNR(im, x); dtcwtpsnr=10*log10(255.^2/MSE);