Pattern Recognition Letters 28 (2007) 1688–1704 www.elsevier.com/locate/patrec
Image analysis by discrete orthogonal dual Hahn moments
Hongqing Zhu a, Huazhong Shu a,c,*, Jian Zhou a, Limin Luo a,c, J.L. Coatrieux b,c
a Laboratory of Image Science and Technology, Department of Computer Science and Engineering, Southeast University, 210096 Nanjing, People’s Republic of China b INSERM U642, Laboratoire Traitement du Signal et de l’Image, Universite´ de Rennes I, 35042 Rennes, France c Centre de Recherche en Information Biome´dicale Sino-franc¸ais (CRIBs), France
Received 29 May 2005; received in revised form 26 February 2007 Available online 29 April 2007
Communicated by L. Younes
Abstract
In this paper, we introduce a set of discrete orthogonal functions known as dual Hahn polynomials. The Tchebichef and Krawtchouk polynomials are special cases of dual Hahn polynomials. The dual Hahn polynomials are scaled to ensure the numerical stability, thus creating a set of weighted orthonormal dual Hahn polynomials. They are allowed to define a new type of discrete orthogonal moments. The discrete orthogonality of the proposed dual Hahn moments not only ensures the minimal information redundancy, but also elim- inates the need for numerical approximations. The paper also discusses the computational aspects of dual Hahn moments, including the recurrence relation and symmetry properties. Experimental results show that the dual Hahn moments perform better than the Legendre moments, Tchebichef moments, and Krawtchouk moments in terms of image reconstruction capability in both noise-free and noisy con- ditions. The dual Hahn moment invariants are derived using a linear combination of geometric moments. An example of using the dual Hahn moment invariants as pattern features for a pattern classification application is given. 2007 Elsevier B.V. All rights reserved.
Keywords: Discrete orthogonal moments; Dual Hahn polynomials; Image reconstruction; Pattern classification
1. Introduction ric moments are not orthogonal. The lack of orthogonality leads to a certain degree of information redundancy and Since Hu (1962) introduced moment invariants, causes the recovery of an image from its geometric moments and functions of moments due to their ability moments strongly ill-posed. Indeed, the fundamental rea- to represent global features of an image have found wide son for the ill-posedness of the inverse moment problem applications in the fields of image processing and pattern is the serious lack of orthogonality of the moment sequence recognition such as noisy signal and image reconstruction (Talenti, 1987). To surmount this shortcoming, Teague (Yin et al., 2002; Mukundan et al., 2001b), image indexing (1980) suggested the use of the continuous orthogonal (Mandal et al., 1996), robust line fitting (Kiryati et al., moments defined in terms of Legendre and Zernike polyno- 2000), and image recognition (Qing et al., 2004). Among mials. Since the Legendre polynomials are orthogonal over the different types of moments, the Cartesian geometric the interval [ 1,1], and Zernike polynomials are defined moments are most extensively used. However, the geomet- inside the unit circle, the computation of these moments requires a suitable transformation of the image co-ordinate space and an appropriate approximation of the integrals * Corresponding author. Present address: Laboratory of Image Science (Shu et al., 2000; Chong et al., 2003). and Technology, Department of Computer Science and Engineering, Southeast University, 210096 Nanjing, People’s Republic of China. Tel.: The study of classical special functions and in particular +86 25 83794249; fax: +86 25 83794298. the discrete orthogonal polynomials has recently received an E-mail address: [email protected] (H. Shu). increasing interest (Mukundan et al., 2001a; Arvesu´ et al.,
0167-8655/$ - see front matter 2007 Elsevier B.V. All rights reserved. doi:10.1016/j.patrec.2007.04.013 H. Zhu et al. / Pattern Recognition Letters 28 (2007) 1688–1704 1689
2003; Foupouagnigni and Ronveaux, 2003; Ronveaux et al., dratic lattices x(s)=s(s + 1)). The dual Hahn polynomials 2000). The use of discrete orthogonal polynomials in image are scaled, to ensure that all the computed moments have analysis was first introduced by Mukundan et al. (2001b) equal weights, and are used to define a new type of discrete who proposed a set of discrete orthogonal moment func- orthogonal moments known as dual Hahn moments. Simi- tions based on the discrete Tchebichef polynomials. An effi- lar to Tchebichef and Krawtchouk moments, there is no cient method for computing the discrete Tchebichef need for spatial normalization; hence, the error in the com- moments was also developed (Mukundan, 2004). Another puted dual Hahn moments due to discretization does not new set of discrete orthogonal moment functions based on exist. However, our new moments are more general because the discrete Krawtchouk polynomials was presented by both the Tchebichef and Krawtchouk polynomials are spe- Yap et al. (2003). It was shown (Mukundan et al., 2001b; cial cases of the dual Hahn polynomials (Nikiforov and Yap et al., 2003) that the discrete orthogonal moments per- Uvarov, 1988). Since the dual Hahn moments contain more form better than the conventional continuous orthogonal parameters (due to the fact that the dual Hahn polynomials moments in terms of image representation capability. are defined on the non-uniform lattice) than the discrete The classical orthogonal polynomials can be character- Tchebichef and Krawtchouk moments, these extra parame- ized by the existence of a differential equation. In fact, ters give more flexibility in describing the image, the the classical continuous polynomials (e.g., Hermite, improvement in performance could thus be expected. Laguerre, Jacobi, Bessel) satisfy a differential equation of It is worth mentioning that although the dual Hahn the form (Nikiforov and Uvarov, 1988) polynomials are orthogonal on a non-uniform lattice, the r~ðxÞy00ðxÞþ~sðxÞy0ðxÞþkyðxÞ¼0 ð1Þ discrete dual Hahn moments defined in this paper are still applied to uniform pixel grid image. The difference between where r~ðxÞ and ~sðxÞ are polynomials of at most second and the dual Hahn moments and the discrete moments based first degree, and k is an appropriate constant. It is possible on the polynomials that are orthogonal on uniform lattice to expand polynomial solutions of partial differential equa- (e.g., discrete Tchebichef moments and Krawtchouk tions in any basis of classical orthogonal polynomials. moments) is that the latter is directly defined on the image When the differential equation (1) is replaced by a differ- grid but, for the former, we should introduce an intermedi- ence equation, we can find the main properties of classical ate, non-uniform lattice, x(s)=s(s + 1). orthogonal polynomials of a discrete variable. Consider the The rest of the paper is organized as follows. In Section simplest case, when (1) is replaced by the following differ- 2, we introduce the dual Hahn polynomials of a discrete ence equation: variable. This section also provides the derivation of 1 yðx þ hÞ yðxÞ yðxÞ yðx hÞ weighted dual Hahn polynomials and the dual Hahn r~ðxÞ moments. Section 3 discusses the computational aspects h h h ~sðxÞ yðx þ hÞ yðxÞ yðxÞ yðx hÞ of dual Hahn moments. It is shown how the recurrence for- þ þ þ kyðxÞ¼0 2 h h mulae and symmetry property of the dual Hahn polynomi- als can be used to facilitate the moment computation. The ð2Þ dual Hahn moment invariants are also derived in this sec- which approximates (1) on a lattice of constant mesh tion. The experimental results are provided and discussed Dx = h. The classical discrete polynomials such as Charlier, in Section 4. Section 5 concludes the paper. Meixner, Tchebichef, Krawtchouk, and Hahn polynomials are all the polynomial solutions of (2). 2. Dual Hahn moments After a change of independent variable x = x(s)in(2),we can obtain a further generalization when (1) is replaced by a 2.1. Discrete orthogonal polynomials on the non-uniform difference equation on a class of lattices with variable mesh lattice Dx = x(s + h) x(s)(Nikiforov and Uvarov, 1988). The discrete orthogonal polynomials on the non-uniform lattice Let us first review some general properties of orthogonal are of great importance for applications in quantum integral polynomials of a discrete variable on a non-uniform lattice systems, quantum field theory and statistical physics (Vega, (Nikiforov and Uvarov, 1988; A´ lvarez-Nodarse and Smir- 1989). In recent years, much attention has been paid to the nov, 1996). As previously indicated, the discrete orthogo- study of this class of polynomials (Koepf and Schmersau, nal polynomials on the non-uniform lattice can be 2001; Kupershmidt, 2003; A´ lvarez-Nodarse et al., 2001a; constructed using a variable mesh in (1) when it is replaced A´ lvarez-Nodarse and Costas-Santos, 2001b; A´ lvarez-Nod- by a difference equation. Let arse and Smirnov, 1996; Temme and Lo´pez, 2000). How- ever, to the best of our knowledge, until now, no discrete D ryðsÞ ~sðxðsÞÞ DyðsÞ ryðsÞ r~ðxðsÞÞ 1 þ þ þ kyðsÞ¼0 orthogonal polynomials defined on a non-uniform lattice Dxs 2 rxðsÞ 2 DxðsÞ rxðsÞ have been used in the field of image analysis. In this paper, ð3Þ we address this problem by introducing a new set of discrete orthogonal polynomials, namely the dual Hahn polynomi- be the second order difference equation of hypergeometric als, which are orthogonal on a non-uniform lattice (qua- type for some lattice function x(s), where 1690 H. Zhu et al. / Pattern Recognition Letters 28 (2007) 1688–1704
2 rgðsÞ¼gðsÞ gðs 1Þ; DgðsÞ¼gðs þ 1Þ gðsÞð4Þ where dn denotes the square of the norm of the correspond- ing orthogonal polynomials, and q(s) is a non-negative denote respectively the backward and forward finite differ- function (weighting function), i.e., ence quotients. r~ðxÞ and ~sðxÞ are polynomials in x(s)of 1 degree at most two and one, respectively, and k is a con- qðsÞ Dxs > 0; a 6 s 6 b 1 ð11Þ stant. It is convenient to rewrite (3) in the following equiv- 2 ´ alent form (Nikiforov and Uvarov, 1988; Alvarez-Nodarse supported in a countable set of the real line (a,b) and such and Smirnov, 1996) that D ryðsÞ DyðsÞ D rðsÞ þ sðsÞ þ kyðsÞ¼0 ð5Þ ½rðxÞqðxÞ ¼ sðxÞqðxÞð12Þ Dxs 1 rxðsÞ DxðsÞ 1 2 Dxs 2 where Some important discrete orthogonal polynomials on the non-uniform lattice are listed in Table 1. Among them, the 1 1 rðsÞ¼r~ðxðsÞÞ ~sðxðsÞÞ Dxs ð6Þ dual Hahn polynomials and Racah polynomials are rela- 2 2 tively simple in terms of the lattice form and weighting sðsÞ¼~sðxðsÞÞ ð7Þ function, moreover, both polynomials have a finite domain of definition that is suited for square images of size N · N The polynomial solutions of equation (5), denoted by pixels. We choose here the dual Hahn polynomials to y (x(s)) P (s), are uniquely determined, up to a normaliz- n n define a new type of moments. The Racah polynomials ing factor B , by the difference analogue of the Rodrigues n have already been used by the authors in another paper formula (Nikiforov and Uvarov, 1988) (Zhu et al., 2007). A comparison of these two moments is B P ðsÞ¼ n rðnÞ½q ðsÞ ; given in Table 2. n qðsÞ n n
ðnÞ r r r 2.2. Dual Hahn polynomials rn ½qnðsÞ ¼ ½qnðsÞ ð8Þ rx1ðsÞ rxn 1ðsÞ rxnðsÞ ðcÞ The classical dual Hahn polynomials wn ðs; a; bÞ, where n =0,1,...,N 1, defined on a non-uniform lattice n Yn x(s)=s(s + 1), are solutions of (5) corresponding to (Nik- x ðsÞ¼xsþ ; q ðsÞ¼qðn þ sÞ rðs þ kÞð9Þ n 2 n iforov and Uvarov, 1988) k¼1 rðsÞ¼ðs aÞðs þ bÞðs cÞ It is known that for some special kind of lattices, solutions s s ab ac bc a b c 1 x s 13 of (5) are orthogonal polynomials of a discrete variable, ð Þ¼ þ þ ð Þ ð Þ i.e., they satisfy the following orthogonality property (Nik- k ¼ n iforov and Uvarov, 1988) and the weighting function q(s) is given by Xb 1 1 2 Cða þ s þ 1ÞCðc þ s þ 1Þ P nðsÞP mðsÞqðsÞ Dxs ¼ dnmdn ð10Þ qðsÞ¼ ð14Þ s¼a 2 Cðs a þ 1ÞCðb sÞCðb þ s þ 1ÞCðs c þ 1Þ
Table 1 Some important discrete orthogonal polynomials on the non-uniform lattice (p, q, a, b, c, c, l, a, b, and x are parameters attached to the respective polynomials, n denoting the order)
Name Notation Lattice smin smax q(s) ðcÞ Cðaþsþ1ÞCðcþsþ1Þ Dual Hahn wn ðs; a; bÞ x(s)=s(s +1) ab Cðs aþ1ÞCðb sÞCðbþsþ1ÞCðs cþ1Þ ða;bÞ Cðaþsþ1ÞCðs aþbþ1ÞCðbþa sÞCðbþaþsþ1Þ Racah un ðs; a; bÞ x(s)=s(s+1) ab Cða bþsþ1ÞCðs aþ1ÞCðb sÞCðbþsþ1Þ sðs 1Þ s q-Krawtchouk kðpÞðs; N; qÞ x(s)=q2s 0 N q l ½N ! ; l ¼ p=ð1 pÞ n Cqðsþ1ÞCqðNþ1 sÞ s q-Meixner mðr;lÞðs; qÞ x(s)=q2s 0 1 l CqðrþsÞ n CqðrÞCqðsþ1Þ sðsþ1Þ ðcÞ q ½sþa q!½sþc q! q Hahn w ðs; a; bÞ x(s)=[s]q[s +1]q ab n q ½s a q!½s c q!½sþb q!½b s 1 q!
Table 2 Comparison of dual Hahn moments and Racah moments Moment Polynomials Symmetry Number of Reconstruction Compression parameters capability capability
Dual Hahn Relatively simple: 2F3 Orthogonal About n Three: (a,b,c) Very good Good
Racah More complex: 3F4 Orthogonal About n + a and s only if Four: (a,b,a,b) Good Very good a = a = b H. Zhu et al. / Pattern Recognition Letters 28 (2007) 1688–1704 1691 where the parameters a, b and c are restricted to In this case, the orthogonality condition given by (19) 1 becomes < a < b; jcj < 1 þ a; b ¼ a þ N ð15Þ Xb 1 2 ðcÞ ðcÞ w^ n ðs; a; bÞw^ m ðs; a; bÞ¼dnm; n; m ¼ 0; 1; ...; N 1 Note that if the uniform lattice, i.e., x(s)=s, is used in s¼a (5), and the parameters a, b, and c are defined as ð22Þ a =(a + b)/2, b = a + N, c =(b a)/2, the dual Hahn ða;bÞ The values of the weighted dual Hahn polynomials are polynomials become the Hahn polynomials hn ðx; NÞ (Nikiforov and Uvarov, 1988). Setting a = 0 and b =0, thus confined within the range of [ 1,1]. Fig. 1 shows the the Hahn polynomials reduce to the Tchebichef polynomi- plots for the first few orders of the weighted dual Hahn polynomials with the parameters a = c = 0 and b = N = 40. als. If we take b = pt and a =(1 p)t in the Hahn polyno- mials and let t !1, we obtain the Krawtchouk polynomials Kn(x;p,N)(Koekoek and Swarttouw, 1998). 2.4. Dual Hahn moments The nth order dual Hahn polynomials are defined as (A´ lvarez-Nodarse and Smirnov, 1996) The dual Hahn moments are a set of moments formed ða b þ 1Þ ða þ c þ 1Þ by using the weighted dual Hahn polynomials. Given a uni- wðcÞðs; a; bÞ¼ n n f s t N N n n! form pixel lattices image ( , ) with size · , the (n + m)th order dual Hahn moment is defined as 3F 2ð n; a s; a þ s þ 1; a b þ 1; a þ c þ 1; 1Þ Xb 1 Xb 1 n ¼ 0; 1; ...; N 1; s ¼ a; a þ 1; ...; b 1 ð16Þ ðcÞ ðcÞ W nm ¼ w^ n ðs; a; bÞw^m ðt; a; bÞf ðs; tÞ; s¼a t¼a where (u)k is the Pochhammer symbol defined as n; m ¼ 0; 1; ...; N 1 ð23Þ Cðu þ kÞ ðuÞ ¼ uðu þ 1Þ ðu þ k 1Þ¼ ð17Þ k CðuÞ The orthogonality property of the dual Hahn polynomials helps in expressing the image intensity function f(s, t) in terms and 3F2(Æ) is the generalized hypergeometric function given of its dual Hahn moments. The reconstructed image can be by obtained by using the following inverse moment transform. X1 k ða1Þ ða2Þ ða3Þ z XN 1 XN 1 F ða ; a ; a ; b ; b ; zÞ¼ k k k ð18Þ 3 2 1 2 3 1 2 b b k! f s; t W w^ðcÞ s; a; b w^ðcÞ t; a; b ; k¼0 ð 1Þkð 2Þk ð Þ¼ nm n ð Þ m ð Þ n¼0 m¼0 The dual Hahn polynomials satisfy the following ortho- s; t ¼ a; a þ 1; ...; b 1 ð24Þ gonality property: Xb 1 1 In (24), s and t represent horizontal and vertical directions wðcÞðs; a; bÞwðcÞðs; a; bÞqðsÞ Dxs ¼ d d2; of reconstructed image with uniform pixel grid. If only the n m 2 nm n s¼a dual Hahn moments of order up to M are used, (24) is n; m ¼ 0; 1; ...; N 1 ð19Þ approximated by where the weighting function q(s) is given by (14) and XM XM ðcÞ ðcÞ f ðs; tÞ¼ W nmw^n ðs; a; bÞw^m ðt; a; bÞð25Þ 2 Cða þ c þ n þ 1Þ d ¼ ; n¼0 m¼0 n n!ðb a n 1Þ!Cðb c nÞ n ¼ 0; 1; ...; N 1 ð20Þ 0.4
The set of dual Hahn polynomials is not suitable for 0.3 defining the moments because the range of values of the polynomials expands rapidly with the order. To surmount 0.2 this shortcoming, we introduce the weighted dual Hahn polynomials in the following subsection. 0.1 w 2.3. Weighted dual Hahn polynomials 0
To avoid numerical instability in polynomial computa- –0.1 tion, the dual Hahn polynomials are scaled by utilizing n=0 –0.2 n=1 the square norm and the weighting function. The set of n=2 the weighted dual Hahn polynomials is defined as n=3 sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi n=4 –0.3 5 10 15 20 25 30 35 40 ðcÞ ðcÞ qðsÞ 1 w^n ðs; a; bÞ¼wn ðs; a; bÞ 2 Dxs ; s dn 2 ðcÞ Fig. 1. Plot of scaled of dual Hahn polynomials w ¼ w^n ðs; a; bÞ for n ¼ 0; 1; ...; N 1 ð21Þ N = 40 with a = c = 0, and b = 40. 1692 H. Zhu et al. / Pattern Recognition Letters 28 (2007) 1688–1704
3. Computational aspects of dual Hahn moments Eq. (26) is used to compute the values of the weighted dual Hahn polynomials. The algorithm for computing the dual ðcÞ In this section, we discuss the computational aspects of Hahn polynomial values w^n ðs; a; bÞ is shown in Fig. 2. dual Hahn moments. We present some properties of dual Fig. 3 shows the reconstruction algorithm based on Eq. (24). Hahn polynomials and show how they can be used to facil- itate the computation of moments. 3.2. Recurrence relation with respect to s 3.1. Recurrence relation with respect to n The recurrence relation of discrete dual Hahn polynomi- als with respect to s is as follows (The detailed derivation is In order to decrease the computational cost in the calcu- given in Appendix B): lation of moments, the recurrence relation can be used to ð2s 1Þ½rðs 1Þþðs 1Þsðs 1Þ 2k sðs 1Þ avoid the overflowing for mathematical functions like the wðcÞðs; a; bÞ¼ wðcÞðs 1; a; bÞ hypergeometric and gamma functions. The weighted dual n ðs 1Þ½rðs 1Þþð2s 1Þsðs 1Þ n s rðs 1Þ Hahn polynomials obey the following recurrence relation wðcÞðs 2; a; bÞ ð31Þ ðs 1Þ½rðs 1Þþð2s 1Þsðs 1Þ n (the derivation of the relation is given in Appendix A)
ðcÞ dn 1 ðcÞ dn 2 ðcÞ To obtain the values for s = 0 and s = 1, we derive w^n ðs; a; bÞ¼A w^n 1ðs; a; bÞþB w^ n 2ðs; a; bÞð26Þ ðnÞ dn dn rn ½qnðsÞ defined by (8) for dual Hahn polynomials as fol- lows (A´ lvarez-Nodarse et al., 2001a): where Xn l n! 1 2 rðnÞ½q ðsÞ ¼ ð 1Þ A ¼ ½sðsþ 1Þ ab þ ac bc ðb a c 1Þð2n 1Þþ2ðn 1Þ n n l! n l ! n l¼0 ð Þ ð27Þ rx ðs l þ 1=2Þ Q n q ðs lÞ 1 n x s m l 1 =2 n B ¼ ða þ c þ n 1Þðb a n þ 1Þðb c n þ 1Þð28Þ m¼0 r nð ð þ Þ Þ n ð32Þ with sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi Using (32) and x(s)=s(s + 1), we have qðsÞ 1 w^ðcÞðs;a;bÞ¼ Dxs ð29Þ 0 2 2 1 nþ1 d0 sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi rxn rx rðnÞq ð0Þ¼Q 2 q ð0Þ¼Q 2 1 q ðsÞ q ðs 1Þ qðsÞ 1 n n n rx ð m 1Þ n n rx n mþ1 w^ðcÞðs;a;bÞ¼ 1 1 Dxs m¼0 n 2 m¼0 2 1 qðsÞ xðs þ 1=2Þ xðs 1=2Þ 2 2 d1 Q n þ 1 qnð0Þ ¼ n qnð0Þ¼ ð33Þ ð30Þ m¼0ðn m þ 1Þ n!
Fig. 2. Algorithm for computing the weighted dual Hahn polynomial values. H. Zhu et al. / Pattern Recognition Letters 28 (2007) 1688–1704 1693
Fig. 3. Algorithm for reconstruction of the original image using Eq. (24).
ðnÞ Q n þ 3 Q nðn þ 1Þ Similarly, we can obtain the recurrence relation for the rn qnð1Þ¼ n qnð1Þ n qnð0Þ m¼0ðn þ 3 mÞ m¼0ðn þ 2 mÞ weighted dual Hahn polynomials with respect to s 2 nðn þ 1Þ c ð2s 1Þ½rðs 1Þþðs 1Þsðs 1Þ 2k sðs 1Þ ¼ qnð1Þ qnð0Þ ð34Þ w^ ð Þðs; a; bÞ¼ ðn þ 2Þ! ðn þ 2Þ! n ðs 1Þ½rðs 1Þþð2s 1Þsðs 1Þ sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiÂÃ we deduce from (8), (33), and (34) that qðsÞ Dxs 1 2 w^ ðcÞðs 1; a; bÞ qðs 1Þð2s 1Þ n ðcÞ 1 wn ð0;a;bÞ¼ ða þ 1Þ ðc þ 1Þ ðb nÞ ð35Þ ðn!Þ2 n n n s rðs 1Þ ðs 1Þ½rðs 1Þþð2s 1Þsðs 1Þ 2 qð0Þ qnð1Þ nðn þ 1Þ sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi wðcÞð1;a;bÞ¼ wðcÞð0;a;bÞ ÂÃ n n 2 n 1 q 1 q 0 2 n 1 ð þ Þð þ Þ ð Þ nð Þ qðsÞ Dxs 2 ðcÞ w^ n ðs 2; a; bÞ; ð36Þ qðs 2Þð2s 3Þ Cða þ s þ n þ 1ÞCðc þ s þ n þ 1Þ n ¼ 0; 1; ...; N 1; s ¼ 2; 3; ...; b 1 q ðsÞ¼ n Cðs a þ 1ÞCðb s nÞCðb þ s þ 1ÞCðs c þ 1Þ ð38Þ ð37Þ with 1694 H. Zhu et al. / Pattern Recognition Letters 28 (2007) 1688–1704
ab0.4 0.5 0.3 0.4
0.3 0.2 0.2 0.1 0.1 w
0 w 0 –0.1 –0.1 –0.2
–0.3 n=0 –0.2 n=0 n=1 n=1 –0.4 n=2 n=2 n=3 –0.3 n=3 –0.5 n=4 n=4 –0.4 5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 50 s s cd0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1
w 0 w 0
–0.1 –0.1
–0.2 –0.2 n=0 n=0 n=1 n=1 n=2 –0.3 –0.3 n=2 n=3 n=3 n=4 n=4 –0.4 –0.4 5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 50 s s
e 0.4
0.3
0.2
0.1
w 0
–0.1
–0.2 n=0 n=1 n=2 –0.3 n=3 n=4 –0.4 5 10 15 20 25 30 35 40 45 50 s
ðcÞ Fig. 4. The influence of parameter c on the weighted dual Hahn polynomials w ¼ w^n ðs; a; bÞ, a =8,b = 48. (a) c = 8, (b) c = 4, (c) c = 0, (d) c = 4, and (e) c =8. sffiffiffiffiffiffiffiffiffi 1 qð0Þ ðcÞ q q n n w^n ð0; a; bÞ¼ 2 ða þ 1Þnðc þ 1Þnðb nÞn 2 ðcÞ 2 ð0Þ nð1Þ ð þ 1Þ n! d w^ n ð1; a; bÞ¼ ð Þ n ðn þ 2Þðn þ 1Þ qð1Þ q ð0Þ 2 1 d sffiffiffiffiffiffiffiffiffiffiffiffi n ¼ ða þ nÞðc þ nÞðb nÞ n 1 w^ ðcÞ ð0; a; bÞ 2 n 1 3qð1Þ n dn w^ðcÞð0; a; bÞð40Þ ð39Þ qð0Þ n H. Zhu et al. / Pattern Recognition Letters 28 (2007) 1688–1704 1695
Table 3 Image reconstruction of the letter ‘‘F’’ of size (40 · 40) without noises, a =8,b =48 Original image (40 · 40)
Reconstructed image Iterative number c= 8 c = 4 c =0 c =4 c =8 1
5
15
39
0.12 The symmetry property is also useful in minimizing the a = 8, b = 48, c = –8 a = 8, b = 48, c = –4 memory requirements for storing the dual Hahn polyno- 0.1 a = 8, b = 48, c = 0 mial values. In fact, Eq. (23) can be rewritten as a = 8, b = 48, c = 4 a = 8, b = 48, c = 8 8 > bP 1 bP 1 0.08 > w^ ðcÞðs; a; bÞw^ ðcÞðt; a; bÞf ðs; tÞ > n m > s¼a t¼a > > if n < N=2 and m < N=2 0.06 > > > bP 1 bP 1 > 1 sw^ ðcÞ s; a; b w^ ðcÞ t; a; b f s; t 0.04 > ð Þ N 1 nð Þ m ð Þ ð Þ Reconstruction Error > > s¼a t¼a <> if n > N=2 and m < N=2 0.02 W nm ¼ > bP 1 bP 1 > c > ð 1Þtw^ðcÞðs; a; bÞw^ ð Þ ðt; a; bÞf ðs; tÞ > n N 1 m 0 > s¼a t¼a 0 5 10 15 20 25 30 35 > > if n < N=2 and m > N=2 Moment Order > > > bP 1 bP 1 Fig. 5. Comparative analysis of reconstruction error of dual Hahn > sþt ðcÞ ðcÞ moment with different coefficients c. > ð 1Þ w^ N 1 nðs; a; bÞw^N 1 mðt; a; bÞf ðs; tÞ > s¼a t¼a :> if n > N=2 and m > N=2 Eqs. (38)–(40) can be used to effectively calculate the ð42Þ weighted dual Hahn polynomial values. Using (41), the inverse transformation (24) can be modified as 3.3. Symmetry property NX=2 1 NX=2 1 ðcÞ ðcÞ The symmetry property can be used to reduce the time f ðs; tÞ¼ w^ n ðs; a; bÞw^m ðt; a; bÞ required for computing the dual Hahn moments. The n¼Â0 m¼0 W þð 1ÞsW þð 1ÞtW weighted dual Hahn polynomials have the symmetry prop- nm N 1 nÃ;m n;N 1 m sþt erty with respect to n þð 1Þ W N 1 n;N 1 m ð43Þ
ðcÞ s ðcÞ where f(s,t) represents image intensity with uniform pixel w^n ðs; a; bÞ¼ð 1Þ w^N 1 nðs; a; bÞð41Þ grid. Such decomposition permits decreasing the computa- The above relation shows that only the values of the dual tional complexity in the reconstruction process. In fact, the Hahn polynomials up to order n = N/2 1 need to be cal- number of multiplication required in (43) is N2/4 assuming culated, thus the computational cost can be reduced. that all the polynomial values have been already calculated, 1696 H. Zhu et al. / Pattern Recognition Letters 28 (2007) 1688–1704
a 0.4
0.3
0.2
0.1 w 0
–0.1 n=0 n=1 –0.2 n=2 n=3 n=4 –0.3 10 20 30 40 50 60 70 80 s
bc0.4 0.4
0.3 0.3
0.2 0.2
0.1 0.1 w w 0 0
–0.1 –0.1 n=0 n=0 n=1 n=1 –0.2 n=2 –0.2 n=2 n=3 n=3 n=4 n=4 –0.3 –0.3 10 20 30 40 50 60 70 80 10 20 30 40 50 60 70 80 s s
ðcÞ Fig. 6. Plot of weighted dual Hahn polynomials w ¼ w^n ðs; a; bÞ for different choices parameters. (a) a = c = 0 and b = 60; (b) a = c = 7, and b = 67; (c) a = c = 18, and b = 78. and the number of addition in (43) is 3N2/4. On the other XN 1 XN 1 c p hand, the multiplication number and addition number in vpq ¼ m00 ½ðs sÞ cos h þðt tÞ sin h (24) are both N2. s¼0 t¼0 ½ðt tÞ cos h ðs sÞ sin h qf ðs; tÞð45Þ 3.4. Invariant pattern recognition using geometric moments where p þ q To obtain the translation, scale and rotation invariants c ¼ þ 1 ð46Þ 2 of dual Hahn moments, we adopt the same strategy used by Yap et al. (2003) for Krawtchouk moments. That is, m10 m01 s ¼ ; t ¼ ð47Þ we derive the dual Hahn moment invariants through the m00 m00 geometric moments. If the geometric moments of an image 2l f(s,t) are expressed using the discrete sum approximation h ¼ 0:5 tan 1 11 ð48Þ as l20 u02
and lpq are the central moments defined by XN 1 XN 1 Z Z m ¼ sptqf ðs; tÞð44Þ 1 1 pq p q s¼0 t¼0 lpq ¼ ðs sÞ ðt tÞ f ðs; tÞdsdt ð49Þ 1 1 then the set of geometric moment invariants, which are The dual Hahn moments of f~ðs; tÞ¼½qðsÞqðtÞð2s þ 1Þ independent to rotation, scaling and translation can be ð2t þ 1Þ 1=2f ðs; tÞ can be written according to the geomet- written as (Hu, 1962) ric moments as H. Zhu et al. / Pattern Recognition Letters 28 (2007) 1688–1704 1697
Table 4 Image reconstruction of the Chinese character of size 60 · 60 without noises Original image (60 · 60)
Reconstructed image Iterative number a = c =0,b =60 a = c =7,b =67 a = c = 18, b =78 10
20
30
50