20 parser = argparse.ArgumentParser(description=
"Resample A Scalar Image.")
21 parser.add_argument(
"input_image")
22 parser.add_argument(
"output_image")
23 parser.add_argument(
"size_x", type=int)
24 parser.add_argument(
"size_y", type=int)
25 args = parser.parse_args()
27 output_size = [args.size_x, args.size_y]
29 input_image = itk.imread(args.input_image)
31 input_spacing = itk.spacing(input_image)
32 input_size = itk.size(input_image)
33 dimension = input_image.GetImageDimension()
35 input_spacing[dim] * input_size[dim] / output_size[dim]
for dim
in range(dimension)
40 output_image = itk.resample_image_filter(
43 output_spacing=output_spacing,
44 output_origin=itk.origin(input_image),
48 itk.imwrite(output_image, args.output_image)