How to combine all CT files to make a single array/matrix?

%load and read CT images
% Preallocate the 512-by-512-by-1-by-263 image array
CT = zeros(512, 512, 263);
 for p = 1:263
   CTfilename = sprintf('CT-0003-%03d.IMA', p);
   CT = dicomread(CTfilename);     
   %imshow(X,[])
 end
% Now I want to make a 3D array of all series
3DArray = cat(3, CT{:}); % how to write a code here so that I can get a 3D array with X*Y*Z dims?

[Geoff Hayes]

blues - your comment says

% Preallocate the 512-by-512-by-1-by-263 image array

but your code to pre-allocate the array removes the 1 and creates a 512x512x263 matrix. Which is correct? Or should your code to save the images be

CT = zeros(512, 512, 263);
for p = 1:263
    CTfilename = sprintf('CT-0003-%03d.IMA', p);
    CT(:,:,p) = dicomread(CTfilename);     
    %imshow(X,[])
end

assuming that dicomread returns a 512x512 matrix?

[blues]

Thank you for pointing that out.

CT without CT(:,:,p) returned 512*512.

Now, with CT(:,:,p), I will not need concatenate (cat) thing, right?

Answer 1

[Sean de Wolski]

Look at dicomreadVolume!

https://www.mathworks.com/help/images/ref/dicomreadvolume.html