Invariant Image Representation Using Novel Fractional-Order Polar Harmonic Fourier Moments

Continuous orthogonal moments, for which continuous functions are used as kernel functions, are invariant to rotation and scaling, and they have been greatly developed over the recent years. Among continuous orthogonal moments, polar harmonic Fourier moments (PHFMs) have superior performance and strong image description ability. In order to improve the performance of PHFMs in noise resistance and image reconstruction, PHFMs, which can only take integer numbers, are extended to fractional-order polar harmonic Fourier moments (FrPHFMs) in this paper. Firstly, the radial polynomials of integer-order PHFMs are modified to obtain fractional-order radial polynomials, and FrPHFMs are constructed based on the fractional-order radial polynomials; subsequently, the strong reconstruction ability, orthogonality, and geometric invariance of the proposed FrPHFMs are proven; and, finally, the performance of the proposed FrPHFMs is compared with that of integer-order PHFMs, fractional-order radial harmonic Fourier moments (FrRHFMs), fractional-order polar harmonic transforms (FrPHTs), and fractional-order Zernike moments (FrZMs). The experimental results show that the FrPHFMs constructed in this paper are superior to integer-order PHFMs and other fractional-order continuous orthogonal moments in terms of performance in image reconstruction and object recognition, as well as that the proposed FrPHFMs have strong image description ability and good stability.


Introduction
The rapid development of information and network technologies has brought about great changes to human life and production. As more and more digital images are conveniently transmitted and downloaded online, negative impacts also arise from such trends. Each digital image can be accessed freely, and such free access provides the opportunity for launching various attacks on digital images, such as rotation, translation, and scaling, making the application of digital watermarking [1,2] and pattern recognition [3] more difficult. To address these problems, researchers have begun to look for a kind of feature vectors that can represent the objective information contained in images, so as to figure out how to use invariant features to describe images and use a very small number of datasets to represent more image information. Image moments are a type of highly concentrated image features, which serve as a powerful tool to characterize images, and are invariant to rotation, translation, and scaling. Image moments have been widely used in various of their radial polynomials. (2) The relationship between the changes in radial polynomials and the reconstructed images is identified by analyzing the rate of change of radial polynomials, and it is experimentally verified that the FrPHFMs constructed using the proposed algorithm have good performance in image reconstruction and noise resistance. (3) The constructed FrPHFMs are used for object recognition and compared with integer-order PHFMs and other fractional-order continuous orthogonal moments from the perspective of performance. The results of comparison show that the proposed FrPHFMs have better performance in image reconstruction and object recognition.
Other sections of this paper are organized as follows: Section 2 introduces the Fr-PHFMs construction process in detail and analyzes the geometric invariance of FrPHFMs; Section 3 mainly analyzes the properties of FrPHFMs from two perspectives, namely the changes in, and the rate of change of, their radial polynomials; Section 4 describes in detail the experiments and discussions with respect to image reconstruction, geometric invariance, and object recognition; and Section 5 draws a conclusion of this study.

FrPHFMs
In this section, the construction and properties of FrPHFMs are described in detail. Firstly, the traditional integer-order PHFMs are introduced, then the definition of FrPHFMs is given, and, finally, the geometric invariance of FrPHFMs is discussed.

