<div dir="ltr">Hi Matt,<div><br></div><div style>Thanks a lot for this explanation :)</div><div style><br></div><div style>AtB</div><div style><br></div><div style>Anja</div><div style><br></div><div class="gmail_extra"><br>
<br><div class="gmail_quote">On 10 April 2013 20:28, Matt McCormick <span dir="ltr">&lt;<a href="mailto:matt.mccormick@kitware.com" target="_blank">matt.mccormick@kitware.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Anja,<br>
<br>
Yes, sampling and the grid play an important role.<br>
<br>
The sampling locations of the output must be chosen, and they are not<br>
necessarily the sampling locations of the input.<br>
<br>
This is what the ResampleImageFilter and its Transform may operator in<br>
the &quot;opposite&quot; direction what is expected -- what is always computed<br>
is the output pixel location.<br>
<br>
HTH,<br>
Matt<br>
<div class="HOEnZb"><div class="h5"><br>
On Wed, Apr 10, 2013 at 7:17 PM, Anja Ende &lt;<a href="mailto:anja.ende@googlemail.com">anja.ende@googlemail.com</a>&gt; wrote:<br>
&gt; Oh yes, of course! That is a subtle point.<br>
&gt;<br>
&gt; May I ask a follow up question then? In inverting a deformation field then,<br>
&gt; how do we decide the sampling of this grid of the inverse field?<br>
&gt;<br>
&gt; Say our original deformation field has a displacement at every voxel and<br>
&gt; these displacement vectors can point to any real location. So, is the<br>
&gt; inverse deformation field going to be calculated at these real locations?<br>
&gt; However, we need to have a uniform grid and this is not guaranteed by this<br>
&gt; scheme now. Can someone explain this to me on how this is done?<br>
&gt;<br>
&gt; Sorry for asking these naive questions. I could use the ITK filters out of<br>
&gt; the box but I would like to know some of these basics.<br>
&gt;<br>
&gt; Thanks,<br>
&gt; Anja<br>
&gt;<br>
&gt;<br>
&gt; On 10 April 2013 20:06, Matt McCormick &lt;<a href="mailto:matt.mccormick@kitware.com">matt.mccormick@kitware.com</a>&gt; wrote:<br>
&gt;&gt;<br>
&gt;&gt; Hi Anja,<br>
&gt;&gt;<br>
&gt;&gt; If displacement vector field gives the displacement at location x, but<br>
&gt;&gt; the negative of the displacement at location x is not the inverse;<br>
&gt;&gt; what is required is negative of the displacement vector at x + \delta<br>
&gt;&gt; x.<br>
&gt;&gt;<br>
&gt;&gt; HTH,<br>
&gt;&gt; Matt<br>
&gt;&gt;<br>
&gt;&gt; On Wed, Apr 10, 2013 at 5:47 PM, Anja Ende &lt;<a href="mailto:anja.ende@googlemail.com">anja.ende@googlemail.com</a>&gt;<br>
&gt;&gt; wrote:<br>
&gt;&gt; &gt; Hello all,<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; This is a very naive question and I would be very grateful if someone<br>
&gt;&gt; &gt; can<br>
&gt;&gt; &gt; clarify this doubt for me.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; I see that there are quite a few filters and some very sophisticated<br>
&gt;&gt; &gt; methods<br>
&gt;&gt; &gt; for inverting deformation fields.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; If I have a dense deformation field, why is it not enough to simply<br>
&gt;&gt; &gt; negate<br>
&gt;&gt; &gt; the displacement vector to compute the inverse? I am sure there is a<br>
&gt;&gt; &gt; simple,<br>
&gt;&gt; &gt; logical explanation for this but I am having trouble figuring it out.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; Thanks a lot for your help!<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; Anja<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; _____________________________________<br>
&gt;&gt; &gt; Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; Visit other Kitware open-source projects at<br>
&gt;&gt; &gt; <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; Kitware offers ITK Training Courses, for more information visit:<br>
&gt;&gt; &gt; <a href="http://www.kitware.com/products/protraining.php" target="_blank">http://www.kitware.com/products/protraining.php</a><br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; Please keep messages on-topic and check the ITK FAQ at:<br>
&gt;&gt; &gt; <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; Follow this link to subscribe/unsubscribe:<br>
&gt;&gt; &gt; <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
&gt;&gt; &gt;<br>
&gt;<br>
&gt;<br>
&gt;</div></div></blockquote></div></div></div>