[Insight-users] Deformation field initialization with a rigid transform

[Torsten Rohlfing]

Hi Moti:

This is actually quite simple, at least as far as the demons 
registration filter is concerned. You can set an initial deformation 
field using the 
itk::DemonsRegistrationFilter::SetInitialDeformationField member function.

The only thing you have to do is to create a deformation field that 
represents your initial rigid transformation. For that, iterate over all 
image pixels, and at each pixel location v set the initial deformation 
field to R(v)-v, where R(.) is your rigid transformation.

It would probably be nice to have this functionality in ITK as a generic 
filter for arbitrary initial transformations, so maybe if you could 
implement one, you can submit it to the Insight Journal? :)

Best,
  Torsten

-- 

Torsten Rohlfing, PhD          SRI International, Neuroscience Program
 Research Scientist             333 Ravenswood Ave, Menlo Park, CA 94025
  Phone: ++1 (650) 859-3379      Fax: ++1 (650) 859-2743
   torsten at synapse.sri.com        http://www.stanford.edu/~rohlfing/

     "Though this be madness, yet there is a method in't"