Next Article in Journal
Exploiting Multiple Detections for Person Re-Identification
Next Article in Special Issue
Open Datasets and Tools for Arabic Text Detection and Recognition in News Video Frames
Previous Article in Journal
Studies of Ancient Russian Cultural Objects Using the Neutron Tomography Method
Previous Article in Special Issue
Transcription of Spanish Historical Handwritten Documents with Deep Neural Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Binarization Algorithm for Historical Documents

1
Departamento de Eletrônica e Sistemas, Centro de Tecnologia, Universidade Federal de Pernambuco, Recife-PE 50670-901, Brazil
2
Centro de Informática, Universidade Federal de Pernambuco, Recife-PE 50740-560, Brazil
3
Departamento de Estatística e Informática, Universidade Federal Rural de Pernambuco, Recife-PE 52171-900, Brazil
4
Programa de Pós-Graduação em Engenharia Elétrica, Universidade Federal de Pernambuco, Recife-PE 50670-901, Brazil
*
Author to whom correspondence should be addressed.
J. Imaging 2018, 4(2), 27; https://doi.org/10.3390/jimaging4020027
Submission received: 31 October 2017 / Revised: 16 January 2018 / Accepted: 16 January 2018 / Published: 23 January 2018
(This article belongs to the Special Issue Document Image Processing)

Abstract

:
Monochromatic documents claim for much less computer bandwidth for network transmission and storage space than their color or even grayscale equivalent. The binarization of historical documents is far more complex than recent ones as paper aging, color, texture, translucidity, stains, back-to-front interference, kind and color of ink used in handwriting, printing process, digitalization process, etc. are some of the factors that affect binarization. This article presents a new binarization algorithm for historical documents. The new global filter proposed is performed in four steps: filtering the image using a bilateral filter, splitting image into the RGB components, decision-making for each RGB channel based on an adaptive binarization method inspired by Otsu’s method with a choice of the threshold level, and classification of the binarized images to decide which of the RGB components best preserved the document information in the foreground. The quantitative and qualitative assessment made with 23 binarization algorithms in three sets of “real world” documents showed very good results.

1. Introduction

Document image binarization plays an important role in the document image analysis, compression, transcription, and recognition pipeline [1]. Binary documents claim for far less storage space and computer bandwidth for network transmission than color or grayscale documents. Historical documents drastically increase the degree of difficulty for binarization algorithms. Physical noises [2] such as stains and paper aging affect the performance of binarization algorithms. Besides that, historical documents were often typed, printed or written on both sides of sheets of paper and the opacity of the paper is often such as to allow the back printing or writing to be visualized on the front side. This kind of “noise”, first called back-to-front interference [3], was later known as bleeding or show-through [4]. Figure 1 presents three examples of documents with such a noise extracted from the three different datasets used in this paper in the assessment of the proposed algorithm. If the document is exhibited either in true-color or gray-scale, the human brain is able to filter out that sort of noise keeping its readability. The strength of the interference present varies with the opacity of the paper, its permeability, the kind and degree of fluidity of the ink used, its storage, age, etc. Thus, the difficulty for obtaining a good binarization performance capable of filtering-out such a noise increases enormously, as a new set of hues of paper and printing colors appear. The direct application of binarization algorithms may yield a completely unreadable document, as the interfering ink of the backside of the paper overlaps with the binary one in the foreground. Several document image compression schemes for color images are based on “adding color” to a binary image. Such compression strategy is unable to handle documents with back-to-front interference [5]. Optical Character Recognizers (OCRs) are also unable to work properly for such documents. Several algorithms were developed specifically to binarize documents with back-to-front interference [3,4,6,7,8,9]. There is no binarization technique to be an all case winner as many parameters may interfere in the quality of the resulting image [9]. The development of new binarization algorithms is still an important research topic. International competitions on binarization algorithms, such as DIBCO - Document Image Binarization Competition [10], are an evidence of the relevance of this area.
This paper presents a new global filter [1] to binarize documents, which is able to remove the back-to-front noise in a wide range of documents. Quantitative and qualitative assessments made in a wide variety of documents from three different “real-world” datasets (typed, printed and handwritten, using different kinds of paper, ink, etc.) allow to witness the efficiency of the proposed scheme.

2. The New Algorithm

The algorithm proposed here is performed in four steps: 1. decision-making for finding the vector of parameters of the image to be filtered, 2. filtering the image using a bilateral filter, 3. splitting the image into the RGB components, and performing their binarization using a method inspired by Otsu’s algorithm for each RGB channel, and 4. choice of which of the RGB components best preserved the document information in the foreground, which is considered the final output of the algorithm. Figure 2 presents the block diagram of the proposed algorithm. The functionality of each block is detailed as follows.

2.1. The Decision Making Block

The decision making block takes as input the image to be binarized and outputs a vector with four parameters: the value of the kernel (kernel) for the bilateral filter and three threshold values (tR, tG, tB) that will be later used in the modified Otsu filtering.
The training of the binarization process proposed here is made with synthetic images which were generated as explained in Section 2.2. After filtering, the matrix of co-occurrence probabilities between the original image and of the binary image was calculated for each of the images in the document training set, whose generation is explained below.
The probabilistic structure applied in the analysis to each of the images in the training set is similar to the transmission of binary data in a Binary Asymmetric Channel, as shown in Figure 3. The probabilities P(f/b) and P(b/f) represent an additive noise in communication channels in information theory, here it represents the inability of the algorithm to correct the back-to-front interference of the image tested in the binarization process. The probabilities P(b/b) and P(f/f) are calculated from the pixel-to-pixel comparison of the binarized image generated by the proposed algorithm with the ground-truth image.
The background-background probability is a function that needs to be optimized in the decision-making block, mapping background pixels (paper) from the original image onto white pixels of the binary image. It depends of all the parameters of the original image texture, strength of the back to front interference (simulated by the coefficient α), paper translucidity, etc. for each RGB channel. Thus, one can represent this dependence as:
P(b/b) = f(α,R,G,B).
The optimal threshold tc* for each channel is calculated in the decision-making block, the index c can be R, G or B, maximizing P(b/b):
tc* = MaxP(b/b),
subject to a given criterion P(f/f) ≥ M. The criterion used here was M = 97%, that is at most 3% of the foreground pixels may be incorrectly mapped. During the training phase, the best tc* will be chosen from the three channels, which best maximizes the P(b/b) for each of the images in the training set. The matrix of co-occurrence probability is calculated and the decision maker chooses the best binary image. The decision-making block was trained with 32,000 synthetic images in such a way to, given a real image to be binarized, it finds the optimal threshold parameters.

2.2. Generating Synthetic Images

The Decision-Making Block needs training to “learn” about the optimal threshold parameters and the value of the kernel to be used in the bilateral filter. Such training must be done using controlled images which are synthesized to mimic the different degrees of back-to-front interference, paper aging, paper translucidity, etc. Figure 4 presents the block diagram for the generation of synthetic images. Two binary images of documents of different nature (typed, handwritten with different pens, printed, etc.) are taken: F—front and V—verso (back). The front image is blurred with a weak Gaussian filter to simulate the digitalization noise [1], the hues that appear in after document scanning.
The verso image is “blurred” by passing through two different Gaussian filters that simulate the low-pass effect of the translucidity of the verso as seen in the front part of the paper. Two different parameters were used to simulate two different classes of paper translucidity. The “blurred” verso image is now faded with a coefficient α varying between 0 and 1 in steps of 0.01. Then, a circular shift of the lines of the document is made of either 5 or 10 pixels, to minimize the chances of the front and verso lines coincide entirely. Finally, the two images are overlapped by performing a “darker” operation pixel-by-pixel in the images. Paper texture is added to the image to simulate the effect of document aging. The texture pattern was extracted from document from late 19th century to the year 2000. The analysis of 3450 documents representative of a wide variety of documents of such a period was analyzed yielding 100 different clusters of textures. The synthetic texture to be applied to the image to simulate paper aging is generated using those 100 clusters by image quilting [11] and randomly, as explained in reference [9]. The training performed in the current version of the presented algorithm was made with 16 of those 200 synthetic textures. The total number of images used for training here was thus 16 (textures), times 10 (0 < α < 1 in steps of 0.10), times 2 blur parameters for the Gaussian filters, times 100 different binary images, totaling 32,000 images. Details of the full generation process of the synthetic image database are out of the scope of this paper and may be found in reference [9].

2.3. The Bilateral Filter

