1 """=========================================================================
5 ' Licensed under the Apache License, Version 2.0 (the "License");
6 ' you may not use this file except in compliance with the License.
7 ' You may obtain a copy of the License at
9 ' http://www.apache.org/licenses/LICENSE-2.0.txt
11 ' Unless required by applicable law or agreed to in writing, software
12 ' distributed under the License is distributed on an "AS IS" BASIS,
13 ' WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14 ' See the License for the specific language governing permissions and
15 ' limitations under the License.
17 '========================================================================="""
19 import SimpleITK
as sitk
28 "Usage: ConnectedThresholdImageFilter inputImage outputImage",
29 "lowerThreshold upperThreshold seedX seedY [seed2X seed2Y ... ]",
37 reader.SetFileName(sys.argv[1])
38 image = reader.Execute()
44 blurFilter.SetNumberOfIterations(5)
45 blurFilter.SetTimeStep(0.125)
46 image = blurFilter.Execute(image)
52 segmentationFilter.SetLower(float(sys.argv[3]))
53 segmentationFilter.SetUpper(float(sys.argv[4]))
54 segmentationFilter.SetReplaceValue(255)
56 for i
in range(5, len(sys.argv) - 1, 2):
57 seed = [int(sys.argv[i]), int(sys.argv[i + 1])]
58 segmentationFilter.AddSeed(seed)
59 print(
"Adding seed at: ", seed,
" with intensity: ", image.GetPixel(*seed))
62 image = segmentationFilter.Execute(image)
69 writer.SetFileName(sys.argv[2])
72 if "SITK_NOSHOW" not in os.environ: