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
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
|
package require vtk
package require vtkinteraction
# create pipeline
#
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 output [[pl3d GetOutput] GetBlock 0]
# create three line probes
vtkLineSource line
line SetResolution 30
vtkTransform transL1
transL1 Translate 3.7 0.0 28.37
transL1 Scale 5 5 5
transL1 RotateY 90
vtkTransformPolyDataFilter tf
tf SetInputConnection [line GetOutputPort]
tf SetTransform transL1
vtkProbeFilter probe
probe SetInputConnection [tf GetOutputPort]
probe SetSourceData $output
probe Update
vtkTransform transL2
transL2 Translate 9.2 0.0 31.20
transL2 Scale 5 5 5
transL2 RotateY 90
vtkTransformPolyDataFilter tf2
tf2 SetInputConnection [line GetOutputPort]
tf2 SetTransform transL2
vtkProbeFilter probe2
probe2 SetInputConnection [tf2 GetOutputPort]
probe2 SetSourceData $output
probe2 Update
vtkTransform transL3
transL3 Translate 13.27 0.0 33.40
transL3 Scale 4.5 4.5 4.5
transL3 RotateY 90
vtkTransformPolyDataFilter tf3
tf3 SetInputConnection [line GetOutputPort]
tf3 SetTransform transL3
vtkProbeFilter probe3
probe3 SetInputConnection [tf3 GetOutputPort]
probe3 SetSourceData $output
probe3 Update
vtkAppendPolyData appendF
appendF AddInputConnection [probe GetOutputPort]
appendF AddInputConnection [probe2 GetOutputPort]
appendF AddInputConnection [probe3 GetOutputPort]
vtkTubeFilter tuber
tuber SetInputConnection [appendF GetOutputPort]
tuber SetRadius 0.1
vtkPolyDataMapper lineMapper
lineMapper SetInputConnection [tuber GetOutputPort]
vtkActor lineActor
lineActor SetMapper lineMapper
# probe the line and plot it
vtkGlyphSource2D triangle
triangle SetGlyphTypeToTriangle
triangle Update
vtkGlyphSource2D cross
cross SetGlyphTypeToCross
cross Update
vtkXYPlotActor xyplot
xyplot AddDataSetInputConnection [probe GetOutputPort]
xyplot AddDataSetInputConnection [probe2 GetOutputPort]
xyplot AddDataSetInputConnection [probe3 GetOutputPort]
[xyplot GetPositionCoordinate] SetValue 0.0 0.67 0
[xyplot GetPosition2Coordinate] SetValue 1.0 0.33 0;#relative to Position
xyplot SetXValuesToArcLength
xyplot SetNumberOfXLabels 6
xyplot SetTitle "Pressure vs. Arc Length (Zoomed View)"
xyplot SetXTitle ""
xyplot SetYTitle "P"
xyplot SetXRange .1 .35
xyplot SetYRange .2 .4
[xyplot GetProperty] SetColor 0 0 0
xyplot PlotLinesOn
[xyplot GetProperty] SetLineWidth 1
xyplot PlotPointsOn
[xyplot GetProperty] SetPointSize 3
xyplot LegendOn
xyplot SetPlotSymbol 2 [triangle GetOutput]
xyplot SetPlotColor 2 0 0 1
xyplot SetPlotColor 1 0 0 0
xyplot SetPlotColor 0 0 0 0
xyplot SetGlyphSize 0.025
# Set text prop color (same color for backward compat with test)
# Assign same object to all text props
set tprop [xyplot GetTitleTextProperty]
eval $tprop SetColor [[xyplot GetProperty] GetColor]
xyplot SetAxisTitleTextProperty $tprop
xyplot SetAxisLabelTextProperty $tprop
xyplot SetLabelFormat "%-#6.2f"
vtkSphereSource vertexGlyph
vertexGlyph Update
vtkXYPlotActor xyplot2
xyplot2 AddDataSetInputConnection [probe GetOutputPort]
xyplot2 AddDataSetInputConnection [probe2 GetOutputPort]
xyplot2 AddDataSetInputConnection [probe3 GetOutputPort]
[xyplot2 GetPositionCoordinate] SetValue 0.00 0.33 0
[xyplot2 GetPosition2Coordinate] SetValue 1.0 0.33 0;#relative to Position
xyplot2 SetXValuesToNormalizedArcLength
xyplot2 SetNumberOfXLabels 6
xyplot2 SetXRange 0.2 1.0
xyplot2 SetTitle "VTK Quality vs. Alcohol Consumption"
xyplot2 SetXTitle ""
xyplot2 SetYTitle "A"
xyplot2 PlotPointsOn
xyplot2 PlotLinesOff
xyplot2 LegendOn
xyplot2 SetLegendPosition 0.4 0.6
xyplot2 SetLegendPosition2 0.40 0.25
[xyplot2 GetProperty] SetColor 1 0 0
[xyplot2 GetProperty] SetPointSize 2
xyplot2 SetPlotSymbol 0 [vertexGlyph GetOutput]
xyplot2 SetPlotLabel 0 "Ken's Mudslide Consumption"
xyplot2 SetPlotColor 0 1 0 0
xyplot2 SetPlotSymbol 1 [cross GetOutput]
xyplot2 SetPlotColor 1 1 0 1
xyplot2 SetPlotLabel 1 "Bill's Beer Consumption"
xyplot2 SetPlotSymbol 2 [triangle GetOutput]
xyplot2 SetPlotColor 2 0 0 1
xyplot2 SetPlotLabel 2 "VTK Quality"
# Set text prop color (same color for backward compat with test)
# Assign same object to all text props
set tprop [xyplot2 GetTitleTextProperty]
eval $tprop SetColor [[xyplot2 GetProperty] GetColor]
xyplot2 SetAxisTitleTextProperty $tprop
xyplot2 SetAxisLabelTextProperty $tprop
xyplot2 SetLabelFormat [xyplot GetLabelFormat]
#Okay exercise data object stuff
vtkDataSetToDataObjectFilter ds2do
ds2do SetInputConnection [probe GetOutputPort]
vtkDataSetToDataObjectFilter ds2do2
ds2do2 SetInputConnection [probe GetOutputPort]
ds2do2 Update
vtkDataSetToDataObjectFilter ds2do3
ds2do3 SetInputConnection [probe GetOutputPort]
ds2do3 Update
vtkXYPlotActor xyplot3
xyplot3 AddDataObjectInputConnection [ds2do GetOutputPort]
xyplot3 SetDataObjectXComponent 0 2
xyplot3 SetDataObjectYComponent 0 5
xyplot3 SetPlotColor 0 1 0 0
xyplot3 SetPlotLabel 0 "Mx"
xyplot3 AddDataObjectInput [ds2do2 GetOutput]
xyplot3 SetDataObjectXComponent 1 2
xyplot3 SetDataObjectYComponent 1 6
xyplot3 SetPlotColor 1 0 1 0
xyplot3 SetPlotLabel 1 "My"
xyplot3 AddDataObjectInput [ds2do3 GetOutput]
xyplot3 SetDataObjectXComponent 2 2
xyplot3 SetDataObjectYComponent 2 7
xyplot3 SetPlotColor 2 0 0 1
xyplot3 SetPlotLabel 2 "Mz"
[xyplot3 GetPositionCoordinate] SetValue 0.0 0.0 0
[xyplot3 GetPosition2Coordinate] SetValue 1.0 0.33 0;#relative to Position
xyplot3 SetXValuesToIndex
xyplot3 SetNumberOfXLabels 6
xyplot3 SetTitle "Momentum Component vs. Point Id"
xyplot3 SetXTitle "Point Id"
xyplot3 SetYTitle "M"
xyplot3 PlotPointsOn
[xyplot3 GetProperty] SetColor 0 0 1
[xyplot3 GetProperty] SetPointSize 3
xyplot3 LegendOn
xyplot3 SetLegendPosition 0.8 0.28
xyplot3 SetLegendPosition2 0.20 0.20
# Set text prop color (same color for backward compat with test)
# Assign same object to all text props
set tprop [xyplot3 GetTitleTextProperty]
eval $tprop SetColor [[xyplot3 GetProperty] GetColor]
xyplot3 SetAxisTitleTextProperty $tprop
xyplot3 SetAxisLabelTextProperty $tprop
xyplot3 SetLabelFormat "%4.f"
# draw an outline
vtkStructuredGridOutlineFilter outline
outline SetInputData $output
vtkPolyDataMapper outlineMapper
outlineMapper SetInputConnection [outline GetOutputPort]
vtkActor outlineActor
outlineActor SetMapper outlineMapper
[outlineActor GetProperty] SetColor 0 0 0
# Create graphics stuff
#
vtkRenderer ren1
vtkRenderer ren2
vtkRenderWindow renWin
renWin SetMultiSamples 0
renWin AddRenderer ren1
renWin AddRenderer ren2
vtkRenderWindowInteractor iren
iren SetRenderWindow renWin
ren1 SetBackground 0.6784 0.8471 0.9020
ren1 SetViewport 0 0 .5 1
ren1 AddActor outlineActor
ren1 AddActor lineActor
ren2 SetBackground 1 1 1
ren2 SetViewport 0.5 0.0 1.0 1.0
ren2 AddActor2D xyplot
ren2 AddActor2D xyplot2
ren2 AddActor2D xyplot3
renWin SetSize 790 400
set cam1 [ren1 GetActiveCamera]
$cam1 SetClippingRange 3.95297 100
$cam1 SetFocalPoint 8.88908 0.595038 29.3342
$cam1 SetPosition -12.3332 31.7479 41.2387
$cam1 SetViewUp 0.060772 -0.319905 0.945498
iren Initialize
renWin Render
# render the image
#
iren AddObserver UserEvent {wm deiconify .vtkInteract}
# prevent the tk window from showing up then start the event loop
wm withdraw .
|
