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
|
#!/usr/bin/env python
import os
import vtk
from vtk.test import Testing
from vtk.util.misc import vtkGetDataRoot
VTK_DATA_ROOT = vtkGetDataRoot()
# Image pipeline
reader = vtk.vtkMINCImageReader()
reader.SetFileName(VTK_DATA_ROOT + "/Data/t3_grid_0.mnc")
reader.RescaleRealValuesOn()
attributes = vtk.vtkMINCImageAttributes()
image = reader
# The current directory must be writeable.
#
try:
channel = open("minc1.mnc", "wb")
channel.close()
minc1 = vtk.vtkMINCImageWriter()
minc1.SetInputConnection(reader.GetOutputPort())
minc1.SetFileName("minc1.mnc")
attributes.ShallowCopy(reader.GetImageAttributes())
attributes.SetAttributeValueAsString(
"patient", "full_name", "DOE^JOHN DAVID")
minc2 = vtk.vtkMINCImageWriter()
minc2.SetImageAttributes(attributes)
minc2.SetInputConnection(reader.GetOutputPort())
minc2.SetFileName("minc2.mnc")
minc3 = vtk.vtkMINCImageWriter()
minc3.SetImageAttributes(attributes)
minc3.AddInputConnection(reader.GetOutputPort())
minc3.AddInputConnection(reader.GetOutputPort())
minc3.SetFileName("minc3.mnc")
minc1.Write()
minc2.Write()
minc3.Write()
reader2 = vtk.vtkMINCImageReader()
reader2.SetFileName("minc3.mnc")
reader2.RescaleRealValuesOn()
reader2.SetTimeStep(1)
reader2.Update()
image = reader2
# cleanup
#
try:
os.remove("minc1.mnc")
except OSError:
pass
try:
os.remove("minc2.mnc")
except OSError:
pass
try:
os.remove("minc3.mnc")
except OSError:
pass
# Write out the file header for coverage
attributes.PrintFileHeader()
viewer = vtk.vtkImageViewer()
viewer.SetInputConnection(image.GetOutputPort())
viewer.SetColorWindow(100)
viewer.SetColorLevel(0)
viewer.Render()
except IOError:
print "Unable to test the writer/reader."
|
