17 python headsq2dcm.py -D /path/to/VTKData
22 from vtk.util.misc
import vtkGetDataRoot
23 VTK_DATA_ROOT = vtkGetDataRoot()
25 reader = vtk.vtkVolume16Reader()
26 reader.SetDataDimensions(64, 64)
27 reader.SetDataByteOrderToLittleEndian()
28 reader.SetFilePrefix(VTK_DATA_ROOT +
"/Data/headsq/quarter")
29 reader.SetImageRange(1, 93)
30 reader.SetDataSpacing(3.2, 3.2, 1.5)
32 cast = vtk.vtkImageCast()
33 cast.SetInput( reader.GetOutput() )
34 cast.SetOutputScalarTypeToUnsignedChar()
37 writer = vtkgdcm.vtkGDCMImageWriter()
38 writer.SetFileName(
"headsq.dcm" )
39 writer.SetInput( reader.GetOutput() )
42 writer.SetFileDimensionality( 3 )