The bilateral filter was first introduced by Aurich and Weule [12] under the name “nonlinear Gaussian filter”. It was later rediscovered by Tomasi and Manduchi [13] who called it the “bilateral filter” which is now the most commonly used name according to reference [14].
The bilateral filter is a technique to smoothen images while preserving their edges. The filter output at each pixel is a weighted average of its neighbors. The weight assigned to each neighbor decreases with both the distance values among pixels of the image plane (the spatial domain S) and the distance on the intensity axis (the range domain R). The filter applies spatial weighted averaging without smoothing the edges. It combines two Gaussian filters; one filter works in the spatial domain, while the other filter works in the intensity domain. Therefore, not only the spatial distance but also the intensity distance is important for the determination of weights. The bilateral filter combines two stages of filtering. These are the geometric closeness (i.e., filter domain) and the photometric similarity (i.e., filter range) among the pixels in a window of size N × N. Let I(x,y) be a 2D discrete image of size N × N, such that {x,y} {0, 1, ..., N − 1} X {0, 1, ..., N − 1}. Assume that I(x,y) is corrupted by an additive white Gaussian noise of variance σ n 2 . For a pixel (x,y), the output of a bilateral filter can be as described by Equation (1):
I B F ( x , y ) = 1 K i = x d x + d j = y d x + d G s ( i ; x , j ; y ) G r [ I ( i , j ) , I ( x , y ) ] I ( i , j ) ,
where I(x,y) is the pixel intensity in the image before applying the bilateral filter, IBF(x,y) is the resulting pixel intensity after applying the bilateral filter and d is a non-negative integer such that (2d + 1) × (2d + 1) stands for the size of the neighborhood window. Let Gs and Gr be the domain and the range components, respectively, which are defined as:
G s ( i ; x ,   j ; y ) = e | ( i x ) 2 + ( j y ) 2 | 2 σ s 2
and
G r ( I ( i , j ) ; I ( x , y ) ) = e | I ( i , j ) I ( x , y ) | 2 2 σ r 2
The normalization constant K is given as:
K = 1 i = x d x + d j = y d x + d G s ( i ; x , j ; y ) G r [ I ( i , j ) , I ( x , y ) ]
Equations (4) and (5) show that the bilateral filter has three parameters: σ s 2 (the filter domain), σ r 2 (the filter range), and the third parameter is the window size N × N [15].
The geometric spread of the bilateral filter is controlled by σ s 2 . If the value of σ s 2 is increased, more neighbours are combined in the diffusion process yielding a “smoother” image, while σ r 2 represents the photometric spreading. Only pixels with a percentage difference of less than σ r 2 are processed [13].

2.4. Otsu Filtering

After passing through the bilateral filter, the image is split into its original (non-gamma corrected) Red, Green and Blue components, as shown in the block diagram in Figure 2. The kernel of the bilateral filter alters the balance of the colors in the original image in such a way to widen the differences between the color of the front and back-to-front interference. A modified version of Otsu [16] algorithm is applied to each RGB channel using the thresholds determined by the Decision Making Block, which may be considered as the “optimal” threshold for each RGB channel, and then three binary images are generated.

2.5. Image Classification

The image classification block was also trained with the synthetic images in such a way to analyze the three binary images generated in each of the channels and outputs the one that is considered the best one. This decision was also made by a naïve Bayes automatic classifier which was trained using the calculated co-occurrence matrix for each of the 32,000 synthetic images by comparing each of them with the original ground truth image, the Front image.

3. Experiments and Results

As already explained, the enormous variety of kinds of text documents makes extremely improbable that one single algorithm is able to satisfactorily binarize all kinds of documents. Depending on the nature (or degree of complexity) of the image several or no algorithm will be able to provide good results. This paper follows the assessment methodology proposed in reference [9], in which one compares the numbers of background and foreground pixels correctly matched with a ground-truth image. Twenty-three binarization algorithms were tested using the methodology described:
  • Mello-Lins [5]
  • DaSilva-Lins-Rocha [6]
  • Otsu [16]
  • Johannsen-Bille [17]
  • Kapur-Sahoo-Wong [18]
  • RenyEntropy (variation of [18])
  • Li-Tam [19]
  • Mean [20]
  • MinError [21]
  • Mixture-Modeling [22]
  • Moments [23]
  • IsoData [24]
  • Percentile [25]
  • Pun [26]
  • Shanbhag [27]
  • Triangle [28]
  • Wu-Lu [29]
  • Yean-Chang-Chang [30]
  • Intermodes [31]
  • Minimum (variation of [31])
  • Ergina-Local [32]
  • Sauvola [33]
  • Niblack [34]
