www.gusucode.com > IPCV_Eval_Kit_R2019a_0ce6858工具箱matlab程序源码 > IPCV_Eval_Kit_R2019a_0ce6858/code/demo_files/I1_02_grademo/mixplot.m

    f1=figure(1); clf reset
set(f1,'units','normalized','position',[0.3652 0.3008 0.6016 0.6016])

load wind
vel = sqrt(u.*u + v.*v + w.*w);

disp('iso...')
p = patch(isosurface(x, y, z, vel, 40));
isonormals (x, y, z, vel,p)
set(p, 'facecolor', 'r', 'edgecolor', 'n');
p2 = patch(isocaps(x,y,z,vel,40), 'facec', 'i', 'edgec',  'n');
caxis([40 65])


daspect([1 1 1])

[ff vv] = isosurface(x, y, z, vel, 30);
[f verts] = reducepatch(ff,vv, .2);

disp('cones...')
h=coneplot(x,y,z,u,v,w,verts(:,1),verts(:,2),verts(:,3),2);
set(h, 'facec',  [0 .7 .7], 'edgec', 'n')


disp('streams...')
[sx sy sz] = meshgrid(71, 18:8:59, 0:5:15);
verts = stream3(x,y,z,u,v,w,sx,sy,sz);
  sl=streamline(verts);
  %set(sl, 'color', [.3 .3 .3]);
axis vis3d tight

view(3)
camlight 
lighting p
%set(gcf, 'color', [.3 .3 .3])
%set(gca, 'color', 'k')
%set(gca, 'xcolor', 'w')
%set(gca, 'ycolor', 'w')
%set(gca, 'zcolor', 'w')


box on;
camproj p;
camva(33)
campos([150 -25 60]); camtarget([105 40 0])