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.
Example 1:
[x y z v] = flow;
p = patch(isosurface(x, y, z, v, -3));
isonormals(x,y,z,v, p)
p.FaceColor = 'red';
p.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)
p.FaceColor = 'interp';
p.EdgeColor = 'none';
daspect([1 1 1]); axis tight;
colormap(prism(28))
camup([1 0 0 ]); campos([25 -55 5])
camlight; lighting phong
See also ISONORMALS, ISOCAPS, SMOOTH3, SUBVOLUME, REDUCEVOLUME,
REDUCEPATCH, SHRINKFACES.
Documentation for isosurface
doc isosurface
