diff --git a/ShoulderCase/@Glenoid/measureCenter.m b/ShoulderCase/@Glenoid/measureCenter.m index 2926319..b7f7706 100644 --- a/ShoulderCase/@Glenoid/measureCenter.m +++ b/ShoulderCase/@Glenoid/measureCenter.m @@ -1,7 +1,7 @@ function measureCenter(obj) - % Glenoid center is the closet point between the mean of the extracted - % surface and all the points of the extracted surface. +% Glenoid center is the point of the segmented glenoid's surface which is the +% closest to the mean of all the points of the segmented glenoid's surface surfacePointsFromMeanPoint = obj.surface.points - obj.surface.meanPoint; centerIndex = findShortest3DVector(surfacePointsFromMeanPoint); obj.center = obj.surface.points(centerIndex,:); end diff --git a/ShoulderCase/@Glenoid/measureVersionAndInclination.m b/ShoulderCase/@Glenoid/measureVersionAndInclination.m index 2998195..cf62099 100644 --- a/ShoulderCase/@Glenoid/measureVersionAndInclination.m +++ b/ShoulderCase/@Glenoid/measureVersionAndInclination.m @@ -1,34 +1,34 @@ function measureVersionAndInclination(obj) % Glenoid version 3D, defined by the amplitude and orientation of the angle % between the glenoid centerline and the scapular transverse axis (z axis). % The orientation is the angle between the -x axis and the glenoid % centerline projected in the xy plane. Zero orientation correspond to -x % axis (posterior side), 90 to superior, 180 to anterior orientaion, and % -90 to inferior. Glenoid version is also reported as version2D (>0 for % retroversion) and inclination (>0 for superior tilt). absoluteCenterLine = obj.centerLine + obj.scapula.coordSys.origin; localCenterLine = obj.scapula.coordSys.express(absoluteCenterLine); X = [1 0 0]; Y = [0 1 0]; Z = [0 0 1]; - versionAmplitudeRotation = vrrotvec(Z,localCenterLine); % angle between centerLine and scapular axis + versionAmplitudeRotation = vrrotvec(Z,localCenterLine); % angle between centerLine and scapular medio-lateral axis obj.versionAmplitude = rad2deg(versionAmplitudeRotation(4)); XYProjectedCenterLine = obj.scapula.coordSys.projectOnXYPlane(absoluteCenterLine); localXYProjectedCenterLine = obj.scapula.coordSys.express(XYProjectedCenterLine); versionOrientationRotation = vrrotvec(-X,localXYProjectedCenterLine); obj.versionOrientation = sign(localXYProjectedCenterLine*Y') * rad2deg(versionOrientationRotation(4)); ZXProjectedCenterLine = obj.scapula.coordSys.projectOnZXPlane(absoluteCenterLine); localZXProjectedCenterLine = obj.scapula.coordSys.express(ZXProjectedCenterLine); versionRotation = vrrotvec(Z,localZXProjectedCenterLine); obj.version = sign(localZXProjectedCenterLine*X') * rad2deg(versionRotation(4)); YZProjectedCenterLine = obj.scapula.coordSys.projectOnYZPlane(absoluteCenterLine); localYZProjectedCenterLine = obj.scapula.coordSys.express(YZProjectedCenterLine); inclinationRotation = vrrotvec(Z,localYZProjectedCenterLine); obj.inclination = sign(localYZProjectedCenterLine*Y') * rad2deg(inclinationRotation(4)); end