diff --git a/ShoulderCase/@Glenoid/measureBetaAngle.m b/ShoulderCase/@Glenoid/measureBetaAngle.m index 9a461d7..d2eb517 100644 --- a/ShoulderCase/@Glenoid/measureBetaAngle.m +++ b/ShoulderCase/@Glenoid/measureBetaAngle.m @@ -1,18 +1,18 @@ function measureBetaAngle(obj) % Glenoid BETA angle is the angle between the floor of the supraspinatus fossa % marked by a sclerotic line and the glenoid fossa line. The floor of the % supraspinatus fossa is the line fitted to the scapular groove points. % The glenoid fossa line is the glenoid supero-inferior line. - supraspinatusFossaFloor = Vector(fitLine(obj.scapula.groove)') + supraspinatusFossaFloor = Vector(fitLine(obj.scapula.groove)'); supraspinatusFossaFloor = supraspinatusFossaFloor.orientToward(obj.scapula.coordSys.ML); betaRotation = vrrotvec(supraspinatusFossaFloor.vector, -obj.inferoSuperiorLine.vector); % Unlike glenoid retroversion angle, glenoid beta angle method doesn't set % the sign of this angle which should be compared to the rotation axis alignment % with scapular PA axis. However, this alignement depends on the shoulder side. % It's assumed that glenoid supero-inferior line will always be aligned with % scapular SI axis, and supraspinatusFossaFloor line will always be aligned % with scapular ML axis. Therefore glenoid beta angle will always be a % positive number. obj.beta = rad2deg(betaRotation(4)); end \ No newline at end of file diff --git a/ShoulderCase/@Humerus/measureSubluxationIndex3D.m b/ShoulderCase/@Humerus/measureSubluxationIndex3D.m index 822979a..0859833 100644 --- a/ShoulderCase/@Humerus/measureSubluxationIndex3D.m +++ b/ShoulderCase/@Humerus/measureSubluxationIndex3D.m @@ -1,15 +1,15 @@ -function measureSubluxationindex3D(obj) +function measureSubluxationIndex3D(obj) % - GL: (Glenoid Line) Line that joins the most anterior glenoid rim point and % the most posterior glenoid rim point. % - HHv: (Humeral Head vector) Vector from the glenoid center to the humeral % head’s fitted sphere center. % - HHp: (Humeral Head projection) Projection of HHv on GL. % - HHr: (Humeral Head radius) Humeral head’s fitted sphere radius. % % The 3D humeral subluxation index is the ratio of ||HHp|| + ||HHr|| over 2*||HHr||. GL = obj.shoulder.scapula.glenoid.posteroAnteriorLine; HHv = Vector(obj.shoulder.scapula.glenoid.center, obj.center); HHp = GL.project(HHv); HHr = obj.radius; obj.subluxationIndex3D = (HHp.norm + HHr) / (2*HHr); end \ No newline at end of file