Problem with sprintf and read DICOM images

5 次查看（过去 30 天）

显示更早的评论

Veronika 2017-4-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/333653-problem-with-sprintf-and-read-dicom-images

评论： Walter Roberson 2017-11-17

采纳的回答： Jan

在 MATLAB Online 中打开

Dear all,

I have CT slices and I would like to load (imread) them to the Workspace. But I always read only one (first) slice (image). This is my code:

for p=1:38
         filename = sprintf('C:/Users/ced0005/DP/DP_framework/DICOM/slicesCT/000048.dcm',p);
          X(:,:,1,p) = dicomread(filename);

I don´t know, where mistake is.

Can you advise me? Thank you for your answers.

采纳的回答

Jan 2017-4-4

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/333653-problem-with-sprintf-and-read-dicom-images#answer_261768

编辑：Jan 2017-4-4

在 MATLAB Online 中打开

There is no format specifier in the format string of sprintf:

sprintf('C:/Users/ced0005/DP/DP_framework/DICOM/slicesCT/000048.dcm',p);

This is the same string for all iterations. Perhaps you want:

filename = sprintf('C:/Users/ced0005/DP/DP_framework/DICOM/slicesCT/%06d.dcm', p);

12 个评论
显示 10更早的评论隐藏 10更早的评论

Walter Roberson 2017-11-9

在 MATLAB Online 中打开

first_to_read = 48;
last_to_read = 86;
num_to_read = last_to_read - first_to_read + 1;

Now you can use

for file_idx = 1 : num_to_read
  img_number = file_idx + first_to_read - 1;
  filename = fullfile(Folder, sprintf('%06d.dcm', img_number));
  X(:, : , 1, file_idx) = dicomread(filename);
end

in which case file_idx refers to the relative number to read, which can be quite useful for several different reasons.

You could also use

for file_number = first_to_read : last_to_read
  file_idx = file_number - first_to_read + 1;
  filename = fullfile(Folder, sprintf('%06d.dcm', img_number));
  X(:, : , 1, file_idx) = dicomread(filename);
end

which is sometimes easier to read. But in practice it turns out to be easier to work with relative indices in a for loop: when relative indices are used it is not necessary that the values be consecutive:

files_to_read = [48:63 72:78 85];   %does not need to be consecutive
num_to_read = length(files_to_read);
for file_idx = 1 : num_to_read
  img_number = files_to_read(file_idx);
  filename = fullfile(Folder, sprintf('%06d.dcm', img_number));
  X(:, : , 1, file_idx) = dicomread(filename);
end