MATLAB: Problem with sprintf and read DICOM images

ctdicomfilenameimagesimreadslicessprintf

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.

Best Answer

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

Related Solutions

MATLAB: Location pixel and slice

If you had only a 2D question, you could use find, but for some odd reason, it doesn't support more than 3D, so I've made a FEX submission for a similar goal as yours: findND.

[r,c,slice]=findND(X>10000 & X<32000);
val=X(X(:)>10000 & X(:)<32000);

Alternatively:

X_temp=X;
X_temp(X>10000 & X<32000)=0;
[r,c,slice,val]=findND(X_temp);

MATLAB: Error using dicomread, file not found

folder_path = uigetdir; 
dirOutput = dir(folder_path);
dirOutput([dirOutput.isdir]) = [];
fileNames = fullfile(folder_path, {dirOutput.name});

If a file does not have a .dcm extension, then dicomread() will not be able to read the file if you add the .dcm extension. You should leave out the code

ext = '.dcm';
for p = 1:file_num
    fileNames(p) = strcat(fileNames(p),ext);
end

Note: you have

    fileNames(:,[1, 2]) = []; % first 2 filenames don't have have an actual file so i clear the first 2 columns

which is not correct code because some files can sort before '.' and '..' . In the code above, the equivalent but more robust code is the one involving isdir

Best Answer

Related Solutions

MATLAB: Location pixel and slice

MATLAB: Error using dicomread, file not found

Related Question