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
|
package require vtk
# Image pipeline
vtkMINCImageReader reader
reader SetFileName "$VTK_DATA_ROOT/Data/t3_grid_0.mnc"
reader RescaleRealValuesOn
vtkMINCImageAttributes attributes
set image reader
#
# write to the temp directory if possible, otherwise use .
#
set dir "."
if {[info commands "rtTester"] == "rtTester"} {
set dir [rtTester GetTempDirectory]
}
# make sure it is writeable first
if {[catch {set channel [open "$dir/test.tmp" "w"]}] == 0 } {
close $channel
file delete -force "$dir/test.tmp"
vtkMINCImageWriter minc1
minc1 SetInputConnection [reader GetOutputPort]
minc1 SetFileName "$dir/minc1.mnc"
attributes ShallowCopy [reader GetImageAttributes]
attributes SetAttributeValueAsString "patient" "full_name" "DOE^JOHN DAVID"
vtkMINCImageWriter minc2
minc2 SetImageAttributes attributes
minc2 SetInputConnection [reader GetOutputPort]
minc2 SetFileName "$dir/minc2.mnc"
vtkMINCImageWriter minc3
minc3 SetImageAttributes attributes
minc3 AddInputConnection [reader GetOutputPort]
minc3 AddInputConnection [reader GetOutputPort]
minc3 SetFileName "$dir/minc3.mnc"
minc1 Write
minc2 Write
minc3 Write
vtkMINCImageReader reader2
reader2 SetFileName "$dir/minc3.mnc"
reader2 RescaleRealValuesOn
reader2 SetTimeStep 1
reader2 Update
set image reader2
file delete -force "$dir/minc1.mnc"
file delete -force "$dir/minc2.mnc"
file delete -force "$dir/minc3.mnc"
# write out the file header for coverage
attributes PrintFileHeader
}
vtkImageViewer viewer
viewer SetInputConnection [$image GetOutputPort]
viewer SetColorWindow 100
viewer SetColorLevel 0
viewer Render
|
