 Technical advance
 Open Access
 Published:
paraFaceTest: an ensemble of regression treebased facial features extraction for efficient facial paralysis classification
BMC Medical Imaging volume 19, Article number: 30 (2019)
Abstract
Background
Facial paralysis (FP) is a neuromotor dysfunction that losses voluntary muscles movement in one side of the human face. As the face is the basic means of social interactions and emotional expressions among humans, individuals afflicted can often be introverted and may develop psychological distress, which can be even more severe than the physical disability. This paper addresses the problem of objective facial paralysis evaluation.
Methods
We present a novel approach for objective facial paralysis evaluation and classification, which is crucial for deciding the medical treatment scheme. For FP classification, in particular, we proposed a method based on the ensemble of regression trees to efficiently extract facial salient points and detect iris or sclera boundaries. We also employ 2^{nd} degree polynomial of parabolic function to improve Daugman’s algorithm for detecting occluded iris boundaries, thereby allowing us to efficiently get the area of the iris. The symmetry score of each face is measured by calculating the ratio of both iris area and the distances between the key points in both sides of the face. We build a model by employing hybrid classifier that discriminates healthy from unhealthy subjects and performs FP classification.
Results
Objective analysis was conducted to evaluate the performance of the proposed method. As we explore the effect of data augmentation using publicly available datasets of facial expressions, experiments reveal that the proposed approach demonstrates efficiency.
Conclusions
Extraction of iris and facial salient points on images based on ensemble of regression trees along with our hybrid classifier (classification tree plus regularized logistic regression) provides a more improved way of addressing FP classification problem. It addresses the common limiting factor introduced in the previous works, i.e. having the greater sensitivity to subjects exposed to peculiar facial images, whereby improper identification of initial evolving curve for facial feature segmentation results to inaccurate facial feature extraction. Leveraging ensemble of regression trees provides accurate salient points extraction, which is crucial for revealing the significant difference between the healthy and the palsy side when performing different facial expressions.
Background
Facial paralysis (FP) or facial nerve palsy is a neuromotor dysfunction that losses voluntary muscles movement in one side of the human face. As a result, this leads to the loss of the person’s ability to mimic facial expressions. FP is not an uncommon medical condition. For every sixty people around the world, one of them can be affected by facial paralysis [1]. Facial paralysis often causes patients to be introverted and eventually suffer from social and psychological distress, which can be even more severe than the physical disability [2]. It is usually encountered in clinical practices, which can be classified into peripheral and central facial palsy [3, 4]. These two categories differ from each other according to the behavior of the upper layer of the human face. Peripheral facial palsy is a nerve disturbance in the pons of the brainstem, which affects the upper, middle and lower facial muscles of one side of the face. On the other hand, central facial palsy is a nerve dysfunction in the cortical areas whereby the forehead and eyes are spared, but the lower half of one side of the face is affected, unlike in peripheral FP [3–5]. Such scenario has triggered the interest of researchers and clinicians of this field, and, consequently, led them to the development of objective grading facial functions and methods in monitoring the effect of medical, rehabilitation or surgical treatment.
Many computeraided analysis systems have been introduced to measure dysfunction of one part of the face and the level of severity, but not much on facial paralysis type as the classification problem. Moreover, the efficiency of the method used for it to be universally accepted is still underway [3]. Classification of each case of facial paralysis into central or peripheral plays a critical role in helping physicians to decide for the most appropriate treatment scheme to use [4]. Image processing has been applied in the existing objective facial paralysis assessment, but the processing methods used are mostly laborintensive, if not; suffer from the sensitivity to the extrinsic facial asymmetry caused by orientation, illumination and shadows. As such, creating a clinically usable and reliable method is still very challenging. Wachtman et al. [6, 7] measured facial paralysis by examining the facial asymmetry on static images. Their methods are prone to facial asymmetry even for healthy subjects due to sensitivity to orientation, illumination and shadows [8]. Other previous works [9–11] were also introduced but most of them are solely based on finding salient points on human face with the use of the standard edge detection tool (e.g. Canny, Sobel, SUSAN) for image segmentation.
Canny edge detection algorithm may yield to inaccurate results as it influences connected edge points. This is because it does comparisons of the adjacent pixels on the gradient direction to determine if the current pixel has a local maximum. Improper segmentation will also result to improper generation of key points. Moreover, it may be difficult to find and detect the exact points when the same algorithm is applied to elder patients. Dong et al. [11] proposed the use of salient points for degree of facial paralysis estimation. A method was proposed to detect the salient points that will be the basis for the estimation of the degree of paralysis. As the salient points may include some unnecessary points for describing facial features, edge detection was used to discard these points. Then Kmeans clustering is applied to classify these salient points into six categories: 2 eyebrows, 2 eyes, nose, and mouth. There are about 14 key points are found in the six facial regions, which would represent the points that may be affected when performing some facial expressions. However, this technique falls short when applied to elder patients, in which exact points can be difficult to apply [12].
Another method was proposed that estimates the degree of facial paralysis by comparing multiple regions on human face. In [13], Liu et al. perform comparison of the two sides of the face and compute the four ratios, which are consequently used to represent the severity of the paralysis. Nevertheless, this method suffers from the influence of uneven illumination [8].
In our previous work [4], a technique that generates closed contours for separating outer boundaries of an object from background using Localized Active Contour (LAC) model for feature extraction was employed, which reasonably reduced these drawbacks. However, one limiting factor of the method introduced is that it has a greater sensitivity to subjects exposed to peculiar images where facial features suffer from different occlusions (e.g. eyeglasses, eyebrows occluded with hair, wrinkles, excessive beard, etc.). Moreover, improper setting or identification of parameters, such as the radius of the initial evolving curve (e.g. minimumbounding box of the eyes, eyebrows, etc.) may lead to improper features segmentation, which may in turn gives inaccurate detection of key points as revealed in Fig. 1. Although such limitation was addressed in [4] by applying a window kernel (i.e. generated based on Otsu’s method [14]) that is run through the binary image (e.g. detected eyes, eyebrows, lip, etc.), Otsu’s method has some drawbacks of the assumption of uniform illumination. Additionally, it does not use any object structure or spatial coherence, which may sometimes result to inaccurate generation of kernel paramaters that may in turn result to improper segmentation.
In order to address these drawbacks, we present a novel approach based on ensemble of regression trees (ERT) model. We leverage ERT model to extract salient points as basis for preparing the features or independent variables. Our target was to keep a natural framework that finds 68 key points on edges of facial features (e.g. eyes, eyebrows, mouth, etc.); and regresses the location of the facial landmarks from a sparse subset of intensity values extracted from an input image. As ERT model provides a simple and intuitive way in performing this task and in generating landmarks that separates boundaries of each significant facial feature from the background, without requiring settingup of parameters for initial evolving curves as required by the previous approach [4], we find this technique appropriate to address our problem. Furthermore, empirical results [15] reveal that ERT model does have an appealing quality that performs shape invariant feature selection while minimizing the same loss function during training and test time, which significantly lessens time complexity of extracting features.
In this study, we make three significant contributions. First, we introduce a more robust approach for efficient evaluation of facial paralysis and classification, which is crucial for determining the appropriate medical treatment scheme. Second, we provide an efficient method in feature extraction based on ensemble of regression trees in a computationally efficient way; and finally, we study in depth, the effect of computing the facial landmarks symmetry features on the quality of facial paralysis classifications.
Methods
The study performs objective evaluation of facial paralysis particularly the facial classification and grading in a computationally efficient way. We capture facial images (i.e. still photos) of the patients with a front view face and with reasonable illumination of lights so that each side of the face achieves roughly similar amount of lighting. The photo taking procedure starts with the patient performing the ’at rest’ face position, followed by the three voluntary facial movements that include raising of eyebrows, screwingup of nose, and showing of teeth or smiling.
The framework of the proposed objective facial paralysis evaluation system is presented in Fig. 2. Facial images of a patient, which are taken while being requested to perform some facial expressions, are stored in the image database. In the rest of this section, we describe the details of the form of individual components of the facial landmark detection and how we perform evaluation and classification of facial paralysis. To begin the process, raw images from the image database are extracted or retrieved. Dimension alignment and image resizing are then performed, followed by the preprocessing of images to enhance contrast and remove undesirable image noise.
Facial Image Acquisition and Preprocessing
Existing literatures show that facial data acquisition methods can be classified into two categories depending on the processing methods they used and on the usage of images or videos as their databases. As imagebased systems do have an appealing advantage of ease of use and cost effectiveness [16], we utilize still images as inputs to our classification model.
Face Normalization
Facial feature extractions can be complicated as the face appearance changes, which is caused by illumination and pose variations. Normalizing the acquired facial images prior to objective FP assessment may significantly reduce these complications. It is worth note taking, however, that while face normalization may be a good approach in connection with the methods for objective facial paralysis assessment, it is but optional, so long as extracted feature parameters are normalized prior to the classification task.
Facial Feature Extraction
Deformation Extraction
Facial features deformation can be categorized by the changes of texture and shape that lead to high spatial gradients, which are good indicators for tracking facial actions [17]. In turn, these are good indicators for facial asymmetry and may be analyzed either in image or spatial frequency domain, which can be computed by highpass gradient or Gabor waveletbased filters. Ngo et al. [18] utilized this method by combining it with local binary patterns and claimed to perform well for the task of quantitative assessment of facial paralysis. Deformation extraction approaches can be holistic or local. Holistic imagebased approaches are those where face is processed as a whole, while local methods focused on facial features or areas that are proned to change and are able to detect subtle changes in small areas [17]. In this work, we extract the facial salient points and iris area for geometricbased and regionbased features generation, respectively.
Salient points detection
Feature extraction starts from the preprocessing of the input image and facial region detection. We find frontal human faces in an image and estimate their pose. Such estimated pose takes the facial form of 68 landmarks that will lead to the detection of the key points of the mouth, eyebrows and eyes. The face detector is made using the classic Histogram of Oriented Gradients (HOG) feature [19] combined with a linear classifier, an image pyramid, and sliding window detection scheme. We utilize the pose estimator employed in [15]. In the context of computer vision, a sliding window is a rectangular region of fixed width and height that slides across an image, as described in the subsequent section. Normally, for each window, we take the window region and apply an image classifier to check if the window has an object (i.e. the face as our region of interest). Combined with image pyramids, image classifiers can be created that can recognize objects even if the scales and locations in the image vary. These techniques, in its simplicity, play an absolutely critical role in object detection and image classification.
Features are extracted from the different facial expressions that the patient will perform, which include: (a) at Rest; (b) Raising of eyebrows; (c) frowning or screwing up of nose; and (d) Smiling or showing of teeth. In this study, we consider geometric and regionbased features as our inputs for modelling a classifier. To extract the geometricbased features, the salient points of the facial features (e.g. eyes, eyebrows and lip) are detected using the concept of ensemble of regression trees [15]. The goal is to find the 68 key points on edges of facial features as shown in Fig. 3a. However, for generating the geometricbased symmetry features of each image, we are interested in the following points: inner canthus (IC), mouth angle (MA), infra orbital (IO), upper eyelids (UE), supraorbital (SO), nose tip (NT) and nostrils (see Fig. 3b). In order to detect the salient points of the facial image, we leverage the dlib library [15] that has the capability to detect points of each facial features. These points are inputs for calculating the distance of the key points as identified in Fig. 4, and the corresponding distance ratio for both sides of the face. Additionally, we also extract regionbased features, which involves detection of iris or sclera boundaries by leveraging Daugman’s IntegroDifferential Operation [20].
Histogram of Oriented Gradients(HOG)
In the field of computer vision and image processing, Histogram of Oriented Gradients (HOG) is a feature descriptor used for object or face detection. It is an algorithm, which takes an image and outputs feature vectors or feature descriptors. Feature descriptors encode interesting information into a series of numbers and act as a sort of numerical fingerprint, which can differentiate one feature from another. Ideally, this information is invariant under image transformation; hence, even if the image is transformed in some way, we can still find the feature again. HOG uses locally normalized histogram of gradient orientation features similar to Scale Invariant Feature Transform (SIFT) descriptor in a dense overlapping grid, which gives very good results in face detection. It is similar to SIFT, except that HOG feature descriptors are computed on a dense grid of uniformly spaced cells and they use overlapping local contrast normalization for improved performance [19, 21].
Implementation of the HOG descriptor algorithm is as follows:

1.
Partition the image into small connected regions called cells, then for each cell, calculate the gradient directions histogram for the pixels within the cell.

2.
According to the gradient orientation, discretize each cell into angular bins. Each pixel of the cell contributes weighted gradient to its corresponding angular bin. Adjacent group of cells are considered as spatial regions, referred to as blocks. The grouping of cells into a block is then the basis for histograms grouping and normalization. Normalized group of histograms represents the block histogram. The set of these block histograms represents the descriptor.
The following basic configuration parameters are required for computation of the HOG descriptor:

1.
Masks to compute derivatives and gradients.

2.
Geometry of splitting an image into cells and grouping cells into a block.

3.
Block overlapping

4.
Normalizing parameters
However, the recommended values for the HOG parameters include: (a)1D centered derivative mask [1, 0, +1]; (b) Detection window size of 64x128; (c) Cell size of 8x8; and Block size of 16x16 (2x2 cells) [21].
Ensemble of Regression Trees
An ensemble of regression tree is a predictive model, which composes a weighted combination of multiple regression trees. Generally, combining multiple regression trees increases the predictive performance. It is with the collection of regression trees that makes a bigger and better regression tree [15]. The core of each regression function r_{t} is the treebased regressors fit to the residual targets during the gradient boosting algorithm.
Shape invariant split tests
Based on the approach used by Kazemi et al. [15], we make a decision based on thresholding the difference between the intensities of two pixels at each split node in the regression tree. When defined in the coordinate system of the mean shape, the pixels used in the test are at positions i and k. For facial images having arbitrary shapes, we intend to index the points that have the same position relative to its shape, as i and k have to the mean shape. To accomplish this task, the image can be warped to the mean shape based on the current shape estimate before extracting the features. Warping the location of points is more efficient rather than the whole image, since we only utilize a very sparse representation of the image. In what follows, the details are precisely shown. We let vi be the index of the facial landmark in the mean shape that is nearest to i and its offset from i is defined as:
Then for a shape S_{j} defined in image I_{j}, the position in I_{j} that is qualitatively similar to i in the mean shape image is given by
where S_{j} and R_{j} are the scale and rotation matrix of the similarity transform which transforms S_{j} to \(\overline {S}\), the mean shape. The scale and rotation are found to minimize
the sum of squares between the mean shape’s facial landmark points, \(\overline {Y}_{m}\)’s, and those of the warped shape. k^{′} is similarly defined. Formally each split is a decision involving 3 parameters θ=(τ,i,k) and is applied to each training and test example as
where i^{′} and k^{′} are defined using the scale and rotation matrix which best warp \(\hat {S}^{(t)}_{j}\) to \(\overline {S}\) according to equation (2).
In practice, during the training phase, the assignments and local translations are identified. Computing the similarity transform, at test time, the most computationally expensive part of the process is only done once at every level of the cascade.
This method starts by using the following:

1.
A training set of labeled facial landmarks on an image. These images are manually labeled, specifying specific (x, y)coordinates of regions surrounding each facial structure.

2.
Priors, or more specifically, the probability on distance between pairs of input pixels.
Given this training data, ensemble of regression trees are trained to estimate the positions or locations of facial landmarks directly from the pixel intensities, that is, no feature extraction is taking place. The final result is a detector of facial landmarks that can be utilized to efficiently detect salient points of an image. Fig. 5 presents sample results of our proposed approach for facial features detection based on ensemble of regression trees (ERT) model when applied in JAFFE dataset [22, 23].
Regionbased Feature Extraction
Iris Detection
A person who has a symptom of facial paralysis is likely to have asymmetric distance between the infra orbital (IO) and mouth angle (MA) of both sides of his face while performing facial movements such as smile and frown. They may also have an unbalanced exposure of iris when performing different voluntary muscle movements (e.g. screwing of nose, showing of teeth or smiling, raising of eyebrows with both eyes directed upward) [4]. We utilize the points (e.g. upper eyelid (UE), infra orbital (IO), inner canthus (IC) and outer canthus (OC)) detected by ensemble of regression trees algorithm and from there, we generate the parameters of the eye region as inputs to Daugman’s Integro Differential operator [20] to detect the iris sclera or boundaries. However, as some eye images do have eyelid occlusions, optimization of Daugman’s algorithm was performed to ensure proper iris boundaries detection. In our previous work, LACbased method was employed to optimize Daugman’s algorithm and perform subtraction method. However, LAC model is quite tedious as it may result to improper segmentation if the initial evolving curves and iterations are not properly defined. Moreover, it has a greater sensitivity to the eyes of old patients due to the presence of excessive wrinkles.
In this paper, we implement the curve fitting scheme to ensure proper detection of iris sclera or boundaries, thereby providing better results in detecting occluded iris. By definition, a parabola is a plane curve. It is said to be mirrorsymmetrical and is approximately Ushaped. To fit a parabolic profile through the data, a second degree polynomial is used, defined as: y=ax^{2}+bx+c
This will exactly fit a simple curve to the three constraints, which include a point, angle, or curvature. More often, the angle and curvature constraints are added to the end of a curve. In this case, they are called end conditions. To fit a curve through a data, we first get the coefficients b and c of the line p(x) = bx + c. To evaluate this line, values of x (e.g. given by xp) must be chosen. For example, the curve will be plotted for x belongs to [0, 7] in steps of delta x = 0.1. The generated coefficients will then be used to generate the y values of the polynomial fit at the desired values of x given by xp. This means that the vector (denoted as yp) can be generated, which contains the parabolic fit to the data evaluated at the xvalues xp.
In what follows, we describe our approach for detecting the iris sclera or boundaries:

1.
Detect eye region (i.e. creating a rectangular bounding box from the detected UE,IO, IC, OC as illustrated in Figure 3b)

2.
Detect upper eyelid edge using gray thresh level (e.g. threshold <= 0.9)

3.
Convert generated grayscale image to binary.

4.
Find significant points of the upper eyelid.

Traverse each column and find the first row whose pixel value=1 and whose location variance (i.e. row address of the last 4 pixels) is minimal within the threshold. We call it vector A.


5.
Implement curve fitting through the generated data points using parabolic profile (i.e. second degree of polynomial). We refer this as A’.

6.
Detect iris using Daugman’s integrodifferential operator and convert it to binary form. We call it vector B.

7.
Find intersections of the two vectors A’ and B and take all pixels of vector B below the parabolic curve A’.

8.
Utilize the points of the lower eyelid detected ERT model, we call it vector C.

9.
Finally, we find the the region of the detected iris within the intersection of vector B and C.
A closer look at Figs. 6 and 7 reveal interesting results of this approach.
Facial Paralysis Measurement
Symmetry Measurement by Iris and Key points
In this paper, the symmetry of both sides of the face is measured using the ratio that are obtained from computing the iris area (i.e. generated while the subject performs raising of eyebrows while both eyes are directed upward; and screwing of nose or frown) and the distance between the two identified points on each side of the face while the subject is asked to perform the different facial expressions (e.g. at rest, raising of eyebrows, screwing of nose, and showing of teeth or smile). Table 1 shows the summary of the salient points as basis for extracting features such as the ratio of iris area as well as the distance ratio between two sides of the face.
With ’at rest’ and ’raising of eyebrows’, we calculate the distance between two points: Infra Orbital (IO) and supraorbital (SO); and the Nose Tip (NT) and SO, while with ’smile’ expression, we get the distances between the two identified points: inner canthus (IC) and Mouth Angle (MA); IO and MA; and the NT and MA. Lastly, for ’frown’ expression, we get the distance between the two points: NT and MA; and NT and nostrils. Consequently, the computed distance ratio of both sides of the face are considered as the symmetrical features of each subject.
Computed distances include: P_{20}P_{R_IO}, P_{25}P_{L_IO}, P_{20}P_{31} and P_{25}P_{31} (see Fig. 4a and 4b); P_{31}P_{33}, P_{31}P_{35}, P_{31}P_{49} and P_{31}P_{55} (see Figure 4c); and P_{40}P_{49}, P_{43}P_{55}, P_{31}P_{49}, P_{31}P_{55}, P_{49}P_{R_IO} and P_{55}P_{L_IO} (see Figure 4d). Additionally, we calculate the area of the extracted iris. This is followed by the computation of ratio between two sides. We find the expression below:
where dRatio is the ratio of the computed distance D_{L} and D_{R} of the specified key points of each half of the face.
We also consider the capability of the patients to raise the eyebrows (i.e. rate of movement) as one important features for symmetry measurement by comparing the two facial images, the ‘at rest position’ and ‘raising of eyebrows’ as shown in Fig. 4a and 4b. We compute the distances a_{1} and b_{1} (Fig. 4a) where a_{1} and b_{1} are the distance from P_{20} and \(P_{\texttt {R\_IO}}\) and P_{20} and P_{31} of the right eye, respectively. We then compute the ratio of a_{1} and b_{1}. Similarly, for the second image (Fig. 4b), we get a_{2} and b_{2} as well as the ratio. Finally, we compute the difference of these two ratio (i.e. difference between a_{1}/ b_{1} and a_{2}/ b_{2}) and denote it as right_Movement.
The same procedure is applied to the two images for finding the ratio difference for the left eye (i.e. difference between y_{3}/ x_{3} and y_{4}/ x_{4}) and denote it as left_Movement. The rate of movement can be computed by finding the ratio between right_Movement and left_Movement. Intuitively, the difference of these two ratio values for normal subjects is likely to be higher (usually approaching to 1, which signifies the ability to raise both of his eyebrows) than those of the FP patients.
Facial Palsy Type Classification
Classification of facial paralysis type involves two phases: (1) discrimination of healthy from unhealthy subjects; and (2) proper facial palsy classification. In this context, we model the mapping of symmetry features (as described in the previous subsection) into each phase as binomial classification problem. As such, we employ two classifiers to be trained, one for healthy and unhealthy discrimination (0healthy, 1unhealthy) and another one for facial palsy type classification (0peripheral palsy(PP), 1central palsy(CP)). For each classifier, we consider Random Forest (RF), Regularized Logistic Regression (RLR), Support Vector Machine (SVM), Decision Tree (DT), naïve bayes (NB) and Hybrid classifier as appropriate classification methods as they have been successfully applied to pattern recognition and classification on datasets with realistic size [4, 8, 13, 22]. With hybrid classifier, we apply rulebased approach prior to employing machine learning (ML) method. This is based on empirical results [4, 11, 13], which show that normal subjects are likely to have an average facial measurement ratio closer to 1.0 and central palsy patients are likely to have a distance ratio from Supraorbital (SO) to Infraorbital (IO) approaching to 1.0. Similarly, iris exposure ratio usually results to values nearly close to 1. Hence, we find hybrid classifier (rulebased + ML) appropriate in our work.
This process is presented in Fig. 8. If rule number 1 is satisfied, the algorithm continues to move to the case path (i.e. for the second task), making a test if rule number 2 is also satisfied; otherwise, it performs a machine learning task, such as RF, RLR, SVM, and NB. The rules are generated after fitting the training set to the DT model. For example, rule 1 may have conditions, like if f_{x}<0.95andf_{y}<0.95 (where f_{x}andf_{y} are two of the predictors, the mean result of all parameters and the IO_MA ratio, respectively, based on Table 1), then the subject is most likely to be diagnosed with facial paralysis, and therefore can proceed for rule no. 2 (i.e. to predict the FP type); otherwise, it performs a machine learning task. If the classifier returns 0, the algorithm exits from the entire process as this signifies that the subject is classified as normal/healthy, else it moves to the case path to perform a test on rule number 2 for facial palsy type classification.
If rule number 2 is satisfied, the system gives 1 (i.e. 0PP; 1CP), else the feature set is fed to another classifier, which can yield either 0 or 1. The same with rule number 1, rule number 2 is also generated by DT model. For example, rule 2 may set conditions like, if f_{a}>0.95 and f_{b}>0.95 (where f_{a} and f_{b} are two of the predictors SO_IO ratio and iris area ratio, respectively), then it is most likely to be diagnosed as having central palsy (CP), otherwise, the feature set is fed to another classifier, which could return either 0 or 1 (i.e. 0PP; 1CP).
Results
In our experiments, we used 440 facial images, which were taken from 110 subjects (i.e. 50 patients and 60 healthy subjects). From the 50 unhealthy subjects, 40 of which have peripheral palsy (PP) cases and 10 have central palsy (CP) cases. We used 70% of the dataset as the training set and 30% as the test set. For example, in discriminating healthy from unhealthy subjects, we used 77 subjects (i.e. 35 patients plus 42 normal subjects) as the training set and 33 subjects (i.e. 15 patients plus 18 normal) as the test set. While in FP type classification problem 35 unhealthy cases (i.e. 28 PP and 7 CP) as our training set and 15 (i.e. 12 PP and 3 CP) as our test set.
Each subject was asked to perform 4 facial movements. During image preprocessing, resolutions were converted to 960 x 720 pixels. Facial palsy type of each subject was prelabeled based on the clinicians’ evaluation, which was used during the training stage. This was followed by the feature extraction or the calculation of the area of the extracted iris and the distances between the identified points as presented in Fig. 4. Overall, we utilize 11 features to train the classifier. The samples for healthy and unhealthy classifier were categorized into two labels: 0healthy, 1unhealthy. Similarly, samples for unhealthy subjects were classified into two labels: 0central palsy and 1peripheral palsy. It can be noted that healthy subjects have very minimal asymmetry in both sides of the face resulting to a ratio that approaches to 1.
Facial paralysis type classification
Regularized logistic regression (RLR), Support Vector Machine (SVM), random forest (RF), naïve bayes (NB), and classification tree (DT) were also utilized to compare with our hybrid classifiers. Given the dataset that is not very large, we adopt the kfold crossvalidation test scheme in forming a hybrid model. The procedure involves 2 steps: rule extraction and hybrid model formation, as applied in our previous work [4].
Step 1: rule extraction.
If we have a dataset D = ((x_{1}, y_{1}),..., (x_{n}, y_{n})), we hold out 30% of D and use it as a test set T = ((x_{1}, y_{1}),..., (x_{t}, y_{t})), leaving 70% as our new dataset D’. We adopt kfold crossvalidation test scheme over D’, i.e. with k = 10. For example, if N = 80 samples, each fold has 8 samples. In each iteration, we leave one fold out and consider it as our validation set and utilize the remaining 9 folds as our training set to learn a model (e.g. rules extraction). Since we have 10 folds, we do this process for 10 repetitions. We extract rules by fitting the training set to a DT model.
Step 2: hybrid model formation.
In this step, we form a hybrid model by combining the rules generated in each fold and the ML classifier. Each model is tested out over the validation set using different parameters (e.g. lambda for logistic regression and gamma for SVM). For example, to form the first hybrid model, we combine the rule extracted from the first fold and a regularized logistic regression model (i.e. rule + RLR) and test out its performance over the validation set (leftout fold) while applying it to each of the 10 parameters. Therefore, for each fold, we generate 10 performance measures. We repeat this procedure for the succeeding folds, which means performing the steps for k times (i.e. with k = 10, since we are using 10fold cross validation) will give us 100 performance measures. We calculate the average performance in all folds. This yields 10 average performance measures (i.e. for each parameter n), each of which corresponds to one specific hybrid model. Then we choose the best hybrid model, i.e. a model with lambda that minimizes errors. We retrain the selected model on all of D’, test this out over the hidden test set T = ((x_{1},y_{1}),·,(x_{t},y_{t})), i.e. 30% of the dataset D and get the performance of the hybrid model.
We evaluate the classifiers by their average performance for 20 repetitions of kfold crossvalidation using k = 10. We repeat this process for evaluation of other hybrid model (e.g. rulebased + RF, rulebased + RLR, etc.) and finally choose the hybrid model that performs best. The hybrid classifiers, RF, SVM, RLR, DT and NB were tested and experiments reveal that hybrid classifier rulebased + RLR (hDT_RLR) outperformed other classifiers for discriminating healthy from unhealthy (i.e. with paralysis) subjects. Similarly, for the classification task of facial palsy (PPperipheral and CPcentral palsy), hDT_RLR hybrid classifier is superior among other classifiers used in the experiments.
Table 2 presents a comparison of the average performance of our hybrid classifier, RF, RLR, SVM, DT and NB based on our proposed approach. For FP type classification, our hDT_RLR hybrid classifier achieves a better performance on sensitivity of 5.2% higher than RLR, RF, SVM, DT and NB (see Table 2). Other hybrid classifiers also show good results comparable with hDT_RLR. However, in the context of our study, we are more concerned on designing a classifier that yields a stable results on sensitivity performance measure without necessarily sacrificing the specificity or fallout or the probability of false alarm. Hence, we employ hDT_RLR.
To illustrate the diagnostic ability of our classifier, we create a graphical plot called receiver operating characteristic (ROC) curve, showing the graphical plot while varying the discrimination threshold. The ROC curve is a plot of the true positive rate (TPR) on the yaxis versus the False Positive Rate (FPR) on the xaxis for every possible classification threshold.
In machine learning, the truepositive rate, also known as sensitivity or recall answers the question how does the classifier predict positive when the actual classification is positive (i.e. not healthy). On the other hand, the falsepositive rate, which is also known as the fallout or probability of false alarm answers the question how does the classifier incorrectly predict positive when the actual classification is negative (i.e. healthy).
The ROC curve is therefore the sensitivity as a function of fallout. Figures 9 and 10 present the comparison of the area under ROC curve (AUC) of our hybrid hDT_RLR classifier for healthy and unhealthy discrimination and for FP type classification (i.e. central or peripheral), respectively, using three different feature extraction methods: (a)Localized Active Contourbased method for key points extraction (LACKPFE));(b)Localized Active Contourbased method for geometric and region features extraction (LACGRFE); and (c)Ensemble of Regression Treebased method for geometric and region features extraction (ERTGRFE).
Our proposed approach ERTGRFE with hybrid classifier hDT_RLR achieves better performance in AUC of 2.7%4.9% higher than other two methods in discriminating healthy from unhealthy (i.e. with paralysis) subjects as shown in Fig. 8. Similarly, for the palsy type classification: central palsy (CP) and peripheral palsy (PP), ERTGRFE plus hDT_RLR hybrid classifier outperformed the two (2) feature extraction methods LACGRFE and LACKPFE used in the experiments with at least 2.5%7.7% as in Fig. 9. Experiments reveal that our method yields more stable results.
Tables 3 and 4 present a comparison of the performance of the three methods for discriminating healthy from unhealthy subjects; and for classifying facial palsy type, respectively. Each approach differs according to the features applied; and the corresponding methods used for extracting such features, which include: (a)Localized Active Contourbased method for key points feature extraction (LACKPFE));(b)Localized Active Contourbased method for geometric and regionbased features (LACGRFE); and (c)Ensemble of Regression Treebased method for geometric and regionbased feature extraction (ERTGRFE).
Table 3 shows that in discriminating healthy from unhealthy subjects, our proposed approach outperforms other methods that use key or salient pointsbased features using LAC model, in terms of sensitivity and specificity with the improvement of 9% and 4.1%, respectively. Similarly, experiments show that our approach outperforms the previous method [4] that used geometric and regionbased features (GRFE) using LAC model in terms of sensitivity and specificity with an improvement of 3.1% AND 5.9%, respectively. On the other hand, Table 4 reveals that for central and peripheral palsy classification, our proposed ERTbased GFRE is way better than the previous approach that solely used key pointsbased features with an improvement of around 12% and 9%, of the sensitivity and specificity performance measure, respectively. Furthermore, experiments reveal that our ERTbased GFRE proposed approach yields better performance, particularly in sensitivity and specificity with an improvement of 6.4% and 6.8%, respectively. Thus, our proposed approach is superior among other three methods.
Discussion
Empirical results [15] reveal that ensemble of regression trees (ERT) model does have an appealing quality that it performs shape invariant feature selection while minimizing the same loss function during training and test time, which significantly lessens time complexity of extracting features. True enough, our experiments reveal that ERTbased method for geometric and region features extraction (ERTGRFE) works well. Furthermore, our proposed approach to combine iris segmentation and the ERTbased key point detection for feature extraction provides a better discrimination of central and peripheral palsy most especially in ‘raising of eyebrows’ and ‘screwing of nose’ movements. It shows changes of the structure on edges of the eye, i.e., the significant difference between the normal side and the palsy side for some facial movements (e.g. eyebrow lifting, nose screwing, and showing of teeth). Also, features based on the combination of iris and key points by utilizing ensemble of regression tree technique can model the typical changes in the eye region. A closer look at the performance of our classifier, as shown in Tables 3 and 4 reveal interesting statistics in terms of the specific abilities of the three methods. Our method proves to have significant contribution in discriminating central from peripheral palsy patients and healthy from facial palsy subjects. The combination of iris segmentation and ERTbased key point approach is more suitable for this operation.
Conclusion
In this study, we present a novel approach to address FP classification problem in facial images. Salient points and iris detection based on ensemble of regression trees are employed to extract the key features. Regularized Logistic Regression (RLR) plus Classification Tree (CT) classifier provides an efficient quantitative assessment of the facial paralysis. Our facial palsy type classifier provides a tool essential to physicians for deciding the medical treatment scheme to begin the patient’s rehabilitation process. Furthermore, our approach has several merits that are very essential to real application:

geometric and iris region features based on ensemble of regression trees and optimized Daugman’s theory (using parabolic function, i.e. 2^{nd} degree polynomial), respectively, allows efficient identification of the asymmetry of the human face, which reveals the significant difference between the normal and afflicted side, which localized active contour model fail to track especially for peculiar images (e.g. wrinkles, excessive mustache, occlusions, etc.);

ERT model does have have a very appealing quality that reduces errors in each iteration, which can be very useful in extracting boundaries of the facial features from the background; and

ERT model does not require proper or perfect identification of initial evolving curves to ensure accurate facial feature extraction.

our method significantly lessens the time complexity of extracting features, without sacrificing the level of accuracy making it more suitable for the real application.
Abbreviations
 1D:

1 Dimensional
 AUC:

area under ROC curve
 CP:

Central palsy
 CT:

Classification tree
 DT:

Decision Tree
 ERT:

Ensemble of regression trees
 ERTGRFE:

Ensemble of Regression Treebased method for geometric and region features extraction
 FP:

Facial Paralysis
 FPR:

False Positive Rate
 HOG:

Histogram of Oriented Gradients
 hDT_RLR:

hybrid Decision tree and Regularized Logistic Regression
 IC:

Inner Canthus
 IO:

Infra Orbital
 IRB:

Institution Review Board
 LAC:

Localized Active Contour
 LACGRFE:

Localized Active Contourbased method for geometric and region features extraction
 LACKPFE:

Localized Active Contourbased method for key points extraction
 MA:

Mouth Angle
 ML:

Machine learning
 NB:

Naive bayes
 NT:

Nose tip
 OC:

Outer Canthus
 PP:

Peripheral palsy
 RF:

Random Forest
 RLR:

Regularized Logistic Regression
 ROC:

Receiver operating characteristic curve
 ROI:

Region of interest
 SO:

Supraorbital
 SVM:

Support Vector Machine
 TPR:

True positive rate
 UE:

Upper Eyelids
References
 1
Baugh R, Ishii G, RSchwartz S, Drumheller CM, Burkholder R, Deckard N, Dawson C, Driscoll C, BoydGillespie M. Clinical practice guideline bellspalsy. OtolaryngolHead Neck Surg. 2013; 149:1–27.
 2
Peitersen E. Bell’s palsy: the spontaneous course of 2,500 peripheral facialnerve palsies of different etiologies. Acta OtoLaryngol. 2002; 122(7):4–30.
 3
Kanerva M. Peripheral facial palsy grading, etiology, and melkersson–rosenthal syndrome. PhD thesis. Finland: University of Helsinki; 2008.
 4
