这个例子展示了在MATLAB®中可视化体积数据的几种方法。
一个等值面是一个曲面,其中空间内的所有点的值都是常数。使用等值面
函数生成表面和外部的面和顶点isocaps
函数生成体的端盖的面和顶点。使用补丁
命令绘制音量和它的结束盖。
负载核磁共振成像D%加载数据D =紧缩(D);删除单例维度limit = [NaN NaN NaN NaN 10];[x, y, z, D] = subvolume(D, limits);%提取卷数据的子集(fo,签证官)=等值面(x, y, z D 5);%等值面为外部体积(铁、ve、ce) = isocaps (x, y, z D 5);% iscaps表示卷的结束caps图p1 = patch(“脸”佛,“顶点”,签证官);%画体积的外面p1。FaceColor =“红色”;p1。EdgeColor =“没有”;p2 =补丁(“脸”、铁、“顶点”已经,...%绘制音量的末端caps“FaceVertexCData”、ce);p2。FaceColor =的插值函数;p2。EdgeColor =“没有”;视图(-40,24)daspect([1 1 0.3])设置轴的长宽比colormap(灰色(100))的盒子在camlight(40、40)%创建两个灯camlight(-20、-10)照明高洛德
的coneplot
命令将速度矢量绘制成锥在x,y,z体积中的点。圆锥表示每一点上矢量场的大小和方向。
班清除当前坐标轴负载风uvwxyz%加载数据(m, n, p) = (u)大小;[Cx, Cy, Cz] = meshgrid(1:4:m,1:4:n,1:4:p);%计算圆锥的位置h = coneplot (u, v, w,残雪,Cy, Cz, y, 4);画圆锥图集(h,“EdgeColor”,“没有”)轴紧平等的37岁的视图(32)在colormap (hsv)光
的简化
函数绘制速度向量的流线x,y,z来说明三维向量场的流动。
class [m,n,p] = size(u);[Sx, Sy, Sz] = meshgrid(1,1:5:n,1:5:p);计算流线的起始点简化(u, v, w, Sx, Sy,深圳)画流线轴紧平等的37岁的视图(32)在
的廉管道
函数绘制速度矢量在的流管x,y,z体积中的点。管的宽度与向量场在每一点的归一化散度成正比。
class [m,n,p] = size(u);[Sx, Sy, Sz] = meshgrid(1,1:5:n,1:5:p);计算流线的起始点h =廉管道(u, v, w, Sx, Sy, Sz);%绘制流管并返回一个曲面数组集(h,“FaceColor”,“青色”)%使用“set”更改对象数组的属性集(h,“EdgeColor”,“没有”)轴紧平等的37岁的视图(32)在光
将体积可视化结合在一个单一的图中,以获得一个更全面的速度场在一个体积。
class SPD =√(u。*u + v *v + w.*w);体积中每一点的风速%(fo,签证官)=等值面(x, y, z,社民党,40);%等值面为外部体积(铁、ve、ce) = isocaps (x, y, z,社民党,40);% iscaps表示卷的结束capsp1 =补丁(“脸”佛,“顶点”,签证官);%画出体积的等值面p1。FaceColor =“红色”;p1。EdgeColor =“没有”;p2 =补丁(“脸”、铁、“顶点”已经,...%绘制音量的末端caps“FaceVertexCData”、ce);p2。FaceColor =的插值函数;p2。EdgeColor =“没有”;[fc, vc] =等值面(x, y, z, spd, 30);圆锥的等值面[fc, vc] = reducepatch(fc, vc, 0.2);减少面和顶点的数量h1 = coneplot (x, y, z, u, v, w, vc (: 1), vc (:, 2), vc (:, 3), 3);画等高线h1。FaceColor =“青色”;h1。EdgeColor =“没有”;[sx, sy, sz] = meshgrid(80, 20:10:50, 0:5:15);%起始点流线h2 =简化(x, y, z, u, v, w, sx, sy, sz);画流线集(h2,“颜色”,(。4 1 4])轴紧平等的37岁的视图(32)在光