[Insight-users] Mattes MI: "Too many samples map outside..."

[Christoph Palm]

Hi Xabi,

you can try to increase the number of spatial samples, but I cannot
believe, this will help you. If the overlap is better than in the
first iteration, this is really a strange thing. Other idea: what about
setting a region and/or a mask? Did you define anything like this or
do you simply use the images as they are? The overlap is not the overlap
of the images but the masked overlap of the regions set (if they are
set). And: regard the spacing, because the registration framework works
on real world measures.

Sorry, but without seeing the images, this seems to be my last idea :-(

Regards

Christoph

On Fri, 2007-02-09 at 12:30 +0100, Xabier Artaechevarria Artieda wrote:
> Hallo Christoph,
> Thanks for the answer.
> I checked after 174 iterations and the overlap seems in fact good, and
> it is considerably better than at the first iteration. Do you think I
> could solve the problem by using a higher number of spatial samples?
> Now I am using one tenth of the total voxels.
> Regards,
> Xabi
> 
> 
> 
> Christoph Palm <christoph.palm at web.de> ha escrito:
> 
> > Hi Xabi,
> >
> > do they overlap well also after the 174 iterations?
> > Perhaps you try to apply the last transformation paramters onto your
> > moving image and look what happens.
> >
> > Regards
> >
> > Christoph
> >
> > On Thu, 2007-02-08 at 11:30 +0100, Xabier Artaechevarria Artieda wrote:
> >> Hi all,
> >> I am registering multiple images using Bspline transform. The program
> >> is working fine with all images except with one couple. I get the
> >> following error after 174 iterations:
> >>
> >> 170   -1.07817   4.07948e-05
> >> 171   -1.07818   4.09635e-05
> >> 172   -1.07818   5.61295e-05
> >> 173   -1.07818   4.65465e-05
> >> 174   -1.07818   4.64441e-05
> >> ExceptionObject caught !
> >>
> >> itk::ExceptionObject (0x6c81a0)
> >> Location: "void
> >> itk::MattesMutualInformationImageToImageMetric<TFixedImage,
> >> TMovingImage>::GetValueAndDerivative(const typename
> >> itk::ImageToImageMetric<TFixedImage, TMovingImage>::ParametersType&,
> >> typename itk::ImageToImageMetric<TFixedImage,
> >> TMovingImage>::MeasureType&, typename
> >> itk::ImageToImageMetric<TFixedImage, TMovingImage>::DerivativeType&)
> >> const [with TFixedImage = main(int, char**)::FixedImageType,
> >> TMovingImage = main(int, char**)::MovingImageType]"
> >> File:
> >> /usr/local/etc/ITK/Code/Algorithms/itkMattesMutualInformationImageToImageMetric.txx
> >> Line: 995
> >> Description: itk::ERROR:
> >> MattesMutualInformationImageToImageMetric(0x6ab920): Too many samples
> >> map outside moving image buffer: 0 / 184251
> >>
> >> I have seen a previous thread where it was mentioned that this might
> >> be due to a lack of overlap between the images. However, at least
> >> visually, images seem to overlap quite well.
> >> Is there any other possible reason? And ways to solve it?
> >> Thanks in advance,
> >> Xabi
> >>