Barbosa J, Lee K, Lee S, Lodhi B, Cho J, Seo W, Kang J. Efficient quantitative assessment of facial paralysis using iris segmentation and active contourbased key points detection with hybrid classifier. BMC Med Imaging. 2016; 16:23–40.
 5
May M. Anatomy for the clinician In: May M, Schaitkin B, editors. The Facial Nerve, May’s Second Edition. 2nd. New York: Thieme, New York: 2000. p. 19–56.
 6
Wachtman GS, Liu Y, Zhao T, Cohn J, Schmidt K, Henkelmann TC, VanSwearingen JM, Manders EK. Measurement of asymmetry in persons with facial paralysis. In: Combined Annu Conf. Robert H. Ivy. Ohio Valley: Societies of Plastic and Reconstructive Surgeons: 2002.
 7
Liu Y, Schmidt K, Cohn J, Mitra S. Facial asymmetry quantification for expression invariant human identification. Comput Vis Image Underst J. 2003; 91:138–59.
 8
He S, Soraghan J, O’Reilly B, Xing D. Quantitative analysis of facial paralysis using local binary patterns in biomedical videos. IEEE Trans Biomed Eng. 2009; 56:1864–70.
 9
Wang S, Li H, Qi F, Zhao Y. Objective facial paralysis grading based on pface and eigenflow. Med Biol Eng Comput. 2004; 42:598–603.
 10
