例9-17 .m

来自「matlab基础与提高」· M 代码 · 共 18 行

I=imread(‘cameraman.itf’)
I=im2double(I);
T=dctmtx(8);
B=blkproc(I,[8 8],’P1*x*p2’,T,T’);
mask=
[1 1 1 1 0 0 0 0
1 1 1 0 0 0 0 0
1 1 0 0 0 0 0 0
1 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0];
B2=blkproc(B,[8 8],’P1.*x’,mask);
I2=blkproc(B2,[8 8],’P1*x*P2’,T’,T);
imshow(I);
figure,imshow(I2)