MATLAB 3D卷可视化 - DICOM文件

Question

我期待将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])

谢谢您的帮助。

Answer 1

解决您的问题的两个选项：

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

2. 为什么isosurface(V,isovalue)不为您解决这个问题？或者，如果您想从（-n：step_size：m）绕过Isovalues，以获取所需的动态范围（使用一些alpha(0.2)查看所需的层）