Research Article  Open Access
Facial Feature Extraction Using Frequency Map Series in PCNN
Abstract
Pulse coupled neural network (PCNN) has been widely used in image processing. The 3D binary map series (BMS) generated by PCNN effectively describes image feature information such as edges and regional distribution, so BMS can be treated as the basis of extracting 1D oscillation time series (OTS) for an image. However, the traditional methods using BMS did not consider the correlation of the binary sequence in BMS and the space structure for every map. By further processing for BMS, a novel facial feature extraction method is proposed. Firstly, consider the correlation among maps in BMS; a method is put forward to transform BMS into frequency map series (FMS), and the method lessens the influence of noncontinuous feature regions in binary images on OTSBMS. Then, by computing the 2D entropy for every map in FMS, the 3D FMS is transformed into 1D OTS (OTSFMS), which has good geometry invariance for the facial image, and contains the space structure information of the image. Finally, by analyzing the OTSFMS, the standard Euclidean distance is used to measure the distances for OTSFMS. Experimental results verify the effectiveness of OTSFMS in facial recognition, and it shows better recognition performance than other feature extraction methods.
1. Introduction
Face recognition is an important research field of pattern recognition; it has good potential applications in biological recognition technology, security system, and so on. At present, there are many face images through sensors, so we need a good algorithm to deal with these images. In the process, because of face data space caused by the problem of dimension disaster, facial feature extraction method with spatial dimension reduction effect is becoming the key technology of face recognition. In the past several decades, many researchers proposed a lot of methods to extract facial features such as geometric characteristics, subspace analysis [1–5], and neural network method [6].
The method of geometric characteristics [7] uses the calculation of geometric parameters as the face features; it has good adaptability to illumination changes, but poor adaptability to the more obvious facial expression, posture, and rotation changes, and so forth. At present, the subspace analytical method is a popular face recognition method, it employs a transform method of linear or nonlinear to make the data in the cast shadow space embodying explicit feature pattern, so as to extract the key features such as the method of PCA [1], LDA [2], and ICA [3] based on linear transformation and the method of KPCA [4], KFDA [5], and KICA based on nonlinear transformation. However, these methods have poor adaptability to the changes of rotation and distortion in the image. Neural network is based on the nonlinear transform ability of the network structure and uses the learning of the training sample to get nonlinear transforming space of the data and then to obtain the facial features according to the nonlinear transforming space, for instance, using SOM neural network and fuzzy RBF neural network, and so forth. But neural network method will likely cause overfitting in the learning process for the samples of empirical risk minimization principle. Besides, the subspace analysis method and the face feature extraction method of traditional neural network need face samples to learn; if the training samples are changed, the projection transformation space of the data also wants to change; thus, the calculated amount in largescale human face feature extraction is too large; its application in the realtime demand higher occasion is limited.
In 1993, Johnson and Ritter proposed pulse coupled neural network (PCNN) [8] based on Eckhorn research in cat’s visual cortex. It has widely been used in image segmentation [9, 10], image fusion [11–13], image retrieval [14], and so forth. In PCNN, the similarity group neurons will issue synchronous pulses under the effect of mutual coupling pulses. These pulses information constitute a 3D binary map series (BMS), which effectively describes the information of the edge and regional distribution of the image. But the data size of the BMS is larger and cannot be directly used as the image features; for this reason, by calculating the area of the binary image, Johnson [15] transform BMS into 1D oscillation time sequences (OTS), and the OTS has a good invariance in geometric transforming such as rotation, translation, and zoom. Based on the pulse mechanism of the neuron in PCNN, the pulses are divided into the capture pulses and selfexcitation pulses, and the OTS is divided into COTS and SOTS, and serve as the facial features extraction in face recognition [16]. Reference [14] extracted 1D OTS of the BMS by calculating the normalization rotational inertia (NRI) of the binary image and applying it to the image retrieval. However, the OTS feature extraction method of the image based on BMS did not fully consider the correlation between BMS in the binary images; those discontinuous features regions will cause influence for the pattern classification capability of OTS features. In addition, the OTS of Johnson’s form is statistical characteristics in the sense of whole situation of the binary images; these did not consider the spatial structure of the image, and the spatial structure information often plays an important role in pattern classification.
In view of the above analysis, this paper proposed a novel face feature extraction method of OTS based on the BMS of PCNN output, and, compared with the traditional subspace analysis and neural network method, the results will not change with the sample space change.
2. Pulse Coupled Neural Network
The PCNN model consists of the receptive field, the modulation field, and the pulse generator. In the receptive field, the neuron, respectively, receives the coupling pulse input and external stimulus input of neighboring neurons and consists of and channels, which is described by (1). In and channels of the neuron, the neuron links with its neighborhood neurons via the synaptic linking weights and , respectively; the two channels accumulate input and exponential decay changes at the same time; the decay exponentials of and channels are and , while the channel amplitudes are and :
In the modulation field, the linking input is added a constant positive bias; then, it is multiplied by the feeding input; the bias is unitary, is the linking strength, and the total internal activity is the result of modulation, which is described by
Pulse generator consists of a threshold adjuster, a comparison organ, and a pulse generator. Its function is to generate the pulse output , and it adjusts threshold value ; is threshold range coefficient, which is described by (3). When the internal state is larger than the threshold , the neuron generates a pulse, which is described by
In the above equations, the subscripts and denote the neuron location in a PCNN and denotes the current iteration (discrete time step), where varies from 1 to ( is the total number of iterations).
The PCNN used for image processing is a single layer twodimensional array of laterally linked pulse coupled neurons as shown in Figure 1, and all neurons are identical. The number of neurons in the network is equal to the number of pixels in the input image. There exists onetoone correspondence between the image pixels and network neurons, and the gray value of a pixel is taken as the external input stimulus of the neuron in channel; namely, . The output of each neuron results is two states, namely, pulse (1 state) and nonpulse (0 state), so the output states of neurons comprise a binary map.
3. Face Feature Extraction Using PCNN
In PCNN, each neuron is connected with neighboring neurons in linking range, so a neuron will receive some pulse inputs from its adjacent neurons. If adjacent neurons have a similar external stimulus with the current neuron, the neuron will issue pulse because of pulse coupled action; that is to say, the neuron and its similar adjacent neurons will emit synchronous pulses. So the similar group of the PCNN neurons possesses the pulse synchronization characteristic and the characteristic benefit features clustering for an image.
3.1. Feature Extraction Based on Binary Map Series of PCNN
The PCNN will output a binary map in each iteration, so a binary map series (BMS) which contains binary maps will be generated after iterations, which is recorded and is shown in Figure 2.
Figure 2 shows a face image in the upper left corner, and others are the face image of BMS generated by PCNN. The results in Figure 2 show that the BMS effectively reflects the edge details and regional features distribution, and on the time axis of the image sequence, it well demonstrates that the process produces changes of the image features by neighborhood neurons pulse coupling. In addition, some characteristics of the regional component cycle are repeated with certain cycle in BMS, but some feature components are contracting on the time axis, while others are expanding. The phenomenon actually shows that the PCNN operates the features region clustering in smooth form.
In PCNN, an image of size will generate a BMS with the size ; this means the data of the BMS is times greater than that of the original image, so the BMS cannot be directly applied to pattern classification. If each 2D binary map in BMS could be translated into 0D data points by some means or another, then the 3D BMS could be translated into a 1D oscillation time sequence (OTS) of the size ; thus, the OTS realizes the feature extraction and data dimension reduction of the image. Because this kind of OTS is generated by BMS, it can be denoted as OTSBMS. Johnson [15] put forward a method of translating a binary image into 0D data points by calculating the area of the binary map; thus, this kind of OTSBMS can be described by (5); an example is shown in Figure 3:
3.2. Frequency Map Series and Feature Extraction
OTSBMS defined by (5) is a global statistical result of each binary image, so in a binary map, those pulses (1 state) at different location and not in same feature region will play the same role in the feature description with OTSBMS. Thus, this king of OTSBMS not only failed to properly describe the correlation among binary maps in chronological order but also cannot well describe the spatial structure information of a binary map, but the information is very important for improving the effectiveness of pattern classification using BMS. In view of this, we define a frequency map series (FMS) based on the above BMS, which is described by (6), and an example is shown in Figure 4:
In FMS, if each 2D frequency map can be transformed into 0D feature points, then 1D time sequence signature OTSFMS can be extracted from 3D FMS. In addition, to make the 0D feature points as possible as including the space structure information of a frequency map, here we use twodimensional entropy of an image to extract 0D feature information for a frequency map. The twodimensional entropy of an image is defined by (8):
In (8), denotes gray levels of the image. In (9), represents pixel gray value, and represents neighbor mean value of this pixel. Obviously, twodimensional entropy not only transforms the image of the size into 0D data points, because of reflecting the relation between pixel gray value and neighbor mean value, but also contains the spatial structure information of the image. Through calculating the twodimensional entropy of the image, the FMS of the size can be transformed into 1D OTSFMS of the size ; that is, ; as a result, this method greatly reduces the feature dimension. In addition, OTSFMS also has good ability of pattern classification for the face image; as shown in Figure 5, OTSFMS is supremely similar with the same kind of face S1_1 and S1_7 and is different for S3_9.
(a) S1_1
(b) S1_7
(c) S3_9
(d) OTSFMS
Because the calculation of the twodimensional entropy of an image is independent of the rotation and translation of the image and has good adaptability to image scaling, therefore, this kind of OTSFMS has good invariance for the image rotation and translation and has invariance for the image scaling with small error. As shown in Figure 6, the face image is rotated 30 degrees, zoomed 0.6, and moved 15 pixels vertically and horizontally, respectively. Obviously, the image occurs geometric changes, but its OTSFMS characteristic curves have little difference.
4. Distance Measure Selection of OTSFMS
After feature extraction and dimensionality reduction of the face image, we can use common Euclidean distance (ED) as distance measure in pattern classification for the face recognition. But ED requires that the distribution of each dimension component of data should be consistent. For the OTSBMS feature data of an image with the size of , because it is global pulse statistics in data space, the mapping space of each dimension component is , so ED can be used to measure the distance among OTSBMS.
However, for OTSFMS feature data of the image, because it is monotone increasing totally, as shown in Figure 5, namely, the frequency rank of each dimension component of data also is increasing in the overall. Therefore, even if we use the same entropy method to calculate the frequency map, mapping spaces of the frequency maps are inconsistent; that is to say, their operation results of twodimensional entropy are not consistent in their dimension. So the distance measure of OTSFMS feature data is inappropriate to select ED, and it should select a distance measure which is independent of dimension. Here, standardized Euclidean distance (SED) is a suitable choice. For eigenvectors and , denotes the standard variance, and then SED is defined by
5. Experimental Results and Analysis
To evaluate the performance of the proposed facial feature extraction in face recognition, experimental studies are carried out on the ORL face database or MITCBCL face database. The ORL database contains 400 images of 40 individuals (4 females and 36 males). Each individual has 10 images varying in position, rotation, zoom, and expression. The MITCBCL database comprises 2000 images from 200 people; each person has 10 images with different illumination, pose, and rotation. The parameters of PCNN are set as , , , , , , , , , the default iterations , and the default statistics starting point .
First of all, to test the classification performance of the OTSBMS or OTSFMS when using ED and SED, we constructed experiments based on ORL database by randomly leaving out image per person each time for testing, the remainder being images per person for training. This was repeated 30 times by leaving out images per person each time. The experimental results listed in Table 1 are the average of 30 times’ results each time.

It can be seen from Table 1, under different training sets, that the recognition rate of SED measure is significantly higher than ED measure for OTSFMS, and the recognition rate of ED measure is also slightly higher than SED measure for OTSBMS. This shows that the OTSFMS by calculating twodimensional entropy of each frequency map in FMS is suitable for using SED measure, while the OTSBMS by calculating the area of each binary map in BMS is suitable for selecting ED measure. In addition, the recognition rate of OTSFMS under SED measure is significantly higher than OTSBMS under ED measure. That shows the proposed OTSFMS in this paper can effectively extract facial image features.
Then, to investigate the influences on the facial recognition performance for different iterations in PCNN or the different statistical starting point in FMS, experiments were carried out on ORL face database, by randomly choosing 6 images per person each time for testing, the remainder being 4 images per person for training. This was repeated 30 times with different and each time. The average face recognition rates in these experiments are shown in Table 2.

In Table 2, we can see that the recognition rate firstly increases with the increase of iterations under definite , but soon afterwards it is gradually stable; hence, can be set as to decrease the amount of calculation for PCNN. Under definite , the recognition rate decreases with the increase of , and to increase the recognition rate, is set as .
On the ORL database, selecting the frontal images of each person as training set, and the rest of images as test set, we make a comparison of the recognition rate of our method (OTSFMS + SED) with other methods such as common subspace face recognition algorithms (PCA, 2DPCA, and KPCA) and PCNN [16], as shown in Table 3. Obviously, under the condition , the recognition rate of the proposed method is significantly higher than other methods, but under the condition , the recognition rate in this paper is lower than other subspace methods, but it is higher than PCNN method [16]. Combined with the experimental results of Table 1 (, the average recognition rate is 92.67%), we can know that there are a few influential samples in the training sets.

In order to verify the validity of this method for face recognition in the bigger data space, on MITCBCL, we randomly select face images of each person at rate to compose training set, the remaining images as test set; then, 30 times experiments were carried out, and the average recognition rate is shown in Table 4.

As shown in Table 4, under the same training set, the recognition rate of OTSFMS in this paper is higher than [16] based on BMS, which illustrates the validity of this paper. In addition, the experimental results in MITCBCL Library are the same as ORL Library; the recognition rate using SED measure is higher than that using ED measure. This again shows the correctness of the analysis of the characteristics of the OTSFMS in this paper.
On ORL or MITCBCL face database above, the face has been tailored manually from original image, but on a practical face recognition system, obtaining the face from an image needs to use a certain face detection algorithm such as the method in [17], and the facial recognition performance is easily affected by the detection accuracy. To test the adaptability of both OTSFMS and OTSBMS for face detection precision, we set up a custom face database, which includes 208 images of 13 individuals; each individual has 16 images varying in facial expression, illumination, pose, background, and facial details (with/without glasses), and the part of the database was shown in Figure 7.
The face has a larger proportion of background on the custom database; first, we use the face detection algorithm [17] for processing (as shown in Figure 8); then, we compare the face recognition performance before and after face detection for both OTSFMS and OTSBMS, and the results are shown in Table 5. Here, experimental repetitions are 30, . As we can see in Table 5, after face detection, the average recognition rates are improved for OTSFMS and OTSBMS; this suggests that the efficient face detection method has important significance for improving identification performance in the proposed method. In addition, both before and after the face detection, the average recognition rate for the OTSFMS is better than that for the OTSBMS; it also suggests that the improvement based on OTSBMS for the facial feature extraction method (i.e., OTSFMS) is effective.

