example_isosurfaces.m

基础程序

% example_isosurfaces.m
% 不同参数情况下的结果
% 生成绘制图形数据
[x,y,z,v]=flow(13);
% 通过isosurface函数寻找等值面
fv=isosurface(x,y,z,v,-2);
% figure1：原始图形;
subplot(2,2,1);
p=patch(fv);
set(p,'facecolor',[.5 .5 .5],'edgecolor','black');
view(-45,30);axis equal tight;grid on;
title('figure a: original');
% figure2：reduce volume;
subplot(2,2,2);
[xr,yr,zr,vr]=reducevolume(x,y,z,v,[3,2,2]);
fvr=isosurface(xr,yr,zr,vr,-2);
p=patch(fvr);
np=size(get(p,'faces'),1);
set(p,'facecolor',[.5 .5 .5],'edgecolor','black');
view(-45,30);axis equal tight;grid on;
title('figure b: reduce volume');
% figure3：reduce patch;
subplot(2,2,3);
p=patch(fv);
set(p,'facecolor',[.5 .5 .5],'edgecolor','black');
view(-45,30);axis equal tight;grid on;
reducepatch(p,0.15);
np=size(get(p,'faces'),1);
title('figure c: reduce patch');
% figure4：smooth3;
subplot(2,2,4);
p=patch(fv,'facecolor','blue','edgecolor','none');
view(-45,30);axis equal tight;grid on;
title('figure c: reduce patch');