rot.m

图像对齐（inage registration）工具箱

function [I]=rot(a,ord)
% Computing rotational moment invariants based on complex moments from
% image a up to the order ord
 
n=size(a,1);
[M]=kmr(a,ord);  % matrix of complex moments up to the ord order

% p ~ rows

[pom1,pom2]=meshgrid(0:ord,0:ord);
pom3 = (pom1+pom2<=ord).*(pom2>=pom1);
idx = find(pom3==1);


dt = M.*pom3.*(M(1,3).^((pom2-pom1)/2));
I = dt(idx)';