FACE RECOGNITION ALGORITHM BASED ON ORIENTATION HISTOGRAM OF HOUGH PEAKS

International Journal of Artificial Intelligence & Applications (IJAIA), Vol. 5, No. 5, September 2014
FACE RECOGNITION ALGORITHM BASED
ON ORIENTATION HISTOGRAM OF
HOUGH PEAKS
Amrita Biswas1 and M.K.Ghose2
1Department of Electronics & Communication Engineering, SMIT, Majitar, Rangpo
2Dean Academics, SMIT, Majitar, Rangpo
ABSTRACT
In this paper we propose a novel face recognition algorithm based on orientation histogram of Hough
Transform Peaks. The novelty of the approach lies in utilizing Hough Transform peaks for determining the
orientation angles and computing the histogram from it. For extraction of feature vectors first the images
are divided into non overlapping blocks of equal size. Then for each of the blocks the orientation
histograms are computed. The obtained histograms are combined to form the final feature vector set.
Classification is done using k nearest neighbor classifier. The algorithm has been tested on the ORL
database, Yale B Database & the Essex Grimace Database.97% Recognition rates have been obtained for
ORL database, 100% for Yale B and 100% for Essex Grimace database
KEYWORDS
Face Recognition, Hough Transform Peaks, Orientation Histograms
1. INTRODUCTION
Face Recognition is a challenging field in pattern recognition and computer vision. Several
researchers have come forward with various methods for face recognition [1-2]. However
recognition under challenging conditions like illumination variation, pose variation, aging etc.
still remains an issue. Feature extraction is an essential part of face recognition. Usually a face
image contains lot of redundant information which are not essential for recognition purposes.
Only specific features like eyes, structure of nose, lips, hair line are more than enough to identify
a face image. Also the features extracted should be robust, inavariant to illumination & pose
variations, easy to compute and compact so that storage space required is minimal. Transform
based feature extraction processes aim to decorrelate the pixels and pack as much information as
possible into the smallest number of transform coefficients[3].Various transform based
recognition algorithms have been proposed by researchers[4-7].Transformations like Fourier
Transform, Discrete Wavelet Transform, Radon Transform, Fourier Mellin Tranform etc. have
been investigated for face recognition. Hough Transform has been widely used by researchers for
determining lines and other shapes in pattern recognition problems. Several researchers have
utilized Hough Transform for face detection. Arindam Kar et al [8] used Hough Transform for
investigating its performance in face recognition. However in their work they computed the
Hough Transform of the binary gradient image and in their work selection of significant blocks
requires several iterations. In this paper we investigate the Hough Transform based histogram
DOI : 10.5121/ijaia.2014.5509 107

