1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
|
# This example demonstrates the generation of a streamsurface.
#
# First we include the VTK Tcl packages which will make available
# all of the vtk commands from Tcl. The vtkinteraction package defines
# a simple Tcl/Tk interactor widget.
#
package require vtk
package require vtkinteraction
package require vtktesting
# Read the data and specify which scalars and vectors to read.
#
vtkMultiBlockPLOT3DReader pl3d
pl3d SetXYZFileName "$VTK_DATA_ROOT/Data/combxyz.bin"
pl3d SetQFileName "$VTK_DATA_ROOT/Data/combq.bin"
pl3d SetScalarFunctionNumber 100
pl3d SetVectorFunctionNumber 202
pl3d Update
set pl3dOutput [[pl3d GetOutput] GetBlock 0]
# We use a rake to generate a series of streamline starting points
# scattered along a line. Each point will generate a streamline. These
# streamlines are then fed to the vtkRuledSurfaceFilter which stitches
# the lines together to form a surface.
#
vtkLineSource rake
rake SetPoint1 15 -5 32
rake SetPoint2 15 5 32
rake SetResolution 21
vtkPolyDataMapper rakeMapper
rakeMapper SetInputConnection [rake GetOutputPort]
vtkActor rakeActor
rakeActor SetMapper rakeMapper
vtkRungeKutta4 integ
vtkStreamLine sl
sl SetInputData $pl3dOutput
sl SetSourceConnection [rake GetOutputPort]
sl SetIntegrator integ
sl SetMaximumPropagationTime 0.1
sl SetIntegrationStepLength 0.1
sl SetIntegrationDirectionToBackward
sl SetStepLength 0.001
#
# The ruled surface stiches together lines with triangle strips.
# Note the SetOnRatio method. It turns on every other strip that
# the filter generates (only when multiple lines are input).
#
vtkRuledSurfaceFilter scalarSurface
scalarSurface SetInputConnection [sl GetOutputPort]
scalarSurface SetOffset 0
scalarSurface SetOnRatio 2
scalarSurface PassLinesOn
scalarSurface SetRuledModeToPointWalk
scalarSurface SetDistanceFactor 30
vtkPolyDataMapper mapper
mapper SetInputConnection [scalarSurface GetOutputPort]
eval mapper SetScalarRange [$pl3dOutput GetScalarRange]
vtkActor actor
actor SetMapper mapper
# Put an outline around for context.
#
vtkStructuredGridOutlineFilter outline
outline SetInputData $pl3dOutput
vtkPolyDataMapper outlineMapper
outlineMapper SetInputConnection [outline GetOutputPort]
vtkActor outlineActor
outlineActor SetMapper outlineMapper
[outlineActor GetProperty] SetColor 0 0 0
# Now create the usual graphics stuff.
vtkRenderer ren
vtkRenderWindow renWin
renWin AddRenderer ren
vtkRenderWindowInteractor iren
iren SetRenderWindow renWin
ren AddActor rakeActor
ren AddActor actor
ren AddActor outlineActor
ren SetBackground 1 1 1
renWin SetSize 300 300
iren AddObserver UserEvent {wm deiconify .vtkInteract}
iren Initialize
# interact with data
wm withdraw .
iren Start
|
