# Complex Zernike moments

Posted 2012.05.21 01:56

The Zernike polynomials were first proposed in 1934 by Zernike [22]. Their moment formulation appears to be one of the most popular, outperforming the alternatives [19] (in terms of noise resilience, information redundancy and reconstruction capability). The pseudo-Zernike formulation proposed by Bhatia and Wolf [3] further improved these characteristics. However, here we study the original formulation of these orthogonal invariant moments.

Complex Zernike moments [18] are constructed using a set of complex polynomials which form a complete orthogonal basis set defined on the unit disc . They are expressed asTwo dimensional Zernike moment:

 (49)

where  and defines the order,  is the function being described and  denotes the complex conjugate. While  is an integer (that can be positive or negative) depicting the angular dependence, or rotation, subject to the conditions:
 (50)

and  is true. The Zernike polynomials [22] Zernike polynomial expressed in polar coordinates are:
 (51)

where  are defined over the unit disc,  and  is the orthogonal radial polynomial, defined asOrthogonal radial polynomial:
 (52)

where:
 (53)

where  and it must be noted that if the conditions in Equation 1.50 are not met, then . The first six orthogonal radial polynomials are:
 (54)

Figure 1.6 shows eight such radial responses, where it can been seen that the polynomials become more grouped, as they approach the edge of the unit disc ( approaches unity). (Care must be taken with regard to the accuracy of these polynomial calculations as the factorial operations can quickly produce large integer values, even at relatively low order ). The difference between these orthogonal polynomials and the non-orthogonal monomials can be seen by comparing Figure 1.6 with Figure 1.2.

So for a discrete image, if  is the current pixel then Equation 1.49 becomes:
 (55)

To calculate the Zernike moments, the image (or region of interest) is first mapped to the unit disc using polar coordinates, where the centre of the image is the origin of the unit disc. Those pixels falling outside the unit disc are not used in the calculation. The coordinates are then described by the length of the vector from the origin to the coordinate point, , and the angle from the  axis to the vector Polar co-ordinate radius, Polar co-ordinate angle, by convention measured from the positive  axis in a counter clockwise direction. The mapping from Cartesian to polar coordinates is:
 (56)

where
 (57)

However,  in practice is often defined over the interval , so care must be taken as to which quadrant the Cartesian coordinates appear in. Translation and scale invariance can be achieved by normalising the image using the Cartesian moments prior to calculation of the Zernike moments [7]. Translation invariance is achieved by moving the origin to the image's COM, causing . Following this, scale invariance is produced by altering each object so that its area (or pixel count for a binary image) is , where  is a predetermined value. Both invariance properties (for a binary image) can be achieved using :
 (58)

and  is the new translated and scaled function. The error involved in the discrete implementation can be reduced by interpolation. If the coordinate calculated by Equation 1.58 does not coincide with an actual grid location, the pixel value associated with it is interpolated from the four surrounding pixels. As a result of the normalisation, the Zernike moments  and are set to known values.  is set to zero, due to the translation of the shape to the centre of the coordinate system. This however will be affected by a discrete implementation where the error in the mapping will decrease as the shape (being mapped) size (or pixel-resolution) increases.  is dependent on , and thus on :
 (59)

Further, the absolute value of a Zernike moment is rotation invariant as reflected in the mapping of the image to the unit disc. The rotation of the shape around the unit disc is expressed as a phase change, if  is the angle of rotation,  is the Zernike moment of the rotated image and  is the Zernike moment of the original image then:
 (60)

Subsections

Next: A note on image Up: Orthogonal moments Previous: Legendre moments
Jamie Shutler 2002-08-15

#### 'Enginius > Machine Learning' 카테고리의 다른 글

 ANN로 headPoseEstimate하기  (0) 2012.06.05 2012.05.31 2012.05.21 2012.05.01 2012.04.26 2012.04.26
« PREV : 1 : ··· : 112 : 113 : 114 : 115 : 116 : 117 : 118 : 119 : 120 : ··· : 143 : NEXT »