[Insight-users] Multiresolution registration with different timestep on each level

[Danielle Pace]

Hi all,

I have a deformable registration algorithm implementing by deriving a filter
from itk::PDEDeformableRegistrationFilter and a function from
itk::PDEDeformableRegistrationFunction.  I use it alongside
itk::MultiResolutionPDEDeformableRegistration and
itk::RecursiveMultiResolutionPyramidImageFilter to accomplish
multi-resolution registration with three levels.

The numerical stability of my deformable registration method's
regularization is dependent on the image spacing.  For images with higher
spacing, I can use larger timesteps without losing accuracy.  So, I would
like to make the registration run faster by using larger timesteps at the
first registration levels.

Is there a particular reason why there is no way to specify different
timesteps for each level in itk::MultiResolutionPDEDeformableRegistration?

The API could be similar to that for the number of iterations:
virtual void itk::MultiResolutionPDEDeformableRegistration<http://www.itk.org/Doxygen318/html/classitk_1_1MultiResolutionPDEDeformableRegistration.html><
TFixedImage, TMovingImage, TDeformationField, TRealType
>::SetNumberOfIterations(unsigned int *data*[] ) [virtual]
In addition, we could add a function SetTimesteps that would input an array
of timesteps, one for each level.

Is this something that people would find useful?

Thanks,

Danielle


-- 
Danielle Pace, M.ESc.
Research and Development Engineer
Kitware Inc.,
North Carolina Office

www.kitware.com
919-969-6990 X 319
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20110713/53bdebbc/attachment.htm>