Can we read dicom files without extension using the dicomread function?
显示 更早的评论
I have 71 image slices and the idea is to make an image volume using these slices. I tried using the following:
close all; clear; clc
%%
Folder = '/etc/';
filenames = fullfile(Folder,'Z*'); % files names are Z01, Z02, ....., Z75
findDir = dir(filenames);
total_Im = numel(filenames);
%%
for n = 1:total_Im
f = fullfile(Folder,filenames(n).name);
dm_images = dicomread(f);
figure(n)
imshow(dm_images)
end
% How can I read these files without extension and make a single image
% volume?
采纳的回答
Rik
2022-6-28
Store the read images in a 3D array. That way you will have the entire volume in one variable, allowing you to do further processing as you need to.
clear dm_images
for n = 1:total_Im
f = fullfile(Folder,filenames(n).name);
dm_images(:, :, n) = dicomread(f);
end
4 个评论
blues
2022-6-28
Rik
2022-6-28
I didn't read your code carefully enough. I had assumed (based on how you treat it) that filenames would be output of the dir call, instead of the input to it. Simply exchange it with findDir to prevent the error.
You can also remove that first line of code.
blues
2022-6-29
Rik
2022-6-29
I would suggest using the debugger to go through your code line by line. Then you would discover your first loop doesn't do anything. You would also find out that your file name variable doesn't actually contain a file name, but is instead a struct. That makes dicomread assume it is a dicominfo struct.
Folder = '/etc/';
findDir = fullfile(Folder,'*');
filenames = dir(findDir);
total_Im = length(filenames);
clear dm_images
for n = 1:total_Im
f = fullfile(Folder,filenames(n).name);
dm_images(:, :, n) = dicomread(f);
end
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Read and Write Image Data from Files 的更多信息
Translated by 
