Diffeomorphometry

From HandWiki
Revision as of 23:05, 6 March 2023 by S.Timg (talk | contribs) (url)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Short description: Metric study of shape and form in computational anatomy

Diffeomorphometry is the metric study of imagery, shape and form in the discipline of computational anatomy (CA) in medical imaging. The study of images in computational anatomy rely on high-dimensional diffeomorphism groups [math]\displaystyle{ \varphi \in \operatorname{Diff}_V }[/math] which generate orbits of the form [math]\displaystyle{ \mathcal{I} \doteq \{ \varphi \cdot I \mid \varphi \in \operatorname{Diff}_V \} }[/math], in which images [math]\displaystyle{ I \in \mathcal{I} }[/math] can be dense scalar magnetic resonance or computed axial tomography images. For deformable shapes these are the collection of manifolds [math]\displaystyle{ \mathcal{M} \doteq \{ \varphi \cdot M \mid \varphi \in \operatorname{Diff}_V \} }[/math], points, curves and surfaces. The diffeomorphisms move the images and shapes through the orbit according to [math]\displaystyle{ (\varphi,I)\mapsto \varphi \cdot I }[/math] which are defined as the group actions of computational anatomy.

The orbit of shapes and forms is made into a metric space by inducing a metric on the group of diffeomorphisms. The study of metrics on groups of diffeomorphisms and the study of metrics between manifolds and surfaces has been an area of significant investigation.[1][2][3][4][5][6][7][8][9] In Computational anatomy, the diffeomorphometry metric measures how close and far two shapes or images are from each other. Informally, the metric is constructed by defining a flow of diffeomorphisms [math]\displaystyle{ \dot \phi_t , t \in [0,1], \phi_t \in \operatorname{Diff}_V }[/math] which connect the group elements from one to another, so for [math]\displaystyle{ \varphi,\psi \in \operatorname{Diff}_V }[/math] then [math]\displaystyle{ \phi_0 = \varphi , \phi_1=\psi }[/math]. The metric between two coordinate systems or diffeomorphisms is then the shortest length or geodesic flow connecting them. The metric on the space associated to the geodesics is given by[math]\displaystyle{ \rho(\varphi,\psi) = \inf_{\phi: \phi_0=\varphi,\phi_1 = \psi} \int_0^1 \| \dot \phi_t \|_{\phi_t} \, dt }[/math]. The metrics on the orbits [math]\displaystyle{ \mathcal{I},\mathcal{M} }[/math] are inherited from the metric induced on the diffeomorphism group.

The group [math]\displaystyle{ \varphi \in \operatorname{Diff}_V }[/math] is thusly made into a smooth Riemannian manifold with Riemannian metric [math]\displaystyle{ \| \cdot \|_\varphi }[/math] associated to the tangent spaces at all [math]\displaystyle{ \varphi \in\operatorname{Diff}_V }[/math]. The Riemannian metric satisfies at every point of the manifold [math]\displaystyle{ \phi \in \operatorname{Diff}_V }[/math] there is an inner product inducing the norm on the tangent space [math]\displaystyle{ \| \dot \phi_t \|_{\phi_t} }[/math] that varies smoothly across [math]\displaystyle{ \operatorname{Diff}_V }[/math].

Oftentimes, the familiar Euclidean metric is not directly applicable because the patterns of shapes and images don't form a vector space. In the Riemannian orbit model of Computational anatomy, diffeomorphisms acting on the forms [math]\displaystyle{ \varphi \cdot I \in \mathcal {I}, \varphi \in \operatorname{Diff}_V, M \in \mathcal{M} }[/math] don't act linearly. There are many ways to define metrics, and for the sets associated to shapes the Hausdorff metric is another. The method used to induce the Riemannian metric is to induce the metric on the orbit of shapes by defining it in terms of the metric length between diffeomorphic coordinate system transformations of the flows. Measuring the lengths of the geodesic flow between coordinates systems in the orbit of shapes is called diffeomorphometry.

The diffeomorphisms group generated as Lagrangian and Eulerian flows

The diffeomorphisms in computational anatomy are generated to satisfy the Lagrangian and Eulerian specification of the flow fields, [math]\displaystyle{ \varphi_t, t \in [0,1] }[/math], generated via the ordinary differential equation

