How to combine all CT files to make a single array/matrix?
1 次查看(过去 30 天)
显示 更早的评论
%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?
2 个评论
Geoff Hayes
2020-1-31
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?
回答(1 个)
Sean de Wolski
2020-1-31
Look at dicomreadVolume!
https://www.mathworks.com/help/images/ref/dicomreadvolume.html
0 个评论
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!