A ground-truth image for each “real” world one is needed to allow a quantitative assessment of the quality of the final binary image. Only the DIBCO dataset [10] had ground-truth images available. This makes the assessment task of real-world images extremely difficult [35]. All care must be taken to guarantee the fairness of the process. The ground-truth images for the other datasets were generated by applying the 23 algorithms above and the bilateral algorithm to all the test images in the Nabuco [7] and LiveMemory [36] datasets. Visual inspection was made to choose the best binary image in a blind process, a process in which the people who selected the best image did not know which algorithm generated it. To increase the degree of fairness and the number of filtering possibilities, the three component images produced by the Decision Making block were all analyzed. The binary images chosen using the methodology above went through salt-and-pepper filtering and were used as ground-truth image for the assessment below. All the processing time figures presented in this paper are from Intel i7-4510U@ 2.00 GHzx2, 8 GB RAM, running Linux Mint 18.2 64-bit. All algorithms were coded in Java, possibly by their authors.

3.1. The Nabuco Dataset

The Nabuco bequest encompasses about 6500 letters and postcards written and typed by Joaquim Nabuco [7], totaling about 30,000 pages. Such documents are of great interest to whoever studies the history of the Americas, as Nabuco was one of the key figures in the freedom of black slaves, and was the first Brazilian Ambassador to the U.S.A. The documents of Nabuco were digitalized by the second author of this paper and the historians of the Joaquim Nabuco Foundation using a table scanner in 200 dpi resolution in true color (24 bits per pixel), back in 1992 to 1994. Due to serious storage limitations then, images were saved in the jpeg format with 1% loss. The historians in the project concluded that 150 dpi resolution would suffice to represent all the graphical elements in the documents, but choice of the 200-dpi resolution was made to be compatible with the FAX devices widely used then. About 200 of the documents in the Nabuco bequest exhibited back-to-front interference. The 15 document images used in this dataset were chosen for being representative of the diversity of documents in such a universe.
Table 1 presents the quantitative results obtained for all the documents in this dataset. P(f/f) stands for the ratio between the number of foreground pixels in the original image mapped onto black pixels and the number of black pixels in the ground-truth image. Similarly, P(b/b) is proportion between the number of background pixels in the original image mapped onto white pixels of the binary image and the number of white pixels in the ground-truth image. The figures for P(b/b) and P(f/f) are followed by “±” and the value of the standard deviation. The time corresponds to the mean processing time elapsed by the algorithm to process the images in this dataset. The results were ranked in P(b/b) decreasing order.
The results presented in Table 1 shows the bilateral filter in third place for this dataset in terms of image quality, however the standard deviation is much lower than the two first. That implies that its quality is more stable for the various document images in this dataset. Figure 5 presents the document for which the bilateral filter presented the best and the worst results in terms of image quality with two zoomed areas from the original and the binarized document.

3.2. The LiveMemory Dataset

This dataset encompasses 15 documents with 200 dpi resolution selected from the over 8000 documents from the LiveMemory project that created a digital library with all the proceedings of technical events from the Brazilian Telecommunications Society. The original proceedings were offset printed from documents either typed or electronically produced. Table 2 presents the performance results for the 12 best ranked algorithms. The bilateral filter obtained the best results in terms of image filtering. It is worth observing that in the case of the worst quality image (Figure 6, right) the performance degraded for all the algorithms. This behavior is due to the shaded area in the hard-bound spine of the volumes of the proceedings.

3.3. The DIBCO Dataset

This dataset has all the 86 images from the Digital Image Binarization Contest from 2009 to 2016. Table 3 presents the results obtained. The performance of the bilateral filter in this set may be considered good, in general. The overall performance of the bilateral filter was strongly degraded by the single image shown in Figure 7 (right) in which the P(f/f) of 25.93 drastically dropped the average result of the algorithm in this test set. It is important to remark that such an image is almost unreadable even for humans and that it degraded the performance of all the best algorithms.

4. Conclusions

