Novel geometric features for off-line writer identification

Writer identification is an important field in forensic document examination. Typically, a writer identification system consists of two main steps: feature extraction and matching and the performance depends significantly on the feature extraction step. In this paper, we propose a set of novel geometrical features that are able to characterize different writers. These features include direction, curvature, and tortuosity. We also propose an improvement of the edge-based directional and chain code-based features. The proposed methods are applicable to Arabic and English handwriting. We have also studied several methods for computing the distance between feature vectors when comparing two writers. Evaluation of the methods is performed using both the IAM handwriting database and the QUWI database for each individual feature reaching Top1 identification rates of 82 and 87 % in those two datasets, respectively. The accuracies achieved by Kernel Discriminant Analysis (KDA) are significantly higher than those observed before feature-level writer identification was implemented. The results demonstrate the effectiveness of the improved versions of both chain-code features and edge-based directional features.


Introduction
Automatic writer identification is important in forensic document examination.Numerous cases have dealt with evidence provided by handwritten documents, such as wills and ransom notes [1].Moreover, writer identification can be used in handwriting recognition when adapting the recognizers to a specific type of writers [2] and in handwriting synthesis when generating a text as it would have been written by a specific writer [3].Writer identification methods generally consist of two main steps.The first step is feature extraction, in which discriminating features are extracted from the handwritten documents to be compared while the second step involves matching or classification in which a comparison between the features is computed, and a decision regarding the authorship is made according to the distance between the extracted features.What makes a system powerful and robust is closely related to a strong feature extraction step because the extraction of discriminative features helps to distinguish between writers.
Automatic methods for writer identification can be classified into two main categories: codebook-based and feature-based approaches.In codebook-based approaches, the writer is assumed to act as a stochastic generator of graphemes.The probability distribution of a grapheme is a characteristic of each writer and can efficiently be used to distinguish between different writers.The methods in this category mainly differ in how the handwriting is segmented into graphemes and how the graphemes are clustered.On the other hand, feature-based approaches compare the handwriting samples according to geometrical [4], structural [5], or textural features [6,7].Feature-based approaches are proven to be efficient and are generally preferred when only a limited amount of handwriting data is available.
In the remainder of this section, an overview of the works done in the field of off-line writer identification is presented.
Srihari et al. proposed a set of macro-features that are extracted from a document, paragraph, or word level (i.e., entropy of gray values; gray-level threshold; number of black pixels; number of interior and exterior contours; number of vertical, horizontal, negative, and positive slope components; and slant) and micro-features that are extracted at the word or character level (i.e., gradient, structural, and concavity features).This approach has been validated on two datasets of 711 writers using the same letter [1] achieving a classification performance of 89 and 87 %, respectively.Said et al. described a text-independent writer identification method based on Gabor filtering and grayscale cooccurrence matrices.The authors obtained 95 and 88 % recognition rate, but evaluated on a small set of 10 and 15 writers [7].
Marti et al. used text line-based features for text-independent writer identification [5].The authors used features related to the position of the top line, the bottom line, the upper baseline, and the lower baseline; the width; and the slant of the writing.They also introduced a set of features based on fractal geometry to distinguish between badly formed and legible writing.The recognition rate obtained with their method reached 90 % on a subset of 50 writers from the IAM database [8].The IAM database will also be used in the present study though more writers will be considered.
A further study by Hertel and Bunke [9] resulted in novel features based on connected components, enclosed regions, lower and upper contours, fractal features, and basic features (including writing skew and slant, the height of the three main writing zones, and the width of the writing).This study has been validated on a subset of 50 writers of the IAM database.A recognition rate of 99 % was obtained for this small data set.Connected component feature alone obtained 31 % in the same dataset.
Schlapbach and Bunke used the results of the HMM text recognizer for both the identification and verification of writers [10].The authors proposed recognizers that are trained for each individual writer.Subsequently, a handwritten sample in question is provided as input to all the recognizers, which the corresponding output likelihood ratio values.By sorting these values, the identity of the most probable writer then is determined.This study has validated their approach using 100 writers in the IAM database with a recognition rate of 96 %.In a follow-up study, the authors proposed an improvement by 2 % to the system by deploying Gaussian mixture models [11].
Bensefia et al. used a textual-based information retrieval model for the writer identification stage [12].This makes it possible to use a particular feature space based on feature frequencies of occurrence.Image queries are handwritten documents projected in this feature space.This approach achieved an 86 % identification rate on a subset of 150 documents of the IAM database.
Siddiqi and Vincent proposed the use of the redundancy of graphemes to characterize writer individuality [13].The authors also proposed some chain code-based features extracted from the handwriting contours.Their best performing feature (local stroke direction histogram) achieved a 77 % identification rate with the IAM database [14].The two categories of features were combined in a follow-up study [15] and validated with the IAM database.Their best performing feature is a distribution of chain code with achievement of 79 %.The chain-code feature initially introduced by these authors will be further improved in this study.
Bulacu and Schomaker used edge-based directional probability distribution functions (PDFs) as features for text-independent writer identification.The joint PDF of ''hinged'' edge-angle combinations outperformed all the other evaluated features (including Contour-direction PDFs, Direction co-occurrence PDFs, Grapheme emission PDFs, Run length on background PDFs, and Autocorrelation during horizontal raster scanning) [16].By using Contour-hinge PDF features a recognition rate of 81 % was achieved.
A study by Imdad et al. used steered Hermite features combined with the SVM classifier [17].The proposed method achieved an 83 % identification rate on a subset of 30 writers from the IAM database.
Muda et al. showed that the discretization of features can significantly improve the identification rates [18].Discretization is performed by exploring the partitioning of features into intervals and unifying the values for each interval.The method achieved around 99 % but 60 writers only were used in the evaluation.
Bozekova ´combined grapheme features obtained using a Kohonen Self-Organizing Map and specified structural features to achieve a 96.5 % identification rate on a subset of 40 writers from the IAM database [19].Dolega et al. used derived canonical stroke frequency descriptors from handwritten text to identify writers [20].The authors reached an 88 % identification rate on a subset of 50 writers in the IAM database.
Steinke et al. combined local features that use different mathematical procedures, such as the reproduction of the write line of individual characters by Legendre polynomials, and global textural features [21].The proposed method achieved a 99.5 % identification rate on a subset of 93 writers.Jain and Doermann extracted codebooks of K-adjacent segments from handwriting text to characterize writer individuality [22].The method achieved a 93 % identification rate on 300 writers from the IAM dataset.
From the review and analysis of existing works available in the literature it can be concluded that existing methods are usually validated using a small number of writers which are usually less than 100, whereas the IAM dataset contains more than 650 writers.In addition, one can note that there has been a consistent amount of research on the combination of features; however the performances of individual features have not been studied and need further analysis with a view to determine their discriminative power to maximize their recognition performances and also to aid as to how one can combine them efficiently.Another motivation of the work relates to the development of novel features especially to capture the peculiar characteristics of the writers' handwriting relating to direction, curvature and also tortuosity.To address the above issues, this paper has the following contributions: • New features based on direction (f1), curvature (f2), tortuosity (f3) quadruple-order chain code (f7) are proposed including an implementation and evaluation on much recent large databases.To the best of the authors' knowledge, the proposed technique is the first attempt the field of handwriting identification.• The paper also proposes an improvement of state-ofthe-art edge-based directional features (by using a filled moving window instead of an edge moving window) and the chain code-based features (by using a fourthorder chain-code list) in terms of the discriminative power (f18-f26).We demonstrate that these improvements lead to much enhanced identification rates.• A kernel discriminant analysis has been used in order to combine the features using several metrics including X 2 , L 2 , and L 1 distances to evaluate the performances.An analytical study is then given showing improvement in the performances.
The remainder of the paper is organized as follows: Sect. 2 introduces the proposed methodology and the system's components including the description of the features used and the matching strategy used.In this section a set of novel features are also discussed.Section 3 discusses the datasets used in our experimentation and evaluation.Section 4 presents the results and their analysis while Sect. 5 concludes this paper.

