[ITK] [ITK-users] Regularization of itkDemonsRegistration

alexlovesmath atayyyy at gmail.com
Fri Mar 21 10:30:32 EDT 2014


All,

I appear to be having troubles with my ITK-users login. I apologize in
advance if I have mistakenly sent multiple copies of this post to the group.

Here is my question:

I am using ITK 4.2 and have a question about the implementation of
itkDemonsRegistrationFilter. 

My question relates to the order of regularization of the accumulated field
in each iteration of Demons. 

What I conjecture is happening in itkDemonsRegistrationFilter is that
regularization of the accumulated field must happen at the beginning of each
iteration. I believe this implies that in the last iteration of the filter,
the displacement field returned by filter->GetOutput() is the result of the
composition of a regularized accumulated field from all past iterations and
the update field, but that a final regularization is not performed. 

I have tested my hypothesis by running the DeminsRegistrationFilter for one
iteration starting with the identity transformation (all zeros) as the
initial condition for the displacement field, once with
SmoothDisplacementFieldOff() and once with SmoothDisplacementFieldOn() (I'm
using the default where SmoothUpdateField is off). In both cases, I obtain
the same result. I think this has to imply that the regularization of the
accumulated field is done at the start of each iteration. 

My questions: 

1) Am I correct that at the end of each iteration, the previously
regularized accumulated field is composed with the update field and that
this is the result returned by filter->GetOutput(). 

2) Is there some reason not to regularize the accumulated field "on the way
out" at the end of each iteration? Is there some algorithmic advantage to
this? I do understand that if a non-identity displacement field was
specified as the initial condition, that the current architecture has the
advantage of regularizing the initial condition. However, it is possible to
regularize the initial condition outside the iterative loop and still
achieve regularization of the accumulated field at the end of the last
iteration. 

Thanks! 

- Alex



--
View this message in context: http://itk-insight-users.2283740.n2.nabble.com/Regularization-of-itkDemonsRegistration-tp7585217.html
Sent from the ITK Insight Users mailing list archive at Nabble.com.
_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.php

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users


More information about the Community mailing list