In addition, after face detection on custom face database, we add Gaussian noise to each kind of faces detected as shown in Figure 8 with the mean value 0 and the variance ; here, the variance satisfies with uniform distribution, and then we investigate the adaptation of both OTSFMS and OTSBMS for noise; the results are shown in Table 6. Here, experimental repetitions are 30, . As we can see in Tables 5 and 6, after the faces are polluted by noise, both for OTSFMS and OTSBMS, the average recognition rate is decreased, but the influence of the noise for the OTSFMS identification is more serious than OTSBMS. This is mainly due to the OTSBMS which is a global statistics on BMS of the facial image, but the OTSFMS is more emphasized on the facial local detail information. Therefore, when a face is polluted by Gaussian noise, the noise damage to the facial local details is fatal, but to the global statistics it is lighter. The fact shows that the face recognition algorithm (OTSFMS + SED) is sensitive to noise.

6. Conclusions
In this paper, a novel method was proposed to extract the facial feature based on PCNN. Through the analysis of the limitations of BMS extracted onedimensional oscillation time sequence (OTS) and considering the correlation between binary images in BMS, we proposed a method to transform BMS into frequency map series (FMS), which reduces the influences of the discontinuous value of binary image on the effectiveness of OTS. Then, the paper considers the method that 2D frequency map is transformed into 0D data points; in order to reflect the spatial structure information in the frequency map, the paper employs the computing method of twodimensional entropy of the image; it will transform 3D FMS into 1D OTSFMS feature, and the feature has a good invariance for face image geometric changes. Finally, based on the analysis of characteristics of OTSFMS data features, the paper we proposed uses the standard Euclidean distance measure as the distance measure of OTSFMS features. The experimental results show that the recognition rate of OTSFMS is significantly higher than PCA and KPCA, and so forth, and also better than method [16] based on BMS extracted OTSBMS features. And compared to traditional subspace analysis and neural network method, face features extracted by this method do not change with the sample space change.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research was financially supported by the National Natural Science Foundation of China (nos. 61365001 and 61463052), the Natural Science Foundation of Yunnan Province (no. 2012FD003), and the Science and Technology Plan of Yunnan Province (no. 2014AB016).
References
 G.C. Luh and C.Y. Lin, “PCA based immune networks for human face recognition,” Applied Soft Computing, vol. 11, no. 2, pp. 1743–1752, 2011. View at: Publisher Site  Google Scholar
 G.F. Lu, J. Zou, and Y. Wang, “Incremental complete LDA for face recognition,” Pattern Recognition, vol. 45, no. 7, pp. 2510–2521, 2012. View at: Publisher Site  Google Scholar
 I. Marques and M. Graña, “Face recognition with lattice independent component analysis and extreme learning machines,” Soft Computing, vol. 16, no. 9, pp. 1525–1537, 2012. View at: Publisher Site  Google Scholar
 R. K. Yadav and R. Upadhayay, “Applications expansion in kernel PCA approach for face recognition,” International Journal of Artificial Intelligence and Knowledge Discovery, vol. 3, no. 4, pp. 1–4, 2013. View at: Google Scholar
 Z. Sun, J. Li, and C. Sun, “Kernel inverse Fisher discriminant analysis for face recognition,” Neurocomputing, vol. 134, pp. 46–52, 2014. View at: Publisher Site  Google Scholar
 J. Haddadnia, K. Faez, and M. Ahmadi, “A fuzzy hybrid learning algorithm for radial basis function neural network with application in human face recognition,” Pattern Recognition, vol. 36, no. 5, pp. 1187–1202, 2003. View at: Publisher Site  Google Scholar
 F. Beverina, G. Palmas, M. Anisetti, and V. Bellandi, “Tracking based face identification: a way to manage occlusions, and illumination, posture and expression changes,” in Proceedings of the 2nd IET International Conference on Intelligent Environments (IE '06), pp. 161–166, July 2006. View at: Publisher Site  Google Scholar
 J. L. Johnson and D. Ritter, “Observation of periodic waves in a pulsecoupled neural network,” Optics Letters, vol. 18, no. 15, pp. 1253–1255, 1993. View at: Publisher Site  Google Scholar
 C. Gao, D. Zhou, and Y. Guo, “Automatic iterative algorithm for image segmentation using a modified pulsecoupled neural network,” Neurocomputing, vol. 119, pp. 332–338, 2013. View at: Publisher Site  Google Scholar
 J. Li, B. Zou, L. Ding, and X. Gao, “Image segmentation with PCNN model and immune algorithm,” Journal of Computers, vol. 8, no. 9, pp. 2429–2436, 2013. View at: Publisher Site  Google Scholar
 I. De and B. Chanda, “Multifocus image fusion using a morphologybased focus measure in a quadtree structure,” Information Fusion, vol. 14, no. 2, pp. 136–146, 2013. View at: Publisher Site  Google Scholar
 J. Lang and Z. Hao, “Novel image fusion method based on adaptive pulse coupled neural network and discrete multiparameter fractional random transform,” Optics and Lasers in Engineering, vol. 52, no. 1, pp. 91–98, 2014. View at: Publisher Site  Google Scholar
 N. Wang, Y. Ma, K. Zhan, and M. Yuan, “Multimodal medical image fusion framework based on simplified PCNN in nonsubsampled contourlet transform domain,” Journal of Multimedia, vol. 8, no. 3, pp. 270–276, 2013. View at: Publisher Site  Google Scholar
 Q. Liu, L.P. Xu, Y.D. Ma, and Y. Wang, “Image NMI feature extraction and retrieval method based on pulse coupled neural networks,” Acta Automatica Sinica, vol. 36, no. 7, pp. 931–938, 2010. View at: Publisher Site  Google Scholar
 J. L. Johnson, “Pulsecoupled neural nets: translation, rotation, scale, distortion, and intensity signal invariance for images,” Applied Optics, vol. 33, no. 26, pp. 6239–6253, 1994. View at: Publisher Site  Google Scholar
 R. C. Nie, S. W. Yao, and D. M. Zhou, “Face recognition using simplified pulse coupled neural network,” Computer Science, vol. 41, no. 2, pp. 297–301, 2014. View at: Google Scholar
 P. Viola and M. J. Jones, “Robust realtime face detection,” International Journal of Computer Vision, vol. 57, no. 2, pp. 137–154, 2004. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Rencan Nie et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.