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
|
package require vtk
package require vtkinteraction
# Now create the RenderWindow, Renderer and Interactor
#
vtkRenderer ren1
vtkRenderWindow renWin
renWin AddRenderer ren1
vtkRenderWindowInteractor iren
iren SetRenderWindow renWin
vtkMath math
set numberOfInputPoints 30
vtkCardinalSpline aSplineX
vtkCardinalSpline aSplineY
vtkCardinalSpline aSplineZ
# generate random points
vtkPoints inputPoints
for {set i 0} {$i < $numberOfInputPoints} {incr i 1} {
set x [math Random 0 1]
set y [math Random 0 1]
set z [math Random 0 1]
aSplineX AddPoint $i $x
aSplineY AddPoint $i $y
aSplineZ AddPoint $i $z
inputPoints InsertPoint $i $x $y $z
}
vtkPolyData inputData
inputData SetPoints inputPoints
vtkSphereSource balls
balls SetRadius .01
balls SetPhiResolution 10
balls SetThetaResolution 10
vtkGlyph3D glyphPoints
glyphPoints SetInput inputData
glyphPoints SetSource [balls GetOutput]
vtkPolyDataMapper glyphMapper
glyphMapper SetInputConnection [glyphPoints GetOutputPort]
vtkActor glyph
glyph SetMapper glyphMapper
eval [glyph GetProperty] SetDiffuseColor 1 0.4 0.4
[glyph GetProperty] SetSpecular .3
[glyph GetProperty] SetSpecularPower 30
ren1 AddActor glyph
# create a polyline
vtkPoints points
vtkPolyData profileData
set numberOfOutputPoints 400
set offset 1.0
proc fit {} {
global numberOfInputPoints numberOfOutputPoints offset
points Reset
for {set i 0} {$i< $numberOfOutputPoints} {incr i 1} {
set t [expr ( $numberOfInputPoints - $offset ) / ( $numberOfOutputPoints - 1) * $i]
points InsertPoint $i [aSplineX Evaluate $t] [aSplineY Evaluate $t] [aSplineZ Evaluate $t]
}
profileData Modified
}
fit
vtkCellArray lines
lines InsertNextCell $numberOfOutputPoints
for {set i 0} {$i < $numberOfOutputPoints} {incr i 1} {
lines InsertCellPoint $i
}
profileData SetPoints points
profileData SetLines lines
vtkTubeFilter profileTubes
profileTubes SetNumberOfSides 8
profileTubes SetInput profileData
profileTubes SetRadius .005
vtkPolyDataMapper profileMapper
profileMapper SetInputConnection [profileTubes GetOutputPort]
vtkActor profile
profile SetMapper profileMapper
eval [profile GetProperty] SetDiffuseColor 1 1 0.6
[profile GetProperty] SetSpecular .3
[profile GetProperty] SetSpecularPower 30
ren1 AddActor profile
ren1 ResetCamera
[ren1 GetActiveCamera] Dolly 1.5
ren1 ResetCameraClippingRange
renWin SetSize 400 400
# render the image
#
iren Initialize
iren AddObserver UserEvent {wm deiconify .vtkInteract}
proc opened {} {
global offset
set offset 1.0
aSplineX ClosedOff
aSplineY ClosedOff
aSplineZ ClosedOff
fit
renWin Render
}
proc varyLeft {} {
for {set left -1} { $left <= 1 } {set left [expr $left + .05]} {
aSplineX SetLeftValue $left
aSplineY SetLeftValue $left
aSplineZ SetLeftValue $left
fit
renWin Render
}
}
proc varyRight {} {
for {set right -1} { $right <= 1 } {set right [expr $right + .05]} {
aSplineX SetRightValue $right
aSplineY SetRightValue $right
aSplineZ SetRightValue $right
fit
renWin Render
}
}
proc constraint {value } {
aSplineX SetLeftConstraint $value
aSplineY SetLeftConstraint $value
aSplineZ SetLeftConstraint $value
aSplineX SetRightConstraint $value
aSplineY SetRightConstraint $value
aSplineZ SetRightConstraint $value
}
proc closed {} {
global offset
set offset 0.0
aSplineX ClosedOn
aSplineY ClosedOn
aSplineZ ClosedOn
fit
renWin Render
}
# prevent the tk window from showing up then start the event loop
wm withdraw .
|
