Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (aterrier-plotUpdate) Browse

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

Last Change
Raw File

Humerus.m
View Options

classdef Humerus < handle
%HUMERUS Properties and methods associated to the humerus
% Detailed explanation goes here
% Author: Alexandre Terrier, EPFL-LBO
% Creation date: 2018-07-01
% Revision date: 2019-06-29
%
% TO DO:
% Local coordinate system
properties
landmarks % 5 3D points
insertions
center % Center of the humeral head (sphere fit on 5 points
radius % Radius of the humeral head (sphere fit on 5 points
jointRadius % Radius of cuvature of the articular surface (todo)
SHSAngle % Scapulo-humeral subluxation angle
SHSPA % Scapulo-humeral subluxation angle in the postero-anterior direction (posterior is negative, as for glenoid version)
SHSIS % Scapulo-humeral subluxation angle in the infero-superior direction
SHSAmpl % Scapulo-humeral subluxation (center ofset / radius)
SHSOrient % Scapulo-humral subluxation orientation (0 degree is posterior)
GHSAmpl % Gleno-humeral subluxation (center ofset / radius)
GHSOrient % Gleno-humral subluxation orientation (0 degree is posterior)
subluxationIndex3D
shoulder
end
methods (Access = ?Shoulder) % Only Shoulder is allowed to construct Humerus
function obj = Humerus(shoulder)
%HUMERUS Construct an instance of this class
% Detailed explanation goes here
obj.landmarks = [];
obj.center = [];
obj.radius = [];
obj.jointRadius = [];
obj.SHSAngle = [];
obj.SHSPA = [];
obj.SHSIS = [];
obj.SHSAmpl = [];
obj.SHSOrient = [];
obj.GHSAmpl = [];
obj.GHSOrient = [];
obj.shoulder = shoulder;
end
end
methods
function output = loadData(obj)
% Call methods that can be run after the ShoulderCase object has
% been instanciated.
SCase = obj.shoulder.SCase;
if exist(...
fullfile(SCase.dataSlicerPath, "HH_landmarks_"+obj.shoulder.side+".mrk.json"),...
"file")
success = Logger.timeLogExecution(...
"Humerus load landmarks (slicer): ",...
@(obj) obj.loadSlicerLandmarks, obj);
elseif not(isempty(dir(fullfile(SCase.dataAmiraPath, "*HH*.landmarksAscii"))))
success = Logger.timeLogExecution(...
"Humerus load landmarks (slicer): ",...
@(obj) obj.loadSlicerLandmarks, obj);
else
success = Logger.timeLogExecution(...
"Humerus load landmarks: ",...
@(message) error(message), "No landmarks file found");
end
output = success;
end
function output = morphology(obj)
% Call methods that can be run after loadData() methods has been run
% by all ShoulderCase objects.
success = Logger.timeLogExecution(...
"Humerus center and radius: ",...
@(obj) obj.measureCenterAndRadius, obj);
output = success;
end
function output = measureFirst(obj)
% Call methods that can be run after morphology() methods has been run
% by all ShoulderCase objects.
success = Logger.timeLogExecution(...
"Humerus scapulo-humeral subluxation: ",...
@(obj) obj.measureScapuloHumeralSubluxation, obj);
output = success;
end
function output = measureSecond(obj)
% Call methods that can be run after measureFirst() methods has been run
% by all ShoulderCase objects.
success = Logger.timeLogExecution(...
"Humerus gleno-humeral subluxation: ",...
@(obj) obj.measureGlenoHumeralSubluxation, obj);
success = success & Logger.timeLogExecution(...
"Humerus 3D subluxation index: ",...
@(obj) obj.measureSubluxationIndex3D, obj);
output = success;
end
end
end
c4science · Help