%DCT变换测试程序DCTTest.m
Cls;
I=imread('lena.bmp');
%s=double(I(1:128,1:128));
s=[200 202 189 188 189 175 175 175
200 203 198 188 189 182 178 175
203 200 200 195 200 187 185 175
200 200 200 200 197 187 187 187
200 205 200 200 195 188 187 175
200 200 200 200 200 190 187 175
205 200 199 200 191 187 187 175
210 200 200 200 188 185 187 186]
%F=Image2DCT(s,128,128);
F=Image2DCT(s,8,8);
%MF=abs(log(abs(F)+1));%将幅值动态范围缩小,便于显示
%figure;imshow(MF,[]);
figure;imshow(F,[]);