How to convert 3D to 4D images

Hi all

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

anyone can help me?

Answer 1

[Simon Chan]

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

Answer 2

[DGM]

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]);

[mohd akmal masud]

TQ Sir DGM and Simon Chan, Its work!!