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
|
package require vtk
# This script is for testing the 3D flood fill filter.
# Image pipeline
vtkRenderWindow imgWin
imgWin SetSize 192 256
vtkImageReader reader
reader ReleaseDataFlagOff
reader SetDataByteOrderToLittleEndian
reader SetDataExtent 0 63 0 63 2 5
reader SetDataSpacing 3.2 3.2 1.5
reader SetFilePrefix "$VTK_DATA_ROOT/Data/headsq/quarter"
reader SetDataMask 0x7fff
vtkPoints seeds
seeds InsertNextPoint 0 0 0
seeds InsertNextPoint 100.8 100.8 0
set replacein "ReplaceInOn ReplaceInOff"
set replaceout "ReplaceOutOn ReplaceOutOff"
set thresholds [list "ThresholdByLower 800" "ThresholdByUpper 1200" "ThresholdBetween 800 1200"]
set k 1
foreach rin $replacein {
foreach rout $replaceout {
foreach t $thresholds {
vtkImageThresholdConnectivity thresh$k
thresh$k SetSeedPoints seeds
thresh$k SetInValue 2000
thresh$k SetOutValue 0
thresh$k $rin
thresh$k $rout
thresh$k SetNeighborhoodRadius 2 2 1
thresh$k SetNeighborhoodFraction 0.5
thresh$k SetInputConnection [reader GetOutputPort]
eval thresh$k $t
vtkImageMapper map$k
map$k SetInputConnection [thresh$k GetOutputPort]
if {$k < 3} {
map$k SetColorWindow 255
map$k SetColorLevel 127.5
} else {
map$k SetColorWindow 2000
map$k SetColorLevel 1000
}
vtkActor2D act$k
act$k SetMapper map$k
vtkRenderer ren$k
ren$k AddActor2D act$k
imgWin AddRenderer ren$k
incr k
}
}
}
ren1 SetViewport 0 0 .33333 .25
ren2 SetViewport .33333 0 .66667 .25
ren3 SetViewport .66667 0 1 .25
ren4 SetViewport 0 .25 .33333 .5
ren5 SetViewport .33333 .25 .66667 .5
ren6 SetViewport .66667 .25 1 .5
ren7 SetViewport 0 .5 .33333 .75
ren8 SetViewport .33333 .5 .66667 .75
ren9 SetViewport .66667 .5 1 .75
ren10 SetViewport 0 .75 .33333 1
ren11 SetViewport .33333 .75 .66667 1
ren12 SetViewport .66667 .75 1 1
imgWin Render
|
