9. 图像调整:转换图像内容#
9.1. 颜色操作#
用于操作颜色通道的大多数函数都可以在子模块 skimage.color
中找到。
9.1.1. 颜色模型之间的转换#
彩色图像可以使用不同的颜色空间来表示。最常见的颜色空间之一是RGB 空间,其中图像具有红色、绿色和蓝色通道。然而,其他颜色模型也被广泛使用,例如HSV 颜色模型,其中色调、饱和度和值是独立的通道,或者用于打印的CMYK 模型。
skimage.color
提供了实用函数,用于将图像转换为不同的颜色空间以及从不同的颜色空间转换。可以通过转换操作将整数类型的数组转换为浮点类型。
>>> # bright saturated red
>>> red_pixel_rgb = np.array([[[255, 0, 0]]], dtype=np.uint8)
>>> color.rgb2hsv(red_pixel_rgb)
array([[[ 0., 1., 1.]]])
>>> # darker saturated blue
>>> dark_blue_pixel_rgb = np.array([[[0, 0, 100]]], dtype=np.uint8)
>>> color.rgb2hsv(dark_blue_pixel_rgb)
array([[[ 0.66666667, 1. , 0.39215686]]])
>>> # less saturated pink
>>> pink_pixel_rgb = np.array([[[255, 100, 255]]], dtype=np.uint8)
>>> color.rgb2hsv(pink_pixel_rgb)
array([[[ 0.83333333, 0.60784314, 1. ]]])
9.1.2. 从 RGBA 到 RGB 的转换 - 通过 alpha 混合移除 alpha 通道#
通过将 RGBA 图像与背景进行 alpha 混合,将其转换为 RGB 图像,可以使用 rgba2rgb()
实现
>>> import skimage as ski
>>> img_rgba = ski.data.logo()
>>> img_rgb = ski.color.rgba2rgb(img_rgba)
9.1.3. 颜色值和灰度值之间的转换#
可以使用 rgb2gray()
将 RGB 图像转换为灰度图像
>>> img = ski.data.astronaut()
>>> img_gray = ski.color.rgb2gray(img)
rgb2gray()
使用颜色通道的非均匀加权,因为人眼对不同颜色的敏感度不同。因此,这种加权确保了从 RGB 到灰度的亮度保留。
>>> red_pixel = np.array([[[255, 0, 0]]], dtype=np.uint8)
>>> ski.color.rgb2gray(red_pixel)
array([[ 0.2125]])
>>> green_pixel = np.array([[[0, 255, 0]]], dtype=np.uint8)
>>> ski.color.rgb2gray(green_pixel)
array([[ 0.7154]])
使用 gray2rgb()
将灰度图像转换为 RGB 图像时,只是将灰度值复制到三个颜色通道上。
9.1.4. 图像反转#
反转图像也称为互补图像。对于二值图像,True 值变为 False,反之亦然。对于灰度图像,像素值将替换为数据类型的最大值和实际值的差值。对于 RGB 图像,每个通道都执行相同的操作。此操作可以使用 skimage.util.invert()
实现。
>>> import skimage as ski
>>> img = ski.data.camera()
>>> inverted_img = ski.util.invert(img)
9.1.5. 用标签绘制图像#
可以使用 label2rgb()
在灰度图像上叠加颜色,使用标签数组来编码要用相同颜色表示的区域。
9.2. 对比度和曝光#
图像像素可以采用由图像的 dtype
确定的值(请参阅图像数据类型及其含义),例如 uint8
图像的 0 到 255,或者浮点图像的 [0, 1]
。然而,大多数图像要么具有较窄的值范围(因为对比度较差),要么大多数像素值集中在可访问值的子范围内。skimage.exposure
提供了将强度值扩展到更大范围的函数。
第一类方法计算强度的非线性函数,该函数独立于特定图像的像素值。这些方法通常用于校正传感器或接收器(如人眼)的已知非线性。一个著名的例子是伽马校正,在 adjust_gamma()
中实现。
其他方法根据图像的直方图重新分配像素值。像素值的直方图使用 skimage.exposure.histogram()
计算
>>> import numpy as np
>>> import skimage as ski
>>> image = np.array([[1, 3], [1, 1]])
>>> ski.exposure.histogram(image)
(array([3, 0, 1]), array([1, 2, 3]))
histogram()
返回每个值箱的像素数和箱的中心。因此,histogram()
的行为与 numpy.histogram()
的行为略有不同,后者返回箱的边界。
最简单的对比度增强 rescale_intensity()
包括使用线性变换将像素值拉伸到整个允许范围。
>>> import skimage as ski
>>> text = ski.data.text()
>>> text.min(), text.max()
(10, 197)
>>> better_contrast = ski.exposure.rescale_intensity(text)
>>> better_contrast.min(), better_contrast.max()
(0, 255)
即使图像使用整个值范围,有时在值范围的末端也只有很少的权重。在这种情况下,使用图像的百分位数裁剪像素值可以提高对比度(以损失一些信息为代价,因为某些像素会因该操作而饱和)
>>> moon = ski.data.moon()
>>> v_min, v_max = np.percentile(moon, (0.2, 99.8))
>>> v_min, v_max
(10.0, 186.0)
>>> better_contrast = ski.exposure.rescale_intensity(moon, in_range=(v_min, v_max))
函数 equalize_hist()
将像素值的累积分布函数 (cdf) 映射到线性 cdf,确保值范围的所有部分在图像中得到同等表示。因此,在对比度较差的较大区域中,细节会得到增强。作为进一步的改进,可以使用 equalize_adapthist()
在图像的子区域中执行直方图均衡化,以校正整个图像的曝光梯度。请参阅示例 直方图均衡化。