hgicf2im.m

图像对齐（inage registration）工具箱

function [H,C]=hgicf2im(X,k)
% HGICF2IM - computes how good is constant fit to image
%               [H,C]=hgicf2im(X,k)
%
% X: image
% C: averages over squares size 2k+1
% H: fit goodness over squares
% k: half size of square to fit locally    
%     
% Copyright Zitova B., Kautsky J., Peters G., Flusser J.
% UTIA, 1998

N=-k:k;
[m,n]=size(X);
C=zeros(m-2*k,n-2*k);
H=C;
y=k+1:n-k; x=k+1:m-k;

for jx=N
   for jy=N
      C=C+X(x+jx,y+jy);
   end
end
C=C/(2*k+1)^2;
for jx=N
   for jy=N
      H=H+(X(x+jx,y+jy)-C).^2;
   end
end