25 if len(sys.argv) != 4:
26 print(
"Usage: " + sys.argv[0] +
" <inputImage> <outputImage> <radius>")
29 inputImage = sys.argv[1]
30 outputImage = sys.argv[2]
31 radiusValue = int(sys.argv[3])
36 ImageType =
itk.Image[PixelType, Dimension]
39 reader = ReaderType.New()
40 reader.SetFileName(inputImage)
43 structuringElement = StructuringElementType.Ball(radiusValue)
46 ImageType, ImageType, StructuringElementType]
47 grayscaleFilter = GrayscaleFilterType.New()
48 grayscaleFilter.SetInput(reader.GetOutput())
49 grayscaleFilter.SetKernel(structuringElement)
52 writer = WriterType.New()
53 writer.SetFileName(outputImage)
54 writer.SetInput(grayscaleFilter.GetOutput())