Stream tubes are similar to stream lines, except the tubes have width, providing another dimension that you can use to represent information.
By default, MATLAB® graphics display the divergence of the vector field by the width of the tube. You can also define widths for each tube vertex and thereby map other data to width.
This example uses the following techniques:
Stream tubes to indicate flow direction and divergence of the vector field
wind data set
Slice planes colored to indicate the speed of the wind currents overlaid with contour line to enhance visibility
Inputs include the coordinates of the volume, vector field components, and starting locations for the stream tubes.
Load the data and calculate values needed to make the plots. These values include:
The location of the slice planes (maximum
y, and a value for the altitude)
x value for the start of the stream
The speed of the wind (magnitude of the vector field)
load wind xmin = min(x(:)); xmax = max(x(:)); ymin = min(y(:)); alt = 7.356; % z value for slice and streamtube plane wind_speed = sqrt(u.^2 + v.^2 + w.^2);
hslice = slice(x,y,z,wind_speed,xmax,ymin,alt); set(hslice,'FaceColor','interp','EdgeColor','none') colormap hsv(16)
Add contour lines (
contourslice) to the slice planes.
Adjust the contour interval so the lines match the color boundaries in the slice
caxis to get the current
Set the interpolation method used by
linear to match the default used by
color_lim = caxis; cont_intervals = linspace(color_lim(1),color_lim(2),17); hcont = contourslice(x,y,z,wind_speed,xmax,ymin,... alt,cont_intervals,'linear'); set(hcont,'EdgeColor',[.4 .4 .4],'LineWidth',1)
meshgrid to create arrays for the
starting points for the stream tubes, which begin at the minimum
x value, range from 20 to 50 in
y, and lie
in a single plane in
z (corresponding to one of the slice
The stream tubes (
streamtube) are drawn at the
specified locations and scaled to be 1.25 times the default width to emphasize the
variation in divergence (width). The second element in the vector [1.25 30]
specifies the number of points along the circumference of the tube (the default is
20). You might want to increase this value as the tube size increases, to maintain a
Set the data aspect ratio (
daspect) before calling
Stream tubes are surface objects, therefore you can control their appearance by setting surface properties. This example sets surface properties to give a brightly lit, red surface.
[sx,sy,sz] = meshgrid(xmin,20:3:50,alt); daspect([1,1,1]) % set DAR before calling streamtube htubes = streamtube(x,y,z,u,v,w,sx,sy,sz,[1.25 30]); set(htubes,'EdgeColor','none','FaceColor','r',... 'AmbientStrength',.5)
view(-100,30) axis(volumebounds(x,y,z,wind_speed)) set(gca,'Projection','perspective') camlight left