7-9.m

来自「MATLAB数字图象处理的主要源码.包括图象去噪变换增强恢复识别压缩的各种算法」· M 代码 · 共 19 行

load wind                                       
%打开保存的数据
lims = [100.64 116.67 17.25 28.75 -0.02 6.86];         
%定义坐标范围
[x,y,z,u,v,w] = subvolume(x,y,z,u,v,w,lims);            
%lims来定义数据子集
[sx sy sz] = meshgrid(110,20:5:30,1:5);               
%定义网格点
verts = stream3(x,y,z,u,v,w,sx,sy,sz,.5);               
%计算彩带顶点
cav = curl(x,y,z,u,v,w);                            
%计算卷曲角速度
wind_speed = sqrt(u.^2 + v.^2 + w.^2);               
%计算流速
h = streamribbon(verts,x,y,z,cav,wind_speed,2);        
%绘制流彩带图
view(3)                                        
%三视图