Image Analysis and Approximation via Generalized Polyharmonic Local Trigonometric Transforms
|
|
- Mark Holt
- 5 years ago
- Views:
Transcription
1 Image Analysis and Approximation via Generalized Polyharmonic Local Trigonometric Transforms Naoki Saito with Jean-François Remy, Noel Smith, Ernest Woei, Katsu Yamatani, & Jim Zhao saito/ Department of Mathematics University of California, Davis Partially supported by ONR N , PECASE, and NSF DMS MGAWS1: Sep. 21, 2004 p.1
2 Outline Motivations: Issues of Localized Fourier Analysis Polyharmonic Local Sine Transform (PHLST) for Rectangular Domains Extension to 3D Extension to General Shape Domain Summary and Future Plan MGAWS1: Sep. 21, 2004 p.2
3 Motivations Want sparse data representation (compression) = Fast decay of the expansion coefficients Want a data representation as less statistically dependent among blocks as possible = No overlaps among blocks Want to develop a local signal analysis tool that can distinguish intrinsic singularities from the artificial discontinuities created by local windowing = many potential applications Want to efficiently represent regions of more general shapes other than rectangular blocks MGAWS1: Sep. 21, 2004 p.3
4 Life on an Interval (or a Domain) If data are periodic and smooth, the Fourier basis is very efficient: Fast decay of the Fourier coefficients Smoothness estimate (e.g., the Lipschitz/Hölder exponents) However, most signals and images of interest have compact supports, and are neither periodic nor smooth. What to do? MGAWS1: Sep. 21, 2004 p.4
5 Life on an Interval (or a Domain)... Design special wavelets ( wavelets on intervals ); Use multiwavelets of Alpert and Rokhlin (segmented orthogonal polynomials); Use Prolate Spheroidal Wave Functions ; Use the Continuous Boundary Local Fourier Transform (CBLFT) that periodizes the signal nicely and expand it into a periodic basis; or Use the Polyharmonic Local Trigonometric Transform (PHLTT) and its generalizations MGAWS1: Sep. 21, 2004 p.5
6 Basic Ideas of PHLTT Split the domain Ω into a set of subregions (often rectangular blocks) {Ω j } j and cut the data/function f into pieces using the characteristic functions, f j = χ Ωj f MGAWS1: Sep. 21, 2004 p.6
7 Basic Ideas of PHLTT Split the domain Ω into a set of subregions (often rectangular blocks) {Ω j } j and cut the data/function f into pieces using the characteristic functions, f j = χ Ωj f f j is decomposed into two parts u j + v j MGAWS1: Sep. 21, 2004 p.6
8 Basic Ideas of PHLTT Split the domain Ω into a set of subregions (often rectangular blocks) {Ω j } j and cut the data/function f into pieces using the characteristic functions, f j = χ Ωj f f j is decomposed into two parts u j + v j u j is a solution of the polyharmonic equations m u = 0 with Dirichlet/Neumann boundary conditions MGAWS1: Sep. 21, 2004 p.6
9 Basic Ideas of PHLTT Split the domain Ω into a set of subregions (often rectangular blocks) {Ω j } j and cut the data/function f into pieces using the characteristic functions, f j = χ Ωj f f j is decomposed into two parts u j + v j u j is a solution of the polyharmonic equations m u = 0 with Dirichlet/Neumann boundary conditions v j is a residual f j u j, which are expanded into the multiple Fourier series (either complex exponentials, sines, or cosines) MGAWS1: Sep. 21, 2004 p.6
10 Basic Ideas of PHLTT Split the domain Ω into a set of subregions (often rectangular blocks) {Ω j } j and cut the data/function f into pieces using the characteristic functions, f j = χ Ωj f f j is decomposed into two parts u j + v j u j is a solution of the polyharmonic equations m u = 0 with Dirichlet/Neumann boundary conditions v j is a residual f j u j, which are expanded into the multiple Fourier series (either complex exponentials, sines, or cosines) Sines and cosines are eigenfunctions of the Laplacian on a rectangular box MGAWS1: Sep. 21, 2004 p.6
11 1D Case In 1D, each u j is a low order algebraic polynomial (e.g., line, cubic poly.) MGAWS1: Sep. 21, 2004 p.7
12 1D Case In 1D, each u j is a low order algebraic polynomial (e.g., line, cubic poly.) v j is a trigonometric polynomial MGAWS1: Sep. 21, 2004 p.7
13 1D Case In 1D, each u j is a low order algebraic polynomial (e.g., line, cubic poly.) v j is a trigonometric polynomial These two compensate the shortcomings of each other: MGAWS1: Sep. 21, 2004 p.7
14 1D Case In 1D, each u j is a low order algebraic polynomial (e.g., line, cubic poly.) v j is a trigonometric polynomial These two compensate the shortcomings of each other: High order algebraic polynomial = Runge phenomenon MGAWS1: Sep. 21, 2004 p.7
15 1D Case In 1D, each u j is a low order algebraic polynomial (e.g., line, cubic poly.) v j is a trigonometric polynomial These two compensate the shortcomings of each other: High order algebraic polynomial = Runge phenomenon Trigonometric polynomial on an interval = Gibbs phenomenon MGAWS1: Sep. 21, 2004 p.7
16 1D Example: Compression Ratio MGAWS1: Sep. 21, 2004 p.8
17 Polyharmonic Global Sine Transform on a Rectangle Consider a function f C 2m (Ω), where m = 1, 2,..., and Ω R n (e.g., Ω = [0, 1] n ), but not periodic. MGAWS1: Sep. 21, 2004 p.9
18 Polyharmonic Global Sine Transform on a Rectangle Consider a function f C 2m (Ω), where m = 1, 2,..., and Ω R n (e.g., Ω = [0, 1] n ), but not periodic. Decompose this function into the following two components: f(x) = u(x) + v(x), MGAWS1: Sep. 21, 2004 p.9
19 Polyharmonic Global Sine Transform on a Rectangle Consider a function f C 2m (Ω), where m = 1, 2,..., and Ω R n (e.g., Ω = [0, 1] n ), but not periodic. Decompose this function into the following two components: f(x) = u(x) + v(x), u(x) satisfies the following polyharmonic equation: m u = 0 in Ω. MGAWS1: Sep. 21, 2004 p.9
20 PHGST on a Rectangle... The boundary condition for u is: p l u ν p l = p l f ν p l on Γ = Ω, l = 0,..., m 1, where p l is the order of the normal derivatives to be specified (p 0 0 u = f on Γ). MGAWS1: Sep. 21, 2004 p.10
21 PHGST on a Rectangle... The boundary condition for u is: p l u ν p l = p l f ν p l on Γ = Ω, l = 0,..., m 1, where p l is the order of the normal derivatives to be specified (p 0 0 u = f on Γ). Now set v(x) = f(x) u(x), which we will call the residual component with p l v ν p l = 0 on Γ, l = 0,..., m 1. MGAWS1: Sep. 21, 2004 p.10
22 A Specific Example: m = 1 (Laplace) Case u = 0 in Ω, u = f on Γ. Variational formulation = minimum gradient interpolation: u 2 dx subject to the above boundary condition. min u H 1 (Ω) Ω Note that in 1D, this is simply a line. MGAWS1: Sep. 21, 2004 p.11
23 A Specific Example: m = 2 (Biharmonic) Case 2 u = 0 in Ω, 2 u u = f, ν = 2 f on Γ. 2 ν 2 Variational formulation = minimum curvature interpolation: min u H 2 (Ω) Ω ( u + 2 j k subject to the above boundary condition. j k u) 2 dx, Note that in 1D, this is simply a cubic polynomial. MGAWS1: Sep. 21, 2004 p.12
24 PHGST on a Rectangle... The polyharmonic component u is smooth inside the domain Ω, and its values (and possibly its normal derivatives) matches those of data. MGAWS1: Sep. 21, 2004 p.13
25 PHGST on a Rectangle... The polyharmonic component u is smooth inside the domain Ω, and its values (and possibly its normal derivatives) matches those of data. The u component can be represented only by the boundary values f Γ. = No need to store the whole u. MGAWS1: Sep. 21, 2004 p.13
26 PHGST on a Rectangle... The polyharmonic component u is smooth inside the domain Ω, and its values (and possibly its normal derivatives) matches those of data. The u component can be represented only by the boundary values f Γ. = No need to store the whole u. The residual component v becomes 0 at the boundary Γ. MGAWS1: Sep. 21, 2004 p.13
27 PHGST on a Rectangle... Therefore, the v component is suitable for Fourier analysis. In fact, if Ω = [0, 1] n and p l = 2l, l = 0,..., m 1, then the Fourier sine analysis should be used to get the matching normal derivatives up to order 2m 1 by odd reflection at the boundaries. MGAWS1: Sep. 21, 2004 p.14
28 PHGST on a Rectangle... Therefore, the v component is suitable for Fourier analysis. In fact, if Ω = [0, 1] n and p l = 2l, l = 0,..., m 1, then the Fourier sine analysis should be used to get the matching normal derivatives up to order 2m 1 by odd reflection at the boundaries. The frequency content (in particular, mid to high frequency range) of the original is retained in the residual = textures remain in v; shading is captured by u. MGAWS1: Sep. 21, 2004 p.14
29 PHGST on a Rectangle... Therefore, the v component is suitable for Fourier analysis. In fact, if Ω = [0, 1] n and p l = 2l, l = 0,..., m 1, then the Fourier sine analysis should be used to get the matching normal derivatives up to order 2m 1 by odd reflection at the boundaries. The frequency content (in particular, mid to high frequency range) of the original is retained in the residual = textures remain in v; shading is captured by u. We can get the decay rate as follows: MGAWS1: Sep. 21, 2004 p.14
30 Theorem Let Ω = [0, 1] n, and f C 2m (Ω), but non-periodic. Assume further that i 2m+1 f, i = 1,..., n, exist and are of bounded variation. Furthermore, let f = u + v be the PHLST representation where the polyharmonic component u is the solution of the polyharmonic equation of order m with the boundary condition 2l u ν 2l = 2l f ν 2l on Γ, l = 0,..., m 1. Then, the Fourier sine coefficient b k of the residual v is of O ( k 2m 1) for all k 0, where k = (k 1,..., k n ), and k is the usual Euclidean (i.e., l 2 ) norm of k. MGAWS1: Sep. 21, 2004 p.15
31 Polyharmonic Local Sine Transform Now, consider a decomposition of Ω into a disjoint set of subdomains {Ω j }, i.e., Ω = J j=1ω j. A typical example is Ω = (0, 1) n, and Ω j is a dyadic subcube. Then, restrict f on Ω j, i.e., for each j, we decompose f locally as follows: fχ Ωj = f j = u j + v j, where we follow the same recipe locally as in the global case. We call this decomposition of f into {u j, v j } Polyharmonic Local Sine Transform (PHLST). For m = 1: Laplace Local Sine Transform (LLST); For m = 2: Biharmonic Local Sine Transform (BLST). MGAWS1: Sep. 21, 2004 p.16
32 Polyharmonic Local Sine Transform (a) Original MGAWS1: Sep. 21, 2004 p.17
33 Polyharmonic Local Sine Transform (a) Original (b) j u j MGAWS1: Sep. 21, 2004 p.17
34 Polyharmonic Local Sine Transform (a) Original (b) j u j (c) j v j MGAWS1: Sep. 21, 2004 p.17
35 Polyharmonic Local Sine Transform... No spatial overlaps MGAWS1: Sep. 21, 2004 p.18
36 Polyharmonic Local Sine Transform... No spatial overlaps Decay of the Fourier sine coefficients are fast if Ω j does not contain any singularity MGAWS1: Sep. 21, 2004 p.18
37 Polyharmonic Local Sine Transform... No spatial overlaps Decay of the Fourier sine coefficients are fast if Ω j does not contain any singularity Can distinguish intrinsic singularities from the artificial discontinuities at Γ j imposed by local windowing MGAWS1: Sep. 21, 2004 p.18
38 Remarks on the Laplace Solver A fast and accurate (no finite difference approximation of the Laplace operator) algorithm (for both 2D and 3D) exists via FFT if Ω and Ω j are dyadic cubes (Averbuch, Braverman, Israeli, and Vozovoi, 1998) MGAWS1: Sep. 21, 2004 p.19
39 Remarks on the Laplace Solver A fast and accurate (no finite difference approximation of the Laplace operator) algorithm (for both 2D and 3D) exists via FFT if Ω and Ω j are dyadic cubes (Averbuch, Braverman, Israeli, and Vozovoi, 1998) The u component can be computed as: u(x, y) = p(x, y) + k 1 where... {b (1) k h k(x, 1 y) + b (2) k h k(y, 1 x) + b (3) k h k(x, y) + b (4) k h k(y, x)}, MGAWS1: Sep. 21, 2004 p.19
40 Remarks on the Laplace Solver p(x, y) is a harmonic polynomial that agrees with f(x, y) at the four corner points of the domain, e.g., p(x, y) = a 3 xy + a 2 x + a 1 y + a 0, MGAWS1: Sep. 21, 2004 p.20
41 Remarks on the Laplace Solver p(x, y) is a harmonic polynomial that agrees with f(x, y) at the four corner points of the domain, e.g., h k (x, y) is defined as: p(x, y) = a 3 xy + a 2 x + a 1 y + a 0, h k (x, y) = sin(πkx) sinh(πky) sinh(πk), where b (j) k are the kth 1D Fourier sine coefficients of the boundary data of f p. MGAWS1: Sep. 21, 2004 p.20
42 Remarks on the Laplace Solver p(x, y) is a harmonic polynomial that agrees with f(x, y) at the four corner points of the domain, e.g., h k (x, y) is defined as: p(x, y) = a 3 xy + a 2 x + a 1 y + a 0, h k (x, y) = sin(πkx) sinh(πky) sinh(πk), where b (j) k are the kth 1D Fourier sine coefficients of the boundary data of f p. The computational cost of u is O(4N 2 log N). MGAWS1: Sep. 21, 2004 p.20
43 Remarks on PHLST Need to store the boundary values = can compress them using the lower dimensional version of PHLST MGAWS1: Sep. 21, 2004 p.21
44 Remarks on PHLST Need to store the boundary values = can compress them using the lower dimensional version of PHLST Can use complex exponentials, wavelets, etc., instead of sines with potentially slower decay MGAWS1: Sep. 21, 2004 p.21
45 Remarks on PHLST Need to store the boundary values = can compress them using the lower dimensional version of PHLST Can use complex exponentials, wavelets, etc., instead of sines with potentially slower decay Can do in the frequency domain = better wavelet packets MGAWS1: Sep. 21, 2004 p.21
46 Remarks on PHLST Need to store the boundary values = can compress them using the lower dimensional version of PHLST Can use complex exponentials, wavelets, etc., instead of sines with potentially slower decay Can do in the frequency domain = better wavelet packets Can be generalized to other geometries (e.g., discs, spheres, star shapes, general smooth boundaries) MGAWS1: Sep. 21, 2004 p.21
47 Remarks on PHLST Need to store the boundary values = can compress them using the lower dimensional version of PHLST Can use complex exponentials, wavelets, etc., instead of sines with potentially slower decay Can do in the frequency domain = better wavelet packets Can be generalized to other geometries (e.g., discs, spheres, star shapes, general smooth boundaries) Useful for interpolation and local feature computation (e.g., gradients, directional derivatives, etc.) MGAWS1: Sep. 21, 2004 p.21
48 f = u + v Models: Interpretation of PHLTT f = u + v models of Yves Meyer (2001): MGAWS1: Sep. 21, 2004 p.22
49 f = u + v Models: Interpretation of PHLTT f = u + v models of Yves Meyer (2001): u objects (e.g., BV (Ω)) MGAWS1: Sep. 21, 2004 p.22
50 f = u + v Models: Interpretation of PHLTT f = u + v models of Yves Meyer (2001): u objects (e.g., BV (Ω)) v textures and noise MGAWS1: Sep. 21, 2004 p.22
51 f = u + v Models: Interpretation of PHLTT f = u + v models of Yves Meyer (2001): u objects (e.g., BV (Ω)) v textures and noise In PHLTT, the polyharmonic component u is what is predictable from only using the boundary behavior of a function. MGAWS1: Sep. 21, 2004 p.22
52 f = u + v Models: Interpretation of PHLTT f = u + v models of Yves Meyer (2001): u objects (e.g., BV (Ω)) v textures and noise In PHLTT, the polyharmonic component u is what is predictable from only using the boundary behavior of a function. The residual component v is what is unpredictable using the boundary information. MGAWS1: Sep. 21, 2004 p.22
53 f = u + v Models: Interpretation of PHLTT f = u + v models of Yves Meyer (2001): u objects (e.g., BV (Ω)) v textures and noise In PHLTT, the polyharmonic component u is what is predictable from only using the boundary behavior of a function. The residual component v is what is unpredictable using the boundary information. Yet another interpretation is: u: trend, v: fluctuation. MGAWS1: Sep. 21, 2004 p.22
54 f = u + v Models: Interpretation of PHLTT... Many possible ideas here once we view this as predicition from the available data: MGAWS1: Sep. 21, 2004 p.23
55 f = u + v Models: Interpretation of PHLTT... Many possible ideas here once we view this as predicition from the available data: A. Cohen et al. uses inpainting algorithm for predicting u components MGAWS1: Sep. 21, 2004 p.23
56 f = u + v Models: Interpretation of PHLTT... Many possible ideas here once we view this as predicition from the available data: A. Cohen et al. uses inpainting algorithm for predicting u components Relationship with geometric harmonic analysis of Coifman & Lafon MGAWS1: Sep. 21, 2004 p.23
57 f = u + v Models: Interpretation of PHLTT... Many possible ideas here once we view this as predicition from the available data: A. Cohen et al. uses inpainting algorithm for predicting u components Relationship with geometric harmonic analysis of Coifman & Lafon Other inverse boundary problems... MGAWS1: Sep. 21, 2004 p.23
58 LLST Approximation Experiment Comparison with DCT, the fixed folding local cosine (LCT), and Wavelet (Coiflet 12). MGAWS1: Sep. 21, 2004 p.24
59 LLST Approximation Experiment Comparison with DCT, the fixed folding local cosine (LCT), and Wavelet (Coiflet 12). Split into a set of squares of pixels or pixels MGAWS1: Sep. 21, 2004 p.24
60 LLST Approximation Experiment Comparison with DCT, the fixed folding local cosine (LCT), and Wavelet (Coiflet 12). Split into a set of squares of pixels or pixels Select the largest 1% and 5% of coefficients. MGAWS1: Sep. 21, 2004 p.24
61 LLST Approximation Experiment Comparison with DCT, the fixed folding local cosine (LCT), and Wavelet (Coiflet 12). Split into a set of squares of pixels or pixels Select the largest 1% and 5% of coefficients. Compression using quantization is in progress. MGAWS1: Sep. 21, 2004 p.24
62 LLST Approximation Experiment: Top 1% (a) Original (b) DCT (c) LCT (d) C12 (e) LLST MGAWS1: Sep. 21, 2004 p.25
63 Zoomed up images: Top 1% (a) Original (b) DCT (c) LCT (d) C12 (e) LLST MGAWS1: Sep. 21, 2004 p.26
64 LLST Approximation Experiment: Top 5% (a) Original (b) DCT (c) LCT (d) C12 (e) LLST MGAWS1: Sep. 21, 2004 p.27
65 Zoomed up images: Top 5% (a) Original (b) DCT (c) LCT (d) C12 (e) LLST MGAWS1: Sep. 21, 2004 p.28
66 Approximation Test: Smooth Function PSNR (db) LLST: J=1 DCT: J=1 LCT: J=1 Daub6: J=3 Symm8: J= (a) Original Ratio of retained coefficients (%) (b) PSNR MGAWS1: Sep. 21, 2004 p.29
67 Approximation Test: Piecewise Smooth Function PSNR (db) LLST: J=1 DCT: J=1 LCT: J=1 Daub6: J=3 Symm8: J= (a) Original Ratio of retained coefficients (%) (b) PSNR MGAWS1: Sep. 21, 2004 p.30
68 Approximation Test: Oscillatory Function LLST: J=1 DCT: J=1 LCT: J=1 Daub6: J=3 Symm8: J= PSNR (db) (a) Original Ratio of retained coefficients (%) (b) PSNR MGAWS1: Sep. 21, 2004 p.31
69 Approximation Test: Oscillatory Function with Discontinuity LLST: J=1 DCT: J=1 LCT: J=1 Daub6: J=3 Symm8: J= PSNR (db) (a) Original Ratio of retained coefficients (%) (b) PSNR MGAWS1: Sep. 21, 2004 p.32
70 Extension to 3D Boundary data contains more information compared to 1D and 2D MGAWS1: Sep. 21, 2004 p.33
71 Extension to 3D Boundary data contains more information compared to 1D and 2D Fast 3D Laplace solver exists (the ABIV algorithm with our new organization of data): O(12N 3 log N) MGAWS1: Sep. 21, 2004 p.33
72 Extension to 3D Boundary data contains more information compared to 1D and 2D Fast 3D Laplace solver exists (the ABIV algorithm with our new organization of data): O(12N 3 log N) Nicely recursive: corners; edges; faces; body contents MGAWS1: Sep. 21, 2004 p.33
73 Extension to 3D Boundary data contains more information compared to 1D and 2D Fast 3D Laplace solver exists (the ABIV algorithm with our new organization of data): O(12N 3 log N) Nicely recursive: corners; edges; faces; body contents Many possibilities: hierarchical decomposition, parallelization, 3D feature extraction,... MGAWS1: Sep. 21, 2004 p.33
74 Extension to 3D Boundary data contains more information compared to 1D and 2D Fast 3D Laplace solver exists (the ABIV algorithm with our new organization of data): O(12N 3 log N) Nicely recursive: corners; edges; faces; body contents Many possibilities: hierarchical decomposition, parallelization, 3D feature extraction,... Promising applications include analysis and compression of 3D seismic, medical, and video data,... MGAWS1: Sep. 21, 2004 p.33
75 Extension to General Shape Domains Consider a general-shaped domain Ω R n. MGAWS1: Sep. 21, 2004 p.34
76 Extension to General Shape Domains Consider a general-shaped domain Ω R n. Want to analyze the frequency information inside of the object. MGAWS1: Sep. 21, 2004 p.34
77 Extension to General Shape Domains Consider a general-shaped domain Ω R n. Want to analyze the frequency information inside of the object. Want to avoid the Gibbs phenomenon due to Γ. MGAWS1: Sep. 21, 2004 p.34
78 Extension to General Shape Domains Consider a general-shaped domain Ω R n. Want to analyze the frequency information inside of the object. Want to avoid the Gibbs phenomenon due to Γ. = Smoothly extend the function defined on Ω to the outside! MGAWS1: Sep. 21, 2004 p.34
79 (c) Object (d) Anomalies MGAWS1: Sep. 21, 2004 p.35 Object-Oriented Image Analysis (a) Original (b) Background
80 Object-Oriented Image Analysis... Suppose its boundary Γ is already detected explicitly. MGAWS1: Sep. 21, 2004 p.36
81 Object-Oriented Image Analysis... Suppose its boundary Γ is already detected explicitly. Let f C 2 (Ω) C(Ω). MGAWS1: Sep. 21, 2004 p.36
82 Object-Oriented Image Analysis... Suppose its boundary Γ is already detected explicitly. Let f C 2 (Ω) C(Ω). Let S be a rectangular domain containing Ω. MGAWS1: Sep. 21, 2004 p.36
83 Object-Oriented Image Analysis... Suppose its boundary Γ is already detected explicitly. Let f C 2 (Ω) C(Ω). Let S be a rectangular domain containing Ω. Then, we will find a smooth extension ṽ to the entire S from Ω and the harmonic function u in Ω: f in Ω, ṽ = 0 in S\Ω, ṽ = 0 on S. MGAWS1: Sep. 21, 2004 p.36
84 Motivations for ṽ ṽ is essentially the same as f in Ω modulo some harmonic function (in fact, this harmonic function becomes the u component); MGAWS1: Sep. 21, 2004 p.37
85 Motivations for ṽ ṽ is essentially the same as f in Ω modulo some harmonic function (in fact, this harmonic function becomes the u component); ṽ is harmonic (i.e., very smooth) in S\Ω (outside of Ω); MGAWS1: Sep. 21, 2004 p.37
86 Motivations for ṽ ṽ is essentially the same as f in Ω modulo some harmonic function (in fact, this harmonic function becomes the u component); ṽ is harmonic (i.e., very smooth) in S\Ω (outside of Ω); can show ṽ C 1 (S) MGAWS1: Sep. 21, 2004 p.37
87 Motivations for ṽ ṽ is essentially the same as f in Ω modulo some harmonic function (in fact, this harmonic function becomes the u component); ṽ is harmonic (i.e., very smooth) in S\Ω (outside of Ω); can show ṽ C 1 (S) can also show ṽ S\Γ C 2 (S\Γ) MGAWS1: Sep. 21, 2004 p.37
88 Motivations for ṽ ṽ is essentially the same as f in Ω modulo some harmonic function (in fact, this harmonic function becomes the u component); ṽ is harmonic (i.e., very smooth) in S\Ω (outside of Ω); can show ṽ C 1 (S) can also show ṽ S\Γ C 2 (S\Γ) ṽ S = 0 = suitable for Fourier sine series expansion, with decay rate O( k 3 ). MGAWS1: Sep. 21, 2004 p.37
89 How to compute ṽ Introduce the potential function: P (x) = f(y)φ(x, y) dy for x R n. Ω where Φ(x, y) = 1 2π log x y if n = 2, x y 2 n (2 n)ω n if n > 2, is the so-called fundamental solution of the Laplace equation. MGAWS1: Sep. 21, 2004 p.38
90 How to compute ṽ... Can show that P satisfies the above Poisson equation: f in Ω, P = 0 in S\Ω, MGAWS1: Sep. 21, 2004 p.39
91 How to compute ṽ... Can show that P satisfies the above Poisson equation: f in Ω, P = 0 in S\Ω, But the boundary condition at the outside box is not satisfied, i.e., P S 0. MGAWS1: Sep. 21, 2004 p.39
92 How to compute ṽ... Introduce now a boundary correction function Q: Q = 0 in S, Q = P on S; MGAWS1: Sep. 21, 2004 p.40
93 How to compute ṽ... Introduce now a boundary correction function Q: Q = 0 in S, Q = P on S; Finally set ṽ(x) = P (x) + Q(x), which satisfies the desired conditions. MGAWS1: Sep. 21, 2004 p.40
94 How to compute ṽ... Can write the solution via single layer and double layer potentials as: P (x) = χ Ω (x)f(x)+ Γ ( f ν ) y) (y)φ(x, y) f(y) Φ(x, ν y ds(y), where χ Ω (x) = 1 if x Ω; 1/2 if x Γ; 0 if x R n \Ω. MGAWS1: Sep. 21, 2004 p.41
95 Object-Oriented Image Analysis... The u component satisfies: u = 0 in Ω, u = f ṽ on Γ. MGAWS1: Sep. 21, 2004 p.42
96 Object-Oriented Image Analysis... The u component satisfies: u = 0 in Ω, u = f ṽ on Γ. u can be recovered from f Γ, f, and ṽ via Fast ν Γ Multipole Method (just potential evaluation, no need to solve the Laplace equation). MGAWS1: Sep. 21, 2004 p.42
97 1D Example 4 3 f=u+χ I v f u v x MGAWS1: Sep. 21, 2004 p.43
98 2D Example (a) Orig.1 MGAWS1: Sep. 21, 2004 p.44
99 2D Example (a) Orig.1 (b) Orig.2 MGAWS1: Sep. 21, 2004 p.44
100 2D Example (a) Orig.1 (b) Orig (c) Ext.1 MGAWS1: Sep. 21, 2004 p.44
101 2D Example (a) Orig.1 (b) Orig (c) Ext.1 (d) Ext.2 MGAWS1: Sep. 21, 2004 p.44
102 2D Example (a) Original MGAWS1: Sep. 21, 2004 p.45
103 2D Example (a) Original (b) u MGAWS1: Sep. 21, 2004 p.45
104 2D Example (a) Original (b) u (c) ṽ MGAWS1: Sep. 21, 2004 p.45
105 The Fourier Transform Magnitudes (log-scale) (a) Original MGAWS1: Sep. 21, 2004 p.46
106 The Fourier Transform Magnitudes (log-scale) (a) Original (b) ṽ MGAWS1: Sep. 21, 2004 p.46
107 Another example MGAWS1: Sep. 21, 2004 p.47
108 Another example (a) Extracted shell (b) Extended (c) u component (d) ṽ component MGAWS1: Sep. 21, 2004 p.48
109 Related Previous Work Whitney (1934): Smooth extensions of C m (Ω) functions in R n, and the Whitney decomposition Zygmund (1935): Removable discontinuities Lanczos (1938,1966): Trigonometric interpolation with linear/polynomial component removal Kantrovich & Krylov (1958): Rapidly convergent trigonometric series by polynomial removal Briggs (1974): Minimum curvature surface interpolation Grimson (1981): Visual surface reconstruction Terzopoulos (1983): Multilevel visual surface reconstruction MGAWS1: Sep. 21, 2004 p.49
110 Related Previous Work... Leclerc (1989): Image segmentation using MDL and local polynomials Wahba (1990): Variational formulation/multivariate splines Dimitrov (1996): Polyharmonic functions for quadratures Kounchev ( ): Theory of polysplines Casas & Torres (1996): Two-stage coding Caselles, Morel, & Sbert (1998): Axiomatic interpolation All of them except Lanczos mainly focus on the u-component. Lanczos explored neither higher dimensions nor multiscale setting. MGAWS1: Sep. 21, 2004 p.50
111 Conclusions PHLTT and its generalizations: Provide a compact image/data representation scheme MGAWS1: Sep. 21, 2004 p.51
112 Conclusions PHLTT and its generalizations: Provide a compact image/data representation scheme Proactive use of boundary information MGAWS1: Sep. 21, 2004 p.51
113 Conclusions PHLTT and its generalizations: Provide a compact image/data representation scheme Proactive use of boundary information Can get faster decay expansion coefficients MGAWS1: Sep. 21, 2004 p.51
114 Conclusions PHLTT and its generalizations: Provide a compact image/data representation scheme Proactive use of boundary information Can get faster decay expansion coefficients Allow object-oriented image analysis & synthesis MGAWS1: Sep. 21, 2004 p.51
115 Conclusions PHLTT and its generalizations: Provide a compact image/data representation scheme Proactive use of boundary information Can get faster decay expansion coefficients Allow object-oriented image analysis & synthesis Variety of applications: segmentation, compression, interpolation, local feature computation, 3D,... MGAWS1: Sep. 21, 2004 p.51
116 Comments This general area is a meeting ground of image processing, harmonic analysis, PDEs, and shape optimization MGAWS1: Sep. 21, 2004 p.52
117 Comments This general area is a meeting ground of image processing, harmonic analysis, PDEs, and shape optimization A large possibility for prediction operators from boundary MGAWS1: Sep. 21, 2004 p.52
118 Comments This general area is a meeting ground of image processing, harmonic analysis, PDEs, and shape optimization A large possibility for prediction operators from boundary Anisotropy and many interesting issues MGAWS1: Sep. 21, 2004 p.52
119 Comments This general area is a meeting ground of image processing, harmonic analysis, PDEs, and shape optimization A large possibility for prediction operators from boundary Anisotropy and many interesting issues Reliability of boundary information and noise is an issue MGAWS1: Sep. 21, 2004 p.52
120 References Check saito/publications/ periodically for more preprints on this work. For preliminary work, see: N. Saito & J.-F. Remy: A new local sine transform without overlaps: A combination of computational harmonic analysis and PDE, in Wavelets X (M. A. Unser, A. Aldroubi, & A. F. Laine, eds.), Proc. SPIE 5207, pp , See also references therein. Thank you very much for your attention! MGAWS1: Sep. 21, 2004 p.53
The Polyharmonic Local Sine Transform: A New Tool for Local Image Analysis and Synthesis without Edge Effect
The Polyharmonic Local Sine Transform: A New Tool for Local Image Analysis and Synthesis without Edge Effect Naoki Saito Jean-François Remy Department of Mathematics University of California Davis, CA
More informationImprovement of DCT-based Compression Algorithms Using Poisson s Equation
1 Improvement of DCT-based Compression Algorithms Using Poisson s Equation Katsu Yamatani* and aoki Saito, Senior Member, IEEE Abstract We propose two new image compression-decompression methods that reproduce
More informationHarmonic Wavelet Transform and Image Approximation
Harmonic Wavelet Transform and Image Approximation Zhihua Zhang and Naoki Saito Dept of Math, Univ of California, Davis, California, 95616, USA Email: zhangzh@mathucdavisedu saito@mathucdavisedu Abstract
More informationSIMULTANEOUS SEGMENTATION, COMPRESSION, AND DENOISING OF SIGNALS USING POLYHARMONIC LOCAL SINE TRANSFORM AND MINIMUM DESCRIPTION LENGTH CRITERION
SIMULTANEOUS SEGMENTATION, COMPRESSION, AND DENOISING OF SIGNALS USING POLYHARMONIC LOCAL SINE TRANSFORM AND MINIMUM DESCRIPTION LENGTH CRITERION Naoki Saito and Ernest Woei Department of Mathematics University
More informationARTICLE IN PRESS Journal of Computational and Applied Mathematics ( )
Journal of Computational and Applied Mathematics ( ) Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam The polynomial
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
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 informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationGEOMETRIC HARMONICS AS A STATISTICAL IMAGE PROCESSING TOOL FOR IMAGES ON IRREGULARLY-SHAPED DOMAINS. Naoki Saito
GEOMETRIC HARMONICS AS A STATISTICAL IMAGE PROCESSING TOOL FOR IMAGES ON IRREGULARLY-SHAPED DOMAINS Naoki Saito Department of Mathematics University of California Davis, CA 95616 USA Email: saito@math.ucdavis.edu
More information1 Motivation. Lecturer: Naoki Saito Scribe: Lawrence Austria/Allen Xue. May 23, This lecture is based on [1] and [2].
MAT 80: Laplacian Eigenfunctions: Theory, Applications, and Computations Lecture 15: Use of Laplacian Eigenfunctions and Eigenvalues for Analyzing Data on a Domain of Complicated Shape Lecturer: Naoki
More informationApplications of Polyspline Wavelets to Astronomical Image Analysis
VIRTUAL OBSERVATORY: Plate Content Digitization, Archive Mining & Image Sequence Processing edited by M. Tsvetkov, V. Golev, F. Murtagh, and R. Molina, Heron Press, Sofia, 25 Applications of Polyspline
More informationUse of Laplacian Eigenfunctions and Eigenvalues for Analyzing Data on a Domain of Complicated Shape
Use of Laplacian Eigenfunctions and Eigenvalues for Analyzing Data on a Domain of Complicated Shape Naoki Saito 1 Department of Mathematics University of California, Davis May 23, 2007 1 Partially supported
More informationSeismic compression. François G. Meyer Department of Electrical Engineering University of Colorado at Boulder
Seismic compression François G. Meyer Department of Electrical Engineering University of Colorado at Boulder francois.meyer@colorado.edu http://ece-www.colorado.edu/ fmeyer IPAM, MGA 2004 Seismic Compression
More informationDiffusion Wavelets and Applications
Diffusion Wavelets and Applications J.C. Bremer, R.R. Coifman, P.W. Jones, S. Lafon, M. Mohlenkamp, MM, R. Schul, A.D. Szlam Demos, web pages and preprints available at: S.Lafon: www.math.yale.edu/~sl349
More informationA HIERARCHICAL 3-D DIRECT HELMHOLTZ SOLVER BY DOMAIN DECOMPOSITION AND MODIFIED FOURIER METHOD
SIAM J. SCI. COMPUT. Vol. 26, No. 5, pp. 1504 1524 c 2005 Society for Industrial and Applied Mathematics A HIERARCHICAL 3-D DIRECT HELMHOLTZ SOLVER BY DOMAIN DECOMPOSITION AND MODIFIED FOURIER METHOD E.
More informationSpectral Processing. Misha Kazhdan
Spectral Processing Misha Kazhdan [Taubin, 1995] A Signal Processing Approach to Fair Surface Design [Desbrun, et al., 1999] Implicit Fairing of Arbitrary Meshes [Vallet and Levy, 2008] Spectral Geometry
More informationJim Lambers ENERGY 281 Spring Quarter Lecture 5 Notes
Jim ambers ENERGY 28 Spring Quarter 27-8 ecture 5 Notes These notes are based on Rosalind Archer s PE28 lecture notes, with some revisions by Jim ambers. Fourier Series Recall that in ecture 2, when we
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 12 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted for noncommercial,
More information524 Jan-Olov Stromberg practice, this imposes strong restrictions both on N and on d. The current state of approximation theory is essentially useless
Doc. Math. J. DMV 523 Computation with Wavelets in Higher Dimensions Jan-Olov Stromberg 1 Abstract. In dimension d, a lattice grid of size N has N d points. The representation of a function by, for instance,
More informationWavelets and Image Compression Augusta State University April, 27, Joe Lakey. Department of Mathematical Sciences. New Mexico State University
Wavelets and Image Compression Augusta State University April, 27, 6 Joe Lakey Department of Mathematical Sciences New Mexico State University 1 Signals and Images Goal Reduce image complexity with little
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 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 informationPolynomial Approximation: The Fourier System
Polynomial Approximation: The Fourier System Charles B. I. Chilaka CASA Seminar 17th October, 2007 Outline 1 Introduction and problem formulation 2 The continuous Fourier expansion 3 The discrete Fourier
More informationA NEW BASIS SELECTION PARADIGM FOR WAVELET PACKET IMAGE CODING
A NEW BASIS SELECTION PARADIGM FOR WAVELET PACKET IMAGE CODING Nasir M. Rajpoot, Roland G. Wilson, François G. Meyer, Ronald R. Coifman Corresponding Author: nasir@dcs.warwick.ac.uk ABSTRACT In this paper,
More informationA WAVELET BASED CODING SCHEME VIA ATOMIC APPROXIMATION AND ADAPTIVE SAMPLING OF THE LOWEST FREQUENCY BAND
A WAVELET BASED CODING SCHEME VIA ATOMIC APPROXIMATION AND ADAPTIVE SAMPLING OF THE LOWEST FREQUENCY BAND V. Bruni, D. Vitulano Istituto per le Applicazioni del Calcolo M. Picone, C. N. R. Viale del Policlinico
More informationFourier series on triangles and tetrahedra
Fourier series on triangles and tetrahedra Daan Huybrechs K.U.Leuven Cambridge, Sep 13, 2010 Outline Introduction The unit interval Triangles, tetrahedra and beyond Outline Introduction The topic of this
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 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 informationResearch Statement. James Bremer Department of Mathematics, University of California, Davis
Research Statement James Bremer Department of Mathematics, University of California, Davis Email: bremer@math.ucdavis.edu Webpage: https.math.ucdavis.edu/ bremer I work in the field of numerical analysis,
More informationMysteries Around Graph Laplacian Eigenvalue 4
Mysteries Around Graph Laplacian Eigenvalue 4 Yuji Nakatsukasa, Naoki Saito, Ernest Woei Department of Mathematics University of California, Davis ICIAM 2011, Vancouver, BC, Canada July 22, 2011 saito@math.ucdavis.edu
More informationIntroduction to Wavelets and Wavelet Transforms
Introduction to Wavelets and Wavelet Transforms A Primer C. Sidney Burrus, Ramesh A. Gopinath, and Haitao Guo with additional material and programs by Jan E. Odegard and Ivan W. Selesnick Electrical and
More informationRepresentation: Fractional Splines, Wavelets and related Basis Function Expansions. Felix Herrmann and Jonathan Kane, ERL-MIT
Representation: Fractional Splines, Wavelets and related Basis Function Expansions Felix Herrmann and Jonathan Kane, ERL-MIT Objective: Build a representation with a regularity that is consistent with
More informationNumerical Approximation Methods for Non-Uniform Fourier Data
Numerical Approximation Methods for Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu 2014 Joint Mathematics Meetings January 18 2014 0 / 16 Joint work with Anne Gelb (Arizona State) Guohui
More informationQuarkonial frames of wavelet type - Stability, approximation and compression properties
Quarkonial frames of wavelet type - Stability, approximation and compression properties Stephan Dahlke 1 Peter Oswald 2 Thorsten Raasch 3 ESI Workshop Wavelet methods in scientific computing Vienna, November
More informationAdapted Feature Extraction and Its Applications
saito@math.ucdavis.edu 1 Adapted Feature Extraction and Its Applications Naoki Saito Department of Mathematics University of California Davis, CA 95616 email: saito@math.ucdavis.edu URL: http://www.math.ucdavis.edu/
More informationInverse Gravimetry Problem
Inverse Gravimetry Problem Victor Isakov September 21, 2010 Department of M athematics and Statistics W ichita State U niversity W ichita, KS 67260 0033, U.S.A. e mail : victor.isakov@wichita.edu 1 Formulation.
More informationH 2 -matrices with adaptive bases
1 H 2 -matrices with adaptive bases Steffen Börm MPI für Mathematik in den Naturwissenschaften Inselstraße 22 26, 04103 Leipzig http://www.mis.mpg.de/ Problem 2 Goal: Treat certain large dense matrices
More informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationPhysics 6303 Lecture 11 September 24, LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation
Physics 6303 Lecture September 24, 208 LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation, l l l l l l. Consider problems that are no axisymmetric; i.e., the potential depends
More informationOn the Phase Transition Phenomenon of Graph Laplacian Eigenfunctions on Trees
On the Phase Transition Phenomenon of Graph Laplacian Eigenfunctions on Trees Naoki Saito and Ernest Woei Department of Mathematics University of California Davis, CA 9566 USA Email: saito@math.ucdavis.edu;
More informationDiffusion Geometries, Diffusion Wavelets and Harmonic Analysis of large data sets.
Diffusion Geometries, Diffusion Wavelets and Harmonic Analysis of large data sets. R.R. Coifman, S. Lafon, MM Mathematics Department Program of Applied Mathematics. Yale University Motivations The main
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationFast Algorithms for the Computation of Oscillatory Integrals
Fast Algorithms for the Computation of Oscillatory Integrals Emmanuel Candès California Institute of Technology EPSRC Symposium Capstone Conference Warwick Mathematics Institute, July 2009 Collaborators
More informationData Analysis and Representation on a General. Domain using Eigenfunctions of Laplacian
Data Analysis and Representation on a General Domain using Eigenfunctions of Laplacian Naoki Saito Department of Mathematics University of California Davis, CA 95616 USA Abstract We propose a new method
More informationConstructing Approximation Kernels for Non-Harmonic Fourier Data
Constructing Approximation Kernels for Non-Harmonic Fourier Data Aditya Viswanathan aditya.v@caltech.edu California Institute of Technology SIAM Annual Meeting 2013 July 10 2013 0 / 19 Joint work with
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationToward a Realization of Marr s Theory of Primal Sketches via Autocorrelation Wavelets: Image Representation using Multiscale Edge Information
Toward a Realization of Marr s Theory of Primal Sketches via Autocorrelation Wavelets: Image Representation using Multiscale Edge Information Naoki Saito 1 Department of Mathematics University of California,
More informationMultiscale Image Transforms
Multiscale Image Transforms Goal: Develop filter-based representations to decompose images into component parts, to extract features/structures of interest, and to attenuate noise. Motivation: extract
More informationNumerical Methods I Orthogonal Polynomials
Numerical Methods I Orthogonal Polynomials Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Nov. 4th and 11th, 2010 A. Donev (Courant Institute)
More informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationFILTERING IN THE FREQUENCY DOMAIN
1 FILTERING IN THE FREQUENCY DOMAIN Lecture 4 Spatial Vs Frequency domain 2 Spatial Domain (I) Normal image space Changes in pixel positions correspond to changes in the scene Distances in I correspond
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationFast and accurate methods for the discretization of singular integral operators given on surfaces
Fast and accurate methods for the discretization of singular integral operators given on surfaces James Bremer University of California, Davis March 15, 2018 This is joint work with Zydrunas Gimbutas (NIST
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationHighly Oscillatory Problems: Computation, Theory and Applications
Highly Oscillatory Problems: Computation, Theory and Applications Isaac Newton Institute for Mathematical Sciences Edited by 1 Rapid function approximation by modified Fourier series Daan Huybrechs and
More informationMarch 13, Paper: R.R. Coifman, S. Lafon, Diffusion maps ([Coifman06]) Seminar: Learning with Graphs, Prof. Hein, Saarland University
Kernels March 13, 2008 Paper: R.R. Coifman, S. Lafon, maps ([Coifman06]) Seminar: Learning with Graphs, Prof. Hein, Saarland University Kernels Figure: Example Application from [LafonWWW] meaningful geometric
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationClassical Fourier Analysis
Loukas Grafakos Classical Fourier Analysis Second Edition 4y Springer 1 IP Spaces and Interpolation 1 1.1 V and Weak IP 1 1.1.1 The Distribution Function 2 1.1.2 Convergence in Measure 5 1.1.3 A First
More informationWavelet Footprints: Theory, Algorithms, and Applications
1306 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 51, NO. 5, MAY 2003 Wavelet Footprints: Theory, Algorithms, and Applications Pier Luigi Dragotti, Member, IEEE, and Martin Vetterli, Fellow, IEEE Abstract
More informationError Analysis for H 1 Based Wavelet Interpolations
Error Analysis for H 1 Based Wavelet Interpolations Tony F. Chan Hao-Min Zhou Tie Zhou Abstract We rigorously study the error bound for the H 1 wavelet interpolation problem, which aims to recover missing
More informationSplines which are piecewise solutions of polyharmonic equation
Splines which are piecewise solutions of polyharmonic equation Ognyan Kounchev March 25, 2006 Abstract This paper appeared in Proceedings of the Conference Curves and Surfaces, Chamonix, 1993 1 Introduction
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 05 Image Processing Basics 13/02/04 http://www.ee.unlv.edu/~b1morris/ecg782/
More informationIntroduction to Fourier Analysis
Lecture Introduction to Fourier Analysis Jan 7, 2005 Lecturer: Nati Linial Notes: Atri Rudra & Ashish Sabharwal. ext he main text for the first part of this course would be. W. Körner, Fourier Analysis
More informationAn Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems
An Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems Justin Romberg Georgia Tech, School of ECE ENS Winter School January 9, 2012 Lyon, France Applied and Computational
More informationDECOUPLING INEQUALITIES IN HARMONIC ANALYSIS AND APPLICATIONS
DECOUPLING INEQUALITIES IN HARMONIC ANALYSIS AND APPLICATIONS 1 NOTATION AND STATEMENT S compact smooth hypersurface in R n with positive definite second fundamental form ( Examples: sphere S d 1, truncated
More informationMA3232 Summary 5. d y1 dy1. MATLAB has a number of built-in functions for solving stiff systems of ODEs. There are ode15s, ode23s, ode23t, ode23tb.
umerical solutions of higher order ODE We can convert a high order ODE into a system of first order ODEs and then apply RK method to solve it. Stiff ODEs Stiffness is a special problem that can arise in
More information( nonlinear constraints)
Wavelet Design & Applications Basic requirements: Admissibility (single constraint) Orthogonality ( nonlinear constraints) Sparse Representation Smooth functions well approx. by Fourier High-frequency
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 informationComputer Problems for Fourier Series and Transforms
Computer Problems for Fourier Series and Transforms 1. Square waves are frequently used in electronics and signal processing. An example is shown below. 1 π < x < 0 1 0 < x < π y(x) = 1 π < x < 2π... and
More informationWavelets in Scattering Calculations
Wavelets in Scattering Calculations W. P., Brian M. Kessler, Gerald L. Payne polyzou@uiowa.edu The University of Iowa Wavelets in Scattering Calculations p.1/43 What are Wavelets? Orthonormal basis functions.
More informationIPAM MGA Tutorial on Feature Extraction and Denoising: A Saga of u + v Models
IPAM MGA Tutorial on Feature Extraction and Denoising: A Saga of u + v Models Naoki Saito saito@math.ucdavis.edu http://www.math.ucdavis.edu/ saito/ Department of Mathematics University of California,
More informationPeriodogram of a sinusoid + spike Single high value is sum of cosine curves all in phase at time t 0 :
Periodogram of a sinusoid + spike Single high value is sum of cosine curves all in phase at time t 0 : X(t) = µ + Asin(ω 0 t)+ Δ δ ( t t 0 ) ±σ N =100 Δ =100 χ ( ω ) Raises the amplitude uniformly at all
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationImage Filtering, Edges and Image Representation
Image Filtering, Edges and Image Representation Capturing what s important Req reading: Chapter 7, 9 F&P Adelson, Simoncelli and Freeman (handout online) Opt reading: Horn 7 & 8 FP 8 February 19, 8 A nice
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationIntegration. Tuesday, December 3, 13
4 Integration 4.3 Riemann Sums and Definite Integrals Objectives n Understand the definition of a Riemann sum. n Evaluate a definite integral using properties of definite integrals. 3 Riemann Sums 4 Riemann
More informationTransform methods. and its inverse can be used to analyze certain time-dependent PDEs. f(x) sin(sxπ/(n + 1))
AMSC/CMSC 661 Scientific Computing II Spring 2010 Transforms and Wavelets Dianne P. O Leary c 2005,2010 Some motivations: Transform methods The Fourier transform Fv(ξ) = ˆv(ξ) = v(x)e ix ξ dx, R d and
More informationESTIMATES FOR MAXIMAL SINGULAR INTEGRALS
ESTIMATES FOR MAXIMAL SINGULAR INTEGRALS LOUKAS GRAFAKOS Abstract. It is shown that maximal truncations of nonconvolution L -bounded singular integral operators with kernels satisfying Hörmander s condition
More information1 Introduction to Wavelet Analysis
Jim Lambers ENERGY 281 Spring Quarter 2007-08 Lecture 9 Notes 1 Introduction to Wavelet Analysis Wavelets were developed in the 80 s and 90 s as an alternative to Fourier analysis of signals. Some of the
More informationHOSVD Based Image Processing Techniques
HOSVD Based Image Processing Techniques András Rövid, Imre J. Rudas, Szabolcs Sergyán Óbuda University John von Neumann Faculty of Informatics Bécsi út 96/b, 1034 Budapest Hungary rovid.andras@nik.uni-obuda.hu,
More informationPhotometric Stereo: Three recent contributions. Dipartimento di Matematica, La Sapienza
Photometric Stereo: Three recent contributions Dipartimento di Matematica, La Sapienza Jean-Denis DUROU IRIT, Toulouse Jean-Denis DUROU (IRIT, Toulouse) 17 December 2013 1 / 32 Outline 1 Shape-from-X techniques
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationWeighted Nonlocal Laplacian on Interpolation from Sparse Data
Noname manuscript No. (will be inserted by the editor) Weighted Nonlocal Laplacian on Interpolation from Sparse Data Zuoqiang Shi Stanley Osher Wei Zhu Received: date / Accepted: date Abstract Inspired
More informationINTRODUCTION TO. Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR)
INTRODUCTION TO WAVELETS Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR) CRITICISM OF FOURIER SPECTRUM It gives us the spectrum of the
More informationWhich wavelet bases are the best for image denoising?
Which wavelet bases are the best for image denoising? Florian Luisier a, Thierry Blu a, Brigitte Forster b and Michael Unser a a Biomedical Imaging Group (BIG), Ecole Polytechnique Fédérale de Lausanne
More informationA Fast N-Body Solver for the Poisson(-Boltzmann) Equation
A Fast N-Body Solver for the Poisson(-Boltzmann) Equation Robert D. Skeel Departments of Computer Science (and Mathematics) Purdue University http://bionum.cs.purdue.edu/2008december.pdf 1 Thesis Poisson(-Boltzmann)
More informationFourier Series and Recent Developments in Analysis
Fourier Series and Recent Developments in Analysis Karlstad, June 2003 Javier Soria (U. Barcelona) 1 Jean Baptiste Joseph Fourier (1768-1830) It was around 1804 that Fourier did his important mathematical
More informationarxiv: v1 [math.ca] 29 Dec 2018
A QUANTITATIVE WEIGHTED WEAK-TYPE ESTIMATE FOR CALDERÓN-ZYGMUND OPERATORS CODY B. STOCKDALE arxiv:82.392v [math.ca] 29 Dec 208 Abstract. The purpose of this article is to provide an alternative proof of
More information3.5 Quadratic Approximation and Convexity/Concavity
3.5 Quadratic Approximation and Convexity/Concavity 55 3.5 Quadratic Approximation and Convexity/Concavity Overview: Second derivatives are useful for understanding how the linear approximation varies
More informationSpectral Algorithms I. Slides based on Spectral Mesh Processing Siggraph 2010 course
Spectral Algorithms I Slides based on Spectral Mesh Processing Siggraph 2010 course Why Spectral? A different way to look at functions on a domain Why Spectral? Better representations lead to simpler solutions
More informationDiffusion Geometries, Global and Multiscale
Diffusion Geometries, Global and Multiscale R.R. Coifman, S. Lafon, MM, J.C. Bremer Jr., A.D. Szlam, P.W. Jones, R.Schul Papers, talks, other materials available at: www.math.yale.edu/~mmm82 Data and functions
More informationA SHORT PROOF OF THE COIFMAN-MEYER MULTILINEAR THEOREM
A SHORT PROOF OF THE COIFMAN-MEYER MULTILINEAR THEOREM CAMIL MUSCALU, JILL PIPHER, TERENCE TAO, AND CHRISTOPH THIELE Abstract. We give a short proof of the well known Coifman-Meyer theorem on multilinear
More informationTime-Frequency Analysis
Time-Frequency Analysis Basics of Fourier Series Philippe B. aval KSU Fall 015 Philippe B. aval (KSU) Fourier Series Fall 015 1 / 0 Introduction We first review how to derive the Fourier series of a function.
More informationUniversal Scalable Robust Solvers from Computational Information Games and fast eigenspace adapted Multiresolution Analysis.
Universal Scalable Robust Solvers from Computational Information Games and fast eigenspace adapted Multiresolution Analysis Houman Owhadi Joint work with Clint Scovel IPAM Apr 3, 2017 DARPA EQUiPS / AFOSR
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 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 informationELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS EXERCISES I (HARMONIC FUNCTIONS)
ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS EXERCISES I (HARMONIC FUNCTIONS) MATANIA BEN-ARTZI. BOOKS [CH] R. Courant and D. Hilbert, Methods of Mathematical Physics, Vol. Interscience Publ. 962. II, [E] L.
More informationRapid evaluation of electrostatic interactions in multi-phase media
Rapid evaluation of electrostatic interactions in multi-phase media P.G. Martinsson, The University of Colorado at Boulder Acknowledgements: Some of the work presented is joint work with Mark Tygert and
More information