An Uncertainty Principle for Quaternion Fourier Transform

Transcription

1 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, Japan b Division of Mathematical Sciences, Osaka Kyoiku University, Osaka , Japan Abstract We review the quaternionic Fourier transform (QFT). Using the properties of the QFT we establish an uncertainty principle for the right-sided QFT. This uncertainty principle prescribes a lower bound on the product of the effective widths of quaternion-valued signals in the spatial and frequency domains. It is shown that only a Gaussian quaternion signal minimizes the uncertainty. Key words: Quaternion algebra, Quaternionic Fourier transform, Uncertainty principle, Gaussian quaternion signal, Hypercomplex functions 99 MSC: 3G35, B, 9A, R5 Introduction Recently it has become popular to generalize the Fourier transform (FT) from real and complex numbers [] to quaternion algebra. In these constructions many FT properties still hold, others are modified. The quaternionic Fourier transform (QFT) plays a vital role in the representation of signals. It transforms a real (or quaternionic) two-dimensional signal into a Current address: School of Mathematical Sciences, Universiti Sains Malaysia, 8 Penang, Malaysia. Corresponding author. addresses: mawardibahri@gmail.com (Mawardi Bahri), hitzer@mech.fukui-u.ac.jp (Eckhard S. M. Hitzer), hayashi@soliton.fukui-u.ac.jp (Akihisa Hayashi), ashino@cc.osaka-kyoiku.ac.jp (Ryuichi Ashino). Preprint submitted to Elsevier June 3

2 quaternion-valued frequency domain signal. The four QFT components separate four cases of symmetry in real signals instead of only two in the complex FT [,3]. The QFT was first proposed by Ell []. He proposed a two-sided QFT and demonstrated some important properties of this type of QFT. He also introduced the use of the QFT in the analysis of two dimensional linear time invariant dynamic systems. Later, Bülow [] made a more extended investigation of important properties of the two-sided QFT mainly for real signals and applied it to signal and image processing. He obtained a local D quaternionic phase. Pei et al. [5] discussed and optimized the implementation of different types of the QFT with applications to linear quaternion filters. Hitzer [6] described in detail properties of different types of QFT applied to fully quaternionic signals and then generalized the QFT to a volume time, as well as to a space time algebra Fourier transform. Sangwine and Ell [7] proposed the QFT application to color image analysis. Bas, Le Bihan and Chassery [8] used the QFT to design a digital color image watermarking scheme. Bayro et al. [9] applied the QFT in image pre-processing and neural computing techniques for speech recognition. It is well known that the uncertainty principle for the FT relates the variances of a function and its Fourier transform which cannot both be simultaneously sharply localized [,]. In signal processing an uncertainty principle states that the product of the variances of the signal in the time and frequency domains has a lower bound. Yet Felsberg [3] notes for two dimensions: In D however, the uncertainty relation is still an open problem. In [] it is stated that there is no straightforward formulation for the D uncertainty relation. A first straightforward directional D uncertainty principle was formulated by Hitzer and Mawardi [3], in Clifford algebras Cl n, with n = (mod ). Now we attempt another formulation for quaternions H = Cl, using the right-sided QFT. This paper briefly reviews the QFT and provides alternative proofs for some of its properties. The QFT considered in this paper enables us to extend the Heisenberg type uncertainty principle from the complex FT to the QFT. The organization of the present paper is as follows. In section, we briefly establish our notation for quaternion algebra and its relationship with the Clifford geometric algebra Cl 3,. In section 3, we demonstrate some important properties of the QFT, which are necessary to prove the uncertainty principle for the QFT. In section, the classical Heisenberg uncertainty principle is generalized for the QFT.

