Accurate Computation of Fractional-Order ExponentialMoments

Qilu University of Technology (Shandong Academy of Sciences), Shandong Computer Science Center (National Supercomputer Center in Jinan), Shandong Provincial Key Laboratory of Computer Networks, Jinan, China School of Computer Science and Technology (School of Cyber Security), Qilu University of Technology (Shandong Academy of Sciences), Jinan 250353, China Shandong Provincial Key Laboratory for Distributed Computer Software Novel Technology, Jinan 250358, China


Introduction
e research on image retrieval has been started since the middle and late last century. At that time, it was mainly text-based image retrieval technology, and the description of image features included text-related information. Later, the image retrieval technology was extended to cloud retrieval, i.e., content-based image retrieval technology, which analyzes the color, texture, and layout of images. Shape is the basic image feature used in content-based image retrieval systems. e image moments are robust and effective shape features. e image moments are excellent image descriptors, and they have strong geometric invariance and global feature description ability. erefore, image moments have also been widely used in the field of image processing [1], including object recognition, image reconstruction, image encryption, and information hiding [2]. e existing image moments are mainly divided into nonorthogonal moments and orthogonal moments. Nonorthogonal moments such as Hu moments [3] and complex moments [4] project images to a set of the nonorthogonal functional polynomial. e translation and rotation of the image and the scale change invariant can be constructed based on nonorthogonal moments. While because their basis function does not have the orthogonal relationship and nonorthogonal moments have large redundancy, it is difficult to realize the image reconstruction, and they are more sensitive to the image noise. Orthogonal moments are projection coefficients which project the image to a set of orthogonal polynomials. An image can be reconstructed based on orthogonal moments. e research work shows that it has high robustness to the image noise, image blur, and other related operations [5]. Orthogonal moments are divided into discrete orthogonal moments and continuous orthogonal moments. e continuous orthogonal moments use the continuous function as the basis function, and it has the rotation, scaling, and the translation invariance, which have been greatly developed in recent years, including Legendre moments (LMs) [6], Zernike moments (ZMs) [7], pseudo-Zernike moments (PZMs) [7], orthogonal Fourier-Mellin moments (OFMMs) [8], Chebyshev-Fourier moments (CHFMs) [9], radial harmonic Fourier moments (RHFMs) [10], Bessel-Fourier moments (BFMs) [11], polar harmonic transforms (PHTs) [12], and exponential moments (EMs) [13], Among them, EMs have good antinoise performances and less information redundancy, and their basis functions have the simple form, low computational complexity, and good image description performance [14]. However, EMs have various errors and numerical instabilities at high orders, which affects the accuracy of EMs [15]. e ubiquitous errors have a very negative impact on the image analysis and reconstruction [16] so that when the order of EMs reaches a critical value, the reconstruction errors are too large to be imaged [17]. e promotion of fractional-order exponential moments (FrEMs) afterward compensated for the numerical instability of EMs effectively [18] and improved the reconstruction and antinoise performance of EMs. In the study of fractional moments, scholars first define the fractional parameter of t(t > 0) and then use r t to replace r in the radial basis function of the orthogonal moment. e radial basis function is further modified to maintain the orthogonality of the moment [19]. e orthogonal moment promoted to the fractional order can adjust the gradient of the radial basis function by assigning different values to the fractional variable of t to further alleviate the problem of numerical instability [20]. e existing fractional moments include fractional-order Legendre-Fourier moments (FrOLFMs) [21], orthogonal fractional-order Fourier-Mellin moments (FrOFMMs) [22], fractional-order Zernike moments (FrZMs) [23], fractionalorder polar harmonic transforms (FrPHTs) [24], fractionalorder orthogonal Chebyshev-Fourier moments (FrCFMs) [25], and fractional-order radial harmonic Fourier moments (FrRHFMs) [26].
Although FrEMs have the excellent image description ability, various errors generated by the traditional calculation method still affect the accuracy of FrEMs. And the calculation accuracy restricts the development and application of continuous orthogonal moments in the fields of pattern recognition and image processing. Among the various errors, the numerical integral error is especially prominent. Since the digital image is stored in the form of Cartesian coordinates in computers and other devices, in the Cartesian coordinate system, the direct calculation of the continuous orthogonal moments of the image cannot obtain the integral value of the polynomial correctly, and it can only be replaced with the estimated value. e numerical integration error occurs during this process [27]. e numerical integration error is more distinct when the high order moment is calculated. To solve the problem of the numerical integration error of the image moment, Liao and Pawlak propose a method based on the numerical integration to reduce the numerical integration error. eir method is to use a unit disk with the radius of r � ����������������� � (1 − 1)/(N − 0.0001) [28]. e reduction of the radius is to ensure that the sampling points used in the numerical integration do not cross the boundary of the unit disk to avoid the radial basis function from becoming unbounded. Because the disk area with a reduced radius is further affected by the reduced radius, this method will cause the geometric error. erefore, they conclude that the geometric error and the numerical integration error cannot be reduced at the same time. Later, Singh et al. proposed a technique that can simultaneously reduce the geometric error and the numerical integration error based on the Gaussian numerical integration (GNI) [29].
is paper uses GNI to propose an accurate calculation method of FrEMs based on this idea.
is method provides very accurate FrEMs and reduces the reconstruction error.
From the above description of image moments, we have summarized two problems: (1) in the traditional algorithm, the calculation of image moments mainly uses zeroth order approximation method, which will produce numerical integration errors and affect the calculation accuracy of moments. (2) e numerical instability of continuous orthogonal moments is common at high order, which affects the accuracy of continuous orthogonal moments. e goal of this paper is to take EMs as an example to solve the above two problems of EMs. e experiments prove that the accuracy of EMs is improved after using the new method. e main innovations of this paper are as follows: (1) an accurate calculation method of FrEMs is proposed, and the nature and comparison of the GNI method and the traditional calculation method are analyzed in depth; (2) the experimental result shows that FrEMs using the GNI accurate calculation method have stronger image reconstruction performance and the antinoising attack performances than FrEMs using the traditional method. e rest of this paper is described as follows: in Section 2, we introduce the construction process of FrEMs in detail; Section 3 mainly introduces the traditional calculation methods and GNI methods of FrEMs; Section 4 conducts detailed experiments and discussions on image reconstruction, antinoising attack, antirotation attack, antiscaling attack, antifiltering attack, and anti-JPEG compression attack performance; and Section 5 summarizes the full text.

