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
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
|
package require vtk
package require vtkinteraction
package require vtktesting
# Contour every quadratic cell type
# Create a scene with one of each cell type.
# QuadraticEdge
vtkPoints edgePoints
edgePoints SetNumberOfPoints 3
edgePoints InsertPoint 0 0 0 0
edgePoints InsertPoint 1 1.0 0 0
edgePoints InsertPoint 2 0.5 0.25 0
vtkFloatArray edgeScalars
edgeScalars SetNumberOfTuples 3
edgeScalars InsertValue 0 0.0
edgeScalars InsertValue 1 0.0
edgeScalars InsertValue 2 0.9
vtkQuadraticEdge aEdge
[aEdge GetPointIds] SetId 0 0
[aEdge GetPointIds] SetId 1 1
[aEdge GetPointIds] SetId 2 2
vtkUnstructuredGrid aEdgeGrid
aEdgeGrid Allocate 1 1
aEdgeGrid InsertNextCell [aEdge GetCellType] [aEdge GetPointIds]
aEdgeGrid SetPoints edgePoints
[aEdgeGrid GetPointData] SetScalars edgeScalars
vtkDataSetMapper aEdgeMapper
aEdgeMapper SetInput aEdgeGrid
aEdgeMapper ScalarVisibilityOff
vtkActor aEdgeActor
aEdgeActor SetMapper aEdgeMapper
[aEdgeActor GetProperty] SetRepresentationToWireframe
[aEdgeActor GetProperty] SetAmbient 1.0
# Quadratic triangle
vtkPoints triPoints
triPoints SetNumberOfPoints 6
triPoints InsertPoint 0 0.0 0.0 0.0
triPoints InsertPoint 1 1.0 0.0 0.0
triPoints InsertPoint 2 0.5 0.8 0.0
triPoints InsertPoint 3 0.5 0.0 0.0
triPoints InsertPoint 4 0.75 0.4 0.0
triPoints InsertPoint 5 0.25 0.4 0.0
vtkFloatArray triScalars
triScalars SetNumberOfTuples 6
triScalars InsertValue 0 0.0
triScalars InsertValue 1 0.0
triScalars InsertValue 2 0.0
triScalars InsertValue 3 1.0
triScalars InsertValue 4 0.0
triScalars InsertValue 5 0.0
vtkQuadraticTriangle aTri
[aTri GetPointIds] SetId 0 0
[aTri GetPointIds] SetId 1 1
[aTri GetPointIds] SetId 2 2
[aTri GetPointIds] SetId 3 3
[aTri GetPointIds] SetId 4 4
[aTri GetPointIds] SetId 5 5
vtkUnstructuredGrid aTriGrid
aTriGrid Allocate 1 1
aTriGrid InsertNextCell [aTri GetCellType] [aTri GetPointIds]
aTriGrid SetPoints triPoints
[aTriGrid GetPointData] SetScalars triScalars
vtkDataSetMapper aTriMapper
aTriMapper SetInput aTriGrid
aTriMapper ScalarVisibilityOff
vtkActor aTriActor
aTriActor SetMapper aTriMapper
[aTriActor GetProperty] SetRepresentationToWireframe
[aTriActor GetProperty] SetAmbient 1.0
# Quadratic quadrilateral
vtkPoints quadPoints
quadPoints SetNumberOfPoints 8
quadPoints InsertPoint 0 0.0 0.0 0.0
quadPoints InsertPoint 1 1.0 0.0 0.0
quadPoints InsertPoint 2 1.0 1.0 0.0
quadPoints InsertPoint 3 0.0 1.0 0.0
quadPoints InsertPoint 4 0.5 0.0 0.0
quadPoints InsertPoint 5 1.0 0.5 0.0
quadPoints InsertPoint 6 0.5 1.0 0.0
quadPoints InsertPoint 7 0.0 0.5 0.0
vtkFloatArray quadScalars
quadScalars SetNumberOfTuples 8
quadScalars InsertValue 0 0.0
quadScalars InsertValue 1 0.0
quadScalars InsertValue 2 1.0
quadScalars InsertValue 3 1.0
quadScalars InsertValue 4 1.0
quadScalars InsertValue 5 0.0
quadScalars InsertValue 6 0.0
quadScalars InsertValue 7 0.0
vtkQuadraticQuad aQuad
[aQuad GetPointIds] SetId 0 0
[aQuad GetPointIds] SetId 1 1
[aQuad GetPointIds] SetId 2 2
[aQuad GetPointIds] SetId 3 3
[aQuad GetPointIds] SetId 4 4
[aQuad GetPointIds] SetId 5 5
[aQuad GetPointIds] SetId 6 6
[aQuad GetPointIds] SetId 7 7
vtkUnstructuredGrid aQuadGrid
aQuadGrid Allocate 1 1
aQuadGrid InsertNextCell [aQuad GetCellType] [aQuad GetPointIds]
aQuadGrid SetPoints quadPoints
[aQuadGrid GetPointData] SetScalars quadScalars
vtkDataSetMapper aQuadMapper
aQuadMapper SetInput aQuadGrid
aQuadMapper ScalarVisibilityOff
vtkActor aQuadActor
aQuadActor SetMapper aQuadMapper
[aQuadActor GetProperty] SetRepresentationToWireframe
[aQuadActor GetProperty] SetAmbient 1.0
# Quadratic tetrahedron
vtkPoints tetPoints
tetPoints SetNumberOfPoints 10
tetPoints InsertPoint 0 0.0 0.0 0.0
tetPoints InsertPoint 1 1.0 0.0 0.0
tetPoints InsertPoint 2 0.5 0.8 0.0
tetPoints InsertPoint 3 0.5 0.4 1.0
tetPoints InsertPoint 4 0.5 0.0 0.0
tetPoints InsertPoint 5 0.75 0.4 0.0
tetPoints InsertPoint 6 0.25 0.4 0.0
tetPoints InsertPoint 7 0.25 0.2 0.5
tetPoints InsertPoint 8 0.75 0.2 0.5
tetPoints InsertPoint 9 0.50 0.6 0.5
vtkFloatArray tetScalars
tetScalars SetNumberOfTuples 10
tetScalars InsertValue 0 1.0
tetScalars InsertValue 1 1.0
tetScalars InsertValue 2 1.0
tetScalars InsertValue 3 1.0
tetScalars InsertValue 4 0.0
tetScalars InsertValue 5 0.0
tetScalars InsertValue 6 0.0
tetScalars InsertValue 7 0.0
tetScalars InsertValue 8 0.0
tetScalars InsertValue 9 0.0
vtkQuadraticTetra aTet
[aTet GetPointIds] SetId 0 0
[aTet GetPointIds] SetId 1 1
[aTet GetPointIds] SetId 2 2
[aTet GetPointIds] SetId 3 3
[aTet GetPointIds] SetId 4 4
[aTet GetPointIds] SetId 5 5
[aTet GetPointIds] SetId 6 6
[aTet GetPointIds] SetId 7 7
[aTet GetPointIds] SetId 8 8
[aTet GetPointIds] SetId 9 9
vtkUnstructuredGrid aTetGrid
aTetGrid Allocate 1 1
aTetGrid InsertNextCell [aTet GetCellType] [aTet GetPointIds]
aTetGrid SetPoints tetPoints
[aTetGrid GetPointData] SetScalars tetScalars
vtkDataSetMapper aTetMapper
aTetMapper SetInput aTetGrid
aTetMapper ScalarVisibilityOff
vtkActor aTetActor
aTetActor SetMapper aTetMapper
[aTetActor GetProperty] SetRepresentationToWireframe
[aTetActor GetProperty] SetAmbient 1.0
# Quadratic hexahedron
vtkPoints hexPoints
hexPoints SetNumberOfPoints 20
hexPoints InsertPoint 0 0 0 0
hexPoints InsertPoint 1 1 0 0
hexPoints InsertPoint 2 1 1 0
hexPoints InsertPoint 3 0 1 0
hexPoints InsertPoint 4 0 0 1
hexPoints InsertPoint 5 1 0 1
hexPoints InsertPoint 6 1 1 1
hexPoints InsertPoint 7 0 1 1
hexPoints InsertPoint 8 0.5 0 0
hexPoints InsertPoint 9 1 0.5 0
hexPoints InsertPoint 10 0.5 1 0
hexPoints InsertPoint 11 0 0.5 0
hexPoints InsertPoint 12 0.5 0 1
hexPoints InsertPoint 13 1 0.5 1
hexPoints InsertPoint 14 0.5 1 1
hexPoints InsertPoint 15 0 0.5 1
hexPoints InsertPoint 16 0 0 0.5
hexPoints InsertPoint 17 1 0 0.5
hexPoints InsertPoint 18 1 1 0.5
hexPoints InsertPoint 19 0 1 0.5
vtkFloatArray hexScalars
hexScalars SetNumberOfTuples 20
hexScalars InsertValue 0 1.0
hexScalars InsertValue 1 1.0
hexScalars InsertValue 2 1.0
hexScalars InsertValue 3 1.0
hexScalars InsertValue 4 1.0
hexScalars InsertValue 5 1.0
hexScalars InsertValue 6 1.0
hexScalars InsertValue 7 1.0
hexScalars InsertValue 8 0.0
hexScalars InsertValue 9 0.0
hexScalars InsertValue 10 0.0
hexScalars InsertValue 11 0.0
hexScalars InsertValue 12 0.0
hexScalars InsertValue 13 0.0
hexScalars InsertValue 14 0.0
hexScalars InsertValue 15 0.0
hexScalars InsertValue 16 0.0
hexScalars InsertValue 17 0.0
hexScalars InsertValue 18 0.0
hexScalars InsertValue 19 0.0
vtkQuadraticHexahedron aHex
[aHex GetPointIds] SetId 0 0
[aHex GetPointIds] SetId 1 1
[aHex GetPointIds] SetId 2 2
[aHex GetPointIds] SetId 3 3
[aHex GetPointIds] SetId 4 4
[aHex GetPointIds] SetId 5 5
[aHex GetPointIds] SetId 6 6
[aHex GetPointIds] SetId 7 7
[aHex GetPointIds] SetId 8 8
[aHex GetPointIds] SetId 9 9
[aHex GetPointIds] SetId 10 10
[aHex GetPointIds] SetId 11 11
[aHex GetPointIds] SetId 12 12
[aHex GetPointIds] SetId 13 13
[aHex GetPointIds] SetId 14 14
[aHex GetPointIds] SetId 15 15
[aHex GetPointIds] SetId 16 16
[aHex GetPointIds] SetId 17 17
[aHex GetPointIds] SetId 18 18
[aHex GetPointIds] SetId 19 19
vtkUnstructuredGrid aHexGrid
aHexGrid Allocate 1 1
aHexGrid InsertNextCell [aHex GetCellType] [aHex GetPointIds]
aHexGrid SetPoints hexPoints
[aHexGrid GetPointData] SetScalars hexScalars
vtkDataSetMapper aHexMapper
aHexMapper SetInput aHexGrid
aHexMapper ScalarVisibilityOff
vtkActor aHexActor
aHexActor SetMapper aHexMapper
[aHexActor GetProperty] SetRepresentationToWireframe
[aHexActor GetProperty] SetAmbient 1.0
# Quadratic wedge
vtkPoints wedgePoints
wedgePoints SetNumberOfPoints 15
wedgePoints InsertPoint 0 0 0 0
wedgePoints InsertPoint 1 1 0 0
wedgePoints InsertPoint 2 0 1 0
wedgePoints InsertPoint 3 0 0 1
wedgePoints InsertPoint 4 1 0 1
wedgePoints InsertPoint 5 0 1 1
wedgePoints InsertPoint 6 0.5 0 0
wedgePoints InsertPoint 7 0.5 0.5 0
wedgePoints InsertPoint 8 0 0.5 0
wedgePoints InsertPoint 9 0.5 0 1
wedgePoints InsertPoint 10 0.5 0.5 1
wedgePoints InsertPoint 11 0 0.5 1
wedgePoints InsertPoint 12 0 0 0.5
wedgePoints InsertPoint 13 1 0 0.5
wedgePoints InsertPoint 14 0 1 0.5
vtkFloatArray wedgeScalars
wedgeScalars SetNumberOfTuples 15
wedgeScalars InsertValue 0 1.0
wedgeScalars InsertValue 1 1.0
wedgeScalars InsertValue 2 1.0
wedgeScalars InsertValue 3 1.0
wedgeScalars InsertValue 4 1.0
wedgeScalars InsertValue 5 1.0
wedgeScalars InsertValue 6 1.0
wedgeScalars InsertValue 7 1.0
wedgeScalars InsertValue 8 0.0
wedgeScalars InsertValue 9 0.0
wedgeScalars InsertValue 10 0.0
wedgeScalars InsertValue 11 0.0
wedgeScalars InsertValue 12 0.0
wedgeScalars InsertValue 13 0.0
wedgeScalars InsertValue 14 0.0
vtkQuadraticWedge aWedge
[aWedge GetPointIds] SetId 0 0
[aWedge GetPointIds] SetId 1 1
[aWedge GetPointIds] SetId 2 2
[aWedge GetPointIds] SetId 3 3
[aWedge GetPointIds] SetId 4 4
[aWedge GetPointIds] SetId 5 5
[aWedge GetPointIds] SetId 6 6
[aWedge GetPointIds] SetId 7 7
[aWedge GetPointIds] SetId 8 8
[aWedge GetPointIds] SetId 9 9
[aWedge GetPointIds] SetId 10 10
[aWedge GetPointIds] SetId 11 11
[aWedge GetPointIds] SetId 12 12
[aWedge GetPointIds] SetId 13 13
[aWedge GetPointIds] SetId 14 14
vtkUnstructuredGrid aWedgeGrid
aWedgeGrid Allocate 1 1
aWedgeGrid InsertNextCell [aWedge GetCellType] [aWedge GetPointIds]
aWedgeGrid SetPoints wedgePoints
[aWedgeGrid GetPointData] SetScalars wedgeScalars
vtkClipDataSet wedgeContours
wedgeContours SetInput aWedgeGrid
wedgeContours SetValue 0.5
vtkDataSetMapper aWedgeContourMapper
aWedgeContourMapper SetInputConnection [wedgeContours GetOutputPort]
aWedgeContourMapper ScalarVisibilityOff
vtkDataSetMapper aWedgeMapper
aWedgeMapper SetInput aWedgeGrid
aWedgeMapper ScalarVisibilityOff
vtkActor aWedgeActor
aWedgeActor SetMapper aWedgeMapper
[aWedgeActor GetProperty] SetRepresentationToWireframe
[aWedgeActor GetProperty] SetAmbient 1.0
vtkActor aWedgeContourActor
aWedgeContourActor SetMapper aWedgeContourMapper
[aWedgeContourActor GetProperty] SetAmbient 1.0
# Quadratic pyramid
vtkPoints pyraPoints
pyraPoints SetNumberOfPoints 13
pyraPoints InsertPoint 0 0 0 0
pyraPoints InsertPoint 1 1 0 0
pyraPoints InsertPoint 2 1 1 0
pyraPoints InsertPoint 3 0 1 0
pyraPoints InsertPoint 4 0 0 1
pyraPoints InsertPoint 5 0.5 0 0
pyraPoints InsertPoint 6 1 0.5 0
pyraPoints InsertPoint 7 0.5 1 0
pyraPoints InsertPoint 8 0 0.5 0
pyraPoints InsertPoint 9 0 0 0.5
pyraPoints InsertPoint 10 0.5 0 0.5
pyraPoints InsertPoint 11 0.5 0.5 0.5
pyraPoints InsertPoint 12 0 0.5 0.5
vtkFloatArray pyraScalars
pyraScalars SetNumberOfTuples 13
pyraScalars InsertValue 0 1.0
pyraScalars InsertValue 1 1.0
pyraScalars InsertValue 2 1.0
pyraScalars InsertValue 3 1.0
pyraScalars InsertValue 4 1.0
pyraScalars InsertValue 5 1.0
pyraScalars InsertValue 6 1.0
pyraScalars InsertValue 7 1.0
pyraScalars InsertValue 8 0.0
pyraScalars InsertValue 9 0.0
pyraScalars InsertValue 10 0.0
pyraScalars InsertValue 11 0.0
pyraScalars InsertValue 12 0.0
vtkQuadraticPyramid aPyramid
[aPyramid GetPointIds] SetId 0 0
[aPyramid GetPointIds] SetId 1 1
[aPyramid GetPointIds] SetId 2 2
[aPyramid GetPointIds] SetId 3 3
[aPyramid GetPointIds] SetId 4 4
[aPyramid GetPointIds] SetId 5 5
[aPyramid GetPointIds] SetId 6 6
[aPyramid GetPointIds] SetId 7 7
[aPyramid GetPointIds] SetId 8 8
[aPyramid GetPointIds] SetId 9 9
[aPyramid GetPointIds] SetId 10 10
[aPyramid GetPointIds] SetId 11 11
[aPyramid GetPointIds] SetId 12 12
vtkUnstructuredGrid aPyramidGrid
aPyramidGrid Allocate 1 1
aPyramidGrid InsertNextCell [aPyramid GetCellType] [aPyramid GetPointIds]
aPyramidGrid SetPoints pyraPoints
[aPyramidGrid GetPointData] SetScalars pyraScalars
vtkClipDataSet pyraContours
pyraContours SetInput aPyramidGrid
pyraContours SetValue 0.5
vtkDataSetMapper aPyramidContourMapper
aPyramidContourMapper SetInputConnection [pyraContours GetOutputPort]
aPyramidContourMapper ScalarVisibilityOff
vtkDataSetMapper aPyramidMapper
aPyramidMapper SetInput aPyramidGrid
aPyramidMapper ScalarVisibilityOff
vtkActor aPyramidActor
aPyramidActor SetMapper aPyramidMapper
[aPyramidActor GetProperty] SetRepresentationToWireframe
[aPyramidActor GetProperty] SetAmbient 1.0
vtkActor aPyramidContourActor
aPyramidContourActor SetMapper aPyramidContourMapper
[aPyramidContourActor GetProperty] SetAmbient 1.0
# Create the rendering related stuff.
# Since some of our actors are a single vertex, we need to remove all
# cullers so the single vertex actors will render
vtkRenderer ren1
[ren1 GetCullers] RemoveAllItems
vtkRenderWindow renWin
renWin AddRenderer ren1
vtkRenderWindowInteractor iren
iren SetRenderWindow renWin
ren1 SetBackground .1 .2 .3
renWin SetSize 400 200
# specify properties
ren1 AddActor aEdgeActor
ren1 AddActor aTriActor
ren1 AddActor aQuadActor
ren1 AddActor aTetActor
ren1 AddActor aHexActor
ren1 AddActor aWedgeActor
ren1 AddActor aPyramidActor
# places everyone!!
aTriActor AddPosition 2 0 0
aQuadActor AddPosition 4 0 0
aTetActor AddPosition 6 0 0
aHexActor AddPosition 8 0 0
aWedgeActor AddPosition 10 0 0
aPyramidActor AddPosition 12 0 0
BuildBackdrop -1 15 -1 4 -1 2 .1
ren1 AddActor base
[base GetProperty] SetDiffuseColor .2 .2 .2
ren1 AddActor left
[left GetProperty] SetDiffuseColor .2 .2 .2
ren1 AddActor back
[back GetProperty] SetDiffuseColor .2 .2 .2
ren1 ResetCamera
[ren1 GetActiveCamera] Dolly 2.5
ren1 ResetCameraClippingRange
renWin Render
# create a little scorecard above each of the cells. These are displayed
# if a ray cast hits the cell, otherwise they are not shown.
vtkPlaneSource pm
pm SetXResolution 1
pm SetYResolution 1
vtkPolyDataMapper pmapper
pmapper SetInputConnection [pm GetOutputPort]
# now try intersecting rays with the cell
vtkCellPicker cellPicker
vtkActor edgeCheck
edgeCheck SetMapper pmapper
edgeCheck AddPosition 0.5 2.5 0
cellPicker Pick 87 71 0 ren1
if {[cellPicker GetCellId] != "-1" } {
ren1 AddActor edgeCheck
}
vtkActor triCheck
triCheck SetMapper pmapper
triCheck AddPosition 2.5 2.5 0
cellPicker Pick 139 72 0 ren1
if {[cellPicker GetCellId] != "-1" } {
ren1 AddActor triCheck
}
vtkActor quadCheck
quadCheck SetMapper pmapper
quadCheck AddPosition 4.5 2.5 0
cellPicker Pick 192 78 0 ren1
if {[cellPicker GetCellId] != "-1" } {
ren1 AddActor quadCheck
}
vtkActor tetCheck
tetCheck SetMapper pmapper
tetCheck AddPosition 6.5 2.5 0
cellPicker Pick 233 70 0 ren1
if {[cellPicker GetCellId] != "-1" } {
ren1 AddActor tetCheck
}
vtkActor hexCheck
hexCheck SetMapper pmapper
hexCheck AddPosition 8.5 2.5 0
cellPicker Pick 287 80 0 ren1
if {[cellPicker GetCellId] != "-1" } {
ren1 AddActor hexCheck
}
vtkActor wedgeCheck
wedgeCheck SetMapper pmapper
wedgeCheck AddPosition 10.5 2.5 0
cellPicker Pick 287 80 0 ren1
if {[cellPicker GetCellId] != "-1" } {
ren1 AddActor wedgeCheck
}
vtkActor pyraCheck
pyraCheck SetMapper pmapper
pyraCheck AddPosition 12.5 2.5 0
cellPicker Pick 287 80 0 ren1
if {[cellPicker GetCellId] != "-1" } {
ren1 AddActor pyraCheck
}
# render the image
#
iren AddObserver UserEvent {wm deiconify .vtkInteract}
iren Initialize
wm withdraw .
|
