Next Article in Journal
Recomputing Causality Assignments on Lumped Process Models When Adding New Simplification Assumptions
Next Article in Special Issue
Dynamics of Trapped Solitary Waves for the Forced KdV Equation
Previous Article in Journal
Nonclassical Symmetries of a Power Law Harry Dym Equation
Previous Article in Special Issue
Three Classes of Fractional Oscillators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Framework for Circular Multilevel Systems in the Frequency Domain

1
School of Mathematical Sciences, University of Electronic Science and Technology of China, Chengdu 610054, China
2
Engineering School, DEIM, University of Tuscia, 01100 Viterbo, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2018, 10(4), 101; https://doi.org/10.3390/sym10040101
Submission received: 9 February 2018 / Revised: 29 March 2018 / Accepted: 31 March 2018 / Published: 8 April 2018
(This article belongs to the Special Issue Symmetry and Complexity)

Abstract

:
In this paper, we will construct a new multilevel system in the Fourier domain using the harmonic wavelet. The main advantages of harmonic wavelet are that its frequency spectrum is confined exactly to an octave band, and its simple definition just as Haar wavelet. The constructed multilevel system has the circular shape, which forms a partition of the frequency domain by shifting and scaling the basic wavelet functions. To possess the circular shape, a new type of sampling grid, the circular-polar grid (CPG), is defined and also the corresponding modified Fourier transform. The CPG consists of equal space along rays, where different rays are equally angled. The main difference between the classic polar grid and CPG is the even sampling on polar coordinates. Another obvious difference is that the modified Fourier transform has a circular shape in the frequency domain while the polar transform has a square shape. The proposed sampling grid and the new defined Fourier transform constitute a completely Fourier transform system, more importantly, the harmonic wavelet based multilevel system defined on the proposed sampling grid is more suitable for the distribution of general images in the Fourier domain.

1. Introduction

Wavelet multiresolution representations are one of the effective techniques for analyzing signals and images. The wavelet multiresolution analysis (MRA) technology has been widely used in signal and image processing. It was first given by Mallat [1], and the authors study the difference of information between approximation of a signal at the resolution 2 j + 1 and 2 j , by decomposing this signal on a wavelet basis of L 2 ( R ) . The 2D general MRA technique possesses a square shape in the frequency domain [2,3,4]. To design the filter of circular-shape in the Fourier domain, the classical polar Fourier transformation is considered. However, the classical polar Fourier transform retains the same shape as in the space domain, so new approaches are investigated. One way is to redefine the sampling grid in the Fourier domain. In [5], the authors introduce a pseudo-polar Fourier transform that samples the Fourier transform on the pseudo-polar grid, also known as the concentric squares grid. We will give more details in Section 3. In addition, [6] samples on points that are equally spaced on an arbitrary arc of the unit circle, which brings about the Fractional Fourier transform; and, in [7], the sampling is on spirals of the form A W k , with A , W C . Using this type of sampling, the authors develops a computation algorithm for numerically evaluating the z t r a n s f o r m . Our goal is to obtain the sampling grid in a circular shape; therefore, we hope to design a new type of sampling that ensures the sampling points concentrated in a circular region. Then, the sampling grid has a circular shape in the Fourier domain. Inspired by the pseudo-polar Fourier transform in [5], we will also redefine the Fourier transform on circular sampling grid.
In recent years, many kinds of directional wavelets filters have been designed, in order to further efficiently capture the details of signals. The most widely used directional multilevel system includes curvelets [8], contourlets [9] and shearlets [10,11]. What these wavelets have in common is that they have compact support multiscale structure in the space domain. In the Fourier domain, the support of a multilevel system constitutes a high redundant partition. To reduce the redundancy, we consider designing the multilevel system in the frequency domain directly. We also must ensure that the multilevel system constructs the basis of L 2 ( R ) in the space domain.
To design the multilevel system in the Fourier domain, wavelets with compact support in frequency are needed. According to the definition of the harmonic wavelet [12,13,14,15], it is suitable to construct a directional multilevel structure with harmonic wavelets whose Fourier transforms are compact and are constructed from simple functions like Haar wavelets [16] in the space domain. We will review the basic definition and property of harmonic wavelet in Section 2.
In this work, by defining the circular-shape Fourier transform (CFT), we will construct the circular-shape directional multilevel system (CMS) in the Fourier domain due to the compact support of harmonic wavelets [12,13,14,15,17]. The specific structure is totally different from the general Cartesian system. By introducing the CFT, we plan to give a parallel analogy with the general classical Descartes Fourier transform, and the corresponding circular-shape directional multilevel system is constructed naturally, which is suitable for the circular shape of images in the Fourier domain. More details will be given in Section 4.
This paper is organized as follows: Section 2 reviews the basic definition and property of harmonic wavelets. The design of CFT is given in Section 3. Then, in Section 4, the multilevel system in the frequency domain based harmonic wavelet is constructed. The quantitative test measures and test results are displayed in Section 5 and Section 6.

