Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (Bashardoust_SystCoord) Browse

R8218/plotSCase.m
cef7dcff8c61Bashardoust_SystCoord

Last Change
Raw File

plotSCase.m
View Options

function SCase = plotSCase(SCaseId)
%PLOTSCASE Plot a SCase.
% Input: id char of shoulder case (e.g. 'P315')
%
% Output: Corresponding ShoulderCase object
%
% Example: SCase = plotSCase('P315');
%
% Author: Alexandre Terrier, EPFL-LBO
% Creation date: 2018-07-01
% Revision date: 2018-12-30
%
% TODO
% Chech that the SCase exist and can be ploted
%% Open log file
% Check if log dir exists in working directory
logDir = 'log';
if ~exist('log', 'dir') % If not create dir log
mkdir('log')
end
logFile = [logDir '/plotSCase.log'];
logFileID = fopen(logFile, 'w');
fprintf(logFileID, 'plotSCase log file');
fprintf(logFileID, '\nDate: ');
fprintf(logFileID, datestr(datetime));
%% Set the data directory from the configuration file config.txt
configFile = 'config.txt';
if ~exist(configFile, 'file')
error([configFile, ' required']);
end
fileID = fopen(configFile,'r');
dataDir = fscanf(fileID, '%s'); % Read config file without spaces
fclose(fileID);
k = strfind(dataDir, 'dataDir='); % Find definition of dataDir
k = k + 8; % Remove 'dataDir='
dataDir = dataDir(k:end); % Assume that config file ends with dataDir content
% Check if dataDir exists
if ~exist(dataDir, 'dir')
fprintf(logFileID, ['Data directory not found, check ' configFile]);
error(['Data directory not found, check ' configFile]);
end
%% Add path of ShoulderCase classes
addpath('ShoulderCase');
%% Instance of a ShoulderCase object
if (exist('SCase','var') > 0)
clear SCase; % Check for delete method
end
SCase = ShoulderCase; % Instanciate a ShoulderCase object
% Load the SCase
% This should be a method of ShoulderCase class
% The validity of the format should be checked Pnnn or Nnnn.
if (numel(regexp(SCaseId,'^[PN]\d{1,3}$')) == 0)
error('Invalid format of SCaseId argument. SCaseID must start with "P" or "N" and be followed by 1 to 3 digits.');
end
% Directory of SCaseId
SCaseDirLevel0 = SCaseId(1); % Either 'P' or 'N'
strLengthSCaseId = strlength(SCaseId(2:end));
if (strLengthSCaseId < 2)
SCaseDirLevel1 = '0'; % Hunderets
SCaseDirLevel2 = '0'; % Dozent
elseif (strLengthSCaseId < 3)
SCaseDirLevel1 = '0'; % Hunderets
SCaseDirLevel2 = SCaseId(2); % Dozent
else
SCaseDirLevel1 = SCaseId(2); % Hunderets
SCaseDirLevel2 = SCaseId(3); % Dozent
end
% Check if a (!unique! to be done) directory exists for this SCaseId
FindSCaseIdFolder = dir([dataDir '/' SCaseDirLevel0 '/' SCaseDirLevel1 '/' SCaseDirLevel2 '/' SCaseId '*']);
if (isempty(FindSCaseIdFolder)) % No directory for this SCaseId
error(['Missing directory for SCaseId: ' SCaseId]);
end
SCaseDirLevel3 = FindSCaseIdFolder.name;
SCaseDir = [dataDir '/' SCaseDirLevel0 '/' SCaseDirLevel1 '/' SCaseDirLevel2 '/' SCaseDirLevel3];
% Check if this SCaseId has a CT directory with '-1' postfix (preoperative
FindCTFolder = dir([SCaseDir '/' 'CT*-1']);
if (isempty(FindCTFolder)) % No CT directory for this SCaseId
error('Missing CT directory for SCaseId');
% ! should exist script here !
end
CTDir = [SCaseDir '/' FindCTFolder.name];
matlabDir = [CTDir '/matlab']; % we should check for existance
matlabFile = [matlabDir '/SCase.mat'];
if exist(matlabFile, 'file')
% SCaseL = load(matlabFile); % Load SCase dat file
% SCase.id = SCaseL.sCase.id;
% SCase.shoulder = SCaseL.sCase.shoulder;
load(matlabFile);
%SCase=sCase;
SCase.plot; % Plot SCase
else
error('Not matlab file for this SCase');
end
end
c4science · Help