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)';