[Insight-users] 3 image registration

[Luis Ibanez]

Hi Niels,

I agree with Torsten in that probably the best way to go is for you
to create a new ImageToImageMetric that internally contains two
instances of the MutualInformationImageToImageMetric, and computes
the GetValue() as the sum of the GetValues returned by each one of
these internal instances.

Your intuition is probably correct in that a 3D Joint Histogram is
the theoretically correct way of implementing this problem, but that
will require a major programming effort on your part.

If you decide to go that way, you probably will find easier to manually
create the Joint histogram by following the examples in


         Insight/Examples/Statistics/
                     ImageMutualInformation1.cxx
                     ImageToListAdaptor.cxx



Keep in mind that Torsten's suggestion will be much easier to implement.



    Regards,



       Luis



----------------------
Niels Zelders wrote:
> Hello itk-users,
> 
> I would really appreciate some help in solving the following problem:
> I'm trying to solve a multiple image registration problem and wondered 
> if the ITK already contained some components for this.
> To be more precise, I am trying to register a (fixed) 3D ultrasound 
> image to a couple of (moving) 3D MR images. The MR images have the same 
> domain, so it's not a general 3 image registration problem, but more 
> like a scalar-image-to-double-valued -'vector'-image-problem. I do not 
> know how to adapt the 'itkMattesMutualInformationMetric' (I prefer to 
> use) in order to accept two moving images or a moving vector image 
> (probably due to the fact that the neccesary 3D histogram construction 
> is rather tricky to implement in a meaningful way??)
> Does anyone know how to achieve this or have some hints on how to proceed?
> 
> Thanks in advance
> 
> Niels
> 
> 
> 
>