Matrix-valued kernels for shape deformation analysis Mario Micheli University of San Francisco Infinite-dimensional Riemannian geometry with applications to image matching and shape analysis Erwin Schr¨ odinger International Institute for Mathematical Physics January 16, 2015 Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
38
Embed
Matrix-valued kernels for shape deformation analysisshape2015/slides/micheli.pdf · Mario Micheli University of San Francisco In nite-dimensional Riemannian geometry with applications
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Infinite-dimensional Riemannian geometry withapplications to image matching and shape analysis
Erwin Schrodinger International Institute for Mathematical Physics
January 16, 2015
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
Summary of talk:
Motivation: shape spaces, LDDMM, manifold of landmarks
Metrics induced my Matrix-valued kernels
Examples
Reference:
MM, J Glaunes, Matrix-valued Kernels for Shape Deformation
Analysis. Geometry, Imaging, and Computing , 1(1):57-139, 2014.
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
Motivation: why measure distance between shapes?
Motivation comes mostly from Computational Anatomy (CA)
E.g.: the hippocampus is deformed in a characteristic wayby Alzheimer’s disease
Idea: if shape and deformation can be describedwe can perform diagnosis from shape
Goals:? build templates,? perform classification on “shape spaces”;? more in general, do statistics on shapes
We need a distance function between shapes:(1) mathematically sound, (2) computable, and(3) relevant for the specific application one has in mind
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
Modern technology allows very accurate identification andvisualization of anatomical structures:
Landmarks
Sulcal ribbons
Diffusion tensor fields
Mathematically, “shape spaces” could be sets of:curves in R2, curves in R3, surfaces in R3, scalar images,diffusion tensor images, landmark points . . .
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
‖x‖2eig(Pr‖x) = a (multiplicity 1) and b (multiplicity d− 1)
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
General form of TRI kernels
Theorem
If ‖ · ‖V is translation- and rotation-invariant then its kernel hasthe form
K(x, y) = k(x− y)
withk(x) = k‖(‖x‖) Pr‖x + k⊥(‖x‖) Pr⊥x
for some functions k‖, k⊥ : R+ → R, scalar valued.
Remark 1: k‖(‖x‖), k⊥(‖x‖) are the eigenvalues of k(x)
Reminder: in the scalar case k(x) = k(‖x‖)Id.
Remark 2:when k‖ = k⊥ =: k we are in the scalar case (Pr
‖x+Pr⊥x =Id)
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
Positive definiteness of the kernels
Theorem (Bochner for matrix-valued functions)
A matrix-valued function k ∈ L1(Rd,Rd×d) is positive definiteiff, ∀ξ, k(ξ) := F [k](ξ) is Hermtian withnon-negative eigenvalues.
Fact: k(x) = k‖(‖x‖) Pr‖x + k⊥(‖x‖) Pr⊥x
⇔ k(ξ) = h‖(‖ξ‖) Pr‖ξ + h⊥(‖ξ‖) Pr⊥ξ
for some pair of functions h‖, h⊥ : R+ → R.
So: k(·) is positive definite iff h‖ ≥ 0 and h⊥ ≥ 0.
Question: how are the pairs (k‖, k⊥) and (h‖, h⊥) related?
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
“Generalized” Hankel Transform
The map T : (k‖, k⊥) 7→ (h‖, h⊥) is given by:
h‖(%) =2π
%µ
∫ ∞0
rµ+1 k‖(r)Jµ(2π%r) dr
− 2µ+ 1
%µ+1
∫ ∞0
rµ(k‖(r)− k⊥(r)
)Jµ+1(2π%r) dr.
h⊥(%) =2π
%µ
∫ ∞0
rµ+1 k⊥(r)Jµ(2π%r) dr
+1
%µ+1
∫ ∞0
rµ(k‖(r)− k⊥(r)
)Jµ+1(2π%r) dr.
where µ := d2 − 1
Remark:
• condition for positive definiteness of k: h‖ ≥ 0 and h⊥ ≥ 0
• when k‖ = k⊥ we retrieve the usual Hankel transform
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
Fact: the Generalized Hankel Transform T is an involution,i.e. T−1 = T
so T−1 : (h‖, h⊥) 7→ (k‖, k⊥) is given by:
k‖(r) =2π
rµ
∫ ∞0
%µ+1 h‖(%)Jµ(2πr%) d%
− 2µ+ 1
rµ+1
∫ ∞0
%µ(h‖(%)− h⊥(%)
)Jµ+1(2πr%) d%.
k⊥(r) =2π
rµ
∫ ∞0
%µ+1 h⊥(%)Jµ(2πr%) d%
+1
rµ+1
∫ ∞0
%µ(h‖(%)− h⊥(%)
)Jµ+1(2πr%) d%.
where µ :=d
2− 1
Remark:
• such formulas can be used to build functions k‖, k⊥
starting from non-negative h‖, h⊥ (generalizes [Schonberg 1938])
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
Main result: meaning of h‖ and h⊥
k(x) = k‖(‖x‖) Pr‖x + k⊥(‖x‖) Pr⊥x
k(ξ) = h‖(‖ξ‖) Pr‖ξ + h⊥(‖ξ‖) Pr⊥ξ
F
Properties of the vector fields in the RKSH V :
Condition on (h‖, h⊥)
divergence-free h‖ = 0
curl-free h⊥ = 0
Remark 1: (h‖, h⊥) correspond to the Hodge decomposition of v.f.’s.Remark 2: provides a technique to build div-free and curl free kernels.
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
Examples of vector fields
scalarkernel
divergence-freekernel
curl-freekernel
v(x) = k(x)p1, with p1 = (1, 0)
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
Examples of vector fields
scalarkernel
divergence-freekernel
curl-freekernel
v(x) =
3∑a=1
k(x− qa)pa
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
A “non-example”: k‖ and k⊥ both Gaussian
Can we have:
k‖(r) = exp(− 1
2
r2
σ21
),
k⊥(r) = exp(− 1
2
r2
σ22
)?
The corresponding kernel k(x)is positive definite if and only if
σ1 = σ2
i.e. if the kernel is scalar!
k||
k!
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
Example #1
k‖(r) = b exp(− 1
2
r2
σ2
),
k⊥(r) = (b− ar2) exp(− 1
2
r2
σ2
).
Positive definiteness of k:(a, b) ∈ D, with:
D ={a ≥ 0 & b ≥ σ2(d−1)a}
a
b
0
D b = σ2(d− 1)a
@Idiv = 0
a = 1, b = 1, σ2 = 12, d = 2
k||
k!
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
a = 0, b = 1, σ2 = 12, d = 2
a = 32, b = 1, σ2 = 1
2, d = 2
a = 12, b = 1, σ2 = 1
2, d = 2
a = 2, b = 1, σ2 = 12, d = 2
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
Example #2
k‖(r) = (b− ar2) exp(− 1
2
r2
σ2
),
k⊥(r) = b exp(− 1
2
r2
σ2
).
Positive definiteness of k:(a, b) ∈ D, with:
D ={a ≥ 0 & b ≥ σ2a}
a
b
0
Db = σ2a
@Icurl = 0
a = 1, b = 1, σ2 = 12, d = 2
k||
k!
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
a = 0, b = 1, σ2 = 12, d = 2
a = 32, b = 1, σ2 = 1
2, d = 2
a = 12, b = 1, σ2 = 1
2, d = 2
a = 2, b = 1, σ2 = 12, d = 2
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
A constructive method for curl-free and div-free kernels
Take a scalar kernel:
k(x) = k(‖x‖)Id, x ∈ Rd
with k at least twice differentiable.
1 If
k‖1(r) := −k′′(r) and k⊥1 (r) := −1
rk′(r)
then k1(x) := k‖1(‖x‖) Pr
‖x + k⊥1 (‖x‖) Pr⊥x is curl-free.
2 If
k‖2(r) := −d− 1
rk′(r) and k⊥2 (r) := −d− 2
rk′(r)−k′′(r)
then k2(x) := k‖2(‖x‖) Pr
‖x + k⊥2 (‖x‖) Pr⊥x is div-free.
Result: all curl-free and div-free kernels can be constructed this way.
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
1 landmark matching (BVP), scalar kernel
a=0, b=0.03125, c=16 (scalar)
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
1 landmark matching, div-free kernel
a=1, b=0.03125, c=16 (div free)
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
1 landmark matching, curl-free kernel
a= 1, b=0.03125, c=16 (curl free)
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
2 landmarks matching (parallel), scalar kernel
a=0, b=0.03125, c=16 (scalar)
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
2 landmarks matching (parallel), div-free kernel
a=1, b=0.03125, c=16 (div free)
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
2 landmarks matching (parallel), curl-free kernel
a= 1, b=0.03125, c=16 (curl free)
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
2 landmarks matching (opposite), scalar kernel
a=0, b=0.03125, c=16 (scalar)
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
2 landmarks matching (opposite), div-free kernel
a=1, b=0.03125, c=16 (div free)
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
2 landmarks matching (opposite), curl-free kernel
a= 1, b=0.03125, c=16 (curl free)
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
Why are we interested in div-free and curl-free kernels?
Divergence-free kernels can be used to model deformationsthat preserve volume (e.g. contractions of the heart’smuscular tissues).
Curl-free kernels are suitable for the longitudinal studies ofdeformations that are purely caused by the growth or loss ofmatter (e.g. during the development of the brain or thedevelopment of neurodegenerative diseases)
Open problem: in general, one has the interest in learningwhat kernels are best suitable for specific applications(challenging, open problem).
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis
Thank you for listening!
Mario Micheli (USF) Matrix-valued kernels for shape deformation analysis