measureGlenoHumeralSubluxation.m
0186d859aeddaterrier-plotUpdate

measureGlenoHumeralSubluxation.m
View Options

			function output = measureGlenoHumeralSubluxation(obj)
			% Gleno-humeral subluxation (GHS)

			scapula = obj.shoulder.scapula;

			xAxis = Vector(...
			scapula.coordSys.origin,...
			scapula.coordSys.origin + scapula.coordSys.PA);
			yAxis = Vector(...
			scapula.coordSys.origin,...
			scapula.coordSys.origin + scapula.coordSys.IS);
			zAxis = Vector(...
			scapula.coordSys.origin,...
			scapula.coordSys.origin + scapula.coordSys.ML);

			glenoidCenterLine = Vector(...
			scapula.glenoid.center,...
			scapula.glenoid.center + scapula.glenoid.centerLine);

			% Vector from glenoid sphere centre to glenoid centerline (perpendicular)
			GHS = glenoidCenterLine.orthogonalComplementTo(obj.center);

			% GHS amplitude (ratio rather than amplitude)
			obj.GHSAmpl = norm(GHS) / obj.radius;

			% GHS orientation
			% Superior is positive
			GHSxy = zAxis.orthogonalComplementTo(GHS);
			obj.GHSOrient = rad2deg(angle(-xAxis, GHSxy));
			obj.GHSOrient = sign(dot(yAxis, GHSxy)) * obj.GHSOrient;
			end