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
|
package require vtk
# Simple volume rendering example.
vtkImageReader reader
reader SetDataByteOrderToLittleEndian
reader SetDataExtent 0 63 0 63 1 93
reader SetFilePrefix "$VTK_DATA_ROOT/Data/headsq/quarter"
reader SetDataMask 0x7fff
reader SetDataSpacing 2 2 1
reader SetDataScalarTypeToUnsignedShort
reader Update
vtkImageChangeInformation changeFilter
changeFilter SetInputConnection [reader GetOutputPort]
changeFilter SetOutputOrigin -63 -63 -46
# Create transfer functions for opacity and color
vtkPiecewiseFunction opacityTransferFunction
opacityTransferFunction AddPoint 600 0.0
opacityTransferFunction AddPoint 2000 1.0
vtkColorTransferFunction colorTransferFunction
colorTransferFunction ClampingOff
colorTransferFunction AddHSVPoint 0.0 0.01 1.0 1.0
colorTransferFunction AddHSVPoint 1000.0 0.50 1.0 1.0
colorTransferFunction AddHSVPoint 2000.0 0.99 1.0 1.0
colorTransferFunction SetColorSpaceToHSV
# Create properties, mappers, volume actors, and ray cast function
vtkVolumeProperty volumeProperty
volumeProperty SetColor colorTransferFunction
volumeProperty SetScalarOpacity opacityTransferFunction
vtkVolumeTextureMapper3D volumeMapper
volumeMapper SetInputConnection [changeFilter GetOutputPort]
volumeMapper SetSampleDistance 0.25
vtkVolume volume
volume SetMapper volumeMapper
volume SetProperty volumeProperty
# Create geometric sphere
vtkSphereSource sphereSource
sphereSource SetRadius 65
sphereSource SetThetaResolution 20
sphereSource SetPhiResolution 40
# Compute random scalars (colors) for each cell
vtkProgrammableAttributeDataFilter randomColors
randomColors SetInputConnection [sphereSource GetOutputPort]
randomColors SetExecuteMethod colorCells
proc colorCells {} {
vtkMath randomColorGenerator
set input [randomColors GetInput]
set output [randomColors GetOutput]
set numCells [$input GetNumberOfCells]
vtkFloatArray colors
colors SetNumberOfTuples $numCells
for {set i 0} {$i < $numCells} {incr i} {
colors SetValue $i [randomColorGenerator Random 0 1]
}
[$output GetCellData] CopyScalarsOff
[$output GetCellData] PassData [$input GetCellData]
[$output GetCellData] SetScalars colors
colors Delete; #reference counting - it's ok
randomColorGenerator Delete
}
# This does not need a hierarchical mapper, but hierarchical
# mapper could use a test that has clipping so we use it here
vtkHierarchicalPolyDataMapper sphereMapper
sphereMapper SetInputConnection [randomColors GetOutputPort 0]
vtkActor sphereActor
sphereActor SetMapper sphereMapper
# Set up the planes
vtkPlane plane1
plane1 SetOrigin 0 0 -10
plane1 SetNormal 0 0 1
vtkPlane plane2
plane2 SetOrigin 0 0 10
plane2 SetNormal 0 0 -1
vtkPlane plane3
plane3 SetOrigin -10 0 0
plane3 SetNormal 1 0 0
vtkPlane plane4
plane4 SetOrigin 10 0 0
plane4 SetNormal -1 0 0
sphereMapper AddClippingPlane plane1
sphereMapper AddClippingPlane plane2
volumeMapper AddClippingPlane plane3
volumeMapper AddClippingPlane plane4
# Okay now the graphics stuff
vtkRenderer ren1
vtkRenderWindow renWin
renWin AddRenderer ren1
renWin SetSize 256 256
vtkRenderWindowInteractor iren
iren SetRenderWindow renWin
[ren1 GetCullers] InitTraversal
set culler [[ren1 GetCullers] GetNextItem]
$culler SetSortingStyleToBackToFront
ren1 AddViewProp sphereActor
ren1 SetBackground 0.1 0.2 0.4
renWin Render
[ren1 GetActiveCamera] Azimuth 45
[ren1 GetActiveCamera] Elevation 15
[ren1 GetActiveCamera] Roll 45
[ren1 GetActiveCamera] Zoom 2.0
set valid [volumeMapper IsRenderSupported volumeProperty ren1]
ren1 AddViewProp volume
if {!$valid} {
ren1 RemoveAllViewProps
vtkTextActor t
t SetInput "Required Extensions Not Supported"
t SetDisplayPosition 128 128
[t GetTextProperty] SetJustificationToCentered
ren1 AddViewProp t
}
wm withdraw .
iren Initialize
for { set i 0 } { $i < 5 } { incr i } {
volume RotateY 17
volume RotateZ 13
sphereActor RotateX 13
sphereActor RotateY 17
renWin Render
}
|
