Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (Bashardoust_SystCoord) Browse

R8218/ShoulderCase/@Glenoid/measureVersionAndInclination.m
b55b6a89fc1aBashardoust_SystCoord

Last Change
Raw File

measureVersionAndInclination.m
View Options

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
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
c4science · Help