How to convert 3D to 4D images

Question

mohd akmal masud 2022-2-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1649870-how-to-convert-3d-to-4d-images

评论： mohd akmal masud 2022-2-15

采纳的回答： Simon Chan

ZubalPhantomDicom.zip

Hi all

I have image dicom as 3D as (attached). How to combine it all the slice to 4D?

anyone can help me?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Simon Chan 2022-2-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1649870-how-to-convert-3d-to-4d-images#answer_896185

在 MATLAB Online 中打开

Try thw following:

[filename,pathname]=uigetfile('*.dcm','Select Files', 'MultiSelect', 'on');
Nz = length(filename);
for k = 1: Nz
    imagedata = dicomread(fullfile(pathname,filename{k}));
    if k == 1
        [Ny,Nx] = size(imagedata);
        rawdata = zeros(Ny,Nx,Nz);
    end
    rawdata(:,:,k) = imagedata;                             % Convert to 3D
end
Dimage = permute(rawdata,[1 2 4 3]);                        % Convert to 4D

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

DGM 2022-2-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1649870-how-to-convert-3d-to-4d-images#answer_896180

编辑：DGM 2022-2-15

在 MATLAB Online 中打开

If you have a volumetric image represented in a 3D array and you want to slice it on dim 3 and arrange each slice into a frame in a 4D image, you can do

B = permute(A,[1 2 4 3]);

If you want to slice on another dimension, just rearrange the first,second, and fourth elements of that vector in the call to permute().

EDIT:

A concrete example:

dirname = 'ZubalPhantomDicom';
dicomlist = dir(fullfile(pwd,dirname,'*.dcm'));
imstack = cell(numel(dicomlist),1);
for cnt = 1:numel(dicomlist)
    imstack{cnt} = dicomread(fullfile(pwd,dirname,dicomlist(cnt).name));  
end
% the image is currently a cell array of pages
% say we arrange the pages on dim3
imstack3 = cat(3,imstack{:});
% let's say we want to arrange the pages on dim4 instead
imstack4 = cat(4,imstack{:});
% let's say we wanted to take imstack3 and turn it into imstack4
imstack34 = permute(imstack3,[1 2 4 3]);