Matlab 3D 剂量阵列可视化

Question

我在 Matlab 中有一个 3D 矩阵，它对应于能量剂量传递矩阵。我正在尝试可视化和表示矩阵。目前我正在执行以下操作 ( thanks to another post )。

diff = double(squeeze(diff));
h = slice(diff, [], [], 1:size(diff,3));
set(h, 'EdgeColor','none', 'FaceColor','interp')
alpha(.1)

这给了我一个很好的 3D 图，但是仍然很难正确地看到它，我必须不断旋转图形才能正确地可视化它。我也用过:

isosurface(diff,'isovalue')

但还是很难看到任何东西。

我想知道是否有办法去除实际剂量表示周围的蓝色区域，因为蓝色区域对应于 0 值。也许摆脱它可以帮助我看到更清晰的画面。

Answer 1

您可以将零设置为nan，这样它们就不会被绘制

diff = double(squeeze(diff));
diff(diff==0)=nan;                              % added line
h = slice(diff, [], [], 1:size(diff,3));
set(h, 'EdgeColor','none', 'FaceColor','interp')
alpha(.1)

例子

使用来自上一个问题的 MRI 数据

load mri
D = double(squeeze(D));
D(D==0)=nan;
h = slice(D, [], [], 1:size(D,3));
set(h, 'EdgeColor','none', 'FaceColor','interp')
alpha(.1)

输出

用 nan 替换零
没有用 nan 替换零