


Xov = OVA.OV(i).Xov;
Yov = OVA.OV(i).Yov;
%[ Xov Yov]

Xov = Xov-mean(Xov);
Yov = Yov-mean(Yov);

C = cov(Xov,Yov);
[L,D] = eig(C);

theta0 = asind(L(1,1));
aspect0 = sqrt(D(1,1)/D(2,2));
