20 parser = argparse.ArgumentParser(description=
"Compute Laplacian.")
21 parser.add_argument(
"input_image")
22 parser.add_argument(
"output_image")
23 args = parser.parse_args()
25 InputPixelType = itk.F
26 OutputPixelType = itk.UC
29 InputImageType =
itk.Image[InputPixelType, Dimension]
30 OutputImageType =
itk.Image[OutputPixelType, Dimension]
33 reader = ReaderType.New()
34 reader.SetFileName(args.input_image)
37 laplacianFilter = FilterType.New()
38 laplacianFilter.SetInput(reader.GetOutput())
41 rescaler = RescaleFilterType.New()
42 rescaler.SetInput(laplacianFilter.GetOutput())
47 rescaler.SetOutputMinimum(outputPixelTypeMinimum)
48 rescaler.SetOutputMaximum(outputPixelTypeMaximum)
51 writer = WriterType.New()
52 writer.SetFileName(args.output_image)
53 writer.SetInput(rescaler.GetOutput())