基于MATLAB的图像变换技术.docVIP

傅立叶级数变换： i=zeros(35,35); f(5:24,13:17)=1; f=fft2(f,256,256); f2=log(abs(f)); Warning: Log of zero. imshow(f2,[-1,5],notruesize); colormap(jet); colorbar; i=zeros(256,256); i(124:132,120:136)=1; subplot(1,3,1);imshow(i); i=im2double(i); f1=fft2(i); subplot(1,3,2);imshow(log(1+abs(f1)),[]); fc1=fftshift(f1); subplot(1,3,3);imshow(fc1); Warning: Displaying real part of complex input In imshowValidateCdata at 497 In imshowParseInputs at 378 In imshow at 85 i1=imread(rice.png); i1=im2double(i1); f1=fft2(i1); fc1=fftshift(f1); subplot(1,2,1);imshow(i1); subplot(1,2,2);imshow(log(1+abs(fc1))); 求线性滤波器的频率响应： h=fspecial(gaussian); freqz2(h); f=zeros(1000,1000); f(350:649,475:524)=1; subplot(2,2,1);imshow(f,notruesize); subplot(2,2,2);f=fftshift(abs(fft2(f))); imshow(f,[-1,5],notruesize); f=zeros(1000,1000); f(350:649,475:524)=1; f=imrotate(f,45,bilinear,crop); subplot(2,2,3);imshow(f,notruesize); f=fftshift(abs(fft2(f))); subplot(2,2,4);imshow(f,[-1,5],notruesize); a=magic(3); b=ones(3); a(8,8)=0; b(8,8)=0; c=fft2(a).*fft2(b); d=ifft2(c); d=d(1:5,1:5); d=real(d) d = 8.0000 9.0000 15.0000 7.0000 6.0000 11.0000 17.0000 30.0000 19.0000 13.0000 15.0000 30.0000 45.0000 30.0000 15.0000 7.0000 21.0000 30.0000 23.0000 9.0000 4.0000 13.0000 15.0000 11.0000 2.0000 离散余弦变换： rgb=imread(autumn.tif); i=rgb2gray(rgb); j=dct2(i); subplot(1,2,1);imshow(rgb); subplot(1,2,2);imshow(log(abs(j)),[]); colormap(jet(64));colorbar; rgb=imread(autumn.tif); i=rgb2gray(rgb); j=dct2(i); subplot(1,2,1);imshow(rgb); subplot(1,2,2);imshow(log(abs(j)),[]); colormap(jet(64));colorbar; clear; rgb=imread(autumn.tif); i=rgb2gray(rgb); j=dct2(i); subplot(1,2,1);imshow(rgb); j(abs(j)10)=0; k=idct2(j); subplot(1,2,2);imshow(k,[0,255]); i=imread(rice.png); a=im2double(i); subplot(1,2,1);imshow(a); d=dctmtx(size(a,1)); dct=d*a*d; subplot(1,2,2);imshow(dct); Walsh--Hadamard变换： sq=[1