1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20""" A SimpleITK example demonstrating the classic Demons image registration. """
21
22import sys
23import os
24import SimpleITK as sitk
25
26
27def command_iteration(sitk_filter):
28 """ Callback invoked when the filter processes an iteration. """
29 print(f"{sitk_filter.GetElapsedIterations():3} = {sitk_filter.GetMetric():10.5f}")
30
31
32if len(sys.argv) < 4:
33 print(
34 f"Usage: {sys.argv[0]}"
35 + " <fixedImageFilter> <movingImageFile> <outputTransformFile>"
36 )
37 sys.exit(1)
38
40
42
44matcher.SetNumberOfHistogramLevels(1024)
45matcher.SetNumberOfMatchPoints(7)
46matcher.ThresholdAtMeanIntensityOn()
47moving = matcher.Execute(moving, fixed)
48
49
50
51
53demons.SetNumberOfIterations(50)
54
55demons.SetStandardDeviations(1.0)
56
57demons.AddCommand(sitk.sitkIterationEvent, lambda: command_iteration(demons))
58
59displacementField = demons.Execute(fixed, moving)
60
61print("-------")
62print(f"Number Of Iterations: {demons.GetElapsedIterations()}")
63print(f" RMS: {demons.GetRMSChange()}")
64
66
68
69if "SITK_NOSHOW" not in os.environ:
71 resampler.SetReferenceImage(fixed)
72 resampler.SetInterpolator(sitk.sitkLinear)
73 resampler.SetDefaultPixelValue(100)
74 resampler.SetTransform(outTx)
75
76 out = resampler.Execute(moving)
79
80
81
82 cimg =
sitk.Compose(simg1, simg2, simg1 // 2.0 + simg2 // 2.0)
83 sitk.Show(cimg,
"DeformableRegistration1 Composition")
Deformably register two images using the demons algorithm.
Normalize the grayscale values for a source image by matching the shape of the source image histogram...
Resample an image via a coordinate transform.
SITKIO_EXPORT Image ReadImage(const PathType &filename, PixelIDValueEnum outputPixelType=sitkUnknown, const std::string &imageIO="")
ReadImage is a procedural interface to the ImageFileReader class which is convenient for most image r...
Image Compose(const std::vector< Image > &images)
ComposeImageFilter combine several scalar images into a multicomponent image.
void SITKIO_EXPORT Show(const Image &image, const std::string &title="", const bool debugOn=ProcessObject::GetGlobalDefaultDebug())
SITKCommon_EXPORT void WriteTransform(const Transform &transform, const PathType &filename)
Image RescaleIntensity(Image &&image1, double outputMinimum=0, double outputMaximum=255)
Applies a linear transformation to the intensity levels of the input Image .
Image Cast(const Image &image, PixelIDValueEnum pixelID)