function [imageAveFFT,phase,sps] = TransformImages(image1Comp,image2Comp);
image1FFT= fftshift(fft2(image1Comp));
image2FFT= fftshift(fft2(image2Comp));
imageAveFFT= (image1FFT + image2FFT)/2;
phase= angle(imageAveFFT);
aveSPS= (image1FFT.*conj(image2FFT) + image2FFT.*conj(image1FFT))/2;
sps= max(0,aveSPS);