Otsu thresholdingΒΆ

This example illustrates automatic Otsu thresholding.

../../../_images/plot_threshold_1.png

Python source code: plot_threshold.py

import matplotlib.pyplot as plt
from skimage import data
try:
from skimage import filters
except ImportError:
from skimage import filter as filters
from skimage import exposure
camera = data.camera()
val = filters.threshold_otsu(camera)
hist, bins_center = exposure.histogram(camera)
plt.figure(figsize=(9, 4))
plt.subplot(131)
plt.imshow(camera, cmap='gray', interpolation='nearest')
plt.axis('off')
plt.subplot(132)
plt.imshow(camera < val, cmap='gray', interpolation='nearest')
plt.axis('off')
plt.subplot(133)
plt.plot(bins_center, hist, lw=2)
plt.axvline(val, color='k', ls='--')
plt.tight_layout()
plt.show()

Total running time of the example: 0.30 seconds ( 0 minutes 0.30 seconds)