Journal of Signal and Information Processing
Vol.09 No.03(2018), Article ID:86760,15 pages
10.4236/jsip.2018.93012

The Properties and Fast Algorithm of Quaternion Linear Canonical Transform

Ye Zhang1, Guanlei Xu2*

1Nanjing Changjiang Electronics Group Co. Ltd., Nanjing, China

2Dalian Navy Academy, Dalian, China

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: June 28, 2018; Accepted: August 18, 2018; Published: August 21, 2018

ABSTRACT

The quaternion linear canonical transform (QLCT) is defined in this paper, with proofs given for its reversibility property, its linear property, its odd-even invariant property and additivity property. Meanwhile, the quaternion convolution (QCV), quaternion correlation (QCR) and product theorem of LCT are deduced. Their physical interpretation is given as classical convolution, correlation and product theorem. Moreover, the fast algorithm of QLCT (FQLCT) is obtained, whose calculation complexity for different signals is similar to FFT. In addition, the paper presents the relationship between the convolution and correlation in LCT domains, and the convolution and correlation can be calculated via product theorem in Fourier transform domain using FFT.

Keywords:

Quaternion Signals (Hyper-Complex Signals), LCT, Convolution, Correlation

1. Introduction

The linear canonical transform (LCT) is a new tool that comes into being in signal processing [1]- [32]. The LCT is the generalization of the FRFT and so on [2][3][12]. Up till now there have been a lot of papers involving the FRFT and the LCT, such as papers [1]- [10]. However, none of them has involved the LCT of quaternion signals (or Hyper-complex signals) even if there has been similar work on FRFT [5]. Quaternion signals can be taken as the generalization of scalar, complex signals and vector, and after the introduction of quaternion signals by Hamilton in 1843 [11]it has become one basic tool for multi-channel and multi-dimensional space. For example, grey image [30]can be taken as scalar, and the analytic signal after Hilbert transformation [12][13][14][15][16][29]is complex signal. The color image can be taken as one vector [17][18], a quaternion number whose real part is zero. In [19], the transform, convolution and correlation have been addressed in fractional Fourier transform (FRFT) domain. In this paper we first propose the definition of the QLCT, QCV and QCR in the LCT domain for quaternion signals, which are the generalization of those in [5]. Meanwhile, some properties and the fast algorithm of QLCT are discussed. We also discover the relationship of QCV and QCR in the LCT domain for quaternion signals. We found that QCV and QCR can be implemented via product theorem in the QLCT domain. Thus we not only yield the generalized frame for scalar, complex signal, vector and quaternion signal [17][20]in the QLCT domain, but also give one new idea and one theoretical base for future engineering use.

In the rest of this paper, we will introduce the definition of QLCT in Section II. We will show the properties in Section III. In Section IV, FRQCV and FRQCR will be addressed. Section V is the fast algorithm. The last section concludes our paper.

2. Definitions of QLCT

For convenience of discussion, we first give some notations used in the following of this paper. f ( x , y ) denotes 2D signal in time domain; F is classical Fourier transform operator; F L ( a , b , c , d ) ( F L in short) is 1D LCT operator, and F L ( u ) is the 1D LCT of f ( x , y ) ; F L 1 , L 2 is the 2D LCT operator of f ( x , y ) ; F Q is classical quaternion Fourier transform operator, and F Q ( u , v ) is quaternion Fourier transform of f ( x , y ) ; I is equivalence operator; P is odd-even operator; “*” is classical convolution operator; “ ” is conjugation operator. “N” is integer set; “R” is real set. Define the product operator of two LCTs’ transform parameter systems:

L 1 L 2 = L 1 ( a 1 , b 1 , c 1 , d 1 ) L 2 ( a 2 , b 2 , c 2 , d 2 ) = L ( a , b , c , d )

where [ a b c d ] = [ a 1 b 1 c 1 d 1 ] [ a 2 b 2 c 2 d 2 ] . Quaternion signals are also called Hypercomplex signals, which are the generalization of complex signals. Complex signals have two components: the real part and the imaginary part. However, one quaternion signal has four parts, one real component and three imaginary parts:

q = q r + i q i + j q j + k q k (1)

where q r , q i , q j , q k R , i , j , k are three imaginary units, which satisfy the following relations: i 2 = j 2 = k 2 = 1 , i j = j i = k , j k = k j = i , k i = i k = j . q a = q r + i q i , q b = q i + i q k . If q r = 0 , then q = i q i + j q j + k q k is called vector, and q r is called scalar. q a and q b are complex signals. Since the sequences of i, j and k will affect the result, the definition of QLCT would take them into account.

Definition 1: For any quaternion signal f ( x , y ) = f r ( x , y ) + i f i ( x , y ) + j f j ( x , y ) + k f k ( x , y ) ( f r ( x , y ) , f i ( x , y ) , f j ( x , y ) , f k ( x , y ) are real ones), the QLCT of f ( x , y ) is F i , j L 1 , L 2 ( u , v )

F i , j L 1 , L 2 ( u , v ) = F i , j L 1 , L 2 { f ( x , y ) } = F i , j L 1 , L 2 { f ( x , y ) } ( u , v ) = + + K L 1 , i ( x , u ) f ( x , y ) K L 2 , j ( y , v ) d x d y (2-1)

where, K L 1 , i ( x , u ) = 1 2 π b 1 i exp ( i ( a 1 + d 1 ) u 2 2 b 1 i u x b 1 ) , K L 2 , j ( y , v ) = 1 2 π b 2 j exp ( j ( a 2 + d 2 ) v 2 2 b 2 j v y b 2 ) . Meanwhile, in the following of this paper we assume a 1 d 1 b 1 c 1 = 1 , a 2 d 2 b 2 c 2 = 1 and b 1 , b 2 0 .

The reversibility transform is defined as

F i , j L 1 , L 2 ( u , v ) = F i , j L 1 , L 2 { f ( x , y ) } = F i , j L 1 , L 2 { f ( x , y ) } ( u , v ) = + + K L 1 , i ( x , u ) f ( x , y ) K L 2 , j ( y , v ) d x d y (2-2)

where, K L 1 , i ( x , u ) = 1 2 π b 1 i exp ( i ( a 1 + d 1 ) u 2 2 b 1 + i u x b 1 ) , K L 2 , j ( y , v ) = 1 2 π b 2 j exp ( j ( a 2 + d 2 ) v 2 2 b 2 + j v y b 2 ) .

If ( a 1 , b 1 , c 1 , d 1 ) = ( a 2 , b 2 , c 2 , d 2 ) = ( 0 , 1 , 1 , 0 ) , definition 1 is quaternion Fourier transform; if ( a 1 , b 1 , c 1 , d 1 ) = ( 0 , 1 , 1 , 0 ) , ( a 2 , b 2 , c 2 , d 2 ) = ( 1 , 0 , 0 , 1 ) , definition 1 is classical 1D Fourier transform of f ( x , y ) for variable x; if ( a 1 , b 1 , c 1 , d 1 ) = ( 1 , 0 , 0 , 1 ) , ( a 2 , b 2 , c 2 , d 2 ) = ( 0 , 1 , 1 , 0 ) , definition 1 is classical 1D Fourier transform of f ( x , y ) for variable y; if ( a 1 , b 1 , c 1 , d 1 ) = ( a 2 , b 2 , c 2 , d 2 ) = ( 1 , 0 , 0 , 1 ) , definition 1 is equivalence transform of f ( x , y ) . As shown above, definition 1 is the generalization of the fractional quaternion Fourier transform and the quaternion Fourier transform [18][19][20][21]. The reversibility (or reconstruction) is one important property for one transform, especially for the processing in another domain. The following gives the proof of the reversibility property.

Theorem 1: One quaternion f ( x , y ) can be reconstructed from F i , j L 1 , L 2 ( u , v ) via QLCT.

Proof: The proof is trivial and omitted here.

3. The Properties of QLCT

In the following section we list the properties and present the proof.

Property 1: For any one quaternion signal f n ( x , y ) ( n ) , the following relationship is true: F i , j L 1 , L 2 { a n f n ( x , y ) } = a n F i , j L 1 , L 2 { f n ( x , y ) } ( a n ) .

Proof: Since QLCT is one linear transform, property 1 can be easily obtained from definition 1.

Property 2: F i , j L 3 , L 4 F i , j L 1 , L 2 = F i , j L 1 , L 2 F i , j L 3 , L 4 = F i , j L 1 L 3 , L 2 L 4 . Proof: For any one quaternion signal f ( x , y ) , from definition 1 we can obtain

F i , j L 3 , L 4 F i , j L 1 , L 2 { f ( x , y ) } = + + K L 3 , i ( u , s ) { + + K L 1 , i ( x , u ) f ( x , y ) K L 2 , j ( y , v ) d x d y } K L 4 , j ( v , w ) d u d v = + + { + + K L 3 , i ( u , s ) K L 1 , i ( x , u ) f ( x , y ) K L 2 , j ( y , v ) K L 4 , j ( v , w ) d u d v } d x d y = + + { [ + K L 3 , i ( u , s ) K L 1 , i ( x , u ) d u ] f ( x , y ) [ + K L 2 , j ( y , v ) K L 4 , j ( v , w ) d v ] } d x d y (3)

For 1D signal the right formula is true [2]:

+ K L 2 ( u , u ) K L 1 ( u , u ) d u = K L 2 L 1 ( u , u ) (4)

Substitute (4) into (3):

F i , j L 3 , L 4 F i , j L 1 , L 2 { f ( x , y ) } = + + { K L 3 L 1 , i ( x , s ) f ( x , y ) K L 4 L 2 , j ( y , w ) } d x d y = F i , j L 1 L 3 , L 2 L 4 { f ( x , y ) }

Therefore,

F i , j L 3 , L 4 F i , j L 1 , L 2 = F i , j L 1 L 3 , L 2 L 4 (5)

The result can be obtained similarly:

F i , j L 1 , L 2 F i , j L 3 , L 4 = F i , j L 1 L 3 , L 2 L 4 (6)

From (5) (6): F i , j L 3 , L 4 F i , j L 1 , L 2 = F i , j L 1 , L 2 F i , j L 3 , L 4 = F i , j L 1 L 3 , L 2 L 4

Property 3: F i , j L 3 , L 4 F i , j L 1 , L 2 = F i , j L 1 , L 2 F i , j L 3 , L 4 , F i , j L 5 , L 6 ( F i , j L 3 , L 4 F i , j L 1 , L 2 ) = ( F i , j L 5 , L 6 F i , j L 3 , L 4 ) F i , j L 1 , L 2 .

Proof: This property can be obtained from property 2.

Property 4: If F i , j L 1 , L 2 { f ( x , y ) } = F i , j L 1 , L 2 ( u , v ) , then F i , j L 1 , L 2 { f ( x , y ) } = F i , j L 1 , L 2 ( u , v ) , F i , j L 1 , L 2 { f ( x , y ) } = F i , j L 1 , L 2 ( u , v ) , F i , j L 1 , L 2 { f ( x , y ) } = F i , j L 1 , L 2 ( u , v ) .

Proof: Let A 1 = 2 π b 1 i , A 2 = 2 π b 2 j , C 1 = d 1 2 b 1 , C 2 = d 2 2 b 2 , and insert them into (2):

F i , j L 1 , L 2 { f ( x , y ) } = A 1 e i u 2 C 1 + + e i u x b 1 e i x 2 C 1 f ( x , y ) e j y 2 C 1 e j v y b 2 d x d y A 2 e j v 2 C 2 (7)

Let s = x , z = y , and substitute them in (7):

