# Integrating Globality and Locality for Robust Representation Based Classification

Authors: Zheng Zhang; Zhengming Li; Binglei Xie; Long Wang; Yan Chen
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

The representation based classification method (RBCM) has shown huge potential for face recognition since it first emerged. Linear regression classification (LRC) method and collaborative representation classification (CRC) method are two well-known RBCMs. LRC and CRC exploit training samples of each class and all the training samples to represent the testing sample, respectively, and subsequently conduct classification on the basis of the representation residual. LRC method can be viewed as a “locality representation” method because it just uses the training samples of each class to represent the testing sample and it cannot embody the effectiveness of the “globality representation.” On the contrary, it seems that CRC method cannot own the benefit of locality of the general RBCM. Thus we propose to integrate CRC and LRC to perform more robust representation based classification. The experimental results on benchmark face databases substantially demonstrate that the proposed method achieves high classification accuracy.


---

## Body

## 1. Introduction

Face recognition has become a popular technique as one of the most promising branches of the pattern recognition [1–4]. Recently, the representation based classification method (RBCM) was proposed for face recognition by the specialists and scholars in the fields of computer vision and pattern recognition [5–7]. The RBCM are completely different from the conventional classification methods, such as principle component analysis (PCA) [8] and kernel PCA [9], linear discriminant analysis (LDA) [10], kernel LDA [11], and Gaussian maximum likelihood [12]. The common procedure of these traditional methods is to first utilize all the training samples to obtain the transform axes and then exploit the predetermined transform axes to convert all the samples into a new lower-dimensional space to perform recognition [4, 8–12]. However, the RBCM generally uses the linear combination of the training samples to represent the test sample and then utilizes the representation results to classify the test sample. The RBCM has been widely applied to the face recognition [13–15], image categorization [16, 17], and image superresolution [18, 19].Generally speaking, the representation based classification consists of two categories, the global representation method and the local representation method. Hereafter, the global representation method refers to use all the training samples to represent the test sample, whereas the local representation method just uses the fraction of the training samples to represent the test sample. What is more, from the viewpoint of the image restoration, the linear combination of training samples can be viewed as a new image with the same size of the samples. If the linear combination of training samples from a class can perfectly reconstruct the probe image, the probe image can be classified as a member of this class. The corresponding residual can be regarded as the distance between the test sample and each class. Both of the global and local representation methods make use of the residual to perform final classification.Since the sparse representation based classification (SRC) method [5] was proposed, the sparse representation based classification has attracted tremendous attention and plenty of representation based classification methods were designed to improve the effectiveness and efficiency of SRC [1, 2, 20, 21]. SRC can be somewhat viewed as the global representation method and SRC exploits the linear combination of all the training samples with a sparse optimizer to represent the given test sample. Moreover, the sparse optimizer is usually implemented by the l0-norm, lp-norm (0<p<1), l1-norm, and even l2-norm minimization. Pati et al. [22] proposed the orthogonal matching pursuit (OMP) method exploiting the l0-norm minimization to obtain the sparse representation solution. Xu et al. [23] exploited the l1/2-norm minimization to restrain the sparsity of representation coefficients. The RBCM with l1-norm or l2-norm minimization constraint seems to be the mainstream approaches [24–26]. The literature [27–29] has proposed some powerful l1-norm and l2-norm minimization based algorithms for SRC. For example, Beck and Teboulle [27] proposed a fast iterative shrinkage-thresholding algorithm (FISTA) using Nesterov’s optimal gradient-based algorithm for convex optimization problems. Kim et al. [29] proposed a novel method exploiting conjugate gradients algorithm to address the convex quadratic problem. However, SRC suffers from distinct time complexity on account of solving the minimization problem iteratively. Thus abundant research communities convert their attention to develop RBCM with loose sparse restriction, l2-norm minimization constraint. We briefly review some remarkable literature resources in the following. Zhang et al. [30] presented an efficient representation method called collaborative representation classification (CRC) method. CRC method introduces l2-norm minimization, rather than l1-norm minimization, to acquire limited-sparse representation solution. Furthermore, it is confirmed that the computational efficiency of CRC is extremely higher than SRC on the basis of analogical classification accuracy [30]. Xu et al. [4] put forward a two-step sparse representation method which first utilized CRC method to figure out M nearest samples and then implemented CRC on the predetermined M samples to make final classification. Liu et al. [31] used reconstruction error and normalized distance of RBCM to perform palmprint recognition method. Xu et al. [32] applied the model of RBCM to improve the nearest neighbor classifier. Yang et al. [26] presented a relaxed collaborative representation classification using the similar and distinctive features. Xu et al. [33] proposed to utilize the bimodal biometric traits to perform sparse representation method for biometrics and palmprint classification. The literature [25, 26, 30–34] also certifies the robustness and efficiency of CRC method in contrast to SRC. Furthermore, CRC also endows with the characteristic of the limited sparsity [30, 34].Local recognition methods have been proposed to improve the classification accuracy and have extraordinary priority in comparison with the global recognition methods [35–37]. The local recognition methods only exploit part of the training samples rather than the whole of training samples to obtain intermediate result and then use the result to make classification. For example, Xu et al. [38] proposed a novel and efficient solution scheme for locality preserving projection (LPP) method and this method can effectively conquer the small training sample size problem. Sugiyama [39] proposed a local LDA method that combined the advantages of LDA and the preponderance of LPP method to perform supervised dimensionality reduction. Liu et al. [40] demonstrated that the local principle component analysis (PCA) had great superiorities compared to the global PCA. In order to obtain better face recognition results, local RBCM method recently was proposed to exploit each class of the training samples to represent the test sample and classify the test sample based on the representation residual. Linear regression classification (LRC) method was proposed to utilize each class of the training samples to represent the test sample and classify the test sample into the class which generates the minimum representation residual [35]. Naseem et al. [36] developed a robust regression method using a linear combination of specific class of training sample to represent the test sample. Tahir et al. [37] fused the histogram and face features using LRC classifier to constitute a texture descriptor, named multiscale local phase quantization histogram (MLPQH) feature, to solve the illumination problem for face recognition. The literature [35–37] certified the significance of locality in RBCM and also demonstrated that the local representation methods could greatly lower the recognition error rate.In order to obtain the bilateral dominant characteristics of the global and local RBCM, we integrate CRC and LRC together to perform face recognition. CRC method possesses the property of limited sparsity and can properly harness the nature of low-complex sparsity to perform global representation based classification. The LRC method is a local representation method and is able to employ the local structure of patterns to perform classification. Section4 verifies the rationality of the proposed method and the substantial experiments in Section 5 also confirm its effectiveness.The rest of this paper is organized as follows: we first briefly review the LRC and CRC methods in Section2. In Section 3, we present the proposed optimizing method. In Section 4, we discuss the potential rationale of the proposed method. In Section 5, the proposed method is verified by conducting extensive experiments on several well-known benchmark face databases. Finally, Section 6 concludes the paper.


