Python 3.7 back compatible with 2.7.
Research code for testing and comparing numerical integrations of vector fields in in the infinite dimensional Lie group / Lie algebra of diffeomorphisms setting.
An axial slice from the BrainWeb dataset - Subject BW 38 is transformed according to:
- Translation
- Rotation of
- Unstable node (linear transformation with real positive eigenvalues)
- Inward spiral (linear transformation, complex conjugate eigenvectors with negative real part)
- Random homographic transformation
All transformations are parametrised with stationary velocity field in Lagrangian coordinates (in red), whose integral or flow field between 0 and 1 represents the actual transformation (in blue).
A stationary velocity field is an elements of the infinite dimensional Lie algebra of diffeomorphisms that parametrise the transformation in an Euclidean space. A diffeomorphism is an element of the infinite dimensional Lie group over , subset of . In the figure above, a stationary velocity field is represented with an arrow of the tangent space of the group of diffeomorphisms. Lie exponential and Lie logarithm map the vector field in the corresponding flow and vice versa.
The code is licenced under BSD 3-Clause.
- This repository is developed within the GIFT-surg research project.
- This work was supported by Wellcome / Engineering and Physical Sciences Research Council (EPSRC) [WT101957; NS/A000027/1; 203145Z/16/Z]. Sebastiano Ferraris was supported by the EPSRC-funded UCL Centre for Doctoral Training in Medical Imaging (EP/L016478/1) and Doctoral Training Grant (EP/M506448/1).