1 August 2005 Efficient Fourier shape descriptor for industrial defect images using wavelets
Author Affiliations +
Optical Engineering, 44(8), 080503 (2005). doi:10.1117/1.1993687
The use of image retrieval and classification has several applications in industrial imaging systems, which typically use large image archives. In these applications, the matter of computational efficiency is essential and therefore compact visual descriptors are necessary to describe image content. A novel approach to contour-based shape description using wavelet transform combined with Fourier transform is presented. The proposed method outperforms ordinary Fourier descriptors in the retrieval of complicated industrial shapes without increasing descriptor dimensionality.
Kunttu, Lepistö, and Visa: Efficient Fourier shape descriptor for industrial defect images using wavelets



The recognition and classification of objects based on their visual similarity has become a central task in current industrial imaging systems. With increasing amounts of real-world image data to be processed and stored, the development of powerful retrieval tools also has become necessary in machine vision applications. Along with texture and color, shape is an essential feature used to describe the objects in the images. Therefore, effective shape description is essential in retrieval systems.

Due to the increasing number of on-line solutions, computational lightness is nowadays considered equally important as classification accuracy. In retrieval, computational efficiency of a particular descriptor is generally dependent on two matters, descriptor dimensionality and matching procedure.

The Fourier descriptor (FD)1 is probably the best-known boundary-based shape descriptor. It has been proven to outperform most other boundary-based methods in terms of retrieval accuracy and efficiency.2 In addition to good retrieval and classification performance, the main advantages of FDs are that (1) they are compact and computationally light, (2) they are easy to implement, (3) their matching is straightforward, (4) they are very easy to normalize to be scale and rotation invariant, and (5) their sensitivity to noise is low.

Wavelet transforms3 have been widely used in multiscale image analysis and also have a few applications in shape description. In Ref. 4, the wavelet descriptors (WDs) are based on zero-crossing points of wavelet approximation of the shape and hence the similarity measurement is dependent on the shape complexity. In Ref. 5, moment invariants are employed in shape description using wavelets. It is also possible to combine wavelets with Fourier descriptors, which yields to rotation and scale invariance. This can be made based on polar coordinates of a shape6 or by Fourier transforming the wavelet coefficients obtained from the complex-valued boundary function.7 On the other hand, when WDs are formed using several scales, the resulting feature vector is typically high dimensional due to spatial information caused by multiple scales.

In this paper, we present an effective approach to wavelet-based shape representation at single scale. We show that it is possible to form rotation and translational invariant WDs, whose matching is as simple and fast as that of FDs. The proposed approach is applied to a practical industrial image retrieval and classification problem.


Shape Description

The contour-based shape description is based on one-dimensional boundary function (shape signature). Let (xk,yk) , k=0,1,2,,N1 represent the object boundary coordinates, in which N is the boundary length. Complex coordinate function z(k) (Ref. 2) expresses the boundary points in an object centered coordinate system:


in which (xc,yc) is the object centroid.


Fourier Descriptors

Fourier descriptors can be formed for the boundary function z(k) using the discrete Fourier transform (DFT):


for n=0,1,2,,N1 and Fn are the transform coefficients of z(k) . The descriptors can be made rotation invariant using the magnitudes of the transform coefficients, Fn . The scale can be normalized by dividing the magnitudes of the coefficients by F1 .

The general shape of the object is represented by the low-frequency coefficients, which are usually selected to be the descriptor. In the contour Fourier method,2 the feature vector of length L is formed as:




Wavelet Shape Descriptor Using Fourier Transform

In the wavelet-based approach, the boundary function z(k) is transformed using some wavelet Ψ .3 The complex wavelet transform8 is based on the continuous wavelet transform (CWT). The CWT of the boundary z(k) is defined as:


In the case of CWT, a set of coefficients Ca(b) of scale a are obtained. The coefficients are defined for all positions b=0,1,2,,N1 .

The problem with the CWT coefficients is that they are dependent on the starting point of the object boundary. Hence, the obtained descriptor is not rotation invariant. Also the dimensionality of the feature vector depends on the boundary length. Therefore, the coefficient vectors of different shapes cannot be directly matched. The proposed solution for this problem is to apply the Fourier transform to the whole set of wavelet coefficients. This way the normalization and matching are straightforward operations. The proposed descriptor is formed by applying the DFT to the coefficients Ca(b) :


In this paper, we use the wavelet shape descriptor at a single scale to keep the descriptor dimensionality as low as in the case of Fourier descriptors. The feature vector of this new descriptor is equal to that of the contour Fourier descriptor presented in Eq. 3.


Experiments with Industrial Defect Shapes

The validation presented in this section is twofold. Simple classification experiments are first carried out to show the influence of scale selection on the shape description. The second part of the validation, the retrieval accuracy of the proposed methods, is compared to that of an ordinary FD (contour Fourier). In all the experiments, Euclidean distance and the “leave one out” validation principle are used.


Testing Database

