clear
clc

%%?Read?image?from?file
inImg?=?im2double（rgb2gray（imread（‘sheep.jpg‘）））;
%%inImg?=?imresize（inImg?64/size（inImg?2））;

%%?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‘?[10?10]?2.5）））;
imshow（saliencyMap）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????359040??2013-04-03?18:10??SR\sheep.jpg

?????文件????????536??2014-05-02?19:16??SR\SR.m

?????目錄??????????0??2014-05-02?19:27??SR

-----------?---------??----------?-----??----

???????????????359576????????????????????3