3 Quaternion Algebra. The Quaternion Algebra H The quaternion algebra [] was first invented by Sir W. R. Hamilton in 83 and is denoted by H in his honor. It is an extension of complex numbers to a fourdimensional algebra. Every element of H is a linear combination of a real scalar and three orthogonal imaginary units (denoted i, j, and k) with real coefficients H = {q = q + iq + jq + kq 3 q, q, q, q 3 R}, (.) where the elements i, j, and k obey Hamilton s multiplication rules i = j = k =, ij = ji = k, jk = kj = i, ki = ik = j. (.) Because H is according to (.) non-commutative, one cannot directly extend various results on complex numbers to quaternions. For simplicity, we express a quaternion q as sum of a scalar q, and a pure 3D quaternion q q = q + q = q + iq + jq + kq 3, (.3) where the scalar part is also denoted Sc(q) = q. The conjugate of a quaternion q is obtained by changing the sign of the pure part, i. e. q = q q = q q i q j q 3 k. (.) The quaternion conjugation (.) is a linear anti-involution p = p, p + q = p + q, pq = q p, p, q H. (.5) Given a quaternion q and its conjugate, we can easily check that the following properties are correct q = (q + q), q = (q q), q = q q = q. (.6) Using (.) the multiplication of the two quaternions q = q + q and p = p + p can be expressed as qp = q p + q p + q p + p q + q p, (.7) where we recognize the scalar product q p = (q p + q p + q 3 p 3 ) and the antisymmetric cross type product q p = i(q p 3 q 3 p ) + j(q 3 p q p 3 ) + k(q p 3

4 q p ). The scalar part of the product is Sc(qp) = q p + q p, (.8) and the pure part is q p + p q + q p. (.9) Especially, if both q and p are pure quaternions (.7) reduces to qp = q p + q p. (.) According to (.7) the multiplication of a quaternion q and its conjugate can be expressed as q q = q q q q + q ( q) + q q + q ( q) = q o + q + q + q 3. (.) Equation (.) leads to the modulus q of a quaternion q defined as q = q q = q o + q + q + q 3. (.) It is straightforward to see that with (.5) and (.) the following modulus properties hold pq = p q, p = p, p, q H. (.3) Using the conjugate (.) and the modulus of a quaternion q, we can define the inverse of q H \ {} as q = q (.) q which shows that H is a normed division algebra. For unit quaternions with q = equation (.) simplifies to and for pure unit quaternions equation (.) becomes q = q, (.5) q = q. (.6) It is important to note that with (.6), we have for two quaternion-valued functions f, g (independent of their domain space) (gf + fg) = g f g f = Sc(gf). (.7)

5 . Quaternion Module According to (.) a quaternion-valued function f : R H can be expressed as f(x) = f + if (x) + jf (x) + kf 3 (x), f, f, f, f 3 R. (.8) We introduce an inner product of functions f, g defined on R with values in H as follows f, g = f(x)g(x) d x, (.9) R and its associated scalar norm f by defining f = f, f = f(x) d x. (.) R The quaternion module L (R ; H) is then defined as L (R ; H) = {f f : R H, f < }. (.).3 Connection Between H and Clifford Algebras Cl 3, and Cl, Quaternions are isomorphic to the even subalgebra Cl + 3, of scalars and bivectors (see [5,6]) of the real associative 8-dimensional Clifford geometric algebra Cl 3,. The latter has the basis of scalar, e, e, e 3 vectors, e e, e 3 e, e e 3 i 3 = e e e 3 bivectors, pseudoscalar. (.) In equation (.) the set {e, e, e 3 } is an orthonormal vector basis of the real 3D Euclidean vector space R 3. The isomorphism means that any quaternion q can be expanded in the form q Cl + 3, = Cl 3, + Cl 3,, i.e. q = α + i 3 b, (.3) where α R and vector b R 3. Equation (.3) tells us that the elements of Cl + 3, form a four-dimensional linear space with one scalar and three bivector dimensions. Quaternions are also isomorphic to Cl,. For this we identify i, j with vectors e, e with square, respectively, and k as their product e e. This fact is helpful for defining quaternionic Fourier and wavelet transforms and to compare them with other Clifford Fourier and wavelet transformations [6,7]. 5