For testing purposes, we use defect images that are collected from an industrial process using a paper inspection system.9 A reason for collecting defect image databases in process industry is a practical need for controlling the quality of production.9 When retrieving images from a database, the defect shape is one essential property describing the defect class. Therefore, effective methods for the shape representation are necessary. The test set consisted of 1204 paper defect shapes, which represented 14 defect classes with each class consisting of 27–103 images (Fig. 1).

Fig. 1

Example contours of each 14-paper defect class in the testing database.



Classification and Retrieval

The feature extraction in the testing database was carried out by calculating the descriptors for the images in the database. The dimensionality (L) was 8 with all the descriptors [Eq. 3]. In the case of the wavelet-based approach, the selected wavelets ψ were first and second order complex Gaussian wavelets that have been implemented in the Matlab wavelet toolbox.8 To compare different scales, we made preliminary k -nearest neighbor ( k -NN) classification experiments. Figure 2a presents the average classification rates of the proposed wavelet descriptors at different scales using a 5-NN classifier. In this figure, the classification rate of the contour Fourier descriptor (41.87%) is also presented. The scales that produce the highest classification rates were compared to contour Fourier in the retrieval experiment by calculating average precision versus recall curves for the queries [Fig. 2b].

Fig. 2

(a) The average 5-NN classification rates of the proposed methods using different scales of the first and second order complex Gaussian wavelets. (b) Average precision/recall curves of the queries using proposed descriptors that employ the first and second order complex Gaussian wavelets the scales 16 and 2, respectively.




In this paper, we showed that it is possible to overcome the difficulties with shape description using wavelet coefficients (rotational variance and complicated matching) by Fourier transforming the coefficients. The results of the classification and retrieval experiments reveal that the proposed wavelet-based shape description approach clearly outperforms ordinary FDs in defect shape description. It is also essential to note that the proposed descriptors have the same dimensionality and matching procedure as FDs. The computational cost of the feature extraction is somewhat higher than that of FDs due to the wavelet transform. However, the dimensionality of the descriptors is more essential than the feature extraction time, because in retrieval applications the feature extraction is usually an off-line operation. If the computational efficiency of feature extraction is critical, the cost of wavelet transform can be decreased using the algorithm presented in Ref. 10.


The authors wish to thank ABB Oy (Mr. Juhani Rauhamaa) for the paper defect image database used in the experiments.


1.  C. T. Zahn and R. Z. Roskies, “Fourier descriptors for plane closed curves,” IEEE Trans. Comput.0018-9340 C-21(3), 269–281 (1972). Google Scholar

2.  H. Kauppinen, T. Seppänen, and M. Pietikäinen, “An experimental comparison of autoregressive and Fourier-based descriptors in 2D shape classification,” IEEE Trans. Pattern Anal. Mach. Intell.0162-8828 17(2), 201–207 (1995). Google Scholar

3.  C. K. Chui, An Introduction to Wavelets, Academic Press, San Diego (1992). Google Scholar

4.  Q. M. Tieng and W. W. Boles, “Recognition of 2D object contours using the wavelet transform zero-crossing representation,” IEEE Trans. Pattern Anal. Mach. Intell.0162-8828 19(8), 910–916 (1997). Google Scholar

5.  D. Shen and H. Ip, “Discriminative wavelet shape descriptors for recognition of 2-D patterns,” Pattern Recogn.0031-3203 10.1016/S0031-3203(98)00137-X 32, 151–165 (1999). Google Scholar

6.  G. Chen and T. Bui, “Invariant Fourier-wavelet descriptors for pattern recognition,” Pattern Recogn.0031-3203 32, 1083–1088 (1999). Google Scholar

7.  I. Kunttu, L. Lepistö, J. Rauhamaa, and A. Visa, “Multiscale Fourier descriptor for shape-based image retrieval,” Proc. 17th Intl. Conf. Patt. Recog., Vol. 2, pp. 765–768 (2004). Google Scholar

8.  M. Misiti, Y. Misiti, G. Oppenheim, and J.-M. Poggi, Wavelet Toolbox for Use with Matlab, Mathworks Inc. (2001). Google Scholar

9.  J. Rauhamaa and R. Reinius, “Paper web imaging with advanced defect classification,” Proc. 2002 TAPPI Technology Summit, (2002). Google Scholar

10.  F. Nicolier, O. Laligant, and F. Truchetet, “Discrete wavelet transform implementation in Fourier domain for multidimensional signal,” J. Electron. Imaging1017-9909 10.1117/1.1479701 11(3), 338–346 (2002). Google Scholar

Iivari Kunttu, Leena Lepistö, Ari J.E. Visa, "Efficient Fourier shape descriptor for industrial defect images using wavelets," Optical Engineering 44(8), 080503 (1 August 2005). http://dx.doi.org/10.1117/1.1993687


Shape analysis


Image processing

Wavelet transforms

Continuous wavelet transforms

Feature extraction


Back to Top