Retracted: Gabor Weber Local Descriptor for Bovine Iris Recognition

Iris recognition is a robust biometric technology. This paper proposes a novel local descriptor for bovine iris recognition, named Gabor Weber local descriptor (GWLD). We first compute the Gabor magnitude maps for the input bovine iris image, and then calculate the differential excitation and orientation for each pixel over each Gabor magnitude map. After that, we use these differential excitations and orientations to construct the GWLD histogram representation. Finally, histogram intersection is adopted to measure the similarity between different GWLD histograms. The experimental results on the SEU bovine iris database verify the representation power of our proposed local descriptor.


Introduction
The food safety and security has become one of the hot topics noticed widely by the society in recent years due to its harm to the health.Traceability system is considered to be an effective solution to meet heightened consumer expectations.Identification technology is crucial for the traceability system of meat supply chain.In existing traceability systems, different identification technologies have been adopted, such as ear tag [1], tattoo, and RFID [2].Although these technologies have achieved some results, they still have deficiency to support the traceability system, as these technologies rely on manual code.Recently, a tracking and traceability system for meat supply chain based on iris recognition was proposed by Zhao et al. [3] to ensure food safety.
Iris recognition has gained popularity as a robust and reliable biometric technology for its characteristics of uniqueness, stability, and noninvasiveness of iris [4].The iris's complex texture and its apparent stability hold tremendous promise for leveraging iris recognition in diverse application scenarios, such as border control, forensic investigations, and cryptosystems [5].One of the key procedures for iris recognition systems is feature representation.Feature representation aims to describe the unique iris texture pattern, which should be highly distinctive among iris images of different eyes while maintaining invariance among iris images of the same eyes [6].
The most common practice is to generate the feature representation by extracting the spatial frequency information from the iris texture with different filters.The wellstudied filter is 2D Gabor filter, proposed by Daugman [4], as the Gabor filter has shown to be optimal in the sense of minimizing the joint two-dimensional uncertainty in space and frequency.Some other researchers also considered the application of Gabor filter, its variant, or other filters to extract texture information.Considering the ordinal Gabor filter would under represent the high frequency components in natural images, Yao et al. [7] used modified Log-Gabor filters to generate iris code.Wildes [8] applied Laplacian of Gaussian filters to encode the iris texture.Ma et al. [9] used a variant of the Gabor filter at two scales to analyze the iris texture, and generated the feature vector with Fisher's linear discriminant.Sun et al. [10] used a Gaussian filter to compute the local orientation at each pixel in the unwrapped template and then used the binary code quantized from the angle to represent iris image.Monro et al. [11] applied discrete cosine transform (DCT) to generate a binary code of iris texture by calculating the differences between the DCT coefficients.Miyazawa et al. [12] proposed to use Fourier transform to compute the representation of iris texture.
In the past years, local descriptors, which represent the image using local characteristic directly in the spatial domain, have been proven to be successfully in image retrieval, image classification, and object recognition.Several popular local descriptors have been introduced to iris recognition, such as local binary pattern (LBP) [13], scale invariant feature transform (SIFT) [14], and DAISY [15].Unlike the spatial frequency methods, local descriptors extract image properties directly in the spatial domain.Sun et al. [16] proposed an efficient graph-matching scheme for iris recognition, by using LBP to generate the histogram representation of iris texture.Belcher and Du [17] applied the SIFT to the iris recognition.They detected keypoints in the iris images, and represented these keypoints using descriptor of orientation histogram.Zhang et al. [18] designed the Deformable DAISY Matcher (DDM) for iris recognition, based on the DAISY descriptor.
Recently, several literatures suggested to combine local descriptor with Gabor filter to improve its performance, as the Gabor filter has shown its high discriminant ability.Zhang et al. [19] proposed the so-called local Gabor binary pattern histogram sequence (LGBPHS) for face representation, by combining Gabor filter with LBP, and Takahashi et al. [20] applied spatiotemporal Gabor filtering to bicycle detection together with Histogram of Oriented Gradients (HOG) [21].
While the texture of bovine iris images is not as salient as human iris images and there exist noise and illumination in the bovine iris images, it is difficult for traditional feature representation algorithms to achieve high performance, which further affects the application of bovine iris recognition technology in food safety.Therefore, a more accurate bovine iris representation is required.Motivated by Weber local descriptor (WLD) [22] and the combination of Gabor filter and local descriptor [20,21], we propose a novel local descriptor for bovine iris recognition, named Gabor Weber local descriptor (GWLD).The novelty of the GWLD lies in generating the histogram representation by computing the differential excitation and gradient orientation over these Gabor magnitude maps of the bovine iris image.The combination of Gabor and WLD enhances the representation power of the spatial histogram for bovine iris recognition greatly.
The remaining part of the paper is organized as follows.Section 2 introduces the preprocessing of bovine iris images.Section 3 describes the computation scheme of GWLD and bovine iris matching algorithm.And experimental results are shown in Section 4. At last, the conclusion is given in Section 5.

Bovine Iris Image Preprocessing
Similar to human iris recognition, bovine iris recognition usually involves several crucial steps: image acquisition, image preprocessing, feature extraction, and iris matching.Image preprocessing includes iris localization, iris normalization, and enhancement.In this paper, we use the image preprocessing procedure mentioned in our previous work [23].The preprocessing method is briefly described as follows.Firstly, the inner and outer iris edge curves are detected by region-based active contour model, which uses the global and local intensity information to guide the contour evolution.Secondly, the inner and outer iris edge curves are used to fit ellipses with least square method, and then elliptic iris boundaries are obtained.Thirdly, to eliminate the influence of eyelid and lash, which usually exists in the original bovine iris image, [−/4, /4] and [3/4, 5/4] are selected to generate the normalized rectangle iris block with polar coordinate transformation.Finally, the normalized bovine iris block is enhanced with histogram equalization.
Figure 1(a) shows the original bovine iris image with initial contour, Figure 1(b) illustrates the final evolution curves after 50 iterations by using the region-based active contour model, and Figures 1(c), 1(d), and 1(e) show the elliptic iris boundaries, normalized bovine iris image, and enhanced iris block, respectively.

Bovine Iris Representation and Matching
We first describe how to compute the Gabor Weber local descriptor in Section 3.1 and then introduce the bovine irismatching algorithm using GWLD.[22], inspired by the Weber's Law, was proposed as a simple and powerful image descriptor and considered to be able to combine the advantages of LBP [13] in computing the intensity difference pixel-wise and SIFT [14] in using gradient information.The WLD describes the intensity relationship between a pixel and its neighborhoods with two components: differential excitation and orientation.

Computation of GWLD. Weber local descriptor
The overall framework of the proposed bovine iris representation approach based on GWLD is illustrated in Figure 2. In this approach, the bovine iris image is modeled as a GWLD histogram by the following procedures: (a) an input bovine iris image is transformed to obtain multiple Gabor magnitude maps in frequency domain by applying multiscale and multiorientation Gabor filters; (b) compute the differential excitation and orientation for each pixel over each Gabor magnitude map; (c) compute the 2D GWLD histogram for bovine iris image with these differential excitations and orientations of each pixel of each Gabor magnitude map, and then generate the more discriminative 1D GWLD histogram.The following subsections will describe these procedures in detail.

Gabor Magnitude Maps.
Considering the advantages of the Gabor filter in object recognition [4,24], we exploit the multiscale and multiorientation Gabor filters to de-composite the input bovine iris images for feature extraction.The Gabor representation of a bovine iris image is derived by convolving the iris image with the Gabor filters.Let (, ) be the enhanced bovine iris image; its convolution with a Gabor filter is  ,] () defined as where * denotes the convolution operator.The Gabor filters we used are defined as follows: where  and ] denote the orientation and scale of the Gabor filters, respectively.‖ ⋅ ‖ defines the norm operator, and the wave vector is defined as where  ] =  max / V and   = /8. max is the maximum frequency, and  is the spacing factor between different filters in the frequency domain.

Differential Excitation and Orientation.
Differential excitation and orientation are the two components to generate the WLD histogram of input image.Notice that there exist relatively flat areas in the input normalized bovine iris image.We apply the modified WLD operator [25] over each Gabor magnitude map to compute the differential excitation and orientation.
Differential excitation (  ) represents the intensity difference between a current pixel   and its neighbors.To compute (  ), we first calculate the difference V 00  between a current pixel and its neighbors using the following filter  00 : Here   is the corresponding neighbor pixel of the Chebyshev distance where   denotes the th neighbor of   and  is the number of neighbors.Then, we compute the ratio of the difference V 00  to the intensity of the current pixel V 01  , which is the output of filter  01 (as shown in Figure 3): Finally, the differential excitation of the current pixel   is computed as Furthermore, we use V 10  and V 11  , which are the outputs of filters  10   and  11  (as shown in Figure 3), to compute the orientation (  ) for each pixel   : where V 10  =  5 −  1 and V 11  =  7 −  3 .the quantized differential excitation along with the dominant orientation [22], and the length of the 1D GWLD histogram representation is

GWLD Histogram
3.2.Iris Matching.Several similarity measurement approaches have been presented for histogram matching [26].In this paper, we use the histogram intersection as the similarity measurement of two GWLD histogram representations.Suppose GWLD  and GWLD  are two GWLD histogram representations, and the histogram intersection is defined as where  is the length of the histogram.The intuitive motivation for the histogram intersection measurement is the calculation of the common part of two histograms.It is clear that larger value of the histogram intersection means higher similarity of the two input bovine iris images.

Experiments
We apply the proposed descriptor to bovine iris to test its performance.We select 18 subjects, 90 original images from  the SEU bovine iris database.The size of the original iris images is 320 × 240, and the size of normalized iris image is 253 × 61.For each input bovine iris image, there are 40 Gabor magnitude maps after convolution.Therefore, 40 × 253 × 61 pairs of quantized differential excitations and orientations are calculated to generate the 1D GWLD histogram.In addition, accurate recognition rate (ARR) is adopted to evaluate the performance of these experiments.

Parameters Assessment.
According to these applications of WLD in iris recognition [25], face recognition [27], and pedestrian detection [28], the length of histogram representations affects the performance of the local descriptor.Furthermore, according to the aforementioned analysis, the first step of the proposed local descriptor GWLD is to generate these Gabor magnitude maps by convolving the input image with the Gabor filters.As shown in Figure 4, the results of different convolution mask sizes are quite different.Therefore, the mask size may affect the performance of GWLD.
Figure 5 illustrates the performance variation against the length of GWLD histogram and the convolution mask size.The ARR increases considerably with the rise in the length of GWLD histogram and reaches the maximum at  = 32000.Nevertheless, it decreases when the length rises further.That is because the GWLD histogram gets much sparser while the number of the differential excitations and orientations is limited.In additional, large length  increases the computation requirements of quantization of differential excitations and orientations.

Performance Comparison.
Feature descriptor plays an important role in the iris recognition system.In this paper, we compare our proposed GWLD with LBP, center-epsilon LBP, WLD, modified WLD and Gabor filter on the same dataset to evaluate the performance of our proposed descriptor.As can be seen from Table 1, the LBP achieves the lowest accurate recognition rate, and the center-epsilon LBP gets a better performance (while epsilon equals 3), as it takes the central pixel information into consideration to generate the feature representation [29].Meanwhile, the WLD operator shows a more discriminative ability than LBP, for its representation is generated from the intensity difference and the orientation information.And by applying a different computation scheme for the difference excitation, the modified WLD achieves a higher ARR [25].Furthermore, as our proposed GWLD integrates the advantages of Gabor filter and modified WLD, it outperforms other feature descriptors with the highest ARR.In other words, our proposed feature descriptor is quite effective for bovine iris recognition.

Conclusion and Future Work
This paper has proposed a novel feature descriptor GWLD for bovine iris recognition.The GWLD constructs the histogram representation by computing the differential excitation and orientation over the Gabor magnitude maps of the input bovine iris image.Experiments on the SEU bovine iris database show that the proposed feature descriptor has encouraging performance, and the length of the GWLD histogram and the mask size of Gabor filters affect the recognition accuracy.One problem of our approach might be that its computation cost is high, because the differential excitations and orientations are calculated over each Gabor magnitude map for each pixel.Hence, one of our future efforts is to reduce the computational cost with some feature selection strategies.Other problems might be its relatively high-dimensional histogram representation and the similarity measurement.The high dimensional histogram representation can be alleviated by applying some analysis approach (e.g., PCA) to the histogram, and its performance can be improved with more suitable similarity measurement (e.g., ADTW [30]).Furthermore, the applicability of the bovine iris recognition in food safety might be increased by other enhancements, for example, taking the motion images into consideration [31].
(a) Input bovine iris image (b) Gabor magnitude maps with mask size 5 × 5 (c) Gabor magnitude maps with mask size 31 × 31

Figure 4 :
Figure 4: Different convolution masks of the Gabor filter.
[22]esentation.After calculating the differential excitation and orientation, we get the WLD features {, } for each pixel of these 40 Gabor magnitude maps.Note that both ranges of differential excitation and orientation are [−/2, /2].Same to WLD[22], the orientation of each pixel  is mapping to the range [0, 2], and then it is quantized to  dominant orientations Φ  .Then, we compute the 2D GWLD histogram with the differential excitation and the quantized orientation of each pixel of 40 Gabor magnitude maps.And in the 2D histogram, each column corresponds to a dominant orientation Φ  .After computing the 2D histogram, we further quantize the differential excitation .The range of differential excitation  is divided into  subranges, and each subrange includes  dominant differential excitations.Therefore, the 2D GWLD histogram can be encoded into a 1D histogram by regrouping

Table 1 :
Performance comparison among different feature descriptors.