Question

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

不同之处在于我的矩阵不是由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

解决问题的两个选项：

天真地，你为什么不设置v=v+min(v(:))所以图像会从零扩展到不同的最大值？
为什么带有负isovlaue的isosurface(V,isovalue)不会为您解决此问题？或者，如果你想要从（-n：step_size：m）的isovalues循环获得你想要的动态范围（用一些alpha(0.2)来查看你想要的图层）

Matlab 3D体积可视化 - dicom文件

1 个答案: