[Python] B-Spline deformable registration for 2 CT images

I have two Dicom CT images and I want to perform B-Spline deformable registration. But I'm not very familiar with registration.

The below code is part of the code related to SimpleITK Examples. It seems that this code is written for images with different modalities. But I want to register two CT images together. can you help me how can I modify it for mono-modality registration?

transformDomainMeshSize = [10] * moving.GetDimension()
tx = sitk.BSplineTransformInitializer(fixed, transformDomainMeshSize)
R = sitk.ImageRegistrationMethod()
R.SetMetricAsMattesMutualInformation(50)
R.SetOptimizerAsGradientDescentLineSearch(
5.0, 100, convergenceMinimumValue=1e-4, convergenceWindowSize=5)
R.SetOptimizerScalesFromPhysicalShift()
R.SetInitialTransform(tx)

R.SetInterpolator(sitk.sitkLinear)

R.SetShrinkFactorsPerLevel([6, 2, 1])
R.SetSmoothingSigmasPerLevel([6, 2, 1])


The similarity metrics and optimizers I was able to find are the following, but I don't know which one I can use.

R.SetMetricAsANTSNeighborhoodCorrelation
R.SetMetricAsCorrelation
R.SetMetricAsDemons
R.SetMetricAsJointHistogramMutualInformation
R.SetMetricAsMattesMutualInformation
R.SetMetricAsMeanSquares


R.SetOptimizerAsAmoeba
R.SetOptimizerAsConjugateGradientLineSearch
R.SetOptimizerAsExhaustive
R.SetOptimizerAsGradientDescent
R.SetOptimizerAsGradientDescentLineSearch
R.SetOptimizerAsLBFGS2
R.SetOptimizerAsLBFGSB
R.SetOptimizerAsOnePlusOneEvolutionary
R.SetOptimizerAsPowell
R.SetOptimizerAsRegularStepGradientDescent


Thanks a lot

Continue reading...