An Improved Algorithm for Pseduo-Jacobi-Fourier Moments

Image moments have been used in many research fields of the engineering. However, the related computation of invariant moments mostly adopted the polar coordinate system, which not only increase the computational load, but also cause large quantized error. To solve this problem, an improved algorithm to compute Pseudo-Jacobi-Fourier moments in the Cartesian coordinate system is proposed in this paper. The experimental results show that the reconstructed image with improved PJFM’s has more advantages than polar coordinate system, such as more information, fewer moments, less time consuming. And the recognition rate of the microscopic images of 8 helminth eggs was also higher than in polar coordinate system.


Introduction
Image moments, because of their powerful description of the image content, have been used in many research fields of the engineering, such as image processing [1] [2], pattern recognition and machine vision [3].The first introduction of image moments for classification purposes was performed by Hu [4], and then developed into other families, such as Legendre [5], Zernike [6], Pseudo-Zernike [7], Fourier-Mellin [8], Tchebichef [9], Krawtchouk [10], Pseudo-Jacobi-Fourier [11] moment.These moments can be used as image descriptors after an appropriate normalization procedure in order to achieve translation, scale and rotation invariance.However, the related computation for these moments mostly adopted the polar coordinate system, which not only increase the computational load, but also cause large quantized error [12] [13].To solve this G.Amu DOI: 10.4236/opj.2017.78B00969 Optics and Photonics Journal problem, an improved algorithm to compute Pseudo-Jacobi-Fourier moments in the Cartesian coordinate system is proposed in this paper.This improved algorithm is applied to classify the microscopic images of helminth eggs by using Euclidean distance classifier [14], and the recognition rate is 92.2%.

Definition of Pseduo-Jacobi-Fourier Moments
Bhatia and Wolf have shown [15] that a polynomial that is invariant in form for any rotation of axes about the origin must be of the form ( cos , sin ) ( ) exp( ) where the function where nk δ is Kronecker symbol, and ( ) ( ) n n n q n p n p q b p q n p n p The radial polynomials of OFMM's, CHM's and ZM's belong to Jacobi polynomials with 2, 1 2 p q p q = = = = and 1 p q m = = + respectively.When 4, 3 p q = = , radial Jacobi polynomials become ( ) ( ) ( ) So a new set of orthogonal polynomial function as an infinite series expansion in terms of ( ) where then nm φ is defined as Pseudo-Jacobi-Fourier Moments (PJFM's), and 1 r = is the maximum size of the objects that can be encountered in a particular application.
Actually, most of the images are defined in the Cartesian coordinate system.
When calculating PJFM's, the images need to be converted to polar coordinates, so may cause some problems: 1) Quantization error is introduced into the calculation of PJFM's, which increases calculating amount and additional noise; 2) In polar coordinates, the discrete points near from the origin are more than those in the Cartesian coordinate system, which lead to information redundancy; 3) Far from the origin, the scattered points in the polar coordinate system are less than those in the Cartesian coordinate system, which lead to information loss.Therefore, an improved algorithm for computing PJFM's Cartesian coordinate system is developed in this paper.

Improved Algorithm for Pseduo-Jacobi-Fourier Moments
When calculating PJFM's in the Cartesian coordinates, the image should first be normalized into the unit circle, and the integral region of Equation ( 11 Figure 1 shows that s1 -s2 is not the part of the image and there is no pixel over it, Equation ( 15) can be written as So Equation ( 16) is the general formula for calculating PJFM's in Cartesian coordinate system.

Image Reconstruction Using Improved Algorithm
Image reconstruction can be used as an effective means to evaluate the quality of feature extraction.The more PJFM's used to reconstruct images, the closer to the original image.An image of capital E, shown in Figure 2, is reconstructed by Equation (17).

ˆ( , ) ( )
The four corners of the reconstructed image are black and some pixels have been lost in the Figure 3.But the black area of the four corners of the reconstructed image is shrinking with the increase of N and M, and the edge information is well preserved in the Figure 4. Figure 4 also shows that E is differentiated well when N = M = 6 instead of N = M = 10 as in Figure 3.
As for computation time, reconstruction time in the Cartesian coordinate system is much shorter than the polar coordinates system, shown in Table 1.

Image Recognition Using Improved Algorithm
In order to test the feature extraction performance of the improved PJFM's, a The training sets consist of 20 different versions from each kind of helminth eggs, including 160 images.Testing sets consist of 307 untrained images from different version of 8 kinds of helminth eggs.Figure 5 gives a multi-distorted image of the microscopic image of some helminth eggs in a testing set.
Choosing N = M = 8, the image feature extracted respectively by using Equation (16) in Cartesian coordinate system and Equation (11) where nm φ is the PJFM of the testing object, and ( ) φ is the PJFM of the reference object of class i.
As can be seen from Table 2, the total recognition rate of the microscopic images of 8 helminth eggs was 92.2% with improved PJFM's, 7.2% higher than PJFM's in polar coordinate system.

Conclusion
An improved algorithm to compute Pseudo-Jacobi-Fourier moments in the Cartesian coordinate system is proposed in this paper.The experimental results show that the reconstructed image with improved PJFM's has more advantages than polar coordinate system, such as more information, fewer moments, less time consuming.And the recognition rate of the microscopic images of 8 helminth eggs was also higher than in polar coordinate system.

9 )
According to the orthogonal theory, the image function ( , ) f r θ can be written G. Amu DOI: 10.4236/opj.2017.78B00970 Optics and Photonics Journal

Figure 2 .
Figure 2. The original image of E.
in polar coordinate system, then the target objects were recognized by the minimum average distance rules .Euclidean distances are calculated by Equation (18), and the result is shown in Table2.

Table 2 .
Recognition result of parasite egg microscopic images