6 3 Quaternionic Fourier Transform (QFT) It is natural to extend the Fourier transform to quaternion algebra. These extensions are broadly called the quaternionic Fourier transform (QFT). Due to the noncommutative properties of quaternions, there are three different types of QFT: a left-sided QFT, a right-sided QFT and a two-sided QFT [5]. By reasons explained in more detail below we choose to apply the right-sided QFT of a D quaternionvalued signal. This version of the QFT defined here is also known as the D Clifford FT of Delanghe, Sommen and Brackx [5]. 3. Definition of QFT Definition 3.. The QFT of f L (R ; H) is the function F q {f}: R H given by F q {f}(ω) = f(x)e iω x e jω x d x, (3.) R where x = x e + x e, ω = ω e + ω e, and the quaternion exponential product e iω x e jω x is the quaternion Fourier kernel. Remark 3.. Apart from the convention used in Definition 3. with in the inverse QFT (3.5), there are two other common conventions: One is obtained by sub- (π) stituting (3.) ω πω. The other is obtained by evenly distributing the π factors between the transformation and the inverse transformation F q = π... d x, Fq = π... d ω. All calculations in this paper can easily be converted to these other conventions. Using the Euler formula for the quaternion Fourier kernel we can rewrite (3.) in the following form F q {f}(ω) = f(x) cos(ω x ) cos(ω x ) d x R f(x) i sin(ω x ) cos(ω x ) d x R f(x) j cos(ω x ) sin(ω x ) d x R + f(x) k sin(ω x ) sin(ω x ) d x. (3.) R Equation (3.) clearly shows how the QFT separates real signals into four quaternion components, i. e. the even-even, odd-even, even-odd and odd-odd components of f. Let us now take an example to illustrate this expression. Example 3.3. Consider the quaternionic distribution signal (see Fig. ), i. e. the 6

7 3 3 x x x 3 3 x 3 3 x x x 3 3 x Fig.. Quaternionic signal of Example 3.3 in the spatial domain (u = v = ). The resulting patterns are identical, apart from π/ phase shifts along x and x [8]. QFT kernel of (3.) f(x) = e jv x e iu x. (3.3) It is easy to see that the QFT of f is a Dirac quaternion function, i. e. F q {f}(ω) = (π) δ(ω ω ), ω = u e + v e. (3.) The following theorem tells us that the QFT is invertible, that is, the original signal f can be recovered by simply taking the inverse of the quaternionic Fourier transform (3.). Theorem 3.. Suppose that f L (R ; H) and F q {f} L (R ; H). Then the QFT F q {f} of f is an invertible transform and its inverse is given by Fq [F q {f}](x) = f(x) = F (π) q {f}(ω)e jω x e iω x d ω. (3.5) R 7

8 3. Major Properties of the QFT This subsection describes important properties of the QFT which will be used to establish a new uncertainty principle for the QFT. For detailed discussions of the properties of the QFT and their proofs, see e.g. [6,5,8]. We now first establish a Plancherel theorem, specific to the right-sided QFT. Theorem 3.5 (QFT Plancherel). The inner product (.9) of two quaternion module functions f, g L (R ; H) and their QFT is related by f, g L (R ;H) = (π) F q{f}, F q {g} L (R ;H). (3.6) In particular, with f = g, we get Parseval s theorem, i. e. f L (R ;H) = (π) F q{f} L (R ;H). (3.7) This shows that the total signal energy computed in the spatial domain is equal to the total signal energy computed in the quaternionic domain. The Parseval theorem allows the energy of a quaternion-valued signal to be considered in either the spatial domain or the quaternionic domain and the change of domains for convenience of computation. In following we give an alternative proof of Plancherel s theorem (compare to Hitzer [6]). = (π) R R f, g L (R ;H) = [ F (π) q {f}(ω)e jω x e iω x d ω R R ] e iω x e jω x F q {g}(ω )d ω d x R = (π) F q {f}(ω)e jω x e ix (ω ω ) e jω x F q {g}(ω )d xd ω d ω R F q {f}(ω)δ(ω ω )F q {g}(ω )d ωd ω R R = F (π) q {f}(ω)f q {g}(ω)d ω. (3.8) R This completes the proof of theorem 3.5. Due to the non-commutativity of the quaternion exponential product factors we only have a left linearity property for general linear combinations with quaternionic constants and a special shift property. Different from the QFT Plancherel Theorem 3.5, the Parseval theorem (3.7) can be established for all three variants of the QFT: right-sided, left-sided and two-sided. 8