Historical documents are far more difficult to binarize as several factors such as paper texture, aging, thickness, translucidity, permeability, the kind of ink, its fluidity, color, aging, etc. all may influence the performance of the algorithms. Besides all that, many historical documents were written or printed on both sides of translucent paper, giving rise to the back-to-front interference.
This paper presents a new binarization scheme based on the bilateral filter. Experiments performed in three datasets of “real world” historical documents with twenty-three other binarization algorithms. Image quality and processing time figures were provided, at least for the top 10 algorithms assessed. The results obtained showed that the proposed algorithm yields good quality monochromatic images that may compensate its high computational cost. This paper provides evidence that no binarization algorithm is an “all-kind-of-document” winner, as the performance of the algorithms varied depending of the specific features of each document. A much larger test set of synthetic about 250,000 images is currently under development, such a test set will allow much better training of the Decision Making and Image Classifier blocks of the bilateral algorithm presented. The authors are currently attempting to integrate the Decision Making and Image Classifier blocks in such a way to anticipate the choice of the best component image. This would highly improve the time performance of the proposed algorithm.
The authors of this paper are promoting a paramount research effort to assess the largest possible number of binarization algorithms for scanned documents using over 5.4 million synthetic images in the DIB-Document Image Binarization platform. An image matcher, a more general and complex version of the Decision Making block, is also being developed and trained with that large set of images, in order to whenever fed with a real world image, to be able to match with the most similar synthetic one. Once that match is made, the most suitable binarization algorithms are immediately known. If this paper were accepted, all the test images and algorithms will be included in the DIB platform. The preliminary version of the DIB-Document Image Binarization platform and website is publicly available at https://dib.cin.ufpe.br/.

Acknowledgments

The authors of this paper are grateful for the referees whose comments much helped in improving the current version of this paper and to those researchers who made the code of their algorithms publicly available for testing and performance analysis and to the DIBCO team from making their images publicly available. The authors also acknowledge the partial financial support of to CNPq and CAPES—Brazilian Government.

Author Contributions

