function metric=MergeMetrics(metric1,metric2) M1xx=metric1(:,1); M1xy=metric1(:,2); M1yy=metric1(:,3); M1lambda1=0.5*((M1xx+M1yy)+sqrt(4*M1xy.^2+(M1xx-M1yy).^2)); M1lambda2=0.5*((M1xx+M1yy)-sqrt(4*M1xy.^2+(M1xx-M1yy).^2)); lambda1=min(M1lambda1,M1lambda2); M2xx=metric1(:,1); M2xy=metric1(:,2); M2yy=metric1(:,3); M2lambda1=0.5*((M2xx+M2yy)+sqrt(4*M2xy.^2+(M2xx-M2yy).^2)); M2lambda2=0.5*((M2xx+M2yy)-sqrt(4*M2xy.^2+(M2xx-M2yy).^2)); lambda2=min(M2lambda2,M2lambda2); metric=metric1; pos=find(lambda2