how to form the volume from the single dicom image using isosurface in matlab

1 次查看（过去 30 天）

andhavarapu lokesh 2016-11-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/314665-how-to-form-the-volume-from-the-single-dicom-image-using-isosurface-in-matlab

评论： Walter Roberson 2016-12-1

want to do volume rendering for the data containing dicom images . Each dicom contains again 72 frames can any one help me how to get the volume from the below code

rojectdir = 'E:\SHIVA BACKUP\THYROID\P1\newcodes\data1\13002';
% y = length(projectdir);
y=72;
X = zeros(128, 128, 1, 72, y, 'uint8');
% Read the series of images.
for p=1:1:y
     thisfile = sprintf('IM_%d.dcm', p);
     filename = fullfile( projectdir, thisfile );
     imdata = dicomread(filename);
     imsize = size(imdata);
     if ~isequal( imsize, [128 128 1 72] )
       fprintf('file is unexpected size %s instead of [128 128 1 72], skipping "%s"\n', mat2str(imsize), filename);
     else
       X(:, :, :, :, p) = imdata;
     end
  end
isoval=0.5;
hiso=patch(isosurface(X,isoval),...
       'FaceColor',[1,0.75,0.65],'EdgeColor','none');%this is for volume display 
    set(hiso,'FaceAlpha',0.74);
    lighting phong;
    lightangle(45,30);
    rotate3d on;

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

回答（2 个）

Walter Roberson 2016-11-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/314665-how-to-form-the-volume-from-the-single-dicom-image-using-isosurface-in-matlab#answer_245299

编辑：Walter Roberson 2016-11-30

在 MATLAB Online 中打开

volumes_of_objects = squeeze( sum(sum(sum(X,1),2),3) );

Note: this assumes, though, that voxels are cubic. If your distance between z slices is different than the size of the pixels in the x and y directions, then you should multiply the above by the ratio of the z slice distance compared to the x distance. Or multiply by the physical X pixel width, physical Y pixel width, and physical Z distance: you can extract those parameters from the DICOM information structure.

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

andhavarapu lokesh 2016-12-1

编辑：Walter Roberson 2016-12-1

在 MATLAB Online 中打开

sir still iam getting error iam telling frankly i don,t know any coding please help me sir

projectdir = 'E:\SHIVA BACKUP\THYROID\P1\newcodes\data1\13002';
% y = length(projectdir);
y=72;
X = zeros(128, 128, 1, 72, y, 'uint8');
% Read the series of images.
for p=1:1:y
     thisfile = sprintf('IM_%d.dcm', p);
     filename = fullfile( projectdir, thisfile );
     imdata = dicomread(filename);
     imsize = size(imdata);
     if ~isequal( imsize, [128 128 1 72] )
       fprintf('file is unexpected size %s instead of [128 128 1 72], skipping "%s"\n', mat2str(imsize), filename);
     else
       X(:, :, :, :, p) = imdata;
     end
  end
for p = 1 : y
  fig = figure();
  ax = axes('Parent', fig);
  v = squeeze(X(:,:,:,:,p));
  hiso = patch( isosurface( v, isoval), ...
     'Parent', ax, 'FaceAlpha', 0.74, ...
     'FaceColor', [1,0.75,0.65],' EdgeColor', 'none');
  lighting(ax, 'phong');
  lightangle(ax, 45, 30);
  rotate3d(ax, 'on');
  title( sprintf('p = %d', p) )
end

and also the error is

Undefined function or variable 'isoval'.
Error in Untitled3 (line 24)
  hiso = patch( isosurface( v, isoval), ...
 There is no  EdgeColor property on the Patch class.
Error in Untitled3 (line 25)
  hiso = patch( isosurface( v, isoval), ..