multi-scale differential geometry and applications - casa phd-day€¦ · 13-11-2008 ·...
TRANSCRIPT
![Page 1: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/1.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Multi-scale Differential Geometry andApplicationsCASA PhD-day
Laura Astola
13 November 2008
![Page 2: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/2.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Outline
1 Introduction
2 Two ways of fitting a HOT to HARDI data
3 Finsler Geometry on a HOT field
4 Applications in the analysis of HARDI data
5 Goals in near future
![Page 3: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/3.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Outline
1 Introduction
2 Two ways of fitting a HOT to HARDI data
3 Finsler Geometry on a HOT field
4 Applications in the analysis of HARDI data
5 Goals in near future
![Page 4: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/4.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
Main application: The human brain white matter architecture.
![Page 5: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/5.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
MRI: Diagnosis of brain tumor, Multiple Sclerosis,Schizophrenia(a scalar per voxel)
DTI: Promising in diagnosis of Alzheimer, pre-operativetract detectionOne dominant fiber direction, may be found(3× 3 two-tensor per voxel)
HARDI:Several fiber directions, may be found(4.− 8.th order tensor per voxel)
![Page 6: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/6.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
MRI: Diagnosis of brain tumor, Multiple Sclerosis,Schizophrenia(a scalar per voxel)
DTI: Promising in diagnosis of Alzheimer, pre-operativetract detectionOne dominant fiber direction, may be found(3× 3 two-tensor per voxel)
HARDI:Several fiber directions, may be found(4.− 8.th order tensor per voxel)
![Page 7: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/7.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
MRI: Diagnosis of brain tumor, Multiple Sclerosis,Schizophrenia(a scalar per voxel)
DTI: Promising in diagnosis of Alzheimer, pre-operativetract detectionOne dominant fiber direction, may be found(3× 3 two-tensor per voxel)
HARDI:Several fiber directions, may be found(4.− 8.th order tensor per voxel)
![Page 8: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/8.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
MRI: Diagnosis of brain tumor, Multiple Sclerosis,Schizophrenia(a scalar per voxel)
DTI: Promising in diagnosis of Alzheimer, pre-operativetract detectionOne dominant fiber direction, may be found(3× 3 two-tensor per voxel)
HARDI:Several fiber directions, may be found(4.− 8.th order tensor per voxel)
![Page 9: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/9.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
MRI: Diagnosis of brain tumor, Multiple Sclerosis,Schizophrenia(a scalar per voxel)
DTI: Promising in diagnosis of Alzheimer, pre-operativetract detectionOne dominant fiber direction, may be found(3× 3 two-tensor per voxel)
HARDI:Several fiber directions, may be found(4.− 8.th order tensor per voxel)
![Page 10: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/10.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
MRI: Diagnosis of brain tumor, Multiple Sclerosis,Schizophrenia(a scalar per voxel)
DTI: Promising in diagnosis of Alzheimer, pre-operativetract detectionOne dominant fiber direction, may be found(3× 3 two-tensor per voxel)
HARDI:Several fiber directions, may be found(4.− 8.th order tensor per voxel)
![Page 11: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/11.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
MRI: Diagnosis of brain tumor, Multiple Sclerosis,Schizophrenia(a scalar per voxel)
DTI: Promising in diagnosis of Alzheimer, pre-operativetract detectionOne dominant fiber direction, may be found(3× 3 two-tensor per voxel)
HARDI:Several fiber directions, may be found(4.− 8.th order tensor per voxel)
![Page 12: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/12.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
To find the diffusion tensor D, in
S(y) = S0e−b·yT Dy ,
D =
(d11 d12d12 d22
), m(ϕ) = −1
b ln(SϕS0
),
we solve
(cosϕ sinϕ
) (d11 d12d12 d22
)(cosϕsinϕ
)= m(ϕ) ,
i.e. cos2(0) 2 cos(0) sin(0) sin2(0)
cos2(π4 ) 2 cos(π4 ) sin(π4 ) sin2(π4 )
cos2(π2 ) 2 cos(π2 ) sin(π2 ) sin2(π2 )
d11
d12d22
=
m(0)m(π4 )m(π2 )
![Page 13: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/13.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
To find the diffusion tensor D, in
S(y) = S0e−b·yT Dy ,
D =
(d11 d12d12 d22
), m(ϕ) = −1
b ln(SϕS0
),
we solve
(cosϕ sinϕ
) (d11 d12d12 d22
)(cosϕsinϕ
)= m(ϕ) ,
i.e. cos2(0) 2 cos(0) sin(0) sin2(0)
cos2(π4 ) 2 cos(π4 ) sin(π4 ) sin2(π4 )
cos2(π2 ) 2 cos(π2 ) sin(π2 ) sin2(π2 )
d11
d12d22
=
m(0)m(π4 )m(π2 )
![Page 14: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/14.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
To find the diffusion tensor D, in
S(y) = S0e−b·yT Dy ,
D =
(d11 d12d12 d22
), m(ϕ) = −1
b ln(SϕS0
),
we solve
(cosϕ sinϕ
) (d11 d12d12 d22
)(cosϕsinϕ
)= m(ϕ) ,
i.e. cos2(0) 2 cos(0) sin(0) sin2(0)
cos2(π4 ) 2 cos(π4 ) sin(π4 ) sin2(π4 )
cos2(π2 ) 2 cos(π2 ) sin(π2 ) sin2(π2 )
d11
d12d22
=
m(0)m(π4 )m(π2 )
![Page 15: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/15.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
To find the diffusion tensor D, in
S(y) = S0e−b·yT Dy ,
D =
(d11 d12d12 d22
), m(ϕ) = −1
b ln(SϕS0
),
we solve
(cosϕ sinϕ
) (d11 d12d12 d22
)(cosϕsinϕ
)= m(ϕ) ,
i.e. cos2(0) 2 cos(0) sin(0) sin2(0)
cos2(π4 ) 2 cos(π4 ) sin(π4 ) sin2(π4 )
cos2(π2 ) 2 cos(π2 ) sin(π2 ) sin2(π2 )
d11
d12d22
=
m(0)m(π4 )m(π2 )
![Page 16: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/16.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
More than 5 measurements ?
Solve a four-tensor dijkl , best (in L2) approximating theODF-data :
cos4 ϕ41 4 cos3 ϕ1 sinϕ1 6 cos2 ϕ1 sin2 ϕ1 4 cosϕ1 sin3 ϕ1 sin4 ϕ1
cos4 ϕ42 4 cos3 ϕ2 sinϕ2 6 cos2 ϕ2 sin2 ϕ2 4 cosϕ2 sin3 ϕ2 sin4 ϕ2
.
.
.
.
.
.cos4 ϕ4
5 4 cos3 ϕ5 sinϕ5 6 cos2 ϕ5 sin2 ϕ5 4 cosϕ5 sin3 ϕ5 sin4 ϕ5
d1111d1112d1122d1222d2222
=
ODF (ϕ1)ODF (ϕ2)ODF (ϕ3)ODF (ϕ4)ODF (ϕ5)
More than seven measurements ?
Solve a sixth order tensor dijklmn and so forth . . .
![Page 17: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/17.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
More than 5 measurements ?
Solve a four-tensor dijkl , best (in L2) approximating theODF-data :
cos4 ϕ41 4 cos3 ϕ1 sinϕ1 6 cos2 ϕ1 sin2 ϕ1 4 cosϕ1 sin3 ϕ1 sin4 ϕ1
cos4 ϕ42 4 cos3 ϕ2 sinϕ2 6 cos2 ϕ2 sin2 ϕ2 4 cosϕ2 sin3 ϕ2 sin4 ϕ2
.
.
.
.
.
.cos4 ϕ4
5 4 cos3 ϕ5 sinϕ5 6 cos2 ϕ5 sin2 ϕ5 4 cosϕ5 sin3 ϕ5 sin4 ϕ5
d1111d1112d1122d1222d2222
=
ODF (ϕ1)ODF (ϕ2)ODF (ϕ3)ODF (ϕ4)ODF (ϕ5)
More than seven measurements ?
Solve a sixth order tensor dijklmn and so forth . . .
![Page 18: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/18.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
Order of the tensor ≤ 2, Order of the tensor > 2,
![Page 19: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/19.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Introduction
A sixth order totallysymmetric tensor Dijklmn,i , j , k , l ,m,n = 1,2,3 incolors.
![Page 20: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/20.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Outline
1 Introduction
2 Two ways of fitting a HOT to HARDI data
3 Finsler Geometry on a HOT field
4 Applications in the analysis of HARDI data
5 Goals in near future
![Page 21: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/21.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Direct fitting vs. hierarchical fitting
Previous examples were direct n.th order monomial fittings
Fitting can also be done hierarchically
As functions on the sphere, these two are equivalent.
![Page 22: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/22.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Direct fitting vs. hierarchical fitting
Previous examples were direct n.th order monomial fittings
Fitting can also be done hierarchically
As functions on the sphere, these two are equivalent.
![Page 23: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/23.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Direct fitting vs. hierarchical fitting
Previous examples were direct n.th order monomial fittings
Fitting can also be done hierarchically
As functions on the sphere, these two are equivalent.
![Page 24: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/24.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Direct fitting vs. hierarchical fitting
Hierarchical D,
D(y) =∞∑
n=0
Di1...iny i1 · · · y in n ∈ 2N .
(1)y = (sin θ cosϕ, sin θ sinϕ, cos θ) = (y1, y2, y3)
Direct D,
D(y) = Di1···iny i1 · · · y in . (2)
cos(6ϕ− π7 ) + 3.4
![Page 25: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/25.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Direct fitting vs. hierarchical fitting
Hierarchical D, Laplace-Beltrami smoothing is easy !
Dτ (y) = e−τ(2+1)2Dijy iy j + e−τ(4+1)4Dijkly iy jyky l + · · ·
=∞∑
k=0
e−τk(k+1)Di1...iny i1 · · · y in (3)
![Page 26: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/26.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Direct fitting vs. hierarchical fitting
Direct D, assigning a Finslernorm F is easy !
F (y) = (Di1···iny i1 · · · y in)(1/n)
(4)
![Page 27: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/27.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Direct fitting vs. hierarchical fitting
(y1)2 + (y2)2 + (y3)2 = 1.
Example, a fourth order tensor.
Dijkly iy jyky l = D4ijklyiy jyky l + D2ijy
iy j + D0 =⇒ (5)
Dijkl = Dijkl
+∑
i,j,k ,(i≤j)
1µ(ijkk)
(µ∑
N=1
DiσN jσN kσN kσN
)|k=l
+∑i,k
1µ(iikk)
(µ∑
N=1
DiσN iσN kσN kσN
)|k=l,i=j ,
(6)
where µ(ijkl) is the number of permutations of (i , j , k , l) and σN theN.th permutation.
![Page 28: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/28.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Direct fitting vs. hierarchical fitting
(y1)2 + (y2)2 + (y3)2 = 1.
Example, a fourth order tensor.
Dijkly iy jyky l = D4ijklyiy jyky l + D2ijy
iy j + D0 =⇒ (5)
Dijkl = Dijkl
+∑
i,j,k ,(i≤j)
1µ(ijkk)
(µ∑
N=1
DiσN jσN kσN kσN
)|k=l
+∑i,k
1µ(iikk)
(µ∑
N=1
DiσN iσN kσN kσN
)|k=l,i=j ,
(6)
where µ(ijkl) is the number of permutations of (i , j , k , l) and σN theN.th permutation.
![Page 29: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/29.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Direct fitting vs. hierarchical fitting
![Page 30: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/30.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Outline
1 Introduction
2 Two ways of fitting a HOT to HARDI data
3 Finsler Geometry on a HOT field
4 Applications in the analysis of HARDI data
5 Goals in near future
![Page 31: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/31.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Finsler geometry, roughly
Vector space Rn
Inner product 〈u,v〉, in Rn
Norm ||y || in Rn
Manifold M
Riemann metricgx = 〈u, v〉x , x ∈ M
Finsler metric F (x , y),x ∈ M, y ∈ TxM
![Page 32: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/32.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Finsler geometry, roughly
Vector space Rn
Inner product 〈u,v〉, in Rn
Norm ||y || in Rn
Manifold M
Riemann metricgx = 〈u, v〉x , x ∈ M
Finsler metric F (x , y),x ∈ M, y ∈ TxM
![Page 33: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/33.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Finsler geometry, roughly
Vector space Rn
Inner product 〈u,v〉, in Rn
Norm ||y || in Rn
Manifold M
Riemann metricgx = 〈u, v〉x , x ∈ M
Finsler metric F (x , y),x ∈ M, y ∈ TxM
![Page 34: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/34.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Finsler geometry, in applications
Optics
Seismology
Mathematical ecology
Relativity theory
Medical field
![Page 35: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/35.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Finsler structure, definition
Let M be n-dimensional C∞ manifold, TxM the tangent spaceat x ∈ M, TM = {(x , y) | x ∈ M, y ∈ TxM} .
Finsler structure on M is a function F : TM → [0,∞) satisfying:
Differentiability: F is C∞ in TM.
Homogeneity: F (λy) = λF (y), ∀λ > 0.
Strong convexity: 12∂2F 2
∂y i∂y j positive definite.
![Page 36: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/36.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Finsler geometry, with a tensor-valued F
Differentiability, OK.
Homogeneity, OK.
Strong convexity, OK with a condition on unit level set.
![Page 37: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/37.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Finsler geometry, with tensor-valued F
Let g(θ, ϕ) satisfy F (g) = 1. Define matrices
m =
g1 g2 g3
g1θ g2
θ g3θ
g1ϕ g2
ϕ g3ϕ
,mθ =
g1θ g2
θ g3θ
g1θ g2
θ g3θ
g1ϕ g2
ϕ g3ϕ
,m =
g1ϕ g2
ϕ g3ϕ
g1θ g2
θ g3θ
g1ϕ g2
ϕ g3ϕ
.
Then
det(mθ)
det(m)> 0 and
det(mϕ)
det(m)>
(gij y iθy
jϕ)2
gij y iθy
jθ
(7)
implies strong convexity.
![Page 38: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/38.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Outline
1 Introduction
2 Two ways of fitting a HOT to HARDI data
3 Finsler Geometry on a HOT field
4 Applications in the analysis of HARDI data
5 Goals in near future
![Page 39: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/39.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Whenever strongconvexity is satisfied, weobtain directional metrictensors gij(y) = ∂2F 2
∂y i∂y j .
A sixth order tensor withone of its metric tensors.Blue line is principaleigenvector. Red lineparameter y .
![Page 40: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/40.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Bunch of metric tensors in a higher order tensor
![Page 41: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/41.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
”If 〈c(t),vg(x(t), y(t))〉 > αproceed,else stop”
Tracking in subthalamicnucleus-area of a rat brain, 20steps with step-size 0.2/voxelheight.
![Page 42: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/42.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Outline
1 Introduction
2 Two ways of fitting a HOT to HARDI data
3 Finsler Geometry on a HOT field
4 Applications in the analysis of HARDI data
5 Goals in near future
![Page 43: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/43.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
soon Track Finsler-fibers in HARDI data
next Compute various Finsler curvatures insome interesting data
. . . later ? Study Laplace-Beltrami smoothing vs. Ricci flow.
![Page 44: Multi-scale Differential Geometry and Applications - CASA PhD-day€¦ · 13-11-2008 · Multi-scale Differential Geometry and Applications CASA PhD-day Laura Astola 13 November](https://reader033.vdocument.in/reader033/viewer/2022042313/5edc05b7ad6a402d6666825d/html5/thumbnails/44.jpg)
/centre for analysis, scientific computing and applications
Intro Tensor fitting Finsler Geometry Applications Future work
Questions?
?