Application to Glioblatomas
Anton François
suppervised by
Joan Glaunès & Pietro Gori
About Me :
2012-2015
2014-2016
2016-2019
2019-2022
Statistical Altas : [Roux, 2019] Frequence of apparition by voxels.
Anatomical Atlas : [Beg and Khan, 2005] Average shape.
Topological Altas : Collection of homeomorphisms covering a topological space (ex: manifold)
Normalisation of medical images is a important step for data acquisition
Before
linear registration :
After
Diffeomorphic matching => Keeps the Topology
Constant Vector fields
Temporal vector field
LDDMM
Comparing the amount of deformation in between two images.
LDDMM
Image I deformed by a temporal vector field v=(vt)t∈[0,1].
I˙t=vt⋅It
LDDMM
Energy of the deformation generated by v
E(v)≐21∫01∥vt∥V2dt=21∫01⟨vt,Kσ⋆vt⟩2dt
Kσ⋆ being the convolution with a Gaussian kernel.
Distance between two images
d(I,J)=minvE(v),s.t. I0=I; I1=∫01vt⋅∇Itdt=J
Geodesic := shortest path on a Manifold
To find the minimum of this cost :
E(v,I)=21∥I1−T∥22+21∫01∥vt∥V2 dt
We have to integrate over this set of geodesic equations :
⎩⎨⎧vt∂tzt∂tIt=−μρKσ⋆ (zt∇It)=−∇⋅(ztvt) =−⟨vt,∇It⟩
Advection equation
Continuity equation
Given a group of image I1,…,In and the Riemannian metric d we compute the template I by minimising:
M(I)=2n1k=1∑nd(I,I1k)2
Metamorphosis := Diffeomorphic registration + Intensity changes
[Trouvé, 2005]