F i , j L 1 , L 2 { f ( x , y ) } = F i , j L 1 , L 2 { f ( s , z ) } = A 1 e i u 2 C 1 + + e i u ( s ) b 1 e i s 2 C 1 f ( s , z ) e j z 2 C 2 e j v ( z ) b 2 d s d z A 2 e j v 2 C 2 = A 1 e i ( u ) 2 C 1 + + e i ( u ) s b 1 e i s 2 C 1 f ( s , z ) e j z 2 C 2 e j ( v ) z b 2 d s d z A 2 e j ( v ) 2 C 2 = F i , j L 1 , L 2 ( u , v )

It can be obtained as well: F i , j L 1 , L 2 { f ( x , y ) } = F i , j L 1 , L 2 ( u , v ) and F i , j L 1 , L 2 { f ( x , y ) } = F i , j L 1 , L 2 ( u , v ) .

We can draw the conclusion that transformed signal of the odd is odd, and even is even.

Property 5: If n , then ( F i , j L 1 , L 2 ) n = F i , j ( L 1 ) n , ( L 2 ) n .

Proof: From property 2,

F i , j L 1 , L 2 F i , j L 1 , L 2 F i , j L 1 , L 2 n = F i , j L 1 × × L 1 n , L 2 × × L 2 n = F i , j ( p 1 ) n , ( p 2 ) n

then

( F i , j L 1 , L 2 ) n = F i , j ( L 1 ) n , ( L 2 ) n .

QLCT doesn’t satisfy Parseval’s principle. Meanwhile, it is hard to find one obvious relationship between QLCT and Wigner-Ville time-frequency plane. Some other properties [2]cannot find physical interpretation in QLCT domains.

4. FRQCV and FRQCR

Convolution and correlation play an important role in signal processing, especially for linear system design and filter design, etc. The convolution in time domain is to the product in Fourier transform domain, that is to say, the classical convolution in time domain can be implemented in Fourier transform domain via FFT, which is beneficial for real-time engineering use. In classical time-frequency analysis correlation is special convolution in that the original signals are implemented via conjugation and so on. This is very important for engineering use [17][20][24]. The key to this paper is to discover the relationships in fractional quaternion Fourier transform domain between them so that we can find the physical interpretation as that of the classical Fourier transform. Paper [26]yielded fractional convolution and product theorem for 1D signals first, however, it didn’t give the similar physical interpretation as that of the classical theorem. Later papers [27][28][29]obtained similar result as the classical theorems. However, they are only for 1D signals. In this section the QCV and QCR of the LCT would be discussed, and can be implemented via FFT.

4.1. Fractional Convolution and Product Theorem

In the following, four theorems are yielded, and theorem 2 and 3 are suitable for scalar and complex signals, and theorem 4 and 5 are suitable for scalar, complex signals, vector and quaternion signals.

Theorem 2: For any real scalar or complex signal f ( x , y ) and convolution kernel h ( x , y ) ,

g ( x , y ) = f ( x , y ) ¯ h ( x , y ) = ( f ¯ h ) ( x , y ) A 1 , 2 e i ( x 2 C 1 + y 2 C 2 ) ( e i ( x 2 C 1 + y 2 C 2 ) f ( x , y ) ) ( h ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) )

where B 1 = 1 / b 1 , B 2 = 1 / b 2 , A 1 , 2 = 1 2 π i b 1 b 2 , C 1 = a 1 2 b 1 , C 2 = a 2 2 b 2 , then:

F L 1 , L 2 { g ( x , y ) } = e i ( u 2 C 1 + v 2 C 2 ) F L 1 , L 2 { f ( x , y ) } F L 1 , L 2 { h ( x , y ) } (8)

Proof:

F L 1 , L 2 { g ( x , y ) } = + + K L 1 , L 2 ( x , y , u , v ) g ( x , y ) d x d y = A 1 , 2 e i ( u 2 C 1 + v 2 C 2 ) + + e i ( x 2 C 1 + y 2 C 2 ) A 1 , 2 e i ( x 2 C 1 + y 2 C 2 ) e i ( x u B 1 + y v B 2 ) ( e i ( x 2 C 1 + y 2 C 2 ) f ( x , y ) ) ( h ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) ) d x d y (9)

Substitute (9) with s = x τ , z = y η :

F L 1 , L 2 { g ( x , y ) } = A 1 , 2 2 e i ( u 2 C 1 + v 2 C 2 ) + + e i ( τ 2 C 1 + η 2 C 2 ) f ( τ , η ) e i ( τ u B 1 + η v B 2 ) { + + h ( s , z ) e i ( s 2 C 1 + z 2 C 2 ) e i ( s u B 1 + z v B 2 ) d s d z } d τ d η = A 1 , 2 + + e i ( τ 2 C 1 + η 2 C 2 ) f ( τ , η ) e i ( τ u B 1 + η v B 2 ) d τ d η F L 1 , L 2 { h ( x , y ) } = e i ( u 2 C 1 + v 2 C 2 ) F L 1 , L 2 { f ( x , y ) } F L 1 , L 2 { h ( x , y ) }

From theorem 2 it can be concluded that the convolution of scalar or complex signal is to the product, frequency-modulated by a chirp, of them in linear canonical transform.

Theorem 3: For any real scalar or complex signal f ( x , y ) and convolution kernel h ( x , y ) ,

g ( x , y ) = ( f ¯ ¯ h ) ( x , y ) ( e i ( x 2 C 1 + y 2 C 2 ) ) / ( e i ( x 2 C 1 + y 2 C 2 ) f ( x , y ) ) ( h ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) )

where B 1 = 1 / b 1 , B 2 = 1 / b 2 , A 1 , 2 = 1 2 π i b 1 b 2 , A 1 , 2 1 , 1 = 1 2 π i ( b 1 ) ( b 2 ) , C 1 = a 1 2 b 1 , C 2 = a 2 2 b 2 , then

