Question

Matlab 中如何得到用isosurface 形成的曲面的点集

rt

Answer 1

//isosurface函数用于由体积数据中提取等值数据。
//使用方法如下
fv =

//使用方法如下
fv = isosurface(X,Y,Z,V,isovalue)
fv = isosurface(V,isovalue)
fvc = isosurface(...,colors)
[f,v] = isosurface(...)
[f,v,c] = isosurface(...)
isosurface(...)
//参数意义为，从由X,Y,Z构造的体积V数据中提取由isovalue指定的等值数据，返回结果fv为一个结构体，包含了等值面的顶点和面（顶点的次序），这些参数可以直接传给patch命令画出图形来。
//另外，可以在最后加上colors，返回值中也会得到每个面的颜色插值，这样能够使你用不同于计算等值面的数据控制等值面的颜色映射。

//示例：

[x y z v] = flow;
[fc vt] = isosurface(x, y, z, v, -3);%vt就是曲面点集
p=patch('Faces',fc,'Vertices',vt);
       isonormals(x,y,z,v, p)
       set(p, 'FaceColor', 'red', 'EdgeColor', 'none');
       daspect([1 1 1])
       view(3)
       camlight; lighting phong
figure
plot3(vt(:,1),vt(:,2),vt(:,3),'.')

Answer 2

这个命令的使用要根据程序需要而定的，它的输出参数是一个结构，可直接传给patch函数进行画图，如
[x, y, z, v] = flow;
q = z./x.*y.^3;
k=isosurface(x, y, z, q, -.08, v);
p = patch(k);
isonormals(x,y,z,q, p)
set(p, 'FaceColor', 'interp', 'EdgeColor', 'none');
daspect([1 1 1]); axis tight;
>> k
k =
vertices: [2318x3 double]
faces: [4344x3 double]
facevertexcdata: [2318x1 double]
其中的vertices是顶点集合，facevertexcdata是面集合
上面程序加上：
figure,plot(k.vertices)
figure,plot(k. facevertexcdata)
看看效果。

Answer 3

举个例子
[x y z v] = flow;
[fc vt] = isosurface(x, y, z, v, -3);%vt就是曲面点集
p=patch('Faces',fc,'Vertices',vt);
isonormals(x,y,z,v, p)
set(p, 'FaceColor', 'red', 'EdgeColor', 'none');
daspect([1 1 1])
view(3)
camlight; lighting phong
figure
plot3(vt(:,1),vt(:,2),vt(:,3),'.')