<div dir="ltr">see line 142 here:<div><br></div><div><a href="https://github.com/stnava/ANTs/blob/d99d9c9e23ae17c116377a248cd370ddbd8c2599/Examples/antsApplyTransformsToPoints.cxx">https://github.com/stnava/ANTs/blob/d99d9c9e23ae17c116377a248cd370ddbd8c2599/Examples/antsApplyTransformsToPoints.cxx</a><br>
</div><div><br></div><div>you can do the same thing in your application but note that the deformation field , when applied to a point , should be the inverse of the deformation field applied to an image</div><div><br></div>
<div>see my previous email regarding SyN and inverse warps.</div></div><div class="gmail_extra"><br clear="all"><div><div><br></div>brian<br><div><br></div><div><br></div></div>
<br><br><div class="gmail_quote">On Wed, Dec 4, 2013 at 2:33 PM, Yago Diez <span dir="ltr"><<a href="mailto:yagodiezdonoso@gmail.com" target="_blank">yagodiezdonoso@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div><div><div><div>Hi everyone,<br><br></div>I have a problem regarding itk deformation fields:<br><br></div>Specifically, I want to know where a particular point ends up after registration. Said registration has produced a deformation field. I have already used this deformation field to deform the source image and produce the output image for registration (this was done using itkWarpImageFilter). <br>
<br></div>My problem is that at this moment I only want to transform one point. The obvious solution was to read the corresponding voxel in the deformation field and add it to the point, but the problem is that the deformation field is not "dense enough", so in most cases, the read deformation for a point ends up being zero (see attached image for a screen capture of a 2D slice of my 3D deformation field).<br>
<br> I have noticed how when I deform an image with itkWarpImageFilter, there is an "extra" interpolation step which would avoid the problem I am having, but so far I have not been able to interpolate a deformation field. Any ideas on how to do that?<br>
<br></div>Thank you<br><div><br><br><div><div><div><div><br><br clear="all"><div><div dir="ltr">Yago Diez Donoso (PhD)<br>Computer Vision and Robotics Group <a href="http://vicorob.udg.es/" target="_blank">http://vicorob.udg.es/</a><br>
E-mail: <a href="mailto:yago@eia.udg.es" target="_blank">yago@eia.udg.es</a>; <a href="mailto:yagodiezdonoso@gmail.com" target="_blank">yagodiezdonoso@gmail.com</a><br>Phone: (int. code) 34 972418013<br>University of Girona </div>
</div>
</div></div></div></div></div></div>
<br>_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.php" target="_blank">http://www.kitware.com/products/protraining.php</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
<br></blockquote></div><br></div>