请教如何画四维图(4D-plot)

Hello_Mr.li

        想用matlab建一个地下地质体模型，大致是一个长方体。这个长方体由很多个单位立方体组成，并且每个单元都有一个数值，用颜色显示数值的大小，数值都是非负数，0值显示无色，即这个单元不显示。
        也就是在matlab里可以用什么方法画这样的四维图，变量为[x,y,z,v],[x,y,z]确定小单元的位置，v表示大小并用颜色显示，0值是无色（这个若不好处理显示一种背景色也可以）。
         我用了slice试一下，但没达到想要的效果（理想的模型，就像魔方从上表面第一层第二层拿掉几个小方块）。   求指点！

• %用切片函数构出长方体
• [x,y,z]=meshgrid(-1:.1:1,-1:.1:1,-3:.1:3);
• %画图范围
• v=x.*exp(-x.^2-y.^2-z.^2);
• %计算函数值
• slice(x,y,z,v,[],[],-2:.01:2,'linear');
• axis([-3 3 -3 3 -3 3]);
• shading flat
  

ExxNEN

下面不是有代码了吗

SsaaM7

你这个要求有点厉害了啊

haidaowang

因为slice画图的时候，默认的 EdgeColor 是黑色，所以，尽管你将颜色数据设置为了NAN, 仍会显示黑色的edge。下面的设置去掉了edge颜色。另外，为了便于观察，调整了一下视角：
- j5 Y3 C/ P& }1 s) |1 u: Z

• [x,y,z] = meshgrid(0:.5:10,0:.5:10,0:.5:10);
• c = x.^2+y.^2+z.^2;
• xs = 0:0.5:10;
• ys = xs;
• zs = xs;
• c(7:15,7:15,13:21)=NaN;
• h = slice(x,y,z,c,xs,ys,zs);
• set(h,'FaceColor','interp',...
•     'EdgeColor','none')
• camproj perspective
• box on
• view(-70,70)
• colormap hsv
• colorbar
  ! B7 b* ^% q/ \

效果如下：
& Q/ {+ U! ^# L* q5 Q
# `! G8 r4 M2 r1 a/ w; a
5 Z/ P' [2 b! i5 P, ]; ]9 e