Ultrasound and MRI Rigid registration (RESECT Dataset)

102 views
Skip to first unread message

David Sixela

unread,
Oct 12, 2021, 4:59:24 AM10/12/21
to elastix-imageregistration
Hello everyone, 

For some experiments, I would like to use Elastix to register Ultrasound and MRI images ( RESECT Dataset) using Advanced Mutual information and adaptive stochastic gradient descent.

So far, the target registration error isn't changing or is getting worse (larger). 
Knowing how difficult registering ultrasound and MRI is, does anyone have a suggestion of parameters I could use to get a first "good enough" registration? 

Thanks in advance.

Please find below the parameter file I use for an Euler transform:

(FixedInternalImagePixelType "float")
(MovingInternalImagePixelType "float")
(UseDirectionCosines "true")

// **************** Main Components **************************

(Registration "MultiResolutionRegistration")
(Interpolator "BSplineInterpolator")
(ResampleInterpolator "FinalBSplineInterpolator")
(Resampler "DefaultResampler")

(FixedImagePyramid "FixedRecursiveImagePyramid")
(MovingImagePyramid "MovingRecursiveImagePyramid")

(Optimizer "AdaptiveStochasticGradientDescent")
(Transform "EulerTransform")
(Metric "AdvancedMattesMutualInformation")

// ***************** Transformation **************************

(AutomaticScalesEstimation "true")
(AutomaticTransformInitialization "false")
(HowToCombineTransforms "Compose")

// ******************* Similarity measure *********************

(NumberOfHistogramBins 32)
(ErodeMask "false")

// ******************** Multiresolution **********************

(NumberOfResolutions 2)

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

(MaximumNumberOfIterations 50)
(MaximumStepLength 0.1)

// **************** Image sampling **********************
(NumberOfSpatialSamples 2048)
(NewSamplesEveryIteration "true")
(ImageSampler "Random")

// ************* Interpolation and Resampling ****************

(BSplineInterpolationOrder 1)
(FinalBSplineInterpolationOrder 3)
(DefaultPixelValue 0)
(WriteResultImage "true")
(ResultImagePixelType "float")
(ResultImageFormat "nii.gz")
(CompressResultImage "true")

m.st...@lumc.nl

unread,
Oct 22, 2021, 4:10:35 AM10/22/21
to elastix-imag...@googlegroups.com

Dear David,


Did you check the model zoo for something that may be a good starting point?

https://elastix.lumc.nl/modelzoo/


Best, Marius



From: elastix-imag...@googlegroups.com <elastix-imag...@googlegroups.com> on behalf of David Sixela <sixel...@gmail.com>
Sent: 12 October 2021 10:59 AM
To: elastix-imageregistration
Subject: [elastix-imageregistration] Ultrasound and MRI Rigid registration (RESECT Dataset)
 
--
You received this message because you are subscribed to the Google Groups "elastix-imageregistration" group.
To unsubscribe from this group and stop receiving emails from it, send an email to elastix-imageregis...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/elastix-imageregistration/f91195d5-e972-4777-9926-09082bf1e931n%40googlegroups.com.
Reply all
Reply to author
Forward
0 new messages