measureWidthAndHeightWithPCSA.m
8170d98ab94fmaster

measureWidthAndHeightWithPCSA.m
View Options

			function measureWidthAndHeightWithPCSA(obj)
			% Calculation of the PCA coeficient (vectors heights and width)
			[coef,~]=pca(obj.surface.points);

			% Base matrix in the "PCA" referential, in the order
			% height-width-depth
			px = [];% empty vector used to store the dot products between the PCA axes and the scapula axes
			py = [];
			pz = [];
			for i = 1:3
			px = [px abs(dot(coef(:,i),xAxis))];
			py = [py abs(dot(coef(:,i),yAxis))];
			pz = [pz abs(dot(coef(:,i),zAxis))];
			end
			[~,a] = max(px); % Here, the width is mainly on the xAxis,
			[~,b] = max(py); % height on the yAxis
			[~,c] = max(pz); % and depth on the zAxis

			if a==b && a==c
			a=1;b=2;c=3;
			else
			if a==b
			if max(py)>max(px)
			a=6-b-c;
			else
			b=6-a-c;
			end
			end
			if a==c
			if max(px)>max(pz)
			c=6-a-b;
			else
			a=6-b-c;
			end
			end
			if b==c
			if max(py)>max(pz)
			c=6-a-b;
			else
			b=6-a-c;
			end
			end
			end

			NewPCAAxis = [coef(:,b),coef(:,a),coef(:,c)]; % New basis in PCA referential, columns are then in the order "height-width-depth"
			PCAGlenSurf = [];
			for i=1:size(obj.surface.points(:,1))
			vp = NewPCAAxis^-1 * obj.surface.points(i,:)';
			PCAGlenSurf = [PCAGlenSurf;vp'];
			end

			heightPCA = max(PCAGlenSurf(:,1))-min(PCAGlenSurf(:,1));
			widthPCA = max(PCAGlenSurf(:,2))-min(PCAGlenSurf(:,2));

			obj.height=heightPCA;
			obj.width=widthPCA;
			end