20 if len(sys.argv) != 4:
21 print(
"Usage: " + sys.argv[0] +
22 " <InputFileName> <LabelMap> <OutputFileName>")
25 inputFileName = sys.argv[1]
26 labelFileName = sys.argv[2]
27 outputFileName = sys.argv[3]
29 PixelType = itk.ctype(
'unsigned char')
32 ImageType =
itk.Image[PixelType, Dimension]
36 reader.SetFileName(inputFileName)
39 labelReader.SetFileName(labelFileName)
42 LabelType = itk.ctype(
'unsigned long')
47 converter.SetInput(labelReader)
51 overlayFilter.SetInput(converter.GetOutput())
52 overlayFilter.SetFeatureImage(reader.GetOutput())
53 overlayFilter.SetOpacity(0.5)
56 itk.imwrite(overlayFilter.GetOutput(), outputFileName)