clear all
close all

im_org = imread('prostata2.png'); %obraz oryginalny
im_n = imnoise(im_org, 'gaussian', 0, 0.000001) ; %obraz zaszumiony

wname = 'rbio1.1'; %rodzaj falki
wname2 ='rbio3.5';
wname3 = 'rbio6.8';
level = 5;
[C,S] = wavedec2(im_n,level,wname);
[C2,S2] = wavedec2(im_n,level,wname2);
[C3,S3] = wavedec2(im_n,level,wname3);

thr = wthrmngr('dw2ddenoLVL','penalhi',C,S,3);
thr2 = wthrmngr('dw2ddenoLVL','penalhi',C2,S2,3);
thr3 = wthrmngr('dw2ddenoLVL','penalhi',C3,S3,3);
sorh = 's';
[XDEN,cfsDEN,dimCFS] = wdencmp('lvd',C,S,wname,level,thr,sorh);
[XDEN2,cfsDEN,dimCFS] = wdencmp('lvd',C2,S2,wname2,level,thr2,sorh);
[XDEN3,cfsDEN,dimCFS] = wdencmp('lvd',C3,S3,wname3,level,thr3,sorh);

figure;
subplot(2,3,1);
imagesc(im_org); colormap gray; axis off;
title('Obraz oryginalny');
subplot(2,3,2);
imagesc(im_n); colormap gray; axis off;
title('Obraz zaszumiony');
subplot(2,3,4);
imagesc(XDEN); colormap gray; axis off;
title(' Biorthogonal 1.1');
subplot(2,3,5);
imagesc(XDEN2); colormap gray; axis off;
title(' Biorthogonal 3.5');
subplot(2,3,6);
imagesc(XDEN3); colormap gray; axis off;
title(' Biorthogonal 6.8');
% parametry oceny
[peaksnr, snr] = psnr(XDEN, double(im_org))
[peaksnr2, snr2] = psnr(XDEN2,double(im_org))
[peaksnr3, snr3] = psnr(XDEN3,double(im_org))

err = immse(XDEN,double(im_org))
err2 = immse(XDEN2,double(im_org))
err3 = immse(XDEN3,double(im_org))