20 if len(sys.argv) != 4:
21 print(
"Usage: " + sys.argv[0] +
" <input_image> <output_image> <scale>")
24 input_file_name = sys.argv[1]
25 output_file_name = sys.argv[2]
26 scale = float(sys.argv[3])
28 input_image = itk.imread(input_file_name)
29 input_size = itk.size(input_image)
30 input_spacing = itk.spacing(input_image)
31 input_origin = itk.origin(input_image)
32 Dimension = input_image.GetImageDimension()
48 output_size = input_size
49 output_spacing = input_spacing
50 output_origin = input_origin
52 scale_transform_parameters = scale_transform.GetParameters()
53 for i
in range(len(scale_transform_parameters)):
54 scale_transform_parameters[i] = scale
55 scale_transform_center = [float(int(s / 2))
for s
in input_size]
56 scale_transform.SetParameters(scale_transform_parameters)
57 scale_transform.SetCenter(scale_transform_center)
61 resampled = itk.resample_image_filter(
63 transform=scale_transform,
64 interpolator=interpolator,
66 output_spacing=output_spacing,
67 output_origin=output_origin,
70 itk.imwrite(resampled, output_file_name)