20 if len(sys.argv) != 6:
21 print(
"Usage: " + sys.argv[0] +
" <inputImage> <outputImage> "
22 "<numberOfHistogramBins> <numberOfThresholds> <labelOffset>")
25 inputImage = sys.argv[1]
26 outputImage = sys.argv[2]
27 numberOfHistogramBins = int(sys.argv[3])
28 numberOfThresholds = int(sys.argv[4])
29 labelOffset = int(sys.argv[5])
34 ImageType =
itk.Image[PixelType, Dimension]
37 reader.SetFileName(inputImage)
42 thresholdFilter.SetInput(reader.GetOutput())
44 thresholdFilter.SetNumberOfHistogramBins(numberOfHistogramBins)
45 thresholdFilter.SetNumberOfThresholds(numberOfThresholds)
46 thresholdFilter.SetLabelOffset(labelOffset)
49 rescaler.SetInput(thresholdFilter.GetOutput())
50 rescaler.SetOutputMinimum(0)
51 rescaler.SetOutputMaximum(255)
54 writer.SetFileName(outputImage)
55 writer.SetInput(rescaler.GetOutput())