2

我期待将 3D 矩阵从 dicom 文件可视化到 matlab 中。由于我对matlab不太熟悉，因此我设法从这篇文章中获得了帮助：

不同之处在于我的矩阵不是由 1 和 0 组成，而是由正确红色的负数组成imshow(dicomeread(dicomFile))

如何获得相同的对比度但使用 3D 渲染？

我的代码：

dicomFilesZm = dir(fullfile(myDcmFolder, 'SLZ-*.dcm')); %Get files name
dicomFilesZp = dir(fullfile(myDcmFolder, 'SLZ+*.dcm')); %~

Z = dicomFilesZm(end:-1:1); % sort

dicomFilesZ = [Z ; dicomFilesZp]; % recompose final array with files name

Iz1 = fullfile(myDcmFolder, dicomFilesZ(1).name); 
v = NaN([size(dicomread(Iz1)) numel(dicomFilesZ)]); % creation of empty matrix with the good size    
for i = 1 : numel(dicomFilesZ)
    Iz = fullfile(myDcmFolder, dicomFilesZ(i).name);
    v(:,:,i) = dicomread(Iz); % fill the matrix with each image
end

p = patch( isosurface(v,0) );
isonormals(v, p)
set(p, 'FaceColor','r', 'EdgeColor','none')
daspect([1 1 1])

谢谢您的帮助。

4

1 回答 1

3

解决您的问题的两种选择：

  1. 天真地，您为什么不直接设置v=v+min(v(:)) 图像将从零缩放到不同的最大值？

  2. 为什么不isosurface(V,isovalue)使用负 isovlaue 为您解决这个问题？或者，如果您想遍历 (-n:step_size:m) 中的等值以获得所需的动态范围（其中一些alpha(0.2)可以查看所需的图层）

于 2012-10-27T00:06:32.283 回答