## 2. Related Work

### 2.1. Linear Regression Classification Algorithm

Suppose that there arec known pattern classes with N training samples x1,x2,…,xN and the query sample is y. All the samples should be firstly normalized to unit column vectors with the length of 1.LetXi be the matrix composed by the training samples of the ith class; that is, Xi=[xi1,xi2,…,xiMi]∈Rd×Mi where d is the dimension of the sample and Mi is the number of training samples of the ith class. Let us note a matrix X=[X1,X2,…,Xc]∈Rd×N and, apparently, the matrix X denotes the whole of training samples.Linear regression classification (LRC) method [35] first exploits a linear combination of each class of training samples to represent the testing sample y; that is, the following equation is satisfied:
(1)y=a1xi1+a2xi2+⋯+anxiMi,
where xi1,xi2,…,xiMi are the ith class of the training samples and aj (j=1,2,…,n) is the coefficient of xij. For convenient presentation, we can rewrite (1) into the following equation:
(2)y=Xiβi,
where Xi=[xi1,xi2,…,xiMi]∈Rd×Mi, βi=[a1,a2,…,an]T, and y∈Rd×1 are column vectors.According to the least square theory,β can be calculated by the equation below:
(3)βi=(XiTXi)-1XiTy,i=1,2,…,c.We refer toy~=Xiβi as the representation result of the ith class of the training samples using preobtained βi and we can transform y~ into a two-dimension image with the same size as the original sample image.Finally, we calculate the representation residuals of each class of the training samples and then exploit the representation residuals to conduct classification. The residual of representing the testing sample,y, using the training samples of the ith class is
(4)resi(i)=∥y-Xiβi∥2.We noteXiβi as the representation contribution using the training samples from the ith class and the RBCM presumes that the test sample, y, can be adequately represented by only the training samples from the same class [41]. Thus, a smaller representation residual can reflect a larger representation contribution, and then the minimum representation residual, resi(i), means the nearest distance between the testing sample and the ith class. Consequently, we classify the testing sample, y, into the class which achieves the smallest representation residual.


### 2.2. Collaborative Representation Classification Algorithm

Collaborative representation classification (CRC) [30] method first utilizes all the training samples to represent the test sample and exploits the linear combination of the whole of training samples to approximate the test sample. That is, the following approximation is satisfied:
(5)y=a1x1+a2x2+⋯+anxN
or
(6)y=XA,
where y is the testing sample. The ai (i=1,2,…,N) is the coefficient of xi and A=[a1a2⋯an]T. The solution of A can be calculated by solving l2-norm minimization problem:
(7)A^=argmin⁡⁡∥A∥2s.t.y=XAor(8)min⁡(∥y-XA∥2+μ∥A∥2),
where μ is a small positive constant and (8) is to utilize the Lagrangian Method to obtain the solution of (7). We solve (8) by A^=(XTX+μI)-1XTy, where I is the identity matrix.Similarly, we evaluate the representation contributions and representation residuals of each class of the training samples using(9)devi(i)=∥y-∑k=(i-1)×m+1i×mXkA~k∥2,
where ∑k=(i-1)×m+1i×mXkA~k and devi(i) are viewed as the representation contribution and the representation residuals of the training samples from the ith class, respectively, and then classify the testing sample, y, into the class which leads to the minimum representation residual.