Definition of Integer-Order PHFMs
The integer-order PHFMs with order of n ≥ 0 and repetition of |m| ≥ 0 of image f (r, θ) in a polar coordinate system [31] is defined as: where [·] is the conjugate of a complex number, and basis function H nm (r, θ) is composed of radial polynomial T n (r) and angular Fourier factor exp(jmθ): H nm (r, θ) = T n (r) exp(jmθ) (2) where radial polynomial T n (r) is T n (r) =    1/ √ 2 while n is 0 sin((n + 1)πr 2 ) while n is odd cos(nπr 2 ) while n is even (3) T n (r) is orthogonal within the range of 0 ≤ r ≤ 1: From the property of angular Fourier factor exp(jmθ) and the formula above, it can be known that basis function H nm (r, θ) is orthogonal in the unit circle [32]: 2π 0 1 0 H nm (r, θ)H kl (r, θ)rdrdθ = π 2 δ nk δ ml (5) where 0 ≤ r ≤ 1, 0 ≤ θ ≤ 2π, and δ is the Kronecker delta.
According to the theory of complete system of orthogonal functions, original image function f (r, θ) can be approximately reconstructed using a finite number of PHFMs. Given the PHFMs with the maximum order of n max and the maximum repetition of m max , the formula for approximate reconstruction of the original image is

Definition of FrPHFMs
In this paper, integer-order PHFMs are extended to FrPHFMs. The key to the construction of FrPHFMs is to construct the fractional-order radial polynomials; therefore, the extension of orthogonal radial polynomial T n (r) is considered. Letting t > 0, U (t) n (r) be defined as the radial polynomial below: then the basis function of the FrPHFMs is: n (r) complies with the following orthogonal relationship within the range of 0 ≤ r ≤ 1: From the properties of the angular Fourier factor and radial polynomials, it can be known that basis function W (t) nm (r, θ) is orthogonal in the unit circle and complies with the following orthogonal relationship: Given t > 0, the FrPHFMs with an order of n ≥ 0 and a repetition of |m| ≥ 0 is defined as: From the formula above, it can be seen that, when t = 1, the FrPHFMs is an integerorder PHFMs, and thus is an extension to the integer-order PHFMs.
Given the FrPHFMs with the max moment order of n max and the maximum repetition of m max , the original image can be approximately reconstructed using the formula below:

Property 1. Rotation invariance of FrPHFMs
The amplitudes of FrPHFMs are invariant to image rotation. Assuming that f (r, θ) is an image function in a polar coordinate system, its FrPHFMs is FP nm , image f (r, θ + ϕ) is obtained by rotating the original image ϕ degree and the FrPHFMs thereof is FP nm , then according to the calculation formula, the FrPHFMs in the polar coordinate system is: The amplitudes are taken on both sides of the equation above: From the formula above, it can be known that the amplitudes of the FrPHFMs of the image obtained by rotating the original image are equal to those of the FrPHFMs of the original image, indicating that FrPHFMs are invariant to image rotation. In this way, the angle of rotation, ϕ, can also be estimated by comparing the moments of the two images.

Property 2. Scaling invariance of FrPHFMs
When calculating the scaled FrPHFMs, for a given image function g(r , θ), find the k of the image radius, the range of variation of r will be 0 ≤ r ≤ k, and the normalized image function will be: where the variation range of r = r k is 0 ≤ r ≤ 1. f (r, θ) is the normalized image function, and the FrPHFMs calculated by using the normalized function f (r, θ) has scale invariance.
Because any image f ( r k , θ) obtained by scaling the same image function f (r, θ) is finally normalized to the same function f (r, θ) according to Formula (15), the normalized image FrPHFMs has scaling invariance.

