第四章 图像增强
本科目的重点在于研究图像处理的办法,因此首先将图像简化,用没有颜色的单通道的灰度图来表示图像,仅在最后一节讨论了彩色图像的模型。
灰度是一种表示图像亮度或灰暗程度的测量标准。在数字图像中,灰度通常以数字值的形式表示,范围通常在0到255之间。而书本对灰度的定义则采用了概率,首先对灰度进行了归一化处理,每个像素取
[
0
,
1
]
[0,1]
[0,1]区间内的灰度
r
k
r_k
rk,用
P
(
r
k
)
P(r_k)
P(rk)表示灰度的概率分布。
假设一个变换后的图像,则设变换后灰度为
s
s
s,有
s
=
T
(
r
)
,
r
=
T
−
1
(
s
)
s=T(r), r=T^{-1}(s)
s=T(r),r=T−1(s)
取随机变量
η
,
ξ
\eta,\xi
η,ξ分别表示新图像s,原图像p,则根据随机变量的函数可以写出
F
η
(
s
)
=
P
(
η
<
s
)
=
P
[
ξ
<
r
]
F_\eta (s)=P(\eta < s)=P[\xi < r]
Fη(s)=P(η<s)=P[ξ<r]
为了让变换是可逆的,且保持各像素的相对位置不变,要求
T
(
r
)
T(r)
T(r)单调递增,即对
有
因此
F
η
(
s
)
=
P
(
η
<
s
)
=
P
[
ξ
<
r
]
=
∫
−
∞
r
p
r
(
x
)
d
x
F_\eta (s)=P(\eta < s)=P[\xi < r] = \int\limits_{-\infin}^{r}p_r(x)dx
Fη(s)=P(η<s)=P[ξ<r]=−∞∫rpr(x)dx
对
s
s
s求导得密度函数
p
s
(
s
)
p_s(s)
ps(s):
p
s
(
s
)
=
d
r
d
s
⋅
d
[
∫
−
∞
r
p
r
(
x
)
d
x
]
d
r
=
p
r
(
r
)
⋅
d
r
d
s
=
p
r
(
r
)
⋅
d
[
T
−
1
(
s
)
]
d
s
p_s(s)=\frac{dr}{ds}\cdot \frac{d[\int\limits_{-\infin}^{r}p_r(x)dx]}{dr}=p_r(r)\cdot \frac{dr}{ds}=p_r(r)\cdot \frac{d[T^{-1}(s)]}{ds}
ps(s)=dsdr⋅drd[−∞∫rpr(x)dx]=pr(r)⋅dsdr=pr(r)⋅dsd[T−1(s)]
4.1 直方图均衡化处理
直方图均衡化是指使图像的整体对比度更均匀,从而提高图像的视觉质量。因此,需要利用直方图获取图片的统计信息。
4.1.1 直方图均衡化原理
- 变换函数为 s = T ( r ) = ∫ 0 r p r ( ω ) d ω s=T(r) =\int\limits_0^r p_r(\omega)d\omega s=T(r)=0∫rpr(ω)dω
- 对 r r r求导得到 d s d r = p r ( r ) \frac{ds}{dr}=\mathcal{p}_r(r) drds=pr(r) 反转得到 d r d s = 1 p r ( r ) \frac{dr}{ds}=\frac{1}{\mathcal{p}_r(r)} dsdr=pr(r)1
- 代入 p s ( s ) = p r ( r ) ⋅ d r d s p_s(s)=p_r(r)\cdot \frac{dr}{ds} ps(s)=pr(r)⋅dsdr 即 p s ( s ) = p r ( r ) ⋅ 1 p r ( r ) = 1 p_s(s)=p_r(r)\cdot \frac{1}{\mathcal{p}_r(r)} = 1 ps(s)=pr(r)⋅pr(r)1=1
- 得到 s s s是在 [ 0 , 1 ] [0,1] [0,1]上的均匀分布
直方图均衡化就是产生一副灰度分布具有均匀概率密度的图像,是直接利用灰度的累积分布函数,可以看作是灰度直方图的内在性质。
这样处理后的图像是一幅灰度层次较为适中,比原始图像清晰,明快得多的图像,认为这样的图像是质量得到了改善。
例题1
4.1.2 直方图均衡化计算
对于数字图像处理而言,灰度值都是离散值,因此要用
∑
\sum
∑取代
∫
d
x
\int dx
∫dx
即将函数
p
r
(
r
)
p_r(r)
pr(r)改为分布:
P
r
(
r
k
)
=
n
k
N
,
0
≤
r
k
≤
1
P_r(r_k)=\frac{n_k}{N},0\leq r_k \leq 1
Pr(rk)=Nnk,0≤rk≤1
累计分布
∫
0
r
p
r
(
x
)
d
x
\int\limits_{0}^{r} p_r(x) dx
0∫rpr(x)dx改为:
T
(
r
k
)
=
S
k
=
∑
i
=
0
k
P
r
(
r
j
)
=
∑
j
=
0
k
n
j
N
T(r_k) = S_k = \sum\limits_{i=0}^{k} P_r(r_j)=\frac{\sum\limits_{j=0}^{k}n_j}{N}
T(rk)=Sk=i=0∑kPr(rj)=Nj=0∑knj
例题
首先计算灰度的分布列
而后由 S k = ∑ i = 0 k P r ( r j ) = ∑ j = 0 k n j N S_k = \sum\limits_{i=0}^{k} P_r(r_j)=\frac{\sum\limits_{j=0}^{k}n_j}{N} Sk=i=0∑kPr(rj)=Nj=0∑knj计算 S 0 … S 7 S_0\dots S_7 S0…S7,并修正:
r i → s i r_i \to s_i ri→si就表示转换后的效果:
处理前后直方图对比如下:
s的直方图高度几乎在一条水平线上,也即“均衡化”处理
使用均衡化处理之后,灰度级的数量通常会减少,这种现象叫做简并现象,主要来自于求
s
k
s_k
sk时使用了舍入归为近似值(以免出现不允许的灰度值)。
因此,减少简并的方法就是增加像素的比特数,允许的灰度值更多,减少灰度层次的损失。
4.2 图像平滑化处理
无例题,略 无例题,略 无例题,略
4.3 图像尖锐化处理
只介绍微分尖锐化处理。
微分尖锐化属于空域法,用于增强图像边缘和细节,基本思想是通过对图像应用导数运算来突出像素值的变化
无例题,略
无例题,略
无例题,略
4.4 彩色图像处理
人们对图像进行分析时,人眼能区别的灰度层次 大约只有二十几种,但却能够识别成千上万的色
彩,彩色图像处理被划分为三个主要领域,即
- 真彩色(True color processing):被处理的图像一般从全彩色传感器中获得,例如彩色摄影机或彩色扫描仪。
- 假彩色 (False color processing): 假彩色并不使用颜色通道,而是将不同的图像信息映射到伪造的颜色上。例如,将高度信息映射到颜色,从而用颜色来表示高度。
- 伪彩色 (pseudo color processing):伪彩色通过将单通道的灰度图像映射到彩色图像上,以提供更多视觉信息。通常,通过使用渐变的颜色映射将灰度值映射到不同的颜色,以增强图像的可视化。
4.4.1 颜色模型
通常使用的多数彩色模型或者是面向硬件设备(例如彩色监视器或打印机),或者是面向应用的。
通常与硬件有关的模型有:
- RGB 模型:这种模型用在彩色监视器和彩色摄相机等领域。
- CMY 模型:用在彩色打印机上。
- YIQ 模型:用于彩色电视广播。该模型中,Y表示亮度,而I和Q是正交的两个颜色分量。
- HSI 模型:常用于彩色图像模拟。
在图像处理中常用的模型是RGB、YIQ和HSI模型。
4.4.2 RGB模型
RGB 模型使用红色(R)、绿色(G)、蓝色(B)三个通道来表示颜色,书本规定每个通道的强度常在 [ 0 , 1 ] [0,1] [0,1]的范围内。
书本对灰度、颜色等指标的定义都采取的是 [ 0 , 1 ] [0,1] [0,1]的随机变量,而不是现实应用中 [ 0 , 255 ] [0,255] [0,255]的离散值。
可能是考虑到可以采用更高或更低的分辨率或比特数,才将这一维度抽象出来。
优势:
- RGB图像由三个独立的图像平面构成,当原图像本身(如摄像机)就是由三原色(光谱数据)描述时,数字图像处理中运用RGB模型就很有意义。
- 对人脸的彩色图像进行增强处理,部分图像隐藏在阴影中,可以对它进行直方图均衡化,把每个图像平面单独地进行直方图均衡,所有可能隐藏在阴影中的图像部分都将被增强。
劣势:
- 三个图像平面的强度改变的同时也将不同地改变颜色性能(如色调),显示在RGB 监视器上的图片就不再是自然和谐的了。
4.4.3 CMY模型
RGB到CMY的转换表示为:
[
C
M
Y
]
=
[
1
1
1
]
−
[
R
G
B
]
\begin{bmatrix} C\\ M \\ Y \end{bmatrix} = \begin{bmatrix} 1\\ 1 \\ 1 \end{bmatrix} - \begin{bmatrix} R\\ G \\ B \end{bmatrix}
CMY
=
111
−
RGB
4.4.4 YIQ模型
YIQ 模型包括亮度(Y)、色度(I、Q)三个分量。亮度(Y)表示图像的亮度信息,而色度(I、Q)表示颜色信息。
RGB到YIQ的转换表示为:
[
Y
I
Q
]
=
[
0.299
0.587
0.114
0.596
−
0.275
−
0.321
0.212
−
0.523
0.311
]
[
R
G
B
]
\begin{bmatrix} Y\\ I \\ Q \end{bmatrix} = \begin{bmatrix} 0.299 & 0.587 & 0.114\\ 0.596 & -0.275 & -0.321 \\ 0.212 & -0.523 & 0.311 \end{bmatrix} \begin{bmatrix} R\\ G \\ B \end{bmatrix}
YIQ
=
0.2990.5960.2120.587−0.275−0.5230.114−0.3210.311
RGB
优势:
- YIQ 模型利用人的可视系统对亮度变化比对色调和饱和度变化更敏感而设计,表示Y时 给予较大的带宽,表示I 、Q 时赋予较小的带宽。
- 去掉了亮度( Y )和颜 色信息( I 和Q )间的紧密联系,使得处理图像的亮度成分时能在不影响颜色成分的情况下进行。
4.4.5 HSI模型
HSI 模型包括**色相(Hue)、饱和度(Saturation)、亮度(Intensity)**三个分量。色相表示颜色的基本色调,饱和度表示颜色的纯度,亮度表示颜色的明亮度。
RGB到HSI的转换表示为:
H
=
c
o
s
−
1
{
1
2
[
(
R
−
G
)
+
(
R
−
B
)
]
[
(
R
−
G
)
2
+
(
R
−
B
)
(
G
−
B
)
]
1
2
}
H=cos^{-1}\lbrace \frac{\frac{1}{2} [(R-G)+(R-B)]}{[(R-G)^2 + (R-B)(G-B)]^{\frac{1}{2}}} \rbrace
H=cos−1{[(R−G)2+(R−B)(G−B)]2121[(R−G)+(R−B)]}
S
=
1
−
3
⋅
m
i
n
(
R
,
G
,
B
)
(
R
+
G
+
B
)
S = 1-\frac{3\cdot min(R,G,B)}{(R+G+B)}
S=1−(R+G+B)3⋅min(R,G,B)
I
=
1
3
(
R
+
G
+
B
)
I=\frac{1}{3}(R+G+B)
I=31(R+G+B)
色相 H 的范围是
[
0
°
,
360
°
)
[0\degree, 360\degree)
[0°,360°),饱和度 S 和亮度 I 的范围是 [0, 1]。
模型的颜色分量(色调(hue)和饱和度 (saturation))的定义与三原色模型有关。颜色点的
H
H
H是该向量与红色轴的夹角。因此,当
H
=
0
°
H=0\degree
H=0° 时,为红色,
H
=
120
°
H=120\degree
H=120°时,为绿色,240为蓝色,60为黄色。
图片的亮度由垂直于三角形并通过其中心的直线表示,从下到上由暗变亮。