Definition of EMs.
EMs are the mapping of the image on the basis function. e basis function of EMs is mainly composed of the radial basis function and the angular Fourier factor. e definition of the radial basis function of EMs is as follows [30]: where n is the order and the value range is − ∞ < n < ∞. e range of r is 0 ≤ r ≤ 1. e definition expression of EMs is where m is the repetition and the value range is − ∞ < m < ∞.
[·] * is the conjugate of a complex number. f(r, θ) is the image function in polar coordinates. exp(jmθ) is the angular Fourier factor, and θ is the polar angle with the value range of 0 ≤ θ ≤ 2π.

Security and Communication Networks
EMs have good reconstruction performance, and they can use a limited number of EMs to reconstruct the original image. e reconstruction formula is as follows: (3)

Definition of FrEMs.
In order to improve the performance of EMs, this section we extend EMs to fractional order and construct FrEMs. e radial basis function of FrEMs is defined as follows [31]: where the fractional parameter t > 0, and the basis function of FrEMs is defined as follows: e definition of FrEMs is It can be known from formulas (1) and (4) that when t � 1, the radial basis functions of FrEMs will be those of EMs; therefore, EMs can be deemed as a special form of FrEMs. e radial basis function of EMs is orthogonal within the range of 0 ≤ r ≤ 1: where δ is the Kronecker delta. From the properties of angular Fourier factor and radial basis function, it can be known that the basis function H nm (r, θ) of EMs is orthogonal in the unit circle: From the definition of radial basis function FrA (t) n (r) of FrEMs expressed in formula (4), it can be known that and then, e basis function of FrEMs satisfies the following orthogonal relationship [32]: FrEMs have very strong image reconstruction ability, and the reconstruction formula is as follows:

Analysis of Radial Basis Function.
In this section, we analyze the influence of the selection of the fractional parameter t on the radial basis function of the fractional exponential moment. Figure 1 shows that when t is taken as 1, 1.3, 1.6, and 1.9, respectively, and the order is 30, the radial basis function changes from 0 ≤ r ≤ 1. From Section 2.2, we can see that when the fractional parameter t � 1, the fractional radial basis function is equivalent to the traditional radial basis function. As can be seen from Figure 1(a), when the fractional parameter t � 1, the traditional radial basis function has a larger variation range around r � 0; thus, the variation rate of the radial basis function is larger around r � 0, which leads to numerical instability and large errors. From Figure 1, we can see that the rate of change in the radial basis function gradually becomes moderate with the continuous increase in the fractional parameter t. erefore, we can alleviate the numerical instability of the exponential moment by adjusting the fractional parameter t. However, different fractional parameters t also lead to different calculation emphasis areas. erefore, the specific application of fractional exponential moments should also be considered when selecting fractional parameter t.

Traditional Method.
When the image moment of a digital image is calculated by using the computer simulation, the expression of the integral form should firstly be discretized and the integral should be converted to a sum [8].
e discrete integration first needs to discretize the integration area into small areas. In these small areas, the center point is served as the sample point of the function value of the integrand, and then, the area of each small area is multiplied by the integrand value on the sample point. e products for all the small areas are summed, and the result is the approximate integral value [33]. Since FrEMs are defined in the polar coordinate system, while the image is defined in the rectangular coordinate system, the traditional calculation method needs to first convert FrEMs to the rectangular coordinate system, and then, FrEMs are calculated in the rectangular coordinate system. e polar coordinates (r, θ) and the rectangular coordinates (x, y) are converted first here, and the conversion formula is as follows: x � r cos θ, Security and Communication Networks e infinitesimal relationship between rectangular coordinates and polar coordinates is e definition of FrEMs in the rectangular coordinate system is obtained as follows: e range of the integral change is 0 ≤ x 2 + y 2 ≤ 1, so the image needs to be mapped in the unit circle when FrEMs are calculated in the rectangular coordinate system. Since FrEMs calculated by the circumcircle mapping method do not have rotation invariance, this paper uses a calculation method based on the inscribed circle, as shown in Figure 2.
e formula for mapping the inscribed circle portion of a grayscale image with the size of N × N into the unit circle is as follows: e above mapping relationship is shown in Figure 2(b). e image center is mapped to the center of the unit circle. (x q , y p ) represents the center of the small image region of where Δx � Δy � (2/N). e discrete summation form of FrEMs can be obtained as follows: where r p,q � ������ � x 2 q + y 2 p , θ p,q � arctan(y p /x q ).

GNI Method.
e traditional calculation method has the numerical integration problem [34]. In view of this defect, an accurate calculation method of FrEMs is proposed in this section by using GNI.
For the one-dimensional function of f(x), the integral over the interval [a, b] can be expressed as and then, so − 1 ≤ s ≤ 1, and it can be obtained that where w k and s k are the weight and position of the image sampling point, respectively, o is the order of GNI, and the above formula can be transformed into the following form: Similarly, for the two-dimensional function of f(x, y), its expression of the double GNI in the integration area can be expressed as Now, we use the double GNI method to precisely calculate FrEMs. For formula (13), it can be obtained as follows: where e constraint given in formula (21) is an improvement over the constraint x 2 q + y 2 p ≤ 1 used in the zeroth-order approximation for inscribing circular disk. is constraint also allows those grids to take part in computation whose centers fall outside the circle.

Experiment and Result Analysis
In this section, image reconstruction, antinoising attack, antirotation attack, antiscaling attack, antifiltering attack, and anti-JPEG compression attack performance of the accurate calculation method of FrEMs were tested via the experiment. irty grayscale images with size of 128 × 128 were used as the image library. e images shown in Figure 3 are ten images selected randomly from the image library. For convenience, TFrEMs and GFrEMs were used in the Security and Communication Networks 5 following to refer to the traditional method and GNI method, respectively; GFrPHFMs, GFrRHFMs, GFrPCETs, GFrPCTs, and GFrPSTs, respectively, represent FrPHFMs, FrRHFMs, FrPCETs, FrPCTs, and FrPSTs calculated by the GNI method. Here, we chose the order of GNI as o � 5.

