ProgowanieTes

clear all;

T = 60;

I = imread('imgnoise.jpg');

RI = I(:,:,1);
GI = I(:,:,2);
BI = I(:,:,3);

%transforamta
RIt = dct2(RI);
BIt = dct2(BI);
GIt = dct2(GI);
RIf = fwht(double(RI));
BIf = fwht(double(BI));
GIf = fwht(double(GI));
[CR,SR]=wavedec2(RI,2,'haar');
[CG,SG]=wavedec2(GI,2,'haar');
[CB,SB]=wavedec2(BI,2,'haar');

%progowanie
RIt(abs(RIt) < T) = 0;
BIt(abs(BIt) < T) = 0;
GIt(abs(GIt) < T) = 0;

RIf(abs(RIf) < T) = 0;
BIf(abs(BIf) < T) = 0;
GIf(abs(GIf) < T) = 0;

CR(abs(CR) < T) = 0;
CG(abs(CG) < T) = 0;
CB(abs(CB) < T) = 0;


%odwrotana trtanformata
RIti = idct2(RIt);
BIti = idct2(BIt);
GIti = idct2(GIt);
RIfi = ifwht(RIf);
BIfi = ifwht(BIf);
GIfi = ifwht(GIf);
Rwi = waverec2(CR, SR, 'haar');
Bwi = waverec2(CB, SB, 'haar');
Gwi = waverec2(CG, SG, 'haar');

finishedT(:,:,1) = RIti;
finishedT(:,:,2) = GIti;
finishedT(:,:,3) = BIti;
finishedT = uint8(finishedT);

finishedF(:,:,1) = RIfi;
finishedF(:,:,2) = GIfi;
finishedF(:,:,3) = BIfi;
finishedF = uint8(finishedF);

finishedW(:,:,1) = Rwi;
finishedW(:,:,2) = Gwi;
finishedW(:,:,3) = Bwi;
finishedW = uint8(finishedW);

f = figure;
imshowpair(I,finishedT,'montage');
figure, imshow(I-finishedT)
f = figure;
imshowpair(I,finishedF,'montage');
f = figure;
imshowpair(I,finishedW,'montage');