F L 1 , L 2 { e i ( u 2 C 1 + v 2 C 2 ) f ( x , y ) g ( x , y ) } = A 1 , 2 1 , 1 ( f L 1 , L 2 ( u , v ) ¯ ¯ h L 1 , L 2 ( u , v ) ) (10)

Proof:

F L 1 , L 2 ( A 1 , 2 1 , 1 f L 1 , L 2 ( u , v ) ¯ ¯ h L 1 , L 2 ( u , v ) ) = A 1 , 2 1 , 1 e i ( x 2 C 1 + y 2 C 2 ) + + e i ( C 1 u 2 + C 2 v 2 ) { A 1 , 2 1 , 1 f L 1 , L 2 ( u , v ) ¯ ¯ h L 1 , L 2 ( u , v ) } e i ( x u B 1 + y v B 2 ) d u d v

Substitute s = x τ , z = y η in above equation

F L 1 , L 2 { A 1 , 2 1 , 1 f L 1 , L 2 ( u , v ) ¯ ¯ h L 1 , L 2 ( u , v ) } = ( A 1 , 2 1 , 1 ) 2 e i ( x 2 C 1 + y 2 C 2 ) 4 π 2 + + { + + h L 1 , L 2 ( s , z ) e i ( C 1 s 2 + C 2 z 2 ) e i ( x s B 1 + y z B 2 ) d s d z } e i ( x τ B 1 + y η B 2 ) f L 1 , L 2 ( τ , η ) e i ( C 1 τ 2 + C 2 η 2 ) d τ d η = h ( x , y ) A 1 , 2 1 , 1 + + e i ( x τ B 1 + y η B 2 ) f L 1 , L 2 ( τ , η ) e i ( C 1 τ 2 + C 2 η 2 ) d τ d η = f ( x , y ) h ( x , y ) e i ( C 1 x 2 + C 2 y 2 )

Therefore, F L 1 , L 2 { e i ( u 2 C 1 + v 2 C 2 ) f ( x , y ) g ( x , y ) } = A 1 , 2 1 , 1 ( f L 1 , L 2 ( u , v ) ¯ ¯ h L 1 , L 2 ( u , v ) ) .

Theorem 4: For one given quaternion function f ( x , y ) = f a ( x , y ) + f b ( x , y ) j and convolution kernel function

h ( x , y ) = h a ( x , y ) + h b ( x , y ) j

where f a ( x , y ) = f r ( x , y ) + i f i ( x , y ) , f b ( x , y ) = f j ( x , y ) + i f k ( x , y ) , h a ( x , y ) = h r ( x , y ) + i h i ( x , y ) , h b ( x , y ) = h j ( x , y ) + i h k ( x , y ) .

Set A 1 , 2 = 1 2 π i b 1 b 2 , A 1 , 2 1 , 1 = 1 2 π i ( b 1 ) ( b 2 ) , and define

g ( x , y ) = ( f ¯ h ) ( x , y ) A 1 , 2 e i ( x 2 C 1 + y 2 C 2 ) ( e i ( x 2 C 1 + y 2 C 2 ) f ( x , y ) ) ( h ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) )

where, C 1 = a 1 2 b 1 , C 2 = a 2 2 b 2 , B 1 = 1 / b 1 , B 2 = 1 / b 2 , α = arcsin b 1 , β = arcsin b 2 , then

F L 1 , L 2 { g ( x , y ) } = e i ( u 2 C 1 + v 2 C 2 ) { F L 1 , L 2 [ f a ( x , y ) ] F L 1 , L 2 [ h a ( x , y ) ] F L 1 , L 2 [ f b ( x , y ) ] F L 1 , L 2 [ h b ¯ ( x , y ) ] } + e i ( u 2 C 1 + v 2 C 2 α β ) { F L 1 , L 2 [ f a ( x , y ) ] F L 1 , L 2 [ h b ( x , y ) ] + F L 1 , L 2 [ f b ( x , y ) ] F L 1 , L 2 [ h a ¯ ( x , y ) ] } j (11)

Proof:

( e i ( x 2 C 1 + y 2 C 2 ) f ( x , y ) ) ( h ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) ) = ( e i ( x 2 C 1 + y 2 C 2 ) ( f a ( x , y ) + f b ( x , y ) j ) ) ( ( h a ( x , y ) + h b ( x , y ) j ) e i ( x 2 C 1 + y 2 C 2 ) ) = ( e i ( x 2 C 1 + y 2 C 2 ) f a ( x , y ) ) ( h a ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) ) + ( e i ( x 2 C 1 + y 2 C 2 ) f a ( x , y ) ) ( h b ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) ) j + ( e i ( x 2 C 1 + y 2 C 2 ) f b ( x , y ) ) ( h a ¯ ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) ) j ( e i ( x 2 C 1 + y 2 C 2 ) f b ( x , y ) ) ( h b ¯ ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) )

From theorem 2 it can be obtained

F L 1 , L 2 { A 1 , 2 e i ( x 2 C 1 + y 2 C 2 ) ( e i ( x 2 C 1 + y 2 C 2 ) f a ( x , y ) ) ( h a ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) ) } = e i ( u 2 C 1 + v 2 C 2 ) ( F L 1 , L 2 { f a ( x , y ) } F L 1 , L 2 { h a ( x , y ) } )

F L 1 , L 2 { A α , β e i ( x 2 C 1 + y 2 C 2 ) ( e i ( x 2 C 1 + y 2 C 2 ) f b ( x , y ) ) ( h b ¯ ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) ) } = e i ( u 2 C 1 + v 2 C 2 ) ( F L 1 , L 2 { f b ( x , y ) } F L 1 , L 2 { h b ¯ ( x , y ) } )

From the linear property of fractional Fourier transfor

F L 1 , L 2 { g ( x , y ) } = e i ( u 2 C 1 + v 2 C 2 ) { F L 1 , L 2 [ f a ( x , y ) ] F L 1 , L 2 [ h a ( x , y ) ] F L 1 , L 2 [ f b ( x , y ) ] F L 1 , L 2 [ h b ¯ ( x , y ) ] } + e i ( u 2 C 1 + v 2 C 2 α β ) { F L 1 , L 2 [ f b ( x , y ) ] F L 1 , L 2 [ h b ( x , y ) ] + F L 1 , L 2 [ f b ( x , y ) ] F L 1 , L 2 [ h a ¯ ( x , y ) ] } j

From theorem 4 we draw the conclusion that the convolution of two quaternion signals is to the summation of product of their components, conjugated or odd-even operated, and the product is frequency modulated by chirps. Meanwhile, it must be noted that the orders of i and j in cannot be disordered.

Theorem 5: For any two quaternion signals

f ( x , y ) = f a ( x , y ) + f b ( x , y ) j and h ( x , y ) = h a ( x , y ) + h b ( x , y ) j

where f a ( x , y ) = f r ( x , y ) + i f i ( x , y ) , f b ( x , y ) = f j ( x , y ) + i f k ( x , y ) , h a ( x , y ) = h r ( x , y ) + i h i ( x , y ) , h b ( x , y ) = h j ( x , y ) + i h k ( x , y ) , set

A 1 , 2 = 1 2 π i b 1 b 2 and A 1 , 2 1 , 1 = 1 2 π i ( b 1 ) ( b 2 ) ,

g ( x , y ) = ( f ¯ ¯ h ) ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) 2 π ( e i ( x 2 C 1 + y 2 C 2 ) f ( x , y ) ) ( h ( x , y ) e i ( x 2 C 1 + y 2 C 2 ) )

where B 1 = 1 / b 1 , B 2 = 1 / b 2 , C 1 = a 1 2 b 1 , C 2 = a 2 2 b 2 , then

F L 1 , L 2 { e i ( u 2 C 1 + v 2 C 2 ) f ( x , y ) g ( x , y ) } = A 1 , 2 1 , 1 { ( f a ( x , y ) ) L 1 , L 2 ¯ ¯ ( h a ( x , y ) ) L 1 , L 2 ( f b ( x , y ) ) L 1 , L 2 ¯ ¯ ( h b ¯ ( x , y ) ) L 1 , L 2 } ( u , v ) + A 1 , 2 1 , 1 { ( f a ( x , y ) ) L 1 , L 2 ¯ ¯ ( h b ( x , y ) ) L 1 , L 2 + ( f b ( x , y ) ) L 1 , L 2 ¯ ¯ ( h a ¯ ( x , y ) ) L 1 , L 2 } ( u , v ) j (12)

Proof: Since

e i ( u 2 C 1 + v 2 C 2 ) f ( x , y ) g ( x , y ) = e i ( u 2 C 1 + v 2 C 2 ) { f a ( x , y ) h a ( x , y ) + f a ( x , y ) h b ( x , y ) j + f b ( x , y ) h a ¯ ( x , y ) j f b ( x , y ) h b ¯ ( x , y ) }

From theorem 3, it can be obtained:

F L 1 , L 2 { e i ( u 2 C 1 + v 2 C 2 ) f a ( x , y ) h a ( x , y ) } = A 1 , 2 1 , 1 { ( f a ( x , y ) ) L 1 , L 2 ¯ ¯ ( h a ( x , y ) ) L 1 , L 2 } ( u , v )

F L 1 , L 2 { e i ( u 2 C 1 + v 2 C 2 ) f a ( x , y ) h b ( x , y ) j } = A 1 , 2 1 , 1 { ( f a ( x , y ) ) L 1 , L 2 ¯ ¯ ( h b ( x , y ) ) L 1 , L 2 } ( u , v ) j

F L 1 , L 2 { e i ( u 2 C 1 + v 2 C 2 ) f b ( x , y ) h a ¯ ( x , y ) j } = A 1 , 2 1 , 1 { ( f b ( x , y ) ) L 1 , L 2 ¯ ¯ ( h a ¯ ( x , y ) ) L 1 , L 2 } ( u , v ) j

F L 1 , L 2 { e i ( u 2 C 1 + v 2 C 2 ) f b ( x , y ) h b ( x , y ) } = A 1 , 2 1 , 1 { ( f b ( x , y ) ) L 1 , L 2 ¯ ¯ ( h b ¯ ( x , y ) ) L 1 , L 2 } ( u , v )

From the linear property of Fourier transform:

F L 1 , L 2 { e i ( u 2 C 1 + v 2 C 2 ) f ( x , y ) g ( x , y ) } = A 1 , 2 1 , 1 { ( f a ( x , y ) ) L 1 , L 2 ¯ ¯ ( h a ( x , y ) ) L 1 , L 2 ( f b ( x , y ) ) L 1 , L 2 ¯ ¯ ( h b ¯ ( x , y ) ) L 1 , L 2 } ( u , v ) + A 1 , 2 1 , 1 { ( f a ( x , y ) ) L 1 , L 2 ¯ ¯ ( h a ( x , y ) ) L 1 , L 2 + ( f b ( x , y ) ) L 1 , L 2 ¯ ¯ ( h a ¯ ( x , y ) ) L 1 , L 2 } ( u , v ) j

From theorem 5 we draw the conclusion that, the product, frequency modulated by a chirp, of two quaternion signals is to the summation, amplitude modulated, of their pseudo convolution.

4.2. FRQCR

Headings, or heads, are organizational devices that guide the reader through your paper. There are two types: component heads and text heads.

Theorem 6 is suitable for scalar and complex signals, and theorem 7 is suitable for scalar, complex signals, vector and quaternion signals.

