24 if len(sys.argv) != 6:
25 print(
"Usage: " + sys.argv[0] +
" <inputImage> <outputImage> "
26 "<numberOfHistogramBins> <numberOfThresholds> <labelOffset>")
29 inputImage = sys.argv[1]
30 outputImage = sys.argv[2]
31 numberOfHistogramBins = int(sys.argv[3])
32 numberOfThresholds = int(sys.argv[4])
33 labelOffset = int(sys.argv[5])
38 ImageType =
itk.Image[PixelType, Dimension]
41 reader.SetFileName(inputImage)
46 thresholdFilter.SetInput(reader.GetOutput())
48 thresholdFilter.SetNumberOfHistogramBins(numberOfHistogramBins)
49 thresholdFilter.SetNumberOfThresholds(numberOfThresholds)
50 thresholdFilter.SetLabelOffset(labelOffset)
53 rescaler.SetInput( thresholdFilter.GetOutput() )
54 rescaler.SetOutputMinimum( 0 )
55 rescaler.SetOutputMaximum( 255 )
58 writer.SetFileName(outputImage)
59 writer.SetInput(rescaler.GetOutput())