Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import itk, sys
- if len(sys.argv) != 7:
- print("Usage: " + sys.argv[0] + " <inputImagePath> <maskImagePath>"
- " <NumberOfBinsPerAxis> <PixelValueMin> "
- "<PixelValueMax> <NeighborhoodRadius>")
- sys.exit(1)
- im = itk.imread(sys.argv[1])
- maskReader = itk.ImageFileReader.IUC3.New()
- maskReader.SetFileName(sys.argv[2])
- mask = maskReader.GetOutput()
- filtr = itk.CoocurrenceTextureFeaturesImageFilter.New(im)
- filtr.SetMaskImage(mask)
- filtr.SetNumberOfBinsPerAxis(int(sys.argv[3]))
- filtr.SetHistogramMinimum(int(sys.argv[4]))
- filtr.SetHistogramMaximum(int(sys.argv[5]))
- filtr.SetNeighborhoodRadius([int(sys.argv[6]),int(sys.argv[6]),int(sys.argv[6])])
- result = filtr.GetOutput()
- itk.imwrite(result, "result.nrrd")
Add Comment
Please, Sign In to add comment