[math]\displaystyle{ \frac{d}{dt} \varphi_t = v_t \circ \varphi_t , \ \varphi_0 = \operatorname{id}; }[/math]

 

 

 

 

(Lagrangian flow)

with the Eulerian vector fields [math]\displaystyle{ v \doteq (v_1,v_2,v_3) }[/math] in [math]\displaystyle{ {\mathbb R}^3 }[/math] for [math]\displaystyle{ v_t = \dot \varphi_t \circ \varphi_t^{-1}, t \in [0,1] }[/math]. The inverse for the flow is given by [math]\displaystyle{ \frac{d}{dt} \varphi_t^{-1} = -(D \varphi_t^{-1}) v_t, \ \varphi_0^{-1} = \operatorname{id}, }[/math] and the [math]\displaystyle{ 3 \times 3 }[/math] Jacobian matrix for flows in [math]\displaystyle{ \mathbb{R}^3 }[/math] given as [math]\displaystyle{ \ D\varphi \doteq \left(\frac{\partial \varphi_i}{\partial x_j}\right). }[/math]

To ensure smooth flows of diffeomorphisms with inverse, the vector fields [math]\displaystyle{ {\mathbb R}^3 }[/math] must be at least 1-time continuously differentiable in space[10][11] which are modelled as elements of the Hilbert space [math]\displaystyle{ (V, \| \cdot \|_V ) }[/math] using the Sobolev embedding theorems so that each element [math]\displaystyle{ v_i \in H_0^3, i=1,2,3, }[/math] has 3-square-integrable derivatives thusly implies [math]\displaystyle{ (V, \| \cdot \|_V ) }[/math] embeds smoothly in 1-time continuously differentiable functions.[10][11] The diffeomorphism group are flows with vector fields absolutely integrable in Sobolev norm:

[math]\displaystyle{ \operatorname{Diff}_V \doteq \{\varphi=\varphi_1: \dot \varphi_t = v_t \circ \varphi_t , \varphi_0 = \operatorname{id}, \int_0^1 \|v_t \|_V \,dt \lt \infty \} \ . }[/math]

 

 

 

 

(Diffeomorphism Group)

The Riemannian orbit model

Shapes in Computational Anatomy (CA) are studied via the use of diffeomorphic mapping for establishing correspondences between anatomical coordinate systems. In this setting, 3-dimensional medical images are modelled as diffemorphic transformations of some exemplar, termed the template [math]\displaystyle{ I_{temp} }[/math], resulting in the observed images to be elements of the random orbit model of CA. For images these are defined as [math]\displaystyle{ I \in \mathcal {I} \doteq \{ I = I_{temp} \circ \varphi, \varphi \in \operatorname{Diff}_V \} }[/math], with for charts representing sub-manifolds denoted as [math]\displaystyle{ \mathcal{M} \doteq \{ \varphi \cdot M_{temp} : \varphi \in \operatorname{Diff}_V \} }[/math].

The Riemannian metric

The orbit of shapes and forms in Computational Anatomy are generated by the group action [math]\displaystyle{ \mathcal{I} \doteq \{ \varphi \cdot I : \varphi \in \operatorname{Diff}_V \} }[/math] , [math]\displaystyle{ \mathcal{M} \doteq \{ \varphi \cdot M : \varphi \in \operatorname{Diff}_V \} }[/math]. These are made into a Riemannian orbits by introducing a metric associated to each point and associated tangent space. For this a metric is defined on the group which induces the metric on the orbit. Take as the metric for Computational anatomy at each element of the tangent space [math]\displaystyle{ \varphi \in \operatorname{Diff}_V }[/math] in the group of diffeomorphisms

[math]\displaystyle{ \| \dot \varphi \|_\varphi \doteq \| \dot \varphi \circ \varphi^{-1} \|_V=\| v \|_V, }[/math]

with the vector fields modelled to be in a Hilbert space with the norm in the Hilbert space [math]\displaystyle{ (V, \| \cdot \|_V ) }[/math]. We model [math]\displaystyle{ V }[/math] as a reproducing kernel Hilbert space (RKHS) defined by a 1-1, differential operator [math]\displaystyle{ A: V \rightarrow V^* }[/math], where [math]\displaystyle{ V^* }[/math] is the dual-space. In general, [math]\displaystyle{ \sigma \doteq Av \in V^* }[/math] is a generalized function or distribution, the linear form associated to the inner-product and norm for generalized functions are interpreted by integration by parts according to for [math]\displaystyle{ v,w \in V }[/math],

[math]\displaystyle{ \langle v , w \rangle_V \doteq \int_X A v \cdot w \, dx, \ \| v\|_V^2 \doteq \int_X A v \cdot v \, dx, \ v,w \in V \ . }[/math]

When [math]\displaystyle{ Av \doteq \mu \,dx }[/math], a vector density, [math]\displaystyle{ \int Av \cdot v \,dx \doteq \int \mu \cdot v \, dx = \sum_{i=1}^3 \mu_i v_i \, dx. }[/math]

The differential operator is selected so that the Green's kernel associated to the inverse is sufficiently smooth so that the vector fields support 1-continuous derivative. The Sobolev embedding theorem arguments were made in demonstrating that 1-continuous derivative is required for smooth flows. The Green's operator generated from the Green's function(scalar case) associated to the differential operator smooths.

For proper choice of [math]\displaystyle{ A }[/math] then [math]\displaystyle{ (V,\| \cdot \|_V) }[/math] is an RKHS with the operator [math]\displaystyle{ K = A^{-1}: V^* \rightarrow V }[/math]. The Green's kernels associated to the differential operator smooths since for controlling enough derivatives in the square-integral sense the kernel [math]\displaystyle{ k(\cdot,\cdot) }[/math] is continuously differentiable in both variables implying

[math]\displaystyle{ K Av (x)_i \doteq \sum_j \int_{{\mathbb R}^3} k_{ij}(x,y) Av_j(y) \,dy \in V \ . }[/math]

The diffeomorphometry of the space of shapes and forms

The right-invariant metric on diffeomorphisms

The metric on the group of diffeomorphisms is defined by the distance as defined on pairs of elements in the group of diffeomorphisms according to

[math]\displaystyle{ d_{\mathrm{Diff}_V}(\psi, \varphi) = \inf_{v_t} \left(\int_0^1 \int_X Av_t \cdot v_t \, dx \, dt: \phi_0 = \psi, \phi_1 = \varphi, \dot \phi_t = v_t \circ \phi_t \right)^{1/2} \ . }[/math]

 

 

 

 

(metric-diffeomorphisms)

This distance provides a right-invariant metric of diffeomorphometry,[12][13][14] invariant to reparameterization of space since for all [math]\displaystyle{ \phi \in \operatorname{Diff}_V }[/math],

[math]\displaystyle{ d_{\operatorname{Diff}_V}(\psi, \varphi) = d_{\operatorname{Diff}_V}(\psi \circ \phi, \varphi \circ \phi). }[/math]

The metric on shapes and forms

The distance on images,[15] [math]\displaystyle{ d_{\mathcal{I}}:\mathcal{I} \times \mathcal{I}\rightarrow \R^+ }[/math],


[math]\displaystyle{ d_{\mathcal{I}}(I,J)=\inf_{\phi \in \operatorname{Diff}_V: \phi \cdot I = J } d_{\operatorname{Diff}_V}(id,\phi) \ ; }[/math]

 

 

 

 

(metric-shapes-forms)

The distance on shapes and forms,[16] [math]\displaystyle{ d_{\mathcal{M}}:\mathcal{M} \times \mathcal{M}\rightarrow \R^+ }[/math],


[math]\displaystyle{ d_{\mathcal{M}} (M,N) = \inf_{\phi \in \operatorname{Diff}_V: \phi \cdot M = N } d_{\mathrm{Diff}_V}(\operatorname{id},\phi) \ . }[/math]

 

 

 

 

(metric-shapes-forms)

The metric on geodesic flows of landmarks, surfaces, and volumes within the orbit

For calculating the metric, the geodesics are a dynamical system, the flow of coordinates [math]\displaystyle{ t \mapsto \phi_t \in \operatorname{Diff}_V }[/math] and the control the vector field [math]\displaystyle{ t \mapsto v_t \in V }[/math] related via [math]\displaystyle{ \dot \phi_t = v_t \cdot \phi_t,\phi_0=\operatorname{id}. }[/math] The Hamiltonian view [17] [18] [19] [20][21] reparameterizes the momentum distribution [math]\displaystyle{ Av \in V^* }[/math] in terms of the Hamiltonian momentum, a Lagrange multiplier [math]\displaystyle{ p: \dot \phi \mapsto (p\mid\dot \phi) }[/math] constraining the Lagrangian velocity [math]\displaystyle{ \dot \phi_t = v_t \circ \phi_t }[/math].accordingly:

[math]\displaystyle{ H(\phi_t,p_t,v_t)=\int_X p_t \cdot (v_t \circ \phi_t) \, dx-\frac{1}{2}\int_X Av_t \cdot v_t \, dx . }[/math]

The Pontryagin maximum principle[17] gives the Hamiltonian [math]\displaystyle{ H(\phi_t,p_t) \doteq \max_v H( \phi_t, p_t,v) \ . }[/math] The optimizing vector field [math]\displaystyle{ v_t \doteq \operatorname{argmax}_v H(\phi_t,p_t,v) }[/math] with dynamics [math]\displaystyle{ \dot \phi_t = \frac{\partial H( \phi_t, p_t)}{\partial p}, \dot p_t = -\frac{\partial H(\phi_t,p_t)}{\partial \phi} }[/math]. Along the geodesic the Hamiltonian is constant:[22] [math]\displaystyle{ H(\phi_t,p_t) = H(\operatorname{id},p_0)=\frac{1}{2} \int_X p_0 \cdot v_0 \, dx }[/math]. The metric distance between coordinate systems connected via the geodesic determined by the induced distance between identity and group element:

[math]\displaystyle{ d_{\mathrm{Diff}_V}(\operatorname{id},\varphi) =\| v_0 \|_V = \sqrt{2H(\operatorname{id},p_0)} }[/math]

Landmark or pointset geodesics

For landmarks, [math]\displaystyle{ x_i, i=1,\dots,n }[/math], the Hamiltonian momentum

[math]\displaystyle{ p(i), i=1,\dots,n }[/math]

with Hamiltonian dynamics taking the form

[math]\displaystyle{ H(\phi_t,p_t) =\frac{1}{2}\textstyle \sum_j \sum_i \displaystyle p_t(i)\cdot K(\phi_t (x_i),\phi_t (x_j)) p_t(j) }[/math]

with

[math]\displaystyle{ \begin{cases} v_t = \textstyle \sum_i \displaystyle K(\cdot, \phi_t (x_i)) p_t(i) , \ \\ \dot p_t (i) = - (Dv_t)^T_{|_{\phi_t(x_i)}} p_t(i), i=1,2,\dots, n \\ \end{cases} }[/math]

The metric between landmarks [math]\displaystyle{ d^2 =\textstyle \sum_i p_0(i)\cdot \sum_j \displaystyle K(x_i,x_j) p_0(j). }[/math]

The dynamics associated to these geodesics is shown in the accompanying figure.

Surface geodesics

For surfaces, the Hamiltonian momentum is defined across the surface has Hamiltonian

[math]\displaystyle{ H(\phi_t,p_t) =\frac{1}{2} \int_U \int_U p_t(u)\cdot K(\phi_t (m(u)), \phi_t (m(v))) p_t(v) \, du \, dv }[/math]

and dynamics

[math]\displaystyle{ \begin{cases} v_t= \textstyle \int_U \displaystyle K(\cdot, \phi_t ( m(u)))p_t(u)\,du \ , \\ \dot p_t(u) = - (Dv_t)^T_{|_{\phi_t(m(u))} } p_t(u), u \in U \end{cases} }[/math]
The metric between surface coordinates [math]\displaystyle{ d^2 = (p_0 \mid v_0) =\int_U p_0(u) \cdot \int_U K(m(u), m(u^\prime)) p_0(u^\prime) \, du \, du^\prime }[/math]

Volume geodesics

For volumes the Hamiltonian

[math]\displaystyle{ H(\phi_t,p_t) = \frac{1}{2}\int_{{\mathbb R}^3} \int_{{\mathbb R}^3} p_t(x)\cdot K(\phi_t(x),\phi_t(y)) p_t(y) \, dx \, dy \displaystyle }[/math]

with dynamics

