Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (aterrier-plotUpdate) Browse

R8218/ShoulderCase/@Humerus/measureScapuloHumeralSubluxation.m
d32ee1bcc847aterrier-plotUpdate

Last Change
Raw File

measureScapuloHumeralSubluxation.m
View Options

function measureScapuloHumeralSubluxation(obj)
% Scapulo-humeral subluxation (SHS)
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);
glenoidToHH = Vector(scapula.glenoid.center, obj.center);
% SHS amplitude (ratio rather than amplitude)
% glenoidToHHxy is the vector from humeral head center to scapular (z) axis,
% perpendicular to z-axis --> similar to projection in xy plane
glenoidToHHxy = zAxis.orthogonalComplementTo(glenoidToHH);
obj.SHSAmpl = norm(glenoidToHHxy) / obj.radius;
% SHS orientation
% Convention that superior is positive
obj.SHSOrient = rad2deg(angle(-xAxis, glenoidToHHxy));
obj.SHSOrient = sign(dot(yAxis, glenoidToHHxy)) * obj.SHSOrient;
% SHSAngle: angle between glenHumHead and zAxis
obj.SHSAngle = rad2deg(angle(glenoidToHH, zAxis));
% SHSPA is angle between zxProjection and zAxis
% Convention that posterior is negative
glenoidToHHzx = yAxis.orthogonalComplementTo(glenoidToHH);
obj.SHSPA = rad2deg(angle(glenoidToHHzx, zAxis));
obj.SHSPA = sign(dot(xAxis, glenoidToHHzx)) * obj.SHSPA;
% SHSIS is angle between yzProjection and zAxis
% Convention that superior is negative
glenoidToHHyz = xAxis.orthogonalComplementTo(glenoidToHH);
obj.SHSIS = rad2deg(angle(glenoidToHHyz, zAxis));
obj.SHSIS = sign(dot(yAxis, glenoidToHHyz)) * obj.SHSIS;
end
c4science · Help