2. Preliminary

The Basic Definition

Harmonic wavelets are complex wavelets defined in the Fourier domain. It is consists of an even function H e ( ω ) (see Figure 1a) as the real part and an odd function H o ( ω ) (see Figure 1b) as the imaginary part, which are defined by
H e ( ω ) = 1 / 4 π , ω [ 4 π , 2 π ) [ 2 π , 4 π ) , 0 , otherwise . H o ( ω ) = i / 4 π , ω [ 4 π , 2 π ) , i / 4 π , ω [ 2 π , 4 π ) , 0 , otherwise .
Combining H e and H o , we get the harmonic function
H ( ω ) = H e ( ω ) + i H o ( ω ) .
From Label (1), we have
H ( ω ) = 1 / 2 π , ω [ 2 π , 4 π ) , 0 , otherwise .
This is shown in Figure 1c.
The corresponding scaling function S is given in the same way, and the even and odd functions are defined as
S e ( ω ) = 1 / 4 π , ω [ 2 π , 2 π ) , 0 , otherwise . S o ( ω ) = i / 4 π , ω [ 2 π , 0 ) , i / 4 π , ω [ 0 , 2 π ) , 0 , otherwise .
so that, from Label (4),
S ( ω ) = S e ( ω ) + i S o ( ω ) .
Therefore, we have
S ( ω ) = 1 / 2 π , ω [ 0 , 2 π ) , 0 , otherwise ,
shown in Figure 2.
Then, the shifting and scaling of basic functions are denoted as S j , ( ω ) and H j , ( ω ) , which are given as
S j , ( ω ) = 1 / 2 j S ( ω / 2 j ) , H j , ( ω ) = 1 / 2 j H ( ω / 2 j ) ,
where j Z is the scaling parameter, and R is the shifting parameter. According to Label (7), the harmonic wavelet system constructs a basis of L 2 ( R ) in the frequency domain; then, for f L 2 ( R ) , we have
f ( ω ) = j = + = + a j , H ( 2 j ω ) ,
f ( ω ) = = + a S ( ω ) + j = 0 + = + a j , H ( 2 j ω ) ,
where
a = + f ( ω ) S ( x ) d x , a j , = + f ( ω ) H ( 2 j x ) d x .

3. The Circular-Shape Fourier Transform (CFT)

This section describes the circular-shape Fourier transform (CFT). We begin with the sampling grid in the Fourier domain, including the Cartesian coordinates (see Figure 3a) for classical Fourier transform and the pseudo-polar grid in [5] (see Figure 3b). This grid samples points of equally spaced long rays but not equally angles. In order to have a circular structure, the sampling grid in concentric circles (see Figure 3c) is designed, which has equally arc and angle in each circle. This type of sampling is consistent with the distribution of images in the frequency domain.

3.1. The Pseudo-Polar Grid

The pseudo-polar grid Ω R is given as
Ω R = Ω R 1 Ω R 2 ,
where
Ω R 1 = { ( 4 k R N , 2 k R ) : | | N / 2 , | k | R N / 2 } , Ω R 2 = { ( 2 k R , 4 k R N ) : | | N / 2 , | k | R N / 2 } ,
with R = 2 the oversampling parameter. The nod of Ω R 1 is on the solid line in Figure 3b and the nod of Ω R 2 is on dotted line.
For an N × N image u, the general discrete Fourier transform u ^ is evaluated on the N × N Cartesian grid in the form
u ^ ( ω x , ω y ) = x , y = N / 2 N / 2 1 u ( x , y ) e 2 π i N ( x ω x + y ω y ) ,
where { ( ω x , ω y ) : ω x , ω y = N / 2 , , N / 2 . } , and
x , y = N / 2 N / 2 1 | u ( x , y ) | 2 = 1 N 2 ω x , ω y = N / 2 N / 2 1 | u ^ ( ω x , ω y ) | 2 .
Analogously, the pseudo-polar Fourier transform is the same as (13) (see [5]), but { ( ω x , ω y ) Ω R . } . According to the Plancherel theorem, (14) can be modified by introducing the weighting function w
x , y = N / 2 N / 2 1 | u ( x , y ) | 2 = ( ω x , ω y ) Ω R N / 2 1 w ( ω x , ω y ) | u ^ ( ω x , ω y ) | 2 .

3.2. The Circular-Polar Grid (CPG)

In this section, the circular-polar grid (CPG) is designed (see Figure 3c), which is defined as
C R = C R 0 C R ;
where C R 0 = { ( 0 , 0 ) } and C R = C R 1 C R 2
C R 1 = { ( r cos ( π m 0 ) , r sin ( π m 0 ) ) : 1 | r | R , | | m 0 2 } , C R 2 = { ( r sin ( π m 0 ) , r cos ( π m 0 ) ) : 1 | r | R , | | m 0 2 } ,
where m 0 is the sampling number in each circle.
As can be seen from Figure 3c, the nod of C R 1 is on a solid line and the nod on a dotted line belongs to C R 2 . In addition, r in (17) serves as the radius and serves as the parameter of angle. m 0 = 16 is the sampling number. In the CPG coordinates, the nod has the following characteristics, for
C R 1 ( ω x , ω y ) = ( r 1 , θ 1 ) , C R 2 ( ω x , ω y ) = ( r 2 , θ 2 ) ,
where
r 1 = k 1 , r 2 = k 2 , θ 1 = 1 π / m 0 ; θ 2 = 2 π / m 0 .
k i = 0 , R ; i = 1 , 2 and i = m 0 / 2 , m 0 / 2 ; i = 1 , 2 . For each fixed angle θ , the samples of the CPG are equally spaced in the radial direction, and, for each fixed radius r, the grid possesses the same angle. Formally,
Δ r 1 ( k 1 + 1 ) k 1 = 1 ; Δ r 2 ( k 2 + 1 ) k 2 = 1 , Δ θ 1 ( 1 + 1 ) π / m 0 1 π / m 0 = π / m 0 , Δ θ 2 ( 2 + 1 ) π / m 0 2 π / m 0 = π / m 0 ,
where r 1 , r 2 and θ 1 , θ 2 are given by (19).
For an N × N image u, the CFT of u ^ on CPG holds
x , y = N / 2 N / 2 1 | u ( x , y ) | 2 = ( ω x , ω y ) C R w c ( ω x , ω y ) | u ^ ( ω x , ω y ) | 2 ,
and
u ^ C R ( ω x , ω y ) = x , y = N / 2 N / 2 1 u ( x , y ) e 2 π i R m 0 + 1 ( x ω x + y ω y ) .
Using operator notation, we denote the refined C F T of an image u as F p , where
( F p u ) ( r , ) u ^ C R ( r , ) ,
with r = R , , R , = m 0 / 2 , , m 0 / 2 . Now, our goal is to choose weight w c , such that w c satisfies (21), and we have
( ω x , ω y ) C R w c ( ω x , ω y ) | u ^ ( ω x , ω y ) | 2 = ( ω x , ω y ) C R w c ( ω x , ω y ) | x , y = N / 2 N / 2 1 u ( x , y ) E ( x , y ) | 2 = ( ω x , ω y ) C R w c ( ω x , ω y ) [ x , y = N / 2 N / 2 1 x , y = N / 2 N / 2 1 u ( x , y ) E ( x , y ) u ( x , y ) E ( x , y ) ¯ ] = ( ω x , ω y ) C R w c ( ω x , ω y ) x , y = N / 2 N / 2 1 | u ( x , y ) | 2 + ( x , y ) ( x , y ) u ( x , y ) u ( x , y ) ¯ [ ( ω x , ω y ) C R w c ( ω x , ω y ) E ( x , y ) E ( x , y ) ¯ ] ,
where E ( x , y ) e 2 π i R m 0 ( x ω x + y ω y ) . Compared with the left of equation (21), the weights w c holds
( ω x , ω y ) C R w c ( ω x , ω y ) e 2 π i R m + 1 ( x ω x + y ω y ) = δ ( x , y ) ;
with N / 2 x , y N / 2 1 .

