1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20""" An example script showing the use of the ConnectedThresholdImageFilter """
21
22import sys
23import os
24import SimpleITK as sitk
25
26
27
28
29if len(sys.argv) < 7:
30 print(
31 "Usage: ConnectedThresholdImageFilter inputImage outputImage",
32 "lowerThreshold upperThreshold seedX seedY [seed2X seed2Y ... ]",
33 )
34 sys.exit(1)
35
36
37
38
40reader.SetFileName(sys.argv[1])
41image = reader.Execute()
42
43
44
45
47blurFilter.SetNumberOfIterations(5)
48blurFilter.SetTimeStep(0.125)
49image = blurFilter.Execute(image)
50
51
52
53
55segmentationFilter.SetLower(float(sys.argv[3]))
56segmentationFilter.SetUpper(float(sys.argv[4]))
57segmentationFilter.SetReplaceValue(255)
58
59for i in range(5, len(sys.argv) - 1, 2):
60 seed = [int(sys.argv[i]), int(sys.argv[i + 1])]
61 segmentationFilter.AddSeed(seed)
62 print("Adding seed at: ", seed, " with intensity: ", image.GetPixel(*seed))
63
64
65image = segmentationFilter.Execute(image)
66image[seed] = 255
67
68
69
70
72writer.SetFileName(sys.argv[2])
73writer.Execute(image)
74
75if "SITK_NOSHOW" not in os.environ:
Label pixels that are connected to a seed and lie within a range of values.
Denoise an image using curvature driven flow.
Read an image file and return a SimpleITK Image.
Write out a SimpleITK image to the specified file location.
void SITKIO_EXPORT Show(const Image &image, const std::string &title="", const bool debugOn=ProcessObject::GetGlobalDefaultDebug())