[Insight-users] register of MRI and a partial fMRI

[Karthik Krishnan]

Hi Yu,

Print out the parameters of the transform. There just isn't enough
overlap between the fixed and the moving image to start with (or the
transform is driving you too far away too soon). See
itkMattesMutualInformationImageToImageMetric.txx


Line: 702
if( nSamples < m_NumberOfSpatialSamples / 4 )
{
itkExceptionMacro( "Too many samples map outside moving image buffer: "
<< nSamples << " / " << m_NumberOfSpatialSamples << std::endl );
}


Remember that if you are using the centered transform, set the center
(for example with the CenteredTransformInitializer helper class).

See Examples/Registration/ImageRegistration9.cxx

Thanks
kk

Yu Qi wrote:

> Hi Luis,
> Thanks for your instruction.
> I tried to set region to the fixed image by setting
> start[2]=size[2]/2 and size[2]/=2 immediately after reading the image.
> I use the Mattes MI metric and Centered Rigid 3D image transformation.
> But every time I tried to run the .exe, it always told me there are
> too many samples outside the moving image. And I choose a smaller
> number for number of samples. I still got the same exception.
> Why should this happen?
> It works well before I set the region.
> Best Regards,
> Yu Qi
>
> */Luis Ibanez <luis.ibanez at kitware.com>/* wrote:
>
>
>     Hi Yu,
>
>     Yes, you can use ITK for registering a fMRI to an MRI
>     even if they dont fully cover the same anatomical
>     regions.
>
>     You may want to use the MRI image as the Fixed image and
>     set the FixedImageRegion to the rectangular box that covers
>     the brain plus a border of about 10 millimeters. Take the
>     fMRI image as the Moving image.
>
>     You may want to use any of the MutualInformation Metrics
>     available in ITK.
>
>     For a full list of the Image Metrics please look at the
>     Groups in Doxygen:
>
>     http://www.itk.org/Insight/Doxygen/html/group__RegistrationMetrics.html
>
>
>
>     Regards,
>
>
>     Luis
>
>
>     ------------------
>     Yu Qi wrote:
>
>     > Hi,
>     >
>     > I have a MRI volume of the whole head and
>     > a fMRI volume of the brain.
>     >
>     > Can I register these two directly with ITK?
>     >
>     > Do I! need to use the SetRegion() to MRI volume?
>     >
>     > Thanks!
>     >
>     >
>     >
>     >
>     > *Best Regards!*
>     > **
>     > *Yu Qi*
>     >
>     >
>     >
>     >
>     >
>     ------------------------------------------------------------------------
>     > *Do You Yahoo!?*
>     > 150万曲MP3疯狂搜，带您闯入音乐殿堂
>     > 美女明星应有尽有，搜遍美图、艳图和酷图
>     > 1G就是1000兆，雅虎电邮自助扩容！
>     >
>     >
>     >
>     >
>     >
>     ------------------------------------------------------------------------
>     >
>     > _______________________________________________
>     > Insight-users mailing list
>     > Insight-users at itk.org
>     > http://www.itk.org/mailman/listinfo/insight-users
>
>
>
>
> ------------------------------------------------------------------------
> *Do You Yahoo!?*
> 150万曲MP3疯狂搜，带您闯入音乐殿堂 <http://music.yisou.com>
> 美女明星应有尽有，搜遍美图、艳图和酷图 <http://image.yisou.com>
> 1G就是1000兆，雅虎电邮自助扩容！
> <http://cn.rd.yahoo.com/mail_cn/tag/1g/*http://cn.mail.yahoo.com/event/mail_1g/>
>
>
>------------------------------------------------------------------------
>
>_______________________________________________
>Insight-users mailing list
>Insight-users at itk.org
>http://www.itk.org/mailman/listinfo/insight-users
>  
>