sigprocexperts.com Git - sie_pcm/blob - run_example.m

   1 % load images
   2 a = imread('img1.png');
   3 b = imread('img2.png');
   4
   5 % transform image data to Fourier domain
   6 A = fftshift(fft2(double(a)));
   7 B = fftshift(fft2(double(b)));
   8
   9 % call pcm estimation code
  10 [x,y]=siepcm(A,B,1);                    % 3rd arg == 1: display debug images
  11