reconstruct 3D image

Question

Zaitul 2016-11-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/310942-reconstruct-3d-image

编辑： Ahmet Cecen 2016-11-10

Hi, i have 20 of 2D images in dicom format and i want to reconstruct the 3D image. How can i do this in MATLAB? here i attached my code and i hope anyone could help me to fix it. Many thanks

a = dicomread ('brain_001.dcm')
mask = zeros(size(a));
mask(100:end-50,100:end-50) = 1;
bw = activecontour(a,mask,1000);
% for p=1:20
%    filename = sprintf('brain_%03d.dcm', p);
%    a(:,:,1,p) =double(dicomread(filename));
% end
 c = im2double(bw);
 shading flat
 d = imgaussfilt3 (c,2);
 for p=1:20
   filename = sprintf('brain_%03d.dcm', p);
   a(:,:,1,p) =double(dicomread(filename));
end
 colormap(copper)
  e = mesh(d.*100)

3 个评论
显示 1更早的评论隐藏 1更早的评论

Walter Roberson 2016-11-7

implay looks for M-by-N-by-1-by-K for grayscale.

Zaitul 2016-11-10

编辑：Zaitul 2016-11-10

i want to make the 3D image from those 2D images. what should i do for c and d? i'm a bit confused when i'm doing the loop for p. i'm sorry as i'm just realized there's only 20 slices of image actually.

请先登录，再进行评论。

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

Answer 1

Ahmet Cecen 2016-11-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/310942-reconstruct-3d-image#answer_242705

编辑：Ahmet Cecen 2016-11-10

在 MATLAB Online 中打开

Too vague to answer directly, but here are two helpful tips:

1) Whatever you are doing, in the end you need to stack the images on the 3rd dimension not 4th.

NOT

a(:,:,1,p) =double(dicomread(filename));

BUT

a(:,:,p) =double(dicomread(filename));

2) Do not try to visualize 3D data in MATLAB, unless it is a surface. MATLAB figures (as far as I know) do not support pixel shading, so your figures will look very flat. Instead use Paraview. You can use this function to export your 3D matrix once you have that in MATLAB. Below is an example of a 3D reconstruction I visualized following the same procedure.Notice the depth/shading.