Rotation Angle Estimation. Let the rotation angle of the
and then, Real(ϕ) � cos(mφ), where Real(ϕ) referred to the real part of ϕ and Imag(ϕ) referred to the imaginary part of ϕ. According to the above formula, the rotation angle φ could be estimated by the inverse trigonometric function. GFrEMs of any order could be used to estimate the rotation angle [35]. From Section 2.3, we know that the selection of t will lead to different focus areas. After a lot of experiments, we know that when t � 1.5, the angle estimation experiment result of GFrEMs is the best. Here, we chose the maximum moment order for each estimate as N max � 1, 2, . . . , 20 and the fractional parameter t � 1.5. And for each maximum moment order N max , we selected 2N max moments with all the repetition of m � 1, respectively, to estimate to get 2N max estimated angles. Finally, the average values of these 2N max estimated angles were taken as the corresponding final result of the maximum moment order N max . e experiment used a Lena grayscale image with a size of 128 × 128. Denote φ * by the estimated angle and mean relative error (MRE) as the measurement standard, and the experimental result was as follows. Figure 4 shows the estimated rotation angle after the original images were rotated 30°and 60°, respectively. As could be seen from Figure 4, the angle estimated by GFrEMs is relatively accurate at low maximum moment order, and some deviations occur as the maximum moment order increases, but the MRE can still be kept small. It could be seen from Table 1 that when 30°was rotated, the MRE of the estimated rotation angles of the real and imaginary parts of ϕ used was 0.0688 and 0.0425, respectively. When 60°was rotated, the MRE of the estimated rotation angles of the real and imaginary parts of ϕ used was 0.0247 and 0.0399, respectively, which verified that the rotation angle estimation using GFrEMs was relatively accurate.

Rotation Invariance.
e rotation invariance of GFrEMs was tested in this section. e Lena grayscale image with size of 128 × 128 was rotated by 5°, 15°, 25°, 35°, and 45°, respectively. e GFrEMs amplitude of the original image and the GFrEMs amplitude of the rotated image were compared. MRE was used to represent the change rate of the GFrEMs amplitude of the rotated image relative to the original image. e selection of the fractional parameter t will affect the zero distribution of the radial basis function and further affect the rate of change in the radial basis function. After lots of experiments, it was known that the MRE was the smallest when the fractional parameter t � 1.5, so the fractional parameter t � 1.5 was selected here. Figure 5 shows the experimental image after rotating at different angles, and the experimental results obtained are shown in Table 2.
It could be seen from the results that the value of MRE was less than 0.02 under different rotation angles, which indicated that the FrQEMs amplitude after the image was rotated was approximately the same as that of the original image, which verified the rotation invariance of GFrEMs.

Scaling Invariance.
GFrEMs were calculated for a set of scaled images below. In this experiment, the Lena grayscale images with a size of 128 × 128 were scaled by 0.75, 1.25, 1.5, 1.75, and 2 times, and their GFrEMs amplitudes were calculated, respectively, to be compared with GFrEMs amplitudes values of the original images. e fractional parameter t � 1.5 was selected here. Figure 6 is the experimental image after it was scaled different times. e experimental result obtained is shown in Table 3.
In the experiment, the original image was scaled different times, and the moment value of the scaled image was calculated. As could be seen from Figure 6, the original image was blurred to different degrees after being scaled to different degrees. From the above experimental data, it could be seen that the amplitude of the same GFrEMs of each scaled image was approximately equal, which verified the scaling invariance of GFrEMs.

