Niblack 和 Sauvola 阈值化

Niblack 和 Sauvola 阈值是局部阈值化技术，它们对于背景不均匀的图像很有用，特别是对于文本识别 [1], [2]。与为整个图像计算单个全局阈值不同，通过使用考虑局部邻域（由以像素为中心的窗口定义）的均值和标准差的特定公式，为每个像素计算多个阈值。

在这里，我们使用这些算法将图像二值化，并将其与常见的全局阈值化技术进行比较。参数 window_size 确定包含周围像素的窗口的大小。

[1]

Niblack, W (1986), An introduction to Digital Image Processing, Prentice-Hall.

[2]

J. Sauvola 和 M. Pietikainen，“自适应文档图像二值化”，模式识别 33(2)，第 225-236 页，2000 年。 DOI:10.1016/S0031-3203(99)00055-2

import matplotlib
import matplotlib.pyplot as plt

from skimage.data import page
from skimage.filters import threshold_otsu, threshold_niblack, threshold_sauvola


matplotlib.rcParams['font.size'] = 9


image = page()
binary_global = image > threshold_otsu(image)

window_size = 25
thresh_niblack = threshold_niblack(image, window_size=window_size, k=0.8)
thresh_sauvola = threshold_sauvola(image, window_size=window_size)

binary_niblack = image > thresh_niblack
binary_sauvola = image > thresh_sauvola

plt.figure(figsize=(8, 7))
plt.subplot(2, 2, 1)
plt.imshow(image, cmap=plt.cm.gray)
plt.title('Original')
plt.axis('off')

plt.subplot(2, 2, 2)
plt.title('Global Threshold')
plt.imshow(binary_global, cmap=plt.cm.gray)
plt.axis('off')

plt.subplot(2, 2, 3)
plt.imshow(binary_niblack, cmap=plt.cm.gray)
plt.title('Niblack Threshold')
plt.axis('off')

plt.subplot(2, 2, 4)
plt.imshow(binary_sauvola, cmap=plt.cm.gray)
plt.title('Sauvola Threshold')
plt.axis('off')

plt.show()