ami.m

图像对齐（inage registration）工具箱

function [I]=ami(a)
% Computing six AMIs of an image stored in matrix a
 
n=size(a,1);
m=cm(a,4);


 m00=m(1,1);
 m11=m(2,2);
 m02=m(1,3);
 m12=m(2,3);
 m03=m(1,4);
 m21=m(3,2);
 m30=m(4,1);
 m20=m(3,1);
 m13=m(2,4);
 m04=m(1,5);
 m31=m(4,2);
 m40=m(5,1);
 m22=m(3,3);
 
 
I(1)=(m20*m02-m11^2)/m00^4;
I(2)=(m30*m30*m03*m03-6*m30*m21*m12*m03+4*m30*m12^3+4*m03*m21^3-3*(m21*m12)^2)/m00^10;
I(3)=(m20*(m21*m03-m12*m12)-m11*(m30*m03-m21*m12)+m02*(m30*m12-m21*m21))/m00^7;
I(4)=(m20^3*m03^2-6*m20*m20*m11*m12*m03-6*m20*m20*m02*m21*m03+9*m20*m20*m02*m12*m12);
I(4)=I(4)  +12*m20*m11*m11*m21*m03+6*m20*m11*m02*m30*m03-18*m20*m11*m02*m21*m12;
I(4)=I(4)   -8*m11*m11*m11*m30*m03-6*m20*m02*m02*m30*m12+9*m20*m02*m02*m21*m21;
I(4)=(I(4)  +12*m11*m11*m02*m30*m12-6*m11*m02*m02*m30*m21+m02^3*m30^2)/m00^11;
I(5)=(m40*m04-4*m31*m13+3*m22^2)/m00^6;
I(6)=(m40*m04*m22+2*m31*m13*m22-m40*m13^2-m04*m31^2-m22^3)/m00^9;