Filtering Attack.
e filtering attack blurred the edges of the image [36], including the median filtering, Gaussian filtering, and the average filtering. Here, we selected the fractional parameter t � 1.5 and added the filtering attacks of (3 × 3) and (5 × 5) median filtering, Gaussian filtering, and average filtering to the original image, respectively. e Lena grayscale image with a size of 128 × 128 was adopted in the experiment. e images after filtering attack are shown in Figure 7, and the experimental result obtained is shown in Table 4.
As could be seen from Figure 7, filtering attacks did blur the edges of images, and different filtering attacks had different effects. As could be seen from Table 4, as the filtering quality increased, the MRE of the GFrEMs amplitude of the filtered attack image also gradually increased, which showed that the filtering attack did affect the image quality. However, the MRE of the GFrEMs amplitude of the image after (3 × 3) and (5 × 5) Gaussian filtering attack still remained below 0.05, which indicated that GFrEMs could well resist the Gaussian filtering. For median filtering and average filtering, MRE could be kept below 0.1, indicating that GFrEMs had a certain degree of resistance to median filtering and average filtering.     [37]. e purpose of the compression was to compress the amount of data and improve the effectiveness. However, pixels of the image were lost in this process. Here, we selected the fractional parameter t � 1.5 and performed the JPEG compression attack on the original image with a quality factor of 10, 20, . . ., 90 and then compared the amplitude of the GFrEMs of the attacked image with the amplitude of the GFrEMs of the original image. Figure 8 shows the Lena image after the JPEG compression attack with different quality factors, and the experiment result is shown in Table 5.
As could be seen from Figure 8, after JPEG compression attacks with different quality factors are added to the original image, the image quality deteriorates to different degrees. As could be seen from the above Table 5, as the quality factor became larger, the MRE also became smaller and smaller. When the quality factor reached 30, the MRE of the GFrEMs amplitude of the image after the JPEG compression attack was kept below 0.02 compared with the original image, which indicated that GFrEMs had strong resistance to JPEG compression attacks.

Image Reconstruction.
e image reconstruction performance was an important feature of the image orthogonal moment, which reflected the accuracy of the image moment.
For the image f(x, y) with a size of M × N and its reconstructed image f(x, y), the mean square error was used in this paper to measure the reconstruction error [38]: 4.6.1. Experiment 1. GFrEMs have good image reconstruction ability. GFrEMs, TFrEMs, GFrPHFMs, GFrRHFMs, GFrPCETs, GFrPCTs, and GFrPSTs were compared in this section. e Lena grayscale image with a size of 128 × 128 was used in the experiment. e maximum moment order N max � 10, 20, . . . , 50, and the fractional parameter t � 1.9. e experiment result is shown in Table 6. As can be seen from Table 6, the images reconstructed by TFrEMs and GFrEMs kept small errors at low orders, but when the maximum moment order reached to a certain value, the images reconstructed by TFrEMs had a large error; both the edge region and the center region of the image reconstructed by TFrEMs are deteriorated. And the image errors reconstructed by GFrEMs were smaller and smaller.
ere was no deterioration in the edge area and center area of the image reconstructed by GFrEMs, and the image quality got better and better with the increase in the maximum moment order. For a more intuitive explanation, the line chart of the reconstruction error is shown in Figure 9. It can be seen from the line chart that the image reconstructed by GFrEMs always maintained a small error, which indicated that the GNI method further improved the accuracy of FrEMs. From Figure 9, we can see that when t � 1.9, the mean square error of the reconstructed images of GFrEMs and GFrRHFMs was similar, and both GFrEMs and GFrRHFMs could maintain great reconstruction effect. GFrPCETs and GFrPCTs had similar reconstruction effects and could keep the mean square error small, but the overall reconstruction effect was worse than GFrEMs. However, GFrPHFMs and GFrPSTs had the worst reconstruction effect compared with other fractional moments calculated by the GNI method. On the whole, the reconstruction error of GFrEMs is always smaller than other fractional-order moments calculated by the GNI method, which once again verifies the reconstruction performance of GFrEMs.

