<div dir="ltr"><div><div><div><div>Hello Mathew,<br><br></div>To my knowledge there is no way to get the inverse transform directly when using ITKDiffeomorphicDemonsRegistrationFilter. You could try to invert the deformation field using one of the algorithm that are available in ITK for that purpose. It seems that you tried one (InverDisplacementField) but have you tried the remote module [1]<br></div>If you want to get the inverse transformation directly, you can also try to use [2] or use the registration method that has been implemented in this Insight Journal paper [3].<br><br></div>Hope this helps,<br></div>Francois<br><br>[1] <a href="https://github.com/InsightSoftwareConsortium/ITK/blob/master/Modules/Remote/FixedPointInverseDisplacementField.remote.cmake">https://github.com/InsightSoftwareConsortium/ITK/blob/master/Modules/Remote/FixedPointInverseDisplacementField.remote.cmake</a><br>[2] <a href="https://itk.org/Doxygen/html/classitk_1_1VariationalSymmetricDiffeomorphicRegistrationFilter.html">https://itk.org/Doxygen/html/classitk_1_1VariationalSymmetricDiffeomorphicRegistrationFilter.html</a><br>[3] <a href="http://www.insight-journal.org/browse/publication/644">http://www.insight-journal.org/browse/publication/644</a><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, May 22, 2017 at 3:25 PM, Mathew Guilfoyle <span dir="ltr"><<a href="mailto:mrguilfoyle@gmail.com" target="_blank">mrguilfoyle@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I am using (Simple)ITK DiffeomorphicDemonsRegistratio<wbr>nFilter() to perform deformable registration of two 3D images.  This works fine.  Output from the filter is 'outImage' vector image which I then convert to a transform using<br>
<br>
dfTransform = DisplacementFieldTransform(<wbr>outImage).<br>
<br>
I can then resample my moving image using this transform and all looks as expected.<br>
<br>
My reading/understanding of diffeomorphic algorithms is that there should be an exact and readily computable inverse transform?  I need this as I want to inverse warp a template segmentation performed on my reference image into the original space of my moving image.  How do I obtain the inverse DisplacementFieldTransform or the inverse vector image?<br>
<br>
I have tried the InvertDisplacementField filter but the output of this does not appear to correctly invert my warped image back to the original so I do not think this filter generates is a n accurate inverse, at least using default settings.<br>
<br>
Any ideas of if/how I can get an exact inverse of a deformable transform?<br>
<br>
Thanks<br>
<br>
<br>
______________________________<wbr>_______<br>
Powered by <a href="http://www.kitware.com" rel="noreferrer" 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" rel="noreferrer" target="_blank">http://www.kitware.com/<wbr>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" rel="noreferrer" target="_blank">http://www.kitware.com/<wbr>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" rel="noreferrer" target="_blank">http://www.itk.org/Wiki/ITK_<wbr>FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/insight-users" rel="noreferrer" target="_blank">http://public.kitware.com/<wbr>mailman/listinfo/insight-users</a><br>
</blockquote></div><br></div>