Anguraj K, Padma S. Analysis of facial paralysis disease using image processing technique. Int J Comput Appl(0975 – 8887). 2012; 54:1–4.
 11
Dong J, Ma L, Li W, Wang S, Liu L, Lin Y, Jian M. An approach for quantitative evaluation of the degree of facial paralysis based on salient point detection. In: International Symposium on Intelligent Information Technology Application Workshops. New York City: IEEE: 2008.
 12
Liu L, Cheng G, Dong J, Qu H. Evaluation of facial paralysis degree based on regions In: IEEE Computer Science Society Washington U.DC (ed.), editor. Third International Conference onKnowledge Discovery and Data Mining. Washington, DC: IEEE Computer Science Society: 2010. p. 514–7.
 13
Liu L, Cheng G, Dong J, Qu H. Evaluation of facial paralysis degree based on regions. In: Proceedings of the 2010 Third International Conference onKnowledge Discovery and Data Mining. Washington, DC: IEEE Computer ScienceSociety: 2010. p. 514–7.
 14
Otsu N. A threshold selection method from graylevel histogram. IEEE Trans Syst, Man Cybern. 1979; 9:62–6.
 15
Kazemi V, Sullivan J. One millisecond face alignment with an ensemble of regression trees. In: Proc. IEEE Conf Comput Vis Pattern Recog. New York City: IEEE: 2014. p. 1867–1874.
 16
Samsudin W, Sundaraj K. Image processing on facial paralysis for facial rehabilitation system: A review. In: IEEE International Conference on Control System, Computing and Engineering. Malaysia: IEEE: 2012. p. 259–63.
 17
Fasel B, Luttin J. Automatic facial expression analysis: Survey. Pattern Recog. 2003; 36:259–75.
 18
Ngo T, Seo M, Chen Y, Matsushiro N. Quantitative assessment of facial paralysis using local binary patterns and gabor filters. In: Proceedings of the 5th International Symposium on Information and Communication Technology (SoICT). New York: ACM: 2014. p. 155–61.
 19
Déniz O, Bueno J, Salido F, De la Torre F. Face recognition using histograms of oriented gradients. Pattern Recog Lett. 2011; 32:1598–603.
 20
Daugman J. How iris recognition works. New York City: IEEE: 2004. p. 21–30. vol. 14, IEEE Trans Circ Syst Vi Technol.
 21
Dalal N, Triggs B. Histograms of oriented gradients for human detection. In: Proc IEEE Conf Computer Vision and Pattern Recognition. New York City: IEEE: 2005.
 22
Lyons M, Budynek J, Akamatsu S. Automatic classification of single facial images. IEEE Trans Pattern Anal Mach Intell. 1999; 21(12):57–62.
 23
Lyons M, Akamatsu S, Kamachi M, Gyoba J. Coding facial expressions with gabor wavelets. In: Third IEEE Conf. Face and Gesture Recognition. IEEE: 1998. p. 200–5.
 24
Proenca H, Alexandre L. Ubiris: A noisy iris image database. In: Proceed. of ICIAP 2005  Intern. Confer. on Image Analysis and Processing. United States: Springer Link: 2005. p. 970–977.
Acknowledgments
There are no acknowledgments.
Funding
This research was supported by the National Research Foundation of Korea (NRF2017M3C4A7065887) and National IT Industry Promotion Agency grant funded by the Ministry of Science and ICT and Ministry of Health and Welfare (NO. C1202181001, Development Project of The Precision Medicine Hospital Information System (PHIS)); and the scholarship was granted by the Korean Government Scholarship Program  NIIED, Ministry of Education, South Korea.
Availability of data and materials
The datasets generated and/or analysed during the current study are not available due to patient privacy.
Author information
Affiliations
Contributions
Conceived and designed the methodology and and analyzed the methods: JB, JK. Performed experiments and programming: JB. Collected images and performed manual evaluation: WS. All authors wrote the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study was approved by the Institution Review Board (IRB)of Korea University, Guro Hospital (with reference number MD14041002). The board permitted not taking an informed consent due to the retrospective design of this study.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Barbosa, J., Seo, W. & Kang, J. paraFaceTest: an ensemble of regression treebased facial features extraction for efficient facial paralysis classification. BMC Med Imaging 19, 30 (2019). https://doi.org/10.1186/s1288001903308
Received:
Accepted:
Published:
Keywords
 Facial paralysis classification
 Facial paralysis objective evaluation
 Ensemble of regression trees
 Salient point detection
 Iris detection
 Facial paralysis evaluation framework