9 Theorem 3.6 (Left linearity property). The QFT of two quaternion module functions f, g L (R ; H) is a left linear operator, i. e. F q {µf + λg}(ω) = µf q {f}(ω) + λf q {g}(ω), (3.9) where µ and λ H are quaternionic constants. Theorem 3.7 (Shift property). If the argument of f L (R ; H) is offset by a constant vector x = x e + y e, i. e. fx (x) = f(x x ), then [6] Proof. Equation (3.) gives F q {fx }(ω) = F q {fe iω x }(ω) e jω y. (3.) F q {fx }(ω) = f(x x ) e iω x e jω x d x. R (3.) We substitute t for x x in the above expression, and get with d x= d t F q {fx }(ω) = f(t) e iω (t +x ) e jω (t +y ) d t R ) = (f(t)e iω x e iω t e jω t d t e jω y. (3.) This proves (3.). R Dual to Theorem 3.7 the following modulation type formula holds for the inverse QFT. Theorem 3.8. If the argument of F q {f} L (R ; H), F q {f} L (R ; H) is offset by a constant frequency vector ω = ω e + ω e R, then f (x) and F q {f}(ω) are related by f (x) = Fq [F q {f}(ω ω )](x) = Fq {F q {f}(ω) e jω x }(x) e iω x. (3.3) Remark 3.9. Equation (3.) and (3.3) are specific for the right-sided definition of Definition 3.. The usual form of the modulation property of the complex FT does not hold for the QFT. It is obstructed by the non-commutativity of the exponential factors (eq. (38) of [6]) e iω x e jω x e jω x e iω x. (3.) Next we give an explicit proof of the derivative properties stated in Table of [6]. The QFT is also right linear for real constants µ, λ R. 9

10 Theorem 3.. If the QFT of the n-th partial derivative of f L (R ; H) with respect to the variable x exists and is in L (R ; H), then the QFT of n f i n is x n given by F q { n f i n }(ω) = ω x n n F q {f}(ω), n N. (3.5) Proof. We first prove the theorem for n =. Applying integration by parts and using the fact that f tends to zero for x we immediately obtain F q { ( ) f i }(ω) = f(x) i e iω x e jω x d x x R x [ ( ) ] = f(x) i e iω x dx e jω x dx R R x [ = f(x) i e iω x x = x = R f(x) i ] e iω x dx e jω x dx R x = f(x) ω e iω x e jω x d x R = ω F q {f}(ω). (3.6) Using mathematical induction we can finish the proof of Theorem 3.. Theorem 3.. If the QFT of the m-th partial derivative of a quaternion-valued function f L (R ; H) with respect to the variable x exists and is in L (R ; H), then F q { m f } (ω) = F x m q {f}(ω)(jω ) m, m N. (3.7) Proof. Direct calculation gives f(x) x = x (π) = (π) = (π) F q {f}(ω)e jω x e iω x d ω R ( ) F q {f}(ω) e jω x e iω x d ω R x [F q {f}(ω) jω ] e jω x e iω x d ω R = F q [F q {f}(ω)jω ]. (3.8) We therefore get F q { f x } (ω) = F q {f}(ω)jω. (3.9)

11 By successive differentiation with respect to the variable x and with induction we easily obtain This ends the proof of (3.7). F q { m f } (ω) = F x m q {f}(ω)(jω ) m, m N. (3.) As consequence of Theorem 3. we immediately obtain the following corollary. Corollary 3.. Suppose that the QFT of a partial derivative n+m f/ x n x m of a quaternion-valued function f L (R ; H) is in L (R ; H), and that f = f + if, then F q { n+m f x n x m }(ω) = (iω ) n F q {f}(ω)(jω ) m, m, n N. (3.) Proof. For n N and m = multiplication of (3.5) with i n from the left gives F q {( x ) n f}(ω) = (iω ) n F q {f}(ω), n N. (3.) The combination of (3.) and (3.) for f = f + if gives (3.). Another important consequence of Theorems 3. and 3. is formulated in the following lemma. Lemma 3.3. If the QFT of the st partial derivative of f L (R ; H) with respect to the variable x k, k {, } exists and is in L (R ; H), then (π) R f(x) d x = ω x k F q {f}(ω) d ω, k {, }. (3.3) k R Proof. For k = straightforward calculation using Parseval s theorem (3.7) and Theorem 3. gives (π) R f(x) d x (.3) = (π) f(x) i x R d x x (3.7) = F q { f i }(ω) d ω R x (3.5) = For k = we similarly use Theorem 3. to get R ω F q {f}(ω) d ω. (3.)

12 Table Properties of quaternionic Fourier transform of quaternion functions f, g L (R ; H), the constants are α, β H, a R \ {}, x = x e + y e R and n N. Property Quaternion Function Quaternionic Fourier Transform Left linearity αf(x)+βg(x) αf q {f}(ω)+ βf q {g}(ω) Scaling f(ax) a F q{f}( ω a ) x-shift f(x x ) F q {fe iω x }(ω) e jω y Part. deriv. ( ) n x f(x) i n ω nf q{f}(ω), f L (R ; H) ( ) n x f(x) (iω ) n F q {f}(ω), f = f + if ( ) n x f(x) Fq {f}(ω)(jω ) n, f L (R ; H) Plancherel f, f L (R ;H) = (π) F q {f }, F q {f } L (R ;H) Parseval f L (R ;H) = π F q{f} L (R ;H) (π) R f(x) d x (3.7) = F q { f }(ω) d ω x R x (3.7) = (.3) = F q {f}(ω)jω d ω. R R ω F q {f}(ω) d ω. (3.5) Some important properties of the QFT are summarized in Table. For more details we refer to [6]. Example 3.. Consider a two-dimensional Gaussian quaternion function (Fig. ) of the form f(x) = C e (a x +a x ), (3.6) where C = C +ic +jc +kc 3 H is a quaternion constant and a, a R are positive real constants. Then the QFT of f as shown Fig. 3 is given by F q {f}(ω) = C e a x e a x e iω x e jω x dx dx R R ( ) = C e a x e a x e iω x dx e jω x dx R R ( ) = C e a x π e ω /(a ) e jω x dx. R a π = C e ω /(a ) π e ω /(a ) a a π = C e ( ω a a a + ω a ). (3.7)

13 x x x x x x x x Fig.. Quaternion Gaussian function for a = a =, C =, C =, C =, and C 3 = 5 in the spatial domain. Top row: real part and imaginary part i. Bottom row: imaginary parts j and k. This shows that the QFT of the Gaussian quaternion function is another Gaussian quaternion function. Uncertainty principle for QFT In physics the uncertainty principle [] was introduced for the first time 8 years ago by Heisenberg who demonstrated the impossibility of simultaneous precise measurements of a particle s momentum and its position. In a communication theory setting, an uncertainty principle states that a signal cannot be arbitrarily confined in both the spatial and frequency domains. Many efforts have been devoted to extend the uncertainty principle to various types of functions and Fourier transforms. Shinde et al. [9] established an uncertainty principle for fractional Fourier transforms which provides a lower bound on the uncertainty product of signal representations in both time and frequency domains for real signals. Korn [] proposed Heisenberg type uncertainty principles for Cohen transforms which describe lower limits for the time-frequency concentration. In our previous papers [3,6,7,], we established a new directional uncertainty principle for the Clifford Fourier transform which describes how the variances (in arbitrary but fixed 3

14 ω ω ω ω ω ω ω ω Fig. 3. Quaternion Gaussian function in the quaternionic frequency domain. Top row: real part and imaginary part i. Bottom row: imaginary parts j and k. directions) of a multivector-valued function and its Clifford Fourier transform are related. Bülow [] showed a quaternion uncertainty principle, for quaternion valued signals according to which x x ω ω 6π, (.) where x k, k =, is the effective width and ω k, k =, is its effective bandwidth, defined as in Definition., only replacing F q by a two-sided version of the QFT. [The factor results from Bülow s use of the linear substitution π ω πω in (3.), compare Remark 3..] He showed that a Gabor filter can lead to equality in (.). It must be remembered that he applied the two-sided QFT for his uncertainty principle. His uncertainty principle is similar to the uncertainty principle for the conventional two-dimensional Fourier transform. In the following we explicitly generalize and prove the classical uncertainty principle to quaternion module functions. We also give an explicit proof for Gaussian quaternion functions (Gabor filters) to be indeed the only functions that minimize the uncertainty. We further emphasize that our generalization is non-trivial because the multiplication of quaternions and the quaternion Fourier kernel are non-commutative. For this purpose we introduce the following definition.

15 Definition.. Let f L (R ; H) be a quaternion-valued signal such that x k f L (R ; H), k =,, and let F q {f} L (R ; H) be its QFT such that ω k F q {f} L (R ; H), k =,. The effective spatial width or spatial uncertainty x k of f is evaluated by x k = V ar k {f}, k {, }, (.) where V ar k {f} is the variance of the energy distribution of f along the x k axis defined by V ar k {f} = x kf L (R ;H) f L (R ;H) = R f(x) x k d x, k {, }. (.3) R f(x) d x Similarly, in the quaternionic domain we define the effective spectral width as ω k = V ar k {F q {f}}, k {, }, (.) where V ar k {F q {f}} is the variance of the frequency spectrum of f along the ω k frequency axis given by V ar k {F q {f}} = ω kf q {f} L (R ;H) F q {f} L (R ;H) = R F q{f}(ω) ω k d ω R F q{f}(ω) d ω. (.5) Theorem.. Let f L (R ; H) be a quaternion-valued signal such that both ( + x k )f(x) L (R ; H) and f(x) L (R ; H) for k =,. Then two uncertainty relations are fulfilled x ω, and x ω. (.6) The combination of the two spatial uncertainty principles above leads to the uncertainty principle for the two-dimensional quaternion signal f(x) of the form x x ω ω. (.7) Equality holds in (.7) if and only if f is a Gaussian quaternion function, i. e. f(x) = K e ( x + x a a ), (.8) where K is a quaternion constant, and a, a R are positive real constants. Analogous to complex numbers, we will use equation (.7) to derive the following lemma which will be necessary to prove Theorem.. Lemma.3. For two quaternion-valued functions f, g L (R ; H), the Schwarz inequality takes the form [ (g f ] + fḡ)d x f fd x gḡd x = f d x g d x (.9) R R R R R 5

16 Remark.. An alternative form of Lemma.3 is given by equation (.). To prove Schwarz s inequality, let ɛ R be a real constant. Then (f + ɛg)(f + ɛg) d x. R (.) Applying the second part of (.5) and then expanding the above inequality we can rewrite it in the following form (f + ɛg)( f + ɛḡ) d x R = f f d x + ɛ (g f + fḡ) d x + ɛ R R R gḡ d x. (.) The right-hand side of equation (.) is a quadratic expression in ɛ. The discriminant of this quadratic polynomial must be negative or zero and gives therefore [ (g f ] + fḡ)d x f fd x gḡd x, (.) R R R which is equivalent to (.9). This finishes the proof of Lemma.3. Now let us begin the proof of Theorem.. Proof. We prove Theorem. for k,. First, by applying Lemma 3.3 and the Parseval theorem (3.7) we immediately obtain x k ω k = Lemma 3.3 = (3.7) = Lemma.3 = R x k f(x) d x R ω k F q {f}(ω) d ω R f(x) d x R F q{f}(ω) d ω R x k f(x) d x (π) R f(x) d x R f(x) d x R F q{f}(ω) d ω R x kf(x) d x R f(x) d x ( R f(x) d x) [ R ( f(x) x k f(x) + xk f(x) f(x) ) d x ] f L (R ;H) ( R x [ ] k f(x) f(x) d x ). (.3) f L (R ;H) Second, using integration by parts we further get 6

17 x k ω k ( [ R x k f(x) dx l ] x k= x k = R f(x) d x ) = ( R f(x) d x) f L (R ;H) f L (R ;H) where l {, }, l k. This proves (.6). =, Remark.5. Consequently replacing the right-sided QFT (3.) by the left-sided QFT Fq left {f}(ω) = R e jω x e iω x f(x) d x, allows to establish a corresponding Parseval theorem, left-sided QFT formulas for the partial derivatives (analogous to Theorems 3. and 3.), and formulas for the norms k f, k {, } (corresponding to Lemma 3.3). Theorem. applies therefore also to the left-sided QFT. We finally show that the equality in (.6) is satisfied if and only if f is a Gaussian quaternion function. Using (.7) we can rewrite Lemma.3 in the following form (compare to Chui []) [ ] Sc(gh)d x h d x g d x. (.) R R R For k =, we first take g = x k f(x) L (R ; H) and h = Equation (.) can then be expressed as [ ( Sc x k f(x) ) f(x) d x] x k f(x) d x R R Equality in (.5) implies that ( Sc x k f(x) ) f(x) f(x) L (R ; H). R f(x) d x. (.5) = x k f(x) f(x), (.6) and x k f(x) = a k f(x), (.7) where the a k are positive real constants. From equation (.7) we obtain x k f(x) = q a k f(x), k =,, (.8) where q is a unit quaternion. Equation (.6) implies that x k f(x) f(x). (.9) 7

18 Multiplying both sides of (.8) by f(x) we get x k f(x) f(x) = q a k f(x), k =,. (.) Applying (.9) we get q =. Hence, we conclude that f(x) = a k x k f(x), k =,. (.) Solving the equations (.) we further obtain that f must be a Gaussian quaternion function f(x) = K e ( x + x a a ), k =,, (.) where K H is a quaternion constant. Since the Gaussian quaternion function f(x) of (.) achieves the minimum widthbandwidth product, it is theoretically a very good prototype wave form. One can therefore construct a basic wave form using spatially or frequency scaled versions of f(x) to provide multiscale spectral resolution. Such a wavelet basis construction from a Gaussian quaternion function prototype waveform has for example been realized in the quaternion wavelet transforms of []. The optimal localization in space and frequency is also the reason why (algebraically related to quaternions) two-dimensional Clifford Gabor bandpass filters (with Gaussian impulse response) were suggested in [3]. 5 Conclusion Using the basic concepts of quaternion algebra H we introduced the two-dimensional quaternionic Fourier transform (QFT). Important properties of the QFT such as partial derivative, Plancherel and Parseval theorems, specific shift- and modulation properties, and the quaternion function Schwarz inequality were demonstrated. We finally proposed a new uncertainty principle for the right-sided QFT. So far no such uncertainty principle for a one-sided QFT (left or right) had been established. In our previous works on Clifford FT uncertainty principles, we always had the benefit of invariant vector derivatives. As far as we know, in quaternion calculus no suitable analogue to such a vector derivative has been established. Before introducing the QFT Plancherel theorem 3.5, we pointed out that this theorem is specific for the right-sided QFT. With the definition of the inner product (.9) it seems not possible to establish a similar QFT Plancherel theorem for the left-sided QFT. But as explained in Remark.5, the uncertainty principle for the right-sided QFT can be shown to apply to the left-sided QFT as well. 8

19 The case of the two-sided QFT is solved. A Parseval theorem can be shown [,6] and equation (3.) holds for general f L (R ; H). Acknowledgements We do thank O. Yasukura for helpful comments and T. Khairuman for assistance in producing the figures. References [] R. Bracewell, The Fourier Transform and its Applications, third ed., McGraw-Hill Book Company, New York,. [] T. Bülow, Hypercomplex Spectral Signal Representations for the Processing and Analysis of Images, Ph.D. Thesis, Institut für Informatik und Praktische Mathematik, University of Kiel, Germany, 999. [3] M. Felsberg, Low-Level Image Processing with the Structure Multivector, Ph.D. Thesis, Institut für Informatik und Praktische Mathematik, University of Kiel, Germany,. [] T. A. Ell, Quaternion-Fourier Transforms for Analysis of Two-dimensional Linear Time-Invariant Partial Differential Systems, in: Proceeding of the 3nd Conference on Decision and Control, San Antonio, Texas, 993, pp [5] S. C. Pei, J. J. Ding and J. H. Chang, Efficient Implementation of Quaternion Fourier Transform, Convolution, and Correlation by -D Complex FFT, IEEE Transactions on Signal Processing 9() () [6] E. Hitzer, Quaternion Fourier Transform on Quaternion Fields and Generalizations, Advances in Applied Clifford Algebras 7(3) (7) [7] S. J. Sangwine and T. A. Ell, Hypercomplex Fourier Transforms of Color Images, IEEE Transactions on Image Processing 6() (7) 35. [8] P. Bas, N. Le Bihan and J.M. Chassery, Color Image Watermarking using Quaternion Fourier Transform, in: Proceedings of the IEEE International Conference on Acoustics Speech and Signal Processing (ICASSP), Hong-Kong, 3, pp [9] E. Bayro-Corrochano, N. Trujillo and M. Naranjo, Quaternion Fourier Descriptors for Preprocessing and Recognition of Spoken Words Using Images of Spatiotemporal Representations, Journal of Mathematical Imaging and Vision 8() (7) [] H. Weyl, The Theory of Groups and Quantum Mechanics, second ed., Dover, New York, 95. 9

20 [] C. K. Chui, An Introduction to Wavelets, Academic Press, New York, 99. [] G. H. Granlund and H. Knutsson, Signal Processing for Computer Vision, Kluwer, Dordrecht, 995. [3] E. Hitzer and B. Mawardi, Clifford Fourier Transform on Multivector Fields and Uncertainty Principle for Dimensions Cl n,, n = (mod ) and n = 3 (mod ), accepted for P. Angeles (ed.), Proceedings of the Seventh International Conference on Clifford Algebra (ICCA7), Toulouse, France, May 9-9, 5. [] J. B. Kuipers, Quaternions and Rotation Sequences, Princeton University Press, New Jersey, 999. [5] F. Brackx, R. Delanghe and F. Sommen, Clifford Analysis, Vol. 76 of Research Notes in Mathematics, Pitman Advanced Publishing Program, Boston, 98. [6] B. Mawardi and E. Hitzer, Clifford Fourier Transformation and Uncertainty Principle for the Clifford Geometric Algebra Cl 3,, Advances in Applied Clifford Algebras, 6() (6) 6. [7] B. Mawardi and E. Hitzer, Clifford Algebra Cl 3, -valued Wavelet Transformation, Clifford Wavelet Uncertainty Inequality and Clifford Gabor Wavelets, International Journal of Wavelet, Multiresolution and Information Processing 5(6) (7) [8] T. Bülow, M. Felsberg and G. Sommer, Non-commutative Hypercomplex Fourier Transforms of Multidimensional Signals, in: G. Sommer (ed.), Geometric Computing with Clifford Algebras, (Chapter 8). Springer, Heidelberg,, pp [9] S. Shinde and V. M. Gadre, An Uncertainty Principle for Real Signals in the Fractional Fourier Transform Domain, IEEE Transactions on Signal Processing 9() () [] P. Korn, Some Uncertainty Principle for Time-Frequency Transforms for the Cohen Class, IEEE Transactions on Signal Processing 53() (5) [] E. Hitzer and B. Mawardi, Uncertainty Principle for the Clifford Geometric Algebra Cl n,, n = 3(mod ) Based on Clifford Fourier transform, in: T. Qian, M. I. Vai, and Y. Xu (eds.), the Springer (SCI) book series Applied and Numerical Harmonic Analysis, 6, pp [] E. Bayro-Corrochano, The Theory and Use of the Quaternion Wavelet Transform, Journal of Mathematical Imaging and Vision () (6) [3] F. Brackx, N. De Schepper and F. Sommmen, The Two-dimensional Clifford-Fourier Transform, Journal of mathematical Imaging and Vision 6() (6) 5 8.