AAA Abstract and Applied Analysis 1687-0409 1085-3375 Hindawi Publishing Corporation 162769 10.1155/2013/162769 162769 Research Article Convolution Theorems for Quaternion Fourier Transform: Properties and Applications Bahri Mawardi 1 http://orcid.org/0000-0003-4807-4300 Ashino Ryuichi 2 Vaillancourt Rémi 3 Apreutesei Narcisa C. 1 Department of Mathematics Hasanuddin University Makassar 90245 Indonesia unhas.ac.id 2 Division of Mathematical Sciences Osaka Kyoiku University Osaka 582-8582 Japan osaka-kyoiku.ac.jp 3 Department of Mathematics and Statistics University of Ottawa Ottawa ON Canada K1N 6N5 uottawa.ca 2013 3 11 2013 2013 01 06 2013 01 09 2013 07 09 2013 2013 Copyright © 2013 Mawardi Bahri et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

General convolution theorems for two-dimensional quaternion Fourier transforms (QFTs) are presented. It is shown that these theorems are valid not only for real-valued functions but also for quaternion-valued functions. We describe some useful properties of generalized convolutions and compare them with the convolution theorems of the classical Fourier transform. We finally apply the obtained results to study hypoellipticity and to solve the heat equation in quaternion algebra framework.

1. Introduction

Convolution is a mathematical operation with several applications in pure and applied mathematics such as numerical analysis, numerical linear algebra, and the design and implementation of finite impulse response filters in signal processing. In , the authors introduced the Clifford convolution. It is found that some properties of convolution, when generalized to the Clifford Fourier transform (CFT), are very similar to the classical ones.

On the other hand, the quaternion Fourier transform (QFT) is a nontrivial generalization of the classical Fourier transform (FT) using quaternion algebra. The QFT has been shown to be related to the other quaternion signal analysis tools such as quaternion wavelet transform, fractional quaternion Fourier transform, quaternionic windowed Fourier transform, and quaternion Wigner transform . A number of already known and useful properties of this extended transform are generalizations of the corresponding properties of the FT with some modifications, but the generalization of convolution theorems of the QFT is still an open problem. In the recent past, several authors  tried to formulate convolution theorems for the QFT. But they only treated them for real-valued functions which is quite similar to the classical case. In , the authors briefly introduced, without proof, the QFT of the convolution of two-dimensional quaternion signals.

In this paper, we establish general convolutions for QFT. Because quaternion multiplication is not commutative, we find new properties of the QFT of convolution of two quaternion-valued functions. These properties describe closely the relationship between the quaternion convolution and its QFT. The generalization of the convolution theorems of the QFT is mainly motivated by the Clifford convolution of general geometric Fourier transform, which has been recently studied in [15, 16]. We further establish the inverse QFT of the product of the QFT, which is very useful in solving partial differential equations in quaternion algebra framework.

This paper consists of the following sections. Section 2 deals with some results on the real quaternion algebra and the definition of the QFT and its basic properties. We also review some basic properties of QFT, which will be necessary in the next section. Section 3 establishes convolution theorems of QFT and some of their consequences. Section 4 presents an application of QFT to study hypoellipticity and to solve the heat equation in quaternion algebra. Some conclusions are drawn in Section 5.

2. Quaternion Algebra

For convenience, we specify the notation used in this paper. The quaternion algebra over , denoted by , is an associative noncommutative four-dimensional algebra, (1)={q=q0+iq1+jq2+kq3;q0,q1,q2,q3}, which obeys the following multiplication rules: (2)ij=-ji=k,jk=-kj=i,ki=-ik=j,i2=j2=k2=ijk=-1. For a quaternion q=q0+iq1+jq2+kq3, q0 is called the scalar part of q denoted by Sc(q) and iq1+jq2+kq3 is called the vector (or pure) part of q. The vector part of q is conventionally denoted by q. Let p, q and let p, q be their vector parts, respectively. It is common to write for short (3)q·p=q1p1+q2p2+q3p3,q×p=i(q2p3-q3p2)+j(q3p1-q1p3)+k(q1p2-q2p1). Then, (2) yields the quaternionic multiplication qp as (4)qp=q0p0-q·p+q0p+p0q+q×p. The quaternion conjugate of q, given by (5)q-=q0-iq1-jq2-kq3,q0,q1,q2,q3, is an anti-involution; that is, (6)qp¯=p-q-. From (5) we obtain the norm or modulus of q defined as (7)|q|=qq-=q02+q12+q22+q32. It is not difficult to see that (8)|qp|=|q||p|,p,q. Using the conjugate (5) and the modulus of q, we can define the inverse of q{0} as (9)q-1=q-|q|2, which shows that is a normed division algebra. As in the algebra of complex numbers, we can define three nontrivial quaternion involutions : (10)α(q)=-iqi=-i(q0+iq1+jq2+kq3)i=q0+iq1-jq2-kq3,β(q)=-jqj=-j(q0+iq1+jq2+kq3)j=q0-iq1+jq2-kq3,γ(q)=-kqk=-k(q0+iq1+jq2+kq3)k=q0-iq1-jq2+kq3.

Hereinafter, besides the quaternion units i, j, and k and the vector part q of a quaternion q, we will use the real vector notation: (11)x=(x1,x2)2,|x|2=x12+x22,x·y=x1y1+x2y2,f(x)=f(x1,x2), and so on when there is no confusion. This gives the following definition.

Definition 1 (see [<xref ref-type="bibr" rid="B10">10</xref>]).

A function f:2 is called quaternionic Hermitian if, for the involutions α and β, (12)f(-x1,x2)=β(f(x)),f(x1,-x2)=α(f(x)), for each x2.

For any unit quaternion (13)q=q0+q=cos(θ2)+sin(θ2), and for any vector v3 the action of the operator (14)Lq(v)=qvq- on v is equivalent to a rotation of the vector v through an angle θ about u as the axis of rotation.

It is convenient to introduce an inner product for two functions f,g:2 as follows: (15)(f,g)L2(2;)=2f(x)g(x)¯d2x. In particular, for f=g, we obtain the scalar product of the above inner product (15) given by (16)f,f=fL2(2;)=(2|f(x)|2d2x)1/2.

2.1. Multiindices and Derivatives

A couple α=(α1,α2) of nonnegative integers is called a multiindex. We denote (17)|α|=α1+α2,α!=α1!α2!, and for x2, (18)xα=x1α1x2α2. Derivatives are conveniently expressed by multiindices: (19)α=|α|x1α1x2α2. Denote by {e1,e2} the standard basis of 2. The vector differential   a· along the direction a is defined by (20)a·=a11+a22, where =e11+e22.

2.2. QFT and Its Properties Definition 2.

The QFT of fL2(2;) is the transform q{f}L2(2,) given by the integral (21)q{f}(ω)=2e-iω1x1f(x)e-jω2x2dx, where q is called the quaternion Fourier transform operator or the quaternion Fourier transformation.

Using the Euler formula for the quaternion Fourier kernel e-iω1x1e-jω2x2, we can rewrite (21) in the following form: (22)q{f}(ω)=2f(x)cos(ω1x1)cos(ω2x2)dx-2if(x)sin(ω1x1)cos(ω2x2)dx-2f(x)jcos(ω1x1)sin(ω2x2)dx+2if(x)jsin(ω1x1)sin(ω2x2)dx.

Definition 3.

The inverse QFT of gL2(2;) is the transform q-1{g}L2(2,) given by the integral (23)q-1[g(ω)](x)=1(2π)22eiω1x1g(ω)ejω2x2dω.

Some important properties of the QFT are stated in the following lemmas proved in [12, 17].

Lemma 4.

Let fL1(2;)L2(2;). If q{αf}L1(2;), then (24)q{αf}(ω)=(iω1)α1q{f}(ω)(jω2)α2.

In particular, if q{(2,0)f}L1(2;), then (25)q{(2,0)f}(ω)=(iω1)2q{f}(ω). And if q{(0,2)f}L1(2;), then (26)q{(0,2)f}(ω)=q{f}(ω)(jω2)2.

Lemma 5 (scalar QFT Parseval).

The scalar product of f,gL2(2;) and its QFT are related by (27)f,gL2(2;)=1(2π)2q{f},q{g}L2(2;). And in particular, with f=g, the Plancherel theorem indicates that (28)fL2(2;)2=1(2π)2q{f}L2(2;)2.

This shows that the total signal energy computed in the spatial domain is equal to the total signal energy computed in the quaternion domain.

3. Convolution of QFT

In this section, we establish the quaternion convolution of the QFT which extends the classical convolution to quaternion fields. Let us first define the convolution of two quaternion-valued functions.

Definition 6.

The convolution of fL2(2;) and gL2(2;), denoted by fg, is defined by (29)(fg)(x)=2f(y)g(x-y)dy.

Example 7.

To illustrate the general noncommutativity (fg)(gf), let us compute the convolution of f(x)=ke-|x|2 and g(x)=(i+j)  eiω1x1ejω2x2. Although gL(2;), we can still define the convolution of f and g, because f decays rapidly at infinity. A simple calculation gives (30)(fg)(x)=2ke-|y|2(i+j)eiω1(x1-y1)ejω2(x2-y2)dy=2k(i+j)eiω1x1e-iω1y1e-|y|2e-jω2y2dyejω2x2=k(i+j)eiω1x1e-|ω|2/4ejω2x2=(j-i)eiω1x1e-|ω|2/4ejω2x2. On the other hand, we have (31)(gf)(x)=2(i+j)eiω1y1ejω2y2ke-|x-y|2dy=2(i+j)ke-iω1y1e-|x-y|2e-jω2y2dy=2(i+j)ke-iω1y1e-|y-x|2e-jω2y2dy=(i-j)e-iω1x1e-|ω|2/4e-jω2x2.

In the following, we summarize the elementary properties of quaternion convolution as shown in Table 1 (compared to Folland ).

Basic properties of quaternion convolution.

Basic property Quaternion convolution
Linearity ( κ 1 f + κ 2 g ) h = κ 1 ( f h ) + κ 2 ( f h ) , κ 1 , κ 2 ℍ,
h ( κ 1 f + κ 2 g ) = κ 1 ( h f ) + κ 2 ( h g ) , κ 1 , κ 2
Shifting ( τ a f g ) = τ a ( f g ),
( f τ a g ) = τ a ( f g )
Conjugation ( f g ) ¯ = ( g - f - )
Associativity ( f g ) h = f ( g h )
Distributivity f ( g + h ) = ( f g ) + ( f h )
Vector differential a · ( f g ) = ( a · f ) g = f ( a · g )
Impulse convolution f δ = f
Lemma 8 (linearity).

For quaternion functions f,g, and h and quaternion constants κ1 and κ2 one gets (32)(κ1f+κ2g)h=κ1(fh)+κ2(gh). One also gets for real constants κ1 and κ2 (due to the noncommutativity of the quaternion multiplication, (33) does not hold for quaternion constants κ1 and κ2) (33)h(κ1f+κ2g)=κ1(hf)+κ2(hg).

Lemma 9 (shifting).

Given a quaternion function fL2(2;), let τaf(x) denote the shifted (translated) function defined by τaf(x)=f(x-a), where a2. Then one gets (34)(τafg)(x)=τa(fg)(x),(35)(fτag)(x)=τa(fg)(x).

Proof.

For (34), a direct calculation gives (36)(τafg)(x)=2τaf(y)g(x-y)dy=2f(y-a)g(x-y)dy=2f(y-a)g((x-a)-(y-a))dy=2f(s)g((x-a)-s)ds=(fg)(x-a)=τa(fg)(x), which finishes the proof.

Remark 10.

From (34) and (35), it is not difficult to see that (τafg)(x)=(fτag)(x) and (τagf)(x)=(gτaf)(x).

Lemma 11 (conjugation).

For all quaternion functions f,gL2(2;) one has (37)(fg)¯(x)=(g-f-)(x).

Proof.

A straightforward computation gives (38)(fg)¯(x)=2f(y)g(x-y)¯dy=2g-((x-y))f-(y)dy=2g-(z)f-((x-z))dz=(g-f-)(x). This finishes the proof.

Lemma 12 (see [<xref ref-type="bibr" rid="B19">2</xref>, <xref ref-type="bibr" rid="B22">3</xref>]; vector differential).

For all quaternion functions f,gL2(2;) one has (39)a·(fg)=(a·f)g=f(a·g).

Ell and Sangwine  distinguish between right and left discrete quaternion convolution due to the non-commutative property of the quaternion multiplication. Here, we only consider one kind of quaternion convolutions. We come now to the main theorem (generalization of the QFT of the quaternion convolution in general geometric Fourier transform is investigated in . It can easily be seen that the result is closely related to equation (4.30) of ) of this paper. This theorem describes the relationship between the convolution of two quaternion functions and its QFT.

Theorem 13.

Let f(x)=f0(x)+if1(x)+jf2(x)+kf3(x) and g(x)=g0(x)+ig1(x)+jg2(x)+kg3(x) be two quaternion-valued functions, then the QFT of the convolution of fL2(2;) and gL2(2;) is given by (40)q{fg}(ω)=(q{f0}(ω)+iq{f1}(ω))×(q{g0}(ω)+jq{g2}(-ω1,ω2))+(q{f0}(ω1,-ω2)+iq{f1}(ω1,-ω2))×(iq{g1}(ω)+kq{g3}(-ω1,ω2))+(jq{f2}(-ω1,ω2)+kq{f3}(-ω1,ω2))×(q{g0}(-ω1,ω2)+jq{g2}(ω))+(jq{f2}(-ω)+kq{f3}(-ω))×(iq{g1}(-ω1,ω2)+kq{g3}(ω)).

Proof.

In this proof we will use the decomposition of quaternion functions and their QFTs. Let q{f}(ω) and q{g}(ω) denote the QFT of fL2(2;) and gL2(2;), respectively. Expanding the QFT of the left-hand side of (40), we immediately get (41)q{fg}(ω)=(21)2e-iω1x1(fg)(x)e-jω2x2dx=(29)2e-iω1x1[2f(y)g(x-y)dy]e-jω2x2dx=2e-iω1x1f(y)[2g(x-y)e-jω2x2dx]dy.

By the change of variables z=x-y, the above transform can be written as (42)q{fg}(ω)=2e-iω1(y1+z1)f(y)g(z)e-jω2(y2+z2)dzdy=2e-iω1(y1+z1)({f0(y)+if1(y)}+jf2(y)+kf3(y))×2(+kg3(z){g0(z)+jg2(z)}+ig1(z)+kg3(z))g(z)e-jω2(y2+z2)dzdy=2e-iω1(y1+z1)(f0(y)+if1(y))×(g0(z)+jg2(z))e-jω2(y2+z2)dzdy+2e-iω1(y1+z1)(f0(y)+if1(y))×(ig1(z)+kg3(z))e-jω2(y2+z2)dzdy+2e-iω1(y1+z1)(jf2(y)+kf3(y))×(g0(z)+jg2(z))e-jω2(y2+z2)dzdy+2e-iω1(y1+z1)(jf2(y)+kf3(y))×(ig1(z)+kg3(z))e-jω2(y2+z2)dzdy=2e-iω1y1(f0(y)+if1(y))×(q{g0}(ω)+jq{g2}(-ω1,ω2))×e-jω2y2dzdy+2e-iω1y1(f0(y)+if1(y))×(iq{g1}(ω)+kq{g3}(-ω1,ω2))×e-jω2y2dzdy+2e-iω1y1(jf2(y)+kf3(y))×(q{g0}(-ω1,ω2)+jq{g2}(ω))×e-jω2y2dzdy+2e-iω1y1(jf2(y)+kf3(y))×(iq{g1}(-ω1,ω2)+kq{g3}(ω))×e-jω2y2dzdy=2e-iω1y1(f0(y)+if1(y))e-jω2y2dy×(q{g0}(ω)+jq{g2}(-ω1,ω2))+2e-iω1y1(f0(y)+if1(y))ejω2y2dy×(iq{g1}(ω)+kq{g3}(-ω1,ω2))+2e-iω1y1(jf2(y)+kf3(y))e-jω2y2dy×(q{g0}(-ω1,ω2)+jq{g2}(ω))+2e-iω1y1(jf2(y)+kf3(y))ejω2y2dy×(iq{g1}(-ω1,ω2)+kq{g3}(ω)), where the assumption q{gi}L2(2;) for i=1,2,3 is used in the fourth line. This gives the desired result.

The following lemmas are special cases of Theorem 13.

Lemma 14.

Let f,gL2(2;), where (43)f=f0+if1+jf2+kf3,g=g0+ig1+jg2+kg3. If q{g}L2(2;), then (40) takes the form (44)q{fg}(ω)=(q{f0}(ω)+iq{f1}(ω))q{g}(ω)+(jq{f2}(-ω1,ω2)+kq{f3}(-ω1,ω2))q{g}(-ω1,ω2). On the other hand, if q{f}L2(2;), then (45)q{fg}(ω)=q{f}(ω)(q{g0}(ω)+jq{g2}(-ω1,ω2))+q{f}(ω1,-ω2)(iq{g1}(ω)+kq{g3}(-ω1,ω2)).

Proof.

We only prove expression (44) of Lemma 14, with the other being similar. Following the steps of (41) we immediately get (46)q{fg}(ω)=2e-iω1(y1+z1)f(y)g(z)e-jω2(y2+z2)dzdy=2e-iω1(y1+z1)×({f0(y)+if1(y)}+jf2(y)+kf3(y))×g(z)e-jω2(y2+z2)dzdy=2e-iω1(y1+z1)(f0(y)+if1(y))×g(z)e-jω2(y2+z2)dzdy+2e-iω1(y1+z1)(jf2(y)+kf3(y))×g(z)e-jω2(y2+z2)dzdy=2e-iω1y1(f0(y)+if1(y))q{g}(ω)e-jω2y2dy+2e-iω1y1(jf2(y)+kf3(y))×q{g}(-ω1,ω2)e-jω2y2dy=(q{f0}(ω)+iq{f1}(ω))q{g}(ω)+(jq{f2}(-ω1,ω2)+kq{f3}(-ω1,ω2))×q{g}(-ω1,ω2), which was to be proven.

Lemma 15.

Let f,gL2(2;), where (47)f=f0+if1,g=g0+ig1+jg2+kg3. If q{f},q{g}L2(2;), then (48)q{fg}(ω)=q{f}(ω)q{g}(ω)=q{g}(ω)q{f}(ω), which is of the same form as a convolution of the classical Fourier transform .

Remark 16.

It is important to notice that, if f,gL2(2;), where (49)f=jf2+kf3,g=g0+ig1+jg2+kg3, then Lemma 15 reduces to (50)q{fg}(ω)=q{f}(ω)q{g}(-ω1,ω2)=q{g}(-ω1,ω2)q{f}(ω), where q{f},q{g}L2(2;).

Table 2 compares convolution theorems of the QFT and classical FT for f,gL2(2;).

Comparison of convolution theorems of the QFT and classical FT for f,gL2(2;).

Assumptions on quaternion functions QFT of convolution
q { f } , q { g } L 2 ( 2 ; ) QFT classical FT
q { f } L 2 ( 2 ; ) and q{g}L2(2;) QFT classical FT
q { f } L 2 ( 2 ; ) and q{g}L2(2;) QFT classical FT
f = f 0 + i f 1 and q{f},q{g}L2(2;) QFT = classical FT
f = j f 2 + k f 3 and q{f},q{g}L2(2;) QFT = classical FT

The following theorem is useful for solving the heat equation in quaternion algebra.

Theorem 17.

If f,gL2(2;) and q{g}L2(2;), then (51)q-1[q{f}q{g}](x)=(f0g)(x)+(if1g)(x)+(jf2g)(-x1,x2)+(kf3g)(-x1,x2).

Proof.

By the QFT inversion, we get, after some simplification, (52)q-1[q{f}q{g}](x)=(23)1(2π)22eiω1x1e-iω1y1f(y)e-jω2y2dy×q{g}(ω)ejω2x2dω=1(2π)22eiω1(x1-y1)f(y)×q{g}(ω)ejω2(x2-y2)dydω=1(2π)22eiω1(x1-y1)×(f0(y)+if1(y)+jf2(y)+kf3(y))×q{g}(ω)ejω2(x2-y2)dydω=1(2π)22eiω1(x1-y1)(f0(y)+if1(y))+eiω1(x1-y1)(jf2(y)+kf3(y))×q{g}(ω)ejω2(x2-y2)dydω=2(f0(y)+if1(y))g(x-y)dy+2(jf2(y)+kf3(y))g(-(x1-y1),x2-y2)dy=(f0g)(x)+(if1g)(x)+(jf2g)(-x1,x2)+(kf3g)(-x1,x2), where, in the second line, we have used the assumption q{g}L2(2;). This completes the proof of (51).

As an immediate consequence of Theorem 17, we get the following corollaries.

Corollary 18.

If f,gL2(2;) and q{g}L2(2;), where (53)f=f0+if1,g=g0+ig1+jg2+kg3, then (51) reduces to (54)q-1[q{f}q{g}](x)=(fg)(x).

Corollary 19.

Let (55)f(x)={e-(x1+x2),ifx1>0,x2>0,0,otherwise. And consider the quaternionic Gabor filter (56)g(x)=eiu0x1ejv0x2e-(1/2)|x|2.

Then, (57)q{fg}(ω)=e-(1/2)((ω1-u0)2+(ω2-v0)2)(1-iω1-jω2-kω1ω2)(2π)2(1+ω12+ω22+ω12ω22).

Proof.

The QFT of f is given by (58)q{f}(ω)=1-iω1-jω2-kω1ω2(2π)2(1+ω12+ω22+ω12ω22), and the QFT of g is given by (59)q{g}(ω)=e-(1/2)((ω1-u0)2+(ω2-v0)2). Therefore, using Corollary 18, we obtain (57).

4. Applications of QFT

In , the authors proposed to use quaternions in order to define a Fourier transform applicable to color images. Their framework makes it possible to compute a single, holistic, Fourier transform which treats a color image as a vector field. In image processing, taking a given image as the initial value, the forward solution to the heat equation or a diffusion equation in general, produces blurred images and the backward solution produces sharpen images for example, see [21, pages 342–350].

In this section, we present two applications of QFT to partial differential equations in quaternion algebra.

4.1. Hypoellipticity

In this paper, since we only deal with QFT in the L2(2;) framework, we will discuss the hypoellipticity in this framework; that is, we will only deal with L2(2;) solutions for linear partial differential operators with constant quaternion coefficients: (60)P()=0|α|naααbα,aα,bα. The noncommutativity of quaternion gives different aspects of P() with constant complex coefficients aα,bα.

Example 20.

Let P()=x1b, b, and f(x),g(x)C1(2;).

Since x1(f(x)g(x))=(x1f(x))g(x)+f(x)(x1g(x)), we have P()(f(x)g(x))=(P()f(x))g(x)+f(x)(P()g(x)) when b. But, when b, as bf(x)f(x)b in general, we cannot have P()(f(x)g(x))=(P()f(x))g(x)+f(x)(P()g(x)) in general.

We have P()fg=fP()g when b. But, by the same reason as (i), when b, we cannot have P()fg=fP()g in general.

Let us start with the definition of our L2(2;) version of hypoellipticity (compared to [22, page 110]).

Definition 21.

The linear partial differential operator P() in 2 is said to be L2(2;)-hypoelliptic if, given any subset U of 2 and any solution u in L2(2;) such that P()u is a C function in U, then all its components u(ui,i=0,1,2,3) are a C function in U.

Definition 22.

Given a linear partial differential operator P() of (60) with the quaternion constant coefficients. One says that a solution E(x) of P()u=δ, where δ is the delta function, is called a fundamental solution of P().

Let A and B be subsets of 2. Define the sum A+B by A+B={x+y2;xA,yB}.

Theorem 23.

Assume that there is one fundamental solution E(x) of P() which is a C function in 2{0}, and the identities (61)P()(fg)=P()fg=fP()g are satisfied for arbitrary sufficiently smooth quaternion-valued functions f and g such that gf is a compactly supported C quaternion function with aα of P() being quaternion constant coefficients and bα of P() being real constant coefficients. Then, the linear partial differential operator P() is L2(2;)-hypoelliptic in 2.

Proof.

Firstly, let U be an arbitrary open subset of 2 and u a solution in U with values in such that f=P()u is a C function in U. Let x0 be an arbitrary point in U. It will suffice to show that u is a C function in some open neighborhood of x0. Take an open disc Dη(x0)={x2;|x-x0|<η} such that Dη(x0)¯U. There exists a function gC(U;) such that suppgU and g=1 in Dη(x0). Then, we have (62)P()(gu)=gP()u+v=gf+v, where every term of v contains a derivative of g of nonzero order; therefore v=0, where the derivatives of g vanish, especially in Dη(x0) and outside of suppg. For the fundamental solution E(x), the hypothesis (61) implies (63)EP()(gu)={P()E}(gu)=gu. Hence (64)gu=E(gf)+Ev. But gf is a compactly supported C function and the convolution of any function with any compactly supported C function is a C function. Therefore, it suffices to show that Ev is a C function in an open neighborhood of x0, because gu is also a C function in an open neighborhood of x0 and gu=u in Dη(x0).

Finally, we will show that Ev is a C function in an open neighborhood of x0. Let us select ε>0 such that ε<1/2η. Then, the open disc Dε(x0) is a neighborhood of x0. Let ζε(x)C(2;), another cutoff function, be equal to one for |x|<ε/2 and to zero for |x-x0|>ε. We have (65)Ev=(ζεE)v+{(1-ζε)E}v. The hypothesis implies that (1-ζε)EC(2;), and therefore (1-ζε)EvC(2;). Since (66)supp{(ζεE)v}supp(ζεE)+suppv,supp{(ζεE)v} is contained in the ε-neighborhood of suppv. We have already seen that v=0Dη(x0). Hence, (ζεE)v vanishes in Dε(x0), and, therefore, Ev is a C function in Dε(x0).

4.2. Parabolic Initial Value Problem

Let us consider the parabolic initial value problem (67)tu-2u=0,on2×(0,), with (68)u(x,0)=f(x),f𝒮(2;), where 𝒮(2;) is the quaternion Schwartz space. Applying the QFT, we easily obtain (69)q{ut}=(iω1)2q{u}(ω)+q{u}(ω)(jω2)2=-|ω|2q{u}(ω). The general solution of (69) is given by (70)q{u}(ω,t)=Ce-|ω|2t, where C is a quaternion constant. We impose the initial condition q{u}(ω,0)=q{f}(ω) to obtain (71)q{u}(ω,t)=e-|ω|2tq{f}(ω). Notice that the QFT of a Gaussian quaternion function is also a Gaussian quaternion function (compared to Bahri et al. ). Hence (72)14πtq{(e-|x|2/(4t))}=e-|ω|2t. Applying the inverse QFT, we have (73)u(x,t)=q-1[e-|ω|2tq{f}](x)=14πtq-1[q{f}q{e-|x|2/(4t)}](x). Since (74)q{e-|x|2/(4t)}(ω)=4πte-|ω|2tL2(2;), then we can apply the convolution theorem of (51) to get (75)u(x,t)=(f0Kt)(x)+(if1Kt)(x)+(jf2Kt)(-x1,x2)+(kf3Kt)(-x1,x2), where Kt=(1/4πt)e-|x|2/(4t), and fiL2(2;), i=0,1,2,3. By Definition 6 of the convolution, we finally obtain (76)u(x,t)=14πt2f0(y)e-((x1-y1)2+(x2-y2)2)/(4t)dy+i4πt2f1(y)e-((x1-y1)2+(x2-y2)2)/(4t)dy+j4πt2f2(y)e-((-x1-y1)2+(x2-y2)2)/(4t)dy+k4πt2f3(y)e-((-x1-y1)2+(x2-y2)2)/(4t)dy. In an actual application, one often takes the quaternionic Gabor filter (see [6, 10]) as (77)f(x)=eiu0x1e-(x12+x22)ejv0x2. Therefore, the above identity will reduce to (78)u(x,t)=14πt2cos(u0x1)cos(v0x2)e-(x12+x22)×e-((x1-y1)2+(x2-y2)2)/(4t)dy,+i4πt2sin(u0x1)cos(v0x2)e-(x12+x22)×e-((x1-y1)2+(x2-y2)2)/(4t)dy,+j4πt2cos(u0x1)sin(v0x2)e-(x12+x22)×e-((-x1-y1)2+(x2-y2)2)/(4t)dy,+k4πt2sin(u0x1)sin(v0x2)e-(x12+x22)×e-((-x1-y1)2+(x2-y2)2)/(4t)dy.

5. Conclusion

Due to the non-commutative property of quaternion multiplication, there are three different types of two-dimensional QFTs. These three QFTs are the so-called left-sided QFT, right-sided QFT, and a double-sided QFT, respectively. In this work, we have established convolution theorem of the double-sided QFT applied to real fields f:2 and quaternion fields f:2. Some important properties of the QFT convolution are investigated. We have shown that the QFT convolution is useful to study hypoellipticity and to solve the heat equation in quaternion algebra framework. It can easily be seen that the solution of generalized heat equation is extension of solution of the classical heat equation.

The future work will establish the convolution theorems of the right-sided QFT. We compare some properties of the convolution theorems of the two types of QFTs. We will apply the properties to find the solution of partial differential equations in quaternion algebra framework. The solutions of generalized partial differential equations using the properties of the three types of two-dimensional QFTs will be compared too.

Acknowledgments

The authors would like to thank the reviewers whose deep and extensive comments greatly contributed to improve this paper. The first author is partially supported by Hibah Penelitian Kompetisi Internal Tahun 2013 (no. 110/UN4-.42/LK.26/SP-UH/2013) from the Hasanuddin University, Indonesia. The second author is partially supported by JSPS.KAKENHI (C)22540130 and (C)25400202 of Japan and the third author is partially supported by NSERC of Canada.

Ebling J. Scheuermann G. Clifford Fourier transform on vector fields IEEE Transactionson Visualization and Computer Graphics 2005 11 4 469 479 10.1109/TVCG.2005.54 Hitzer E. Mawardi B. Clifford Fourier transform on multivector fields and uncertainty principles for dimensions n=2(mod4) and n=3(mod4) Advances in Applied Clifford Algebras 2008 18 3-4 715 736 10.1007/s00006-008-0098-3 MR2490585 ZBL1177.15029 Mawardi B. Hitzer E. Clifford Fourier transformation and uncertainty principle for the Clifford geometric algebra Cl3,0 Advances in Applied Clifford Algebras 2006 16 1 41 61 10.1007/s00006-006-0003-x MR2226529 ZBL1133.42305 de Bie H. de Schepper N. Fractional Fourier transform of hypercomplex signals Signal, Image and Video Processing 2012 6 3 381 388 10.1007/s11760-012-0315-3 Guanlei X. Xiaotong W. Xiagang X. Fractional quaternion Fourier transform, convolutionand correlation Signal Processing 2008 88 10 2511 2517 10.1016/j.sigpro.2008.04.012 Bahri M. Hitzer E. Ashino R. Vaillancourt R. Windowed Fourier transform of two-dimensional quaternionic signals Applied Mathematics and Computation 2010 216 8 2366 2379 10.1016/j.amc.2010.03.082 MR2647110 ZBL1196.42009 Bahri M. Ashino R. Vaillancourt R. Two-dimensional quaternion wavelet transform Applied Mathematics and Computation 2011 218 1 10 21 10.1016/j.amc.2011.05.030 MR2821445 ZBL1232.65192 Bahri M. Quaternion algebra-valued wavelet transform Applied Mathematical Sciences 2011 5 71 3531 3540 MR2860561 ZBL1247.42032 Sangwine S. J. Biquaternion (complexified quaternion) roots of -1 Advances in Applied Clifford Algebras 2006 16 1 63 68 10.1007/s00006-006-0005-8 MR2226530 ZBL1107.30038 Bülow T. Hypercomplex spectral signal representations for the processing and analysis of images [Ph.D. thesis] 1999 Kiel, Germany University of Kiel Ell T. A. Quaternion-Fourier transformations for analysis of two-dimensional lineartime-invariant partial differential systems 2 Proceedings of the 32nd IEEE Conference on Decision and Control December 1993 San Antonio, Tex, USA 1830 1841 Hitzer E. Quaternion Fourier transform on quaternion fields and generalizations Advances in Applied Clifford Algebras 2007 17 3 497 517 10.1007/s00006-007-0037-8 MR2350594 ZBL1143.42006 Bahri M. Hitzer E. Hayashi A. Ashino R. An uncertainty principle for quaternion Fourier transform Computers & Mathematics with Applications 2008 56 9 2411 2417 10.1016/j.camwa.2008.05.032 MR2466762 ZBL1165.42310 Assefa D. Mansinha L. Tiampo K. F. Rasmussen H. Abdella K. Local quaternion Fourier transform and color images texture analysis Signal Processing 2010 90 6 1825 1835 10.1016/j.sigpro.2009.11.031 Bujack R. Scheuermann G. Hitzer E. A general geometric fourier transform convolution theorem Advances in Applied Clifford Algebras 2013 23 1 15 38 10.1007/s00006-012-0338-4 ZBL1263.15021 Bujack R. Scheuermann G. Hitzer E. Guerlebeck K. A general geometric Fourier transform Proceedings of the 9th International Conference on Clifford Algebras and Their Applications in Mathematical Physics (ICCA '11) July 2011 Weimar, Germany 10.1007/s00006-012-0338-4 ZBL 1263.15021 Hitzer E. Directional uncertainty principle for quaternion Fourier transform Advances in Applied Clifford Algebras 2010 20 2 271 284 10.1007/s00006-009-0175-2 MR2645348 ZBL1198.42006 Folland G. B. Real Analysis: Modern Techniques and Their Applications 1999 2nd New York, NY, USA John Wiley & Sons Pure and Applied Mathematics Ell T. A. Sangwine S. J. Hypercomplex Fourier transforms of color images IEEE Transactions on Image Processing 2007 16 1 22 35 10.1109/TIP.2006.884955 MR2460142 Mallat S. A Wavelet Tour of Signal Processing 1999 2nd San Diego, Calif, USA Academic Press MR2479996 Petrou M. Petrou C. Image Processing: The Fundamentals 2010 2nd West Sussex, UK John Wiley & Sons Hörmander L. The Analysis of Linear Partial Differential Operators I: Distribution Theory and Fourier Analysis 1983 Berlin, Germany Springer 10.1007/978-3-642-96750-4 MR717035