Two Simple Nonlinear Edge Detectors
|
|
- Rosalyn Marsh
- 6 years ago
- Views:
Transcription
1 Two Simple Nonlinear Edge Detectors Mladen Victor Wickerhauser Washington University in St. Louis, Missouri International Conference on Wavelet Analysis and Applications 28 May 2004, Chongqing, China Joint work with Wojciech Czaja Research supported in part by NSF grant DMS
2 Introduction Goal: Detect and characterize singularities in images by local dimensionality. Idea: Model images and functions as random vectors. Use an autocovariance matrix as a substitute for the Jacobian, or derivative. Application: Detect points (local dimension 0) and edges (local dimension 1) in images using the eigenvalues of the autocovariance matrix. 1
3 Prior Similar Ideas Old Observation: if a function is rough along some direction, then its local Fourier transform stays big along that direction. Duda and Hart [1972]: Hough transform for detection of lines and curves in images. Mallat and Zhong [1992]: edge detection from wavelet maxima. Aron and Kurz [1997]: detect lines and edges. linear hypothesis testing of variances in small windows to New View: recognize nonsmooth directions from eigenvalues of an autocovariance matrix accumulated from the localized Fourier transform. Use the eigenvector of the larger eigenvalue as the normal to the edge. 2
4 Probabilistic Motivation (2-D) Function φ : T C on the unit circle T R 2. Random direction vector Φ : T C 2 is Φ(θ) def = φ(θ)(cos θ, sin θ), θ T. Direction probability density is φ 2, normalized. Symmetric 2 2 Autocovariance Matrix 2π E(Φ Φ) ij = T Φ i Φ j = φ(θ) 2 τ i (θ)τ j (θ) dθ, 0 for i, j {1, 2}, τ 1 (θ) = sin θ, and τ 2 (θ) = cos θ. Then E(Φ Φ)v, v = T Φ, v 2 attain sup v =1 E(Φ Φ)v, v at a unit eigenvector v R 2 of the largest eigenvalue of E(Φ Φ). Such an eigenvector v always exists. Use v to approximate the max of φ = Φ. 3
5 Example (2-D) For f : R 2 C and large fixed R > 0, define φ(θ) = ˆf(R cos θ, R sin θ) 2 This samples ˆf 2 far out in the direction θ. Suppose φ is highly concentrated near the point θ 0 T. Then E(Φ Φ) is approximately proportional to ( cos 2 ) θ 0 cos θ 0 sin θ 0 cos θ 0 sin θ 0 sin 2. θ 0 Any nonzero vector in the θ 0 direction is an eigenvector of the largest eigenvalue. 4
6 Alternatives (A-1) Replace circle T with disc B: for i, j {1, 2}, and φ L 2 (B). E(Φ Φ) ij = B φ(ξ) 2 ξ i ξ j dξ, (A-2) Replace disk B with R 2, under additional integrability assumptions. 5
7 Localization Start with f : R 2 R, point of interest x, scale ɛ = 1/R. Localize by f gf, for some nonzero radial Schwartz function g : R 2 R: ( ) y x g ɛ (y) = g, ɛ > 0. ɛ which is: smooth, to avoid introducing new singularities; radial, to avoid introducing directional bias; nonzero and concentrated, to emphasize the point of interest x. 6
8 Dual Local Autocovariance Definition 1 The dual local autocovariance matrix of f at x is the 2 2 matrix whose ij-coefficient is: E ɛ,g (f; x) ij = B(0,1/ɛ) ξ iξ j ĝɛf(ξ) 2 dξ. This is the real, symmetric second moment matrix of the unnormalized probability density function ĝɛf 2. Localization gf is integrable are well defined. ĝf is bounded and continuous matrix coefficients f 0 near x E ɛ,g (f; x) is positive definite. 7
9 Straight Edges Change of variables: E ɛ,g (f; x) ij = B ξ iξ j ( g R 2 y + ɛ 1 ɛ ) x f(ɛy)e 2πiy ξ dy 2 dξ. If f is homogeneous of degree 0, and x = 0, this formula is independent of ɛ > 0 Example: f = 1 L, indicator function of left half-plane L = {(x 1, x 2 ) : x 1 0}. E ɛ,g (1 L ; 0) is a matrix with no ɛ-dependence. Use g(x) = exp( π x 2 ), for ease of computing Fourier transforms. Get dual local autocovariance matrix explicitly for the edge x 1 = 0. 8
10 Straight edges (continued... ) Lemma 1 E ɛ,g (1 L ; 0) has the eigensystem ( ) 1 λ 1, λ 2 0 satisfying λ 1 > λ 2 > 0, independently of ɛ. ( ) 0, 1 Numerical estimate: the matrix is diagonal, so λ i = B ξ2 i 0 2 e πξ2 2 e πx2 2πixξ 1 dx for i = 1, 2. In terms of the Gaussian integral G(ξ) = ξ 0 exp(πx2 ) dx, this becomes: [ ] 1 λ i = B ξ2 i e 2π ξ G2 (ξ 1 ) dξ, for i = 1, 2. dξ, By Mathematica, λ 1 λ
11 Off the Straight Edge Conversely, take g(ξ) to be a smooth radial function supported in B(0, 1), so: for any x / L, and 0 < ɛ < dist (L, x) small enough, since g ɛ 1 L = 0. E ɛ,g (1 L ; x) 0, for each point x IntL, E ɛ,g (1 L ; x) = const I, for all sufficiently small 0 < ɛ < dist (x, L). In either case, E has two equal eigenvalues. 10
12 Rotation and Translation Translation moves x, preserves eigenvectors: Lemma 2 Let T be a translation by x in R 2 : f T (y) = f(y + x). Then E ɛ,g (f; x) = E ɛ,g (f T ; 0). Rotation about x rotates eigenvectors, changes edge direction: Lemma 3 Let U be a rotation about x in R 2. Then E ɛ,g (f U; x) = U E ɛ,g (f; x) U 1. 11
13 Straight Edge Detection Let 1 H be the characteristic function of a half-plane H = {x R 2 : ν x β}, with nonzero normal vector ν R 2, for some constant β R. Its edge is the line H. Theorem 4 x H iff E ɛ,g (1 H ; x) has distinct eigenvalues λ 1 > λ 2 > 0 for every smooth, radial function g 0, and ɛ > 0. In that case, ν will be an eigenvector of the larger eigenvalue. Note; degree-0 homogeneity of 1 H eigenvalues of E ɛ,g (1 H ; x) do not depend on ɛ. For x / H off the edge, and each smooth, compactly-supported radial function g, E ɛ,g (1 H ; x) = const I has eventually constant and equal eigenvalues as ɛ 0. 12
14 Smooth-Curve Edges Idea: Boundaries of smooth domains look like boundaries of half-planes. Fix: Domain D R 2, smooth boundary D. 1 D, the characteristic function of D. g : R 2 R, a smooth radial function, supported in B = B(0, 1). H, half-plane with H D at some x. As ɛ 0, expect: x D distinct eigenvalues for E ɛ,g (1 D, x). x / D equal eigenvalues for E ɛ,g (1 D, x). 13
15 Smooth-Curve Edges (continued... ) Lemma 5 In any matrix norm, as ɛ 0. E ɛ,g (1 H ; x) E ɛ,g (1 D ; x) = O(ɛ), (Note: for Schwartz g, any δ > 0 gives as ɛ 0.) E ɛ,g (1 H ; x) E ɛ,g (1 D ; x) = O(ɛ 1 δ ), Theorem 6 For x D, and eigenvalues λ 1 (ɛ) and λ 2 (ɛ) of E ɛ,g (1 D ; x), we have For x / D, λ 1 (ɛ) λ 2 (ɛ) 0 as ɛ 0. lim inf ɛ 0+ λ 1(ɛ) λ 2 (ɛ) > 0. 14
16 Affine and Differentiable Functions At points x where f is differentiable, the eigenvalues λ 1 (ɛ), λ 2 (ɛ) of E ɛ,g (f; x) are equal in the limit: Lemma 7 Let A : R 2 R be an affine function: A(x) = a x + b for some constants a R 2, b R. Let g : R 2 R be a smooth radial function supported in B = B(0, 1). Then: lim [λ 1(ɛ) λ 2 (ɛ)] = 0. ɛ 0+ Lemma 8 If f : R 2 R is differentiable at x with affine A tangent to f at x, then as ɛ 0. E ɛ,g (f; x) E ɛ,g (A; x) = o(ɛ), (Continuously differentiable f near x satisfies as ɛ 0.) E ɛ,g (f; x) E ɛ,g (A; x) = O(ɛ 2 ), 15
17 Detecting Non-Edge Points Theorem 9 Suppose that f : R 2 R is differentiable at x. radial compactly-supported function g : R 2 R, Then for any smooth E ɛ,g (f; x) const I, as ɛ 0. Converse is false: lim ɛ 0+ E ɛ,g (f; x) = const I does not imply that f is differentiable at x, or even continuous. Symmetry can masquerade as smoothness, as in f = 1 x1,x 2 >0. Fix g(x) = exp( π x 2 ) as before. Then E ɛ,g (f; 0) = const I 0 for every ɛ > 0, so λ 1 (ɛ) = λ 2 (ɛ) = λ > 0 for every ɛ, even though f is discontinuous at 0. 16
18 Higher-Dimensional Theory Study the geometry of f : R p R d in high dimensions p, d. Suppose f = (f 1,..., f d ) is polynomially bounded. Fix x R p, Schwartz g : R p R d with radial components, and g ɛ (y) def = g ( ) y x, as before, for ɛ > 0. Definition 2 The dual local autocovariance matrix of f at x is the p p matrix: E ɛ,g (f; x) ij = for i, j {1,..., p}. ɛ ξ iξ j f g ɛ (ξ) 2 dξ, B(0,1/ɛ) Smooth p p case is like smooth 2 2: Theorem 10 Suppose that f : R p R d is differentiable at x. Then for any smooth radial function g : R p R d of compact support, the matrix E ɛ,g (f; x) const I, as ɛ 0. 17
19 Discretization Discrete Fourier transform on N real samples {f(n) : 0 n < N}: for integer k B N = [ N 2, N 2 ]. ˆf(k) = N 1 n=0 exp ( 2πi kn N ) f(n), If only the first q N samples of f are nonzero, then the sum is over {0, 1,..., q 1}. When f is real-valued, for k B N. ˆf(k) 2 = q 1 n,n =0 exp ( 2πi k(n n ) N ) f(n)f(n ), 18
20 Discretization (continued... ) The r-th moment of ˆf(k) 2 is q 1 n,n =0 k B N k r ˆf(k) 2 = f(n)f(n ) k r exp k B N ( 2πi k(n n ) N ). The innermost sum in k, if normalized with N r+1, is a Riemann sum for the integral evaluated at n (n n ): µ r (n) = 1 2 µ 0 (n) = 1 2 x r exp( 2πinx) dx, { 1, if n = 0, 0, otherwise; µ 1 (n) = µ 2 (n) = { 0, if n = 0, i ( 1)n 2πn, otherwise; { 1/12, if n = 0, ( 1) n 2π 2 n 2, otherwise. 19
21 Discrete Dual Local Autocovariance Image: f on {(m, n) : 0 m < M; 0 n < N}. Bump: g on {(m, n) : 0 m < p; 0 n < q}, for p M; q N. Autocovariance matrix: localize f gf, then compute E 11 = p 1 q 1 m,m =0 n=0 f(m, n)f(m, n)µ 2 (m m ); E 22 = p 1 q 1 m=0 n,n =0 f(m, n)f(m, n )µ 2 (n n ); E 12 = E 21 = p 1 q 1 m,m =0 n,n =0 f(m, n)f(m, n )µ 1 (m m )µ 1 (n n ). 20
22 Implementation in General Around each point z of the image, do: Localization. Extract pixel values on the square subgrid contained in z + [ ɛ, ɛ] 2. Multiply f(z + y) by the weighting bump function g(y). [O(ɛ 2 ) operations per pixel.] Dual autocovariance. Compute E = (E ij ), a real-valued, symmetric, positive semidefinite 2 2 matrix. [O(ɛ 4 ) operations per pixel.] Eigenvalues. For symmetric 2 2 matrices E: λ = 1 ( E 11 + E 22 ± (E 11 E 22 ) 2 + 4E where + gives λ 1, gives λ 2. [O(1) operations per pixel.] ), Edginess. Compute ratio λ 1 /λ 2 or difference λ 1 λ 2, or use reciprocal λ 2 /λ 1 [0, 1] in write-black mode so darker edgier. [O(1) operations per pixel.] 21
23 Special 7x7 Gridpoint Case Let z be of the form (m, n), a gridpoint with integer coordinates. Use ɛ = 3, to localize to 7 7 subgrids. Zero-pad at boundaries. Let g be the Gaussian bump g(y) = exp( π y 2 /ɛ 2 ). Use reciprocal edginess e = λ 2 /λ 1 [0, 1] in write-black mode so darker edgier. 22
24 Special 2x2 Midpoint Case Let z be of the form (m + 1 2, n ), a midpoint. Let ɛ = 1, so E is computed from the pixels immediately NE, SE, SW, and NW of z: (m, n + 1) (m + 1, n + 1) z (m, n) (m + 1, n) def = NW f(m, n + 1) SW def = f(m, n) [z] NE def = f(m + 1, n + 1) SE def = f(m + 1, n) Use eigenvalue difference for edginess: e def = λ 1 (ɛ, z) λ 2 (ɛ, z), All pixels are equidistant from z, so normalize g = 1 at that distance. 23
25 Special 2x2 Midpoint Case (continued... ) Evaluate µ 1 (0) = 0, µ 1 (±1) = 2πi ±1, µ 2(0) = 12 1,and µ 2(±1) = 1 2, so E 11 = (NW 2 + NE 2 )µ 2 (0 0) + (SW 2 + +SE 2 )µ 2 (1 1) + +(NW SW + NE SE) (µ 2 (0 1) + µ 2 (1 0)) = 1 ( NW 2 + NE 2 + SW 2 + SE 2 ) 1 (NW SW + NE SE) ; 12 π2 E 22 = 1 ( NW 2 + NE 2 + SW 2 + SE 2 ) 1 (NW NE + SW SE) ; 12 π2 E 11 E 22 = 1 (NW NE + SW SE NW SW NE SE) π2 = 1 (NW SE) (NE SW ) ; π2 E 12 = E 21 = 1 (NW SE NE SW ). 2π2 Edginess formula is a nonlinear four-tap filter with a square 2 2 pixel mask: e = λ 1 λ 2 = 1 π 2 (NW SE) 2 (NE SW ) 2 + (NW SE NE SW ) 2. 2π 24
26 Special 4+1 Centerpoint Case Let z be of the form (m, n). Let ɛ = 1, so E is computed from the central pixel at z = (m, n) and the 4 pixels N, E, S, and W of z: (m, n+1) N def = f(m, n+1) (m 1, n) z = (m, n) (m+1, n) def = W f(m 1, n) C def = f(m, n) E def = f(m+1, n) (m, n 1) S def = f(m, n 1) Use eigenvalue difference for edginess: e def = λ 1 (ɛ, z) λ 2 (ɛ, z), Use g with g(0) = 1 and g(1) = t > 0, so z gets 1 while N,E,S,W get t. 25
27 Special 4+1 Centerpoint Case (continued... ) Evaluate µ 1 (0) = 0, µ 1 (±1) = 2πi ±1, µ 1(±2) = 4π ±i, µ 2(0) = 12 1, µ 2(±1) = 1 2π2, and µ 2 (±2) = 1 8π 2: D def = 1 12 (C2 + t 2 [N 2 + E 2 + S 2 + W 2 ]); E 11 = D t t2 π2(n C + C S) + (N S) ; 4π2 E 22 = D t t2 π2(e C + C W ) + 4π 2 (E W ) ; E 12 = E 21 = t 2 2π 2 (N W N E + S E S W ) = t2 (N S) (E W ). 2π2 Edginess formule is a nonlinear five-tap filter with a mask shaped like a plus sign: λ 1 λ 2 = t [ π 2 C (E+W N S) + t 4 ] 2 (N S E W ) +t 2 (N S) 2 (E W ) 2. Each weight t (0, 1) tunes a different filter. Use t = 0.75 for no good reason. The choice t = 0 gives e = 0 everywhere. 26
28 Example Images 1. Geometrical figures, piecewise constant functions with jump discontinuities along rectifiable, mostly smooth curves. 2. Toy house demonstration image from the XHoughTool website. 3. Fingerprint, part of the FBI/NIST compliance test suite for WSQ compression. 4. Lena Sjøblum, the famous woman in a hat. 5. Cone, a ray-traced image by Craig Kolb. 6. Truck, one frame of video from Peng Li. 27
29 Example: Geometrical Figures 28
30 Geometrical Figures 7x7 Gridpoint Edginess 29
31 Geometrical Figures 2x2 Midpoint Edginess 30
32 Geometrical Figures 4+1 Centerpoint Edginess 31
33 Example: Toy House Image 32
34 Toy House 7x7 Gridpoint Edginess 33
35 Toy House 2x2 Midpoint Edginess 34
36 Toy House 4+1 Centerpoint Edginess 35
37 Example: Fingerprint Image 36
38 Fingerprint 7x7 Gridpoint Edginess 37
39 Fingerprint 2x2 Midpoint Edginess 38
40 Fingerprint 4+1 Centerpoint Edginess 39
41 Example: Lena Sjøblum Image 40
42 Lena 7x7 Gridpoint Edginess 41
43 Lena 2x2 Midpoint Edginess 42
44 Lena 4+1 Centerpoint Edginess 43
45 Example: Cone Image 44
46 Cone 7x7 Gridpoint Edginess 45
47 Cone 2x2 Midpoint Edginess 46
48 Cone 4+1 Centerpoint Edginess 47
49 Example: Truck Video Frame 48
50 Truck 7x7 Gridpoint Edginess 49
51 Truck 2x2 Midpoint Edginess 50
52 Truck 4+1 Centerpoint Edginess 51
53 Comparison: Toy House Image 7x7, 2x2, 4+1 Edginess 52
54 Comparison: Fingerprint Image 7x7, 2x2, 4+1 Edginess 53
55 Comparison: Lena Sjøblum Image 7x7, 2x2, 4+1 Edginess 54
56 Comparison: Cone Image 7x7, 2x2, 4+1 Edginess 55
57 Comparison: Truck Video Frame 7x7, 2x2, 4+1 Edginess 56
A Simple Nonlinear Filter for Edge Detection in Images
A Simple Nonlinear Filter for Edge Detection in Images laden Victor Wickerhauser Department of athematics, Washington University, St. Louis, 0 63130 USA and Wojciech Czaja Instytut atematyczny, Uniwersytet
More information17 The functional equation
18.785 Number theory I Fall 16 Lecture #17 11/8/16 17 The functional equation In the previous lecture we proved that the iemann zeta function ζ(s) has an Euler product and an analytic continuation to the
More informationFOURIER TRANSFORMS. 1. Fourier series 1.1. The trigonometric system. The sequence of functions
FOURIER TRANSFORMS. Fourier series.. The trigonometric system. The sequence of functions, cos x, sin x,..., cos nx, sin nx,... is called the trigonometric system. These functions have period π. The trigonometric
More informationLinear Independence of Finite Gabor Systems
Linear Independence of Finite Gabor Systems 1 Linear Independence of Finite Gabor Systems School of Mathematics Korea Institute for Advanced Study Linear Independence of Finite Gabor Systems 2 Short trip
More informationTopics in Harmonic Analysis Lecture 1: The Fourier transform
Topics in Harmonic Analysis Lecture 1: The Fourier transform Po-Lam Yung The Chinese University of Hong Kong Outline Fourier series on T: L 2 theory Convolutions The Dirichlet and Fejer kernels Pointwise
More informationMATH 220 solution to homework 5
MATH 220 solution to homework 5 Problem. (i Define E(t = k(t + p(t = then E (t = 2 = 2 = 2 u t u tt + u x u xt dx u 2 t + u 2 x dx, u t u xx + u x u xt dx x [u tu x ] dx. Because f and g are compactly
More information2 Infinite products and existence of compactly supported φ
415 Wavelets 1 Infinite products and existence of compactly supported φ Infinite products.1 Infinite products a n need to be defined via limits. But we do not simply say that a n = lim a n N whenever the
More information1.5 Approximate Identities
38 1 The Fourier Transform on L 1 (R) which are dense subspaces of L p (R). On these domains, P : D P L p (R) and M : D M L p (R). Show, however, that P and M are unbounded even when restricted to these
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationInternational Competition in Mathematics for Universtiy Students in Plovdiv, Bulgaria 1994
International Competition in Mathematics for Universtiy Students in Plovdiv, Bulgaria 1994 1 PROBLEMS AND SOLUTIONS First day July 29, 1994 Problem 1. 13 points a Let A be a n n, n 2, symmetric, invertible
More informationTOOLS FROM HARMONIC ANALYSIS
TOOLS FROM HARMONIC ANALYSIS BRADLY STADIE Abstract. The Fourier transform can be thought of as a map that decomposes a function into oscillatory functions. In this paper, we will apply this decomposition
More informationLecture Note 5: Semidefinite Programming for Stability Analysis
ECE7850: Hybrid Systems:Theory and Applications Lecture Note 5: Semidefinite Programming for Stability Analysis Wei Zhang Assistant Professor Department of Electrical and Computer Engineering Ohio State
More informationMATH 220 solution to homework 4
MATH 22 solution to homework 4 Problem. Define v(t, x) = u(t, x + bt), then v t (t, x) a(x + u bt) 2 (t, x) =, t >, x R, x2 v(, x) = f(x). It suffices to show that v(t, x) F = max y R f(y). We consider
More informationa k 0, then k + 1 = 2 lim 1 + 1
Math 7 - Midterm - Form A - Page From the desk of C. Davis Buenger. https://people.math.osu.edu/buenger.8/ Problem a) [3 pts] If lim a k = then a k converges. False: The divergence test states that if
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationProblem 1A. Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1. (Hint: 1. Solution: The volume is 1. Problem 2A.
Problem 1A Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1 (Hint: 1 1 (something)dz) Solution: The volume is 1 1 4xydz where x = y = 1 z 2 This integral has value 16/3 Problem 2A Let f(x)
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationFRAMES AND TIME-FREQUENCY ANALYSIS
FRAMES AND TIME-FREQUENCY ANALYSIS LECTURE 5: MODULATION SPACES AND APPLICATIONS Christopher Heil Georgia Tech heil@math.gatech.edu http://www.math.gatech.edu/ heil READING For background on Banach spaces,
More informationComplex Analysis Important Concepts
Complex Analysis Important Concepts Travis Askham April 1, 2012 Contents 1 Complex Differentiation 2 1.1 Definition and Characterization.............................. 2 1.2 Examples..........................................
More informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
More information1 Fourier Integrals on L 2 (R) and L 1 (R).
18.103 Fall 2013 1 Fourier Integrals on L 2 () and L 1 (). The first part of these notes cover 3.5 of AG, without proofs. When we get to things not covered in the book, we will start giving proofs. The
More informationReal Analysis Problems
Real Analysis Problems Cristian E. Gutiérrez September 14, 29 1 1 CONTINUITY 1 Continuity Problem 1.1 Let r n be the sequence of rational numbers and Prove that f(x) = 1. f is continuous on the irrationals.
More informationMath Homework 2
Math 73 Homework Due: September 8, 6 Suppose that f is holomorphic in a region Ω, ie an open connected set Prove that in any of the following cases (a) R(f) is constant; (b) I(f) is constant; (c) f is
More informationSolution to Homework 4, Math 7651, Tanveer 1. Show that the function w(z) = 1 2
Solution to Homework 4, Math 7651, Tanveer 1. Show that the function w(z) = 1 (z + 1/z) maps the exterior of a unit circle centered around the origin in the z-plane to the exterior of a straight line cut
More informationLet R be the line parameterized by x. Let f be a complex function on R that is integrable. The Fourier transform ˆf = F f is. e ikx f(x) dx. (1.
Chapter 1 Fourier transforms 1.1 Introduction Let R be the line parameterized by x. Let f be a complex function on R that is integrable. The Fourier transform ˆf = F f is ˆf(k) = e ikx f(x) dx. (1.1) It
More informationSTA 294: Stochastic Processes & Bayesian Nonparametrics
MARKOV CHAINS AND CONVERGENCE CONCEPTS Markov chains are among the simplest stochastic processes, just one step beyond iid sequences of random variables. Traditionally they ve been used in modelling a
More informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
More informationCones of measures. Tatiana Toro. University of Washington. Quantitative and Computational Aspects of Metric Geometry
Cones of measures Tatiana Toro University of Washington Quantitative and Computational Aspects of Metric Geometry Based on joint work with C. Kenig and D. Preiss Tatiana Toro (University of Washington)
More informationCHAPTER VI APPLICATIONS TO ANALYSIS
CHAPTER VI APPLICATIONS TO ANALYSIS We include in this chapter several subjects from classical analysis to which the notions of functional analysis can be applied. Some of these subjects are essential
More informationFourier Series. ,..., e ixn ). Conversely, each 2π-periodic function φ : R n C induces a unique φ : T n C for which φ(e ix 1
Fourier Series Let {e j : 1 j n} be the standard basis in R n. We say f : R n C is π-periodic in each variable if f(x + πe j ) = f(x) x R n, 1 j n. We can identify π-periodic functions with functions on
More informationFFTs in Graphics and Vision. Homogenous Polynomials and Irreducible Representations
FFTs in Graphics and Vision Homogenous Polynomials and Irreducible Representations 1 Outline The 2π Term in Assignment 1 Homogenous Polynomials Representations of Functions on the Unit-Circle Sub-Representations
More informationSome basic elements of Probability Theory
Chapter I Some basic elements of Probability Theory 1 Terminology (and elementary observations Probability theory and the material covered in a basic Real Variables course have much in common. However
More informationSize properties of wavelet packets generated using finite filters
Rev. Mat. Iberoamericana, 18 (2002, 249 265 Size properties of wavelet packets generated using finite filters Morten Nielsen Abstract We show that asymptotic estimates for the growth in L p (R- norm of
More informationMultidimensional partitions of unity and Gaussian terrains
and Gaussian terrains Richard A. Bale, Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureu ABSTRACT Partitions of unity play an important rôle as amplitude-preserving windows for nonstationary
More informationIn particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with
Appendix: Matrix Estimates and the Perron-Frobenius Theorem. This Appendix will first present some well known estimates. For any m n matrix A = [a ij ] over the real or complex numbers, it will be convenient
More informationComplex Analysis, Stein and Shakarchi The Fourier Transform
Complex Analysis, Stein and Shakarchi Chapter 4 The Fourier Transform Yung-Hsiang Huang 2017.11.05 1 Exercises 1. Suppose f L 1 (), and f 0. Show that f 0. emark 1. This proof is observed by Newmann (published
More informationHW1 solutions. 1. α Ef(x) β, where Ef(x) is the expected value of f(x), i.e., Ef(x) = n. i=1 p if(a i ). (The function f : R R is given.
HW1 solutions Exercise 1 (Some sets of probability distributions.) Let x be a real-valued random variable with Prob(x = a i ) = p i, i = 1,..., n, where a 1 < a 2 < < a n. Of course p R n lies in the standard
More informationMath 127C, Spring 2006 Final Exam Solutions. x 2 ), g(y 1, y 2 ) = ( y 1 y 2, y1 2 + y2) 2. (g f) (0) = g (f(0))f (0).
Math 27C, Spring 26 Final Exam Solutions. Define f : R 2 R 2 and g : R 2 R 2 by f(x, x 2 (sin x 2 x, e x x 2, g(y, y 2 ( y y 2, y 2 + y2 2. Use the chain rule to compute the matrix of (g f (,. By the chain
More informationDecay profiles of a linear artificial viscosity system
Decay profiles of a linear artificial viscosity system Gene Wayne, Ryan Goh and Roland Welter Boston University rwelter@bu.edu July 2, 2018 This research was supported by funding from the NSF. Roland Welter
More informationTEST CODE: MMA (Objective type) 2015 SYLLABUS
TEST CODE: MMA (Objective type) 2015 SYLLABUS Analytical Reasoning Algebra Arithmetic, geometric and harmonic progression. Continued fractions. Elementary combinatorics: Permutations and combinations,
More informationWavelets in Image Compression
Wavelets in Image Compression M. Victor WICKERHAUSER Washington University in St. Louis, Missouri victor@math.wustl.edu http://www.math.wustl.edu/~victor THEORY AND APPLICATIONS OF WAVELETS A Workshop
More informationBernstein s inequality and Nikolsky s inequality for R d
Bernstein s inequality and Nikolsky s inequality for d Jordan Bell jordan.bell@gmail.com Department of athematics University of Toronto February 6 25 Complex Borel measures and the Fourier transform Let
More informationComplex Analysis, Stein and Shakarchi Meromorphic Functions and the Logarithm
Complex Analysis, Stein and Shakarchi Chapter 3 Meromorphic Functions and the Logarithm Yung-Hsiang Huang 217.11.5 Exercises 1. From the identity sin πz = eiπz e iπz 2i, it s easy to show its zeros are
More information7: FOURIER SERIES STEVEN HEILMAN
7: FOURIER SERIES STEVE HEILMA Contents 1. Review 1 2. Introduction 1 3. Periodic Functions 2 4. Inner Products on Periodic Functions 3 5. Trigonometric Polynomials 5 6. Periodic Convolutions 7 7. Fourier
More informationSolutions to Complex Analysis Prelims Ben Strasser
Solutions to Complex Analysis Prelims Ben Strasser In preparation for the complex analysis prelim, I typed up solutions to some old exams. This document includes complete solutions to both exams in 23,
More informationx. Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ 2 ).
.8.6 µ =, σ = 1 µ = 1, σ = 1 / µ =, σ =.. 3 1 1 3 x Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ ). The Gaussian distribution Probably the most-important distribution in all of statistics
More information1 Distributions (due January 22, 2009)
Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =
More informationMapping problems and harmonic univalent mappings
Mapping problems and harmonic univalent mappings Antti Rasila Helsinki University of Technology antti.rasila@tkk.fi (Mainly based on P. Duren s book Harmonic mappings in the plane) Helsinki Analysis Seminar,
More informationQUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1)
QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1) Each of the six questions is worth 10 points. 1) Let H be a (real or complex) Hilbert space. We say
More informationProblem 1A. Suppose that f is a continuous real function on [0, 1]. Prove that
Problem 1A. Suppose that f is a continuous real function on [, 1]. Prove that lim α α + x α 1 f(x)dx = f(). Solution: This is obvious for f a constant, so by subtracting f() from both sides we can assume
More informationSOLUTION OF THE DIRICHLET PROBLEM WITH A VARIATIONAL METHOD. 1. Dirichlet integral
SOLUTION OF THE DIRICHLET PROBLEM WITH A VARIATIONAL METHOD CRISTIAN E. GUTIÉRREZ FEBRUARY 3, 29. Dirichlet integral Let f C( ) with open and bounded. Let H = {u C ( ) : u = f on } and D(u) = Du(x) 2 dx.
More informationLecture 7 Local properties of analytic functions Part 1 MATH-GA Complex Variables
Lecture 7 Local properties of analytic functions Part 1 MATH-GA 2451.001 omplex Variables 1 Removable singularities 1.1 Riemann s removable singularity theorem We have said that auchy s integral formula
More informationTheorem [Mean Value Theorem for Harmonic Functions] Let u be harmonic on D(z 0, R). Then for any r (0, R), u(z 0 ) = 1 z z 0 r
2. A harmonic conjugate always exists locally: if u is a harmonic function in an open set U, then for any disk D(z 0, r) U, there is f, which is analytic in D(z 0, r) and satisfies that Re f u. Since such
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More information1 Fourier Integrals of finite measures.
18.103 Fall 2013 1 Fourier Integrals of finite measures. Denote the space of finite, positive, measures on by M + () = {µ : µ is a positive measure on ; µ() < } Proposition 1 For µ M + (), we define the
More informationTopics in Harmonic Analysis, Sparse Representations, and Data Analysis
Topics in Harmonic Analysis, Sparse Representations, and Data Analysis Norbert Wiener Center Department of Mathematics University of Maryland, College Park http://www.norbertwiener.umd.edu Thesis Defense,
More informationLinear Algebra And Its Applications Chapter 6. Positive Definite Matrix
Linear Algebra And Its Applications Chapter 6. Positive Definite Matrix KAIST wit Lab 2012. 07. 10 남성호 Introduction The signs of the eigenvalues can be important. The signs can also be related to the minima,
More information1. Let A be a 2 2 nonzero real matrix. Which of the following is true?
1. Let A be a 2 2 nonzero real matrix. Which of the following is true? (A) A has a nonzero eigenvalue. (B) A 2 has at least one positive entry. (C) trace (A 2 ) is positive. (D) All entries of A 2 cannot
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 information1 Continuity Classes C m (Ω)
0.1 Norms 0.1 Norms A norm on a linear space X is a function : X R with the properties: Positive Definite: x 0 x X (nonnegative) x = 0 x = 0 (strictly positive) λx = λ x x X, λ C(homogeneous) x + y x +
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 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 informationHölder regularity estimation by Hart Smith and Curvelet transforms
Hölder regularity estimation by Hart Smith and Curvelet transforms Jouni Sampo Lappeenranta University Of Technology Department of Mathematics and Physics Finland 18th September 2007 This research is done
More information8 Nonlinear change of variables
Tel Aviv University, 206 Analysis-III 98 8 Nonlinear change of variables 8a Introduction....................... 98 8b Examples......................... 00 8c Measure 0 is preserved................. 02
More informationTangent spaces, normals and extrema
Chapter 3 Tangent spaces, normals and extrema If S is a surface in 3-space, with a point a S where S looks smooth, i.e., without any fold or cusp or self-crossing, we can intuitively define the tangent
More informationSynopsis of Complex Analysis. Ryan D. Reece
Synopsis of Complex Analysis Ryan D. Reece December 7, 2006 Chapter Complex Numbers. The Parts of a Complex Number A complex number, z, is an ordered pair of real numbers similar to the points in the real
More information4 Riesz Kernels. Since the functions k i (ξ) = ξ i. are bounded functions it is clear that R
4 Riesz Kernels. A natural generalization of the Hilbert transform to higher dimension is mutiplication of the Fourier Transform by homogeneous functions of degree 0, the simplest ones being R i f(ξ) =
More informationINTRODUCTION TO REAL ANALYTIC GEOMETRY
INTRODUCTION TO REAL ANALYTIC GEOMETRY KRZYSZTOF KURDYKA 1. Analytic functions in several variables 1.1. Summable families. Let (E, ) be a normed space over the field R or C, dim E
More informationMore on Fourier Series
More on Fourier Series R. C. Trinity University Partial Differential Equations Lecture 6.1 New Fourier series from old Recall: Given a function f (x, we can dilate/translate its graph via multiplication/addition,
More informationX n D X lim n F n (x) = F (x) for all x C F. lim n F n(u) = F (u) for all u C F. (2)
14:17 11/16/2 TOPIC. Convergence in distribution and related notions. This section studies the notion of the so-called convergence in distribution of real random variables. This is the kind of convergence
More informationbe the set of complex valued 2π-periodic functions f on R such that
. Fourier series. Definition.. Given a real number P, we say a complex valued function f on R is P -periodic if f(x + P ) f(x) for all x R. We let be the set of complex valued -periodic functions f on
More informationA quick introduction to regularity structures
A quick introduction to regularity structures Lorenzo Zambotti (LPMA, Univ. Paris 6) Workshop "Singular SPDEs" Foreword The theory of Regularity Structures by Martin Hairer (2014) is a major advance in
More informationLecture Notes on the Gaussian Distribution
Lecture Notes on the Gaussian Distribution Hairong Qi The Gaussian distribution is also referred to as the normal distribution or the bell curve distribution for its bell-shaped density curve. There s
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 informationMath Final Exam.
Math 106 - Final Exam. This is a closed book exam. No calculators are allowed. The exam consists of 8 questions worth 100 points. Good luck! Name: Acknowledgment and acceptance of honor code: Signature:
More informationIntroduction to Mathematical Programming
Introduction to Mathematical Programming Ming Zhong Lecture 25 November 5, 2018 Ming Zhong (JHU) AMS Fall 2018 1 / 19 Table of Contents 1 Ming Zhong (JHU) AMS Fall 2018 2 / 19 Some Preliminaries: Fourier
More informationProblem 1A. Calculus. Problem 3A. Real analysis. f(x) = 0 x = 0.
Problem A. Calculus Find the length of the spiral given in polar coordinates by r = e θ, < θ 0. Solution: The length is 0 θ= ds where by Pythagoras ds = dr 2 + (rdθ) 2 = dθ 2e θ, so the length is 0 2e
More informationFeature extraction: Corners and blobs
Feature extraction: Corners and blobs Review: Linear filtering and edge detection Name two different kinds of image noise Name a non-linear smoothing filter What advantages does median filtering have over
More informationDalal-Schmutz (2002) and Diaconis-Freedman (1999) 1. Random Compositions. Evangelos Kranakis School of Computer Science Carleton University
Dalal-Schmutz (2002) and Diaconis-Freedman (1999) 1 Random Compositions Evangelos Kranakis School of Computer Science Carleton University Dalal-Schmutz (2002) and Diaconis-Freedman (1999) 2 Outline 1.
More informationAn idea how to solve some of the problems. diverges the same must hold for the original series. T 1 p T 1 p + 1 p 1 = 1. dt = lim
An idea how to solve some of the problems 5.2-2. (a) Does not converge: By multiplying across we get Hence 2k 2k 2 /2 k 2k2 k 2 /2 k 2 /2 2k 2k 2 /2 k. As the series diverges the same must hold for the
More informationMathematical Tripos Part III Michaelmas 2017 Distribution Theory & Applications, Example sheet 1 (answers) Dr A.C.L. Ashton
Mathematical Tripos Part III Michaelmas 7 Distribution Theory & Applications, Eample sheet (answers) Dr A.C.L. Ashton Comments and corrections to acla@damtp.cam.ac.uk.. Construct a non-zero element of
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 informationKernel methods, kernel SVM and ridge regression
Kernel methods, kernel SVM and ridge regression Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012 Collaborative Filtering 2 Collaborative Filtering R: rating matrix; U: user factor;
More information= F (b) F (a) F (x i ) F (x i+1 ). a x 0 x 1 x n b i
Real Analysis Problem 1. If F : R R is a monotone function, show that F T V ([a,b]) = F (b) F (a) for any interval [a, b], and that F has bounded variation on R if and only if it is bounded. Here F T V
More informationMeasure and Integration: Solutions of CW2
Measure and Integration: s of CW2 Fall 206 [G. Holzegel] December 9, 206 Problem of Sheet 5 a) Left (f n ) and (g n ) be sequences of integrable functions with f n (x) f (x) and g n (x) g (x) for almost
More informationTEST CODE: MIII (Objective type) 2010 SYLLABUS
TEST CODE: MIII (Objective type) 200 SYLLABUS Algebra Permutations and combinations. Binomial theorem. Theory of equations. Inequalities. Complex numbers and De Moivre s theorem. Elementary set theory.
More informationMATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5
MATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5.. The Arzela-Ascoli Theorem.. The Riemann mapping theorem Let X be a metric space, and let F be a family of continuous complex-valued functions on X. We have
More informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More information13. Fourier transforms
(December 16, 2017) 13. Fourier transforms Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/real/notes 2017-18/13 Fourier transforms.pdf]
More informationReview of Multi-Calculus (Study Guide for Spivak s CHAPTER ONE TO THREE)
Review of Multi-Calculus (Study Guide for Spivak s CHPTER ONE TO THREE) This material is for June 9 to 16 (Monday to Monday) Chapter I: Functions on R n Dot product and norm for vectors in R n : Let X
More informationMA3111S COMPLEX ANALYSIS I
MA3111S COMPLEX ANALYSIS I 1. The Algebra of Complex Numbers A complex number is an expression of the form a + ib, where a and b are real numbers. a is called the real part of a + ib and b the imaginary
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationMADHAVA MATHEMATICS COMPETITION, December 2015 Solutions and Scheme of Marking
MADHAVA MATHEMATICS COMPETITION, December 05 Solutions and Scheme of Marking NB: Part I carries 0 marks, Part II carries 30 marks and Part III carries 50 marks Part I NB Each question in Part I carries
More informationMa 416: Complex Variables Solutions to Homework Assignment 6
Ma 46: omplex Variables Solutions to Homework Assignment 6 Prof. Wickerhauser Due Thursday, October th, 2 Read R. P. Boas, nvitation to omplex Analysis, hapter 2, sections 9A.. Evaluate the definite integral
More informationFunctional Analysis Review
Outline 9.520: Statistical Learning Theory and Applications February 8, 2010 Outline 1 2 3 4 Vector Space Outline A vector space is a set V with binary operations +: V V V and : R V V such that for all
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 informationMobile Robotics 1. A Compact Course on Linear Algebra. Giorgio Grisetti
Mobile Robotics 1 A Compact Course on Linear Algebra Giorgio Grisetti SA-1 Vectors Arrays of numbers They represent a point in a n dimensional space 2 Vectors: Scalar Product Scalar-Vector Product Changes
More informationAPPROXIMATING CONTINUOUS FUNCTIONS: WEIERSTRASS, BERNSTEIN, AND RUNGE
APPROXIMATING CONTINUOUS FUNCTIONS: WEIERSTRASS, BERNSTEIN, AND RUNGE WILLIE WAI-YEUNG WONG. Introduction This set of notes is meant to describe some aspects of polynomial approximations to continuous
More informationDetermine for which real numbers s the series n>1 (log n)s /n converges, giving reasons for your answer.
Problem A. Determine for which real numbers s the series n> (log n)s /n converges, giving reasons for your answer. Solution: It converges for s < and diverges otherwise. To see this use the integral test,
More information