|Home | About | Journals | Submit | Contact Us | Français|
A detailed analysis of the quaternion generic Jacobi-Fourier moments (QGJFMs) for color image description is presented. In order to reach numerical stability, a recursive approach is used during the computation of the generic Jacobi radial polynomials. Moreover, a search criterion is performed to establish the best values for the parameters and of the radial Jacobi polynomial families. Additionally, a polar pixel approach is taken into account to increase the numerical accuracy in the calculation of the QGJFMs. To prove the mathematical theory, some color images from optical microscopy and human retina are used. Experiments and results about color image reconstruction are presented.
In recent years, moments have been successfully used in a variety of research areas such as image registration,1 face recognition,2 angle estimation,3 watermarking,4 pattern reconstruction,5 medical imaging,6–8 focus measures,9 image analysis,10 forensic applications,11 gait phase detection,12 and so forth. In the 1960s, Hu13 introduced a set of invariants based on the low-order geometric moments for pattern recognition tasks. Almost two decades later, Teague14 proposed Zernike and Legendre moments derived from the basis set of orthogonal polynomials. It is known that one of the main advantages of moments with orthogonal basis functions is the ability to represent the global features of an image by a set of mutually independent descriptors, with a minimal amount of information redundancy. As stated by Bhatia and Wolf,15 there are an infinite number of orthogonal sets inside the unit circle. Some orthogonal moments are derived from the basis set of pseudo-Zernike,16 Chebyshev-Fourier,17 orthogonal Fourier-Mellin,18 radial harmonic Fourier,19 and Bessel-Fourier20 polynomials. Additionally, new orthogonal basis sets of circular moments have been proposed from the generic formula of the Jacobi radial polynomials,21–23 where each set can be generated by combinations of two real parameters, which are commonly denoted as and . Also, Jacobi-Fourier moments (JFMs) have been successfully proven in pattern recognition,24 image analysis,25 and machine vision applications.26 Recently, quaternion color moments have been introduced on the basis set of Fourier-Mellin,27 Zernike,28 and Bessel-Fourier polynomials.29 These works have extended the capabilities of moments for describing color or images. Consequently, new invariant descriptors have been proposed for pattern recognition while new algorithms for image reconstruction also have been studied. In the same direction, Karakasis et al.30 presented a general framework for computing accurate quaternion color moments based on different orthogonal polynomials, from continuous to discrete polynomials and from Cartesian to circular symmetries.
On the other hand, in Ref. 23, there exists an analysis of the interconversion between the two historical definitions of the radial Jacobi polynomials; this interesting paper points out and clears the confusing points between both definitions. However, the computations of the polynomials through these formulas take into account the calculation of factorials of high orders. Unfortunately, this method of computing the radial polynomials causes strong numerical instability and inaccuracy for orders higher than . This effect is most evident in the reconstruction of large images as those acquired in the biomedical areas.
In many biomedical laboratories, it is a common task to count some microscopic specimens in order to help the clinical analysis or during the diagnostic of some diseases. A quality control program for clinical laboratories is commonly described by features such as accuracy, precision, tendency to avoid erroneous results, and tendency to avoid systematically different results in different laboratories. All of these situations should be supported by the use of automatic methods of pattern recognition and by the employment of digital image techniques, which must be based on numerical features increasingly more accurate and faster during the data processing.
In another context, retina images have been of interest in many scientific analyses as the biometric identification of people, among other applications. Due to this fact, new numerical approaches to analyze these kinds of images are continuously required. In this paper, we research the capabilities of image description from the quaternion generic Jacobi-Fourier moments (QGJFMs). This investigation makes use of the variations of the Jacobi polynomials by the and parameters, because they can use different sets of quaternion orthogonal moments such as quaternion orthogonal Fourier-Mellin moments (, ), quaternion Chebyshev-Fourier moments (, ), quaternion pseudo Jacobi-Fourier moments (, ), quaternion Legendre-Fourier moments (, ), quaternion pseudo-Zernike moments, and quaternion Zernike moments. The analysis is done in the context of a polar pixel scheme and by means of an algorithm for fast computation and high precision.25,26 Furthermore, this study proposes the search for optimal parameters and of QGJFMs to significantly improve the quality of reconstruction. The analysis is then applied to reconstruct some color biomedical images in circular symmetry as validation tests of the algorithm.
This work is organized as follows. Section 2 introduces QGJFMs. In Sec. 3, the QGJFMs are computed on the geometry of polar pixels and presents an image conversion via bicubic interpolation for the high-precision computing of QGJFMs. In Sec. 4, an inverse quaternion transformation for reconstruction of RGB images is presented. Also, in this section, the normalized image reconstruction error is used as a measurement of image reconstruction. Section 5 takes into account a numerical comparison with other methods in the context of a test RGB image for different orders of reconstruction. Some experimental results are presented in Sec. 6, where an evaluation of accurate and fast computation through some color biomedical images has been used. Finally, the conclusions of this research are presented in Sec. 7.
Recently, quaternion moments have been used as a novel method for describing color images. The quaternion theory is a generalization of the complex theory of numbers. A quaternion is composed of one real part and three imaginary parts; it takes the following form:
where , , , and are real numbers and , , and are three imaginary units. The basic operations among them are in accord with the following rules:
The conjugate and modulus of a quaternion are, respectively, defined by
Using the quaternion representation, an RGB image defined in polar coordinates can be considered as a quaternion-based model for color images, which is given as follows:
where the functions , , and are, respectively, the red, green, and blue channels of the image.
According to the noncommutative property of quaternion multiplication, there are two types of QGJFMs of order and repetition , namely, the right-side moments defined by
and the left-side moments, which are given as follows:
where are the generic Jacobi radial polynomials.23 The real numbers and define a specific radial polynomial family and is a unit pure quaternion. It can be used as . Due to the anti-involution property of quaternion conjugation, the left-side and right-side transforms have the following relationship:
For the case of discrete image functions, Eqs. (6) and (7) cannot be directly applied. The integrals of these equations have to be replaced by summations and the discrete image function has to be normalized inside the unit disk; this approach is known as the zeroth-order approximation or direct method. This procedure is briefly described below.
Let be a discrete image function with spatial dimensions . Its right discrete quaternion moments are given by
where the discrete polar coordinates are expressed by
and they are transformed as follows:
where and .
In a previous analysis,31 it has been demonstrated that two types of errors occur in the direct calculation method of orthogonal moments, namely, geometric error and numerical integration error. To increase the numerical accuracy calculation, Xin et al.32 proposed an algorithm based on changing the shape of the square pixels for a polar pixels scheme. Recently, Camacho-Bello et al.25 proposed a recurrence relation to eliminate the numerical instability of the orthogonal polynomials and to compare with other methods in terms of invariance, reconstruction error, and computation time. In terms of our research in this section, a recurrence relationship for computing the generic radial Jacobi polynomials will be used.26
The recurrence relation with respect to is used to increase the numerical stability for the computation of the shifted radial Jacobi polynomials. The recursive relationship is given by
where , , , , and the coefficients and are computed as follows:
To start with the numerical computation, the Jacobi polynomials of zero and first normalized orders are given by
The advantage of the recurrence relationships is that they can diminish the computation time of the radial polynomials because the calculation of factorials is not taken into account. Moreover, using the recurrence relation, the computation accuracy of the ’th power of () is increased.
The unit disk of polar pixels scheme is uniformly divided along the radial direction into sections, with radial distance of , where . The number of polar pixels in the ’th ring is , with angles , where and is the number of sectors contained in the innermost section. In practice, Xin et al.32 recommended to set the value of and for an image. Figure 1 shows the polar pixels scheme with and . The calculation of the QGJFMs in the framework of polar pixels is performed by the summation of all sectors inside the unit disk. For this calculation, Eq. (9) can be rewritten as
where is an approximation of the function defined over a set of concentric sectors , and the factor is given by
(a) Polar pixel representation of an image and (b) polar pixel sector .
The integrals of Eq. (19) can be separated as
where and denote, respectively, the starting and ending radii of the sector , while and denote the starting and ending angles of the sector. A scheme of a polar pixel is shown in Fig. 1(b). Using the Euler formula for quaternions, Eq. (20) can be seen as follows:
Rewriting Eq. (21) in the next form produces
It is evident from Eq. (22) that three integrals can be obtained, which can be expressed as
In order to solve the integral of Eq. (23), Camacho-Bello et al.25 proposed to use a 10-point Gaussian quadrature rule, along with the recurrence relationship with respect to the order for computing the generic radial Jacobi polynomials.
Typically, the aforementioned procedure is less accurate than the direct method, but is numerically more stable for orders greater than 21. Thus, the composite Gaussian quadrature rule for numerical integration of the generic Jacobi radial polynomials can be stated as
where are weights and are the points where the function is evaluated. The values for and are given in Table 1.
Weights () and location of sampling points () for 10-point Gaussian quadrature.
The integrals of Eqs. (24) and (25) can be analytically calculated as follows:
For the calculation of QGJFMs from the right-side transformation in a scheme of polar pixels, we have to represent , , and as the components of an RGB image. These images can be transformed to a scheme of polar pixels.20,25,26 Therefore, from Eq. (5) can be replaced as
Substituting Eq. (29) into Eq. (6) leads to
Using the basic rules of Eq. (2), we can express the QGJFMs as
As stated earlier, can be represented by Eq. (22), and by virtue of this relation, Eq. (31) takes the following form:
At this point, we can again use the basic rules for quaternions defined in Eq. (2). Then, distributing the integrals for both sine and cosine, Eq. (32) can be written as follows:
Taking into account the proposed value for the quaternion , the QGJFMs of an image in polar pixels are expressed by
The capacity of description of the QGJFMs is done by means of its inverse transformation. According to the inverse quaternion Jacobi-Fourier transform, an original image can be reconstructed by a finite number of QGJFMs, which is given by
Note that , , and represent the red, green, and blue components of the reconstructed color image, respectively, and is the maximum order of QGJFMs used in the reconstruction. Typically, the reconstruction formula has been used to recover the image information up to a certain level of approximation.
In this section, a reconstruction metric based on the normalized image reconstruction error (NIRE) is presented. It is defined as the normalized square error between the input image and its reconstruction . This error is presented in discrete form for each channel of the RGB image, and it is given by
The measure for the reconstruction error of a color image that considers the three channels is given by
Also, we have assumed that the image sizes are .
The zero-order approximation is commonly used because of its easy implementation. However, it requires replacing the double integral by a double sum. Typically, this approximation during the computation of the quaternion moments creates a poor reconstruction in RGB images. For this reason, Karakasis et al.30 proposed an alternative numerical approximation that uses some recurrence relations in order to improve the performance calculation. In this section, we compare the approach of polar pixels with zero-order approximation and the numerical approximation in terms of image reconstruction, reconstruction error, and computation time. As in Ref. 30, we use the standard image Lena with the quaternion orthogonal Fourier-Mellin moments (). Figures 2 and and33 show the Lena image reconstruction and the NIRE, respectively.
Lena image reconstruction of with different methods for , 50, 100.
NIRE of Lena image with size for different computation methods.
Since the computation time does not depend on image content, only one image is considered for the analysis. The code is implemented using MATLAB® with a Sony Electronics Inc. Model VAIO R Computer Notebook PC with Intel Core™i5-2430M CPU 2.40 GHz processor with 4 GB of RAM. Figure 4 shows the computation time of the different methods. The zero-order approximation calculation has better times. However, the scheme of polar pixels has reasonable computing times with minimal reconstruction error.
Computation time of Lena image with size for different computation methods.
In this section, we present the reconstruction results obtained using some color test images from optical microscopy and the human retina. These RGB test images are shown in Fig. 5. Using Eq. (38), we have reconstructed the three cases of test images. The results are shown in Figs. 6(a), 6(b), and 6(c).The results for the NIRE from the test images are shown in Figs. 6(d), 6(e), and 6(f). To obtain these results, we have considered image sizes of and , 2, 3, 4, 5.
RGB test images from (a) yogurt bacteria, (b) human retina, and (c) diatom. Yogurt bacteria and diatom are acquired using a microscope in bright field illumination. In the case of the diatom, a 100× microscope objective ...
Results of the test images: (a), (b), and (c) RGB image reconstruction, (d), (e), and (f) RGB image reconstruction error, (g), (h), and (i) searching of and parameters, and (j), (k), and (l) mean of PSNR.
As mentioned in Sec. 1, Bhatia and Wolf15 pointed out that there is an infinite number of complete sets of orthogonal polynomials, which are invariant to rotation and can be obtained from the generic Jacobi polynomial formula. Each set is obtained by the combination of the parameters and . Therefore, a set of orthogonal moments could be chosen to fit specific applications. We define the mean of Mean_NIRE as metric to quantitatively evaluate the best combinations of and . The mean value of Mean_NIRE is given by
where is the maximum order of QGJFMs used in the reconstruction of the input image and is the cutting order. The cutting order can find the optimal parameters for moments of low and high orders. In general, low-order moments represent the global shape of an image, and high-order moments the detail.5 The search space is restricted by and the sections with minimum values represent the neighborhood of the best combinations of and . Figures 6(g), 6(h), and 6(i) show the exhaustive searches performed for the RGB test images with values of and .
In general, image reconstruction has been commonly used to describe how well an image can be retrieved by a small set of its moments. Here, the interest is to discuss the results obtained in each test case.
As is shown in Fig. 6(g), the best values of the search are located when . According to its NIRE graph of Fig. 6(d), the tendency of error is to go down as the order is increased. Moreover, the RGB image is completely reconstructed before the order . Figure 6(a) shows the RGB image reconstruction when , which is the best value found in the search of Fig. 6(g). It is clear that the reconstruction is not severely affected by the Gaussian noise acquired during the registration of the input image.
In this case, the search for the minimum values for the average of the NIRE yields . The NIRE graph of the best values found shown in Fig. 6(e) presents numerical oscillations and the tendency of error after the order to increase. This behavior is explained by the low contrast presented in the input image of Fig. 5(b). As in the previous case, the best value in the search of Fig. 6(h) is found when ; the image reconstruction results are shown in Fig. 6(b). In spite of these results, QGJFMs are capable of recovering the global information of the image including some high frequencies.
In the third case, the best results of the search neighborhood are when . Moreover, the image reconstruction and NIRE with and best parameters present qualitatively and quantitatively good results, respectively. As shown in Fig. 6(f), NIRE decreases as increases, reaching zero at . Unlike the prior case, the diatom image has a clear background and the object of interest is clearly identified. Finally, the RGB image reconstruction shown in Fig. 6(c) has , which are the best parameters found in Fig. 6(i).
For all RGB test images, the best results are found when , especially when , which would be the case of the quaternion Legendre-Fourier moments.
Sensitivity to noise is a critical issue for quaternion moments. We compare the best parameters of QGJFMs () for peak signal-to-noise ratio (PSNR) and their performance for image reconstruction in the presence of noise, which is the ratio between a signal maximum power and the power of the signal noise, expressed in decibels. The PSNR for each channel of the image is given by
where the mean square error () is defined as
The measure for the PSNR of an RGB image that considers the three channels is given by
Figures 6(j), 6(k), and 6(l) show the Mean_PSNR result from the test images corrupted by an additive Gaussian noise. The results of the graphs show that the QGJFMs when remain unchanged for values , i.e., are resistant to Gaussian noise for quaternion moments of low order.
We have introduced QGJFMs for color image description in the scheme of polar pixels. We have proven a reconstruction algorithm for measuring the capacity of the QGJFMs in retrieving the features of some test RGB images. Also, we have a approach for high precision and fast computation of QGJFMs compared to other methods. Furthermore, the arrangement of polar pixels and the recurrence relationship used for different polynomial families improves the computation time and the reconstruction error in the case of frames with large sizes.
The three test images used here present interesting characteristics. In the yogurt bacteria image, the background and object of interest are near the original colors. However, the QGJFMs have been able to recover all the information. In the case of the human retina image, new QGJFMs have obtained global information and the main high frequencies. Finally, for the diatom image, the background is strongly different from the specimen; this fact helps to practically recover all the RGB image. In general, these results have been found using the and parameters through the searching spaces for each particular case.
C. Camacho-Bello thanks CONACyT for the scholarship with number 423493.
César Camacho-Bello received his BS degree in industrial engineering from the Autonomous University of Hidalgo State in 2006, his MS degree in optical computing from the Polytechnic University of Tulancingo in 2011, and his PhD in optics from the National Institute of Astrophysics, Optics, and Electronics in 2014. He is now a researcher in the Polytechnic University of Tulancingo, Mexico. His research interests include pattern recognition, biometric analysis, computer vision, and digital image processing.
Alfonso Padilla-Vivanco received his physics degree in 1994 from Puebla Autonomous University, Mexico, and his MS degree in 1995 and his PhD in 1999, both in optics, from the National Institute of Astrophysics, Optics, and Electronics. In 2000, he held a postdoctoral position in the Physics Department at the University of Santiago de Compostela, Spain. He is now a researcher at Polytechnic University of Tulancingo, Hidalgo, Mexico. His research interests include optical information processing, image analysis, and computer vision.
Carina Toxqui-Quitl is a researcher at the Polytechnic University of Tulancingo. She received her BS degree from the Autonomous University of Puebla, Mexico, in 2005 and her MS and PhD degrees in optics from the National Institute of Astrophysics, Optics, and Electronics in 2006 and 2010, respectively. Her current research areas include image moments, digital image processing, biomedical optics, and computer vision. She is also a member of SPIE.
José Javier Báez-Rojas received his MSc and PhD degrees from the Tokyo Institute of Technology, Japan, in 1991 and 1994, respectively. He is currently a researcher with the National Institute of Astrophysics, Optics, and Electronics in Mexico. His research interests include digital color, multispectral imagery, and digital image processing.