## 2.1. Linear Regression Classification Algorithm

Suppose that there arec known pattern classes with N training samples x1,x2,…,xN and the query sample is y. All the samples should be firstly normalized to unit column vectors with the length of 1.LetXi be the matrix composed by the training samples of the ith class; that is, Xi=[xi1,xi2,…,xiMi]∈Rd×Mi where d is the dimension of the sample and Mi is the number of training samples of the ith class. Let us note a matrix X=[X1,X2,…,Xc]∈Rd×N and, apparently, the matrix X denotes the whole of training samples.Linear regression classification (LRC) method [35] first exploits a linear combination of each class of training samples to represent the testing sample y; that is, the following equation is satisfied:
(1)y=a1xi1+a2xi2+⋯+anxiMi,
where xi1,xi2,…,xiMi are the ith class of the training samples and aj (j=1,2,…,n) is the coefficient of xij. For convenient presentation, we can rewrite (1) into the following equation:
(2)y=Xiβi,
where Xi=[xi1,xi2,…,xiMi]∈Rd×Mi, βi=[a1,a2,…,an]T, and y∈Rd×1 are column vectors.According to the least square theory,β can be calculated by the equation below:
(3)βi=(XiTXi)-1XiTy,i=1,2,…,c.We refer toy~=Xiβi as the representation result of the ith class of the training samples using preobtained βi and we can transform y~ into a two-dimension image with the same size as the original sample image.Finally, we calculate the representation residuals of each class of the training samples and then exploit the representation residuals to conduct classification. The residual of representing the testing sample,y, using the training samples of the ith class is
(4)resi(i)=∥y-Xiβi∥2.We noteXiβi as the representation contribution using the training samples from the ith class and the RBCM presumes that the test sample, y, can be adequately represented by only the training samples from the same class [41]. Thus, a smaller representation residual can reflect a larger representation contribution, and then the minimum representation residual, resi(i), means the nearest distance between the testing sample and the ith class. Consequently, we classify the testing sample, y, into the class which achieves the smallest representation residual.


## 2.2. Collaborative Representation Classification Algorithm

Collaborative representation classification (CRC) [30] method first utilizes all the training samples to represent the test sample and exploits the linear combination of the whole of training samples to approximate the test sample. That is, the following approximation is satisfied:
(5)y=a1x1+a2x2+⋯+anxN
or
(6)y=XA,
where y is the testing sample. The ai (i=1,2,…,N) is the coefficient of xi and A=[a1a2⋯an]T. The solution of A can be calculated by solving l2-norm minimization problem:
(7)A^=argmin⁡⁡∥A∥2s.t.y=XAor(8)min⁡(∥y-XA∥2+μ∥A∥2),
where μ is a small positive constant and (8) is to utilize the Lagrangian Method to obtain the solution of (7). We solve (8) by A^=(XTX+μI)-1XTy, where I is the identity matrix.Similarly, we evaluate the representation contributions and representation residuals of each class of the training samples using(9)devi(i)=∥y-∑k=(i-1)×m+1i×mXkA~k∥2,
where ∑k=(i-1)×m+1i×mXkA~k and devi(i) are viewed as the representation contribution and the representation residuals of the training samples from the ith class, respectively, and then classify the testing sample, y, into the class which leads to the minimum representation residual.


## 3. The Proposed Method

This section describes the proposed method of integrating CRC and LRC in detail. As presented above, (1) just uses training samples of each class to represent the test sample, whereas (5) uses all the training samples to represent the test sample.The literature [30, 34] has proved that CRC method possesses limited-sparse representation capability to perform an excellent RBCM for face recognition. The literature [35, 36] also has shown the LRC method owns the ability to address the challenges of varying facial expressions, illumination, and occlusion. Thus, the motivation of the proposed method is to combine LRC and CRC to make robust and effective face recognition.Therefore, we propose a novel model by integrating CRC method and LRC method. Then the objective function can be constructed by(10)argmin⁡⁡(∥y-XA∥2+μ∥A∥2+γ∥y-XiAi∥),
where y is the test sample and X=[X1⋯Xc]. Xi is the matrix consisting of the training samples from the ith class. A is the coefficient of X and Ai is the coefficient of Xi. γ and μ are small positive constants.Equation (10) is a joint optimization problem of CRC in (8) and LRC in (2). The main steps of the proposed method to solve (10) are described as follows.Step 1.
Normalizing all the samples including the training samples and test samples to unit vectors of length 1 usingxi=xi/∥xi∥ and yi=yi/∥yi∥.Step 2.
Using CRC method by solving (8) to obtain each class of representation residual, that is, applying (9) to compute the representation residual of the jth class, devi(j).Step 3.
Using LRC method by solving (2) to obtain each class of the representation residual, that is, utilizing (4) to calculate the representation residual of the ith class, resi(i).Step 4.
Normalizing representation residuals,devi and resi. For test sample y, devi and resi are normalized by
(11)devi′=revi-revimin⁡devimax⁡-devimin⁡,resi′=resi-resimin⁡resimax⁡-resimin⁡,
where devimax⁡ and resimax⁡ denote the maximum representation residual of devi and resi, respectively. devimin⁡ and resimin⁡ are the minimum representation residual of devi and resi.Step 5.
Integrating two different residuals obtained by CRC and LRC. For the test sampley, we calculate the fusing score of the test sample with respect to the training samples of ith class using
(12)disti=devii′+γresii′,
where γ is the weight of the representation residual.Step 6.
Outputting the identity of the test sampley as
(13)identity(y)=argmin⁡i⁡(dist⁡).
We believe that a lower representation residual can lead to a more accurate representation and we assign the test sample to the class with the smallest final representation residual.


