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
|
package require vtk
package require vtkinteraction
package require vtktesting
#define a Single Cube
vtkFloatArray Scalars
Scalars InsertNextValue 1.0
Scalars InsertNextValue 0.0
Scalars InsertNextValue 0.0
Scalars InsertNextValue 1.0
vtkPoints Points
Points InsertNextPoint 0 0 0
Points InsertNextPoint 1 0 0
Points InsertNextPoint 0 1 0
Points InsertNextPoint 0 0 1
vtkIdList Ids
Ids InsertNextId 0
Ids InsertNextId 1
Ids InsertNextId 2
Ids InsertNextId 3
vtkUnstructuredGrid Grid
Grid Allocate 10 10
Grid InsertNextCell 10 Ids
Grid SetPoints Points
[Grid GetPointData] SetScalars Scalars
#Clip the tetra
vtkClipDataSet clipper
clipper SetInputData Grid
clipper SetValue 0.5
clipper Update
# build tubes for the triangle edges
#
vtkExtractEdges tetEdges
tetEdges SetInputConnection [clipper GetOutputPort]
vtkTubeFilter tetEdgeTubes
tetEdgeTubes SetInputConnection [tetEdges GetOutputPort]
tetEdgeTubes SetRadius .005
tetEdgeTubes SetNumberOfSides 6
vtkPolyDataMapper tetEdgeMapper
tetEdgeMapper SetInputConnection [tetEdgeTubes GetOutputPort]
tetEdgeMapper ScalarVisibilityOff
vtkActor tetEdgeActor
tetEdgeActor SetMapper tetEdgeMapper
eval [tetEdgeActor GetProperty] SetDiffuseColor $lamp_black
[tetEdgeActor GetProperty] SetSpecular .4
[tetEdgeActor GetProperty] SetSpecularPower 10
#shrink the triangles so we can see each one
vtkShrinkFilter aShrinker
aShrinker SetShrinkFactor 1
aShrinker SetInputConnection [clipper GetOutputPort]
vtkDataSetMapper aMapper
aMapper ScalarVisibilityOff
aMapper SetInputConnection [aShrinker GetOutputPort]
vtkActor Tets
Tets SetMapper aMapper
eval [Tets GetProperty] SetDiffuseColor $banana
#build a model of the cube
vtkExtractEdges Edges
Edges SetInputData Grid
vtkTubeFilter Tubes
Tubes SetInputConnection [Edges GetOutputPort]
Tubes SetRadius .01
Tubes SetNumberOfSides 6
vtkPolyDataMapper TubeMapper
TubeMapper SetInputConnection [Tubes GetOutputPort]
TubeMapper ScalarVisibilityOff
vtkActor CubeEdges
CubeEdges SetMapper TubeMapper
eval [CubeEdges GetProperty] SetDiffuseColor $khaki
[CubeEdges GetProperty] SetSpecular .4
[CubeEdges GetProperty] SetSpecularPower 10
# build the vertices of the cube
#
vtkSphereSource Sphere
Sphere SetRadius 0.04
Sphere SetPhiResolution 20
Sphere SetThetaResolution 20
vtkThresholdPoints ThresholdIn
ThresholdIn SetInputData Grid
ThresholdIn ThresholdByUpper .5
vtkGlyph3D Vertices
Vertices SetInputConnection [ThresholdIn GetOutputPort]
Vertices SetSourceConnection [Sphere GetOutputPort]
vtkPolyDataMapper SphereMapper
SphereMapper SetInputConnection [Vertices GetOutputPort]
SphereMapper ScalarVisibilityOff
vtkActor CubeVertices
CubeVertices SetMapper SphereMapper
eval [CubeVertices GetProperty] SetDiffuseColor $tomato
eval [CubeVertices GetProperty] SetDiffuseColor $tomato
#define the text for the labels
vtkVectorText caseLabel
caseLabel SetText "Case 1"
vtkTransform aLabelTransform
aLabelTransform Identity
aLabelTransform Translate -.2 0 1.25
aLabelTransform Scale .05 .05 .05
vtkTransformPolyDataFilter labelTransform
labelTransform SetTransform aLabelTransform
labelTransform SetInputConnection [caseLabel GetOutputPort]
vtkPolyDataMapper labelMapper
labelMapper SetInputConnection [labelTransform GetOutputPort];
vtkActor labelActor
labelActor SetMapper labelMapper
#define the base
vtkCubeSource baseModel
baseModel SetXLength 1.5
baseModel SetYLength .01
baseModel SetZLength 1.5
vtkPolyDataMapper baseMapper
baseMapper SetInputConnection [baseModel GetOutputPort]
vtkActor base
base SetMapper baseMapper
# Create the RenderWindow, Renderer and both Actors
#
vtkRenderer ren1
vtkRenderWindow renWin
renWin AddRenderer ren1
vtkRenderWindowInteractor iren
iren SetRenderWindow renWin
# position the base
base SetPosition .5 -.09 .5
ren1 AddActor tetEdgeActor
ren1 AddActor base
ren1 AddActor labelActor
ren1 AddActor CubeEdges
ren1 AddActor CubeVertices
ren1 AddActor Tets
eval ren1 SetBackground $slate_grey
iren AddObserver UserEvent {wm deiconify .vtkInteract}
Grid Modified
renWin SetSize 400 400
ren1 ResetCamera
[ren1 GetActiveCamera] Dolly 1.2
[ren1 GetActiveCamera] Azimuth 30
[ren1 GetActiveCamera] Elevation 20
ren1 ResetCameraClippingRange
renWin Render
iren Initialize
set mask "1 2 4 8 16 32"
proc cases {id} {
global mask
for {set i 0} {$i < 4} {incr i} {
set m [lindex $mask $i]
if {[expr $m & $id] == 0} {
Scalars SetValue $i 0
} else {
Scalars SetValue $i 1
}
caseLabel SetText "Case $id"
}
Grid Modified
renWin Render
}
cases 3
wm withdraw .
|