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
|
package require vtk
package require vtkinteraction
set res 6
vtkPlaneSource plane
plane SetResolution $res $res
vtkElevationFilter colors
colors SetInputConnection [plane GetOutputPort]
colors SetLowPoint -0.25 -0.25 -0.25
colors SetHighPoint 0.25 0.25 0.25
vtkPolyDataMapper planeMapper
planeMapper SetInputConnection [colors GetOutputPort]
vtkActor planeActor
planeActor SetMapper planeMapper
[planeActor GetProperty] SetRepresentationToWireframe
# create simple poly data so we can apply glyph
vtkSuperquadricSource squad
vtkElevationFilter squadColors
squadColors SetInputConnection [squad GetOutputPort]
squadColors SetLowPoint -0.25 -0.25 -0.25
squadColors SetHighPoint 0.25 0.25 0.25
vtkCastToConcrete squadCaster
squadCaster SetInputConnection [squadColors GetOutputPort]
vtkTransform squadTransform
vtkTransformPolyDataFilter transformSquad
transformSquad SetInputConnection [squadColors GetOutputPort]
transformSquad SetTransform squadTransform
transformSquad Update
# procedure for generating glyphs
proc Glyph {} {
global res
set ptId [glypher GetPointId]
set pd [glypher GetPointData]
set xyz [glypher GetPoint]
set x [lindex $xyz 0]
set y [lindex $xyz 1]
set length [[glypher GetInput 0] GetLength]
set scale [expr $length / (2.0*$res)]
squadTransform Identity
if { $x == $y } {
squad ToroidalOn
eval squadTransform Translate $xyz
squadTransform RotateX 90
} else {
eval squadTransform Translate $xyz
squad ToroidalOff
}
squadTransform Scale $scale $scale $scale
squad SetPhiRoundness [expr abs($x)*5.0]
squad SetThetaRoundness [expr abs($y)*5.0]
}
vtkProgrammableGlyphFilter glypher
glypher SetInputConnection [colors GetOutputPort]
glypher SetSourceConnection [transformSquad GetOutputPort]
glypher SetGlyphMethod Glyph
glypher SetColorModeToColorBySource
vtkPolyDataMapper glyphMapper
glyphMapper SetInputConnection [glypher GetOutputPort]
vtkActor glyphActor
glyphActor SetMapper glyphMapper
# Create the rendering stuff
#
vtkRenderer ren1
vtkRenderWindow renWin
renWin SetMultiSamples 0
renWin AddRenderer ren1
vtkRenderWindowInteractor iren
iren SetRenderWindow renWin
ren1 AddActor planeActor
ren1 AddActor glyphActor
ren1 SetBackground 1 1 1
renWin SetSize 450 450
renWin Render
[ren1 GetActiveCamera] Zoom 1.3
# Get handles to some useful objects
#
iren AddObserver UserEvent {wm deiconify .vtkInteract}
renWin Render
# prevent the tk window from showing up then start the event loop
wm withdraw .
|
