StackingDicoms.m

function [a] = StackingDicom(filename,short,nFrames)
%Read in file information
info = dicominfo(filename);
nRows = info.Rows;
nColumns = info.Columns;
%nPlanes = info.SamplesPerPixel;
X = repmat(int16(0), [nRows, nColumns, nFrames]); %create empty array
%Stacking slices
for p=0:nFrames-1
    filename = sprintf('%s%05d.DCM', short,p);
    X(:,:,p+1) = dicomread(filename);%add files to array
end
a = size(X);


end