//ImageTypes

(FixedInternalImagePixelType "float")

(FixedImageDimension 3)

(MovingInternalImagePixelType "float")

(MovingImageDimension 3)

//Components

(Registration "MultiResolutionRegistration")

(FixedImagePyramid "FixedSmoothingImagePyramid")

(MovingImagePyramid "MovingSmoothingImagePyramid")

(Interpolator "BSplineInterpolator")

(Metric "AdvancedMattesMutualInformation")

(Optimizer "StandardGradientDescent")

(ResampleInterpolator "FinalBSplineInterpolator")

(Resampler "DefaultResampler")

(Transform "AffineTransform")

// ********** Pyramid

// Total number of resolutions

(NumberOfResolutions 4)

// ********** Transform

(AutomaticTransformInitialization "true")

(AutomaticScalesEstimation "true")

(HowToCombineTransforms "Compose")

// ********** Optimizer

// Maximum number of iterations in each resolution level:

(MaximumNumberOfIterations 1000)

//SP: Param_a in each resolution level. a_k = a/(A+k+1)^alpha

(SP_a 500.0)

//SP: Param_alpha in each resolution level. a_k = a/(A+k+1)^alpha

(SP_alpha 0.602)

//SP: Param_A in each resolution level. a_k = a/(A+k+1)^alpha

(SP_A 50.0)

// ********** Metric

//Number of grey level bins in each resolution level:

(NumberOfHistogramBins 32)

(FixedLimitRangeRatio 0.0)

(MovingLimitRangeRatio 0.0)

(FixedKernelBSplineOrder 1)

(MovingKernelBSplineOrder 3)

// ********** Several

(WriteTransformParametersEachIteration "false")

(WriteTransformParametersEachResolution "false")

(WriteResultImage "false")

(ShowExactMetricValue "false")

(ErodeFixedMask "false")

(ErodeMovingMask "false")

(UseDifferentiableOverlap "false")

// ********** ImageSampler

//Number of spatial samples used to compute the mutual information in each resolution level:

(ImageSampler "RandomSparseMask")

(NumberOfSpatialSamples 2048)

(NewSamplesEveryIteration "true")

// ********** Interpolator and Resampler

//Order of B-Spline interpolation used in each resolution level:

(BSplineInterpolationOrder 1)

//Order of B-Spline interpolation used for applying the final deformation:

(FinalBSplineInterpolationOrder 3)

//Default pixel value for pixels that come from outside the picture:

(DefaultPixelValue 0)