<br><br><div class="gmail_quote">2010/6/17 Luis Ibanez <span dir="ltr"><<a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.com</a>></span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div class="gmail_quote">On Wed, Jun 16, 2010 at 11:36 AM, Simon Warfield <span dir="ltr"><<a href="mailto:simon.warfield@childrens.harvard.edu" target="_blank">simon.warfield@childrens.harvard.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<a href="mailto:insight-users-request@itk.org" target="_blank">insight-users-request@itk.org</a> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Date: Wed, 16 Jun 2010 09:46:09 -0400<br>
From: Luis Ibanez <<a href="mailto:luis.ibanez@kitware.com" target="_blank">luis.ibanez@kitware.com</a>><br>
Subject: Re: [Insight-users] Interpolation and centered pixel<br>
coordinates<br>
To: D?enan Zuki? <<a href="mailto:dzenanz@gmail.com" target="_blank">dzenanz@gmail.com</a>><br>
Cc: <a href="mailto:insight-users@itk.org" target="_blank">insight-users@itk.org</a><br>
Message-ID:<br>
<<a href="mailto:AANLkTikQna2jtbdcwpMTZtWmGC5KZnEMf1dTN2_9fOY-@mail.gmail.com" target="_blank">AANLkTikQna2jtbdcwpMTZtWmGC5KZnEMf1dTN2_9fOY-@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
D?enan,<br>
<br>
We did quite a heavy testing when we fixed<br>
the problem of centered pixel coordinates.<br>
<br>
The Dashboard was green when we released<br>
ITK 3.18, which give us some confidence on<br>
the correctness of the code.<br>
<br>
--<br>
<br>
Please note that it is your responsibility to call<br>
<br>
interpolator->IsInsideBuffer( point )<br>
<br>
<br>
</blockquote>
<br>
There is a bug in the interpolator API. The expectation of the API is that if the interpolator can legally address the coordinates of the point then all of the surrounding discrete indexes can be legally addressed. Of course this is not always true, but last time I looked, none of the interpolators check that the discrete coordinates they look up are still inside the image. It is a problem for interpolation near the boundaries of the image. This is particularly a problem for interpolators with a large window region.<br>
<br>
--Simon<br>
<br></blockquote></div>------------------------------------------------------------------<br><br><br>Simon,<br><br>Thanks for pointing this out.<br><br>This seems to be the same problem that Ivan is reporting.<br><br><br>
Dzenan,<br><br>Please to log a bug entry for this issue:<br><a href="http://public.kitware.com/Bug/my_view_page.php" target="_blank">http://public.kitware.com/Bug/my_view_page.php</a><br><br><br> Thanks<br><font color="#888888"><br>
<br> Luis<br>
<br><br>
</font></blockquote></div><br>Hi again,<br><br>Just to let Simon know, I answered in another post that it was my fault this time, sorry. A bounds check is performed when using ITK_USE_CENTERED_PIXEL_COORDINATES_CONSISTENTLY. This was introduced by Luis in revision 1.42 of itkLinearInterpolateImageFunction.txx. However there seems to be no check when this is not defined. I will check my old failing programs again.<br>
<br>Best regards<br><br>Iván<br><br><br><br><br>