Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (aterrier-plotUpdate) Browse

R8218/degenerationAll.m
da3e00881c47aterrier-plotUpdate

Last Change
Raw File

degenerationAll.m
View Options

function degenerationAll(caseID, musclesList)
%%DEGENERATIONALL fills the Shoulder database with muscle degeneration measurements
%
% This function fills the Shoulder database with the muscle degeneration
% values of the case given as input. The measurements are stored in a XLS
% file.
%
% USE: degenerationAll(caseID, musclesList)
%
% INPUT caseID: Case IDs for which you want to compute muscle degeneration.
% Cell array that stores the case IDs as a char in the form 'P###' or
% 'N###' (starts with "P" or "N" followed by 1 to 3 digits).
%
% musclesList: List of muscles which you want to compute muscle
% degeneration. Cell array that stores each msucle name as a char.
%
% REMARKS The muscle values are saved in the muscles.xls file.
%
% created with MATLAB ver.: 9.2.0.556344 (R2017a) on Windows 7
% Author: EPFL-LBO-VMC
% Date: 12-Jun-2017
%
CTDatabaseLocation = 'Z://data'; % Location of the CT database
XLSShoulderDatabaseLocation = 'C:/Users/vmalfroy/Dropbox/shoulderDatabase/'; % Location of the XLS ShoulderDatabase
HUThresholdMuscle = -29; %#ok<NASGU>
HUThresholdOsteochondroma = 166.25; %#ok<NASGU>
% Check validity of input argument (cell array)
validateattributes(caseID,{'cell'},{'nonempty'});
% Check that imput case IDs are unique
[caseIDNames,~,uniqueCaseID] = unique(caseID);
countOfCaseID = hist(uniqueCaseID,unique(uniqueCaseID))';
freqCaseIDs = struct('name',caseIDNames','freq',num2cell(countOfCaseID));
repeatedValues = find([freqCaseIDs.freq] > 1);
if(~isempty(repeatedValues))
warning('Input contains repeated Case ID: %s. Repeated occurences of that Case ID will be omitted. \n', freqCaseIDs(repeatedValues).name)
caseID = caseIDNames';
end
% open mySQL connection
conn = openSQL();
% Get the list of exisiting cases in XLS database
[~,~,currDatabase] = xlsread([XLSShoulderDatabaseLocation 'Excel/xlsFromMatlab/muscles.xls']);
currDatabaseCaseIDlist = currDatabase(2:end,1);
% Compute muscle degeneration values and fill database
for i = 1:length(caseID)
% Check validity of input arguments (char and format 'P###' or 'N###')
validateattributes(caseID{i},{'char'},{'nonempty'});
if (numel(regexp(caseID{i},'^[PN]\d{1,3}$')) == 0)
error(['Invalid format of CaseID: ' caseID{i} '. CaseID must start with "P" or "N" and be followed by 1 to 3 digits.']);
end
% Find CT folder for the given Case ID
levelDir1 = caseID{i}(1);
if (length(caseID{i}(2:end)) < 2)
levelDir2 = '0';
levelDir3 = '0';
elseif (length(caseID{i}(2:end)) < 3)
levelDir2 = '0';
levelDir3 = caseID{i}(2);
else
levelDir2 = caseID{i}(2);
levelDir3 = caseID{i}(3);
end
FindCaseCTFolder = dir([CTDatabaseLocation '/' levelDir1 '/' levelDir2 '/' levelDir3 '/' caseID{i} '*']);
if (isempty(FindCaseCTFolder))
error(['Missing CT directory for CaseID: ' caseID{i}]);
end
CaseID_IPP = FindCaseCTFolder.name;
CaseCTFolder = [CTDatabaseLocation '/' levelDir1 '/' levelDir2 '/' levelDir3 '/' CaseID_IPP];
% Compute musle degeneration values and store in the XLS database
if exist([CaseCTFolder '/CT-' CaseID_IPP '-1/muscles'],'dir') == 7 %if a folder "muscles" exist, the patient name is added to the "measured" column
muscleDirectory = [CaseCTFolder '/CT-' CaseID_IPP '-1/muscles']; %#ok<NASGU>
muscles_results = struct('CT_id', caseID{i});
for j=1:length(musclesList)
% Compute muscle degeneration values
eval([musclesList{j} '= muscleDegeneration(caseID{i}, musclesList{j},HUThresholdMuscle,HUThresholdOsteochondroma,muscleDirectory);']);
eval(['muscles_results = setfield(muscles_results, musclesList{j}, ' musclesList{j} ');']);
end
muscles = cellstr(muscles_results.CT_id)';
for j = 1:length(musclesList)
eval(['muscles = horzcat(muscles, transpose(struct2cell(muscles_results.' musclesList{j} ')));']);
end
% Write to XLS file
% Replace line if already existing or append
if(any(strcmp(currDatabaseCaseIDlist, caseID{i})))
xlswrite([XLSShoulderDatabaseLocation 'Excel/xlsFromMatlab/muscles.xls'],muscles,'Feuil1',['A' int2str(find(strcmp(currDatabaseCaseIDlist, caseID{i}))+1)]);
else
xlswrite([XLSShoulderDatabaseLocation 'Excel/xlsFromMatlab/muscles.xls'],muscles,'Feuil1',['A' int2str(length(currDatabaseCaseIDlist)+1+i)]);
end
% Write to MySQL database
% Replace line if already existing or issue warning message
sqlquery = ['SELECT CT_id FROM CT WHERE shoulder_id IN (SELECT shoulder_id FROM sCase WHERE folder_name = "' caseID{i} '")'];
curs = exec(conn,sqlquery);
if ~isempty(curs.Message)
fprintf('\nMessage: %s\n', curs.Message);
end
curs=fetch(curs);
sqlresult = curs.Data;
if isnumeric(sqlresult{1})
% Update data in table 'muscle'
data = {sqlresult{1} muscles_results.SS.atrophy muscles_results.SS.infiltration muscles_results.SS.osteochondroma muscles_results.SS.degeneration muscles_results.IS.atrophy muscles_results.IS.infiltration muscles_results.IS.osteochondroma muscles_results.IS.degeneration muscles_results.SC.atrophy muscles_results.SC.infiltration muscles_results.SC.osteochondroma muscles_results.SC.degeneration muscles_results.TM.atrophy muscles_results.TM.infiltration muscles_results.TM.osteochondroma muscles_results.TM.degeneration};
fieldNames = {'CT_id','SSA','SSI','SSO','SSD','ISA','ISI','ISO','ISD','SCA','SCI','SCO','SCD', 'TMA', 'TMI', 'TMO', 'TMD'};
upsert(conn, 'muscle', fieldNames, [1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0], data);
else
error(['No entry for caseID ' caseID{i} ' was found in the ' conn.Instance ' MySQL database.']);
end
else
error(['"muscles" folder missing for CaseID: ' caseID{i}]);
end
end
% Close mySQL database
closeSQL(conn);
end
c4science · Help