Image Offset from Matched Filter

image file: out2.img
image size: 512 x 512
block size: 5
threshold: 0.05
max # peaks: 20

Searching for peaks
peak	   value	hoff	voff
   1	1.000000	  62	 362
-1
   2	0.945098	 511	 350
   3	0.925490	 423	 306
   4	0.925490	 423	 322
   5	0.913725	   0	 190
   6	0.909804	   0	 274
   7	0.905882	 511	 151
   8	0.901961	 424	  88
   9	0.890196	 424	  51
  10	0.886275	   0	 151

>> [m1,map]=tiffread('map1.tif');
>> m1 = ind2gray(m1,map);
>> [m2,map]=tiffread('smap2.tif');
>> m2 = ind2gray[m2,map];

>> sect = m2(1:64,1:64);
>> whos

  Name       Size         Bytes  Class

  m1       583x424      1977536  double array
  m2       545x424      1848640  double array
  map      256x3           6144  double array
  sect      64x64         32768  double array

>> m1c=zeros(512,512);
>> m1c(1:512,1:424) = m1(72:583,1:424);
>> imshow(m1c);

>> fft1= fft2(fftshift(m1c));
>> sectfft = fft2(sect,512,512);
>> prod = fft1.*conj(sectfft);
>> out = abs(ifft2(prod));
>> out = out/max(out(:));
>> out = fftshift(out);
>> imshow(out);

>> bpof = imag(sectfft)>0;
>> bpof = 2*bpof-1;
>> out2 = abs(fftshift(ifft2(fft1.*bpof)));
>> out2 = out2/max(out2(:));
>> imshow(out2);
>> tiffwrite(sect,'sect.tif');
>> imshow(bpof);
>> tiffwrite((bpof+1)/2.0,'bpof.tif');
>> tiffwrite(out2,'out2.tif');
>> tiffwrite(out,'out.tif');
>> tiffwrite(m1c,'m1c.tif');