mnc.m

图像对齐（inage registration）工具箱

function [odp,s,uhel]=mnc(pary)
%computes transformation parameters by means of LSM
% pary - pairs of coordinates of mutually corresponding CPs

body1=pary(:,[1 2]);
body2=pary(:,[3 4]);

sumax1=sum(body1(:,1));
sumay1=sum(body1(:,2));
sumaxy1=sum(body1(:,1).^2+body1(:,2).^2);

sumax2=sum(body2(:,1));
sumay2=sum(body2(:,2));
sumaxy21=sum(body1(:,1).*body2(:,1)+body1(:,2).*(body2(:,2)));
sumaxy22=sum(body1(:,1).*body2(:,2)-body1(:,2).*(body2(:,1)));

n=size(body1(:,1));
matice=[sumaxy1 0 sumax1 sumay1;0 sumaxy1 -1*sumay1 sumax1;sumax1 -1*sumay1 n(1) 0; sumay1 sumax1 0 n(1)];
matice2=[sumaxy21; sumaxy22; sumax2; sumay2];

odp=matice\matice2;

s=sqrt((odp(1)^2+odp(2)^2));
muhel(1)=(acos(odp(1)/s));
muhel(2)=(asin(odp(2)/s));

znameni=sign(muhel(2))<0;
uhel=(~znameni)*muhel(1)+(znameni)*(2*pi-muhel(1));


uhel=uhel/(pi/180);