International Journal of Artificial Intelligence & Applications (IJAIA), Vol. 5, No. 5, September 2014
features for face recognition. Invariance to image scaling, rotation and partial invariance to
illumination is an essential feature of Hough Transform [9]. The results obtained using orientation
histograms of Hough Transform are as good as the results obtained using state of the art methods
used for face recognition.
In the proposed algorithm first the images are divided into non overlapping blocks. Then each
block is binarized and the Hough Transform of each binarized block is computed. Finally the
Hough Peaks have been detected.30 Hough transform peaks were selected from each image block
and then the orientation angle corresponding to the Hough peaks were selected. The orientation
histogram for the selected peaks was computed. The orientation bins range from -90° to
90°.The same procedure is repeated for each of the blocks and to form the final feature vector set
the orientation histograms obtained from each image block are concatenated.
The next sections are divided as follows. In Section 2 we describe the Hough Transform and in
Section 3 we describe the concept of Orientation Histograms. In Section 4 the experiments
performed and results obtained have been discussed and we finally conclude in section 5.
108
2. HOUGH TRANSFORM
Hough Transform is a feature extraction technique widely used in image analysis, computer
vision and digital image processing. Hough Transform performs groupings of edge points into
object candidates by performing an explicit voting procedure over a set of parameterized image
objects and thus helps in selecting the edge features of the face images which play an important
role in face recognition [10]. The standard Hough transform which is used for detecting straight
lines has been used in the proposed algorithm. In the image space, considering a point (xi,yi) in
the xy plane and the general equation of a straight line in slope intercept form yi=axi+b. Plenty of
lines pass through (xi,yi) but they all satisfy the equation yi=axi+b for varying values of a and b.
Let us write this equation as b=-xia+yi and consider the ab plane yields the equation of a single
line for a fixed pair (xi,yi).A second point (xj,yj) also has a line in parameter space associated with
it and unless and until they are parallel this line intersects the line associated with (xi,yi) at some
point (a’,b’) where a’ is the slope and b’ is the intercept of the line containing both (xi,yi) and
(xj,yj) in the xy plane. Precisely, all the points on this line have lines in parameter space that
intersect at (a’,b’).Fig.1 shows this concept. In principle the parameter space lines corresponding
to all points (xk,yk) in the xy plane could be plotted and the principal lines in that plane could be
found by identifying points in parameter space where large numbers of parameter space lines
intersect. However, vertical lines give rise to unbounded values of the parameters a and b.i.e
slope of the line will approach infinity as the line approaches vertical direction. So the other
option is to use the normal representation of the line: x cos+y sin=r
Fig.2(a) illustrates the geometrical interpretation of the parameters r . A horizontal line has
=0°,with r being equal to the positive x intercept. Similarly a vertical line has =90°, with r
being equal to the positive y intercept or =-90°,with r being equal to the negative y intercept.
Each sinusoidal curve in Fig.2(b) represents the family of lines that pass through a particular
point (xk,yk) in the xy plane. The intersection point (r’,’) in Fig.2(b) corresponds to the line that
passes through both (xi,yi) and (xj,yj) in Fig.2(a). The r parameter space is subdivided into
accumulator cells shown in Fig.2(c) where (rmin,rmax) and (min,max) are the expected ranges of the
parameter values : 90 q 90 ° ° − £ £ and −D £ r £ D,where D is the maximum distance between
opposite corners of the image. The cell at coordinates (i,j) with accumulator values A(i,j)

International Journal of Artificial Intelligence Applications (IJAIA), Vol. 5, No. 5, September 2014
,corresponds to the square associated with parameter space coordinates (ri,j).Initially these cells
are set to zero. Then for every non background pixel (xk,yk) in the xy plane is equaled each of
the allowed subdivision values on the axis and solved for the corresponding value of r using the
equation : r=xk cos+yk sin
The resulting r values are then rounded off to the nearest allowed cell value along the r axis. If a
choice of p results in solution rq, then A(p,q)=A(p,q)+1. After the procedure is over a value of P
in A(i,j) means that P points in the xy plane lie on the line x cosj+y sinj=ri. The number of
computations in this method is linear with respect to the number of non background points in the
xy plane.
109
An approach based on Hough Transform is as follows:
• Obtain binary edge image
• Specify subdivisions in the r plane.
• Examine the counts of the accumulator cells for high pixel concentrations.
• Examine the relationship (principally for continuity) between pixels in a chosen cell.
(a) (b)
Fig. 1 (a) xy plane (b) Parameter Space
(a)

110
(b) (c)
Fig.2. (a)(r,) parameterization of line in xy plane (b) Sinusoidal curves in the (r,) plane;the point of
intersection (r’,’) corresponds to the line passing through points (xi,yi) and (xj,yj) in the xy plane (c)
Division of the (r,) plane into accumulator cells
3.ORIENTATION HISTOGRAMS
The histogram of a digital image with gray levels in the range [0, L-1] is a discrete function h (rk)
= nk, where rk is the kth gray level and nk is the number of pixel in the image having gray level rk.
It is common practice to normalize a histogram by dividing each of its values by the total number
of pixels in the image, denoted by n. Thus, a normalized histogram is given by p(rk) = nk/n, for k
= 0, 1……, L-1. In other words, p(rk) gives an estimate of probability of occurrence of gray level
rk. The motivation behind the Orientation Histogram descriptors is that local object appearance
and shape within an image can be described by the distribution of intensity gradients or edge
directions. These descriptors can be computed by dividing the image into small blocks, and for
each block computing a histogram of the orientation .The concatenation of these histograms will
form the descriptor. To improve the efficiency, the local histograms can be normalized by
calculating a measure of the intensity across a larger region of the image, and then using this
value to normalize all blocks within that region. Invariance to changes in illumination or
shadowing improves due to this normalization. The major advantages of this descriptor will be
invariance to geometric and photometric transformations and hence the feature vectors obtained
from these descriptors are well suited for application in face recognition problems.
4. PROPOSED APPROACH
• The images are partitioned into equal sized non overlapping blocks
• The gray level block are converted into binary form
• The Hough Transform of each image block is computed.
• The Hough Peaks are detected and the first 30 Hough peaks are retained
• The orientation angle for the retained peaks are selected
• Based on the orientation histograms are computed for the image blocks
• The histograms of all the blocks were concatenated to form the final feature vector set

