傅里叶变换专题讲座.ppt

a=[1 1 1 1;1 1 1 1;0 0 0 0; 0 0 0 0]; a=repmat(a,64); figure,subplot(1,2,1),imshow(a,[]) subplot(1,2,2),imshow(log(1+abs(fftshift(fft2(a)))),[]) help fftshift a=[1 1 1 1;1 1 1 1;0 0 0 0; 0 0 0 0]; a=repmat(a,64); figure,subplot(1,2,1),imshow(a,[]) subplot(1,2,2),imshow(log(1+abs(fftshift(fft2(a)))),[]) help fftshift a=[1 1 0 0;1 1 0 0;0 0 1 1; 0 0 1 1]; a=repmat(a,64); figure,subplot(1,2,1),imshow(a,[]) subplot(1,2,2),imshow(log(1+abs(fftshift(fft2(a)))),[]) help fftshift I=imread(lenna.bmp); figure,imshow(I,[]) figure,imshow(log(1+fftshift(abs(fft2(I)))),[]) 高频和低频 Saliency Detection: A?Spectral Residual?Approach Saliency Detection: A Spectral Residual Approach %% Read image from file inImg = im2double(rgb2gray(imread(giraffe.jpg))); %% Spectral Residual myFFT = fft2(inImg); myLogAmplitude = log(abs(myFFT)); myPhase = angle(myFFT); mySpectralResidual = myLogAmplitude - imfilter(myLogAmplitude, fspecial(average, 3), replicate); saliencyMap = abs(ifft2(exp(mySpectralResidual + i*myPhase))).^2; %% After Effect saliencyMap = mat2gray(imfilter(saliencyMap, fspecial(gaussian, [3, 3], 1))); figure,imshow(saliencyMap,[]); /blog-493885-589542.html (The PhD Grind) Guo, Chenlei, Qi Ma, and Liming Zhang. “Spatio-temporal saliency detection using phase spectrum of quaternion fourier transform.”?Computer vision and pattern recognition, 2008. cvpr 2008. ieee conference on. IEEE, 2008. The results are good, but the reason is questionable The no-end story %% Read image from file inImg = im2double(rgb2gray(imread(giraffe.jpg))); %% Spectral Residual myFFT = fft2(inImg); myPhase = angle(myFFT); saliencyMap = abs(ifft2(exp(i*myPhase))).^2; %% After Effect saliencyMap = mat2gray(imfilter(saliencyMap, fspecial(gaussian, [3, 3], 1))); figure,imshow(saliencyMap,[]); * Fourier变换 What does the Fourier Transform do? Given a smoothie, it finds the recipe. What does the Fourier transform do? Given a smoothie, it finds the recipe. How? Run the smoothie through filt