I have a dicom file of size 2562563. However, when I read this file using dicomread command in matlab the size of the dicom file is 256*256? How to fix this ?

Question

kitty varghese 2018-3-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/387722-i-have-a-dicom-file-of-size-256-256-3-however-when-i-read-this-file-using-dicomread-command-in-mat

评论： kitty varghese 2018-3-13

 close all;
clear all;
clear all;
bb=8; % block size
RR=4; % redundancy factor
K=RR*bb^2; % number of atoms in the dictionary
%here we are reading the image and adding noise into that image.
sigma = 25;% standard deviation
P ='C:\Users\kitty\Dropbox\denoise_ksvd\ADNI';
D=dir(fullfile(P,'*.dcm'));
C=cell(size(D));
for k=1:numel(D)
    C=dicomread(fullfile(P,D(k).name));
   IMin0(:,k)=C(:);

2 个评论
显示无隐藏无

Stephen23 2018-3-12

编辑：Stephen23 2018-3-12

在 MATLAB Online 中打开

@kitty Varghese: please show the outputs of these commands:

size(C)
info = dicominfo(fullfile(P,D(k).name));
info.ColorType
info.Width
info.Height

Note that this line is totally superfluous, as you redefine C on each loop iteration:

C = cell(size(D)); % this line does nothing!
for k = 1:numel(D)
    C = dicomread(fullfile(P,D(k).name)); % because of this line.
    ...
end

kitty varghese 2018-3-12

在 MATLAB Online 中打开

I updated the code and now my code is as shown below.

   P = 'C:\Users\kitty\Dropbox\denoise_ksvd';
S = dir(fullfile(P,'ADNI','*.dcm'));
I = cell(1,numel(S));
for k = 1:numel(S)
   I{k} = dicomread(fullfile(P,'ADNI',S(k).name));
end
M = cell2mat(cellfun(@(m)m(:),I(:).','uni',0));

then size(I)= 1 21

info = dicominfo(fullfile(P,'ADNI',S(k).name))

info.ColorType= grayscale

info.Width=256

info.Height=256

请先登录，再进行评论。

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

Answer 1

Stephen23 2018-3-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/387722-i-have-a-dicom-file-of-size-256-256-3-however-when-i-read-this-file-using-dicomread-command-in-mat#answer_309665

编辑：Stephen23 2018-3-12

在 MATLAB Online 中打开

According to the dicomread documentation "For single-frame grayscale images, X is an M-by-N array". In your response to my comment you wrote

info.ColorType= grayscale

so therefore I would expect the imported image to have size MxN, just as the dicomread help states. If all of the other images are also grayscale, then all of them will be 2D matrices. You can check this yourself in the loop.

In the original question title you wrote "I have a dicom file of size 256*256*3": how are you checking this? With what tool or command?