Experiment 2.
Since the noise would seriously affect the reconstruction performance of images [39], this experiment tested the comparison of the image reconstruction     errors between GFrEMs, TFrEMs, GFrPHFMs, GFrRHFMs, GFrPCETs, GFrPCTs, and GFrPSTs after adding the salt and pepper noise. is experiment used the average value of thirty images. Maximum moment order N max � 10, 20, . . . , 50, and fractional parameter t � 1.9 was selected. e salt and pepper noise σ � 0.01 was added separately. After the salt and pepper noise was added, the image reconstruction result is shown in Table 7, and the error line chart is shown in Figure 10.
It can be seen from Table 7 that the overall effect of the image reconstruction became worse with the addition of noise, which indicated that the noise did affect the image reconstruction. However, from Figure 10, the image reconstruction effect of GFrEMs was always better than TFrEMs and other fractional-order moments, which indicated that GFrEMs had further improved the antinoise performance of FrEMs.
is verified again that the GNI  method was superior to the traditional method in the image reconstruction performance.

Application of GFrEMs in Medical Images.
In this section, the image reconstruction, antinoising attack, antirotation attack, antiscaling attack, antifiltering attack, and anti-JPEG compression attack performance of GFrEMs applied to medical images were tested via the experiment [40]. Seventy grayscale images with a size of 128 × 128 were used as the image library. e images shown in Figure 11 are ten images selected randomly from the image library.

Experiment 1.
In this section, GFrEMs is applied to the reconstruction of 128 × 128 grayscale medical images, and the salt and pepper noise σ � 0.01, 0.02, 0.03 were added separately to test its antinoising attack performance. e maximum moment order N max � 10, 20, . . . , 50, and the fractional parameter t � 1.9. e experiment result is shown in Table 8.
After adding salt and pepper noise, the reconstruction effect of the image becomes worse and the relative error also increases, which proves once again that salt and pepper noise can seriously affect the reconstruction performance of the image moments. From Table 8, it can be seen that the GFrEMs applied to the reconstruction of medical images can still maintain a small error, which once again verifies the reconstruction performance and antinoising attack performance of GFrEMs.

Experiment 2.
In this section, we test the antigeometric attack capability of GFrEMs applied to medical images. e GFrEMs amplitude of the original image and the GFrEMs amplitude of the image being attacked were     Table 9.
As can be seen from Table 9, as the rotation angle of the original image increased, the value of MRE also increased, but the MRE could always be kept within 0.01. And as the quality factor became larger, the MRE also became smaller and smaller, and the MRE could always be kept within 0.04. However, when scaling the original image, we found that too large or too small scaling will lead to larger errors. When the original image was scaled by 0.05 times, the MRE reached 0.0405, which showed that scaling attacks do affect the image quality. As can be seen from Table 9, GFrEMs had strong resistance to Gaussian filtering attacks and certain resistance to average filtering attacks. However, GFrEMs had relatively weak resistance to median attacks, but it could still keep small errors. On the whole, GFrEMs could still keep small errors after various geometric attacks, which proved that GFrEMs applied to medical images have strong robustness.

Conclusion
EMs have good image description abilities, but they have numerical instability problems, which seriously affects the accuracy of EMs. And the traditional calculation method of EMs will generate the numerical integration errors, which restricts the development and application of EMs in the field of pattern recognition and image processing. e main advantages of this paper are manifested as follows: (1) the promotion of FrEMs effectively solves the numerical instability problem and improves the antinoising attack performance and reconstruction performance of EMs to a certain degree. (2) e accurate calculation method based on GNI proposed in this paper effectively reduces the numerical integration error. e experiments show that FrEMs which use this calculation method have a very strong performance in the image reconstruction, rotation angle estimation,  geometric invariance, and resistance to filtering attacks and JPEG compression attacks. Although the proposed method has the above advantages, it takes a long time to calculate. Because the core of the proposed scheme is GFrEMs, for a N × N sized image, we know from (24) that the number of multiplications in the computation of one GFrEMs is o × o × N × N. Since the GFrEMs with maximum moment order N max contains (2N max + 1) × (2N max + 1) moments, the number of multiplications of the GFrEMs with maximum moment order N max is O(o × o × N × N × (2N max + 1) × (2N max + 1)). Hence, the computational complexity of the proposed scheme is O(o 2 N 2 N 2 max ), which is very high. us, our future work is to propose a fast calculation method for this method.

Data Availability
e data used to support the findings of this study are available at https://download.csdn.net/download/weixin_40394701/ 10228207.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this article.