Proposed methodology
Similar to any writer identification system, our proposed method consists of two steps: feature extraction and matching steps.In the following we will discuss both steps in detail.

Feature extraction
In this stage, the characterizing features are extracted from the handwriting.To develop a pen-independent system, images are first binarized using the Otsu thresholding algorithm [25].It is worth noting that in writer identification, features do not correspond to a single value but to a probability distribution function (PDF) extracted from the handwriting images to characterize a writer's uniqueness [26].Novel sets of features were extracted from different handwritten datasets.In this paper we propose to extract features based on direction (which we refer to as f1), curvature (f2), tortuosity (f3), quadruple-order chain code (f7), and the edge-based directional features using the whole window computed from size 2 (f18, whose PDF size is 12) to size 10 (f26, whose PDF size is 220).We have also considered and used other state-of-the-art features, such as chain-code features (f4-f6) and Edge-based directional features (f8-f17) to compare the results.These features will be explained in the what follows.

Direction (f1)
Direction is known to be useful for characterizing writers [16].The methodology used in this paper is novel and has been used before in writer identification.Its implementation is somehow similar to the method proposed by Matas et al. [27].First, we compute the Zhang skeleton of the binarized image.This skeleton is well known for not producing parasitic branches in contrast to most skeletonization algorithms [28].The skeleton is then segmented at its junction pixels.Then, we move along the pixels of the obtained segments of the skeleton using the predefined order, favoring the four-connectivity neighbors, as shown in Fig. 1a.A result of such an ordering is shown in Fig. 1b.For each pixel p, we consider the 2N ? 1 neighboring pixels centered at position p.The linear regression of these pixels gives a good estimation of the tangent at pixel p (Fig. 1c).The value of N has empirically been set to 5 pixels.This feature is illustrated in Fig. 1d.
The PDF of the resulting directions is computed as a vector of probabilities for which the size has been empirically set to 10.Note that this way of computing directions has never been used before in off-line writer identification.
Here, we have introduced a novel method for computing curvature and its deployment in the field of writer identification for the first time.It is to ne noted that this technique has been used previously for estimating the curvature of the peaks and valleys in optical soundtracks [30].We have modified and adapted this method to handwritings as follows: for each pixel p belonging to the contour, we consider a neighboring window of size t.Inside this window, we compute the number of pixels n 1 and n 2 which belong to the background and the foreground, respectively (Fig. 2a).Therefore, the difference n 1n 2 increases with the local curvature of the contour.Therefore, we estimate the curvature as follows: This value is illustrated in Fig. 2b on a binary shape for which t has been empirically set to 5.
The PDF of curvature is computed in a vector whose size has been empirically set to 100.To the best of our knowledge, this method of computing curvature is also novel in the field of off-line writer identification.

