20 if len(sys.argv) != 4:
21 print(
'Usage: ' + sys.argv[0] +
22 ' <InputFileName> <DisplacementFieldFileName> <OutputFileName>')
25 inputFileName = sys.argv[1]
26 displacementFieldFileName = sys.argv[2]
27 outputFileName = sys.argv[3]
31 VectorComponentType = itk.F
32 VectorPixelType =
itk.Vector[VectorComponentType, Dimension]
34 DisplacementFieldType =
itk.Image[VectorPixelType, Dimension]
37 ImageType =
itk.Image[PixelType, Dimension]
40 reader.SetFileName(inputFileName)
43 fieldReader.SetFileName(displacementFieldFileName)
46 deformationField = fieldReader.GetOutput()
53 warpFilter.SetInterpolator(interpolator)
55 warpFilter.SetOutputSpacing(deformationField.GetSpacing())
56 warpFilter.SetOutputOrigin(deformationField.GetOrigin())
57 warpFilter.SetOutputDirection(deformationField.GetDirection())
59 warpFilter.SetDisplacementField(deformationField)
61 warpFilter.SetInput(reader.GetOutput())
64 writer.SetInput(warpFilter.GetOutput())
65 writer.SetFileName(outputFileName)