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
