DigitalBreastTomosynthesis/Functions/readDicom.m

%% Author: Rodrigo de Barros Vimieiro
% Date: May, 2018
% rodrigo.vimieiro@gmail.com
% =========================================================================
%{
% -------------------------------------------------------------------------
%                         readDicom(imgdir,parameter,answer)
% -------------------------------------------------------------------------
%     DESCRIPTION:
%     This function read a set of Dicom images from a directory.
%   
%     INPUT:
%     - imgdir = Directory for Dicom files 
%     - parameter = Parameter of all geometry 
% 
%     OUTPUT:
%     - dataDicom = Stack of Dicom images.
%     - infoDicom = Stack of Dicom headers.
% 
%     ---------------------------------------------------------------------
%     Copyright (C) <2018>  <Rodrigo de Barros Vimieiro>
% 
%     This program is free software: you can redistribute it and/or modify
%     it under the terms of the GNU General Public License as published by
%     the Free Software Foundation, either version 3 of the License, or
%     (at your option) any later version.
% 
%     This program is distributed in the hope that it will be useful,
%     but WITHOUT ANY WARRANTY; without even the implied warranty of
%     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%     GNU General Public License for more details.
% 
%     You should have received a copy of the GNU General Public License
%     along with this program.  If not, see <http://www.gnu.org/licenses/>.
%}
% =========================================================================
%% Read Dicom images from files
function [dataDicom,infoDicom] = readDicom(imgdir,parameter,answer)

if(strcmp(parameter.type,'vct'))
    dc = 1;
    typecase = 1;
end
if(strcmp(parameter.type,'hologic'))
    dc = 1;
    typecase = 0;
end
if(strcmp(parameter.type,'ge'))
    dc = 0;
    typecase = 1;
end
if(strcmp(answer,'Microcalcifications') || strcmp(answer,'Noised'))
    dc = 0;
    typecase = 0;
end

img_list = dir([imgdir,filesep,'*.dcm']);  % List dicom files
img_count = size(img_list,1);

fprintf('%2d/%2d', 1, img_count)

for i=1:img_count
    fprintf('\b\b\b\b\b%2d/%2d', i, img_count)
    imgAux = single(dicomread([imgdir, filesep, img_list(i,:).name]));  % Read Dicom
    infoAux = dicominfo([imgdir, filesep, img_list(i,:).name]);         % Read Dicom headers
    
    if(typecase == 1)
        nProj = infoAux.InstanceNumber+dc;
    else
        filename = strtrim(img_list(i,:).name);
        filenameSlpited = strsplit(filename, '_');
        filenameSlpited = strsplit(filenameSlpited{end}, '.');
        nProj = str2double(filenameSlpited{1})+dc;
    end 

    dataDicom(:,:,nProj) = imgAux;
    infoDicom(:,nProj) = infoAux;
end

end