Theorem 6: For two scalar (or complex) signals f ( x , y ) and h ( x , y ) , A 1 , 2 = 1 2 π i b 1 b 2 , C 1 = a 1 2 b 1 , C 2 = a 2 2 b 2 , f ( x , y ) , h ( x , y ) = + + f ( τ , η ) h ( x + τ , y + η ) ¯ d τ d η , and set g ( x , y ) = f ( x , y ) h ( x , y ) = A 1 , 2 e i ( C 1 x 2 + C 2 y 2 ) e i ( C 1 x 2 + C 2 y 2 ) f ( x , y ) , e i ( C 1 x 2 + C 2 y 2 ) h ( x , y ) , then:

f ( x , y ) h ( x , y ) = ( f ( x , y ) ) ¯ ( h ( x , y ) ¯ ) (13)

Proof: the proof is similar with that of FRQCV and is omitted here.

From theorem 6 we draw the conclusion that correlation can be implemented by convolution.

Theorem 7: For any two quaternion signals f ( x , y ) and h ( x , y ) , A 1 , 2 = 1 2 π i b 1 b 2 , C 1 = a 1 2 b 1 , C 2 = a 2 2 b 2 , f ( x , y ) , h ( x , y ) = + + f ( τ , η ) h ( x + τ , y + η ) ¯ d τ d η , and let g ( x , y ) = f ( x , y ) h ( x , y ) = A 1 , 2 e i ( C 1 x 2 + C 2 y 2 ) e i ( C 1 x 2 + C 2 y 2 ) f ( x , y ) , e i ( C 1 x 2 + C 2 y 2 ) h ( x , y ) , “ ” is correlation operator, then

f ( x , y ) h ( x , y ) = ( f a ( x , y ) ) ¯ ( h a ( x , y ) ¯ ) + ( f b ( x , y ) ) ¯ ( h b ( x , y ) ¯ ) A 1 , 2 e i ( x 2 C 1 + y 2 C 2 ) { ( e i ( x 2 C 1 + y 2 C 2 ) f a ( x , y ) ) ( h b ( x , y ) ¯ e i ( x 2 C 1 + y 2 C 2 ) ) } j + A 1 , 2 e i ( x 2 C 1 + y 2 C 2 ) { ( e i ( x 2 C 1 + y 2 C 2 ) f b ( x , y ) ) ( h a ( x , y ) ¯ e i ( x 2 C 1 + y 2 C 2 ) ) } j (14)

Proof: The proof is similar with that of FRQCV and is omitted here.

From theorem 7 we draw the conclusion that the correlation of two quaternion signals is to the summation of convolution of their components, conjugated or odd-even operated. It means that correlation can be implemented by convolution via FFT.

5. Fast Algorithm of QLCT

Fast algorithm of QLCT is the key to engineering use. The following discusses the efficient implementation in great detail through the decomposition of quaternion [24]and the definition of the QLCT. For one quaternion function f ( x , y ) , from definition 1 we have

F i , j L 1 , L 2 { f ( x , y ) } = + + K L 1 , i ( x , u ) f ( x , y ) K L 2 , j ( y , v ) d x d y = 1 2 π b 1 i e i d 1 u 2 2 b 1 + + e i u x b 1 e i x 2 a 1 2 b 1 f ( x , y ) e j y 2 a 2 2 b 2 e j v y b 2 d x d y 1 2 π b 2 j e j d 2 v 2 2 b 2 = G i ( u ) + + e i u x b 1 g ( x , y ) e j v y b 2 d x d y G j (v)

where,

G i ( u ) = 1 2 π b 1 i e i d 1 u 2 2 b 1 , G j ( v ) = 1 2 π b 2 j e j d 2 v 2 2 b 2 ,

g ( x , y ) = e i x 2 a 1 2 b 1 f ( x , y ) e j y 2 a 2 2 b 2 = g r ( x , y ) + i g i ( x , y ) + j g j ( x , y ) + k g k ( x , y )

g r ( x , y ) , g i ( x , y ) are real signals.

Let

W ( u , v ) = + + e i u x b 1 g ( x , y ) e j v y b 2 d x d y (17)

Then,

W ( u , v ) + W ( u , v ) 2 = + + e i u x b 1 g ( x , y ) cos ( v y b 2 ) d x d y (18)

W ( u , v ) W ( u , v ) 2 = + + e i u x b 1 g ( x , y ) sin ( v y b 2 ) d x d y ( i ) (19)

Therefore,

W ( u , v ) + W ( u , v ) 2 + W ( u , v ) W ( u , v ) 2 ( k ) = + + e i u x b 1 g ( x , y ) e j v y b 2 d x d y (20)

Therefore,

F i , j L 1 , L 2 ( u , v ) = G i ( u ) W ( u , v ) ( 1 k ) + W ( u , v ) ( 1 + k ) 2 G j ( v ) (21)

Then the following task is to implement W ( u , v ) .

g ( x , y ) can be expressed as

g ( x , y ) = g r ( x , y ) + i g i ( x , y ) + j g j ( x , y ) + k g k ( x , y ) = g a ( x , y ) + g b ( x , y ) j

where, g a ( x , y ) = g r ( x , y ) + i g i ( x , y ) , g b ( x , y ) = g j ( x , y ) + i g k ( x , y ) .

Therefore,

W ( u , v ) = + + e i u x b 1 e i v y b 2 g a ( x , y ) d x d y + ( + + e i u x b 1 e i v y b 2 g b ( x , y ) d x d y ) j = F { g a ( x , y ) } ( u b 1 , v b 2 ) + F { g b ( x , y ) } ( u b 1 , v b 2 ) j (22)

W ( u , v ) can be Calculated by two 2D FFT and some scaling transform. The steps of calculating QLCT:

1) Calculate g ( x , y ) from f ( x , y ) using (16);

2) Calculate W ( u , v ) from g ( x , y ) using (22) and (17);

3) Calculate G i ( u ) and G j ( v ) using (16);

4) At last Calculate F i , j L 1 , L 2 ( u , v ) using (20) and (21).

For one 2D discrete signal with size M × N, one 2D-DFT needs M N log 2 ( M N ) real number multiplications [25]. To implement W ( u , v ) , we require 2 M N log 2 ( M N ) real number multiplications. Therefore, the complexity of quaternion signal f ( x , y ) is O ( 2 M N log 2 ( M N ) ) . And the complexity of scalar, complex signal and vector is: O ( M N log 2 ( M N ) 2 ) , O ( M N log 2 ( M N ) ) , O ( 3 M N log 2 ( M N ) 2 ) .

For any discrete 2D signal f ( m , n ) ( m [ 1 , M ] , n [ 1 , N ] ), it can be expressed as:

f ( m , n ) = f e e ( m , n ) + f e o ( m , n ) + f o e ( m , n ) + f o o ( m , n )

where:

f e e ( m , n ) = f ( m , n ) + f ( n , N n ) + f ( M m , n ) + f ( M m , N n ) 4

f o e ( m , n ) = f ( m , n ) + f ( n , N n ) f ( M m , n ) f ( M m , N n ) 4

f e o ( m , n ) = f ( m , n ) f ( n , N n ) + f ( M m , n ) f ( M m , N n ) 4

f o o ( m , n ) = f ( m , n ) f ( n , N n ) f ( M m , n ) + f ( M m , N n ) 4

If in the right side of f ( m , n ) = f e e ( m , n ) + f e o ( m , n ) + f o e ( m , n ) + f o o ( m , n ) there is only one term, we call f ( m , n ) symmetric;

If f ( M m , n ) = ± f ( m , n ) , we call f ( m , n ) symmetric about x;

If f ( m , N n ) = ± f ( m , n ) , we call f ( m , n ) symmetric about y;

If any above relationship is not true, we call f ( m , n ) asymmetric.

The symmetry is of great importance to greatly decreasing the calculation complexity of them. Table 1 lists the calculation complexity of different types of signals. It gives the conclusion that the symmetry can decrease the calculation complexity by a few times. Meanwhile, the calculation complexity will increase with the number of components by a few times.

Meanwhile, the calculation complexity of QLCT for different signals is multiplications. Also, the complexity of QCV and QCR for the same type of signals is the same and is much less than calculation in time-domain directly.

Figure 1 shows one intuitive result. The QCR of the quaternion signal f ( x , y ) and kernel h ( x , y ) is calculated. We take different signals (scalar, complex, vector and quaternion) as the convolution kernel h ( x , y ) . The red lines denote the complexity of implementing QCR in time domain directly, and the blue lines denote the complexity of implementing QCR via FFT. For example,

Table 1. The calculation complexity of QLCT for different signals.

Figure 1. The comparison of complexity via FFT and calculation directly.

when the size is 60, there is one nearly-ten-times relationship. Moreover, with the increase of size the gap would become bigger and bigger.

6. Conclusion

One contribution of this paper is that the definition of QLCT is obtained, and its properties are given, and its generalization is proved. The reversibility property disclosed the efficiency of QLCT. The linear property indicated that LCT is linear transform. Another contribution of this paper is that the QCV and QCR of LCT are defined and their relationships and physical interpretation are discovered: the fractional convolution of two quaternion signals is to the summation of product of their components, conjugated or odd-even operated, and the product is frequency modulated by chirps; and the product, frequency modulated by a chirp, of two quaternion signals is to the summation, amplitude modulated, of their pseudo convolution; and the correlation of two quaternion signals is to the summation of convolution of their components, which are conjugated or odd-even operated. The last contribution is that the complexity of QLCT, QCV and QCR are given, and its Fast Algorithm is obtained through implementing them via the product theorem in transformed domain whose complexity is similar to FFT, which is of great importance to engineering use [31][32].

Acknowledgements

This work was fully supported by the NSFCs (61471412, 61771020, 61002052, 61250006).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Zhang, Y. and Xu, G.L. (2018) The Properties and Fast Algorithm of Quaternion Linear Canonical Transform. Journal of Signal and Information Processing, 9, 202-216. https://doi.org/10.4236/jsip.2018.93012

References

  1. 1. Exampleh, V., Ozaktas, M. and Aytur, O. (1995) Fractional Fourier Domains. Signal Process, 46, 119-124. https://doi.org/10.1016/0165-1684(95)00076-P

  2. 2. Tao, R., Deng, B. and Wang, Y. (2009) Theory and Application of the Fractional Fourier Transform. Tsinghua University Press, Beijing.

  3. 3. Moshinsky, M. and Quesne, C. (1971) Linear Canonical Transformations and Their Unitary Representation. Journal of Mathematical Physics, 12, 1772-1783. https://doi.org/10.1063/1.1665805

  4. 4. Xu, G.L., Wang, X.T. and Xu, X.G. (2009) The Logarithmic, Heisenberg’s and Windowed Uncertainty Principles in Fractional Fourier Transform Domains. Signal Processing, 89, 339-343. https://doi.org/10.1016/j.sigpro.2008.09.002

  5. 5. Xu, G.L., Wang, X.T. and Xu, X.G. (2009) Generalized Hilbert Transform and Its Properties in 2D LCT. Signal Process, 89, 1395-1402. https://doi.org/10.1016/j.sigpro.2009.01.009

  6. 6. Mendlovic, D. and Ozaktas, H.M. (1993) Fractional Fourier Transforms and Their Optical Implementation: I. Journal of the Optical Society of America A, 10, 1875-1881.

  7. 7. Ozaktas, H.M. and Mendlovic, D. (1993) Fractional Fourier Transforms and Their Optical Implementation. II. Journal of the Optical Society of America A, 10, 2522-2531. https://doi.org/10.1364/JOSAA.10.002522

  8. 8. Pei, S.C. and Yeh, M.H. (1998) Two Dimensional Discrete Fractional Fourier Transform. Signal Processing, 67, 99-108. https://doi.org/10.1016/S0165-1684(98)00024-3

  9. 9. Aytur, O. and Ozaktas, H.M. (1995) Non-Orthogonal Domains in Phase Space of Quantum Optics and Their Relation to Fractional Fourier Transforms. Optics Communications, 120, 166-170. https://doi.org/10.1016/0030-4018(95)00452-E

  10. 10. Ozaktas, H.M., Zalevsky, Z. and Kutay, M.A. (2000) The Fractional Fourier Transform with Applications in Optics and Signal Processing. Wiley, New York.

  11. 11. Hamilton, W.R. (1864) Elements of Quaternions. Longman, London.

  12. 12. Xu, G.L., Wang, X.T. and Xu, X.G. (2010) On Uncertainty Principle for the Linear Canonical Transform of Complex Signals. IEEE Transactions on Signal Processing, 58, 4916-4918. https://doi.org/10.1109/TSP.2010.2050201

  13. 13. Stark, H. (1971) An Extension of the Hilbert Transform Product Theorem. Proceedings of the IEEE, 59, 1359-1360.

  14. 14. Havlicek, J.P., Havlicek, J.W., Ngao, D., et al. (1998) Skewed 2D Hilbert Transforms and Computed AM-FM Models. Proceedings 1998 International Conference on Image Processing, Chicago, 7 October 1998, 602-606.

  15. 15. Xu, G.L., Wang, X.T. and Xu, X.G. (2008) Extended Hilbert Transform for Multidimensional Signals. 2008 5th International Conference on Visual Information Engineering, Xi’an, 29 July-1 August 2008, 292-297.

  16. 16. Hahn, S.L. (1992) Multidimensional Complex Signals with Single-Orthant Spectra. Proceedings of the IEEE, 80, 1287-1300.

  17. 17. Moxey, C.E., Sangwine, S.J. and Ell, T.A. (2003) Hypercomplex Correlation Techniques for Vector Images. IEEE Transactions on Signal Processing, 51, 1941-1953. https://doi.org/10.1109/TSP.2003.812734

  18. 18. Xu, G.L., Wang, X.T. and Xu, X.G. (2007) Neighborhood Limited Empirical Mode Decomposition and Application in Image Processing. 4th International Conference on Image and Graphics, Sichuan, 22-24 August 2007, 149-154.

  19. 19. Xu, G.L., Wang, X.T. and Xu, X.G. (2008) Fractional Quaternion Fourier Transform, Convolution and Correlation. Signal Processing, 88, 2511-2517. https://doi.org/10.1016/j.sigpro.2008.04.012

  20. 20. Sangwine, S.J. and Ell, T.A. (1999) Hypercomplex Auto- and Cross-Correlation of Color Images. Proceedings of the 1999 International Conference on Image Processing, Kobe, 24-28 October 1999, 319-322.

  21. 21. Ell, T.A. (1993) Quaternion-Fourier Transforms for Analysis of Two-Dimensional Linear Time-Invariant Partial Differential Systems. Proceedings of 32nd IEEE Conference on Decision and Control, San Antonio, 15-17 December 1993, 1830-1841. https://doi.org/10.1109/CDC.1993.325510

  22. 22. Xu, G.L., Wang, X.T. and Xu, X.G. (2009) Improved Bi-Dimensional EMD and Hilbert Spectrum for the Analysis of Textures. Pattern Recognition, 42, 718-734. https://doi.org/10.1016/j.patcog.2008.09.017

  23. 23. Ell, T.A. (1992) Hypercomplex Spectral Transforms. Ph.D. Dissertation, University of Minnesota, Minneapolis.

  24. 24. Pei, S.C., Ding, J.J. and Chang, J.H. (2001) Efficient Implementation of Quaternion Fourier Transform, Convolution, and Correlation by 2-D Complex FFT. IEEE Transactions on Signal Processing, 49, 2783-2797. https://doi.org/10.1109/78.960426

  25. 25. Duhamel, P. (1986) Implementation of Split-Radix FFT Algorithms for Complex, Real and Real-Symmetric Data. IEEE Transactions on Acoustics, Speech, and Signal Processing, 34, 285-295. https://doi.org/10.1109/TASSP.1986.1164811

  26. 26. Almeida, L.B. (1997) Product and Convolution Theorems for the Fractional Fourier Transform. IEEE Signal Processing Letters, 4, 15-17. https://doi.org/10.1109/97.551689

  27. 27. Zayed, A.I. (1998) A Convolution and Product Theorem for the Fractional Fourier Transform. IEEE Signal Processing Letters, 5, 101-103. https://doi.org/10.1109/97.664179

  28. 28. Akay, O. and Boudreaus, G.F. (1998) Linear Fractionally Invariant Systems: Fractional Filtering and Correlation via Fractional Operators. Conference Record of the 31st Asilomar Conference on Signals, Systems and Computers, Vol. 2, Pacific Grove, 2-5 November 1997, 1494-1498.

  29. 29. Deng, B., Tao, R. and Wang, Y. (2006) Convolution Theorems for the Linear Canonical Transform and Their Applications. Science in China F, 49, 592-603. https://doi.org/10.1007/s11432-006-2016-4

  30. 30. Xu, G., Wang, X. and Xu, X. (2012) On Analysis of Bi-Dimensional Component Decomposition via BEMD. Pattern Recognition, 45, 1617-1625. https://doi.org/10.1016/j.patcog.2011.11.004

  31. 31. Xu, G., Zhou, L., Wang, X. and Xu, X. (2017) Assisted Signals Based Mode Decomposition. International Conference on Image, Vision and Computing, Chengdu, 2-4 June 2017, 868-874.

  32. 32. Xu, G., Wang, X., Zhou, L. and Xu, X. (2018) Image Decomposition and Texture Analysis via Combined Bi-Dimensional Bedrosian’s Principles. IET Image Processing, 12, 262-273. https://doi.org/10.1049/iet-ipr.2017.0494