20 import SimpleITK
as sitk
25 def command_iteration(method):
27 f
"{method.GetOptimizerIteration():3} "
28 + f
"= {method.GetMetricValue():10.5f} "
29 + f
": {method.GetOptimizerPosition()}"
37 "ImageRegistrationMethod1",
38 "<fixedImageFilter> <movingImageFile>",
39 "<outputTransformFile>",
48 R.SetMetricAsMeanSquares()
49 R.SetOptimizerAsRegularStepGradientDescent(4.0, 0.01, 200)
51 R.SetInterpolator(sitk.sitkLinear)
53 R.AddCommand(sitk.sitkIterationEvent,
lambda: command_iteration(R))
55 outTx = R.Execute(fixed, moving)
59 print(f
"Optimizer stop condition: {R.GetOptimizerStopConditionDescription()}")
60 print(f
" Iteration: {R.GetOptimizerIteration()}")
61 print(f
" Metric value: {R.GetMetricValue()}")
66 resampler.SetReferenceImage(fixed)
67 resampler.SetInterpolator(sitk.sitkLinear)
68 resampler.SetDefaultPixelValue(100)
69 resampler.SetTransform(outTx)
71 out = resampler.Execute(moving)
74 cimg =
sitk.Compose(simg1, simg2, simg1 // 2.0 + simg2 // 2.0)
76 return_images = {
"fixed": fixed,
82 if __name__ ==
"__main__":
83 return_dict = main(sys.argv)
84 if "SITK_NOSHOW" not in os.environ:
85 sitk.Show(return_dict[
"composition"],
"ImageRegistration1 Composition")