Multiscale Frame-based Kernels for Image Registration
|
|
- Kevin Cooper
- 5 years ago
- Views:
Transcription
1 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 Registration 22 July, 16 1 / 31
2 Outline 1 Motivation 2 Frame Theory, Reproducing Kernels and Frames 3 Eg. Spline-based Reproducing Kernels Conditions and Interpretations Relationship across different scales 4 Simple Simulations Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 2 / 31
3 Kernels and Deformation Fields LDDMM formulation: Vector field, v t, is constrained to an RKHS, (V,, V ), such that it has the form: v t( ) = N K(g t(x n), )α n,t (1) n=1 1 K : Ω Ω R d, a matrix-valued function, is the reproducing kernel associated with the RKHS V. (*Here, K : Ω Ω R is a scalar-valued function multiplied with a d d identity matrix.) 2 K determines the properties of v t. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 3 / 31
4 ... from Multiscale Gaussian Kernels 1 Deformable objects deformed with the LDDMM framework are commonly used with Gaussian kernels of chosen bandwidths. 2 multiscale Gaussian: weighted sum of Gaussian kernels of varying bandwidths such that both smooth and sharp variations can be captured within a single deformation map. [Risser et al. ] Figure: Taken from Risser et al. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 4 / 31
5 to Multiscale frame-based Kernels 1 Propose to use a generic multiscale reproducing kernel that is frame-based 1 Exploit known relationships between b-splines, Riesz bases, frames to build RKHS from spline-based frames. 2 Multiscale structure 3 Compact 4 Example 1: compact b-splines (denoted as φ) are dilated iteratively using pre-defined refinement masks M : ˆφ(2 ) ˆφ( ) and combined together like the sum of Gaussians. Example 2: wavelet masks W : ˆφ(2 ) ˆψ( ) that generate compactly supported spline-based functions (wavelet frames) (a) Refinement Mask (b) Wavelet Mask Figure: Refinement masks and wavelet masks; B-spline of order 2 Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 / 31
6 Literature Review Reproducing kernels from Frames Construction of reproducing kernels, and consequently, RKHS, is well-studied in the machine learning literature, within the regularisation theory framework for function approximations. 1 In Rakotomamonjy(0), general conditions for the construction of arbitrary reproducing kernel Hilbert spaces from frame and dual frame elements were given. 2 Wavelet frames with multiresolution properties were used to design kernels for classification or regression problems utilising support vector machines (Zhang, 04; Opfer, 06; Ling, 06; Shang, 08; Zhang, ) 3 In this work, we consider tight frames constructed from b-splines and spline framelets (Ron, 199; Ron, 1997; Daubechies, 03) with multiscale structures inherited from wavelet frames. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 6 / 31
7 Preliminaries - Frame Theory - Operators 1 ϕ := {ϕ i } i I - a family of functions, where I is a finite or infinitely countable index set. 2 H - a Hilbert space with inner product, H. Let ϕ i H, f H and c := {c i } i I l 2. We define a synthesis operator T ϕ and an adjoint decomposition operator Tϕ such that T ϕ : l 2 (ϕ) H : {c i } i I c i ϕ i (2) i I T ϕ : H l 2 (ϕ) : f { f,ϕ i H } i I (3) Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 7 / 31
8 Preliminaries - Frame Theory - Definition of Frames Definition of Frames ϕ is a frame of (H,, H ) if and only if the frame condition C 1 f 2 H f,ϕ i H 2 C 2 f 2 H (4) i I holds for some C 1,C 2 such that 0 < C 1 C 2 <. 1 Tϕ is injective. 2 T ϕ is surjective. 3 When C 1 = C 2, ϕ is known as a tight frame. If C 1 = C 2 = 1, ϕ is supertight/parseval ( ϕ i H = c). Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 8 / 31
9 Preliminaries - Frame Theory - Frame Operator Frame Operator T ϕtϕ : H H : f f,ϕ i H ϕ i () i I 1 It is continuously invertible. 2 We denote its inverse as R := (T ϕt ϕ ) 1, R = R (self-adjoint). 3 T Rϕ is the right inverse of Tϕ and TϕT Rϕ (f ) = i I f,rϕ i H ϕ i = f. 4 Rϕ is known as the dual frame of ϕ and commonly denoted as ϕ. Rϕ helps to reconstruct any f H with the frames ϕ by first decomposing f into a set of specific coefficients { f,rϕ i } i I, known as the canonical frame coefficients. 6 The canonical frame coefficients { f, ϕ i } i I l 2 has the smallest sum of squares amongst all possible coefficients e l 2, where e := {e i } i I,T ϕ(e) = f, i.e. { f, ϕ i } i I 2 l 2 e 2 l 2. 7 Every frame has a dual frame and every dual frame is the frame of some dual frame. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 9 / 31
10 Reproducing kernels and Frames (Rakotomamonjy, 0) Let ϕ be a known set of frames for Hilbert space (H,, H ). Then H is an RKHS with the reproducing kernel K(x,y) = ϕ i (y)ϕ i (x) (6) i I 1 Every finite set of vectors {ϕ i } i I in a vector space equipped with an inner product is a frame for the Hilbert space H := span{ϕ i } i I (Christensen, 08) and every frame has a dual frame. It is possible to build an RKHS with a reproducing kernel of the form (6) (but it s feasible only if the dual frame has an analytical form). 2 Given a (super)tight frame ϕ, i.e. ϕ = ϕ of a Hilbert space (H,, H ), then the corresponding reproducing kernel is of the form K(x,y) = ϕ i (y)ϕ i (x) (7) i I Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 / 31
11 Sum of Reproducing Kernels (Berlinet, 11) Let K 1 and K 2 be reproducing kernels of spaces H 1 and H 2 of functions f 1,f 2 R Ω. Then K = K 1 + K 2 is the reproducing kernel of the space H = H 1 H 2 = {f f = f 1 + f 2,f 1 H 1,f 2 H 2 } with the norm H defined by ) f 2 H = min ( f 1 2 H 1 + f 2 2 H 2 : f = f 1 + f 2,f 1 H 1,f 2 H 2 1 Let {ϕ i } i I and {ϕ i } i I be two sets of basis vectors and H 1 := span{ϕ i } i I, H 2 := span{ϕ i } i I. Then K(x,y) := i (I I ) ϕ i (x)ϕ i (y) is the reproducing kernel of H := H 1 H 2 = span{ϕ i } i (I I ) with the norm f 2 H = min( f 1 2 H 1 + f 2 2 H 2 : f = f 1 + f 2,f 1 H 1,f 2 H 2 ). Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
12 Reproducing kernels and tight frames Let ϕ be a known set of frames for Hilbert space (H,, H ). Then H is an RKHS with the reproducing kernel K(x,y) = ϕ i (y)ϕ i (x) (8) i I if and only if ϕ = ϕ. 1 The inner product, H is equivalent to, l2, since f,f H = f,ϕ i 2 H = ci 2 = c i,c i l2 (9) i I i I Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
13 Reproducing kernels and tight frames (Continued) 1 We have f ( ) = N α nk(x n, ) = c i ϕ i ( ) n=1 where c := {c i },c i = N n=1 αnϕ i (x n). 2 The frame coefficients live in the orthogonal complement of the null space of T ϕ. 3 c are the canonical frame coefficients of ϕ. 4 The reproducing kernel guarantees the l 2 -minimisation of the frame coefficients. i I Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
14 Weighted Sum of Reproducing Kernels Consider a more general form of RK that admits a positive weight parameter λ i > 0 for ϕ i, i I λ i < : K λ (x,y) = λ i ϕ i (y)ϕ i (x) () i I 1 Tunable weights allow us to further control the characteristics or properties of a reproducing kernel constructed via scaled summations. 2 Not unlike the sum of weighted Gaussians, it makes sense to put stronger emphasis on the larger bandwidths to build a smoother function space. 3 Suitable conditions on the weights can produce functions that live in a Sobolev space of arbitrary smoothness. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
15 Weighted Sum of Reproducing Kernels (Continued) Let i I λ i ϕ i (x) 2 <. K λ is admissible and corresponds to (is the RK of) some Hilbert space, i.e. the function space spanned by { λ i ϕ i } i I. Then { λ i ϕ i } i I is a supertight/parseval frame of (H λ,, Hλ ): { } c H λ = c i ϕ i : i 2 < (11) λ i i I i I f,g Hλ : = P (c),p (e) l2,λ (12) { } c l 2,λ : = c = {c i } i I : i 2 < (13) λ i i I c c,e l2,λ = i e i (14) λ i i I P is the projection into the orthogonal complement of the null space of the operator T ϕ,λ, s.t. T ϕ,λ (c) = f, T ϕ,λ (e) = g. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 / 31
16 Eg. Spline-based Reproducing Kernels 1 Certain combinations of b-splines gives us tight frames. 2 B-splines are refinable and compact. 3 Their refinability gives a natural multiscale structure that incorporates a relationship across the different scales. 4 B-splines can be shown to converge asymptotically to Gaussians. B-splines are scaling functions of multiresolution analysis (MRA) and can be used to build spline wavelets with natural multiresolution structure. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
17 Example 1 - Multiscale b-splines (Rakotomamonjy, 0, Opfer, 06, Han, 09) Let φ m j,k := 2jd/2 φ m (2 j k), where φ m is an univariate b-spline of order m, λ j > 0, j Z λ j <. Then { λ j φ m j,k } j,k is a supertight frame of the Hilbert space (B,, B ): B = f : f = c j,k φ m (2 j c j,k 2 k), λ j 2 jd < () j Z k Z d j Z k Z d f 2 c j,k 2 B : = min λ j 2 jd = f 2 V j λ j 2 jd : f j V j := span{φ m (2 j k)} (16) j Z k Z d j Z with the reproducing kernel K S (x,y) = λ j φ m j,k(x)φ m j,k(y) = λ j 2 jd K S,j (x,y), (17) j Z,k Z d j Z where K S,j (x,y) = φ m (2 j x k)φ m (2 j y k). (18) k Z d Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
18 Example 1(Multiscale Interpretation) Using the kernel K S, v t B takes the form N v t = λ j 2 jd α n,tφ m (2 j x n k)φ m (2 j k) = λ j 2 jd v t,j j Z k Z d n=1 j Z The solution can be interpreted as a weighted sum of functions v t,j, where each v t,j belongs to span{φ m (2 j k)}, and can be constructed with the canonical coefficients c j,k = N α n,tφ m (2 j x n k) n=1 Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
19 Example 2 - with Spline Framelets Let φ m be the centred b-spline of order m for d = 1. Then φ is refinable with the refinement mask ˆτ 0 (ξ) := e iξm 2 cos m ( ξ 2 ) and it forms an MRA, {V j} j Z of L 2 (R). Then the set of m functions Ψ m = {ψl m : l = 1,...,m} defined (in the Fourier domain) as ˆψ l m (ξ) := i l e iξm 2 ( m l ) cos m l (ξ/4)sin m+l (ξ/4) (ξ/4) m (19) is a set of generators for a wavelet system X(Ψ m ) := {ψ m j,k,l := 2j/2 ψ m l (2j k) : 1 l m;j,k Z}. 1 X(Ψ m ) L 2 (R) is a tight frame in L 2 (R) that is constructed from MRA (Ron, 1997). 2 We refer to the generators Ψ m as spline framelets. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
20 Example 2 (continued) The reproducing kernel for (L 2 (R),, ) is defined as K F (x,y) = ψj,k,l(x)ψ m j,k,l(y) m () j Z k Z l 1,...,m f,g : = f,ψj,k,l g,ψ m j,k,l m (21) j Z k Z l 1,...,m Since ψ m was created using wavelet masks from the scaling function φ m, the reproducing kernel given in () can be truncated from below: K F (x,y) = φ m j min,k(x)φ m j min,k(y) + ψj,k,l(x)ψ m j,k,l(y) m (22) k Z j Z j min k Z l 1,...,m Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 / 31
21 Example 2 (continued) Extension from L 2 (R) to L 2 (R d ) For simplicity, we consider the case d = 2. Using the tensor products of spline framelets in L 2 (R), the spline framelets for L 2 (R 2 ) are defined as ψ m l,l (x) := ψm l (x 1 )ψ m l (x 2 ), x = (x 1,x 2 ) R 2, 0 l,l m, ψ0 m = φm. Let Ψ m l,l framelets for 0 l,l m, (l,l ) (0,0). Then Ψ m l,l L 2 (R d ) with the reproducing kernel K F,2 (x,y) = j Z k Z 2 (l,l ) L where L := {(l,l ) : 0 l,l m, (l,l ) (0,0)}. be the set of newly defined spline is a tight wavelet frame for ψj,k,(l,l m ) (x)ψm j,k,(l,l )(y) (23) Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
22 Spline-based Kernel Computation - Worked Example Plot of Kernel Splines τ 0 φ τ 1 φ τ 2 φ τ φ 3 τ 4 φ A worked example to Examples 1 and 2 where m = 4. τ 0 = 16 1 [1,4,6,4,1] τ 1 = 16 1 [1, 4,6, 4,1] τ 2 = 1 8 [1, 2,0,2, 1] τ 3 = 6 16 [1,0, 2,0,1] τ 4 = 1 8 [1,2,0, 2, 1] (24) Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
23 Asymptotic properties of b-splines and spline framelets 1 Normalised b-splines are point-wise close to some Gaussian function as m goes to infinity. Specifically, it has been shown that mφ m ( mx) converges uniformly to 6 π e 6x2 as m (Chen, 04). 2 The spline framelets ψl m, 1 l m converges to the lth derivative of some Gaussian function G m,l (Shen, 13) (m G m,l := 1 ) ( ) 6 4 l l π(m l/2) exp 6x 2 m l/2 () Specifically, let G m l dl := dx G l m,l (x m /2). Then we have the inequality max max ψl m (x) Gm l 1 l m x R (x) (lnm)/2 m 3/2 Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
24 Weights and Frames of Sobolev Space Under specific conditions on λ j, { λ j φ m j,k } j,k is also a frame of a Sobolev space of arbitrary smoothness. For s R, we define the Sobolev space H s (R d ) as a space consisting of all tempered distributions f such that f 2 H s (R d ) := 1 (2π) d ˆf 2 (1 + ξ 2 2) s dξ < (26) R d 1 The wavelet system {2 j(d/2 s) φ m (2 j k) : j Z,k Z d } is a wavelet frame in H s (R d ) for all 0 < s < m 1/2 (Han, 09). 2 The upper bound of s is the regularity of the m order b-spline. 3 { λ j φ m j,k } j,k is a frame of H s (R d ) when λ j := 2 js. 4 By varying λ j, a level of smoothness can be imposed on the vector fields. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
25 Relationship Across Scales Let d = 1. Consider the masks {h l } l {0,1,...,m} such that { φ m (x) = k Z h 0,kφ m (2x k) ψ m l (x) = k Z h l,kφ m (2x k) (27) (Relationship between coefficients) Let e j,k,l := N n=1 α(xn)ψm l (2j x n k) and e j,k,0 := N n=1 α(xn)φm (2 j x n k), 1 l m. Then for 0 l m, e j,k,l = h l,k 2ke j+1,k,0 k Z Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 / 31
26 Relationship Across Scales (Relationship between coefficients of derivatives) Using the same masks {h l } l {0,1,...,m}, and the notation D x f = df /dx, we use the following identities to compute D x v. D x ψl m (x) = k Z 2h l,kd x φ m (2x k) e j,k,l = k Z 2h l,k 2ke j+1,k,0 = N n=1 α(xn)dxn ψl m (2j x n k) where, 0 l m, ψ m 0 := φm. (28) Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
27 Test Results (a) (b) (c) (d) (e) Figure: Visualisation of kernels of different scale. The kernel used here is K S, with a b-spline of order 4. (a) Template Image (b) Target Image (c) g 1 (x) for j min = j max = 2 (d) g 1 (x) for j min = j max = 3 (e) g 1 (x) for j min = j max = 4. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
28 Test Results (a) (b) (c) (d) Figure: Examples of K S in LDDMM: Panels (a) Template Image, (b) Target Image, (c) Deformed grid points using K S, j min = 3, j max = 0 (d) Deformed template using K S. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
29 Test Results (a) (b) (c) (d) (e) (f) (g) (h) Figure: A simulation of a deformable diffeomorphic image matching algorithm. Panels (a) Template Image, (b) Target Image, (c,e,g) Deformed grid points using Gaussian kernels, (d,f,h) Deformed template using Gaussian kernels. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
30 Test Results (a) (b) (c) (d) Figure: A simulation of a deformable diffeomorphic image matching algorithm using proposed kernels. Panels (i,j) and (k,l) shows the deformed grid points and deformed template using K S and K F respectively. For K S and K F, s = 3, j min = 3, j max = 3. Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, 16 / 31
31 Thank you for your attention Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image Registration 22 July, / 31
Reproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 11, 2009 About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
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 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 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 informationBAND-LIMITED REFINABLE FUNCTIONS FOR WAVELETS AND FRAMELETS
BAND-LIMITED REFINABLE FUNCTIONS FOR WAVELETS AND FRAMELETS WEIQIANG CHEN AND SAY SONG GOH DEPARTMENT OF MATHEMATICS NATIONAL UNIVERSITY OF SINGAPORE 10 KENT RIDGE CRESCENT, SINGAPORE 119260 REPUBLIC OF
More informationReproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 12, 2007 About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationWAVELETS WITH SHORT SUPPORT
WAVELETS WITH SHORT SUPPORT BIN HAN AND ZUOWEI SHEN Abstract. This paper is to construct Riesz wavelets with short support. Riesz wavelets with short support are of interests in both theory and application.
More informationReproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 9, 2011 About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
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 informationConstruction of Multivariate Compactly Supported Tight Wavelet Frames
Construction of Multivariate Compactly Supported Tight Wavelet Frames Ming-Jun Lai and Joachim Stöckler April 5, 2006 Abstract Two simple constructive methods are presented to compute compactly supported
More informationJournal of Mathematical Analysis and Applications. Properties of oblique dual frames in shift-invariant systems
J. Math. Anal. Appl. 356 (2009) 346 354 Contents lists available at ScienceDirect Journal of Mathematical Analysis and Applications www.elsevier.com/locate/jmaa Properties of oblique dual frames in shift-invariant
More informationRKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets Class 22, 2004 Tomaso Poggio and Sayan Mukherjee
RKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets 9.520 Class 22, 2004 Tomaso Poggio and Sayan Mukherjee About this class Goal To introduce an alternate perspective of RKHS via integral operators
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationAtomic decompositions of square-integrable functions
Atomic decompositions of square-integrable functions Jordy van Velthoven Abstract This report serves as a survey for the discrete expansion of square-integrable functions of one real variable on an interval
More informationDirect Learning: Linear Classification. Donglin Zeng, Department of Biostatistics, University of North Carolina
Direct Learning: Linear Classification Logistic regression models for classification problem We consider two class problem: Y {0, 1}. The Bayes rule for the classification is I(P(Y = 1 X = x) > 1/2) so
More informationConstruction of wavelets. Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam
Construction of wavelets Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam Contents Stability of biorthogonal wavelets. Examples on IR, (0, 1), and (0, 1) n. General domains
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 informationReproducing Kernel Hilbert Spaces Class 03, 15 February 2006 Andrea Caponnetto
Reproducing Kernel Hilbert Spaces 9.520 Class 03, 15 February 2006 Andrea Caponnetto About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationHilbert Spaces: Infinite-Dimensional Vector Spaces
Hilbert Spaces: Infinite-Dimensional Vector Spaces PHYS 500 - Southern Illinois University October 27, 2016 PHYS 500 - Southern Illinois University Hilbert Spaces: Infinite-Dimensional Vector Spaces October
More informationConstruction of Biorthogonal Wavelets from Pseudo-splines
Construction of Biorthogonal Wavelets from Pseudo-splines Bin Dong a, Zuowei Shen b, a Department of Mathematics, National University of Singapore, Science Drive 2, Singapore, 117543. b Department of Mathematics,
More informationCHAPTER VIII HILBERT SPACES
CHAPTER VIII HILBERT SPACES DEFINITION Let X and Y be two complex vector spaces. A map T : X Y is called a conjugate-linear transformation if it is a reallinear transformation from X into Y, and if T (λx)
More informationFrames. Hongkai Xiong 熊红凯 Department of Electronic Engineering Shanghai Jiao Tong University
Frames Hongkai Xiong 熊红凯 http://ivm.sjtu.edu.cn Department of Electronic Engineering Shanghai Jiao Tong University 2/39 Frames 1 2 3 Frames and Riesz Bases Translation-Invariant Dyadic Wavelet Transform
More informationReproducing Kernel Hilbert Spaces
9.520: Statistical Learning Theory and Applications February 10th, 2010 Reproducing Kernel Hilbert Spaces Lecturer: Lorenzo Rosasco Scribe: Greg Durrett 1 Introduction In the previous two lectures, we
More informationUNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES. Tight compactly supported wavelet frames of arbitrarily high smoothness
UNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES Tight compactly supported wavelet frames of arbitrarily high smoothness Karlheinz Gröchenig Amos Ron Department of Mathematics U-9 University
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 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 informationThe Method of Virtual Components in the Multivariate Setting
The Method of Virtual Components in the Multivariate Setting Ming-Jun Lai and Alexander Petukhov May 21, 2007 Abstract We describe the so-called method of virtual components for tight wavelet framelets
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 informationMRA Frame Wavelets with Certain Regularities Associated with the Refinable Generators of Shift Invariant Spaces
Chapter 6 MRA Frame Wavelets with Certain Regularities Associated with the Refinable Generators of Shift Invariant Spaces Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan University,
More informationBANACH FRAMES GENERATED BY COMPACT OPERATORS ASSOCIATED WITH A BOUNDARY VALUE PROBLEM
TWMS J. Pure Appl. Math., V.6, N.2, 205, pp.254-258 BRIEF PAPER BANACH FRAMES GENERATED BY COMPACT OPERATORS ASSOCIATED WITH A BOUNDARY VALUE PROBLEM L.K. VASHISHT Abstract. In this paper we give a type
More informationKernels A Machine Learning Overview
Kernels A Machine Learning Overview S.V.N. Vishy Vishwanathan vishy@axiom.anu.edu.au National ICT of Australia and Australian National University Thanks to Alex Smola, Stéphane Canu, Mike Jordan and Peter
More informationA Banach Gelfand Triple Framework for Regularization and App
A Banach Gelfand Triple Framework for Regularization and Hans G. Feichtinger 1 hans.feichtinger@univie.ac.at December 5, 2008 1 Work partially supported by EUCETIFA and MOHAWI Hans G. Feichtinger hans.feichtinger@univie.ac.at
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 informationA RECONSTRUCTION FORMULA FOR BAND LIMITED FUNCTIONS IN L 2 (R d )
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 127, Number 12, Pages 3593 3600 S 0002-9939(99)04938-2 Article electronically published on May 6, 1999 A RECONSTRUCTION FORMULA FOR AND LIMITED FUNCTIONS
More informationIntroduction to Hilbert Space Frames
to Hilbert Space Frames May 15, 2009 to Hilbert Space Frames What is a frame? Motivation Coefficient Representations The Frame Condition Bases A linearly dependent frame An infinite dimensional frame Reconstructing
More informationSMALL SUPPORT SPLINE RIESZ WAVELETS IN LOW DIMENSIONS BIN HAN, QUN MO, AND ZUOWEI SHEN
SMALL SUPPORT SPLINE RIESZ WAVELETS IN LOW DIMENSIONS BIN HAN, QUN MO, AND ZUOWEI SHEN Abstract. In [B. Han and Z. Shen, SIAM J. Math. Anal., 38 (006), 530 556], a family of univariate short support Riesz
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 informationSemi-orthogonal wavelet frames on positive half-line using the Walsh Fourier transform
NTMSCI 6, No., 175-183 018) 175 New Trends in Mathematical Sciences http://dx.doi.org/10.085/ntmsci.018.83 Semi-orthogonal wavelet frames on positive half-line using the Walsh Fourier transform Abdullah
More informationReproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 9, 2011 About this class Goal In this class we continue our journey in the world of RKHS. We discuss the Mercer theorem which gives
More informationElements of Positive Definite Kernel and Reproducing Kernel Hilbert Space
Elements of Positive Definite Kernel and Reproducing Kernel Hilbert Space Statistical Inference with Reproducing Kernel Hilbert Space Kenji Fukumizu Institute of Statistical Mathematics, ROIS Department
More informationA survey on frame theory. Frames in general Hilbert spaces. Construction of dual Gabor frames. Construction of tight wavelet frames
A survey on frame theory Frames in general Hilbert spaces Construction of dual Gabor frames Construction of tight wavelet frames Ole Christensen Technical University of Denmark Department of Mathematics
More informationWe have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma and (3.55) with j = 1. We can write any f W 1 as
88 CHAPTER 3. WAVELETS AND APPLICATIONS We have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma 3..7 and (3.55) with j =. We can write any f W as (3.58) f(ξ) = p(2ξ)ν(2ξ)
More informationWavelets: Theory and Applications. Somdatt Sharma
Wavelets: Theory and Applications Somdatt Sharma Department of Mathematics, Central University of Jammu, Jammu and Kashmir, India Email:somdattjammu@gmail.com Contents I 1 Representation of Functions 2
More informationKernel Methods. Machine Learning A W VO
Kernel Methods Machine Learning A 708.063 07W VO Outline 1. Dual representation 2. The kernel concept 3. Properties of kernels 4. Examples of kernel machines Kernel PCA Support vector regression (Relevance
More informationStatistical Convergence of Kernel CCA
Statistical Convergence of Kernel CCA Kenji Fukumizu Institute of Statistical Mathematics Tokyo 106-8569 Japan fukumizu@ism.ac.jp Francis R. Bach Centre de Morphologie Mathematique Ecole des Mines de Paris,
More informationMGA Tutorial, September 08, 2004 Construction of Wavelets. Jan-Olov Strömberg
MGA Tutorial, September 08, 2004 Construction of Wavelets Jan-Olov Strömberg Department of Mathematics Royal Institute of Technology (KTH) Stockholm, Sweden Department of Numerical Analysis and Computer
More informationExtremely local MR representations:
Extremely local MR representations: L-CAMP Youngmi Hur 1 & Amos Ron 2 Workshop on sparse representations: UMD, May 2005 1 Math, UW-Madison 2 CS, UW-Madison Wavelet and framelet constructions History bits
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 informationInfinite-dimensional Vector Spaces and Sequences
2 Infinite-dimensional Vector Spaces and Sequences After the introduction to frames in finite-dimensional vector spaces in Chapter 1, the rest of the book will deal with expansions in infinitedimensional
More informationContents. 0.1 Notation... 3
Contents 0.1 Notation........................................ 3 1 A Short Course on Frame Theory 4 1.1 Examples of Signal Expansions............................ 4 1.2 Signal Expansions in Finite-Dimensional
More informationDigital Gabor Filters Do Generate MRA-based Wavelet Tight Frames
Digital Gabor Filters Do Generate MRA-based Wavelet Tight Frames Hui Ji a,,zuoweishen a,yufeizhao a a Department of Mathematics, National University of Singapore, Singapore 7543 Abstract Gabor frames,
More informationThe Learning Problem and Regularization Class 03, 11 February 2004 Tomaso Poggio and Sayan Mukherjee
The Learning Problem and Regularization 9.520 Class 03, 11 February 2004 Tomaso Poggio and Sayan Mukherjee About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing
More informationSymmetric Wavelet Tight Frames with Two Generators
Symmetric Wavelet Tight Frames with Two Generators Ivan W. Selesnick Electrical and Computer Engineering Polytechnic University 6 Metrotech Center, Brooklyn, NY 11201, USA tel: 718 260-3416, fax: 718 260-3906
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationReproducing formulas associated with symbols
Reproducing formulas associated with symbols Filippo De Mari Ernesto De Vito Università di Genova, Italy Modern Methods of Time-Frequency Analysis II Workshop on Applied Coorbit space theory September
More informationNonlinear Approximation Schemes Associated With Nonseparable Wavelet Bi-frames
Nonlinear Approximation Schemes Associated With Nonseparable Wavelet Bi-frames Martin Ehler September 26, 2008 Abstract In the present paper, we study nonlinear approximation properties of multivariate
More informationData fitting by vector (V,f)-reproducing kernels
Data fitting by vector (V,f-reproducing kernels M-N. Benbourhim to appear in ESAIM.Proc 2007 Abstract In this paper we propose a constructive method to build vector reproducing kernels. We define the notion
More informationA Short Course on Frame Theory
A Short Course on Frame Theory Veniamin I. Morgenshtern and Helmut Bölcskei ETH Zurich, 8092 Zurich, Switzerland E-mail: {vmorgens, boelcskei}@nari.ee.ethz.ch April 2, 20 Hilbert spaces [, Def. 3.-] and
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 informationCOMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY
COMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY BIN HAN AND HUI JI Abstract. In this paper, we provide a family of compactly supported orthonormal complex wavelets with dilation
More informationWavelets Marialuce Graziadei
Wavelets Marialuce Graziadei 1. A brief summary 2. Vanishing moments 3. 2D-wavelets 4. Compression 5. De-noising 1 1. A brief summary φ(t): scaling function. For φ the 2-scale relation hold φ(t) = p k
More informationOn the Numerical Evaluation of Fractional Sobolev Norms. Carsten Burstedde. no. 268
On the Numerical Evaluation of Fractional Sobolev Norms Carsten Burstedde no. 268 Diese Arbeit ist mit Unterstützung des von der Deutschen Forschungsgemeinschaft getragenen Sonderforschungsbereiches 611
More informationEXAMPLES OF REFINABLE COMPONENTWISE POLYNOMIALS
EXAMPLES OF REFINABLE COMPONENTWISE POLYNOMIALS NING BI, BIN HAN, AND ZUOWEI SHEN Abstract. This short note presents four examples of compactly supported symmetric refinable componentwise polynomial functions:
More informationReal Variables # 10 : Hilbert Spaces II
randon ehring Real Variables # 0 : Hilbert Spaces II Exercise 20 For any sequence {f n } in H with f n = for all n, there exists f H and a subsequence {f nk } such that for all g H, one has lim (f n k,
More informationWavelets and regularization of the Cauchy problem for the Laplace equation
J. Math. Anal. Appl. 338 008440 1447 www.elsevier.com/locate/jmaa Wavelets and regularization of the Cauchy problem for the Laplace equation Chun-Yu Qiu, Chu-Li Fu School of Mathematics and Statistics,
More informationNEW CONSTRUCTIONS OF PIECEWISE-CONSTANT WAVELETS
NEW CONSTRUCTIONS OF PIECEWISE-CONSTANT WAVELETS YOUNGMI HUR AND AMOS RON Abstract. The classical Haar wavelet system of L 2 (R n ) is commonly considered to be very local in space. We introduce and study
More informationfor compression of Boundary Integral Operators. Steven Paul Nixon B.Sc.
Theory and Applications of the Multiwavelets for compression of Boundary Integral Operators. Steven Paul Nixon B.Sc. Institute for Materials Research School of Computing, Science & Engineering, University
More informationFourier analysis, measures, and distributions. Alan Haynes
Fourier analysis, measures, and distributions Alan Haynes 1 Mathematics of diffraction Physical diffraction As a physical phenomenon, diffraction refers to interference of waves passing through some medium
More informationNonstationary Tight Wavelet Frames, II: Unbounded Intervals 1)
Nonstationary Tight Wavelet Frames, : Unbounded ntervals 1) by Charles K. Chui 2), Wenjie He Department of Mathematics and Computer Science University of Missouri St. Louis St. Louis, MO 63121-4499, USA
More informationWavelets and Multiresolution Processing
Wavelets and Multiresolution Processing Wavelets Fourier transform has it basis functions in sinusoids Wavelets based on small waves of varying frequency and limited duration In addition to frequency,
More informationApproximately dual frames in Hilbert spaces and applications to Gabor frames
Approximately dual frames in Hilbert spaces and applications to Gabor frames Ole Christensen and Richard S. Laugesen October 22, 200 Abstract Approximately dual frames are studied in the Hilbert space
More informationOn Riesz-Fischer sequences and lower frame bounds
On Riesz-Fischer sequences and lower frame bounds P. Casazza, O. Christensen, S. Li, A. Lindner Abstract We investigate the consequences of the lower frame condition and the lower Riesz basis condition
More informationAn Introduction to Kernel Methods 1
An Introduction to Kernel Methods 1 Yuri Kalnishkan Technical Report CLRC TR 09 01 May 2009 Department of Computer Science Egham, Surrey TW20 0EX, England 1 This paper has been written for wiki project
More informationLectures notes. Rheology and Fluid Dynamics
ÉC O L E P O L Y T E C H N IQ U E FÉ DÉR A L E D E L A U S A N N E Christophe Ancey Laboratoire hydraulique environnementale (LHE) École Polytechnique Fédérale de Lausanne Écublens CH-05 Lausanne Lectures
More informationDigital Image Processing
Digital Image Processing Wavelets and Multiresolution Processing () Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science 2 Contents Image pyramids Subband coding
More informationAffine and Quasi-Affine Frames on Positive Half Line
Journal of Mathematical Extension Vol. 10, No. 3, (2016), 47-61 ISSN: 1735-8299 URL: http://www.ijmex.com Affine and Quasi-Affine Frames on Positive Half Line Abdullah Zakir Husain Delhi College-Delhi
More informationCompactly Supported Tight Frames Associated with Refinable Functions 1. Communicated by Guido L. Weiss Received July 27, 1999
Applied and Computational Harmonic Analysis 8, 93 319 (000) doi:10.1006/acha.000.0301, available online at http://www.idealibrary.com on LETTER TO THE EDITOR Compactly Supported Tight Frames Associated
More information5.6 Nonparametric Logistic Regression
5.6 onparametric Logistic Regression Dmitri Dranishnikov University of Florida Statistical Learning onparametric Logistic Regression onparametric? Doesnt mean that there are no parameters. Just means that
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 informationDUALITY PRINCIPLE IN g-frames
Palestine Journal of Mathematics Vol. 6(2)(2017), 403 411 Palestine Polytechnic University-PPU 2017 DUAITY PRINCIPE IN g-frames Amir Khosravi and Farkhondeh Takhteh Communicated by Akram Aldroubi MSC 2010
More informationConstruction of Orthonormal Quasi-Shearlets based on quincunx dilation subsampling
Construction of Orthonormal Quasi-Shearlets based on quincunx dilation subsampling Rujie Yin Department of Mathematics Duke University USA Email: rujie.yin@duke.edu arxiv:1602.04882v1 [math.fa] 16 Feb
More informationarxiv: v3 [math.na] 11 Oct 2017
Indirect Image Registration with Large Diffeomorphic Deformations Chong Chen and Ozan Öktem arxiv:1706.04048v3 [math.na] 11 Oct 2017 Abstract. The paper adapts the large deformation diffeomorphic metric
More informationKernel B Splines and Interpolation
Kernel B Splines and Interpolation M. Bozzini, L. Lenarduzzi and R. Schaback February 6, 5 Abstract This paper applies divided differences to conditionally positive definite kernels in order to generate
More informationDeterminant lines and determinant line bundles
CHAPTER Determinant lines and determinant line bundles This appendix is an exposition of G. Segal s work sketched in [?] on determinant line bundles over the moduli spaces of Riemann surfaces with parametrized
More informationAnalysis of Fractals, Image Compression and Entropy Encoding
Analysis of Fractals, Image Compression and Entropy Encoding Myung-Sin Song Southern Illinois University Edwardsville Jul 10, 2009 Joint work with Palle Jorgensen. Outline 1. Signal and Image processing,
More informationBERGMAN KERNEL ON COMPACT KÄHLER MANIFOLDS
BERGMAN KERNEL ON COMPACT KÄHLER MANIFOLDS SHOO SETO Abstract. These are the notes to an expository talk I plan to give at MGSC on Kähler Geometry aimed for beginning graduate students in hopes to motivate
More informationIntroduction to Signal Spaces
Introduction to Signal Spaces Selin Aviyente Department of Electrical and Computer Engineering Michigan State University January 12, 2010 Motivation Outline 1 Motivation 2 Vector Space 3 Inner Product
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 informationOn Dual Wavelet Tight Frames
On Dual Wavelet Tight Frames Bin Han Department of Mathematical Sciences University of Alberta Edmonton, AB T6G 2G, Canada Email: bhan@vega.math.ualberta.ca Abstract. A characterization of multivariate
More informationSparse Multidimensional Representation using Shearlets
Sparse Multidimensional Representation using Shearlets Demetrio Labate a, Wang-Q Lim b, Gitta Kutyniok c and Guido Weiss b, a Department of Mathematics, North Carolina State University, Campus Box 8205,
More informationKernel Methods. Outline
Kernel Methods Quang Nguyen University of Pittsburgh CS 3750, Fall 2011 Outline Motivation Examples Kernels Definitions Kernel trick Basic properties Mercer condition Constructing feature space Hilbert
More informationSpace-Frequency Atoms
Space-Frequency Atoms FREQUENCY FREQUENCY SPACE SPACE FREQUENCY FREQUENCY SPACE SPACE Figure 1: Space-frequency atoms. Windowed Fourier Transform 1 line 1 0.8 0.6 0.4 0.2 0-0.2-0.4-0.6-0.8-1 0 100 200
More informationKernel metrics on normal cycles for the matching of geometrical structures.
Kernel metrics on normal cycles for the matching of geometrical structures. Pierre Roussillon University Paris Descartes July 18, 2016 1 / 34 Overview Introduction Matching of geometrical structures Kernel
More informationWAVELETS WITH COMPOSITE DILATIONS
ELECTRONIC RESEARCH ANNOUNCEMENTS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Pages 000 000 (Xxxx XX, XXXX S 1079-6762(XX0000-0 WAVELETS WITH COMPOSITE DILATIONS KANGHUI GUO, DEMETRIO LABATE, WANG-Q
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 informationA new class of pseudodifferential operators with mixed homogenities
A new class of pseudodifferential operators with mixed homogenities Po-Lam Yung University of Oxford Jan 20, 2014 Introduction Given a smooth distribution of hyperplanes on R N (or more generally on a
More informationApplied and Computational Harmonic Analysis 11, (2001) doi: /acha , available online at
Applied and Computational Harmonic Analysis 11 305 31 (001 doi:10.1006/acha.001.0355 available online at http://www.idealibrary.com on LETTER TO THE EDITOR Construction of Multivariate Tight Frames via
More informationOptimal series representations of continuous Gaussian random fields
Optimal series representations of continuous Gaussian random fields Antoine AYACHE Université Lille 1 - Laboratoire Paul Painlevé A. Ayache (Lille 1) Optimality of continuous Gaussian series 04/25/2012
More informationFrames and operator representations of frames
Frames and operator representations of frames Ole Christensen Joint work with Marzieh Hasannasab HATA DTU DTU Compute, Technical University of Denmark HATA: Harmonic Analysis - Theory and Applications
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More information