[math]\displaystyle{ \begin{cases} v_t=\textstyle \int_X \displaystyle K(\cdot, \phi_t(x))p_t(x)\,dx \ , \\ \dot p_t(x) = - (Dv_t)^T_{|_{\phi_t(x)} } p_t(x), x \in {\mathbb R}^3 \end{cases} }[/math]
The metric between volumes [math]\displaystyle{ \displaystyle d^2 =(p_0\mid v_0) = \int_{\mathbb R^3} p_0(x)\cdot \int_{{\mathbb R}^3} K(x,y) p_0(y)\,dy \, dx. }[/math]

Software for diffeomorphic mapping

Software suites containing a variety of diffeomorphic mapping algorithms include the following:

Cloud software

References

  1. Miller, M. I.; Younes, L. (2001-01-01). "Group Actions, Homeomorphisms, and Matching: A General Framework" (in en). International Journal of Computer Vision 41 (1–2): 61–84. doi:10.1023/A:1011161132514. ISSN 0920-5691. 
  2. Younes, L. (1998-04-01). "Computable Elastic Distances Between Shapes". SIAM Journal on Applied Mathematics 58 (2): 565–586. doi:10.1137/S0036139995287685. 
  3. Mio, Washington; Srivastava, Anuj; Joshi, Shantanu (2006-09-25). "On Shape of Plane Elastic Curves". International Journal of Computer Vision 73 (3): 307–324. doi:10.1007/s11263-006-9968-0. 
  4. Michor, Peter W.; Mumford, David; Shah, Jayant; Younes, Laurent (2008). "A Metric on Shape Space with Explicit Geodesics". Rend. Lincei Mat. Appl. () 9 (2008): 25–57. Bibcode2007arXiv0706.4299M. 
  5. Michor, Peter W.; Mumford, David (2007). "An overview of the Riemannian metrics on spaces of curves using the Hamiltonian approach". Applied and Computational Harmonic Analysis 23 (1): 74–113. doi:10.1016/j.acha.2006.07.004. 
  6. Kurtek, Sebastian; Klassen, Eric; Gore, John C.; Ding, Zhaohua; Srivastava, Anuj (2012-09-01). "Elastic geodesic paths in shape space of parameterized surfaces". IEEE Transactions on Pattern Analysis and Machine Intelligence 34 (9): 1717–1730. doi:10.1109/TPAMI.2011.233. PMID 22144521. 
  7. Srivastava, Anuj; Klassen, Eric; Joshi, Shantanu H.; Jermyn, Ian H. (2011). "Shape Analysis of Elastic Curves in Euclidean Spaces". IEEE Transactions on Pattern Analysis and Machine Intelligence 33 (7): 1415–1428. doi:10.1109/TPAMI.2010.184. ISSN 1939-3539. PMID 20921581. https://ieeexplore.ieee.org/document/5601739. 
  8. Jermyn, Ian H.; Kurtek, Sebastian; Klassen, Eric; Srivastava, Anuj (2012), Fitzgibbon, Andrew; Lazebnik, Svetlana; Perona, Pietro et al., eds., "Elastic Shape Matching of Parameterized Surfaces Using Square Root Normal Fields" (in en), Computer Vision – ECCV 2012 (Berlin, Heidelberg: Springer Berlin Heidelberg) 7576: pp. 804–817, doi:10.1007/978-3-642-33715-4_58, ISBN 978-3-642-33714-7 
  9. Jermyn, Ian H.; Kurtek, Sebastian; Laga, Hamid; Srivastava, Anuj (2017-09-15). "Elastic Shape Analysis of Three-Dimensional Objects" (in en). Synthesis Lectures on Computer Vision 7 (3): 1–185. doi:10.2200/s00785ed1v01y201707cov012. ISSN 2153-1056. https://doi.org/10.2200/S00785ED1V01Y201707COV012. 
  10. 10.0 10.1 P. Dupuis, U. Grenander, M.I. Miller, Existence of Solutions on Flows of Diffeomorphisms, Quarterly of Applied Math, 1997.
  11. 11.0 11.1 A. Trouvé. Action de groupe de dimension infinie et reconnaissance de formes. C R Acad Sci Paris Sér I Math, 321(8):1031– 1034, 1995.
  12. Miller, M. I.; Younes, L. (2001-01-01). "Group Actions, Homeomorphisms, And Matching: A General Framework". International Journal of Computer Vision 41: 61–84. doi:10.1023/A:1011161132514. 
  13. Miller, M. I; Younes, L; Trouvé, A (2014). "Diffeomorphometry and geodesic positioning systems for human anatomy". Technology 2 (1): 36. doi:10.1142/S2339547814500010. PMID 24904924. 
  14. Miller, Michael I.; Trouvé, Alain; Younes, Laurent (2015-01-01). "Hamiltonian Systems and Optimal Control in Computational Anatomy: 100 Years Since D'Arcy Thompson". Annual Review of Biomedical Engineering 17 (1): 447–509. doi:10.1146/annurev-bioeng-071114-040601. PMID 26643025. 
  15. Miller, M. I.; Younes, L. (2001-01-01). "Group Actions, Homeomorphisms, And Matching: A General Framework". International Journal of Computer Vision 41: 61–84. doi:10.1023/A:1011161132514. 
  16. Miller, Michael I.; Younes, Laurent; Trouvé, Alain (March 2014). "Diffeomorphometry and geodesic positioning systems for human anatomy". Technology 2 (1): 36. doi:10.1142/S2339547814500010. ISSN 2339-5478. PMID 24904924. 
  17. 17.0 17.1 Miller, Michael I.; Trouvé, Alain; Younes, Laurent (2015-01-01). "Hamiltonian Systems and Optimal Control in Computational Anatomy: 100 Years Since D'arcy Thompson". Annual Review of Biomedical Engineering 17 (1): 447–509. doi:10.1146/annurev-bioeng-071114-040601. PMID 26643025. 
  18. Glaunès J, Trouvé A, Younes L. 2006. Modeling planar shape variation via Hamiltonian flows of curves. In Statistics and Analysis of Shapes, ed. H Krim, A Yezzi Jr, pp. 335–61. Model. Simul. Sci. Eng. Technol. Boston: Birkhauser
  19. Arguillère S, Trélat E, Trouvé A, Younes L. 2014. Shape deformation analysis from the optimal control viewpoint. arXiv:1401.0661 [math.OC]
  20. Miller, MI; Younes, L; Trouvé, A (2014). "Diffeomorphometry and geodesic positioning systems for human anatomy". Technology (Singap World Sci) 2 (1): 36. doi:10.1142/S2339547814500010. PMID 24904924. 
  21. Michor, Peter W.; Mumford, David (2007-07-01). "An overview of the Riemannian metrics on spaces of curves using the Hamiltonian approach". Applied and Computational Harmonic Analysis. Special Issue on Mathematical Imaging 23 (1): 74–113. doi:10.1016/j.acha.2006.07.004. 
  22. Miller, Michael I.; Trouvé, Alain; Younes, Laurent (2015-01-01). "Hamiltonian Systems and Optimal Control in Computational Anatomy: 100 Years Since D'Arcy Thompson". Annual Review of Biomedical Engineering 17 (1): 447–509. doi:10.1146/annurev-bioeng-071114-040601. PMID 26643025. 
  23. Software - Stanley Durrleman. https://raweb.inria.fr/rapportsactivite/RA2015/aramis/uid32.html. 
  24. Avants, Brian B.; Tustison, Nicholas J.; Song, Gang; Cook, Philip A.; Klein, Arno; Gee, James C. (2011-02-01). "A Reproducible Evaluation of ANTs Similarity Metric Performance in Brain Image Registration". NeuroImage 54 (3): 2033–2044. doi:10.1016/j.neuroimage.2010.09.025. ISSN 1053-8119. PMID 20851191. 
  25. Ashburner, John (2007-10-15). "A fast diffeomorphic image registration algorithm". NeuroImage 38 (1): 95–113. doi:10.1016/j.neuroimage.2007.07.007. PMID 17761438. 
  26. "Software - Tom Vercauteren". https://sites.google.com/site/tomvercauteren/software. Retrieved 2015-12-11. 
  27. Beg, M. Faisal; Miller, Michael I.; Trouvé, Alain; Younes, Laurent (2005-02-01). "Computing Large Deformation Metric Mappings via Geodesic Flows of Diffeomorphisms" (in en). International Journal of Computer Vision 61 (2): 139–157. doi:10.1023/B:VISI.0000043755.93987.aa. ISSN 0920-5691. 
  28. "Comparing algorithms for diffeomorphic registration: Stationary LDDMM and Diffeomorphic Demons (PDF Download Available)" (in en). https://www.researchgate.net/publication/33419970. 
  29. "MRICloud". The Johns Hopkins University. http://www.mricloud.org.