MATLAB Answers

0

Dicomread images with different x and y sizes into a cellarray.

Asked by GioPapas81 on 18 May 2018
Latest activity Edited by GioPapas81 on 18 May 2018

Hi. This may be easy for some Matlab users. I am trying to dicomread and process images starting my code with the following commands:

[Filename1,Pathname1]=uigetfile('Multiselect','on');
for i=1:length(Filename1)
    images(:,:,i)=dicomread(fullfile(Pathname1,Filename1{i}));
    images=double(images);
end

However, my images appear to change size therefore I get the error: 'Subscripted assignment dimension mismatch.'

Any ideas how I could dicomread and place in order all dicom images independently of their size, would be much appreciated.

Thank you

  0 Comments

Sign in to comment.

1 Answer

Answer by Guillaume
on 18 May 2018
 Accepted Answer

You've answered your question in its title. Use a cell array:

[Filename1,Pathname1]=uigetfile('Multiselect','on');
images = cell(size(Filename1));
for fileidx = 1:numel(Filename1)
  images{fileidx} = double(dicomread(fullfile(Pathname1,Filename1{i})));
end

  4 Comments

Show 1 older comment

I've not changed the way you read the files. It's a bit puzzling why the file would not be found since it's been selected with uigetfile but that's a completely different issue to your question.

Using a cell array, it will read all images. But the output is obviously not a 3D matrix. it is not clear to me whether your requirement to store the images as images(:, :, i), that is a 3D matrix, is one of lack of familiarity with matlab or truly a necessity. Why do you need the images as a 3D matrix? Isn't the cell array good for you?

If you do need a 3D matrix, then yes, you'll either have to drop the images that don't match your required size, or rescale them.

I'm not even sure that you NEED to have all the images in memory simultaneously. Why do you think you do? Usually it's not necessary. Just process each image one at a time in the loop and don't store them all.

I re-sorted my files, tried again and it worked (I am not sure why I got this error in first place).

I think that I have mis-used the term 'cellarray', instead of using 3D matrix initially but your are right, I can replace the 3D matrix with a cellarray. I only need to incorporate changes in a relatively long code, but there is no free lunch in life.

Thank you again both.

Sign in to comment.