图形图像处理[第8章节].ppt

图形图像处理 信息与计算科学系 张宝雷 上节内容回顾 正交矩阵 多元函数极值 积分变换 3.4 变换编码 3.4.4 DCT编码的MATLAB示例 例3-8 cr=0.5; I=imresize(rgb2gray(imread(chinesegirl.jpg)),[512 512]); I1=double(I)/255; figure(1); imshow(I1); fftcoe=blkproc(I1,[8 8],fft2(x)); coevar=im2col(fftcoe,[8 8],distinct); coe=coevar; [y,ind]=sort(coevar); [m,n]=size(coevar); snum=64-64*cr; 例3-8 for i=1:n coe(ind(1:snum),i)=0; end b2=col2im(coe,[8 8],[512 512],distinct); I2=blkproc(b2,[8 8],ifft2(x)); figure(2); imshow(I2); e=double(I1)-double(I2); [m,n]=size(e); erms=sqrt(sum(e(:).^2)/(m*n)) 运行结果 erms = 0.0085 - 0.0000i 为什么看不到呢？ 分析 subplot(121);plot(abs(ccd(1:40))) subplot(122);plot(abs(ccd(40:length(ccd)))) 分析 subplot(121);plot(abs(ccd(1:40))) subplot(122);plot(abs(ccd(40:length(ccd)))) 频谱分析 示例1： subplot(221);plot(fft(ccc(40,:))); subplot(221);plot(ccc(40,:)); subplot(222);plot(abs(fft(ccc(40,:)))); subplot(222);plot(angle(fft(ccc(40,:)))); subplot(222);plot(abs(fft(ccc(40,:)))); subplot(223);plot(abs(fft(ccc(40,:)))); ccd=fft(ccc(40,:); ccd1=[ccd(1:40) zeros(1,length(ccd)-40)]; subplot(224);plot(real(ifft(ccd1))) 结果对比 6.1.4 快速傅立叶变换 计算示例1： x=0:0.1:10; yy=3*cos(2.*x)+6*sin(6.*x); subplot(221);plot(x,yy) subplot(222);plot(x,abs(fft(yy))) subplot(223);plot(x,angle(fft(yy))) subplot(223);plot(x,ifft(fft(yy))) subplot(223);plot(x,angle(fft(yy))) subplot(224);plot(x,ifft(fft(yy))) 示例： 6.2 离散余弦变换 离散余弦矩阵 dctmtx(4) ans = 0.5000 0.5000 0.5000 0.5000 0.6533 0.2706 -0.2706 -0.6533 0.5000 -0.5000 -0.5000 0.5000 0.2706 -0.6533 0.6533 -0.2706 示例2： A = 1 2 3 4 1 2 3 4 2 3 4 5 3 4 5 6 示例2： B=dctmtx(4); cc= B*A*B ans = 13.0000 -4.4609 -0.0000 -0.3170 -3.1543 -0.0000 0.0000 0.0000 1.0000 0.0000 -0.0000 0.0000 0.2242 0.0000 0.0000 -0.0000 示例2 f=zeros(2); ccc=[cc(1:2,1:2) f;f f] ccc = 13.0000 -4.4609