Analysis of Radial Polynomials
Although continuous orthogonal moments have good image description ability, they will be affected by various errors and numerical instability under the condition of high order, and these factors will affect their accuracy. Because such errors have negative effect on image analysis and reconstruction, the image reconstruction performance of continuous orthogonal moments will become very poor when their order reaches the critical value. The properties of continuous orthogonal moments are mainly reflected in their radial polynomials. In this section, the properties of radial polynomials are analyzed. Two groups of test images are shown in the figures below. Figure 1 shows the plots of FrPHFMs radial polynomials versus r with the order n = 85 and different t values, where the value range of r is 0 ≤ r ≤ 1. When t = 1, FrPHFMs will be PHFMs. It can be seen from Figure 1 that the radial polynomial variation of integer-order PHFMs is unstable, while the radial polynomial variation of FrPHFMs is relatively stable. It can also be observed in Figure 1 that the radial polynomials of integerorder PHFMs change rapidly near the point where r = 1, resulting in numerical instability near this point. When t = 0.7, 0.8, 0.9, the radial polynomials of FrPHFMs change in a relatively stable manner, effectively reducing the errors occurring in the edge of images reconstructed with PHFMs. tion results. When the radial polynomials change too fast, the radial polynomials oscillates around r axis − at a higher frequency, which leads to the fact that the radial polynomials cannot be correctly represented by a single value at the pixel center, resulting in poor reconstruction effect and unclear reconstructed image. On the contrary, when the value of t is fractional parameter, the changes in the radial polynomials of FrPHFMs are stable, indicating that FrPHFMs can achieve clear display effects of reconstructed images, reduce experimental errors, and effectively mitigate the deficiency of PHFMs.  To show the changing degree of the radial polynomial versus r more clearly, we calculate the derivative of the radial polynomial to r as the change rate of FrPHFMs radial polynomial. Figure 2 shows the change rate of the corresponding radial polynomials with the increase of r. It can be seen from Figure 2 that the change rate of radial polynomials of integer-order PHFMs is generally higher than the FrPHFMs with the same order near r = 1. And the radial polynomials of PHFMs change rapidly, resulting in numerical instability in the edge regions of images [33] and very poor image reconstruction results. However, when t = 0.7, 0.8, 0.9, the change of radial polynomials near r = 1 is relatively stable, and the change rate is small, thereby, realizing greatly improved image reconstruction results. When the radial polynomials change too fast, the radial polynomials oscillates around r − axis at a higher frequency, which leads to the fact that the radial polynomials cannot be correctly represented by a single value at the pixel center, resulting in poor reconstruction effect and unclear reconstructed image. On the contrary, when the value of t is fractional parameter, the changes in the radial polynomials of FrPHFMs are stable, indicating that FrPHFMs can achieve clear display effects of reconstructed images, reduce experimental errors, and effectively mitigate the deficiency of PHFMs.

Experiments and Analysis of Experimental Results
In this section, the performance of FrPHFMs in image reconstruction and object recognition, as well as their geometric invariance, are tested through a series of experiments. Thirty grayscale images with a size of 128 × 128 were used as the test images in these experiments. Figure 3 shows 10 grayscale images randomly selected from the USC-SIPI Image Database.

Experiments and Analysis of Experimental Results
In this section, the performance of FrPHFMs in image reconstruction and object recognition, as well as their geometric invariance, are tested through a series of experiments. Thirty grayscale images with a size of 128 × 128 were used as the test images in these experiments. Figure 3 shows 10 grayscale images randomly selected from the USC-SIPI Image Database.

Image Reconstruction
Image reconstruction performance is an important feature of the continuous orthogonal moments [34] and reflects the accuracy of such moments in image reconstruction. When FrPHFMs are used to reconstruct images, the computational efficiency can be improved by limiting the number of FrPHFMs involved in image reconstruction. In other words, given , the formula for image reconstruction using FrPHFMs can be written as:

Image Reconstruction
Image reconstruction performance is an important feature of the continuous orthogonal moments [34] and reflects the accuracy of such moments in image reconstruction. When FrPHFMs are used to reconstruct images, the computational efficiency can be improved by limiting the number of FrPHFMs involved in image reconstruction. In other words, given n + |m| ≤ K (0 ≤ n ≤ n max , 0 ≤ |m| ≤ m max ), where K is a constant [31], the formula for image reconstruction using FrPHFMs can be written as: Similarly, other types of moments also are subject to reconstruction constraints. The reconstruction constraints and the numbers of moments used for image reconstruction for three types of moments are listed in Table 1. Table 1. Reconstruction constraints and the numbers of moments used for image reconstruction.

Moment Limit Condition Number of Moments
FrPHFMs

