显著性检测

 1 clear
 2 clc
 3 
 4 %% Read image from file
 5 inImg = im2double(rgb2gray(imread(‘256.png‘)));
 6 %%inImg = imresize(inImg, 64/size(inImg, 2));
 7 
 8 %% Spectral Residual
 9 myFFT = fft2(inImg);
10 myLogAmplitude = log(abs(myFFT));
11 myPhase = angle(myFFT);
12 mySpectralResidual = myLogAmplitude - imfilter(myLogAmplitude, fspecial(‘average‘, 3), ‘replicate‘);
13 saliencyMap = abs(ifft2(exp(mySpectralResidual + i*myPhase))).^2;
14 
15 %% After Effect
16 saliencyMap = mat2gray(imfilter(saliencyMap, fspecial(‘gaussian‘, [10, 10], 2.5)));
17 imshow(saliencyMap);