Tortuosity (f3)
In this work we propose to use tortuosity as a novel feature to allow us to distinguish between fast writers who produce smooth handwriting and slow writers who produce ''tortuous''/twisted handwriting.To estimate tortuosity, for each pixel p of the text, we determine the longest line segment that traverses p and is completely included inside the foreground (Fig. 3a).An example of estimated tortuosity is shown in Fig. 3b.
The PDFs of the angles of the longest traversing segments are stored in a vector whose size has been set to 10.

Chain-code features (f4-f7)
Chain codes are generated by browsing the contour of the text and assigning a number to each pixel according to its location with respect to the previous pixel.Figure 4 shows a contour and its corresponding chain code.
These features make it possible to characterize the detailed distribution of curvature in the handwritings.Chain codes can be applied at different levels: f4: PDF of i patterns in the chain-code list such that i[0, 1,…,7.This PDF has a size of 8. f5: PDF of (i, j) patterns in the chain-code list such that i, j[0, 1,…,7.This PDF has a size of 64.
Similarly, f6 and f7 correspond to PDFs of (i, j, k) and (i, j, k, l) in the chain-code list, and their respective sizes are 512 and 4,096, respectively.
Although, the chain-code features f4, f5, and f6 have previously been applied to writer identification [14], we propose to deploy the quadruple-order chain code f7 for the first time in writer identification.Initially introduced [16], these features provide detailed distributions for the direction and can also be applied at several scales by positioning a window centered at each contour pixel and counting the occurrences of each direction, as shown in Fig. 5a.This feature has been computed from size 1 (f8, whose PDF size is 4) to size 10 (f17, whose PDF size is 40).We have also extended these features to include not only the contour of the moving window but also of the whole window (Fig. 5b).This feature has been computed from size 2 (f18, whose PDF size is 12) to size 10 (f26, whose PDF size is 220).

Matching of feature vectors
When comparing a query document q against any given document i, the difference between their features is computed as shown in Fig. 6.
In this study, three different distances were considered and evaluated: • The v 2 distance is the first method used for comparing two PDFs: • The L 2 distance or Euclidian distance was also tested: • The L 1 distance was also tested: In addition to computing one distance between each pair of feature vectors, we have also computed the list of distances at the element level for each pair of vectors: These distances are combined using a kernel discriminant analysis with spectral regression (SR-KDA).A description of this classifier is given below: Let x i [ Rd, i = 1,… and m be training vectors represented as an m 9 m kernel Matrix K is defined such that K x i ; , where U(x i ) and U(x j ) are the embedded data items x i and x j , respectively.If m denotes a projective function into the kernel feature space, then the objective function for KDA is where C b and C t denote the between-class and total scatter in feature space, respectively.Equation 1 can be solved by the eigenvalue problem C b ¼ kC t .It is proven that Eq. ( 1) is equivalent to the following: where a = [a 1 , a 2 ,…..a m ] T is the eigenvector satisfying KWKa = kKKa.W = (W l ) l=1,…,n is a (m 9 m) block diagonal matrix of labels arranged such that the upper block corresponds to positive examples and the lower one to negative examples of the class.Each eigenvector a gives a projection function m into the feature space.
Instead of solving the eigenvalue problem in KDA, the KDA projections can be obtained by the following two linear equations: where is an eigenvector of W, I is the identity matrix, and d [ 0 is a regularization parameter.W = (W l ) l=1,…,n is an (m 9 m) block diagonal matrix of labels arranged such that the upper block corresponds to positive examples and the lower one to negative examples of the class.Eigenvectors / are obtained directly from the Gram-Schmidt method.Because (K ?dI) is positive definite, the Cholesky decomposition is used to solve the linear equations in (3).Thus, for the resolution of linear system (3), the system becomes i.e., the system is solved first to find vector h and then vector a.In summary, SRKDA only needs to solve a set of regularized regression problems, and there is no eigenvector computation involved.This approach results in substantial improvement in terms of computational cost and makes it possible to handle large kernel matrices.After obtaining a, the decision function for new data items is where Kðx; x iÞ ¼ hUðxÞ; Uðx iÞi.The proposed features have different discriminative levels, which are described below.

Training and testing datasets
For the experimental set up of our system, we have chosen two different datasets: IAM [8] and QUWI [23].The IAM dataset contains English handwritings whereas the QUWI dataset contains both English and Arabic handwritings.The IAM dataset [8] is the most widely used dataset in writer identification.It contains handwritings of 657 writers, each of whom produced two samples of text.The IAM database consists of forms with handwritten English text of variable content that has been scanned at 300 dpi with 8 bits/pixel in grayscale.In addition to writer identity, the images are accompanied by extensive segmentation and ground-truth information at the text line, sentence, and word levels.This dataset includes a variable number of handwritten pages per writer, from one page (350 writers) to 59 pages (one writer).
To have comparable experimental conditions with existing state-of-the-art methods, we have modified the IAM dataset to contain two samples per writer.We kept only the first two documents for the writers who produced more than two documents, and we split the document roughly in half for the writers with a unique page in the original set.This modified IAM is used for testing purposes in a manner similar to that used previously [14][15][16].This modified dataset contains lowercase handwritings from 657 people, with two samples per writer.The amount of ink is roughly equal in the two samples belonging to one writer but varies among writers from three lines up to a full page.Note that the third and fourth documents of the writers who produced at least four documents are used for training purposes.
We also used the QUWI dataset [23], which is a dataset built at Qatar University consisting of handwritings from 1,017 writers.This dataset has been scanned with a spatial resolution of 600 dpi.The writers were asked to produce four pages of text: one similar text and one different text in English and one similar text and one different text in Arabic.Part of the QUWI dataset has been used for organizing a competition for writer identification [24].

Evaluation
The most widely used database for writer identification is the IAM database [8] which was described above.Similar to the method described in [16], we have used only the first two documents from the writers who produced more than two documents, and for the writers who produced only one document, we split the document into two separate documents.Note that the comparison with other systems is to be considered approximate and not exact because the current IAM database contains 657 writers, not 650.
Therefore, each document is compared against all 1,313 other documents with only one possible correct match.If the distance between the document of interest and the correct match is the smallest among all possibilities, then the document is said to be correctly identified.The TOP 10 identification rate considers the matching as correct if the corresponding distance is among the 10 smallest distances.
As previously stated in Sect.3, the other dataset used in this study is QUWI dataset.The pages of this been segmented into three paragraphs; two of them are used for training while the third one is used for testing.The identification rates obtained for the presented features for both datasets where tested using the v 2 distance, L 2 distance, L 1 distance, and KDA classifier.The results obtained for each feature using the three types of distances are given in Table 1.
The results of the KDA classifier on the test datasets are presented in Table 2.
As mentioned previously, in this paper, the proposed features were computed as follows: directional feature is determined using a linear regression technique (f1) while the curvature (f2), tortuosity (f3), quadruple-order chain code (f7) and edge-based directional features are computed using the whole window from a size 2 (f18, whose PDF size is 12) to size 10 (f26, whose PDF size is 220).Other chain-code features (f4 to f6) and direction features (f8 to f17) were reproduced for comparison purposes.
The KDA classifier achieves the best TOP1 recognition rate for most features except the following: f2, f16, and f17 in all datasets; f8 in all QUWI datasets; and f7 and f9 in the QUWI Arabic same-text dataset.The TOP10 recognition rates using kDa classifier are also generally the highest in the IAM dataset.
In the case of the QUWI dataset, the highest recognition rates were obtained using the L 1 distance.It was also found that the first 17 features give high performances for most sub-datasets when using the other distances.The KDA classifier produced the highest recognition rate for the 10 highest-ranked writers for the QUWI dataset for features f18 to f26.
For the IAM dataset, the quadruple-order chain-code feature f7 yielded the highest recognition rate of 82 % for the top writer and of 92 % for the top ten writers using the  KDA classifier.This approach marginally outperformed the approaches described in [15,16] for other individual features using dataset.The proposed new edgebased directional features (f16 to f18) resulted in high recognition rates as well.Other individual features produced reasonable recognition rates for the IAM dataset.
In the case of the QUWI dataset, the highest recognition rate was obtained when using the new edge directional feature for English and Arabic for both text-dependent and text-independent approaches.Recognition rates of 70 % for text-dependent Arabic handwriting for the top writer and 86 % for the top 10 ranked writers were obtained using a edge-based directional features of sizes 7 and 8 when KDA is utilized.In addition, window edge-based directional features for windows of size 8 and 9 have resulted in highest recognition rates for text-dependent Arabic text achieving 87.6 % for the top-ranked writer and 96.25 % for the topranked writers.For English handwriting, recognition rates of 75 % for the top-ranked writer and 90.75 % for the top-ranked writers using a window of size 5 were obtained when edge-based directional features are used for different text samples.In this case, a window of size 6 for edge-based directional feature produced the highest recognition rate for  text-dependent handwritten text.For example, a recognition rate of 83.01 and 93.42 % was achieved for the top-ranked writers and 10 top-ranked writers, respectively.Figure 7 shows the writer recognition rate for selected individual features (f7, f22-f26).These features resulted in the highest recognition rate with other features.Features such as direction (f1), curvature (f2), and tortuosity (f3) can be used with other features to improve the recognition rate.Furthermore, features such as tortuosity (f3) can be used to measure the writing speed.
Table 3 shows the performance of the new features compared to the state-of-the-art counterparts when using the IAM dataset with approximately 650 writers.As shown in the table, the new features result in the highest recognition rate compared with the state-of-the-art features.The new fourth-order chain-code features show a significant improvement with regard to other features especially when combined with the KDA classifier.

Conclusion
This paper has presented several new features for writer identification.These include curvature, direction, and tortuosity features.We have shown through an evaluation and analysis of the recognition performance results obtained using two well-known datasets the usefulness of these features for writer identification.In particular, a comparison of the recognition performances of these features when deployed individually against state-of-the art features was also carried out.The findings have led us to propose methods to improve the discriminative power of both chain-code and edge-based directional features.To further ascertain the discriminative power of the proposed features we have utilized several distance metrics including X 2 , L 2 , and L 1 distances and the distance at the feature level when combined with KDA classifier by computing the differences between the feature vectors.The latter approach outperforms all the other distance measures previously reported in the literature.Work on the use of these features is ongoing for the prediction of demographic categories, including age, nationality, handedness, and gender.Future work will include the validation of the method with other languages and offline signature verification.

Fig. 1 Fig. 2 aFig. 3 Fig. 4 a
Fig. 1 Computing local direction.a The predefined order for traversing shapes.b Example of an ordered shape.c Estimating direction by computing the linear regression of neighboring pixels for three different bold-colored regions.d Binary image and its corresponding Zhang skeleton; the red color corresponds to a p/2 tangent, and the blue color corresponds to a zero tangent (color figure online)

Fig. 5 Fig. 6
Fig. 5 Counting the edge-based directional features when considering a the contour of the moving window and b the whole moving window

Fig. 7
Fig. 7 Writer identification rate for the ten top-ranked writers using different datasets in Arabic and English

Table 1
Top1 identification rates using each category of distance

Table 2
Identification rates of the presented features using the KDA classifier

Table 3
State-of-the-art writer identification performance of individual features compared to the new features on the IAM dataset with 650 writers