isosurface

MATLAB 2013b版本中给出的isosurface的函数说明

isosurface Isosurface extractor.
　　FV = isosurface(X,Y,Z,V,ISOVALUE) computes isosurface geometry for
　　data V at isosurface value ISOVALUE. Arrays (X,Y,Z) specify the points
　　at which the data V is given. The struct FV contains the faces and
　　vertices of the isosurface and can be passed directly to the PATCH
　　command.
　　
　　FV = isosurface(V,ISOVALUE) assumes [X Y Z] = meshgrid(1:N, 1:M, 1:P)
　　where [M,N,P]=SIZE(V).
　　
　　FV = isosurface(X,Y,Z,V) or FV = isosurface(V) selects an
　　isosurface value automatically using the histogram of the
　　data.
　　
　　FVC = isosurface(..., COLORS) interpolates the array COLORS onto
　　the scalar field and returns the interpolated values in
　　facevertexcdata. The size of the COLORS array must be the same
　　as V.
　　
　　FV = isosurface(..., 'noshare') does not attempt to create
　　shared vertices. This is faster, but produces a larger set of
　　vertices.
　　
　　FV = isosurface(..., 'verbose') prints progress messages to the
　　command window as the computation progresses.
　　
　　[F, V] = isosurface(...) or [F, V, C] = isosurface(...)
　　returns the faces and vertices (and facevertexcdata) in
　　separate arrays instead of a struct.
　　
　　isosurface(...) With no output arguments, a patch is created
　　into the current axes with the computed faces and vertices.
　　If no current axes exists, a new axes will be created with
　　a 3-D view and appropriate lighting.

Example 1:
　　[x y z v] = flow;
　　p = patch(isosurface(x, y, z, v, -3));
　　isonormals(x,y,z,v, p)
　　set(p, 'FaceColor', 'red', 'EdgeColor', 'none');
　　daspect([1 1 1])
　　view(3)
　　camlight; lighting phong
　　Example 2:
　　[x y z v] = flow;
　　q = z./x.*y.^3;
　　p = patch(isosurface(x, y, z, q, -.08, v));
　　isonormals(x,y,z,q, p)
　　set(p, 'FaceColor', 'interp', 'EdgeColor', 'none');
　　daspect([1 1 1]); axis tight;
　　colormap(prism(28))
　　camup([1 0 0 ]); campos([25 -55 5])
　　camlight; lighting phong