20 from distutils.version
import StrictVersion
as VS
23 print(
"ITK 5.0.0 or newer is required.")
26 parser = argparse.ArgumentParser(description=
"Segment blood vessels.")
27 parser.add_argument(
"input_image")
28 parser.add_argument(
"output_image")
29 parser.add_argument(
"--sigma", type=float, default=1.0)
30 parser.add_argument(
"--alpha1", type=float, default=0.5)
31 parser.add_argument(
"--alpha2", type=float, default=2.0)
32 args = parser.parse_args()
34 input_image = itk.imread(args.input_image, itk.ctype(
"float"))
36 hessian_image = itk.hessian_recursive_gaussian_image_filter(
37 input_image, sigma=args.sigma
43 vesselness_filter.SetInput(hessian_image)
44 vesselness_filter.SetAlpha1(args.alpha1)
45 vesselness_filter.SetAlpha2(args.alpha2)
47 itk.imwrite(vesselness_filter, args.output_image)