20 if len(sys.argv) != 3:
21 print(
"Usage: " + sys.argv[0] +
" <inputImage> <outputImage> ")
24 inputImage = sys.argv[1]
25 outputImage = sys.argv[2]
27 InputPixelType = itk.F
28 OutputPixelType = itk.UC
31 InputImageType =
itk.Image[InputPixelType, Dimension]
32 OutputImageType =
itk.Image[OutputPixelType, Dimension]
35 reader = ReaderType.New()
36 reader.SetFileName(inputImage)
39 laplacianFilter = FilterType.New()
40 laplacianFilter.SetInput(reader.GetOutput())
43 InputImageType, OutputImageType]
44 rescaler = RescaleFilterType.New()
45 rescaler.SetInput(laplacianFilter.GetOutput())
47 outputPixelTypeMinimum = itk.NumericTraits[OutputPixelType].min()
48 outputPixelTypeMaximum = itk.NumericTraits[OutputPixelType].max()
50 rescaler.SetOutputMinimum(outputPixelTypeMinimum)
51 rescaler.SetOutputMaximum(outputPixelTypeMaximum)
54 writer = WriterType.New()
55 writer.SetFileName(outputImage)
56 writer.SetInput(rescaler.GetOutput())