how to create 3-D volume of selected DICOM images ?

8 次查看(过去 30 天)
Hello , I have written a MATLAB code for the operation of converting selected DICOM images into its corresponded 3D volume . But I am not getting any response for longer time as it should complete its execution in limited time . I checked 'mri.mat' and its variable 'D' for reference . Please tell me where i should do correction in my code ? My code is -
clear all
close all
[filename, pathname] = uigetfile( ...
{'*.dcm','DCM Files (*.dcm,)';
'*.dcm', 'DICOM Images (*.dcm)'; ...
'*.*', 'All Files (*.*)'}, ...
'MultiSelect', 'on',...
'Pick file/(s)');
file_len = length(filename);
for i=1:file_len
I{i} = dicomread(filename{i});
Q{i}=imresize(I{i},0.4);
end
II = cell2mat(Q);
AA=uint8(II);
save('new1.mat','II','AA');
load new1.mat;
K = squeeze(AA);
L = padarray(K,[10 10 10],'both');
% Create an isosurface
Ds = smooth3(L);
surface = isosurface(Ds,5);
% Display the surface
hiso = patch('Vertices',surface.vertices,...
'Faces',surface.faces,...
'FaceColor',[.2,.8,.9],...
'EdgeColor','none');
grid on;
view(27,14)
axis tight
daspect([1,1,.35])
lightangle(40,30);
set(gcf,'Renderer','zbuffer'); lighting phong
isonormals(Ds,hiso)
set(hiso,'SpecularColorReflectance',0,'SpecularExponent',50)
Thank you :)

采纳的回答

Walter Roberson
Walter Roberson 2015-9-12
Your code
for i=1:file_len
I{i} = dicomread(filename{i});
Q{i}=imresize(I{i},0.4);
end
II = cell2mat(Q);
is creating II (and therefor the later AA) by concatenating the data long the second dimension. I suspect you wanted to concenate along the third dimension. If so then right before the above for loop insert
Q = cell(1,1,file_len);
Also I do not understand why you later squeeze(AA). Which dimension might be singleton?
When I am reading other people's code, I am more comfortable with them using permute() or reshape() rather than squeeze(), at least if it is not very obvious that one dimension has just been made singleton and the other dimensions are definitely not singleton. When the reader of the code does not know ahead of time what the dimensions of the data are to be, the reader of the code must assume that squeeze() could be operating upon an unexpected dimension, such as if it turned out that the data only had a single row. permute() and reshape() remove the possibility of changing any dimension not explicitly intended to be changed.
  8 个评论
Ringo Cheung
Ringo Cheung 2018-5-7
In the code, it says that :
Q{i}=imresize(I{i},0.4);
But I am not quite understand why it needs to resize the model by a factor of 0.4, can you briefly explain this?
Thanks, Ringo
Walter Roberson
Walter Roberson 2018-5-7
That imresize() was there in the original, so I presume that the original poster had some good reason for it.

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 DICOM Format 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by