1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23from InsightToolkit import *
24
25from sys import argv
26
27
28readerInput = itkImageFileReaderUC2_New()
29readerPrior = itkImageFileReaderUC2_New()
30
31readerInput.SetFileName(argv[1])
32readerPrior.SetFileName(argv[2])
33
34readerInput.Update()
35readerPrior.Update()
36
37filter = itkVoronoiSegmentationImageFilterUC2UC2UC2_New()
38
39filter.SetInput(readerInput.GetOutput())
40filter.TakeAPrior(readerPrior.GetOutput())
41
42filter.SetMeanPercentError(eval(argv[4]))
43filter.SetSTDPercentError(eval(argv[5]))
44
45writer = itkImageFileWriterUC2_New()
46writer.SetFileName(argv[3])
47writer.SetInput(filter.GetOutput())
48
49writer.Update()