111
4. EXPERIMENTS PERFORMED RESULTS OBTAINED
For our experiments we have tested the ORL, Yale B Essex Grimace Database. The Databases
are described as follows:
ORL (AT and T) database has 10 different images (92 x 112), each of 40 different subjects. All
images were taken against a dark homogenous background with the subjects in upright, frontal
position with some side movement. Sample images of the ORL database have been shown in Fig.
2. Yale Database contains images of 10 subjects (480x 640) under 64 different lighting
conditions. Sample images of this database are shown in Fig. 3. Essex-Grimace database has a set
of 20 different images (180x200), each of 18 individuals. The database contains images of both
male and female subjects, taken with a fixed camera. Sample images of the Essex Grimace
database have been shown in Fig. 4.
All the images from the three databases were resized to 112x112.Then the images were divided
into blocks of size 16x16.So a total of 49 blocks were obtained for each image. Each of the gray
level blocks was converted to binary and the Hough Transform for each image block was
computed. The size of the obtained Hough Transformed matrix for each image block of size
16x16 is 45x1796 .Next the Hough Peaks were detected and a total of 30 peaks were retained for
each image block. The orientation angle for the 30 peaks were selected and the orientation
histogram computed for the image blocks. The orientation bins range from -90° to 90°.Hence the
size of the histogram obtained for each block is 1x181.The orientation histograms thus obtained
were concatenated to form the final feature vector set of the given image. Hence the size of the
feature vector set for a given image is 1x8869
Fig.2.Sample Images of the ORL Database
Fig.3.Sample Images of the Yale B Database

112
Fig.4.Sample Images of the Essex Grimace Database
For each of the databases we have taken different numbers of training images and the remaining
images were selected for testing. k nearest neighbor classifier was used for classification.
Fig.5(a) Fig.5(b) Fig.5(c)
Fig.5(a) Original Image (b) Image Partitioned into blocks of size 16x16 (c) Binarized blocks of
Image
Fig.6(a) Fig.6(b)

113
Fig.6(c) Fig.6(d)
Fig.6(a) One of the image blocks of size 16x16 (b) Corresponding binarized image block (c)
Hough Transform image the 30 Peaks extracted from the same block (d) Orientation Histogram
of Hough Peaks for the same block
Results obtained have been shown in Table.1.
Table 1: Results obtained
Sl.No. Database No.of
Training
Images
Recognition %
1 ORL 5 97
2 ORL 4 96
3 ORL 3 92
4 ORL 2 90
5 ORL 1 80
6 Yale B 4 100
7 Yale B 3 97
8 Yale B 2 96
9 Yale B 1 94
10 Essex-
Grimace
1 100
Table 2.Comparison with other transform based algorithms for single training image
Sl.No Algorithm Database Recog%
1
DWT based Algorithm using coiflets[4] ORL 71.2
Yale 67.5
Essex Grimace 85.3
2 Radon Transform based Algorithm[4] ORL 67.3
Yale 70
Essex Grimace 97
3 FMT based Algorithm[5] ORL 68
Yale 65
Essex Grimace 98
4 Proposed Orientation Histogram based Algorithm ORL 80
Yale 94
Essex Grimace 100