Experiment 1.
A grayscale image named Lena with a size of 128 × 128 pixels was used for this experiment, the max moment order n max was set to varying number from 50 to 90 at an interval of 5, and t = 0.7, 0.8, 0.9, 1.0, 1.1, 1.2. The number of moments used to reconstruct the original image is 2601, 3136, 3721, 4356, 5041, 5776, 6561, 7396 and 8281 respectively. Unless otherwise specified, the value of K in this paper is K = n max .The experimental results are summarized in Table 2. For the purpose of comparison, the results of image reconstruction with integer-order PHFMs (t = 1) are also listed in Table 2, in which ε denotes the mean square reconstruction error (MSRE) [14], which is defined as: where f (x, y) denotes the original image with a size of M × N, and f (x, y) denotes the reconstructed image with the same size.
It can be seen from the data in Table 2 that, within a certain order range, with the increase of n max , the reconstruction effect of FrPHFMs becomes better, the reconstructed image becomes clearer, and the error becomes smaller. When it exceeds this range, the error will gradually increase, and the reconstructed image will become blurred, which Iis caused by the instability of the radial polynomial. When n max is about 70, with the increase of the order, the reconstruction error of FrPHFMs with t = 0.8, 0.9 shows a trend of overall decrease compared with that of integer-order PHFMs. At the same time, with the increase of n max , the edge part of the reconstructed image with the same t value will appear unclear region, which can be explained by the change of radial polynomial. The change of fractional orthogonal moments is mainly realized by changing the radial polynomial. For the same order, the radial polynomial can be adjusted by changing the value of fractional parameter t, which slows down the appearance of white regions and then changes the image reconstruction effect.
Subsequently, the reconstruction errors produced by FrPHFMs were compared with those by PHFMs (t = 1), and the mean value of reconstruction errors was taken for the 30 images. The experimental results are shown in Figure 4, in which reconstruction errors are measured by MSRE. The line chart above clearly shows that, when n max is greater than 70, the reconstruction errors with respect to the images processed with integer-order PHFMs become greater rapidly, and, when t < 1, the reconstruction errors produced by FrPHFMs are remarkably smaller than those by integer-order PHFMs, indicating that FrPHFMs can effectively mitigate the trend of overall increase in image reconstruction errors.   Figure 4, in which reconstruction errors are measured by MSRE. The line chart above clearly shows that, when m ax n is greater than 70, the reconstruction errors with respect to the images processed with integer-order PHFMs become greater rapidly, and, when 1 t < , the reconstruction errors produced by FrPHFMs are remarkably smaller than those by integer-order PHFMs, indicating that FrPHFMs can effectively mitigate the trend of overall increase in image reconstruction errors.

Experiment 2.
Considering the serious effect of noise on image reconstruction [35], the trend of change in reconstruction errors with the increase in order was tested when different values of t were taken after the same salt and pepper noises were added. Thirty images were used for this experiment, the mean value of reconstruction errors was taken for the 30 images, m ax n was set to varying number within the range from 10 to 100, 0.7,0.8,0.9,1.0 t = , and salt and pepper noises at intensity levels of 0.05, 0.1, 0.15, and 0.2 were added, respectively. The line chart for image reconstruction errors is shown in Figure  5.

Experiment 2.
Considering the serious effect of noise on image reconstruction [35], the trend of change in reconstruction errors with the increase in order was tested when different values of t were taken after the same salt and pepper noises were added. Thirty images were used for this experiment, the mean value of reconstruction errors was taken for the 30 images, n max was set to varying number within the range from 10 to 100, t = 0.7, 0.8, 0.9, 1.0, and salt and pepper noises at intensity levels of 0.05, 0.1, 0.15, and 0.2 were added, respectively. The line chart for image reconstruction errors is shown in Figure 5. The results of this experiment show that, as the density of added salt and pepper noise increases, the image reconstruction effects will deteriorate, indicating that noise does affect image reconstruction. When 1 t < , the reconstruction effects of FrPHFMs are always better than those of integer-order PHFMs. When m ax n is greater than 85, the reconstruction errors produced by PHFMs become greater rapidly, while the reconstruction errors produced by FrPHFMs are smaller than those by integer-order PHFMs, indicating that FrPHFMs can effectively mitigate the trend of overall increase in image reconstruction errors.

