Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (aterrier-plotUpdate) Browse

R8218/muscleContours.m
7e5facd3688daterrier-plotUpdate

Last Change
Raw File

muscleContours.m
View Options

%% muscleContours
% This script allows to select the contours on Tiff images of the four
% cuff rotator muscles: Supraspinatus (SS), Infraspinatus (IS), Subscapularis (SC)
% and Teres Minor (TM).
%
% INPUT:
% subject: ### subject number
% type: 'normal' or 'pathologic'
% Author: EPFL-LBO
% Date: 2016-09-05
%%
function [] = muscleContours(subject,type)
directory = '../../../../data';
location = 'Pathologic';
CTn = sprintf('P%03d',subject);
pathologic_list = dir(sprintf('%s/%s/P*', directory, location)); %Create the case list
for i = 1:1:numel(pathologic_list) % Loop over the list of patients and extract the complete name of the specific patient
name = pathologic_list(i).name;
t = strfind(name,CTn);
if t == 1;
patient_name = name;
end
end
finaldirectory = sprintf('%s/%s/%s/CT-%s-1/amira/*tif',directory, location, patient_name, patient_name);
muscleName = cellstr(['SS';'IS';'SC';'TM']);
close all; % Close all figure windows except those created by imtool.
imtool close all; % Close all figure windows created by imtool.
workspace; % Make sure the workspace panel is showing.
fontSize = 16;
imageWindow = [-100 200];
% Choose file
[inputFile, inputPath] = uigetfile(finaldirectory,'Select the image for muscle atrophy measurement');
whereisp = strfind(inputFile, 'P');
if isempty(whereisp)
whereisp = strfind(inputFile, 'N');
if isempty(whereisp)
return
end
end
subject = inputFile(whereisp:whereisp+3);
% Read gray scale image.
inputImage = imread(sprintf('%s%s',inputPath,inputFile));
if size(inputImage,3) == 3
% Calculate the monochrome luminance by combining the RGB values according to the NTSC standard
inputImage = 0.2989*inputImage(:,:,1)...
+0.5870*inputImage(:,:,2)...
+0.1140*inputImage(:,:,3);
end
for i=1:4;
finished = 'Redo';
while strcmp('Redo',finished)
msg = 0; % 0 = no message
% clc; % Clear command window.
close all; % Close all figure windows except those created by imtool.
imtool close all; % Close all figure windows created by imtool.
imshow(inputImage,imageWindow);
imagetitle = sprintf('Select %s ideal section', char(muscleName(i)));
title(imagetitle, 'FontSize', fontSize);
set(gcf, 'Position', get(0,'Screensize')); % Maximize figure.
if msg == 1
message = sprintf('Select the ideal section of your muscle in the image. \nLeft click to place points and draw a polygon. Click again on the first point to close the shape.\nYou can move the points once the line is closed. Double-click or right-click and "Create Mask" to finish.');
msg = 0;
uiwait(msgbox(message));
end
binaryImage1 = roipoly();
if isempty(binaryImage1)
error('Invalid selection or interface closed')
end
binaryImage1 = bwmorph(binaryImage1,'majority');
% Get coordinates of the boundary of the polygonal drawn region.
structBoundaries1 = bwboundaries(binaryImage1);
xy1 = structBoundaries1{1}; % Get n by 2 array of x,y coordinates.
x1 = xy1(:, 2); % Columns.
y1 = xy1(:, 1); % Rows.
hold on
plot(x1, y1, 'color','green');
% check contours
finished = questdlg('Please verify the contours. Press "Save" if they are satisfactory or "Redo" if you would like to redo the selection. ',...
'Contour verification','Save','Redo','Save');
end
%% Report results
% mkdir('./results',subject);
outputPath = sprintf('results/%s/%s/', char(muscleName(i)), subject);
mkdir(outputPath);
outputFile = sprintf('ssContours%s',subject);
% Check for overwriting
overwrite = exist(sprintf('%s%s.mat',outputPath,outputFile), 'file');
switch overwrite
case 2
isoverwrite = questdlg('Overwrite existing file ?', 'Warning', 'No', 'Yes', 'No');
switch isoverwrite
case 'Yes'
save(sprintf('%s%s.mat',outputPath,outputFile), 'inputImage','binaryImage1')
saveas(gcf,sprintf('%s%s.tif',outputPath,outputFile))
fprintf('%s%s.mat\n%s%s.tif\n',outputPath,outputFile,outputPath,outputFile);
end
case 0
save(sprintf('%s%s.mat',outputPath,outputFile), 'inputImage','binaryImage1')
saveas(gcf,sprintf('%s%s.tif',outputPath,outputFile))
fprintf('%s%s.mat\n%s%s.tif\n',outputPath,outputFile,outputPath,outputFile);
end
end
close all;
end
c4science · Help