114
5.CONCLUSION
From the results obtained it can be seen that the orientation histograms obtained from the peaks of
the Hough Transformation matrix does a good job in minimizing the effects of illumination
variance. The algorithm performs well on the Yale B database which consists of images with
major illumination variance. It also takes care of major expression variations which can be
corroborated by the results obtained on the Essex Grimace Database. The orientation histograms
also take care of the pose variations in the ORL database and it can be seen that as we increase
the number of images in the training database the performance of the algorithm improves. The
best performance of the proposed approach for the three databases is almost similar to the
transform based algorithms proposed in [4] and [5].But it can be said that when the number of
training images reduce there is significant reduction in performance of the transformed based
algorithms but the Orientation Histogram based algorithm performs better for less number of
training images. One of the reasons could be the orientation histograms better depict the rotation
variance of the images. Also, while the transform based approaches were holistic, the Orientation
Histogram based approach is local feature based. The histograms contain information about the
distribution of the local orientation over the whole image. But for efficient face representation
spatial information also has to be retained. To retain spatial information the image is has been
partitioned into non overlapping blocks and combination of the histograms of these blocks
provide the spatially enhanced histogram. The results reflected in the comparison table provide
inspiration to investigate local appearance based methods for designing robust algorithms using
single training image.
REFERENCES
[1] Rabia Jafri and Hamid R. Arabnia,2009,A survey of Face Recognition Techniques, Journal of
Information Processing Systems,Vol.5,No.2,pp-41-68
[2] Ripal Patel, Nidhi Rathod, Ami Shah, 2012,Comparitive Analysis of Face Recognition Approaches:A
Survey,Vol:57,No.17,pp-50-69
[3] R. C. Gonzalez and R. E. woods,” Digital Image Processing”, Dorling Kindersley, India, Pearson
Prentice Hall, 2006.
[4] Sambhunath Biswas,Amrita Biswas,2012,Face Recognition Algorithms Based on Transformed Shape
Features, IJCSI International Journal of Computer Science Issues, Vol. 9, Issue 3, No 3,pp-445-451
[5] Sambhunath Biswas and Amrita Biswas,2013,FourierMellin Transform Based Face
Recognition,International Journal of Computer Engineering and Technology(IJCET),vol.4,issue 1,
pp-8-15
[6] Anissa Bouzalmat, Arsalane Zarghili and Jamal Kharroubi,2011, Facial Face Recognition Method
using Fourier Transform Filters Gabor and R_LDA, IJCA Special Issue on Intelligent Systems and
Data Processing, pp-18-24.
[7] Zheng Dezhong Cui Fayi, ‘Face Recognition based on Wavelet Transform and Image
Comparison”,2008,Proc. International Symposium on Computational Intelligence and Design,
Volume: 2, pp-24-29.
[8] Arindam Kar1, Debotosh Bhattacharjee2, Dipak Kumar Basu2, Mita Nasipuri2, Mahantapas
Kundu2,2012,Face Recognition using Hough Peaks extracted from the significant blocks of the
Gradient Image, International Journal of Advanced Research in Computer Science and Software
Engineering,Vol.2,Issuue 1,pp-1-6.
[9] http://homepages.inf.ed.ac.uk/rbf/HIPR2/hough.htm.
[10] Shapiro, Linda and Stockman, George. Computer Vision, Prentice-Hall, Inc. 2001
[11] Use of the Hough Transformation to Detect Lines and Curves in Pictures (http:/ / www. ai. sri. com/
pubs/ files/ tn036-duda71. pdf).
[12]R.C Gonzalez and R.E Woods ,Digital Image Processing,Dorling Kinderlsey,India,Pearson,Prentice
Hall,2006

FACE RECOGNITION ALGORITHM BASED ON ORIENTATION HISTOGRAM OF HOUGH PEAKS

Recommended

More Related Content

What's hot (19)

Similar to FACE RECOGNITION ALGORITHM BASED ON ORIENTATION HISTOGRAM OF HOUGH PEAKS (20)

Recently uploaded (20)

FACE RECOGNITION ALGORITHM BASED ON ORIENTATION HISTOGRAM OF HOUGH PEAKS