3.3. The Choice of Weights w c

In the following, we present the basic condition of weights w c , according to (25), which satisfies that
0 = ( ω x , ω y ) C R w c ( ω x , ω y ) [ cos ( 2 π R m 0 + 1 x ω x ) cos ( 2 π R m 0 + 1 y ω y ) sin ( 2 π R m 0 + 1 x ω x ) sin ( 2 π R m 0 + 1 y ω y ) ] ; 0 = ( ω x , ω y ) C R w c ( ω x , ω y ) [ sin ( 2 π R m 0 + 1 x ω x ) cos ( 2 π R m 0 + 1 y ω y ) + cos ( 2 π R m 0 + 1 x ω x ) sin ( 2 π R m 0 + 1 y ω y ) ] .
According to the symmetry of the C F T , the weighting function w c is assumed to satisfy
w c ( ω x , ω y ) = w c ( ω y , ω x ) , ( ω x , ω y ) C R , w c ( ω x , ω y ) = w c ( ω y , ω x ) , ( ω x , ω y ) C R , w c ( ω x , ω y ) = w c ( ω y , ω x ) , ( ω x , ω y ) C R ,
where four equations of (27) describe the ( ω y = ω x ) -symmetry, ( ω y = ω x ) -symmetry and the o r i g i n -symmetry.
In addition,
( ω x , ω y ) C R w c ( ω x , ω y ) = 1 .
To avoid high complexity, we choose the weight w ( ω x , ω y ) in the form:
w ( ω x , ω y ) = w 0 ( ω x , ω y ) ω x , ω y w 0 ( ω x , ω y ) ,
where
w 0 ( ω x , ω y ) = 1 m 0 r , { ( ω x , ω y ) : ω x 2 + ω y 2 = r 2 } ,
with r [ 1 , R ] , and w ( 0 , 0 ) = 1 .

4. The Construction of Multilevel System in Frequency Domain

In this section, we construct a new type multilevel system on C P G in the frequency domain.

4.1. 2D Basic Harmonic Function