## 4. Potential Rationale of the Proposed Method

In this section, we analyze our method for exploring its underlying characteristics. We will discuss the effectiveness of our method from the following three aspects: the significance of the proposed method, the advantages of the local representation based classification, and the superiorities of the global representation based classification.CRC method and LRC method use different ideas of ideologies to make representation based classification and the correlation coefficients can reflect the difference between the two methods. Figure1 illustrates the correlation coefficients of these two methods. We treated the first three samples of each subject in ORL database [42] as the training samples and regarded the rest of the samples as test samples and then calculated the representation residual of each class using LRC and CRC, respectively. And the representation residuals from all the classes for each testing sample can form a residual vector, and then we compute the correlation coefficients of residual vectors obtained by CRC and LRC. From Figure 1, we can see that the correlation coefficients are from 0.2148 to 0.7765. That means CRC and LRC have weak relationships with each other. Thus, combining the two methods is meaningful and significant. For example, we conducted an experiment on ORL database and used the first sample from each subject as training sample and the rest of the samples as test samples. The testing sample shown in Figure 2(a) from the 21th class is classified into the 36th class, Figure 2(b), using CRC, but LRC assigns the testing sample to the 19th class, Figure 2(c). However, our method can make a correct classification into the 21th class, Figure 2(d). Figure 3 presents the representation residuals using CRC, LRC, and the proposed method, respectively. Moreover, the experimental results in Section 5 also verify that it is practical to exploit the proposed method to make effective classification.Figure 1
The correlation coefficient of CRC and LRC. We exploit the first three samples of each subject as training sample on the ORL database and the remaining samples as test samples. And the figure illustrates the correlation coefficient of representation residuals of each testing sample using CRC and LRC.A test sample that is erroneously and correctly classified by CRC, LRC, and the proposed method, respectively. (a) The test sample. (b), (c) One sample from the class wrongly classified by CRC and LRC. (d) One sample from the class correctly classified by our method.
(a)
(b)
(c)
(d)Figure 3
The representation residuals of the test sample obtained by CRC, LRC, and our method. We use the first sample of each subject as training sample and the remaining samples as test samples on the ORL database. We can see that the minimum representation residuals using CRC and LRC are the 36th and 19th class of training samples, respectively, but the smallest representation residual using our method is the 21th training sample.Representation based classification method (RBCM) usually supposes that only the training samples from the same class can sufficiently represent the test sample [41]. However, if one subject extremely resembles another subject, the above assumption may not stand. For example, we performed an experiment on FERET face database [43] using the first sample of each subject as training sample and the remaining samples for testing samples. Figure 4 shows that for the test sample wrongly and accurately classified by CRC and LRC method, our method certainly can assign the test sample to the right class. Figure 5 presents the coefficients of the representation solution, and it is apparent that the largest two coefficients are the coefficients of training samples from the 110th class and the 72th class. However, CRC assigns the test sample to the 72th class that is the nearest class. On the contrary, LRC and our method can exactly recognize the test sample. Thus, the local representation based classification method can compensate for the imperfection of the global representation based recognition method. From the viewpoint of the proposed method, our method can sufficiently implement the two typical characteristics of the local representation based classification, LRC, and the global representation based classification method, CRC.A test sample that is falsely and accurately classified by CRC and LRC. (a) A test sample from the 110th class. (b) A sample from the 72th class which is the closest to the test sample. (c) One sample from the 110th class precisely classified by LRC method and the proposed method.
(a)
(b)
(c)Figure 5
The coefficients of the training samples obtained by CRC. The first sample of each subject on the FERET database is regarded as training sample and the remaining samples are treated as test samples. This figure illustrates the coefficient of training samples exploiting all the training samples to represent the testing sample, Figure4(a), from the FERET face database.It is obvious that global representation based classification method, which uses all the training samples to represent the test sample, can alleviate the difficulty of the small sample size problem. For example, if we just use one sample of each subject to make classification by exploiting LRC method, LRC will deteriorate to the nearest neighbor (NN) classifier and also lose the preponderance of RBCM. In that regard, score fusion of CRC and LRC is meaningful and can obtain better result.CRC and SRC are the typical examples of global RBCMs and LRC is one of the typical local RBCMs. However, SRC has very high computation complexity and CRC owns the superiority that CRC is extremely faster than SRC with very competitive face recognition accuracy. What is more, CRC and LRC will be almost not able to identify the class to which the test sample truly belongs and will erroneously classify the test sample. However, the proposed method can compensate for the imperfection and shortcomings of each other. Moreover, the correlation coefficients between these two methods are very low. Thus, integrating CRC and LRC can take full advantages of these two methods. At this stage, we can demonstrate that the framework of integrating the advantages of global RBCM and the superiorities of local RBCM is reasonable and meaningful.


## 5. Experimental Results

In this section, we conducted extensive experiments to evaluate the effectiveness and robustness of the proposed method on different face databases including ORL, AR, CMU PIE, and FERET database. In order to present the performance of the proposed method, several competitive face recognition methods are tested as comparison, such as CRC [30], LRC [35], LDA, improved nearest neighbor classification (INNC) method [32], and the sparse representation based classification (SRC) algorithm proposed in [29]. Moreover, we also performed experiments on corrupted face databases and conducted contiguous occlusion experiments. The parameter of the proposed method is set to μ=0.01 for all the experiments.


### 5.1. Experiment on the ORL Face Database

The ORL face database [42] includes 400 face images obtained from 40 distinct subjects and each subject provides ten different face images. For some subjects, the images were taken at different times, varying the lighting conditions, facial expressions (open or closed eyes, smiling, or not smiling), and facial details (glasses or no glasses). All the images were captured against a dark homogeneous background with the subjects in an upright, frontal position. In order to reduce the computational complexity, all the images were resized to 56 by 46 image matrix at first. Figure 6 presents some images from ORL database. When the number of training sample is one, LDA algorithm is not available since the between-class scatter matrix is zero.Figure 6
Some face images from the ORL face database.We took the first 1, 2, 3, 4, and 5 image samples of each subject as original training samples, respectively, and regarded the remaining face images as test samples. Table1 enumerates the variation of the classification error rates of different solutions. And our method with different values of γ also verifies that the proposed method outperforms the other classification methods. For example, when we took the first four face images of each subject as training samples and the rest of the face images as test samples, the classification error rate of our method can greatly outperforms the CRC, LRC, INNC, LDA, and SRC algorithms by a margin of 2.5 percent, 5.42 percent, 9.17 percent, 6.25 percent, and 4.17 percent, respectively. What is more, it is clear that the proposed method can dramatically improve CRC and LRC classification methods and it also reflects the effectiveness of the framework of integrating the global RBCM and local RBCM.Table 1
Rates of classification errors of different methods on the ORL face database.
1
2
3
4
5
Our method (γ=0.9)
0.2944
0.1500
0.1321
0.0833
0.1100
Our method (γ=0.85)
0.2972
0.1469
0.1357
0.0833
0.1100
Our method (γ=0.8)
0.2972
0.1500
0.1321
0.0833
0.1050
CRC
0.3194
0.1656
0.1392
0.1083
0.1150
LRC
0.3250
0.2063
0.1857
0.1375
0.1200
INNC
0.3194
0.2125
0.2178
0.1750
0.1700
LDA
—
0.2531
0.1678
0.1458
0.1750
SRC
0.3177
0.1781
0.1607
0.1250
0.1200


### 5.2. Experiments on the AR Face Database

The AR face database [44] is composed of more than 3000 face images of 126 subjects and 26 facial images for each subject were taken in two separated terms and the resolution of each AR image has been resized to a 40 by 50 image matrix. This face database has much more variations, such as, illumination conditions, facial expressions, and facial disguises. We treated the first 1,2,…,9 face images from each subject as the original training samples and regarded the remaining samples as test samples, respectively. Figure 7 has shown some face images for four subjects from the AR database. Table 2 shows the classification error rates using different classification methods. It is obvious that our method has remarkable priorities compared to the other classification algorithms. For example, when we took the first face image of each subject as training samples and the rest of the face images as test samples, the classification error rates of our method with γ=0.90, γ=0.75, γ=0.65, CRC, LRC, INNC, and SRC are 29.06%, 28.90%, 28.87%, 30.17%, 38.07%, 30.17%, and 31.20%, respectively. That means that 9.2 percent recognition errors can be avoided by using our method instead of LRC. We can see that combining CRC and LRC is a favorable idea to perform face recognition.Table 2
Rates of classification errors of different methods on the AR face database.
1
2
3
4
Our method (γ=0.90)
0.2906
0.3045
0.3061
0.3166
Our method (γ=0.75)
0.2890
0.3059
0.3029
0.3125
Our method (γ=0.65)
0.2887
0.3069
0.2996
0.3140
CRC
0.3017
0.3167
0.3141
0.3246
LRC
0.3807
0.4038
0.4094
0.4186
INNC
0.3017
0.3299
0.3301
0.3413
LDA
—
0.3965
0.3880
0.4102
SRC
0.3120
0.3156
0.3167
0.3220Figure 7
Some face images from the AR face database.


### 5.3. Experiments on FERET Face Database

