Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (aterrier-plotUpdate) Browse

R8218/measureSCase.m
077e11f382f6aterrier-plotUpdate

Last Change
Raw File

measureSCase.m
View Options

function output = measureScase(varargin)
measureSCaseStart = tic;
Logger.start("measureSCase");
config = getConfig();
% Initialise cases
if nargin > 0
casesToMeasure = [string(varargin)];
else
casesToMeasure = config.casesToMeasure;
end
if config.overwriteMeasurements
SCases = createEmptySCase(casesToMeasure);
else
SCases = loadSCase(casesToMeasure);
end
% Run measurements
for i = 1:length(SCases)
SCaseStart = tic;
Logger.newDelimitedSection(SCases(i).id);
if config.standardMeasurementsToRun.loadData
Logger.logn("");
SCases(i).patient.loadData();
end
if config.shouldersToMeasure.rightAuto
Logger.newDelimitedSection("Shoulder Right Auto");
runShoulderMeasurements(SCases(i).shoulders.right.auto);
Logger.closeSection();
end
if config.shouldersToMeasure.rightManual
Logger.newDelimitedSection("Shoulder Right Manual");
runShoulderMeasurements(SCases(i).shoulders.right.manual);
Logger.closeSection();
end
if config.shouldersToMeasure.leftAuto
Logger.newDelimitedSection("Shoulder Left Auto");
runShoulderMeasurements(SCases(i).shoulders.left.auto);
Logger.closeSection();
end
if config.shouldersToMeasure.leftManual
Logger.newDelimitedSection("Shoulder Left Manual");
runShoulderMeasurements(SCases(i).shoulders.left.manual);
Logger.closeSection();
end
Logger.newDelimitedSection(...
sprintf("SCase measured in %.2fs", toc(SCaseStart)));
Logger.closeBlock();
end
% Save results
if config.saveMeasurements
saveMeasurements(SCases);
end
if config.saveAllMeasurementsInOneFile
saveAllMeasurementsInOneFile(SCases, casesToMeasure);
end
Logger.newDelimitedSection(...
sprintf("Whole measureSCase() ran in %.2fs", toc(measureSCaseStart)));
Logger.stop();
output = SCases;
end
function runShoulderMeasurements(shoulder)
shoulderStart = tic;
config = getConfig();
% Basic measurements
if config.standardMeasurementsToRun.loadData
shoulder.loadData();
end
if shoulder.isempty()
Logger.logn("Not enough data loaded to perform further measurements.")
return
end
if config.standardMeasurementsToRun.morphology
shoulder.morphology();
end
if config.standardMeasurementsToRun.measureFirst
shoulder.measureFirst();
end
if config.standardMeasurementsToRun.measureSecond
shoulder.measureSecond();
end
% Special measurements
Logger.newDelimitedSection("Special measurements"); Logger.logn("");
specialMeasurementNames = string(fields(config.specialMeasurementsToRun))';
for specialMeasurementName = specialMeasurementNames
% Arguments parsing
specialMeasurementArguments =...
config.specialMeasurementsArguments.(specialMeasurementName);
if isempty(specialMeasurementArguments)
specialMeasurementArguments = {};
end
% Recursive call a special measurement
if config.specialMeasurementsToRun.(specialMeasurementName)
try
recursiveMethodCall(shoulder, specialMeasurementName, {shoulder.SCase},...
specialMeasurementArguments{:});
end
end
end
Logger.newDelimitedSection(...
sprintf("Shoulder measured in %.2fs", toc(shoulderStart)));
Logger.closeSection();
end
function saveMeasurements(SCases)
savingStart = tic;
Logger.log("Saving new measurements: ")
try
arrayfun(@(SCase) SCase.saveMatlab, SCases);
Logger.logn("Done in %.2fs", toc(savingStart));
catch ME
Logger.logn(ME.message);
end
end
function saveAllMeasurementsInOneFile(SCases, casesToMeasure)
savingStart = tic;
Logger.log("Saving whole database in one file: ")
try
if not(isequal(casesToMeasure, "*"))
SCaseDB = loadSCase("*");
else
SCaseDB = SCases;
end
save(fullfile(getConfig().dataDir, "matlab", "SCaseDB.mat"), "SCaseDB");
Logger.logn("Done in %.2fs", toc(savingStart));
catch ME
Logger.logn(ME.message);
end
end
c4science · Help