kmr.m

图像对齐（inage registration）工具箱

function [M]=kmr(img,p)

% M is a matrix of complex moments up to the p-th order of the image img
% img - image matrix
% \mu{pq} = M(p+1,q+1)

w=cm(img,2*p);       

[n1,n2]=size(img);

for r=1:p+1;
for s=1:p+2-r;
M(r,s)=km(r-1,s-1,w);
end;
end;