Dicom files to one 2d matrix

Hello ,

I have dicom files with 210 slices. I would like to use dicom read to read all slices then I would like to have a matrix with slice spacing in the y axe and all HU numbers in the width for each slice.

any help?

Thank you

Answer 1

[Selva Karna]

clc

clear all;

close all;

warning off

myFolder='your dicom root';

filePattern = fullfile(myFolder,'*.dcm');

dcmfiles = dir(filePattern);

flen=length(dcmfiles);

for k = 1:flen

baseFileName =dcmfiles(k).name;

fullFileName = fullfile(myFolder, baseFileName);

a = dicomread(fullFileName);

dcm_vlm(:,:,k)=a;

end

%% To view 3d volume

volshow(squeeze(dcm_vlm));