Wavelet Frames on the Sphere for Sparse Representations in High Angular Resolution Diusion Imaging. Chen Weiqiang
|
|
- Amberly Stewart
- 5 years ago
- Views:
Transcription
1 Wavelet Frames on the Sphere for Sparse Representations in High Angular Resolution Diusion Imaging Chen Weiqiang
2 Overview 1. Introduction to High Angular Resolution Diusion Imaging (HARDI). 2. Wavelets on the Sphere for L 2 sym(s 2 ). 3. Applications of Wavelets for sparse representation of HARDI Signals. 4. Possible Future Work
3 1. Introduction to HARDI Diusion Magnetic Resonance Imaging to nd directionality of neural bers via water diusion in brain tissues.
4 1. Introduction to HARDI Diusion Magnetic Resonance Imaging to nd directionality of neural bers via water diusion in brain tissues. Directionality of neural bers useful for characterization of neuro-degenerative diseases & surgical planning.
5 1. Introduction to HARDI Diusion Magnetic Resonance Imaging to nd directionality of neural bers via water diusion in brain tissues. Directionality of neural bers useful for characterization of neuro-degenerative diseases & surgical planning. Stejskal & Tanner's model: given MR signal s, reconstruct its diusion probability density function (PDF) P via ˆ P(r) = s(q) exp(iq r) dq, r R 3. R 3
6 1. Introduction to HARDI Diusion Magnetic Resonance Imaging to nd directionality of neural bers via water diusion in brain tissues. Directionality of neural bers useful for characterization of neuro-degenerative diseases & surgical planning. Stejskal & Tanner's model: given MR signal s, reconstruct its diusion probability density function (PDF) P via ˆ P(r) = s(q) exp(iq r) dq, r R 3. R 3 For computational savings, consider High Angular Resolution Diusion Imaging (HARDI), which samples on S 2 (instead of R 3 ).
7 1. Introduction to HARDI Figure: Expected diusion PDFs in (a) isotropic, (b) single ber and (c) 2 crossing bers from a brain coronal slice.
8 1. Introduction to HARDI Compute a diusion orientation distribution function (ODF) ˆ Ψ(u) = P(ru)r 2 dr, u S 2. 0 where Ψ is a 'radial-average' of PDF P with S Ψ = 1. 2
9 1. Introduction to HARDI Compute a diusion orientation distribution function (ODF) ˆ Ψ(u) = P(ru)r 2 dr, u S 2. 0 where Ψ is a 'radial-average' of PDF P with S Ψ = 1. 2 Estimate ODF Ψ from modifying HARDI signal s s(u) = log( log s(u)), where 0 < s(u) < 1, u S 2.
10 1. Introduction to HARDI Compute a diusion orientation distribution function (ODF) ˆ Ψ(u) = P(ru)r 2 dr, u S 2. 0 where Ψ is a 'radial-average' of PDF P with S Ψ = 1. 2 Estimate ODF Ψ from modifying HARDI signal s s(u) = log( log s(u)), where 0 < s(u) < 1, u S 2. Approximate a sharper ODF via Ψ(u) 1 4π π 2 R[ b s](u), u S 2, where b is the spherical Laplacian operator, and R is the Funk Radon Transform (FRT) ˆ R[s](u) := s(q) dσ(q), u S 2, v C(u) with C(u) the great circle with normal u.
11 1. Introduction to HARDI Figure: A 3D-Diusion PDF & its corresponding ODF.
12 1. Introduction to HARDI Analytical computation of ODF Ψ 2 when modied HARDI signal ~s has representation ~s(u) R c j f j (u), u S 2, j=1 where f j are eigenfunctions of R and b, so that R[ b f j ] = λ j f j, j, giving Ψ 2 (u) 1 4π π 2 R λ j c j f j (u), u S 2. j=1
13 1. Introduction to HARDI Analytical computation of ODF Ψ 2 when modied HARDI signal ~s has representation ~s(u) R c j f j (u), u S 2, j=1 where f j are eigenfunctions of R and b, so that R[ b f j ] = λ j f j, j, giving Ψ 2 (u) 1 4π π 2 R λ j c j f j (u), u S 2. j=1 Modied spherical harmonics (modied SH) Ỹj are eigenfunctions of R and b.
14 1. Introduction to HARDI Analytical computation of ODF Ψ 2 when modied HARDI signal ~s has representation ~s(u) R c j f j (u), u S 2, j=1 where f j are eigenfunctions of R and b, so that R[ b f j ] = λ j f j, j, giving Ψ 2 (u) 1 4π π 2 R λ j c j f j (u), u S 2. j=1 Modied spherical harmonics (modied SH) Ỹj are eigenfunctions of R and b. Wavelets are then designed from nite weighted linear combinations of modied SH Ỹ j 's.
15 1. Introduction to HARDI Denote the modied spherical harmonics (SH) index set Λ := {(m, l) Z 2N 0 : m l, l 2N 0 }, where m & l are the degree & order of modied SH Ỹl m (θ, φ), (m, l) Λ, θ [0, π] and φ [0, 2π),
16 1. Introduction to HARDI Denote the modied spherical harmonics (SH) index set Λ := {(m, l) Z 2N 0 : m l, l 2N 0 }, where m & l are the degree & order of modied SH Ỹl m (θ, φ), (m, l) Λ, θ [0, π] and φ [0, 2π), { Ỹl m Θ m l (cos θ) Φ 1,m (φ), if l m 0, (θ, φ) := Θ m l (cos θ) Φ 2,m (φ), if 0 < m l, where Θ m l are dened by the associated Legendre functions Pl m as Θ m l (cos θ) := ( 1) m (2l + 1)(l m)! Pl m (cos θ), 2(l + m)!
17 1. Introduction to HARDI Denote the modied spherical harmonics (SH) index set Λ := {(m, l) Z 2N 0 : m l, l 2N 0 }, where m & l are the degree & order of modied SH Ỹl m (θ, φ), (m, l) Λ, θ [0, π] and φ [0, 2π), { Ỹl m Θ m l (cos θ) Φ 1,m (φ), if l m 0, (θ, φ) := Θ m l (cos θ) Φ 2,m (φ), if 0 < m l, where Θ m l are dened by the associated Legendre functions Pl m as Θ m l (cos θ) := ( 1) m (2l + 1)(l m)! Pl m (cos θ), 2(l + m)! Φ 1,m (φ) := { π 1/2 cos(mφ), if l m < 0, (2π) 1/2 if m = 0, Φ 2,m (φ) := π 1/2 sin(mφ), 0 < m l.
18 1. Introduction to HARDI Figure: Diagram of modied SHs Ỹ m l, m l, l = 0, 2, 4.
19 2. Wavelets on the Sphere for L 2 sym(s 2 ) L 2 sym(s 2 ) is a Hilbert space with inner product given as L 2 sym(s 2 ) := { f L 2 (S 2 ) f ( u) = f (u), u S 2 a.e. }, ˆ f, g := f g dσ, f, g L 2 sym(s 2 ). S 2
20 2. Wavelets on the Sphere for L 2 sym(s 2 ) L 2 sym(s 2 ) is a Hilbert space with inner product given as L 2 sym(s 2 ) := { f L 2 (S 2 ) f ( u) = f (u), u S 2 a.e. }, ˆ f, g := f g dσ, f, g L 2 sym(s 2 ). S 2 Modied SH {Ỹ m l is an orthonormal basis for }(m,l) Λ L 2 sym(s 2 ) with L 2 sym(s 2 ) =H 1 H 2, H i := span {Ỹ m l : (m, l) Λ i }, Λ 1 := {(m, l) Λ : l m 0}, Λ 2 := {(m, l) Λ : 1 m l}.
21 2. Wavelet Frames on the Sphere for L 2 sym(s 2 ) A countable collection {f ν } ν I in a Hilbert space H is a frame for H if constants C & D s.t. C f 2 H ν I f, f ν H 2 D f 2 H, f H. When C = D = 1, {f ν } ν I is a tight frame for H.
22 2. Wavelet Frames on the Sphere for L 2 sym(s 2 ) A countable collection {f ν } ν I in a Hilbert space H is a frame for H if constants C & D s.t. C f 2 H ν I f, f ν H 2 D f 2 H, f H. When C = D = 1, {f ν } ν I is a tight frame for H. Frames give overcomplete, yet stable representations & superior denoising compared to bases.
23 2. Wavelet Frames on the Sphere for L 2 sym(s 2 ) A countable collection {f ν } ν I in a Hilbert space H is a frame for H if constants C & D s.t. C f 2 H ν I f, f ν H 2 D f 2 H, f H. When C = D = 1, {f ν } ν I is a tight frame for H. Frames give overcomplete, yet stable representations & superior denoising compared to bases. { To build } tight (wavelet) frames for L 2 sym(s 2 ), construct f 1 ν & { f 2 } ν I 1 ν forming tight frames for H 1 & H 2. ν I 2
24 2. Wavelet Frames on the Sphere for L 2 sym(s 2 ) A countable collection {f ν } ν I in a Hilbert space H is a frame for H if constants C & D s.t. C f 2 H ν I f, f ν H 2 D f 2 H, f H. When C = D = 1, {f ν } ν I is a tight frame for H. Frames give overcomplete, yet stable representations & superior denoising compared to bases. { To build } tight (wavelet) frames for L 2 sym(s 2 ), construct f 1 ν & { f 2 } ν I 1 ν forming tight frames for H 1 & H 2. ν I 2 Construct tight (wavelet) frames for H i from its orthonormal basis {Ỹ m l : (m, l) Λ i }, i = 1, 2.
25 2. Wavelets on the Sphere for L 2 sym(s 2 ) Denote Ỹ j := Ỹ m l, where j := (m, l) Λ. For i = 1, 2, k 1, wavelets in H i take the form ψk i ( ; ηi k,ν ) = µ k,ν Ỹ j (ηk,ν i )Ỹj, ν J k, j J i k where J i k Λi, η i k,ν S2, µ k,ν R\ {0},ν J k & J k is a nite set s.t. J i k J k.
26 2. Wavelets on the Sphere for L 2 sym(s 2 ) Denote Ỹ j := Ỹ m l, where j := (m, l) Λ. For i = 1, 2, k 1, wavelets in H i take the form ψk i ( ; ηi k,ν ) = µ k,ν Ỹ j (ηk,ν i )Ỹj, ν J k, j J i k where J i k Λi, η i k,ν S2, µ k,ν R\ {0},ν J k & J k is a nite set s.t. J i k J k. Regard ψ i k ( ; ηi k,ν ) as generalized translates of ψ i k = j J i k Ỹ j under T η i k,ν, i.e., ψ i k ( ; ηi k,ν ) = T η i k,ν ψi k.
27 Motivation of Generalized Translation Operator For instance, take H := L 2 [0, 2π), S := [0, 2π), Λ := Z, e j (x) := exp(ijx), x S, j Z.
28 Motivation of Generalized Translation Operator For instance, take H := L 2 [0, 2π), S := [0, 2π), Λ := Z, e j (x) := exp(ijx), x S, j Z. For ζ S, set f (x) := N 1 j=0 α je j (x), then we have T ζ f (x) = = N 1 j=0 N 1 j=0 α j e j (ζ)e j (x) α j exp(ij(x ζ)).
29 Motivation of Generalized Translation Operator For instance, take H := L 2 [0, 2π), S := [0, 2π), Λ := Z, e j (x) := exp(ijx), x S, j Z. For ζ S, set f (x) := N 1 j=0 α je j (x), then we have T ζ f (x) = = N 1 j=0 N 1 j=0 α j e j (ζ)e j (x) α j exp(ij(x ζ)). If we choose ζ := ζ l := 2πl N, l = 0,..., N 1, then, A := [e j (ζ l )] N 1 j,l=0 = [exp(i2πjl/n)]n 1 j,l=0 satises ADA = I with D = 1 N I.
30 2. Wavelets on the Sphere for L 2 sym(s 2 ) Key Construction Idea For i = 1, 2, k N 0, dene I i k := { (m, l) Λ i : 0 l L k }, Lk 2N 0.
31 2. Wavelets on the Sphere for L 2 sym(s 2 ) Key Construction Idea For i = 1, 2, k N 0, dene I i k := { (m, l) Λ i : 0 l L k }, Lk 2N 0. Seek 1. nodes ζ i k,ν S2, ν I k with Ii k I k <, 2. weight matrices D k := diag(d k,ν ) ν I k 0, such that A i k := [Ỹj(ζ k,ν i )] j Ik i,ν I k orthogonality condition satises a discrete weighted A i k D ka i k = I. (1)
32 2. Wavelets on the Sphere for L 2 sym(s 2 ) Key Construction Idea For i = 1, 2, k N 0, dene I i k := { (m, l) Λ i : 0 l L k }, Lk 2N 0. Seek 1. nodes ζ i k,ν S2, ν I k with Ii k I k <, 2. weight matrices D k := diag(d k,ν ) ν I k 0, such that A i k := [Ỹj(ζ k,ν i )] j Ik i,ν I k orthogonality condition satises a discrete weighted A i k D ka i k = I. (1) Dene Fk i := Ai k D1/2 k, so by (1), Fk i F k i = I.
33 2. Wavelets on the Sphere for L 2 sym(s 2 ) Set [ϕ i k ( ; ζi k,ν )] ν I k := F k [Ỹ j ] j I i k, i.e., ϕ i k ( ; ζi k,ν ) = d k,ν Ỹ j (ζk,ν i ) Ỹj, ν I k. j I i k
34 2. Wavelets on the Sphere for L 2 sym(s 2 ) Set [ϕ i k ( ; ζi k,ν )] ν I k := F k [Ỹ j ] j I i k, i.e., Then, ϕ i k ( ; ζi k,ν ) = d k,ν Ỹ j (ζk,ν i ) Ỹj, ν I k. j I i k { } ϕ i k ( ; ζi k,ν ) : ν I k forms a tight frame for V i k := span {Ỹj : j I i k }, k N 0, i = 1, 2. (2)
35 2. Wavelets on the Sphere for L 2 sym(s 2 ) Set [ϕ i k ( ; ζi k,ν )] ν I k := F k [Ỹ j ] j I i k, i.e., Then, ϕ i k ( ; ζi k,ν ) = d k,ν Ỹ j (ζk,ν i ) Ỹj, ν I k. j I i k { } ϕ i k ( ; ζi k,ν ) : ν I k forms a tight frame for V i k := span {Ỹj : j I i k }, k N 0, i = 1, 2. (2) Each ϕ i k ( ; ζi k,ν ) is localized around ζi k,ν, ν I k, i.e., ϕ i k ( ; ζi k,ν ) ϕ i k (ζi k,ν ; ζi k,ν ) = min { f : f V i k, f (ζi k,ν ) = 1}.
36 2. Wavelets on the Sphere for L 2 sym(s 2 ) Lemma: { } For i = 1, 2, k N 0, dene modied Gauss-Legendre quadrature nodes ζk,ν i : ν I k ζ i k,ν := (θ k,α, φ i k,β ), ν := (α, β) I k, (3) where θ k,α := cos 1 (x k,α ), α = 0, 1,..., L k, with {x α} L k α=0 as zeros of the Legendre polynomial P Lk +1 & the zenith points are φ 1 k,β π(β + 1/2) =, φ 2 π(β + 1) k,β = L k + 1 L k + 2, β = 0, 1,..., L k.
37 2. Wavelets on the Sphere for L 2 sym(s 2 ) Lemma: { } For i = 1, 2, k N 0, dene modied Gauss-Legendre quadrature nodes ζk,ν i : ν I k ζ i k,ν := (θ k,α, φ i k,β ), ν := (α, β) I k, (3) where θ k,α := cos 1 (x k,α ), α = 0, 1,..., L k, with {x α} L k α=0 as zeros of the Legendre polynomial P Lk +1 & the zenith points are φ 1 k,β π(β + 1/2) =, φ 2 π(β + 1) k,β = L k + 1 L k + 2, β = 0, 1,..., L k. Enumerate the nodes as { } {(θ k,α, φ i k,β }(α,β) ) := (θ k,0, φ i k,0 ),..., (θ k,0, φ i k,l k ),..., (θ k,lk, φ i k,0 ),..., (θ k,l k, φ i k,l k ), & dene the diagonal weight matrix D k = 2π L k + 1 diag(w k,0,..., w k,0,..., w k,lk,..., w k,lk ), (4) where { w k,α } Lk α=0 are the Gauss-quadrature weights of P L k +1, with each unique w k,α appearing L k + 1 times in D k.
38 2. Wavelets on the Sphere for L 2 sym(s 2 ) Then, for i = 1, 2, k N 0, A i k := [Ỹ j (ζk,ν i )] j Ik i,ν I k satises the weighted orthogonality condition (1), with diag(d k,ν ) ν I k := D k. Furthermore, for i = 1, 2, k N 0, { ϕ i k ( ; ζi k,ν ) : ν I k } yields a tight frame for V i k dened in (2), where ϕ i k ( ; ζi k,ν ) := d k,ν Ỹ j (ζk,ν i )Ỹ j, ν I k. (5) j I i k
39 2. Wavelets on the Sphere for L 2 sym(s 2 ) A sequence of nite-dimensional subspaces {V k } k 0 in H is a multiresolution analysis (MRA) of H if 1. k 0, functions ϕ k ( ; ζ k,ν ), ν I k, s.t. {ϕ k ( ; ζ k,ν ) : ν I k } spans V k. 2. k 0, V k V k k 0 V k = H. We call the functions ϕ k ( ; ζ k,ν ), ν I k, k 0, scaling functions.
40 2. Wavelets on the Sphere for L 2 sym(s 2 ) A sequence of nite-dimensional subspaces {V k } k 0 in H is a multiresolution analysis (MRA) of H if 1. k 0, functions ϕ k ( ; ζ k,ν ), ν I k, s.t. {ϕ k ( ; ζ k,ν ) : ν I k } spans V k. 2. k 0, V k V k k 0 V k = H. We call the functions ϕ k ( ; ζ k,ν ), ν I k, k 0, scaling functions. Proposition: For i = 1, 2, k N 0, let Vk i in (2) & be the subspace dened Ik i := { (m, l) Λ i } : 0 l L k, Lk 2N 0, where {L k } k N0. Then, for i = 1, 2, { Vk i } forms an k N 0 MRA for H i with scaling functions ϕ i k ( ; ζi k,ν ), ν I k, given by (5).
41 2. Wavelets on the Sphere for L 2 sym(s 2 ) Figure: Scaling functions of V 1 0,V 2 0, where V i 0 = { (m, l) Λ i : 0 l 4 }, i = 1, 2.
42 2. Wavelets on the Sphere for L 2 sym(s 2 ) Figure: Scaling functions of V 1 1, V 2 1,where V i 1 = { (m, l) Λ i : 0 l 6 }, i = 1, 2.
43 2. Wavelets on the Sphere for L 2 sym(s 2 ) For i = 1, 2, k 0, let W i k := V i k+1 \V i k & J i k := Ii k+1 \Ii k. Then where W i k := span {Ỹj : j J i k }, (6) J i k := { (m, l) Λ i : L k < l L k+1 }. (7)
44 2. Wavelets on the Sphere for L 2 sym(s 2 ) For i = 1, 2, k 0, let W i k := V i k+1 \V i k & J i k := Ii k+1 \Ii k. Then where W i k := span {Ỹj : j J i k }, (6) J i k := { (m, l) Λ i : L k < l L k+1 }. (7) Thus, we have an orthogonal decomposition H i = W 1 i W0 i W1 i..., where W 1 i := V 0 i, i = 1, 2.
45 2. Wavelets on the Sphere for L 2 sym(s 2 ) Theorem: For i = 1, 2, k 1, dene ψ i k ( ; ηi k,ν ) := ˆd k,ν Ỹl m (ηk,ν i )Ỹ l m, ηk,ν i J k (m,l) J i k where J i k is given by (7), ηi k,ν := ζi k+1,ν, ν J k := I k+1 & diag( ˆd k,ν ) := ˆD k := D k+1 given by (3) & (4).
46 2. Wavelets on the Sphere for L 2 sym(s 2 ) Theorem: For i = 1, 2, k 1, dene ψ i k ( ; ηi k,ν ) := ˆd k,ν Ỹl m (ηk,ν i )Ỹ l m, ηk,ν i J k (m,l) J i k where J i k is given by (7), ηi k,ν := ζi k+1,ν, ν J k := I k+1 & diag( ˆd k,ν ) := ˆD k := D k+1 given by (3) & (4). Denote W i 1 := V i 0, J i 1 := Ii 0, ψi 1 ( ; ηi 1,ν ) := ϕi 0 ( ; ζi 0,ν ), i = 1, 2, and Γ(Ψ) := { ψ i k ( ; ηi k,ν ) : ν J k, k 1, i = 1, 2}. Then Γ(Ψ) forms a normalized tight wavelet frame for L 2 sym(s 2 ).
47 2. Wavelets on the Sphere for L 2 sym(s 2 ) Figure: Wavelet functions of W 1 0, W 2 0, where W i 0 = { (m, l) Λ i : l = 6 }, i = 1, 2.
48 3. Applications of Wavelets to HARDI Signals Key idea: ODF Ψ 2 is localized along a few radial directions, so represent ODF Ψ 2 with localized wavelets in Γ(Ψ).
49 3. Applications of Wavelets to HARDI Signals Key idea: ODF Ψ 2 is localized along a few radial directions, so represent ODF Ψ 2 with localized wavelets in Γ(Ψ). However, only HARDI signal s is known, so seek another wavelet system Γ( Ψ) to rst represent it.
50 3. Applications of Wavelets to HARDI Signals Key idea: ODF Ψ 2 is localized along a few radial directions, so represent ODF Ψ 2 with localized wavelets in Γ(Ψ). However, only HARDI signal s is known, so seek another wavelet system Γ( Ψ) to rst represent it. To reconstruct unknown ODF Ψ 2, given modied HARDI signal s(u), apply R and b as follows: Ψ 2 (u) 1 4π π 2 R[ b s](u), u S 2.
51 3. Applications of Wavelets to HARDI Signals Key idea: ODF Ψ 2 is localized along a few radial directions, so represent ODF Ψ 2 with localized wavelets in Γ(Ψ). However, only HARDI signal s is known, so seek another wavelet system Γ( Ψ) to rst represent it. To reconstruct unknown ODF Ψ 2, given modied HARDI signal s(u), apply R and b as follows: Ψ 2 (u) 1 4π π 2 R[ b s](u), u S 2. { } Thus, nd Γ( Ψ) := ψ k i ( ; ηi k,ν ) : ν J k, k 1, i = 1, 2 to represent ~s, where ψ i k ( ; ηi k,ν ) is a `pre-image' of ψi k ( ; ηi k,ν ) in Γ(Ψ) under b and R.
52 3. Applications of Wavelets to HARDI Signals To nd Γ( Ψ), recall ψk i ( ; ηi k,ν ) are nite combinations of modied SH {Ỹl m } (m,l) Λ, i.e., ψ i k ( ; ηi k,ν ) := ˆd k,ν (m,l) J i k Ỹl m (ηk,ν i )Ỹ l m, ηk,ν i J k.
53 3. Applications of Wavelets to HARDI Signals To nd Γ( Ψ), recall ψk i ( ; ηi k,ν ) are nite combinations of modied SH {Ỹl m } (m,l) Λ, i.e., ψ i k ( ; ηi k,ν ) := ˆd k,ν (m,l) J i k Ỹl m (ηk,ν i )Ỹ l m, ηk,ν i J k. {Ỹ l m } (m,l) Λ are eigenfunctions of b and R, i.e., b Ỹ m l = l(l + 1)Ỹl m, R[Ỹl m ] = 2πP l (0)Ỹ m l, (m, l) Λ.
54 3. Applications of Wavelets to HARDI Signals To nd Γ( Ψ), recall ψk i ( ; ηi k,ν ) are nite combinations of modied SH {Ỹl m } (m,l) Λ, i.e., ψ i k ( ; ηi k,ν ) := ˆd k,ν (m,l) J i k Ỹl m (ηk,ν i )Ỹ l m, ηk,ν i J k. {Ỹ l m } (m,l) Λ are eigenfunctions of b and R, i.e., b Ỹ m l = l(l + 1)Ỹl m, R[Ỹl m ] = 2πP l (0)Ỹ m l, (m, l) Λ. To obtain ψ i k ( ; ηi k,ν ) as a 'pre-image' of ψi k ( ; ηi k,ν ) under b and R, for i = 1, 2, k 1, ν J k, dene ψ i k ( ; ηi k,ν ) : = ˆd k,ν (m,l) J i k δ 1 l where δ l = l(l + 1), l 2N, and δ 0 = 1. [2πP l (0)] 1 Ỹl m (ηk,ν i )Ỹ l m,
55 3. Applications of Wavelets to HARDI Signals Then, approximate modied HARDI signal s(u) R c j ψ j (u). j=1
56 3. Applications of Wavelets to HARDI Signals Then, approximate modied HARDI signal s(u) R c j ψ j (u). j=1 ODF Ψ 2 can be reconstructed analytically: Ψ 2 (u) 1 4π π 2 = 1 4π π 2 R c j R[ b ψj ](u) j=1 R c j ψ j (u). j=1
57 3. Applications of Wavelets to HARDI Signals Then, approximate modied HARDI signal s(u) R c j ψ j (u). j=1 ODF Ψ 2 can be reconstructed analytically: Ψ 2 (u) 1 4π π 2 = 1 4π π 2 R c j R[ b ψj ](u) j=1 R c j ψ j (u). j=1 In practice, ~s is discretized with gradient directions {g k } K k=1, giving f = [~s(g k )] K k=1 & discrete wavelet matrix B where B j,k = ψ j (g k ).
58 3. Applications of Wavelets to HARDI Signals For sparse reconstruction of a single discrete ODF, 1 min c R R 2 f Bc 2 F + λ c 1, where c := [c j ] R j=1 is the coecient vector.
59 3. Applications of Wavelets to HARDI Signals For sparse reconstruction of a single discrete ODF, 1 min c R R 2 f Bc 2 F + λ c 1, where c := [c j ] R j=1 is the coecient vector. Use the fast iterative soft-thresholding algorithm (FISTA) to obtain an approximate solution c.
60 3. Applications of Wavelets to HARDI Signals For sparse reconstruction of a single discrete ODF, 1 min c R R 2 f Bc 2 F + λ c 1, where c := [c j ] R j=1 is the coecient vector. Use the fast iterative soft-thresholding algorithm (FISTA) to obtain an approximate solution c. Reconstruct a single discrete ODF Ψ 2,d by Ψ 2,d 1 4π π 2 BLPc, where L := diag[ l j (l j + 1)] R j=1, P := diag[2πp l j (0)] R j=1.
61 3. Applications of Wavelets to HARDI Signals Figure: A 2D image of HARDI ODFs.
62 3. Applications of Wavelets to HARDI Signals Neigbouring HARDI ODFs tend to be similar in diusion directionality.
63 3. Applications of Wavelets to HARDI Signals Neigbouring HARDI ODFs tend to be similar in diusion directionality. Improve HARDI reconstruction by modifying optimization models with HARDI spatial regularization.
64 3. Applications of Wavelets to HARDI Signals Neigbouring HARDI ODFs tend to be similar in diusion directionality. Improve HARDI reconstruction by modifying optimization models with HARDI spatial regularization. For each sampling gradient in {g k } K k=1, we have an N x N y Nx,Ny image of HARDI signals S(g k ):=[S(g k ) i,j ] i,j=1.
65 3. Applications of Wavelets to HARDI Signals Neigbouring HARDI ODFs tend to be similar in diusion directionality. Improve HARDI reconstruction by modifying optimization models with HARDI spatial regularization. For each sampling gradient in {g k } K k=1, we have an N x N y Nx,Ny image of HARDI signals S(g k ):=[S(g k ) i,j ] i,j=1. K Nx Ny From there, extract an R HARDI data matrix F and impose spatial regularization by adding diag(µ)wc 1 : min C 1 2 F BC 2 F + λ C 1 + diag(µ)wc 1, 4Nx Ny 1 where C : coe matrix, λ R +, µ R 0 are parameters that balance coecient sparsity & spatial regularization by the 2D-tensor Haar tight framelet lter matrix W.
66 3. Applications of Wavelets to HARDI Signals Solved with split Bregman method by setting Q = WC and P = C, giving 1 min F C,P,Q 2 BC 2 F +λ P 1 + diag(µ)q 1 s.t. Q = WC & P = C.
67 3. Applications of Wavelets to HARDI Signals Solved with split Bregman method by setting Q = WC and P = C, giving 1 min F C,P,Q 2 BC 2 F +λ P 1 + diag(µ)q 1 s.t. Q = WC & P = C. Denoting Lagrange multipliers as Λ 1, Λ 2 & ρ 1, ρ 2 as positive parameters, dene an augmented Lagrangian L ρ1,ρ 2 (C, Q, P; Λ 1, Λ 2 ) := 1 2 F BC 2 F + λ P 1 + µ diag(µ)q 1 ρ 1 Λ 1, Q WC + ρ 1 2 Q WC 2 F ρ 2 Λ 2, P C + ρ 2 2 P C 2 F.
68 3. Applications of Wavelets to HARDI Signals Then, for a xed (C k, Q k, P k, Λ k 1, Λk 2 ), we update (C k+1, Q k+1, P k+1, Λ k+1 1, Λ2 k+1 ) by a Gauss-Seidel scheme on the above augmented Lagrangian as follows: 1. C k+1 = arg min C L ρ1,ρ 2 (C, Q k, P k ; Λ k 1, Λk 2 ); 2. Q k+1 = arg min Q L ρ1,ρ 2 (C k+1, Q, P k ; Λ k 1, Λk 2 ); 3. P k+1 = arg min P L ρ1,ρ 2 (C k+1, Q k+1, P; Λ k 1, Λk 2 ); 4. Λ k+1 1 = Λ k 1 (Qk+1 W (C k+1 ) ); 5. Λ k+1 2 = Λ k 2 (Pk+1 C k+1 ).
69 3. Applications of Wavelets to HARDI Signals Spatially regularized wavelet frame method for CSA QBI Input: HARDI data matrix F and tolerance level ɛ. Initialization: set C 0 = [B B + λi ] 1 B F, Q 0 = W (C 0 ), P 0 = C 0, Λ 1, Λ 2 = 0. For k = 0, 1,..., while C k+1 C k > ɛ, perform: 1. C k+1 = [B B +(ρ 1 +ρ 2 )I ] 1 [B F +ρ 1 (Q k Λ k 1 ) W +ρ 2 (P k Λ 2 )]; 2. Q k+1 = T 1 µ/ρ 1 [W (C k+1 ) + Λ k 1 ]; 3. P k+1 = T 1 λ/ρ 2 [C k+1 + Λ k 2 ]; 4. Λ k+1 1 = Λ k 1 (Qk+1 W (C k+1 ) ); 5. Λ k+1 2 = Λ k 2 (Pk+1 C k+1 ). Output: Coecient matrix C := C K to reconstruct multi-voxel version of discrete ODF Ψ 2,d.
70 3. Applications of Wavelets to HARDI Signals Figure: Structured Field Testing Dataset: voxel set Ω d : Data set simulates a realistic 3D conguration of tracts: it comprises 5 dierent ber bundles (Sub-plots A-E). In each voxel, the directions are color-coded based on their orientation (x-axis, y -axis, z-axis). Subplot F shows a representative slice with the ODF orientations.
71 3. Applications of Wavelets to HARDI Signals Method / SNR (db) Spherical Harmonics (SH) Spherical Ridgelets (SR) Wavelet Frames (WF) Table: NMSE values of SH, SR, WF based methods (with spatial regularization) for HARDI testing data for SNR levels 5, 10, 20 db. Normalized mean square error (NMSE): 1 2 Ψ est(r) Ψ true(r) 2 NMSE :=. N xn y N z Ψ r Ω true(r) 2 2 d
72 Possible Future Work Extend the above wavelet frame representation to handle HARDI deconvolution problems. Design other wavelet frame systems for other cases of Hilbert spaces for practical applications. Consider integration of above techniques for feature extraction with deep learning.
Fast and Accurate HARDI and its Application to Neurological Diagnosis
Fast and Accurate HARDI and its Application to Neurological Diagnosis Dr. Oleg Michailovich Department of Electrical and Computer Engineering University of Waterloo June 21, 2011 Outline 1 Diffusion imaging
More informationOrdinary Least Squares and its applications
Ordinary Least Squares and its applications Dr. Mauro Zucchelli University Of Verona December 5, 2016 Dr. Mauro Zucchelli Ordinary Least Squares and its applications December 5, 2016 1 / 48 Contents 1
More informationMultiresolution Analysis
Multiresolution Analysis DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Frames Short-time Fourier transform
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationChapter 7 Wavelets and Multiresolution Processing
Chapter 7 Wavelets and Multiresolution Processing Background Multiresolution Expansions Wavelet Transforms in One Dimension Wavelet Transforms in Two Dimensions Image Pyramids Subband Coding The Haar
More informationSpectral transforms. Contents. ARPEGE-Climat Version 5.1. September Introduction 2
Spectral transforms ARPEGE-Climat Version 5.1 September 2008 Contents 1 Introduction 2 2 Spectral representation 2 2.1 Spherical harmonics....................... 2 2.2 Collocation grid.........................
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More informationA Riemannian Framework for Denoising Diffusion Tensor Images
A Riemannian Framework for Denoising Diffusion Tensor Images Manasi Datar No Institute Given Abstract. Diffusion Tensor Imaging (DTI) is a relatively new imaging modality that has been extensively used
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More informationThe University of Texas at Austin Department of Electrical and Computer Engineering. EE381V: Large Scale Learning Spring 2013.
The University of Texas at Austin Department of Electrical and Computer Engineering EE381V: Large Scale Learning Spring 2013 Assignment 1 Caramanis/Sanghavi Due: Thursday, Feb. 7, 2013. (Problems 1 and
More informationLinear Algebra and Dirac Notation, Pt. 3
Linear Algebra and Dirac Notation, Pt. 3 PHYS 500 - Southern Illinois University February 1, 2017 PHYS 500 - Southern Illinois University Linear Algebra and Dirac Notation, Pt. 3 February 1, 2017 1 / 16
More informationCompressed Sensing and Neural Networks
and Jan Vybíral (Charles University & Czech Technical University Prague, Czech Republic) NOMAD Summer Berlin, September 25-29, 2017 1 / 31 Outline Lasso & Introduction Notation Training the network Applications
More informationVariational Image Restoration
Variational Image Restoration Yuling Jiao yljiaostatistics@znufe.edu.cn School of and Statistics and Mathematics ZNUFE Dec 30, 2014 Outline 1 1 Classical Variational Restoration Models and Algorithms 1.1
More informationGroup Representations
Group Representations Alex Alemi November 5, 2012 Group Theory You ve been using it this whole time. Things I hope to cover And Introduction to Groups Representation theory Crystallagraphic Groups Continuous
More informationIntroduction to Discrete-Time Wavelet Transform
Introduction to Discrete-Time Wavelet Transform Selin Aviyente Department of Electrical and Computer Engineering Michigan State University February 9, 2010 Definition of a Wavelet A wave is usually defined
More informationOle Christensen 3. October 20, Abstract. We point out some connections between the existing theories for
Frames and pseudo-inverses. Ole Christensen 3 October 20, 1994 Abstract We point out some connections between the existing theories for frames and pseudo-inverses. In particular, using the pseudo-inverse
More informationAnouncements. Assignment 3 has been posted!
Anouncements Assignment 3 has been posted! FFTs in Graphics and Vision Correlation of Spherical Functions Outline Math Review Spherical Correlation Review Dimensionality: Given a complex n-dimensional
More informationOrientation Distribution Function for Diffusion MRI
Orientation Distribution Function for Diffusion MRI Evgeniya Balmashnova 28 October 2009 Diffusion Tensor Imaging Diffusion MRI Diffusion MRI P(r, t) = 1 (4πDt) 3/2 e 1 4t r 2 D 1 t Diffusion time D Diffusion
More informationEstimation of Non-Negative ODFs using the Eigenvalue Distribution of Spherical Functions
Estimation of Non-Negative ODFs using the Eigenvalue Distribution of Spherical Functions Evan Schwab, Bijan Afsari, and René Vidal Center for Imaging Science, Johns Hopkins University Abstract. Current
More informationWavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing
Wavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing Qingtang Jiang Abstract This paper is about the construction of univariate wavelet bi-frames with each framelet being symmetric.
More informationAnisotropy of HARDI Diffusion Profiles Based on the L 2 -Norm
Anisotropy of HARDI Diffusion Profiles Based on the L 2 -Norm Philipp Landgraf 1, Dorit Merhof 1, Mirco Richter 1 1 Institute of Computer Science, Visual Computing Group, University of Konstanz philipp.landgraf@uni-konstanz.de
More informationInexact Alternating Direction Method of Multipliers for Separable Convex Optimization
Inexact Alternating Direction Method of Multipliers for Separable Convex Optimization Hongchao Zhang hozhang@math.lsu.edu Department of Mathematics Center for Computation and Technology Louisiana State
More informationA Tutorial on Wavelets and their Applications. Martin J. Mohlenkamp
A Tutorial on Wavelets and their Applications Martin J. Mohlenkamp University of Colorado at Boulder Department of Applied Mathematics mjm@colorado.edu This tutorial is designed for people with little
More information2D Wavelets. Hints on advanced Concepts
2D Wavelets Hints on advanced Concepts 1 Advanced concepts Wavelet packets Laplacian pyramid Overcomplete bases Discrete wavelet frames (DWF) Algorithme à trous Discrete dyadic wavelet frames (DDWF) Overview
More informationDISCRETE CDF 9/7 WAVELET TRANSFORM FOR FINITE-LENGTH SIGNALS
DISCRETE CDF 9/7 WAVELET TRANSFORM FOR FINITE-LENGTH SIGNALS D. Černá, V. Finěk Department of Mathematics and Didactics of Mathematics, Technical University in Liberec Abstract Wavelets and a discrete
More informationIMAGE RESTORATION: TOTAL VARIATION, WAVELET FRAMES, AND BEYOND
IMAGE RESTORATION: TOTAL VARIATION, WAVELET FRAMES, AND BEYOND JIAN-FENG CAI, BIN DONG, STANLEY OSHER, AND ZUOWEI SHEN Abstract. The variational techniques (e.g., the total variation based method []) are
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More information2 Tikhonov Regularization and ERM
Introduction Here we discusses how a class of regularization methods originally designed to solve ill-posed inverse problems give rise to regularized learning algorithms. These algorithms are kernel methods
More informationNew Coherence and RIP Analysis for Weak. Orthogonal Matching Pursuit
New Coherence and RIP Analysis for Wea 1 Orthogonal Matching Pursuit Mingrui Yang, Member, IEEE, and Fran de Hoog arxiv:1405.3354v1 [cs.it] 14 May 2014 Abstract In this paper we define a new coherence
More informationSpherical Harmonics on S 2
7 August 00 1 Spherical Harmonics on 1 The Laplace-Beltrami Operator In what follows, we describe points on using the parametrization x = cos ϕ sin θ, y = sin ϕ sin θ, z = cos θ, where θ is the colatitude
More informationChapter 5 Fast Multipole Methods
Computational Electromagnetics; Chapter 1 1 Chapter 5 Fast Multipole Methods 5.1 Near-field and far-field expansions Like the panel clustering, the Fast Multipole Method (FMM) is a technique for the fast
More informationMultiscale Frame-based Kernels for Image Registration
Multiscale Frame-based Kernels for Image Registration Ming Zhen, Tan National University of Singapore 22 July, 16 Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image
More informationDeep Learning: Approximation of Functions by Composition
Deep Learning: Approximation of Functions by Composition Zuowei Shen Department of Mathematics National University of Singapore Outline 1 A brief introduction of approximation theory 2 Deep learning: approximation
More informationMultiresolution image processing
Multiresolution image processing Laplacian pyramids Some applications of Laplacian pyramids Discrete Wavelet Transform (DWT) Wavelet theory Wavelet image compression Bernd Girod: EE368 Digital Image Processing
More informationRecovering overcomplete sparse representations from structured sensing
Recovering overcomplete sparse representations from structured sensing Deanna Needell Claremont McKenna College Feb. 2015 Support: Alfred P. Sloan Foundation and NSF CAREER #1348721. Joint work with Felix
More informationLinear Algebra, 4th day, Thursday 7/1/04 REU Info:
Linear Algebra, 4th day, Thursday 7/1/04 REU 004. Info http//people.cs.uchicago.edu/laci/reu04. Instructor Laszlo Babai Scribe Nick Gurski 1 Linear maps We shall study the notion of maps between vector
More information2 Regularized Image Reconstruction for Compressive Imaging and Beyond
EE 367 / CS 448I Computational Imaging and Display Notes: Compressive Imaging and Regularized Image Reconstruction (lecture ) Gordon Wetzstein gordon.wetzstein@stanford.edu This document serves as a supplement
More informationOctober 7, :8 WSPC/WS-IJWMIP paper. Polynomial functions are renable
International Journal of Wavelets, Multiresolution and Information Processing c World Scientic Publishing Company Polynomial functions are renable Henning Thielemann Institut für Informatik Martin-Luther-Universität
More informationPhysics 70007, Fall 2009 Answers to Final Exam
Physics 70007, Fall 009 Answers to Final Exam December 17, 009 1. Quantum mechanical pictures a Demonstrate that if the commutation relation [A, B] ic is valid in any of the three Schrodinger, Heisenberg,
More informationWavelets and multiresolution representations. Time meets frequency
Wavelets and multiresolution representations Time meets frequency Time-Frequency resolution Depends on the time-frequency spread of the wavelet atoms Assuming that ψ is centred in t=0 Signal domain + t
More informationSome new method of approximation and estimation on sphere
Some new method of approximation and estimation on sphere N. Jarzebkowska B. Cmiel K. Dziedziul September 17-23, 2017, in Bedlewo (Poland) [BD] Bownik, M.; D.K. Smooth orthogonal projections on sphere.
More informationELEMENTARY LINEAR ALGEBRA WITH APPLICATIONS. 1. Linear Equations and Matrices
ELEMENTARY LINEAR ALGEBRA WITH APPLICATIONS KOLMAN & HILL NOTES BY OTTO MUTZBAUER 11 Systems of Linear Equations 1 Linear Equations and Matrices Numbers in our context are either real numbers or complex
More informationHere is an example of a block diagonal matrix with Jordan Blocks on the diagonal: J
Class Notes 4: THE SPECTRAL RADIUS, NORM CONVERGENCE AND SOR. Math 639d Due Date: Feb. 7 (updated: February 5, 2018) In the first part of this week s reading, we will prove Theorem 2 of the previous class.
More informationWavelets and Multiresolution Processing. Thinh Nguyen
Wavelets and Multiresolution Processing Thinh Nguyen Multiresolution Analysis (MRA) A scaling function is used to create a series of approximations of a function or image, each differing by a factor of
More informationA Novel Fast Computing Method for Framelet Coefficients
American Journal of Applied Sciences 5 (11): 15-157, 008 ISSN 1546-939 008 Science Publications A Novel Fast Computing Method for Framelet Coefficients Hadeel N. Al-Taai Department of Electrical and Electronic
More informationAn Adaptive Pseudo-Wavelet Approach for Solving Nonlinear Partial Differential Equations
An Adaptive Pseudo-Wavelet Approach for Solving Nonlinear Partial Differential Equations Gregory Beylkin and James M. Keiser Wavelet Analysis and Applications, v.6, 1997, Academic Press. Contents 1 Introduction
More informationQuantum Computing Lecture 2. Review of Linear Algebra
Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces
More information7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim
69 Appendix F Molecular Dynamics F. Introduction In this chapter, we deal with the theories and techniques used in molecular dynamics simulation. The fundamental dynamics equations of any system is the
More informationBand-limited Wavelets and Framelets in Low Dimensions
Band-limited Wavelets and Framelets in Low Dimensions Likun Hou a, Hui Ji a, a Department of Mathematics, National University of Singapore, 10 Lower Kent Ridge Road, Singapore 119076 Abstract In this paper,
More informationFirst Year Seminar. Dario Rosa Milano, Thursday, September 27th, 2012
dario.rosa@mib.infn.it Dipartimento di Fisica, Università degli studi di Milano Bicocca Milano, Thursday, September 27th, 2012 1 Holomorphic Chern-Simons theory (HCS) Strategy of solution and results 2
More informationElectromagnetic Theory I
Electromagnetic Theory I Final Examination 18 December 2009, 12:30-2:30 pm Instructions: Answer the following 10 questions, each of which is worth 10 points. Explain your reasoning in each case. Use SI
More informationWavelets, wavelet networks and the conformal group
Wavelets, wavelet networks and the conformal group R. Vilela Mendes CMAF, University of Lisbon http://label2.ist.utl.pt/vilela/ April 2016 () April 2016 1 / 32 Contents Wavelets: Continuous and discrete
More informationWavelets in abstract Hilbert space
Wavelets in abstract Hilbert space Mathieu Sablik Mathematiques, ENS Lyon, 46 allee d Italie, F69364 LYON Cedex 07, France June-July 2000 Introduction The purpose of my training period has been to provide
More informationStrengthened Sobolev inequalities for a random subspace of functions
Strengthened Sobolev inequalities for a random subspace of functions Rachel Ward University of Texas at Austin April 2013 2 Discrete Sobolev inequalities Proposition (Sobolev inequality for discrete images)
More informationSpin one matter elds. November 2015
Spin one matter elds M. Napsuciale, S. Rodriguez, R.Ferro-Hernández, S. Gomez-Ávila Universidad de Guanajuato Mexican Workshop on Particles and Fields November 2015 M. Napsuciale, S. Rodriguez, R.Ferro-Hernández,
More informationMODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione
MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione, Univ. di Roma Tor Vergata, via di Tor Vergata 11,
More informationDigital Affine Shear Filter Banks with 2-Layer Structure
Digital Affine Shear Filter Banks with -Layer Structure Zhihua Che and Xiaosheng Zhuang Department of Mathematics, City University of Hong Kong, Tat Chee Avenue, Kowloon Tong, Hong Kong Email: zhihuache-c@my.cityu.edu.hk,
More informationDiffusion Tensor Imaging I. Jennifer Campbell
Diffusion Tensor Imaging I Jennifer Campbell Diffusion Imaging Molecular diffusion The diffusion tensor Diffusion weighting in MRI Alternatives to the tensor Overview of applications Diffusion Imaging
More informationSparsifying Transform Learning for Compressed Sensing MRI
Sparsifying Transform Learning for Compressed Sensing MRI Saiprasad Ravishankar and Yoram Bresler Department of Electrical and Computer Engineering and Coordinated Science Laborarory University of Illinois
More information1 First and second variational formulas for area
1 First and second variational formulas for area In this chapter, we will derive the first and second variational formulas for the area of a submanifold. This will be useful in our later discussion on
More informationPreconditioning techniques in frame theory and probabilistic frames
Preconditioning techniques in frame theory and probabilistic frames Department of Mathematics & Norbert Wiener Center University of Maryland, College Park AMS Short Course on Finite Frame Theory: A Complete
More informationExpansion of 1/r potential in Legendre polynomials
Expansion of 1/r potential in Legendre polynomials In electrostatics and gravitation, we see scalar potentials of the form V = K d Take d = R r = R 2 2Rr cos θ + r 2 = R 1 2 r R cos θ + r R )2 Use h =
More informationNumerical Methods - Numerical Linear Algebra
Numerical Methods - Numerical Linear Algebra Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Numerical Linear Algebra I 2013 1 / 62 Outline 1 Motivation 2 Solving Linear
More informationAn Introduction to Wavelets and some Applications
An Introduction to Wavelets and some Applications Milan, May 2003 Anestis Antoniadis Laboratoire IMAG-LMC University Joseph Fourier Grenoble, France An Introduction to Wavelets and some Applications p.1/54
More informationPolynomial functions over nite commutative rings
Polynomial functions over nite commutative rings Balázs Bulyovszky a, Gábor Horváth a, a Institute of Mathematics, University of Debrecen, Pf. 400, Debrecen, 4002, Hungary Abstract We prove a necessary
More information446 SCIENCE IN CHINA (Series F) Vol. 46 introduced in refs. [6, ]. Based on this inequality, we add normalization condition, symmetric conditions and
Vol. 46 No. 6 SCIENCE IN CHINA (Series F) December 003 Construction for a class of smooth wavelet tight frames PENG Lizhong (Λ Π) & WANG Haihui (Ξ ) LMAM, School of Mathematical Sciences, Peking University,
More informationV. SUBSPACES AND ORTHOGONAL PROJECTION
V. SUBSPACES AND ORTHOGONAL PROJECTION In this chapter we will discuss the concept of subspace of Hilbert space, introduce a series of subspaces related to Haar wavelet, explore the orthogonal projection
More informationLecture 10: Dimension Reduction Techniques
Lecture 10: Dimension Reduction Techniques Radu Balan Department of Mathematics, AMSC, CSCAMM and NWC University of Maryland, College Park, MD April 17, 2018 Input Data It is assumed that there is a set
More information1 Introduction. Aurobrata Ghosh 1, Elias Tsigaridas 2, Maxime Descoteaux 3, Pierre Comon 4, Bernard Mourrain 2, and Rachid Deriche 1
A polynomial based approach to extract the maxima of an antipodally symmetric spherical function and its application to extract fiber directions from the Orientation Distribution Function in Diffusion
More informationEigenfunctions on the surface of a sphere. In spherical coordinates, the Laplacian is. u = u rr + 2 r u r + 1 r 2. sin 2 (θ) + 1
Eigenfunctions on the surface of a sphere In spherical coordinates, the Laplacian is u = u rr + 2 r u r + 1 r 2 [ uφφ sin 2 (θ) + 1 sin θ (sin θ u θ) θ ]. Eigenfunctions on the surface of a sphere In spherical
More informationRepresentation theory and quantum mechanics tutorial Spin and the hydrogen atom
Representation theory and quantum mechanics tutorial Spin and the hydrogen atom Justin Campbell August 3, 2017 1 Representations of SU 2 and SO 3 (R) 1.1 The following observation is long overdue. Proposition
More informationJoint Optimum Bitwise Decomposition of any. Memoryless Source to be Sent over a BSC. Ecole Nationale Superieure des Telecommunications URA CNRS 820
Joint Optimum Bitwise Decomposition of any Memoryless Source to be Sent over a BSC Seyed Bahram Zahir Azami, Pierre Duhamel 2 and Olivier Rioul 3 cole Nationale Superieure des Telecommunications URA CNRS
More informationLecture 4 Quantum mechanics in more than one-dimension
Lecture 4 Quantum mechanics in more than one-dimension Background Previously, we have addressed quantum mechanics of 1d systems and explored bound and unbound (scattering) states. Although general concepts
More informationWavelets on Two-Dimensional Manifolds
Technical University of Cluj-Napoca, Romania HABILITATION THESIS Wavelets on Two-Dimensional Manifolds author Daniela Roşca 0 Contents Introduction iii Wavelets on domains of R. Piecewise constant wavelet
More informationMIT 9.520/6.860, Fall 2017 Statistical Learning Theory and Applications. Class 19: Data Representation by Design
MIT 9.520/6.860, Fall 2017 Statistical Learning Theory and Applications Class 19: Data Representation by Design What is data representation? Let X be a data-space X M (M) F (M) X A data representation
More informationQuantum Mechanics for Mathematicians: Energy, Momentum, and the Quantum Free Particle
Quantum Mechanics for Mathematicians: Energy, Momentum, and the Quantum Free Particle Peter Woit Department of Mathematics, Columbia University woit@math.columbia.edu November 28, 2012 We ll now turn to
More informationRecent developments on sparse representation
Recent developments on sparse representation Zeng Tieyong Department of Mathematics, Hong Kong Baptist University Email: zeng@hkbu.edu.hk Hong Kong Baptist University Dec. 8, 2008 First Previous Next Last
More informationTHE GEOMETRY AND THE ANALYTIC PROPERTIES OF ISOTROPIC MULTIRESOLUTION ANALYSIS
THE GEOMETRY AND THE ANALYTIC PROPERTIES OF ISOTROPIC MULTIRESOLUTION ANALYSIS JUAN R. ROMERO, SIMON K. ALEXANDER, SHIKHA BAID, SAURABH JAIN, AND MANOS PAPADAKIS ABSTRACT. In this paper we investigate
More informationA Bregman alternating direction method of multipliers for sparse probabilistic Boolean network problem
A Bregman alternating direction method of multipliers for sparse probabilistic Boolean network problem Kangkang Deng, Zheng Peng Abstract: The main task of genetic regulatory networks is to construct a
More informationMultiscale Geometric Analysis: Thoughts and Applications (a summary)
Multiscale Geometric Analysis: Thoughts and Applications (a summary) Anestis Antoniadis, University Joseph Fourier Assimage 2005,Chamrousse, February 2005 Classical Multiscale Analysis Wavelets: Enormous
More informationIMA Preprint Series # 2256
MULTIPLE q-shell ODF RECONSTRUCTION IN q-ball IMAGING By Iman Aganj Christophe Lenglet Guillermo Sapiro Essa Yacoub Kamil Ugurbil and Noam Harel IMA Preprint Series # 2256 ( May 2009 ) INSTITUTE FOR MATHEMATICS
More informationLinear Algebra Notes. Lecture Notes, University of Toronto, Fall 2016
Linear Algebra Notes Lecture Notes, University of Toronto, Fall 2016 (Ctd ) 11 Isomorphisms 1 Linear maps Definition 11 An invertible linear map T : V W is called a linear isomorphism from V to W Etymology:
More informationORTHOGONAL WAVELET FRAMES AND VECTOR-VALUED WAVELET TRANSFORMS
ORTHOGONAL WAVELET FRAMES AND VECTOR-VALUED WAVELET TRANSFORMS GHANSHYAM BHATT, BRODY DYLAN JOHNSON, AND ERIC WEBER Abstract Motivated by the notion of orthogonal frames, we describe sufficient conditions
More informationCONSTRUCTION OF AN ORTHONORMAL COMPLEX MULTIRESOLUTION ANALYSIS. Liying Wei and Thierry Blu
CONSTRUCTION OF AN ORTHONORMAL COMPLEX MULTIRESOLUTION ANALYSIS Liying Wei and Thierry Blu Department of Electronic Engineering, The Chinese University of Hong Kong, Shatin, N.T., Hong Kong ABSTRACT We
More informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods
More information1 The Continuous Wavelet Transform The continuous wavelet transform (CWT) Discretisation of the CWT... 2
Contents 1 The Continuous Wavelet Transform 1 1.1 The continuous wavelet transform (CWT)............. 1 1. Discretisation of the CWT...................... Stationary wavelet transform or redundant wavelet
More informationSupplementary Material for Nonparametric Operator-Regularized Covariance Function Estimation for Functional Data
Supplementary Material for Nonparametric Operator-Regularized Covariance Function Estimation for Functional Data Raymond K. W. Wong Department of Statistics, Texas A&M University Xiaoke Zhang Department
More informationChapter 7 Wavelets and Multiresolution Processing. Subband coding Quadrature mirror filtering Pyramid image processing
Chapter 7 Wavelets and Multiresolution Processing Wavelet transform vs Fourier transform Basis functions are small waves called wavelet with different frequency and limited duration Multiresolution theory:
More informationThe Generalized Haar-Walsh Transform (GHWT) for Data Analysis on Graphs and Networks
The Generalized Haar-Walsh Transform (GHWT) for Data Analysis on Graphs and Networks Jeff Irion & Naoki Saito Department of Mathematics University of California, Davis SIAM Annual Meeting 2014 Chicago,
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationTutorial: Sparse Signal Processing Part 1: Sparse Signal Representation. Pier Luigi Dragotti Imperial College London
Tutorial: Sparse Signal Processing Part 1: Sparse Signal Representation Pier Luigi Dragotti Imperial College London Outline Part 1: Sparse Signal Representation ~90min Part 2: Sparse Sampling ~90min 2
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationStatistical Learning & Applications. f w (x) =< f w, K x > H = w T x. α i α j < x i x T i, x j x T j. = < α i x i x T i, α j x j x T j > F
CR2: Statistical Learning & Applications Examples of Kernels and Unsupervised Learning Lecturer: Julien Mairal Scribes: Rémi De Joannis de Verclos & Karthik Srikanta Kernel Inventory Linear Kernel The
More informationClass Meeting # 12: Kirchhoff s Formula and Minkowskian Geometry
MATH 8.52 COURSE NOTES - CLASS MEETING # 2 8.52 Introduction to PDEs, Spring 207 Professor: Jared Speck Class Meeting # 2: Kirchhoff s Formula and Minkowskian Geometry. Kirchhoff s Formula We are now ready
More informationIsotropic Multiresolution Analysis: Theory and Applications
Isotropic Multiresolution Analysis: Theory and Applications Saurabh Jain Department of Mathematics University of Houston March 17th 2009 Banff International Research Station Workshop on Frames from first
More informationBoundary Value Problems - Solving 3-D Finite-Difference problems Jacob White
Introduction to Simulation - Lecture 2 Boundary Value Problems - Solving 3-D Finite-Difference problems Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Reminder about
More informationSparsity- and continuity-promoting seismic image recovery. with curvelet frames
Sparsity- and continuity-promoting seismic image recovery with curvelet frames Felix J. Herrmann, Peyman Moghaddam and Chris Stolk Seismic Laboratory for Imaging and Modeling, Department of Earth and Ocean
More informationPAIRS OF DUAL PERIODIC FRAMES
PAIRS OF DUAL PERIODIC FRAMES OLE CHRISTENSEN AND SAY SONG GOH Abstract. The time-frequency analysis of a signal is often performed via a series expansion arising from well-localized building blocks. Typically,
More informationL. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS
Rend. Sem. Mat. Univ. Pol. Torino Vol. 57, 1999) L. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS Abstract. We use an abstract framework to obtain a multilevel decomposition of a variety
More information