The FERET database is one of the largest publicly available facial images databases. We just selected a subset made up of 1400 images from 200 individuals and each one contains seven images to conduct experiments [43]. The images were collected in a semicontrolled environment. We resized each image to a 40 by 40 image matrix using the same downsampling algorithm. For some individuals, there were over two years between their first and last time to be photographed. Figure 8 presents some cropped face images from the FERET face database. The parameter of our method is set to γ=1.0 for this experiment.Figure 8
Some face images from the FERET face database.We perform an experiment on FERET face database using the first six samples of each subject as training sample and the remaining samples for testing. The classification error rates exploiting our method, CRC, LRC, INNC, LDA, and SRC are 15.00%, 25.50%, 21.50%, 25.50%, 41.50%, and 18.50%, respectively. In other words, this suggests that the proposed method can greatly outperform the CRC, LRC INNC, LDA, and SRC algorithms by a margin of 10.5 percent, 6.5 percent, 10.5 percent, 26.5 percent, and 2.5 percent, respectively.The first image of each subject was used for training sample, while the remaining six were served as test samples. Also Figure9 shows the detailed experimental results compared to other classification methods for a large variation of feature dimensions. From Figure 9, we can see that our method achieves the lowest classification error rate compared to the other classification methods although the dimensions of the samples are extremely low. The classification error rates under different feature dimensions suggest that the proposed method is superior to other classifiers.Figure 9
The recognition error rates of different method under different feature dimensions. We treat the first image of each subject as training sample from FERET face database and the remaining images are used for test samples. The figure shows the detailed experimental results compared to other classification methods for a large variation of feature dimensions.


### 5.4. Experiments on CMU PIE Database

We chose a subset of the CMU PIE database to perform experiments. For each subject in this face database, we randomly selected 21 face images under different lighting conditions and expressions including 1428 images of 68 subjects for experiments [45]. The size of each image is cropped to 32 by 32 pixels. We just took the first 1, 2, and 3 face images of each subject as the training samples and the rest of the samples as test sample. The experimental results are shown in Table 3 and we can draw a conclusion that the proposed method can achieve the minimum classification error rate. For example, when the number of training sample is one, the classification error rates using our method with γ=1.0, γ=0.85, γ=0.70, CRC, LRC, INNC, and SRC are 10.80%, 11.17%, 11.32%, 16.10%, 13.82%, 16.10%, and 16.10%, respectively. Thus the proposed method can extensively improve the classification accuracy compared to CRC and LRC methods and our method obtains the lowest classification error rates.Table 3
Rates of classification errors of different methods on the illumination variation subset of CMU PIE face database.
1
2
3
The proposed method (γ=1.0)
0.1080
0.0542
0.0049
The proposed method (γ=0.85)
0.1117
0.0503
0.0049
The proposed method (γ=0.70)
0.1132
0.0479
0.0032
CRC
0.1610
0.0975
0.0080
LRC
0.1382
0.1060
0.0849
INNC
0.1610
0.0588
0.0049
LDA
—
0.1262
0.0131
SRC
0.1610
0.1370
0.0179


### 5.5. Recognition on Random Pixel Corruption Databases

We performed experiments on random pixel corrupted face images to test the robustness of the methods. The datasets of training samples and the test samples are the same as those in Sections5.1 and 5.2, respectively. In order to obtain random pixel corruption face images, we add Gaussian white noise of zero mean and variance of 0.01 to the ORL database and SaltPepper noise of density of 0.1 to AR database, respectively. Figure 10 shows some corrupted face images on ORL database and AR database. Table 4 and Figure 11 show the classification error rate results on ORL database and AR database, respectively. The LDA algorithm is not available when the number of the training sample for each class is one, so we set the classification error rate exploiting the LDA algorithm to 0 for corrupted AR database in Figure 11. We see again that our method possesses the absolute dominance compared to all the other methods. For example, when each subject provides 3 corrupted ORL face images as training samples, the classification error rates using our method with γ=0.65, CRC, LRC, INNC, LDA, and SRC are 15.00%, 17.14%, 22.86%, 21.78%, 29.64%, and 16.79%, respectively. Thus the experimental results also demonstrate the robustness and effectiveness of the proposed method.Table 4
Rates of classification errors of different methods on corrupted ORL face database.
2
3
4
5
Our method (γ=0.65)
0.1656
0.1500
0.0916
0.0900
Our method (γ=0.60)
0.1656
0.1500
0.0875
0.0950
CRC
0.1844
0.1714
0.1083
0.1100
LRC
0.2594
0.2286
0.1625
0.1450
INNC
0.2312
0.2178
0.1583
0.1600
LDA
0.4625
0.2964
0.2125
0.1600
SRC
0.1813
0.1679
0.1083
0.1000Figure 10
Some corrupted face images in ORL database and AR database. The first two rows show the corrupted face images used in the ORL database. The last two rows show the corrupted face images used in the AR face database.Figure 11
Rates of classification error rates of different methods on the corrupted AR face database. We exploit the first1,2,…,9 face images of each subject as training samples and take the remaining samples as test samples. The figure shows the recognition error rates using our method with the parameter γ=0.95 and γ=0.85, CRC, LRC, INNC, LDA, and SRC, respectively.


## 5.1. Experiment on the ORL Face Database

The ORL face database [42] includes 400 face images obtained from 40 distinct subjects and each subject provides ten different face images. For some subjects, the images were taken at different times, varying the lighting conditions, facial expressions (open or closed eyes, smiling, or not smiling), and facial details (glasses or no glasses). All the images were captured against a dark homogeneous background with the subjects in an upright, frontal position. In order to reduce the computational complexity, all the images were resized to 56 by 46 image matrix at first. Figure 6 presents some images from ORL database. When the number of training sample is one, LDA algorithm is not available since the between-class scatter matrix is zero.Figure 6
Some face images from the ORL face database.We took the first 1, 2, 3, 4, and 5 image samples of each subject as original training samples, respectively, and regarded the remaining face images as test samples. Table1 enumerates the variation of the classification error rates of different solutions. And our method with different values of γ also verifies that the proposed method outperforms the other classification methods. For example, when we took the first four face images of each subject as training samples and the rest of the face images as test samples, the classification error rate of our method can greatly outperforms the CRC, LRC, INNC, LDA, and SRC algorithms by a margin of 2.5 percent, 5.42 percent, 9.17 percent, 6.25 percent, and 4.17 percent, respectively. What is more, it is clear that the proposed method can dramatically improve CRC and LRC classification methods and it also reflects the effectiveness of the framework of integrating the global RBCM and local RBCM.Table 1
Rates of classification errors of different methods on the ORL face database.
1
2
3
4
5
Our method (γ=0.9)
0.2944
0.1500
0.1321
0.0833
0.1100
Our method (γ=0.85)
0.2972
0.1469
0.1357
0.0833
0.1100
Our method (γ=0.8)
0.2972
0.1500
0.1321
0.0833
0.1050
CRC
0.3194
0.1656
0.1392
0.1083
0.1150
LRC
0.3250
0.2063
0.1857
0.1375
0.1200
INNC
0.3194
0.2125
0.2178
0.1750
0.1700
LDA
—
0.2531
0.1678
0.1458
0.1750
SRC
0.3177
0.1781
0.1607
0.1250
0.1200


## 5.2. Experiments on the AR Face Database

The AR face database [44] is composed of more than 3000 face images of 126 subjects and 26 facial images for each subject were taken in two separated terms and the resolution of each AR image has been resized to a 40 by 50 image matrix. This face database has much more variations, such as, illumination conditions, facial expressions, and facial disguises. We treated the first 1,2,…,9 face images from each subject as the original training samples and regarded the remaining samples as test samples, respectively. Figure 7 has shown some face images for four subjects from the AR database. Table 2 shows the classification error rates using different classification methods. It is obvious that our method has remarkable priorities compared to the other classification algorithms. For example, when we took the first face image of each subject as training samples and the rest of the face images as test samples, the classification error rates of our method with γ=0.90, γ=0.75, γ=0.65, CRC, LRC, INNC, and SRC are 29.06%, 28.90%, 28.87%, 30.17%, 38.07%, 30.17%, and 31.20%, respectively. That means that 9.2 percent recognition errors can be avoided by using our method instead of LRC. We can see that combining CRC and LRC is a favorable idea to perform face recognition.Table 2
Rates of classification errors of different methods on the AR face database.
1
2
3
4
Our method (γ=0.90)
0.2906
0.3045
0.3061
0.3166
Our method (γ=0.75)
0.2890
0.3059
0.3029
0.3125
Our method (γ=0.65)
0.2887
0.3069
0.2996
0.3140
CRC
0.3017
0.3167
0.3141
0.3246
LRC
0.3807
0.4038
0.4094
0.4186
INNC
0.3017
0.3299
0.3301
0.3413
LDA
—
0.3965
0.3880
0.4102
SRC
0.3120
0.3156
0.3167
0.3220Figure 7
Some face images from the AR face database.


## 5.3. Experiments on FERET Face Database

The FERET database is one of the largest publicly available facial images databases. We just selected a subset made up of 1400 images from 200 individuals and each one contains seven images to conduct experiments [43]. The images were collected in a semicontrolled environment. We resized each image to a 40 by 40 image matrix using the same downsampling algorithm. For some individuals, there were over two years between their first and last time to be photographed. Figure 8 presents some cropped face images from the FERET face database. The parameter of our method is set to γ=1.0 for this experiment.Figure 8
Some face images from the FERET face database.We perform an experiment on FERET face database using the first six samples of each subject as training sample and the remaining samples for testing. The classification error rates exploiting our method, CRC, LRC, INNC, LDA, and SRC are 15.00%, 25.50%, 21.50%, 25.50%, 41.50%, and 18.50%, respectively. In other words, this suggests that the proposed method can greatly outperform the CRC, LRC INNC, LDA, and SRC algorithms by a margin of 10.5 percent, 6.5 percent, 10.5 percent, 26.5 percent, and 2.5 percent, respectively.The first image of each subject was used for training sample, while the remaining six were served as test samples. Also Figure9 shows the detailed experimental results compared to other classification methods for a large variation of feature dimensions. From Figure 9, we can see that our method achieves the lowest classification error rate compared to the other classification methods although the dimensions of the samples are extremely low. The classification error rates under different feature dimensions suggest that the proposed method is superior to other classifiers.Figure 9
The recognition error rates of different method under different feature dimensions. We treat the first image of each subject as training sample from FERET face database and the remaining images are used for test samples. The figure shows the detailed experimental results compared to other classification methods for a large variation of feature dimensions.


