clear
clc

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

%%?Spectral?Residual
myFFT?=?fft2（inImg）;
myLogAmplitude?=?log（abs（myFFT））;
myPhase?=?angle（myFFT）;
mySpectralResidual?=?myLogA