km.m

图像对齐（inage registration）工具箱

function m=km(p,q,w)
% auxilliary function for kmr
% m is a complex moment p,q of image 
% \mu{pq} = M(p+1,q+1)
% w = cm(img,p+q);

m=0;
h1=p+q;

for r=1:p+1;
for s=1:q+1;
	h2=r+s-2;
 	m=m+nchoosek(p,r-1)*nchoosek(q,s-1)*(i^(h1-h2))*((-1)^(q-s+1))*w(h2+1,h1-h2+1);

end;
end;