MRI-CT phantom registration

I am trying to register (rigid) a MRI image of a phantom to its CT image. However, I am not getting the results I need. I am using MattesMutual Metric, and regularstepgradient, Euler3DTransform, etc. Has anyone encountered similar problem and solved it?

Here is the code:

Thank you,

It looks like repetitive grid pattern quickly puts the registration metric into a local minimum. This could benefit from multi-resolution approach, especially with a lot of initial downsampling.

I just took a quick look at your source code. It looks like you are using some multi-resolution. Try using about 10 levels instead of 3.

Is the initialization good? It looks from your code that you’ve scaled the CT to be positive, which is a good start, but it’s worth doing a visual check of where the optimization is starting from.

1 Like