Comparison between FrPHFMs and Other Fractional-Order Continuous Orthogonal Moments
In this subsection, FrPHFMs are compared with FrRHFMs and FrEMs from the perspective of image reconstruction effects. Different values of t are required to allow different fractional-order moments to achieve good image reconstruction effects. The results of many experiments indicate that the optimal reconstruction effects achievable by three different types of moments correspond to different values of t . Therefore, the image reconstruction performance of FrPHFMs (  Table 3.  The results of this experiment show that, as the density of added salt and pepper noise increases, the image reconstruction effects will deteriorate, indicating that noise does affect image reconstruction. When t < 1, the reconstruction effects of FrPHFMs are always better than those of integer-order PHFMs. When n max is greater than 85, the reconstruction errors produced by PHFMs become greater rapidly, while the reconstruction errors produced by FrPHFMs are smaller than those by integer-order PHFMs, indicating that FrPHFMs can effectively mitigate the trend of overall increase in image reconstruction errors.

Comparison between FrPHFMs and Other Fractional-Order Continuous Orthogonal Moments
In this subsection, FrPHFMs are compared with FrRHFMs and FrEMs from the perspective of image reconstruction effects. Different values of t are required to allow different fractional-order moments to achieve good image reconstruction effects. The results of many experiments indicate that the optimal reconstruction effects achievable by three different types of moments correspond to different values of t. Therefore, the image reconstruction performance of FrPHFMs (t = 0.7, 0.8, 0.9, 1.0) was compared with that of FrRHFMs (t = 1.2, 1.3, 1.4, 1.0) and FrEMs (t = 1.3, 1.4, 1.5, 1.0), and the max moment order was set to n max = 50, 55, . . . , 75, respectively. The results of comparison of images reconstructed with different moments are summarized in Table 3. In order to demonstrate image reconstruction errors more visually, a line chart was drawn to allow for comparison of the reconstruction errors produced by the three different moments. Curves in different colors represent the reconstruction errors produced by different moments for certain values of t, as shown in Figure 6. From the line chart above, it can be seen that the reconstruction errors produced by FrPHFMs are very small and are much smaller than those produced by FrEMs; and, as n max increases continuously, the reconstruction errors produced by the three different moments become greater gradually, but those produced by FrPHFMs become greater more slowly, indicating that FrPHFMs can effectively mitigate the trend of increase in image reconstruction errors. These results show that FrPHFMs are superior to FrRHFMs and FrEMs in terms of image reconstruction performance. In order to demonstrate image reconstruction errors more visually, a line ch drawn to allow for comparison of the reconstruction errors produced by the thre ent moments. Curves in different colors represent the reconstruction errors prod different moments for certain values of t , as shown in Figure 6. From the line char it can be seen that the reconstruction errors produced by FrPHFMs are very small much smaller than those produced by FrEMs; and, as m ax n increases continuou reconstruction errors produced by the three different moments become greater gr but those produced by FrPHFMs become greater more slowly, indicating that F can effectively mitigate the trend of increase in image reconstruction errors. Thes show that FrPHFMs are superior to FrRHFMs and FrEMs in terms of image recons performance.

Experiments on Geometric Invariance
The amplitudes of FrPHFMs are invariant to rotation and scaling. In other wo amplitudes of the FrPHFMs of the rotated and scaled images will be approximate to those of the FrPHFMs of the original images. The rotation invariance and scalin

Experiments on Geometric Invariance
The amplitudes of FrPHFMs are invariant to rotation and scaling. In other words, the amplitudes of the FrPHFMs of the rotated and scaled images will be approximately equal to those of the FrPHFMs of the original images. The rotation invariance and scaling invariance of FrPHFMs are experimentally demonstrated in this section.
(1) Rotation invariance The amplitudes of the FrPHFMs of a group of rotated images were calculated. A grayscale image named Lena with a size of 128 × 128 pixels was used for this experiment. After the image was rotated 5 • , 15 • , 25 • , 35 • , and 45 • , respectively, when t = 0.7, 0.8, 0.9, the amplitudes of the FrPHFMs of the rotated images were calculated, respectively, and compared with those of the FrPHFMs of the original image. Figure 7 shows the rotated versions of image Lena, and the experimental results are summarized in Tables 4-6 below.    The amplitudes of 9 moments of the rotated images and the original image when different values of t are taken are given in the tables above. Through comparison of the listed data, it is found that the amplitudes of the same FrPHFMs of all rotated images are approximately equal, indicating that the amplitudes of FrPHFMs are invariant to rotation [36].
(2) Scaling invariance The FrPHFMs of a group of scaled images were calculated to verify the scaling invariance thereof. A grayscale image named Lena with a size of 128 × 128 pixels was used for this experiment. After the image was scaled by 0.5, 0.7, 1.25, and 1.5 times, respectively, the amplitudes of the FrPHFMs of the scaled images were calculated, respectively, and compared with those of the FrPHFMs of the original image. Figure 8 shows the scaled versions of image Lena, and the experimental results are summarized in Tables 7-9 below.  The amplitudes of 9 moments of the scaled images and the original image when different values of t are taken are given in the tables above. Through comparison of the listed data, it is found that the amplitudes of the same FrPHFMs of all scaled images are approximately equal, indicating that the amplitudes of FrPHFMs are invariant to scaling.

Object Recognition
In this section, the proposed FrPHFMs are experimentally compared with other fractional-order continuous orthogonal moments, including fractional-order polar complex exponential transforms (FrPCETs), fractional polar cosine transforms (FrPCTs), fractional polar sine transforms (FrPSTs), and FrZMs, from the perspective of object recognition. Firstly, relevant operations were performed on the test images, the moments of tested images were calculated, and these moments were used as image features. Because errors will occur when orthogonal moments are calculated, accurate moments need to be selected and used for comparative experiments on object recognition in order to ensure accurate object recognition. The sets of different fractional-order orthogonal moments and the number of moments in each set used for object recognition are listed in Table 10. It is to be noted that low-order moments are used as image features for comparative experiments on object recognition because such moments can represent more image information using less data. Table 10. Moment sets used in object recognition and the number of moments in moment sets.

Moment Conjugated Moment Final Moment Set Number of Moments
FrPHFMs In this experiment, the binary images of 26 capitalized English letters were used for object recognition, and the relationship between the number of moments used for object recognition and correct classification percentage (CCP) was tested. The size of each image was 64 × 64 pixels, as shown in Figure 9. Firstly, each of the 26 letters was rotated 0 − 180 • at an interval of 5 • , producing 37 rotated images, from which 19 images were selected randomly as training images, and the remaining 18 images were used as test images. Nine hundred and sixty-two test images were produced in total, including 494 images in the training set and 468 images in the test set. Gaussian noise with average intensity of 0 and variance of σ 2 = 0.20 was added to the images in the test set, which were classified using the K-Nearest Neighbors (KNN) algorithm [37]. The changes in the CCP of different fractional-order continuous orthogonal moments with changing number of moments are shown in Figure 10 below.  From the figure above, it can be seen that the CCP of FrPHFMs is higher than those of other fractional-order continuous orthogonal moments, indicating that FrPHFMs are superior to other fractional-order continuous orthogonal moments in terms of classification performance. Additionally, it can also be observed that, as the number of moments increases, the CCP will increase gradually, and it will start to decrease after reaching a certain value, which is known as the overfitting problem in machine learning. Experiment 2. One hundred grayscale images with an equal size of 128 × 128 pixels selected from database COIL-100 were used as the test images for this experiment. Figure 11 shows some images randomly selected from this database. Firstly, each image was rotated 0 − 90 • at an interval of 5 • , producing 19 images, from which 10 images were randomly selected as training images, and the remaining 9 images were used as test images. One thousand nine hundred test images were produced in total, including 1000 images in the training set and 900 images in the test set. Salt and pepper noises with varying density of σ 2 = 0.00, 0.05, · · · , 0.25 were added to the images in the test set, which were classified using the multi-class naive Bayes classifier [38]. Six, 18, 45 and 96 low-order moments were selected and used, respectively, as the image features for this experiment, and the CCPs of these different moments were compared. The experimental results are shown in Figure 12.  From Figure 12, it can be clearly observed that the CCP of FrPHFMs is higher than those of other fractional-order continuous orthogonal moments, indicating that FrPHFMs have the best classification performance. Additionally, it can also be observed that, when nose is not added to the images, the CCP of each moment is 100% but will decline with the increase in noise density. However, the CCP of FrPHFMs decreases very slowly, indicating that FrPHFMs can effectively mitigate the trend of rapid decline of object recognition rate. Experiment 3. The 100 grayscale images used in experiment 2 were used as the training images for this experiment. Various conventional attacks were applied to the images, and the object recognition rates of different moments of images subjected to such attacks were compared. Firstly, all images were rotated 0 • , 5 • , . . . , 90 • , respectively, producing 1900 images in total, which were included into the test set. Then different attacks were applied to the test images. The applied attacks included JPEG compression with quality factor of 10, 20, 50, and 80, Wiener filtering with a window size of 2 × 2 and Gaussian filtering with a window size of 4 × 4, contrast enhancement filtering with alpha factor of 0.2 and 0.4, and circular averaging filtering at radius of 2 and 4. The images in the test set were then classified using a multi-class support vector machine (MSVM) [39]. Five and 13 low-order moments were selected and used, respectively, as the image features for this experiment, and the CCPs of these different moments were compared. The experimental results are shown in Table 11. From the Table 11, it can be seen that the CCPs of all fractional-order continuous orthogonal moments are 100% under various attacks because such moments are highly robust to attacks. It can also be observed that, regardless of which attack is applied to the images, the CCP of FrPHFMs is always higher than those of other fractional-order continuous orthogonal moments, indicating that FrPHFMs are highly robust and superior to other fractional-order continuous orthogonal moments in terms of resistance to attacks.

Conclusions
In this paper, in order to improve the anti-noise and reconstruction performance of PHFMs, the traditional PHFMs, which can only take integer-order, are extended to FrPHFMs. By modifying the radial polynomial of integer-order PHFMs, FrPHFMs are constructed according to fractional radial polynomial, and the properties of FrPHFMs are introduced and detailed experiments are carried out according to their properties. Firstly, the traditional integer-order PHFMs are introduced, and then FrPHFMs are constructed by using fractional radial polynomial, and their properties are described in detail. FrPHFMs have good orthogonality, rotation invariance, and scaling invariance and are superior to integer-order PHFMs. Secondly, the change of radial polynomial is analyzed in detail. Finally, the constructed FrPHFMs are applied to image reconstruction, geometric invariance, and object recognition experiments, which further verifies their good geometric invariance and image description ability. From the numerical and experimental analysis, the following conclusions can be drawn: FrPHFMs not only maintain the orthogonality, rotation invariance, and scaling invariance of integer-order PHFMs, but they also have good image description ability. Their performance in image reconstruction, anti-noise performance, and object recognition is better than integer-order PHFMs and other fractional-order continuous orthogonal moments. In the future, the improvement of FrPHFMs performance will be made an important area of research.