First, we define the 2D basic harmonic wavelet functions in the Fourier domain. For deriving convenience, the wavelet function H and scaling function S can be normalized in the form given by Definition 1 .
Definition 1.
The 2D harmonic basic functions are defined as
H ( r , θ ) : H ( 2 π | r | cos ( | θ | ) ) S ( 2 π | r | sin ( | θ | ) ) , S ( r , θ ) : S ( 2 π | r | cos ( | θ | ) ) S ( 2 π | r | sin ( | θ | ) ) ,
with r [ R , R ] , θ [ 0 , π ] .
Then, the support of H ( r , θ ) and S ( r , θ ) are investigated, according to (3) and (6),
H ( 2 π | r | cos ( | θ | ) ) 0 , S ( 2 π | r | sin ( | θ | ) ) 0
hold simultaneously; therefore,
1 | r | cos ( | θ | ) 2 , 0 | r | sin ( | θ | ) 1 .
Thus, the support of H ( r , θ ) is given as
2 | r | 2 , | θ | π / 4 .
Similarly,
supp S ( r , θ ) = { ( r , θ ) : 0 | r | 1 , | θ | π / 4 } .
Next, the 2D scaling and shifting of H ( r , θ ) and S ( r , θ ) are defined.
Definition 2.
The 2D scaling and shifting of harmonic basic functions in the frequency domain are defined as
H j , ( r , θ ) : H ( 2 π 2 j | r | cos ( 2 j | θ | ) ) S ( 2 π 2 j | r | sin ( 2 j | θ | ) ) , S j , ( r , θ ) : S ( 2 π 2 j | r | cos ( 2 j | θ | ) ) S ( 2 π 2 j | r | sin ( 2 j | θ | ) ) , H j , ( r , θ ) : H ( 2 π 2 j | r | sin ( 2 j | θ | ) S ( 2 π 2 j | r | cos ( 2 j | θ | ) ) , S j , ( r , θ ) : S ( 2 π 2 j | r | sin ( 2 j | θ | π 2 ) ) S ( 2 π 2 j | r | sin ( 2 j | θ | ) ) ,
with j , R .

4.2. The Polar Harmonic Multilevel System in the Frequency Domain (PHMS) on CPG

In this section, we give the definition of the polar harmonic multilevel system (PHMS) defined on CPG.
Definition 3.
The 2D PHMS on CPG is defined as
P H M S : H j , j ( r , θ ) S j , j ( r , θ ) H j , j ( r , θ ) S j , j ( r , θ ) ,
where H j , , H j , , S j , and S j , are given in (36). The level parameter j [ log 2 R ] , the shifting parameter is related to j, and we defined the j = 2 j π 4 with | | 2 j , Z .
From | | 2 j , we have 2 ( 2 j + 1 + 1 ) subbands in each level j, in order to reduce the overlap, we choose 2 ( 2 j + 1 ) subbands; then, the PHMS constructs a partition of the Fourier domain. We displayed the PHMS structure in Figure 4.
For a signal or image u, the corresponding PHMS transform P ( u ) in the frequency domain can be defined as
P ( u ) < u ^ C R , P H M S > = j = 0 J j = ( 2 j π 4 ) ( 2 j π 4 ) = 2 j 2 j ( H j , j . u ^ C R + S j , j . u ^ C R + H j , j . u ^ C R + S j , j . u ^ C R ) ,
with j = 2 j π 4 , | | 2 j , Z , where u ^ C R is the C P F T of u, defined in (22). In addition, . is the dot product, and the matrix M 1 . M 2 is defined as
( M 1 . M 2 ) i , j = ( M 1 ) i , j ( M 2 ) i , j .
Theorem 1.
The discrete polar harmonic multilevel system PHMS defined on CPG forms a framelet of L 2 ( R 2 ) .
According to the framelet defined in [18], for the signal U in (38),
U 2 < U , P H M S > 2 c U 2 ,
where c < + is the constant; therefore, PHMS forms a framelet of L 2 ( R 2 ) .
Then two denoising reconstruction tests of PHMS are shown in Figure 5.

5. Quantitative Test Measures

In the following, several performance measures are introduced to test the quality of the P H M S . The quality measure is the Monte Carlo estimate for the different operator norm by generating a sequence of five random images u i , i = 1 , , 5 on C P G for R = 256 , = 8 with standard normally distributed entries.
  • Isometry of C F T :
    (a)
    Closeness to tight: M c l o = max i = 1 , , 5 F p F p u i u i 2 u i 2 ,
    (b)
    Quality of preconditioning. M q u a = λ max ( F p F p ) λ min ( F p w F p ) .
  • Tight Frame Property: The operator norm P P I o p , which is defined as M t i g = max i = 1 , , 5 P P u i u i 2 u i 2 .
  • Robustness:
    (a)
    Thresholding: Let u be the regular sampling of a Gaussian function with mean 0 and variance 512 on [ 257 ] 2 generating an 512 × 512 image. Two types of robustness are considered, for k = 1 , 2 , and M p k = P T p k P u u 2 u 2 .
    T p 1 :
    T p 1 discards 100 ( 1 2 p 1 ) percent of coefficient, with p 1 = [ 2 : 2 : 20 ] .
    T p 2 :
    T p 2 keeps the absolute value of coefficients bigger than m / 2 p 2 with m is the maximal absolute value of all coefficients, where p 2 = [ 0.5 : 0.5 : 5 ] .
    (b)
    Quantization: The quality measure is given as M p = P Q q P u u 2 u 2 , where Q q ( c ) = r o u n d ( c / ( m / 2 q ) ) · ( m / 2 q ) , and q [ 5 : 0.5 : 0.5 ] .

6. Test Results

In this section, the test results of P H M S on C P G for quantitative measure in (1)–(3) are shown.
First the performance with respect to quantitative measures in ( 1 ) , ( 2 ) are presented in Table 1.
From Table 1, the error of tightness M t i g for the P H M S transform is about 0.1, which confirms that the multilevel system is indeed not a tight frame. The main reason is the redundancy of the multilevel structure due to the radius in (36). The quantity M c l o 0.0095 and M q u a 1.934 suggests that the circular-polar Fourier transform provides good properties in terms of isometry, which allows us to employ the conjugate gradient method to compute the inverse of F p . The weight w should also be chosen carefully.
Second, the robustness measurements in ( 3 ) displayed in Table 2:
Table 2 shows the robustness of P H M S . Even discarding 100 ( 1 2 10 ) 99.9 % of the coefficients, the image still can be recovered with error M P 1 = 0.9 × 10 2 . The second row suggests that just the coefficient greater than thresholding value m ( 1 1 / 2 0.001 ) 0.1 % can give a good reconstruction with M p 2 = 0.009 . In the third row, the quantization of robustness M p is displayed.

7. Conclusions

In this work, we developed and implemented a polar harmonic multilevel system on the circular-polar grid based on the multiscale theory, and testified to the performance of the P H M S in four different quantitative measures, which suggests that the P H M S transform is suitable for the circular-shape Fourier transform. Another advantage of the P H M S is that the circular low frequency region in the frequency domain is consistent with the image spectrum distribution, which can process the multilevel structure more effectively.

Acknowledgments

The first author Guomin Sun is very grateful to the China Scholarship Council for funding the first author visiting the University of Tuscia. This work was supported by the National Natural Science Foundation of China under Grant 11271001, Grant 61370147, and Grant 61573085.

Author Contributions

Guomin Sun and Carlo Cattani designed the circular multilevel system and performed the experiments, Jinsong Leng analyzed the results. Guomin Sun and Carlo Cattani wrote the paper. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CPGCircular-polar grid
MRAMultiresolution analysis
CFTCircular-shape Fourier transform
CMSCircular-shape directional multilevel system
PHMSPolar harmonic multilevel system

References

  1. Mallat, S.G. A Theory for Multiresolution Signal Decomposition: The Wavelet Representation. IEEE Trans. Pattern Anal. Mach. Intell. 1989, 11, 674–693. [Google Scholar] [CrossRef]
  2. Leng, J.; Huang, T.; Jing, Y.; Jiang, W. A Study on Conjugate Quadrature Filters. EURASIP J. Adv. Signal Process. 2011, 2011, 2317–2329. [Google Scholar] [CrossRef]
  3. Leng, J.; Cheng, Z.; Huang, T.; Lai, C. Construction and properties of multiwavelet packets with arbitrary scale and the related algorithms of decomposition and reconstruction. Comput. Math. Appl. 2006, 51, 1663–1676. [Google Scholar] [CrossRef]
  4. Leng, J.; Huang, T.; Fu, Y. Construction of bivariate nonseparable compactly supported biorthogonal wavelets. In Proceedings of the International Conference on Machine Learning and Cybernetics IEEE, Kunming, China, 12–15 July 2008. [Google Scholar]
  5. Averbuch, A.; Coifman, R.R.; Donoho, D.L.; Israeli, M.; Shkolnisky, Y. A Framework for Discrete Integral Transformations I-The Pseudopolar Fourier Transform. Siam J. Sci. Comput. 2007, 30, 764–784. [Google Scholar] [CrossRef]
  6. Bailey, D.H.; Swarztrauber, P.N. The Fractional Fourier Transform and Applications. SIAM Rev. 1991, 33, 389–404. [Google Scholar] [CrossRef]
  7. Rabiner, L.R. Chirp Z-transform algorithm. IEEE Trans. Audio Electroacoust. 1969, 17, 86–92. [Google Scholar] [CrossRef]
  8. Candes, E.J.; Donoho, D.L. New tight frames of curvelets and optimal representations of objects with smooth singularities. Commun. Pure Appl. Math. 2004, 57, 219–266. [Google Scholar] [CrossRef]
  9. Do, M.N.; Vetterli, M. The contourlet transform: an efficient directional multiresolution image representation. IEEE Trans. Image Process. 2005, 14, 2091–2106. [Google Scholar] [CrossRef] [PubMed]
  10. Kutyniok, G.; Shahram, M.; Zhuang, X. ShearLab: A Rational Design of a Digital Parabolic Scaling Algorithm. Siam J. Imaging Sci. 2011, 5, 1291–1332. [Google Scholar] [CrossRef]
  11. Kutyniok, G. Shearlets. In Applied and Numerical Harmonic Analysis; Springer: Birkhauser Basel, Berlin, Germany, 2012. [Google Scholar]
  12. Newland, D.E. Harmonic Wavelet Analysis. Proc. R. Soc. Math. Phys. Eng. Sci. 1993, 443, 203–225. [Google Scholar] [CrossRef]
  13. Cattani, C. Harmonic wavelets towards solution of nonlinear PDE. Comput. Math. Appl. 2005, 50, 1191–1210. [Google Scholar] [CrossRef]
  14. Cattani, C. Fractals Based on Harmonic Wavelets. In Proceedings of the International Conference on Computational Science and Its Applications, Seoul, South Korea, 29 June–2 July 2009; Springer: Berlin, Heidelberg, 2009; pp. 729–744. [Google Scholar]
  15. Cattani, C. Harmonic wavelet approximation of random, fractal and high frequency signals. Telecommun. Syst. 2010, 43, 207–217. [Google Scholar] [CrossRef]
  16. Haar, A. Zur Theorie der orthogonalen Funktionensysteme. Math. Ann. 1911, 71, 38–53. [Google Scholar] [CrossRef]
  17. Duan, X.; Leng, J.; Cattani, C.; Li, C. A Shannon-Runge-Kutta-Gill Method for Convection-Diffusion Equations. Math. Probl. Eng. 2013, 46, 532–546. [Google Scholar] [CrossRef]
  18. Cai, J.; Ji, H.; Liu, C.; Shen, Z. Framelet. IEEE Trans. Image Process. 2012, 21, 562–572. [Google Scholar] [PubMed]
Figure 1. The harmonic wavelet function. (a) the even part H e ( ω ) ; (b) the odd part H o ( ω ) ; (c) the harmonic wavelet function H ( ω ) .
Figure 1. The harmonic wavelet function. (a) the even part H e ( ω ) ; (b) the odd part H o ( ω ) ; (c) the harmonic wavelet function H ( ω ) .
Symmetry 10 00101 g001
Figure 2. The harmonic scaling function. (a) the even part S e ( ω ) ; (b) the odd part S o ( ω ) ; (c) the scaling function S ( ω ) .
Figure 2. The harmonic scaling function. (a) the even part S e ( ω ) ; (b) the odd part S o ( ω ) ; (c) the scaling function S ( ω ) .
Symmetry 10 00101 g002
Figure 3. Three different grids.
Figure 3. Three different grids.
Symmetry 10 00101 g003
Figure 4. The polar harmonic multilevel system (PHMS) in the Fourier domain j 2 , j Z and 2 j j < 2 j .
Figure 4. The polar harmonic multilevel system (PHMS) in the Fourier domain j 2 , j Z and 2 j j < 2 j .
Symmetry 10 00101 g004
Figure 5. Recovery results by PHMS with scale j = 4 , the random noise level σ = 20 .
Figure 5. Recovery results by PHMS with scale j = 4 , the random noise level σ = 20 .
Symmetry 10 00101 g005
Table 1. Results of Labels (1)–(2).
Table 1. Results of Labels (1)–(2).
M clo M qua M tig
0.009460921.93426270.1036731
Table 2. Results of Label ( 3 ) .
Table 2. Results of Label ( 3 ) .
M p 1 3.6 × 10 6 1.7 × 10 5 5.9 × 10 3 2.1 × 10 2 0.9 × 10 2
M p 2 0.0090.0630.1030.1720.197
M p 0.0510.0650.0830.1260.143

Share and Cite

MDPI and ACS Style

Sun, G.; Leng, J.; Cattani, C. A Framework for Circular Multilevel Systems in the Frequency Domain. Symmetry 2018, 10, 101. https://doi.org/10.3390/sym10040101

AMA Style

Sun G, Leng J, Cattani C. A Framework for Circular Multilevel Systems in the Frequency Domain. Symmetry. 2018; 10(4):101. https://doi.org/10.3390/sym10040101

Chicago/Turabian Style

Sun, Guomin, Jinsong Leng, and Carlo Cattani. 2018. "A Framework for Circular Multilevel Systems in the Frequency Domain" Symmetry 10, no. 4: 101. https://doi.org/10.3390/sym10040101

APA Style

Sun, G., Leng, J., & Cattani, C. (2018). A Framework for Circular Multilevel Systems in the Frequency Domain. Symmetry, 10(4), 101. https://doi.org/10.3390/sym10040101

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop