A fly-through is an effect created by moving the camera through three-dimensional space, giving the impression that you are flying along with the camera as if in an aircraft. You can fly through regions of a scene that might be otherwise obscured by objects in the scene or you can fly by a scene by keeping the camera focused on a particular point.
To accomplish these effects you move the camera along a particular path, the x-axis for example, in a series of steps. To produce a fly-through, move both the camera position and the camera target at the same time.
The following example makes use of the fly-though effect to view the interior of an isosurface drawn within a volume defined by a vector field of wind velocities. This data represents air currents over North America.
This example employs a number of visualization techniques. It uses
Isosurfaces and cone plots to illustrate the flow through the volume
Lighting to illuminate the isosurface and cones in the volume
Stream lines to define a path for the camera through the volume
Coordinated motion of the camera position, camera target, and light
The first step is to draw the isosurface and plot the air flow using cone plots.
Setting the data aspect ratio (
[1,1,1] before drawing the cone plot enables MATLAB® software to calculate the size of the cones correctly for the final
load wind wind_speed = sqrt(u.^2 + v.^2 + w.^2); figure p = patch(isosurface(x,y,z,wind_speed,35)); isonormals(x,y,z,wind_speed,p) p.FaceColor = [0.75,0.25,0.25]; p.EdgeColor = [0.6,0.4,0.4]; [f,vt] = reducepatch(isosurface(x,y,z,wind_speed,45),0.05); daspect([1,1,1]); hcone = coneplot(x,y,z,u,v,w,vt(:,1),vt(:,2),vt(:,3),2); hcone.FaceColor = 'blue'; hcone.EdgeColor = 'none';
You need to define viewing parameters to ensure the scene is displayed correctly:
Selecting a perspective projection provides the perception of depth as the
camera passes through the interior of the isosurface (
Setting the camera view angle to a fixed value prevents MATLAB from automatically adjusting the angle to encompass the entire
scene as well as zooming in the desired amount (
camproj perspective camva(25)
Positioning the light source at the camera location and modifying the reflectance characteristics of the isosurface and cones enhances the realism of the scene:
Creating a light source at the camera position provides a "headlight" that
moves along with the camera through the isosurface interior (
Setting the reflection properties of the isosurface gives the appearance
of a dark interior (
AmbientStrength set to 0.1) with highly reflective material
DiffuseStrength set to 1).
SpecularStrength of the cones to 1 makes
them highly reflective.
hlight = camlight('headlight'); p.AmbientStrength = 1; p.SpecularStrength = 1; p.DiffuseStrength = 1; hcone.SpecularStrength = 1; set(gcf,'Color','k') set(gca,'Color',[0,0,0.25])
gouraud lighting for smoother lighting:
Stream lines indicate the direction of flow in the vector field. This example uses the x-, y-, and z-coordinate data of a single stream line to map a path through the volume. The camera is then moved along this path. The tasks include
Create a stream line starting at the point
x = 80,
y = 30,
z = 11.
Get the x-, y-, and z-coordinate data of the stream line.
Delete the stream line (you could also use
stream3 to calculate the
stream line data without actually drawing the stream line).
hsline = streamline(x,y,z,u,v,w,80,30,11); xd = hsline.XData; yd = hsline.YData; zd = hsline.ZData; delete(hsline)
To create a fly-through, move the camera position and camera target along the same
path. In this example, the camera target is placed five elements further along the
x-axis than the camera. A small value is added to the
camera target x position to prevent the position of the camera and target from
becoming the same point if the condition
xd(n) = xd(n+5) should
Update the camera position and camera target so that they both move along the coordinates of the stream line.
Move the light along with the camera.
drawnow to display the
results of each move.
for i=1:length(xd)-5 campos([xd(i),yd(i),zd(i)]) camtarget([xd(i+5)+min(xd)/500,yd(i),zd(i)]) camlight(hlight,'headlight') drawnow end
coneplot for a fixed visualization
of the same data.