transforming_mesh_points.m
4730834938aemaster

transforming_mesh_points.m
View Options

			% R_H2_T
			Zh = BLDATA.Original_Points.GH - BLDATA.Original_Points.HU; Zh = Zh/norm(Zh);
			Yh = cross(BLDATA.Original_Points.GH - BLDATA.Original_Points.EL, BLDATA.Original_Points.EL - BLDATA.Original_Points.EM); Yh = Yh/norm(Yh);
			Xh = cross(Yh, Zh); Xh = Xh/norm(Xh);
			R_H2_T = [Xh, Yh, Zh]; % Local -> absolute


			% R_H1_T
			Zh = BLDATA.Original_Points.HH - BLDATA.Original_Points.HU; Zh = Zh/norm(Zh);
			Yh = cross(BLDATA.Original_Points.HH - BLDATA.Original_Points.EL, BLDATA.Original_Points.EL - BLDATA.Original_Points.EM); Yh = Yh/norm(Yh);
			Xh = cross(Yh, Zh); Xh = Xh/norm(Xh);
			R_H1_T = [Xh, Yh, Zh]; % Local -> absolute


			% transform Humerus mesh points to H2 coordinate system
			for i = 1:length(Humerus_Mesh0.points)
			Humerus_Mesh_H2(:, i) = 1e-3(R_H2_T'(BLDATA.Original_Points.HH + R_H1_T1e3Humerus_Mesh0.points(i, :)' - BLDATA.Original_Points.GH));
			end