Tutorial on Fourier and Wavelet Transformations in Geometric Algebra
|
|
- Russell Brendan Hancock
- 6 years ago
- Views:
Transcription
1 Tutorial on Fourier and Wavelet Transformations in Geometric Algebra E. Hitzer Department of Applied Physics University of Fukui Japan 17. August 2008, AGACSE 3 Hotel Kloster Nimbschen, Grimma, Leipzig, Germany
2 Acknowledgements I do believe in God s creative power in having brought it all into being in the first place, I find that studying the natural world is an opportunity to observe the majesty, the elegance, the intricacy of God s creation. Francis Collins, Director of the US National Human Genome Research Institute, in Time Magazine, 5 Nov I thank my wife, my children, my parents. B. Mawardi, G. Sommer, D. Hildenbrand, H. Li, R. Ablamowicz Organizers of AGACSE 3, Leipzig 2008.
3 CONTENTS 1 Geometric Calculus Multivector Functions 2 Clifford GA Fourier Transformations (FT) Definition and properties of GA FTs Applications of GA FTs, discrete and fast versions 3 Quaternion Fourier Transformations (QFT) Basic facts about Quaternions and definition of QFT Properties of QFT and applications 4 Spacetime Fourier Transformation (SFT) GA of spacetime and definition of SFT Applications of SFT 5 Clifford GA Wavelets GA wavelet basics and GA wavelet transformation GA wavelets and uncertainty limits, example of GA Gabor wavelets 6 Conclusion
4 Multivector Functions Multivectors, blades, reverse, scalar product Multivector M G p,q = Cl p,q, p + q = n, has k-vector parts (0 k n) scalar, vector, bi-vector,..., pseudoscalar M = A M A e A = M + M 1 + M M n, (1) blade index A {0, 1, 2, 3, 12, 23, 31, 123,..., n}, M A R. Reverse of M G p,q M = n k=0 replaces complex conjugation/quaternion conjugation ( 1) k(k 1) 2 M k. (2) The scalar product of two multivectors M, Ñ Gp,q is defined as M Ñ = MÑ = MÑ 0 (3) For M, Ñ Gn = G n,0 we get M Ñ = A M AN A.
5 Multivector Functions Modulus, blade subspace, pseudoscalar The modulus M of a multivector M G n = G n,0 is defined as M 2 = M M = n MA 2. (4) A=1 For n = 2(mod 4), n = 3(mod 4) pseudoscalar i n = e 1 e 2... e n A blade B describes a vector subspace Its dual blade i 2 n = 1. (5) V B = {x R p,q x B = 0}. B = Bi 1 n describes the complimentary vector subspace V B. pseudoscalar i n commutes for n = 3(mod 4) i n M = M i n, M G n,0.
6 Multivector Functions Multivector functions Multivector valued function f : R p,q G p,q = Cl p,q, p + q = n, has 2 n blade components n f(x) = f A (x)e A. (6) A=1 We define the inner product of R n Cl n,0 functions f, g by (f, g) = f(x) g(x) d n x = e A ẽ B f A (x)g B (x) d R n A,B R n x, (7) n and the L 2 (R n ; Cl n,0 )-norm f 2 = (f, f) = f(x) 2 d n x. (8) R n
7 Multivector Functions Vector Differential & Vector Derivative Vector differential of f defined (any const. a R p,q ) NB: a scalar. Vector derivative can be expanded as f(x + ɛa) f(x) a f(x) = lim, (9) ɛ 0 ɛ = n k=1 both coordinate independent! e k k, k = e k = x k, (10)
8 Multivector Functions Geometric Calculus Examples MV Functions f 1 = x f 2 = x 2 f 3 = x f 4 = x A k, 0 k n f 5 = log r r = x x 0 r = r V-Differentials a f 1 = a a f 2 = 2a x a f 3 = a x x a f 4 = a A k a f 5 = a r r 2 V-Derivatives f 1 = 3, (n = 3) f 2 = 2x f 3 = x/ x f 4 = k A k f 5 = r 1 = r/r 2. References 1) Hestenes & Sobczyk, Clifford Algebra to Geometric Calculus, ) Hitzer, Vector Differential Calculus, 2002.
9 Multivector Functions Geometric Calculus Examples MV Functions f 1 = x f 2 = x 2 f 3 = x f 4 = x A k, 0 k n f 5 = log r r = x x 0 r = r V-Differentials a f 1 = a a f 2 = 2a x a f 3 = a x x a f 4 = a A k a f 5 = a r r 2 V-Derivatives f 1 = 3, (n = 3) f 2 = 2x f 3 = x/ x f 4 = k A k f 5 = r 1 = r/r 2. References 1) Hestenes & Sobczyk, Clifford Algebra to Geometric Calculus, ) Hitzer, Vector Differential Calculus, 2002.
10 Multivector Functions Geometric Calculus Examples MV Functions f 1 = x f 2 = x 2 f 3 = x f 4 = x A k, 0 k n f 5 = log r r = x x 0 r = r V-Differentials a f 1 = a a f 2 = 2a x a f 3 = a x x a f 4 = a A k a f 5 = a r r 2 V-Derivatives f 1 = 3, (n = 3) f 2 = 2x f 3 = x/ x f 4 = k A k f 5 = r 1 = r/r 2. References 1) Hestenes & Sobczyk, Clifford Algebra to Geometric Calculus, ) Hitzer, Vector Differential Calculus, 2002.
11 Multivector Functions Geometric Calculus Examples MV Functions f 1 = x f 2 = x 2 f 3 = x f 4 = x A k, 0 k n f 5 = log r r = x x 0 r = r V-Differentials a f 1 = a a f 2 = 2a x a f 3 = a x x a f 4 = a A k a f 5 = a r r 2 V-Derivatives f 1 = 3, (n = 3) f 2 = 2x f 3 = x/ x f 4 = k A k f 5 = r 1 = r/r 2. References 1) Hestenes & Sobczyk, Clifford Algebra to Geometric Calculus, ) Hitzer, Vector Differential Calculus, 2002.
12 Multivector Functions Rules for Vector Differential & Derivative Derivative from differential (regard x = const., a = variable) f = a (a f) (11) Sum rules, product rules exist. Modification by non-commutativity: fġ f ġ (fg) = ( f)g + fġ = ( f)g + n e k f( k g). (12) k=1 Vector differential / derivative of f(x) = g(λ(x)), λ(x) R : a f = {a λ(x)} g λ, Example: For a = e k (1 k n) g f = ( λ) λ. (13) e k f = k f = ( k λ) f λ.
13 Classical scalar complex Fourier Transformations Definition (1D Cassical FT) For an integrable function f L 2 (R) L 1 (R), the Fourier transform of f is the function F{f}: R C = G 0,1 F{f}(ω) = f(x) e iωx dx, (14) R i imaginary unit: i 2 = 1. Inverse FT f(x) = F 1 [F{f}(ω)] = 1 F{f}(ω) e iωx dω. (15) 2π R
14 Examples of continuous and discrete 1D FT [Images: Wikipedia]
15 Definition and properties of GA FTs Geometric Algebra Fourier Transformation (GA FT) complex geometric complex unit i Geometric roots of 1, e.g. pseudoscalars i n, n = 2, 3(mod 4). complex f multivector function f L 2 (R n, G n) Definition (with pseudoscalars in dims. 2,3,6,7,10,11,...) The GA Fourier transform F{f}: R n G n = Cl n,0, n = 2, 3(mod 4) is given by F{f}(ω) = f(x) e inω x d n x, (16) R n for multivector functions f: R n G n. NB: Also possible for Cl 0,n, n = 1, 2(mod 4). Possible applications: dimension specific transformations for actual data and signals, with desired subspace structure. Inversion of these GA FTs f(x) = F 1 [F{f}(ω)] = 1 (2π) n R 3 F{f}(ω) e inω x d n ω. (17)
16 Definition and properties of GA FTs Overview of complex, Clifford, quaternion, spacetime FTs Complex Fourier Transformation (FT) F C {f}(ω) = f(x) e iωx dx (18) R Clifford Geometric Algebra (GA) FT in G 3 (i i 3 ) F G3 {f}( ω) = f( x) e i3 ω x d 3 x (19) R 3 GA FT in G n,0, G 0,n (i i n) F Gn {f}(ω) = f(x) e inω x d n x (20) R n QFT in H (i i, j) F H {f}(u) = ˆf(u) = e ixu f(x) e jyv d 2 x (21) R 2 SFT in spacetime algebra G 3,1 (i e t, j i 3 ) F ST A {f}(u) = f(u) = e e t ts f(x) e i3 x u d 3 xdt (22) R 3,1
17 Definition and properties of GA FTs Properties of the GA FT in G n, n = 3 (mod 4) General properties: Linearity, shift, modulation (ω-shift), scaling, convolution, Plancherel theorem and Parseval theorem f = 1 F{f}. (23) (2π) n/2 Table: Unique properties of GA FT Multiv. Funct. f L 2 (R n, G n), a R n, m N. Property Multiv. Funct. GA FT Vec. diff. (a ) m f(x) (a ω) m F{f}(ω) i m n Moments of x (a x) m f(x) i m n (a ω)m F{f}(ω) x m f(x) i m n m ω F{f}(ω) Vec. deriv. f(x)x m m f(x) n i m n F{f}(ω) m ω i m ω m F{f}(ω) dual moments f(x) m i m n F{f}(ω) ωm NB: Standard complex FT only gives coordinate moments and partial derivative properties.
18 Definition and properties of GA FTs Properties of the GA FT in G n, n = 2 (mod 4) Property Multiv. Funct. CFT Left lin. αf(x)+β g(x) αf{f}(ω)+ βf{g}(ω) x-shift f(x a) F{f}(ω) e inω a ω-shift f(x) e inω 0 x F{f}(ω ω 0 ) 1 Scaling f(ax) a n F{f}( ω a ) Convolution (f g)(x) F{f}( ω) F{g odd }(ω) +F{f}(ω) F{g even}(ω) Vec. diff. (a ) m f(x) (a ω) m F{f}(ω) i m n (a x) m f(x) (a ω) m F{f}(ω) i m n Moments of x x m f(x) m ω F{f}(ω) im n Vec. deriv. f(x) x m m f(x) F{f}(( 1) m ω) m ω i m n ω m F{f}(ω) i m n f(x) m F{f}(( 1) m ω) ω m i m n Plancherel R n f 1 (x) f 2 (x) d n 1 x (2π) n R n F{f 1 }(ω) F{f2 }(ω) d n ω 1 sc. Parseval f (2π) n/2 F{f}
19 Definition and properties of GA FTs Applications Sampling Representation of shape Uncertainty LSI filters (smoothing, edge detection) Signal analysis Image processing Fast (multi)vector pattern matching Visual flow analysis (Multi)vector field analysis
20 Definition and properties of GA FTs Relation with complex FT, Discrete GA FT, Fast GA FT Example of multivector signal decomposition f = [f 0 + f 123 i 3 ] + [f 1 + f 23 i 3 ]e 1 + [f 2 + f 31 i 3 ]e 2 + [f 3 + f 12 i 3 ]e 3 (24) corresponds to 4 complex signals. Implementation of GA FT by 4 complex FTs F G3 {f} = F[f 0 + f 123 i 3 ] + F[f 1 + f 23 i 3 ]e 1 + F[f 2 + f 31 i 3 ]e 2 + F[f 3 + f 12 i 3 ]e 3 (25) This form of the GA FT leads to the discrete GA FT using 4 discrete complex FTs. 4 fast FTs (FFT) can be used to implement fast GA FT.
21 Definition and properties of GA FTs Application: Fourier Descriptor Repr. of Shape (B. Rosenhahn et al [6]) Set of 3D sample points is projected along e 1, e 2, e 3 : f l n 1,n 2, l = 1, 2, 3. Use of rotors for 2D discrete rotor FT ( N 1,2 n 1,2, k 1,2 N 1,2, N = 2N + 1) R k 1,n 1 1,l = exp( 2π N 1 k 1 n 1 e l i 3 ), R k 2,n 2 2,l =... (26) Sample points f l n 1,n 2 give surface phase vectors p l k 1,k 2 = 1 N 1 N 2 Surface estimation 3 F (Φ 1, Φ 2 ) = f l (Φ 1, Φ 2 )e l = l=1 n 1,n2 f l n 1,n 2 Rk 1,n 1 1,l 2,l e l (27) R k 2,n 2 3 R k1,φ1 1,l R k 2,Φ 2 2,l p l k 1,k Rk 1,Φ 1 2 1,l R k 2,Φ 2 2,l (28) l=1 k 1,k 2
22 Definition and properties of GA FTs Application: Fourier Descriptor Repr. of Shape (B. Rosenhahn et al [6])
23 Definition and properties of GA FTs GA FT of Convolution Convolution The GA FT of the convolution of f(x) and g(x) (f g)(x) = f(y)g(x y) d 3 y, (29) R 3 equals for n = 3(mod 4) the product of the GA FTs of f(x) and g(x): F{f g}(ω) = F{f}(ω)F{g}(ω). (30) Convolution for n = 2(mod 4) For n = 2(mod 4) we get due to non-commutativity of the pseudoscalar i n G n i n M M i n, for M G n. (31) that F{f g}(ω) = F{f}( ω)f{g odd }(ω) + F{f}(ω)F{g even}(ω) (32)
24 Definition and properties of GA FTs GA FT of Correlation Correlation The GA FT of the correlation of f(x) and g(x) (f g)(x) = f(y)g(x + y) d 3 y, (33) R 3 equals for n = 3(mod 4) the product of the GA FTs of f (x) = f( x) and g(x): F{f g}(ω) = F{f }(ω)f{g}(ω). (34) Correlation for n = 2(mod 4) For n = 2(mod 4) we get due to non-commutativity of the pseudoscalar i n G n i n M M i n, for M G n. (35) that F{f g}(ω) = F{f }( ω)f{g odd }(ω) + F{f }(ω)f{g even}(ω) (36)
25 Definition and properties of GA FTs Application: Vector Pattern Analysis (J. Ebling, G. Scheuermann [7, 8]) rotation pattern and discrete GA FT black: vector components, red: bivector components (shown by normal vectors)
26 Definition and properties of GA FTs Application: Vector Pattern Matching (J. Ebling, G. Scheuermann [7, 8]) Convolution of = 3 3 (red), 5 3 (yellow), 8 3 (green) rotation patterns with gas furnace chamber flow field.
27 Definition and properties of GA FTs Uncertainty Principle in Physics (Image: Wikipedia) Optical Measurement Heisenberg Uncertainty Princ. Standard deviations x, p of position x, momentum p measurements x p 2, f... wave function. f x, F{f} p.
28 Definition and properties of GA FTs Application of GAFT: Uncertainty Principles Directional Uncertainty Principle [Hitzer&Mawardi, Proc. ICCA7] Multivect. funct. f L 2 (R n, G n), n = 2, 3 (mod 4) GA FT F{f}(ω). arbitrary const. vectors a, b R n (a x) 2 f(x) 2 d n x R n 1 f 2 1 Rn (2π) n f 2 (b ω) 2 F{f}(ω) 2 d n ω Equality (minimal bound) for optimal Gaussian multivector functions C 0 G n arbitrary const. multivector, 0 < k R. Uncertainty principle (direction independent) (a b)2 4 f(x) = C 0 e k x2, (37) For f L 2 (R n, G n), n = 2, 3 (mod 4) we obtain from the dir. UP that 1 f 2 x 2 f(x) 2 d n 1 x R n (2π) n f 2 ω 2 F{f}(ω) 2 d n ω n R n 4. (38)
29 Basic facts about Quaternions and definition of QFT Quaternions Gauss, Rodrigues and Hamilton s 4D quaternion algebra H over R with 3 imaginary units: ij = ji = k, jk = kj = i, ki = ik = j, i 2 = j 2 = k 2 = ijk = 1. (39) Quaternion q = q r + q i i + q j j + q k k H, q r, q i, q j, q k R (40) has quaternion conjugate (reversion in Cl + 3,0 ) q = q r q i i q j j q k k, (41) Leads to norm of q H q = q q = qr 2 + q2 i + q2 j + q2 k, pq = p q. (42) Scalar part Sc(q) = q r = 1 (q + q). (43) 2
30 Basic facts about Quaternions and definition of QFT ± Split of quaternions [Hitzer, AACA 2007] Convenient split of quaternions q = q + + q, q ± = 1 (q ± iqj). (44) 2 Explicitly in real components q r, q i, q j, q k R using (39): q ± = {q r ± q k + i(q i q j )} 1 ± k 2 Consequence: modulus identity Property: Scalar part of mixed split product = 1 ± k {q r ± q k + j(q j q i )}. (45) 2 q 2 = q 2 + q + 2. (46) Given two quaternions p, q and applying the ± split we get zero for the scalar part of the mixed products Sc(p + q ) = 0, Sc(p q + ) = 0. (47)
31 Basic facts about Quaternions and definition of QFT Definition of quaternion FT (QFT) F q{f}(ω) = ˆf(ω) = e ix 1ω 1 f(x) e jx 2ω 2 d 2 x. R 2 (48) Linearity leads to F q{f}(ω) = F q{f + f + }(ω) = F q{f }(ω) + F q{f + }(ω). (49) Remark: Other variations exist (see later). Simple complex forms for QFT of f ± The QFT of f ± split parts of quaternion function f L 2 (R 2, H) have simple complex forms ˆf ± = f ± e j(x 2ω 2 x 1 ω 1 ) d 2 x = R 2 e i(x 1ω 1 x 2 ω 2 ) f ± d 2 x. R 2 (50) Free of coordinates: ˆf = f e j x ω d 2 x, R 2 ˆf+ = f + e j x (U e ω) 1 d 2 x, R 2 (51) the reflection Ue 1 ω changes component ω 1 ω 1.
32 Basic facts about Quaternions and definition of QFT 2D complex FT and QFT (Images: T. Bülow, thesis [4])
33 Properties of QFT and applications Applications of QFT, discrete and fast versions partial differential systems colour image processing filtering disparity estimation (two images differ by local translations) texture segmentation wide ranging higher dimensional generalizations
34 Properties of QFT and applications Discrete and fast QFT: Pei, Ding, Chang (2001) [13] Discrete QFT F DQF T {f}(ω) = M 1 N 1 x 1 =0 x 2 =0 e ix 1ω 1 /M f(x) e jx 2ω 2 /N. (52) Fast QFT The simple complex forms ˆf = ˆf + ˆf +, = f e j x ω d 2 x, + f + e j x (U e ω) 1 d 2 x, (53) R 2 R 2 show that the QFT can be implemented with 2 complex M N 2D discrete FTs. This requires 2MN log 2 MN (54) real number multiplications.
35 Properties of QFT and applications Application to image processing Application: Image transformations The QFT of a quaternion function f L 2 (R 2 ; H) with a GL(R 2 ) transformation A (stretches, reflections, rotations) of its vector argument x is f(ax)(u) = det A 1 { ˆf (A 1 u) + ˆf + (Ue 1 A 1 Ue 1 u) }, (55) where A 1 is the adjoint of the inverse of A.
36 Properties of QFT and applications Properties of QFT Standard properties: Linearity, shift, modulation, dilation, Plancherel theorem and Parseval theorem (signal energy) f = 1 Fq{f}. (56) 2π Table: Special QFT properties. Quat. Funct. f L 2 (R 2 ; H), with x, u R 2, m, n N 0. Easy with correct order! Property Quat. Funct. QFT m+n Part. deriv. x m y n f(x) (iu) m ˆf(u)(jv) n Moments of x, y x m y n f(x) i m m+n u m v n j n Powers of i, j i m f(x) j n i m ˆf(u) j n ± split and QFT commute: F q{f ± } = F q{f} ±.
37 Properties of QFT and applications Properties of QFT Integration of parts With the vector differential a = a a 2 2, with arbitrary constant a R 2, g, h L 2 (R 2, H) [ g(x)[a h(x)]d 2 x = R 2 g(x)h(x)dx R ] a x= a x= R 2 [a g(x)]h(x)d 2 x. (57) Modulus identities Due to q 2 = q 2 + q + 2 we get for f : R 2 H the following identities f(x) 2 = f (x) 2 + f + (x) 2, (58) F q{f}(ω) 2 = F q{f }(ω) 2 + F q{f + }(ω) 2. (59)
38 Properties of QFT and applications Properties of QFT QFT of vector differentials Using the split f = f + f + we get the QFTs of the split parts. Let b R 2 be an arbitrary constant vector. F q{b f }(ω) = b ωf q{f }(ω) j = i b ωf q{f }(ω), (60) F q{b f + }(ω) = (b Ue 1 ω) F q{f + }(ω) j = i (b Ue 2 ω) F q{f + }(ω), (61) F q{(ue 1 b )f + }(ω) = b ω F q{f + }(ω) j, (62) F q{(ue 2 b )f + }(ω) = i b ω F q{f + }(ω). (63)
39 Properties of QFT and applications Gaussian quaternion filter (GF) h(x) = g(x) e i2πx 1ω 01 e j2πx 2ω 02 2D Gaussian amplitude: g(x) = ( 1 e 1 x 2 ) 1 2 σ 2 + x2 2 (2π) 3 1 σ σ 1 σ 2
40 Properties of QFT and applications Application of QFT to texture segmentation (T. Bülow, thesis [4])
41 Properties of QFT and applications Application of QFT to disparity estimation (T. Bülow, thesis [4]) QFT method overcomes directional limitations of complex methods.
42 Properties of QFT and applications Uncertainty Principles for QFT (componentwise and directional) Right sided quaternion Fourier transform (QFT) in H F r{f}(ω) = f(x) e ix 1ω 1 e jx 2ω 2 d 2 x (64) R 2 Componentwise (k = 1, 2) uncertainty [Mawardi&Hitzer, 2008] 1 f 2 x 2 R 2 k f(x) 2 d 2 1 x (2π) 2 f 2 ω 2 R 2 k f(ω) 2 d 2 ω 1 4. (65) Equality holds if and only if f is a Gaussian quaternion function. Directional QFT Uncertainty Principle [Hitzer, 2008] For two arbitrary constant vectors a, b R 2 (selecting two directions), and f L 2 (R 2, H), x 1/2 f L 2 (R 2, H), b = b 1 e 1 + b 2 e 2 we obtain (a x) 2 f(x) 2 d 2 x (b ω) 2 F q{f}(ω) 2 d 2 ω R 2 R 2 (2π)2 4 [ (a b) 2 f 4 + (a b ) 2 f + 4], (66)
43 GA of spacetime and definition of SFT Motion in time: Video sequences, flow fields,... GA of spacetime (STA) {e t, e 1, e 2, e 3 }, e 2 1 = e 2 2 = e 2 3 = 1. (67) e 2 t = 1, i 3 = e 1 e 2 e 3, i 2 3 = 1, ist = ete 1e 2 e 3, i 2 st = 1, (68) with time vector e t, 3D space volume i 3, hyper volume of spacetime i st. Isomorphism: Volume-time subalgebra to Quaternions {1, e t, i 3, i st} {1, i, j, k} (69) NB: Does not work for G 2,1 or G 4,1, but works for conformal model of spacetime G 5,2. Split Spacetime Split (e t = eti 1 st = e ti st = e te ti 3 = i 3 ) f ± = 1 2 (f ± etfe t ) (70) Time direction e t determines complimentary 3D space i 3 as well!
44 GA of spacetime and definition of SFT QFT Spacetime Fourier transform (SFT): f F SF T {f} with F SF T {f}(u) = f(u) = e e t ts f(x) e i3 x u d 4 x. (71) R 3,1 spacetime vectors x = te t + x R 3,1, x = xe 1 + ye 2 + ze 3 R 3 spacetime volume d 4 x = dtdxdydz spacetime frequency vectors u = se t + u R 3,1, u = ue 1 + ve 2 + we 3 R 3 maps 16D spacetime algebra functions f : R 3,1 G 3,1 = Cl 3,1 to 16D spacetime spectrum functions f : R 3,1 Cl 3,1 Part f(x) e i 3 x u d 3 x... GA FT of G 3. Reversion replaced by principal involution (reversion and e t e t).
45 Applications of SFT Application: Multivector Wave Packets in Space and Time SFT = Sum of Right/Left Propagating Multivector Wave Packets f = f + + f = f + e i 3( x u ts ) d 4 x + f e i 3( x u+ ts ) d 4 x. (72) R 3,1 R 3,1 Directional 4D Spacetime Uncertainty for Multivector Wave Packets For two arbitrary constant spacetime vectors a, b R 3,1 (selecting two directions), and f L 2 (R 3,1, G 3,1 ), x 1/2 f L 2 (R 3,1, G 3,1 ) we obtain (a tt a x) 2 f(x) 2 d 4 x (b tω t b ω) 2 F{f}(ω) 2 d 4 ω R 3,1 R 3,1 (2π)4 4 [ (a tb t a b) 2 f 4 + (a tb t + a b) 2 f + 4],
46 Applications of SFT Spacetime signal transformations Application: Space & Time Signal Transformations The SFT of a spacetime function f L 2 (R 3,1 ; G 3,1 ) with a GL(R 3,1 ) transformation A (stretches, reflections, rotations, acceleration, boost) of its spacetime vector argument x is F SF T {f(ax)}(u) = det A 1 { f + (Ue t A 1 Ue t u) + f (A 1 u) }, (73) where A 1 is the adjoint of the inverse of A.
47 GA wavelet basics and GA wavelet transformation GA Wavelet Basics [Mawardi&Hitzer, 2007] Transformation group SIM(3) We apply 3D (elements of G = SIM(3)) translations, scaling and rotations to a mother wavelet ψ : R 3 G 3 = Cl 3,0 Wavelet admissibility (condition on ψ) ψ(x) ψ a,θ,b (x) = 1 ψ(r 1 a3/2 θ ( x b )) (74) a ψ(ω) C ψ = ψ(ω) R 3 ω 3 d 3 ω. (75) is an invertible multivector constant and finite at a.e. ω R 3, GA FT ψ = F{ψ}.
48 GA wavelet basics and GA wavelet transformation GA Wavelet Transformation Definition (GA wavelet transformation) T ψ : L 2 (R 3 ; Cl 3,0 ) L 2 (G; Cl 3,0 ) f T ψ f(a, θ, b) = f(x) ψa,θ,b (x) d 3 x, (76) R 3 ψ L 2 (R 3 ; Cl 3,0 )... GA mother wavelet, f... image, multivector field, data, etc. Properties: Locality (different from global FTs), left linearity; translation, dilation, rotation covariance. Inversion of GA wavelet transform Any f L 2 (R 3 ; Cl 3,0 ) can be decomposed as (C 1 ψ... admissibility) f(x) = T ψ f(a, b, θ) ψ a,θ,b C 1 ψ dµd3 b, (77) G with left Haar measure on G = SIM(3): dλ(a, θ, b) = dµ(a, θ)d 3 b
49 GA wavelet basics and GA wavelet transformation Applications (expected) of GA wavelets multi-dimensional image/signal processing JPEG local spherical harmonics (lighting) geological exploration seismology local multivector pattern matching and the like
50 GA wavelet basics and GA wavelet transformation Notation: Inner products and norms We define the inner product of f, g : G Cl 3,0 by (f, g) = f(a, θ, b) g(a, θ, b) dλ(a, θ, b), (78) G and the L 2 (G; Cl 3,0 )-norm f 2 = (f, f) = f(a, θ, b) 2 dλ. (79) G
51 GA wavelets and uncertainty limits, example of GA Gabor wavelets GA Wavelet Uncertainty Generalized GA wavelet uncertainty principle Let ψ be a Clifford algebra wavelet that satisfies the admissibility condition. Then for every f L 2 (R 3 ; Cl 3,0 ), the following inequality holds bt ψ f(a, θ, b) 2 L 2 (G;Cl 3,0 ) C ψ ( ω ˆf, ω ˆf) L 2 (R 3 ;Cl 3,0 ) 3(2π)3 4 [ C ψ (f, f) L 2 (R 3 ;Cl 3,0 )] 2. (80)
52 GA wavelets and uncertainty limits, example of GA Gabor wavelets Uncertainty principle for scalar admissibility constant For scalar admissibility constant C ψ we get Uncertainty principle for GA wavelet Let ψ be a Clifford algebra wavelet with scalar admissibility constant. Then for every f L 2 (R 3 ; Cl 3,0 ), the following inequality holds b T ψ f(a, θ, b) 2 L 2 (G;Cl 3,0 ) ω ˆf 2 L 2 (R 3 ;Cl 3,0 ) 3C (2π) 3 ψ f 4 4 L 2 (R 3 ;Cl 3,0 ). (81) NB1: This shows indeed, that the previous inequality (80) represents a multivector generalization of the uncertainty principle of inequality (81) for Clifford wavelets with scalar admissibility constant. NB2: Compare with (direction independent) uncertainty principle for GA FT.
53 GA wavelets and uncertainty limits, example of GA Gabor wavelets GA Gabor Wavelets Example: GA Gabor Wavelets ψ c (x) = g(x) e i 3ω 0 x e 2 1 (σ2 1 u2 0 +σ2 2 u2 0 +σ2 3 w2 0 }{{} ) constant 3D Gaussian: g(x) = ( 1 e 1 x σ (2π) σ 1 σ 2 σ 3 ) + x2 2 σ x2 3 σ 3 2. Uncertainty principle for GA Gabor wavelet Let ψ c be an admissible GA Gabor wavelet. Assume f L 2 (R 3 ; Cl 3,0 ), then the following inequality holds b T ψ c f(a, θ, b) 2 L 2 (G;Cl 3,0 ) ω ˆf 2 L 2 (R 3 ;Cl 3,0 ) 3C (2π) 3 ψ c f 4. (82) 4
54 Conclusion Geometric calculus allows to construct a variety of GA Fourier transformations. Especially quaternion Fourier transformations are well studied. We studied the space time Fourier transformation giving rise to a multivector wave packet decomposition. Discrete and fast versions exist. Uncertainty limits established (principle bounds to accuracy). A diverse range of applications exists. For the future the application of local GA wavelets may overcome the limitations of global GA Fourier transformations.
55 D. Hestenes, G. Sobczyk, Clifford Algebra to Geometric Calculus, Kluwer, Dordrecht, E. Hitzer, Multivector Differential Calculus, AACA, 12(2), pp (2002). T. A. Ell, Quat.-Four. Transf. for Analysis of Two-dim. Linear Time-Invariant Partial Diff. Systems, Proc. 32nd IEEE Conf. on Decision and Control, Dec. 1993, pp T. Bülow, Hypercomplex Spectral Signal Representations for the Processing and Analysis of Images, Ph. D. thesis, University of Kiel (1999). T. Bülow, M. Felsberg and G. Sommer, Non-commutative Hypercomplex Four. Transf. of Multidim. Signals, in G. Sommer (ed.), Geom. Comp. with Cliff. Alg., Theor. Found. and Appl. in Comp. Vis. and Rob., Springer, (2001), Rosenhahn B. and Sommer G., Pose estim. of free-form obj., Proc. of the Eur. Conf. on Comp. Vis. ECCV 04, Part I, T. Pajdla and J. Matas (Eds.), Springer, Berlin, LNCS 3021, J. Ebling, G. Scheuermann, Clifford Fourier Transform on Vector Fields, IEEE Trans. on Vis. and Comp. Graph., 11 (4), (July/August 2005). J. Ebling, G. Scheuermann, Clifford Convolution And Pattern Matching On Vector Fields, Proc 14th IEEE Vis (VIS 03), p.26, (October 22-24, 2003)
56 B. Mawardi, E. Hitzer. Clifford Fourier Transform and Uncertainty Principle for the Clifford Geometric Algebra Cl 3,0, Adv. App. Cliff. Alg., 16(1), pp , E. Hitzer, B. Mawardi, Cliff. FT on Multivector Fields and Unc. Princ. for Dimensions n = 2 (mod 4) and n = 3 (mod 4), Proc. of ICCA7, Toulouse, B. Mawardi, E. Hitzer, A. Hayashi, R. Ashino, An uncertainty principle for quaternion FT, in press, CAMWA, E. Hitzer, Quaternion Fourier Transform on Quaternion Fields and Generalizations, Adv. in App. Cliff. Alg., 17, pp , S.C. Pei, J.J. Ding, J.H. Chang, Efficient Implementation of Quaternion Fourier Transform, Convolution, and Correlation by 2-D Complex FFT, IEEE TRANS. on SIG. PROC., Vol. 49, No. 11, p (2001). B. Mawardi, E. Hitzer, Cliff. Alg. Cl 3,0 -valued Wavelet Transf., Cliff. Wavelet Uncertainty Inequality and Cliff. Gabor Wavelets, IJWMIP, F. Brackx, R. Delanghe, and F. Sommen, Clifford Analysis, Vol. 76 of Research Notes in Mathematics, Pitman Advanced Publishing Program, C. Li, A. McIntosh and T. Qian, Cliff. Algs., Fourier Transf. and Sing. Conv. Operators On Lipschitz Surfaces, Rev. Mat. Iberoam., 10 (3), (1994), Other publications:
57 News about Clifford Geometric Algebra GA-Net Electronic newsletter: GA-Net Updates (blog) Immediate access to latest GA news: Soli Deo Gloria J. S. Bach (Leipzig)
An Uncertainty Principle for Quaternion Fourier Transform
An Uncertainty Principle for Quaternion Fourier Transform Mawardi Bahri a, Eckhard S. M. Hitzer a Akihisa Hayashi a Ryuichi Ashino b, a Department of Applied Physics, University of Fukui, Fukui 9-857,
More informationLOGARITHMIC UNCERTAINTY PRINCIPLE FOR QUATERNION LINEAR CANONICAL TRANSFORM
Proceedings of the 06 International Conference on Wavelet Analysis and Pattern Recognition Jeju South Korea 0-3 July LOGARITHMIC UNCERTAINTY PRINCIPLE FOR QUATERNION LINEAR CANONICAL TRANSFORM MAWARDI
More informationTwo-Dimensional Clifford Windowed Fourier Transform
Two-Dimensional Clifford Windowed Fourier Transform Mawardi Bahri, Eckhard M. S. Hitzer and Sriwulan Adji Abstract Recently several generalizations to higher dimension of the classical Fourier transform
More informationProduct Theorem for Quaternion Fourier Transform
Int. Journal of Math. Analysis, Vol. 8, 204, no. 2, 8-87 HIKAI Ltd, www.m-hikari.com http://dx.doi.org/0.2988/ijma.204.3290 Product Theorem for Quaternion ourier Transform Mawardi Bahri Department of Mathematics,
More informationDemystification of the Geometric Fourier Transforms
Demystification of the Geometric Fourier Transforms Roxana Buack, Gerik Scheuermann and Eckhard Hitzer Universität Leipzig, Institut für Informatik, Abteilung für Bild- und Signalverarbeitung, Augustuplatz
More informationThe orthogonal planes split of quaternions and its relation to quaternion geometry of rotations 1
The orthogonal planes split of quaternions and its relation to quaternion geometry of rotations 1 Eckhard Hitzer Osawa 3-10-, Mitaka 181-8585, International Christian University, Japan E-mail: hitzer@icu.ac.jp
More informationThe orthogonal planes split of quaternions and its relation to quaternion geometry of rotations
Home Search Collections Journals About Contact us My IOPscience The orthogonal planes split of quaternions and its relation to quaternion geometry of rotations This content has been downloaded from IOPscience.
More informationOperator Exponentials for the Clifford Fourier Transform on Multivector Fields
Operator Exponentials for the Clifford Fourier Transform on Multivector Fields David Eelbode and Echard Hitzer Mathematics Subject Classification 000). Primary 4A38; Secondary 15A66,34L40. Keywords. Fourier
More informationInternational Journal of Wavelets, Multiresolution and Information Processing c World Scientific Publishing Company
October 6, 014 1:7 WSPC/WS-IJWMIP QTF-ijwmip revf International Journal of Wavelets, Multiresolution and Information Processing c World Scientific Publishing Company CONTINUOUS QUATENION FOUIE AND WAVELET
More informationConnecting spatial and frequency domains for the quaternion Fourier transform
Connecting spatial and frequency domains for the quaternion Fourier transform Ghent University (joint work with N. De Schepper, T. Ell, K. Rubrecht and S. Sangwine) MOIMA, Hannover, June, 2016 Direct formulas
More informationShort-time Fourier transform for quaternionic signals
Short-time Fourier transform for quaternionic signals Joint work with Y. Fu and U. Kähler P. Cerejeiras Departamento de Matemática Universidade de Aveiro pceres@ua.pt New Trends and Directions in Harmonic
More informationResearch Article A Simplified Proof of Uncertainty Principle for Quaternion Linear Canonical Transform
Abstract and Applied Analysis Volume 06, Article ID 5874930, pages http://doiorg/055/06/5874930 Research Article A Simplified Proof of Uncertainty Principle for Quaternion Linear Canonical Transform Mawardi
More informationAnalyzing Real Vector Fields with Clifford Convolution and Clifford Fourier Transform
Analyzing Real Vector Fields with Clifford Convolution and Clifford Fourier Transform Wieland Reich and Gerik Scheuermann Abstract Post-processing in computational fluid dynamics as well as processing
More informationKeywords: Vector eld, Cliord's geometric algebra, Geometric calculus
Vector eld computations in Cliord's geometric algebra * Eckhard Hitzer (International Christian University, Japan), Roxana Bujack and Gerik Scheuermann (Leipzig University, Germany) Abstract Exactly 125
More informationVector Field Processing with Clifford Convolution and Clifford Fourier Transform
Vector Field Processing with Clifford Convolution and Clifford Fourier Transform Wieland Reich ¹, Gerik ¹ ¹ Computer Science Institute, University of Leipzig Outline The problem: CFD data and PIV measurements
More informationarxiv:math/ v1 [math.dg] 6 Jul 1998
arxiv:math/9807024v [math.dg] 6 Jul 998 The Fundamental Theorem of Geometric Calculus via a Generalized Riemann Integral Alan Macdonald Department of Mathematics Luther College, Decorah, IA 520, U.S.A.
More informationLocal textures parameters of images in space domain obtained with different analytical approaches
European Journal of Mathematics and Computer Science Vol. 4 No. 1, 017 Local textures parameters of images in space domain obtained with different analytical approaches Robert Goutte University of Lyon,
More informationThe Fueter Theorem and Dirac symmetries
The Fueter Theorem and Dirac symmetries David Eelbode Departement of Mathematics and Computer Science University of Antwerp (partially joint work with V. Souček and P. Van Lancker) General overview of
More informationA Novel Approach to the 2D Analytic Signal? Thomas Bulow and Gerald Sommer. Christian{Albrechts{Universitat zu Kiel
A Novel Approach to the 2D Analytic Signal? Thomas Bulow and Gerald Sommer Christian{Albrechts{Universitat zu Kiel Institute of Computer Science, Cognitive Systems Preuerstrae 1{9, 24105 Kiel Tel:+49 431
More informationAn Introduction to Geometric Algebra
An Introduction to Geometric Algebra Alan Bromborsky Army Research Lab (Retired) brombo@comcast.net November 2, 2005 Typeset by FoilTEX History Geometric algebra is the Clifford algebra of a finite dimensional
More informationA General Geometric Fourier Transform Convolution Theorem
A General Geometric Fourier Transform Convolution Theorem Roxana Bujack, Gerik Scheuermann and Eckhard Hitzer Abstract. The large variety of Fourier transforms in geometric algebras inspired the straight
More informationLie Algebra of Unit Tangent Bundle in Minkowski 3-Space
INTERNATIONAL ELECTRONIC JOURNAL OF GEOMETRY VOLUME 12 NO. 1 PAGE 1 (2019) Lie Algebra of Unit Tangent Bundle in Minkowski 3-Space Murat Bekar (Communicated by Levent Kula ) ABSTRACT In this paper, a one-to-one
More informationClifford Fourier Transform for Color Image Processing
Clifford Fourier Transform for Color Image Processing Thomas Batard, Michel Berthier and Christophe Saint-Jean Abstract The aim of this paper is to define a Clifford Fourier transform that is suitable
More informationNew Developments in Clifford Fourier Transforms
New Developments in Clifford Fourier Transforms Eckhard Hitzer Abstract We show how real and complex Fourier transforms are extended to W.R. Hamiltons algebra of quaternions and to W.K. Clifford s geometric
More informationEckhard Hitzer and Stephen J. Sangwine. 1. Introduction
Construction of multivector inverse for Clifford algebras over 2m+1-dimensional vector spaces from multivector inverse for Clifford algebras over 2m-dimensional vector spaces. Eckhard Hitzer and Stephen
More informationA General Geometric Fourier Transform
A General Geometric Fourier Transform Roxana Bujack, Gerik Scheuermann and Eckhard Hitzer Abstract. The increasing demand for Fourier transforms on geometric algebras has resulted in a large variety. Here
More informationSome results on the lattice parameters of quaternionic Gabor frames
Some results on the lattice parameters of quaternionic Gabor frames S. Hartmann Abstract Gabor frames play a vital role not only modern harmonic analysis but also in several fields of applied mathematics,
More informationOld Wine in New Bottles: A new algebraic framework for computational geometry
Old Wine in New Bottles: A new algebraic framework for computational geometry 1 Introduction David Hestenes My purpose in this chapter is to introduce you to a powerful new algebraic model for Euclidean
More informationClassical Physics. SpaceTime Algebra
Classical Physics with SpaceTime Algebra David Hestenes Arizona State University x x(τ ) x 0 Santalo 2016 Objectives of this talk To introduce SpaceTime Algebra (STA) as a unified, coordinate-free mathematical
More informationAlgorithm for conversion between geometric algebra versor notation and conventional crystallographic symmetry-operation symbols
Algorithm for conversion between geometric algebra versor notation and conventional crystallographic symmetry-operation symbols Eckhard Hitzer and Christian Perwass June, 2009 Introduction This paper establishes
More informationContinuous quaternion fourier and wavelet transforms
International Journal of Wavelets, Multiresolution and Information Processing Vol., No. 4 (04) 46000 ( pages) c World Scientific Publishing Company DOI: 0.4/S0969460000 Continuous quaternion fourier and
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 informationECS 178 Course Notes QUATERNIONS
ECS 178 Course Notes QUATERNIONS Kenneth I. Joy Institute for Data Analysis and Visualization Department of Computer Science University of California, Davis Overview The quaternion number system was discovered
More informationAn Introduction to Geometric Algebra and Calculus. Alan Bromborsky Army Research Lab (Retired)
An Introduction to Geometric Algebra and Calculus Alan Bromborsky Army Research Lab Retired abrombo@verizon.net April 1, 2011 2 Introduction Geometric algebra is the Clifford algebra of a finite dimensional
More informationand multiplication as q + p = (q 0 + iq 1 + jq 2 + kq 3 ) + (p 0 + ip 1 + j p 2 + kp 3 )
Ben Witten* and Jeff Shragge, Stanford University SUMMAY Hypercomlex numbers are primarily used for pattern recognition, offer many useful applications to geophysics. Image disparity estimation is a hypercomplex,
More informationConic Sections and Meet Intersections in Geometric Algebra
Conic Sections and Meet Intersections in Geometric Algebra Eckhard M.S. Hitzer Department of Physical Engineering, University of Fukui, Japan hitzer@mech.fukui-u.ac.jp Abstract. This paper first gives
More informationMotivation towards Clifford Analysis: The classical Cauchy s Integral Theorem from a Clifford perspective
Motivation towards Clifford Analysis: The classical Cauchy s Integral Theorem from a Clifford perspective Lashi Bandara November 26, 29 Abstract Clifford Algebras generalise complex variables algebraically
More information1 Singular Value Decomposition
1 Singular Value Decomposition Factorisation of rectangular matrix (generalisation of eigenvalue concept / spectral theorem): For every matrix A C m n there exists a factorisation A = UΣV U C m m, V C
More informationGeometry and Kinematics with Uncertain Data
Geometry and Kinematics with Uncertain Data Christian Perwass, Christian Gebken, and Gerald Sommer Institut für Informatik, CAU Kiel Christian-Albrechts-Platz 4, 24118 Kiel, Germany {chp,chg,gs}@ks.informatik.uni-kiel.de
More informationArchiv der Mathematik Holomorphic approximation of L2-functions on the unit sphere in R^3
Manuscript Number: Archiv der Mathematik Holomorphic approximation of L-functions on the unit sphere in R^ --Manuscript Draft-- Full Title: Article Type: Corresponding Author: Holomorphic approximation
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 informationQUATERNIONS AND ROTATIONS
QUATERNIONS AND ROTATIONS SVANTE JANSON 1. Introduction The purpose of this note is to show some well-known relations between quaternions and the Lie groups SO(3) and SO(4) (rotations in R 3 and R 4 )
More informationSpinor Fourier Transform for Image Processing
Spinor Fourier Transform for Image Processing Thomas Batard, Michel Berthier To cite this version: Thomas Batard, Michel Berthier. Spinor Fourier Transform for Image Processing. 2013. HAL
More informationThe Quaternion Domain Fourier Transform and its Application in Mathematical Statistics
The Quaternion Domain Fourier Transform and its Application in Mathematical Statistics Mawardi Bahri, Amir Kamal Amir, Resnawati, and Chrisandi Lande Abstract Recently a generalization of the quaternion
More information1. Fourier Transform (Continuous time) A finite energy signal is a signal f(t) for which. f(t) 2 dt < Scalar product: f(t)g(t)dt
1. Fourier Transform (Continuous time) 1.1. Signals with finite energy A finite energy signal is a signal f(t) for which Scalar product: f(t) 2 dt < f(t), g(t) = 1 2π f(t)g(t)dt The Hilbert space of all
More informationGabor wavelet analysis and the fractional Hilbert transform
Gabor wavelet analysis and the fractional Hilbert transform Kunal Narayan Chaudhury and Michael Unser (presented by Dimitri Van De Ville) Biomedical Imaging Group, Ecole Polytechnique Fédérale de Lausanne
More informationThe Monogenic Wavelet Transform Sofia C. Olhede and Georgios Metikas
3426 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 57, NO. 9, SEPTEMBER 2009 The Monogenic Wavelet Transform Sofia C. Olhede and Georgios Metikas Abstract This paper extends the 1-D analytic wavelet transform
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 informationWavelets, wavelet networks and the conformal group
Wavelets, wavelet networks and the conformal group R. Vilela Mendes CMAF, University of Lisbon http://label2.ist.utl.pt/vilela/ April 2016 () April 2016 1 / 32 Contents Wavelets: Continuous and discrete
More informationMathematical Methods for Computer Science
Mathematical Methods for Computer Science Computer Laboratory Computer Science Tripos, Part IB Michaelmas Term 2016/17 Professor J. Daugman Exercise problems Fourier and related methods 15 JJ Thomson Avenue
More informationMultiresolution analysis
Multiresolution analysis Analisi multirisoluzione G. Menegaz gloria.menegaz@univr.it The Fourier kingdom CTFT Continuous time signals + jωt F( ω) = f( t) e dt + f() t = F( ω) e jωt dt The amplitude F(ω),
More informationWavelets: Theory and Applications. Somdatt Sharma
Wavelets: Theory and Applications Somdatt Sharma Department of Mathematics, Central University of Jammu, Jammu and Kashmir, India Email:somdattjammu@gmail.com Contents I 1 Representation of Functions 2
More informationarxiv:math/ v1 [math.ra] 2 Jun 2005
Quaternion Involutions arxiv:math/0506034v1 [math.ra] 2 Jun 2005 Todd A. Ell 5620 Oak View Court, Savage, MN 55378-4695, USA. Email: T.Ell@IEEE.org. Stephen J. Sangwine Department of Electronic Systems
More informationA Spherical Harmonic Expansion of the Hilbert Transform on S 2. Oliver Fleischmann Cognitive Systems Group Kiel University.
A Spherical Harmonic Expansion of the Hilbert Transform on S 2 Oliver Fleischmann Cognitive Systems Group Kiel University München 2010 Motivation Goal: Extract local structural information from signals
More informationAlgebra of Complex Vectors and Applications in Electromagnetic Theory and Quantum Mechanics
Mathematics 2015, 3, 781-842; doi:10.3390/math3030781 OPEN ACCESS mathematics ISSN 2227-7390 www.mdpi.com/journal/mathematics Article Algebra of Complex Vectors and Applications in Electromagnetic Theory
More informationIntroduction to Modern Quantum Field Theory
Department of Mathematics University of Texas at Arlington Arlington, TX USA Febuary, 2016 Recall Einstein s famous equation, E 2 = (Mc 2 ) 2 + (c p) 2, where c is the speed of light, M is the classical
More informationDecomposition algorithms in Clifford algebras
Decomposition algorithms in Clifford algebras G. Stacey Staples 1 Department of Mathematics and Statistics Southern Illinois University Edwardsville Combinatorics and Computer Algebra Fort Collins, 2015
More informationSerge Ballif January 18, 2008
ballif@math.psu.edu The Pennsylvania State University January 18, 2008 Outline Rings Division Rings Noncommutative Rings s Roots of Rings Definition A ring R is a set toger with two binary operations +
More informationY. Liu and A. Mohammed. L p (R) BOUNDEDNESS AND COMPACTNESS OF LOCALIZATION OPERATORS ASSOCIATED WITH THE STOCKWELL TRANSFORM
Rend. Sem. Mat. Univ. Pol. Torino Vol. 67, 2 (2009), 203 24 Second Conf. Pseudo-Differential Operators Y. Liu and A. Mohammed L p (R) BOUNDEDNESS AND COMPACTNESS OF LOCALIZATION OPERATORS ASSOCIATED WITH
More informationWavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing
Wavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing Qingtang Jiang Abstract This paper is about the construction of univariate wavelet bi-frames with each framelet being symmetric.
More informationHilbert Space Problems
Hilbert Space Problems Prescribed books for problems. ) Hilbert Spaces, Wavelets, Generalized Functions and Modern Quantum Mechanics by Willi-Hans Steeb Kluwer Academic Publishers, 998 ISBN -7923-523-9
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 informationAtomic decompositions of square-integrable functions
Atomic decompositions of square-integrable functions Jordy van Velthoven Abstract This report serves as a survey for the discrete expansion of square-integrable functions of one real variable on an interval
More 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 informationarxiv: v1 [math.ca] 18 Nov 2016
BEURLING S THEOREM FOR THE CLIFFORD-FOURIER TRANSFORM arxiv:1611.06017v1 [math.ca] 18 Nov 016 RIM JDAY AND JAMEL EL KAMEL Abstract. We give a generalization of Beurling s theorem for the Clifford-Fourier
More informationQuaternion Wiener Deconvolution for Noise Robust Color Image Registration
1278 IEEE SIGNAL PROCESSING LETTERS, VOL. 22, NO. 9, SEPTEMBER 2015 Quaternion Wiener Deconvolution for Noise Robust Color Image Registration Matteo Pedone, Eduardo Bayro-Corrochano, Jan Flusser, and Janne
More informationChapter 1 Numerical evaluation of Versors with Clifford Algebra
This is page 1 Printer: Opaque this Chapter 1 Numerical evaluation of Versors with Clifford Algebra Christian B. U. Perwass, Gerald Sommer 1 ABSTRACT This paper has two main parts. In the first part we
More informationGeometric Algebra 2 Quantum Theory
Geometric Algebra 2 Quantum Theory Chris Doran Astrophysics Group Cavendish Laboratory Cambridge, UK Spin Stern-Gerlach tells us that electron wavefunction contains two terms Describe state in terms of
More informationHyperbolic Conformal Geometry with Clifford Algebra 1)
MM Research Preprints, 72 83 No. 18, Dec. 1999. Beijing Hyperbolic Conformal Geometry with Clifford Algebra 1) Hongbo Li Abstract. In this paper we study hyperbolic conformal geometry following a Clifford
More informationMultivector Calculus
In: J. Math. Anal. and Appl., ol. 24, No. 2, c Academic Press (1968) 313 325. Multivector Calculus David Hestenes INTRODUCTION The object of this paper is to show how differential and integral calculus
More informationIr O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )
Section 3.2 Theorem 3.6. Let A be an m n matrix of rank r. Then r m, r n, and, by means of a finite number of elementary row and column operations, A can be transformed into the matrix ( ) Ir O D = 1 O
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 informationClifford Analysis, Homework 1
Clifford Analysis, Homework 1 November 1, 017 1 Let w v 1... v k, for vectors v j. Show that ŵ is the result of negating the vectors: ŵ ( v 1 )... ( v k ). Show that w is the result of changing the order
More informationThe following definition is fundamental.
1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic
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 informationSOME TOPICS ON WAVELETS
2 SOME TOPICS ON WAVELETS RYUICHI ASHINO 1. Introduction Let us consider music. If the music is recorded from a live broadcast onto tape, we have a signal, that is, a function f(x). The time-frequency
More informationTRANSLATION-INVARIANT CLIFFORD OPERATORS
TRANSLATION-INVARIANT CLIFFORD OPERATORS JEFF HOGAN AND ANDREW J. MORRIS Abstract. This paper is concerned with quaternion-valued functions on the plane and operators which act on such functions. In particular,
More informationIntro to harmonic analysis on groups Risi Kondor
Risi Kondor Any (sufficiently smooth) function f on the unit circle (equivalently, any 2π periodic f ) can be decomposed into a sum of sinusoidal waves f(x) = k= c n e ikx c n = 1 2π f(x) e ikx dx 2π 0
More information1 From Fourier Series to Fourier transform
Differential Equations 2 Fall 206 The Fourier Transform From Fourier Series to Fourier transform Recall: If fx is defined and piecewise smooth in the interval [, ] then we can write where fx = a n = b
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 informationThe Fourier Transform Method
The Fourier Transform Method R. C. Trinity University Partial Differential Equations April 22, 2014 Recall The Fourier transform The Fourier transform of a piecewise smooth f L 1 (R) is ˆf(ω) = F(f)(ω)
More informationMatrices over Hyperbolic Split Quaternions
Filomat 30:4 (2016), 913 920 DOI 102298/FIL1604913E Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://wwwpmfniacrs/filomat Matrices over Hyperbolic Split
More informationFourier transforms. R. C. Daileda. Partial Differential Equations April 17, Trinity University
The Fourier Transform R. C. Trinity University Partial Differential Equations April 17, 214 The Fourier series representation For periodic functions Recall: If f is a 2p-periodic (piecewise smooth) function,
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationMATH45061: SOLUTION SHEET 1 II
MATH456: SOLUTION SHEET II. The deformation gradient tensor has Cartesian components given by F IJ R I / r J ; and so F R e x, F R, F 3 R, r r r 3 F R r, F R r, F 3 R r 3, F 3 R 3 r, F 3 R 3 r, F 33 R
More informationConstruction of wavelets. Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam
Construction of wavelets Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam Contents Stability of biorthogonal wavelets. Examples on IR, (0, 1), and (0, 1) n. General domains
More informationOn the Hilbert Transform of Wavelets
On the Hilbert Transform of Wavelets Kunal Narayan Chaudhury and Michael Unser Abstract A wavelet is a localized function having a prescribed number of vanishing moments. In this correspondence, we provide
More informationLectures notes. Rheology and Fluid Dynamics
ÉC O L E P O L Y T E C H N IQ U E FÉ DÉR A L E D E L A U S A N N E Christophe Ancey Laboratoire hydraulique environnementale (LHE) École Polytechnique Fédérale de Lausanne Écublens CH-05 Lausanne Lectures
More informationDecomposition of Riesz frames and wavelets into a finite union of linearly independent sets
Decomposition of Riesz frames and wavelets into a finite union of linearly independent sets Ole Christensen, Alexander M. Lindner Abstract We characterize Riesz frames and prove that every Riesz frame
More informationFourier Series. (Com S 477/577 Notes) Yan-Bin Jia. Nov 29, 2016
Fourier Series (Com S 477/577 otes) Yan-Bin Jia ov 9, 016 1 Introduction Many functions in nature are periodic, that is, f(x+τ) = f(x), for some fixed τ, which is called the period of f. Though function
More informationChapter One. The Calderón-Zygmund Theory I: Ellipticity
Chapter One The Calderón-Zygmund Theory I: Ellipticity Our story begins with a classical situation: convolution with homogeneous, Calderón- Zygmund ( kernels on R n. Let S n 1 R n denote the unit sphere
More informationf k j M k j Among them we recognize the gradient operator on polynomial-valued functions which is part of the complete gradient (see [9]) u j xj.
Advances in Applied Clifford Algebras 4, No. 1 (1994) 65 FUNCTIONS OF TWO VECTOR VARIABLES F. Sommen* and N. Van Acker Department of Mathematical Analysis, University of Gent, Galglaan 2 B-9000 Gent, Belgium
More informationTighter Uncertainty Principles Based on Quaternion Fourier Transform
Adv. Appl. Clifford Algebras 6 (016, 479 497 c 015 Springer Basel 0188-7009/010479-19 published online July, 015 DOI 10.1007/s00006-015-0579-0 Advances in Applied Clifford Algebras Tighter Uncertainty
More informationA COUNTEREXAMPLE TO AN ENDPOINT BILINEAR STRICHARTZ INEQUALITY TERENCE TAO. t L x (R R2 ) f L 2 x (R2 )
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 5, pp. 6. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) A COUNTEREXAMPLE
More informationMORE NOTES FOR MATH 823, FALL 2007
MORE NOTES FOR MATH 83, FALL 007 Prop 1.1 Prop 1. Lemma 1.3 1. The Siegel upper half space 1.1. The Siegel upper half space and its Bergman kernel. The Siegel upper half space is the domain { U n+1 z C
More informationAGEOMETRICALGEBRA APPROACH TO SOME PROBLEMS OF ROBOT VISION
AGEOMETRICALGEBRA APPROACH TO SOME PROBLEMS OF ROBOT VISION Gerald Sommer Institut für Informatik und Praktische Mathematik Christian-Albrechts-Universität zu Kiel, Kiel, Germany gs@ks.informatik.uni-kiel.de
More informationON THE BOUNDEDNESS BEHAVIOR OF THE SPECTRAL FACTORIZATION IN THE WIENER ALGEBRA FOR FIR DATA
ON THE BOUNDEDNESS BEHAVIOR OF THE SPECTRAL FACTORIZATION IN THE WIENER ALGEBRA FOR FIR DATA Holger Boche and Volker Pohl Technische Universität Berlin, Heinrich Hertz Chair for Mobile Communications Werner-von-Siemens
More informationIntroduction to Gabor Analysis
Theoretical and Computational Aspects Numerical Harmonic Group under the supervision of Prof. Dr. Hans Georg Feichtinger 30 Oct 2012 Outline 1 2 3 4 5 6 7 DFT/ idft Discrete Given an input signal f of
More informationShannon-Like Wavelet Frames on a Class of Nilpotent Lie Groups
Bridgewater State University From the SelectedWorks of Vignon Oussa Winter April 15, 2013 Shannon-Like Wavelet Frames on a Class of Nilpotent Lie Groups Vignon Oussa Available at: https://works.bepress.com/vignon_oussa/1/
More informationAlgebra C Numerical Linear Algebra Sample Exam Problems
Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric
More information