20 parser = argparse.ArgumentParser(description=
"Wrap An Image Using A Deformation Field.")
21 parser.add_argument(
"input_image")
22 parser.add_argument(
"displacement_field")
23 parser.add_argument(
"output_image")
24 args = parser.parse_args()
28 VectorComponentType = itk.F
29 VectorPixelType =
itk.Vector[VectorComponentType, Dimension]
31 DisplacementFieldType =
itk.Image[VectorPixelType, Dimension]
34 ImageType =
itk.Image[PixelType, Dimension]
37 reader.SetFileName(args.input_image)
40 fieldReader.SetFileName(args.displacement_field)
43 deformationField = fieldReader.GetOutput()
49 warpFilter.SetInterpolator(interpolator)
51 warpFilter.SetOutputSpacing(deformationField.GetSpacing())
52 warpFilter.SetOutputOrigin(deformationField.GetOrigin())
53 warpFilter.SetOutputDirection(deformationField.GetDirection())
55 warpFilter.SetDisplacementField(deformationField)
57 warpFilter.SetInput(reader.GetOutput())
60 writer.SetInput(warpFilter.GetOutput())
61 writer.SetFileName(args.output_image)