<DIV>Dear Luis,</DIV>
<DIV> </DIV>
<DIV>Thank you for your reply. Sorry that I was not able to feed back promptly. I didn't express the question clearly in my previous email.</DIV>
<DIV> </DIV>
<DIV>During the registration process, the metric will be calculated based on the overlap of the fixed and moving image for each iteration. Since the area of the overlap may change for each iteration, the number of total pixels counted as the overlap may change, so the value of metric may have discontinuity that might cause the failure in finding the global extreme.</DIV>
<DIV> </DIV>
<DIV>To avoid above, different weights could be given to pixels or voxels within the overlap. For example, when calculating the metric, pixels near the boundary are given lower weights, but those ones near the center are given higher weights. I don't know if ITK has done this kind of stuff. Thanks. </DIV>
<DIV> </DIV>
<DIV>Regards,</DIV>
<DIV> </DIV>
<DIV>Jay</DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV><BR><B><I>Luis Ibanez <luis.ibanez@kitware.com></I></B> wrote:</DIV>
<BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #1010ff 2px solid"><BR>Hi Jay,<BR><BR>By apodization do you refer to the attenuation of<BR>the high frequencies in the spectral domain of the<BR>images ?<BR><BR>I assume that you are concerned about this due to the<BR>resampling of the images during the registration process<BR>as well as after the transform has been extracted.<BR><BR>If that's the case, you don't really need to be too<BR>concerned. It practice this is rarely a problem, unless<BR>your images are actually populated with white noise,<BR>or at least blue noise, which are the type of spectra<BR>having large components in the high frequency range and<BR>therefore requiring attenuation in order to avoid aliasing<BR>during the interpolation used for resampling.<BR><BR>For typical medical images, the high frequency components<BR>are very low, and you are usually safe by using simple<BR>interpolations. That is, the Fourier transform of
the<BR>shape functions used by typical interpolation are already<BR>attenuatting high frequencies.<BR><BR>You may want to look at the documentation of the<BR>SincWindowed interpolator where some of this issues<BR>are discussed.<BR><BR><BR>From the practical point of view your best options<BR>for avoiding local minima are :<BR><BR>- Multi-resolution<BR>- Good initialization<BR><BR><BR><BR>Please let us know if you have further questions.<BR><BR><BR><BR>Regards,<BR><BR><BR><BR>Luis<BR><BR><BR><BR>-----------------<BR>Jay Li wrote:<BR><BR>> Dear all,<BR>> <BR>> I'm new to ITK. Could anyone give me a hint if ITK has<BR>> taken care of the apodization related problems in the<BR>> process of image registration. I looked through the<BR>> manual and searched the online archives, but couldn't<BR>> find anything related. Will there be any problems to<BR>> find the global extreme without the apodization?<BR>> Thanks in advance.<BR>> <BR>> Jay<BR>> <BR>>
<BR>> <BR><BR><BR><BR></BLOCKQUOTE><p>
<hr size=1>Do you Yahoo!?<br>
<a href="http://us.rd.yahoo.com/mail_us/taglines/aac/*http://promotions.yahoo.com/new_mail/static/ease.html">Yahoo! Mail Address AutoComplete</a> - You start. We finish.