## 5.4. Experiments on CMU PIE Database

We chose a subset of the CMU PIE database to perform experiments. For each subject in this face database, we randomly selected 21 face images under different lighting conditions and expressions including 1428 images of 68 subjects for experiments [45]. The size of each image is cropped to 32 by 32 pixels. We just took the first 1, 2, and 3 face images of each subject as the training samples and the rest of the samples as test sample. The experimental results are shown in Table 3 and we can draw a conclusion that the proposed method can achieve the minimum classification error rate. For example, when the number of training sample is one, the classification error rates using our method with γ=1.0, γ=0.85, γ=0.70, CRC, LRC, INNC, and SRC are 10.80%, 11.17%, 11.32%, 16.10%, 13.82%, 16.10%, and 16.10%, respectively. Thus the proposed method can extensively improve the classification accuracy compared to CRC and LRC methods and our method obtains the lowest classification error rates.Table 3
Rates of classification errors of different methods on the illumination variation subset of CMU PIE face database.
1
2
3
The proposed method (γ=1.0)
0.1080
0.0542
0.0049
The proposed method (γ=0.85)
0.1117
0.0503
0.0049
The proposed method (γ=0.70)
0.1132
0.0479
0.0032
CRC
0.1610
0.0975
0.0080
LRC
0.1382
0.1060
0.0849
INNC
0.1610
0.0588
0.0049
LDA
—
0.1262
0.0131
SRC
0.1610
0.1370
0.0179


## 5.5. Recognition on Random Pixel Corruption Databases

We performed experiments on random pixel corrupted face images to test the robustness of the methods. The datasets of training samples and the test samples are the same as those in Sections5.1 and 5.2, respectively. In order to obtain random pixel corruption face images, we add Gaussian white noise of zero mean and variance of 0.01 to the ORL database and SaltPepper noise of density of 0.1 to AR database, respectively. Figure 10 shows some corrupted face images on ORL database and AR database. Table 4 and Figure 11 show the classification error rate results on ORL database and AR database, respectively. The LDA algorithm is not available when the number of the training sample for each class is one, so we set the classification error rate exploiting the LDA algorithm to 0 for corrupted AR database in Figure 11. We see again that our method possesses the absolute dominance compared to all the other methods. For example, when each subject provides 3 corrupted ORL face images as training samples, the classification error rates using our method with γ=0.65, CRC, LRC, INNC, LDA, and SRC are 15.00%, 17.14%, 22.86%, 21.78%, 29.64%, and 16.79%, respectively. Thus the experimental results also demonstrate the robustness and effectiveness of the proposed method.Table 4
Rates of classification errors of different methods on corrupted ORL face database.
2
3
4
5
Our method (γ=0.65)
0.1656
0.1500
0.0916
0.0900
Our method (γ=0.60)
0.1656
0.1500
0.0875
0.0950
CRC
0.1844
0.1714
0.1083
0.1100
LRC
0.2594
0.2286
0.1625
0.1450
INNC
0.2312
0.2178
0.1583
0.1600
LDA
0.4625
0.2964
0.2125
0.1600
SRC
0.1813
0.1679
0.1083
0.1000Figure 10
Some corrupted face images in ORL database and AR database. The first two rows show the corrupted face images used in the ORL database. The last two rows show the corrupted face images used in the AR face database.Figure 11
Rates of classification error rates of different methods on the corrupted AR face database. We exploit the first1,2,…,9 face images of each subject as training samples and take the remaining samples as test samples. The figure shows the recognition error rates using our method with the parameter γ=0.95 and γ=0.85, CRC, LRC, INNC, LDA, and SRC, respectively.


## 6. Conclusion

The proposed method combines CRC and LRC methods to perform effective and robust face recognition. The significance of our method is to sufficiently exploit the globality of CRC and the locality of LRC, respectively. The potential rationale of this method is that integrating LRC and CRC is meaningful because they have the low correlation coefficients showed in Figure1. Moreover, it is remarkable that our method obtains a promising performance on face databases with random pixel corruption, and then we can confirm the good efficiency of the proposed method under complex conditions.All the experiments in Section5 can testify that the framework of integrating the globality and locality of training samples can tremendously improve the classification accuracy. CRC, as one of the typical global RBCMs, can completely utilize the globality of all the training samples, whereas LRC, as one of the typical local RBCMs, can make the most of the locality of the training samples. A large number of experiments also demonstrate that the proposed framework is reasonable and practical and it can extensively improve the classification accuracy compared to the state-of-art classification methods on different face databases. Furthermore, the proposed method is simple and computationally efficient. Besides the method description and the experimental analysis, the rationale of our method also has been explicitly interpreted and the superiority of our method also has been visually presented in this paper.



---

*Source: 415856-2014-03-26.xml*