kvetsc
kvetn = kvet/256
imwrite(kvetn,'kvet.bmp');

I=imread('kvet.bmp');

I = im2double(I);

imshow(I,'InitialMagnification','fit');

dctI = dct2(I);
imshow(dctI,'InitialMagnification','fit');

dctI(6:8,6:8) = 0; % omezení frekvencí (filtrace)

% inverzní transformace
invDctI = idct2(dctI);
imshow(invDctI,'InitialMagnification','fit');

% Proveďte diskrétní kosinovou transformaci pomocí blokových procedur.