% load images
a = imread('img1.png');
b = imread('img2.png');

% transform image data to Fourier domain
A = fftshift(fft2(double(a)));
B = fftshift(fft2(double(b)));

% call pcm estimation code
[x,y]=siepcm(A,B,1);                    % 3rd arg == 1: display debug images

