clear
clc

%%?Read?image?from?file
inImg?=?im2double（rgb2gray（imread（‘D:\matlab\test\image\41.jpg‘）））;
%%inImg?=?imresize（inImg?64/size（inImg?2））;

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