Marcos Almeida and Rafael Dueire Lins contributed in equal proportion to the development of the algorithm presented in this paper, which was written by the latter author. Bruno Lima was responsible for the first implementation of the algorithm proposed. Rodrigo Bernardino and Darlisson Jesus re-implemented the algorithm and were also responsible for all the quality and time assessment figures presented here.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chaki, N.; Shaikh, S.H.; Saeed, K. Exploring Image Binarization Techniques; Springer: New Delhi, India, 2014. [Google Scholar]
  2. Lins, R.D. A Taxonomy for Noise in Images of Paper Documents-The Physical Noises. In Proceedings of the International Conference Image Analysis and Recognition, Halifax, NS, Canada, 6–8 July 2009; Volume 5627, pp. 844–854. [Google Scholar]
  3. Lins, R.D. An Environment for Processing Images of Historical Documents. Microprocess. Microprogr. 1995, 40, 939–942. [Google Scholar] [CrossRef]
  4. Sharma, G. Show-through cancellation in scans of duplex printed documents. IEEE Trans. Image Process. 2001, 10, 736–754. [Google Scholar] [CrossRef] [PubMed]
  5. Mello, C.A.B.; Lins, R.D. Generation of Images of Historical Documents by Composition. In Proceedings of the 2002 ACM Symposium on Document Engineering, New York, NY, USA, 8–9 November 2002; pp. 127–133. [Google Scholar]
  6. Silva, M.M.; Lins, R.D.; Rocha, V.C. Binarizing and Filtering Historical Documents with Back-to-Front Interference. In Proceedings of the 2006 ACM Symposium on Applied Computing, New York, NY, USA, 23–27 April 2006; pp. 853–858. [Google Scholar]
  7. Lins, R.D. Nabuco—Two Decades of Processing Historical Documents in Latin America. J. Univers. Comput. Sci. 2011, 17, 151–161. [Google Scholar]
  8. Roe, E.; Mello, C.A.B. Binarization of Color Historical Document Images Using Local Image Equalization and XDoG. In Proceedings of the 12th International Conference on Document Analysis and Recognition, Washington, DC, USA, 25–28 August 2013; pp. 205–209. [Google Scholar]
  9. Lins, R.D.; Almeida, M.A.M.; Bernardino, R.B.; Jesus, D.; Oliveira, J.M. Assessing Binarization Techniques for Document Images. In Proceedings of the ACM Symposium on Document Engineering, Valletta, Malta, 4–7 September 2017. [Google Scholar]
  10. Pratikakis, I.; Zagoris, K.; Barlas, G.; Gatos, B. ICDAR 2017 Competition on Document Image Binarization (DIBCO 2017). In Proceedings of the 14th IAPR International Conference on Document Analysis and Recognition, Kyoto, Japan, 13–15 November 2017; pp. 2140–2379. [Google Scholar]
  11. Efros, A.A.; Freeman, W.T. Image quilting for texture synthesis and transfer. In Proceedings of the 28th Annual Conference on Computer Graphics and Interactive Techniques (SIGGRAPH ’01), New York, NY, USA, 12–17 August 2001; pp. 341–346. [Google Scholar]
  12. Aurich, V.; Weule, J.B. Non-Linear Gaussian Filters Performing Edge Preserving Diffusion. In Proceedings of the DAGM Symposium, London, UK, 13–15 September 1995; pp. 538–545. [Google Scholar]
  13. Tomasi, C.; Manduchi, R. Bilateral Filtering for Gray and Color Images. In Proceedings of the 6th International Conference on Computer Vision, Washington, DC, USA, 4–7 January 1998; pp. 836–846. [Google Scholar]
  14. Paris, P.; Kornprobst, P.; Tumblim, J.; Durand, F. Bilateral Filtering: Theory and Applications. Found. Trends Comput. Graph. Vis. 2008, 4, 1–73. [Google Scholar] [CrossRef]
  15. Shyam Anand, C.; Sahambi, J.S. Pixel Dependent Automatic Parameter Selection for Image Denoising with Bilateral Filter. Int. J. Comput. Appl. 2012, 45, 41–46. [Google Scholar]
  16. Otsu, N. A Threshold Selection Method from Gray-Level Histograms. IEEE Trans. Syst. Man Cybern. 1979, 9, 62–66. [Google Scholar] [CrossRef]
  17. Johannsen, G.; Bille, J.A. A Threshold Selection Method Using Information Measure. In Proceedings of the 6th International Conference on Pattern Recognition (ICPR’82), Munich, Germany, 19–22 October 1982; pp. 140–143. [Google Scholar]
  18. Kapur, N.; Sahoo, P.K.; Wong, A.K.C. A New Method for Gray-Level Picture Thresholding Using the Entropy of the Histogram. Comput. Vis. Graph. Image Process. 1985, 29, 273–285. [Google Scholar] [CrossRef]
  19. Li, C.H.; Tam, P.K.S. An iterative algorithm for minimum cross entropy thresholding. Pattern Recognit. Lett. 1998, 19, 771–776. [Google Scholar] [CrossRef]
  20. Glasbey, C.A. An analysis of histogram-based thresholding algorithms. Graph. Models Image Process. 1993, 55, 532–537. [Google Scholar] [CrossRef]
  21. Kittler, J.; Illingworth, J. Minimum error thresholding. Pattern Recognit. 1986, 19, 41–47. [Google Scholar] [CrossRef]
  22. Mixture Modeling. ImageJ. Available online: http://imagej.nih.gov/ij/plugins/mixture-modeling.html (accessed on 20 January 2018).
  23. Tsai, W.H. Moment-preserving thresholding: A new approach. Comput. Vis. Graph. Image Process. 1985, 29, 377–393. [Google Scholar] [CrossRef]
  24. Doyle, W. Operation useful for similarity-invariant pattern recognition. J. Assoc. Comput. Mach. 1962, 9, 259–267. [Google Scholar] [CrossRef]
  25. Pun, T. Entropic Thresholding, A New Approach. Comput. Vis. Graph. Image Process. 1981, 16, 210–239. [Google Scholar] [CrossRef]
  26. Shanbhag, A.G.G. Utilization of Information Measure as a Means of Image Thresholding. Comput. Vis. Graph. Image Process. 1994, 56, 414–419. [Google Scholar] [CrossRef]
  27. Zack, G.W.; Rogers, W.E.; Latt, S.A. Automatic measurement of sister chromatid exchange frequency. J. Histochem. Cytochem. 1977, 25, 741–753. [Google Scholar] [CrossRef] [PubMed]
  28. Wu, U.L.; Songde, A.; Haqing, L.U.A. An Effective Entropic Thresholding for Ultrasonic Imaging. In Proceedings of the International Conference Pattern Recognition, Brisbane, Australia, 16–20 August 1998; pp. 1522–1524. [Google Scholar]
  29. Yen, J.C.; Chang, F.J.; Chang, S. A New Criterion for Automatic Multilevel Thresholding. IEEE Trans. Image Process. 1995, 4, 370–378. [Google Scholar] [PubMed]
  30. Ridler, T.W.; Calvard, S. Picture Thresholding Using an Iterative Selection Method. IEEE Trans. Syst. Man Cybern. 1978, 8, 630–632. [Google Scholar]
  31. Prewitt, M.S.; Mendelsohn, M.L. The Analysis of Cell Images. Ann. N. Y. Acad. Sci. 1996, 128, 836–846. [Google Scholar] [CrossRef]
  32. Kavallieratou, E.; Stamatatos, S. Adaptive binarization of historical document images. In Proceedings of the 18th International Conference on Pattern ICPR 2006, Hong Kong, China, 20–24 August 2006; Volume 3. [Google Scholar]
  33. Sauvola, J.; Pietikainen, M. Adaptive document image binarization. Pattern Recognit. 2000, 33, 225–236. [Google Scholar] [CrossRef]
  34. Niblack, W. An introduction to Digital Image Processing; Prentice-Hall: Upper Saddle River, NJ, USA, 1986. [Google Scholar]
  35. Ntirogiannis, K.; Gatos, B.; Pratikakis, I. Performance Evaluation Methodology for Historical Document Image Binarization. IEEE Trans. Image Process. 2013, 22, 595–609. [Google Scholar] [CrossRef] [PubMed]
  36. Lins, R.D.; Silva, G.F.P.; Torreão, G.; Alves, N.F. Efficiently Generating Digital Libraries of Proceedings with the LiveMemory Platform. In IEEE International Telecommunications Symposium; IEEE Press: Rio de Janeiro, Brazil, 2010; pp. 119–125. [Google Scholar]
