Medical Image Processing - Digital Signal Processing for Medical Imaging Using Matlab - page 80

Fig. 4.8 Illustration of the original image and the corresponding compressed image obtained using

2D-DCT

4.3.1.1 imagecompusingdct.m

load mri
D=double(D);
V=D(:,:,1,1);
V=double(V);
V1=blkproc(V,[8 8],'dctt(x)');
V2=blkproc(V1,[8 8],'idctt(x)');
figure
colormap(gray)
subplot(1,2,1)
imagesc(abs(V))
subplot(1,2,2)
imagesc(abs(V2))

4.3.1.2 dctt.m

function [res]=dctt(x)
%Function called by the imagecompusingdct.m file
x=x-128;
x=dct2(x);
JPEGQ=[16 11 10 16 24 40 51 61;12 12 14 19 26 58 60 55;14 13