Figure 1. Images with back-to-front interference from the three test sets used in this paper: Nabuco bequest (left), LiveMemory (center) and DIBCO (right).
Figure 1. Images with back-to-front interference from the three test sets used in this paper: Nabuco bequest (left), LiveMemory (center) and DIBCO (right).
Jimaging 04 00027 g001
Figure 2. Block diagram of the proposed algorithm.
Figure 2. Block diagram of the proposed algorithm.
Jimaging 04 00027 g002
Figure 3. Generation of the co-occurrence matrix for each of the images in the training set.
Figure 3. Generation of the co-occurrence matrix for each of the images in the training set.
Jimaging 04 00027 g003
Figure 4. Block diagram of the scheme for the generation of synthetic images for the Decision-Making Block.
Figure 4. Block diagram of the scheme for the generation of synthetic images for the Decision-Making Block.
Jimaging 04 00027 g004
Figure 5. Historical documents from Nabuco bequest with the best ((left)—P(f/f) = 100, P(b/b) = 99.99) and the worst ((right)—P(f/f) = 89.76, P(b/b) = 99.98) binarization results for the bilateral filter with zooms from the original (top) and binary (bottom) parts.
Figure 5. Historical documents from Nabuco bequest with the best ((left)—P(f/f) = 100, P(b/b) = 99.99) and the worst ((right)—P(f/f) = 89.76, P(b/b) = 99.98) binarization results for the bilateral filter with zooms from the original (top) and binary (bottom) parts.
Jimaging 04 00027 g005
Figure 6. Images from LiveMemory with the best ((left)—P(f/f) = 100.00, P(b/b) = 99.99) and the worst ((right)—P(f/f) = 100.00, P(b/b) = 95.97) binarization results for the bilateral filter with zooms from the original (top) and binary (bottom) parts.
Figure 6. Images from LiveMemory with the best ((left)—P(f/f) = 100.00, P(b/b) = 99.99) and the worst ((right)—P(f/f) = 100.00, P(b/b) = 95.97) binarization results for the bilateral filter with zooms from the original (top) and binary (bottom) parts.
Jimaging 04 00027 g006
Figure 7. Two documents from DIBCO dataset: (left-top) original image (left-bottom) binary image obtained using the bilateral filter best result (P(f/f) = 97.05, P(b/b) = 99.88); (right-top) original image. (right-bottom) the worst binarization results for the bilateral filter (P(f/f) = 25.93, P(b/b) = 99.99).
Figure 7. Two documents from DIBCO dataset: (left-top) original image (left-bottom) binary image obtained using the bilateral filter best result (P(f/f) = 97.05, P(b/b) = 99.88); (right-top) original image. (right-bottom) the worst binarization results for the bilateral filter (P(f/f) = 25.93, P(b/b) = 99.99).
Jimaging 04 00027 g007
Table 1. Binarization results for images from Nabuco bequest.
Table 1. Binarization results for images from Nabuco bequest.
AlgNameP(f/f)P(b/b)Time (s)
IsoData98.08 ± 3.3999.38 ± 0.600.0171
Otsu98.08 ± 3.3999.36 ± 0.630.0159
Bilateral99.57 ± 1.2399.29 ± 0.931.0790
Huang99.40 ± 2.1498.69 ± 0.880.0200
Moments99.39 ± 1.3498.40 ± 1.700.0160
Ergina-Local99.99 ± 0.0398.13 ± 0.640.3412
RenyEntropy100.0097.56 ± 1.170.0188
Kapoo-Sahoo-Wong100.0097.51 ± 1.070.0172
Yean-Chang-Chang100.0097.38 ± 1.260.0161
Triangle100.0095.94 ± 1.460.0160
Mello-Lins98.61 ± 5.1489.63 ± 24.430.0160
Mean100.0081.77 ± 5.990.0168
Johannsen-Bille98.87 ± 2.9759.77 ± 48.800.0164
Pun100.0055.44 ± 2.570.0185
Percentile100.0053.21 ± 1.330.0185
Sauvola85.51 ± 12.9399.95 ± 0.111.2977
Niblack99.75 ± 0.3477.06 ± 5.630.2135
Table 2. Binarization results for images from the LiveMemory project.
Table 2. Binarization results for images from the LiveMemory project.
AlgNameP(f/f)P(b/b)Time (s)
Bilateral100.0098.90 ± 1.073.3325
IsoData-ORIG99.56 ± 0.6998.61 ± 1.990.0734
Otsu99.60 ± 0.6898.57 ± 2.080.0735
Moments99.99 ± 0.0397.91 ± 1.870.0716
Ergina-Local98.98 ± 2.8297.62 ± 1.040.9917
Huang99.93 ± 0.2796.42 ± 4.200.0865
Triangle100.0094.24 ± 2.150.0728
Mean100.0083.58 ± 5.590.0747
Niblack99.76 ± 0.7678.31 ± 2.970.6710
Pun100.0055.28 ± 3.600.0800
Percentile100.0053.91 ± 1.960.0795
Kapur-Sahoo-Wong98.62 ± 4.9297.15 ± 1.440.0729
Table 3. Binarization results for images from Document Image Binarization Competition (DIBCO).
Table 3. Binarization results for images from Document Image Binarization Competition (DIBCO).
AlgNameP(f/f)P(b/b)Time (s)
Ergina-local91.37 ± 6.2599.88 ± 1.890.1844
RenyEntropy90.13 ± 14.1996.77 ± 3.500.0125
Yean-Chang-Chang90.61 ± 14.4496.16 ± 4.350.0112
Moments90.75 ± 9.9195.80 ± 5.190.0112
Bilateral92.99 ± 9.0690.78 ± 16.010.6099
Huang95.62 ± 6.3784.22 ± 18.360.0147
Triangle96.40 ± 5.7280.80 ± 23.320.0113
Mean99.35 ± 1.1478.99 ± 9.350.0115
MinError92.79 ± 23.4674.29 ± 19.360.0115
Pun99.68 ± 0.8256.20 ± 6.180.0122
Percentile99.71 ± 0.7255.06 ± 3.580.0121
Sauvola59.75 ± 30.0699.58 ± 0790.6933
Niblack95.91 ± 2.3178.61 ± 5.690.1241

Share and Cite

MDPI and ACS Style

Almeida, M.; Lins, R.D.; Bernardino, R.; Jesus, D.; Lima, B. A New Binarization Algorithm for Historical Documents. J. Imaging 2018, 4, 27. https://doi.org/10.3390/jimaging4020027

AMA Style

Almeida M, Lins RD, Bernardino R, Jesus D, Lima B. A New Binarization Algorithm for Historical Documents. Journal of Imaging. 2018; 4(2):27. https://doi.org/10.3390/jimaging4020027

Chicago/Turabian Style

Almeida, Marcos, Rafael Dueire Lins, Rodrigo Bernardino, Darlisson Jesus, and Bruno Lima. 2018. "A New Binarization Algorithm for Historical Documents" Journal of Imaging 4, no. 2: 27. https://doi.org/10.3390/jimaging4020027

APA Style

Almeida, M., Lins, R. D., Bernardino, R., Jesus, D., & Lima, B. (2018). A New Binarization Algorithm for Historical Documents. Journal of Imaging, 4(2), 27. https://doi.org/10.3390/jimaging4020027

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop