第五章图像复原讲解
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
V2 SNR 10 log10 s V2 n
灰度对比度 C SNR ob 噪声均方差
2
2
信噪比, 狭义来讲是指放大器的输出信号的电压与同时输出的噪声电压的比,常常用 分贝数表示。设备的信噪比越高表明它产生的杂音越少。
行周期扩展为周期为M ≥ A + B – 1。
f ( x) fe ( x) 0
0 ≤ x ≤ A 1 A≤ x ≤ M 1
h( x ) he ( x ) 0
0 ≤ x ≤ B 1 B ≤ x ≤ M 1
那么它们的时域离散卷积可定义为下式：
g e ( x)
此式说明，如果系统H对冲激函数的响应为已知，则对任意输入的响应可 用上式求得，即，线性系统H完全可以由冲激响应来表征。图像中冲激响应
也称为点扩散函数。
在有噪音的情况下：
4.离散函数的退化模型 对和进行均匀取样后，就可引伸出离散函数的退化模型。用一维的来 说明。如果f (x)和h(x)周期分别A和B的序列，为避免卷积周期重叠需要对它们进
2
2 5/ 6
]
其中，C是与湍流性质有关的常数。 （2）光学散焦转移函数：
H (u, v) J1 (d )
d
(u 2 v 2 )1/ 2
其中，d是散焦点扩展函数的直径, J1(•) 是第一类贝塞尔函数。
(3)照相机与景物相对运动 设T为快门时间，x0(t)，y0(t)是位移的x分量和y分量。
H1 f e (0) ne (0) H2 f e (1) ne (1) f ( MN 1) n ( MN 1) H0 e e
H0 H1 g Hf n H M 1
2
ˆ n n g Hf
T
2
ˆ ) T ( g Hf ˆ) ( g Hf
n
2
ˆ gH f
2
fˆ 2 最小，即求 J ( fˆ ) g H × fˆ 2 求 n 2最小等效于求 g H × 的极小值问题。这里选择 fˆ 除了要求 J ( fˆ ) 为最小外，不受任何 其它条件约束，因此称为非约束复原。求 J ( fˆ ) 的极小值方法就 是一般的求极值的方法。把 J ( fˆ ) 对 fˆ 微分，并使结果为 0，即:
如果 如果 其他
za z b
4）噪声的形成 高斯噪声：电子噪声、弱光照/温度条件下的传感器噪声 瑞利分布：深度成像、超声波图像 指数和Gamma分布：激光成像 椒盐噪声：快速瞬变、误切换 周期噪声：图像采集过程中的电子或电磁干扰
图像中噪声的概率密度函数举例1：
原图
图像和其直方图
图像中噪声的概率密度函数举例2：
5.退化参数的估计 退化参数的估计包括噪声估计和点扩展函数的估计。下面主要介绍点扩 散函数的估计。 1）运用先验知识估计 大气湍流、光学系统散焦、照相机与景物相对运动等，根据导致模糊的 物理过程（先验知识）来确定h(x,y)或H(u,v)。 (1) 长时间曝光下大气湍流造成的转移函数：
H (u, v) exp[c u v
原图
图像和其直方图
图像中噪声的概率密度函数举例3：
图像中的周期噪声
5）噪声参数估计
(1) 周期噪声的参数估计 一般可以通过图像的频谱进行估计；特殊情况下可以直接从图像中 噪声分量的周期性进行推断（简单情形）。 (2) 一般噪声参数的估计 可以根据所采用的传感器类型进行噪声分布的部分推断；通常通过 特定的成像安排进行估计当只有已采集到的图像时，一般通过图像中的 平滑区域进行PDF参数的估计。如下图：
2）常见噪声 热噪声：白噪声（频率覆盖整个频谱均匀） 高斯噪声（幅度符合高斯分布） 闪烁噪声：具有反比于频率（1/f）的频谱； 粉色噪声（在对数频率间隔内有相同的能量） 发射噪声：高斯分布（电子运动的随机性）
3）噪声的概率密度函数
噪声为随机变量，用概率密度来刻画。
(1) 高斯噪声
(z )2 p( z ) exp 2 2π 2 1
推广到2-D 扩展：
f ( x, y ) f e ( x, y ) 0 h( x, y ) he ( x, y ) 0
0 ≤ x ≤ A 1 A ≤ x ≤ M 1 0 ≤ x ≤ C 1
和 或 和
0 ≤ y ≤ B 1 B ≤ y ≤ N 1 0 ≤ y ≤ D 1 D ≤ y ≤ N 1
5.2 图像退化的数学模型
1.退化模型示意图
n ( x, y )
f ( x, y )
H
g ( x, y )
其中H为退化过程，n(x, y)为加性噪声（统计特性已知）。 2. 系统H的基本定义 就一般而言，系统是某些元件或部件以某种方式构造而成的整体。系统 本身所具有的某些特性就构成了通过系统的输入信号与输出信号的某种联系。 系统的分类可有：线性系统和非线性系统，时变系统和非时变系统，集总参 数系统和分布参数系统，连续系统和离散系统。 1）线性系统：是具有均匀性和相加性的系统 2）时不变系统：满足各个参数不随时间变化。 3）空间不变系统：满足H [ f ( x a , y b )] g ( x a , y b )
H (u, v) exp j 2 (ux0 (t ) vy0 (t )dt
T 0
点扩散函数示图：
(a) 光脉冲 (b) 退化的光脉冲
5.3 无约束恢复
1.无约束和有约束复原
由退化模型得： n g Hf 最小均方误差准则： n 1) 无约束复原：
在最小二乘方意义上说，希望找到一个 f 使下式达到最小:
M 1 m 0
f e (m)he ( x m)
x 0, 1, , M 1
显然，上式也是具有周期M的序列。 如果用矩阵来表示上述离散退化模型，可写成下式之形式:
用矩阵形式表示
g e (0) he (0) g e (1) he (1) g Hf g ( M 1) h ( M 1) e e he ( 1) he (0) he ( M 2) he ( M 1) he ( M 2) he (0) f e ( 0) f e (1) f ( M 1) e
其处理方法如下：
ˆ ，使下述准则函数为最小 寻找一个 f J( f ) Q × fˆ
2
( g H fˆ
2
n
2
)
式中 为一常数，是拉格朗日系数。加上约束条件后，就可以按一般 求极小值的方法进行求解。将上式对 fˆ 微分，并使结果为零，则有:
ˆ) J ( f ˆ) 0 2H ' ( g Hf ˆ f
ˆ H 'g H ' Hf
ˆ (H ' H ) 1 H ' g f
因为 H 是一方阵，并且设 H
-1
存在，则可求得
f
:
ˆ H 1 ( H ' ) 1 H ' g H 1 g f
这种方法要求知道成象系统的表达式H。 2）有约束复原方法 在最小二乘方复原处理中，为了在数学上更容易处理，常常附加某 种约束条件。 例如，可以令 Q 为 f 的线性算子，那么，最小二乘方复原问题可看 ˆ 2 的函数，服从约束条件 2 × 成是使形式为 Q × ˆ 2 的最小化 f n g H f 问题。而这种有附加条件的极值问题可用拉格朗日乘数法来处理。
实际上，大部分系统是非线性和空间变化的，但以这样的模型处理起 来困难很大，一般都简化为线性的非时变和非空间变化的近似模型进行处 理。这样近似的优点是使线性系统理论中的许多理论可以直接用来解决图 象复原问题。 3.连续函数的退化模型 设系统H对坐标为(a,b)处的冲激函数(x-a,y-b)的冲激响应为h(x,a,y,b)， 则
(2) 均匀噪声
1 /(b a) 如果 a ≤ z ≤ b p( z ) 其他 0
其均值和方差为：
( a b) / 2
2 (b a ) 2 / 12
(3) 脉冲噪声 噪声脉冲可以是正的或负的，一 般假设a和b。都是“饱和”值双极性脉 冲噪声也称椒盐噪声。
Pa p ( z ) Pb 0
数字图像处理
(Digital Image Processing）
第五章
• • • • • 5.1 5.2 5.3 5.4 5.5
图像复原
退化及噪声 图像退化的数学模型 无约束恢复 有约束恢复 交互式恢复
5.1 退化及噪声
1. 图象复原的概念
1）图像复原的定义 图像复原也称图象恢复，是图象处理中的一大类技术。所谓图像复原， 是指去除或减轻在获取数字图像过程中发生的图像质量下降（退化）这些退 化包括由光学系统、运动等等造成图像的模糊，以及源自电路和光度学因素 的噪声。 图像复原的目标是对退化的图像进行处理，使它趋向于复原成没有退化 的理想图像。成像过程的每一个环节（透镜，感光片，数字化等等）都会引 起退化。 在进行图像复原时，既可以用连续数学，也可以用离散数学进行处理。 其次，处理既可在空间域，也可在频域进行。 2）图象恢复与图象增强的异同
由于的h(x，y)的周期性，使得[H]成为一个块循环矩阵。每一 块如下所示：
he (i,0) he (i,1) Hi h (i, N 1) e
he (i, N 1) he (i,0) he (i, N 2)
he (i,1) he (i,2) he (i,0)
H M 1 H0 H M 2
H1 f e (0) ne (0) H2 f e (1) ne (1) f ( MN 1) n ( MN 1) H0 e e
根据周期性：he(x) = he(x+M)
he ( M 1) he (0) h (1) he (0) e H he ( M 1) he ( M 2)
he (1) he (2) he (0)
由于的h(x)周期性，使得[H]成为一个循环矩阵。
聚焦不准（失焦，限制了图象锐度） 模糊（限制频谱宽度）
噪声（是一个统计过程）
抖动（机械、电子） 3. 图象复原方法分类 按采用的技术可分为：无约束和有约束 按采用的策略可分为：自动和交互 按采用的处理所在域可分为：频域和空域
图像退化举例1：
图像退化举例2：
4.噪声及其特性 噪声是最常见的退化因素之一，对信号来说，噪声是一种外部干扰。 但噪声本身也是一种信号（携带了噪声源的信息）。 1）关于噪声的度量 人们常只关心噪声的强度 ，可用信噪比(signal-to-noise ratio,SNR)、 能量比(电压平方比) 等来描述。分别表示为：
相同点：改进输入图象的视觉质量 。
不同点：图象增强目的是取得较好的视觉结果(不考虑退化原因)； 图象 恢复根据相应的退化模型和知识重建或恢复原始的图象(考虑退化原因) 。
2. 图象退化的原因 图象退化指由场景得到的图象没能完全地反映场景的真实内容，产生了 失真等问题。其原因是多方面的。如：
透镜象差/色差
C ≤ x ≤ Mwk.baidu.com1 或
则退化过程为：
g e ( x, y )
M 1 N 1 m 0 n 0
f e (m, n)he ( x m, y n)
H M 1 H0 H M 2
x 0, 1, , M 1 y 0, 1, , N 1
H0 H1 g Hf n H M 1
灰度对比度 C SNR ob 噪声均方差
2
2
信噪比, 狭义来讲是指放大器的输出信号的电压与同时输出的噪声电压的比,常常用 分贝数表示。设备的信噪比越高表明它产生的杂音越少。
行周期扩展为周期为M ≥ A + B – 1。
f ( x) fe ( x) 0
0 ≤ x ≤ A 1 A≤ x ≤ M 1
h( x ) he ( x ) 0
0 ≤ x ≤ B 1 B ≤ x ≤ M 1
那么它们的时域离散卷积可定义为下式：
g e ( x)
此式说明，如果系统H对冲激函数的响应为已知，则对任意输入的响应可 用上式求得，即，线性系统H完全可以由冲激响应来表征。图像中冲激响应
也称为点扩散函数。
在有噪音的情况下：
4.离散函数的退化模型 对和进行均匀取样后，就可引伸出离散函数的退化模型。用一维的来 说明。如果f (x)和h(x)周期分别A和B的序列，为避免卷积周期重叠需要对它们进
2
2 5/ 6
]
其中，C是与湍流性质有关的常数。 （2）光学散焦转移函数：
H (u, v) J1 (d )
d
(u 2 v 2 )1/ 2
其中，d是散焦点扩展函数的直径, J1(•) 是第一类贝塞尔函数。
(3)照相机与景物相对运动 设T为快门时间，x0(t)，y0(t)是位移的x分量和y分量。
H1 f e (0) ne (0) H2 f e (1) ne (1) f ( MN 1) n ( MN 1) H0 e e
H0 H1 g Hf n H M 1
2
ˆ n n g Hf
T
2
ˆ ) T ( g Hf ˆ) ( g Hf
n
2
ˆ gH f
2
fˆ 2 最小，即求 J ( fˆ ) g H × fˆ 2 求 n 2最小等效于求 g H × 的极小值问题。这里选择 fˆ 除了要求 J ( fˆ ) 为最小外，不受任何 其它条件约束，因此称为非约束复原。求 J ( fˆ ) 的极小值方法就 是一般的求极值的方法。把 J ( fˆ ) 对 fˆ 微分，并使结果为 0，即:
如果 如果 其他
za z b
4）噪声的形成 高斯噪声：电子噪声、弱光照/温度条件下的传感器噪声 瑞利分布：深度成像、超声波图像 指数和Gamma分布：激光成像 椒盐噪声：快速瞬变、误切换 周期噪声：图像采集过程中的电子或电磁干扰
图像中噪声的概率密度函数举例1：
原图
图像和其直方图
图像中噪声的概率密度函数举例2：
5.退化参数的估计 退化参数的估计包括噪声估计和点扩展函数的估计。下面主要介绍点扩 散函数的估计。 1）运用先验知识估计 大气湍流、光学系统散焦、照相机与景物相对运动等，根据导致模糊的 物理过程（先验知识）来确定h(x,y)或H(u,v)。 (1) 长时间曝光下大气湍流造成的转移函数：
H (u, v) exp[c u v
原图
图像和其直方图
图像中噪声的概率密度函数举例3：
图像中的周期噪声
5）噪声参数估计
(1) 周期噪声的参数估计 一般可以通过图像的频谱进行估计；特殊情况下可以直接从图像中 噪声分量的周期性进行推断（简单情形）。 (2) 一般噪声参数的估计 可以根据所采用的传感器类型进行噪声分布的部分推断；通常通过 特定的成像安排进行估计当只有已采集到的图像时，一般通过图像中的 平滑区域进行PDF参数的估计。如下图：
2）常见噪声 热噪声：白噪声（频率覆盖整个频谱均匀） 高斯噪声（幅度符合高斯分布） 闪烁噪声：具有反比于频率（1/f）的频谱； 粉色噪声（在对数频率间隔内有相同的能量） 发射噪声：高斯分布（电子运动的随机性）
3）噪声的概率密度函数
噪声为随机变量，用概率密度来刻画。
(1) 高斯噪声
(z )2 p( z ) exp 2 2π 2 1
推广到2-D 扩展：
f ( x, y ) f e ( x, y ) 0 h( x, y ) he ( x, y ) 0
0 ≤ x ≤ A 1 A ≤ x ≤ M 1 0 ≤ x ≤ C 1
和 或 和
0 ≤ y ≤ B 1 B ≤ y ≤ N 1 0 ≤ y ≤ D 1 D ≤ y ≤ N 1
5.2 图像退化的数学模型
1.退化模型示意图
n ( x, y )
f ( x, y )
H
g ( x, y )
其中H为退化过程，n(x, y)为加性噪声（统计特性已知）。 2. 系统H的基本定义 就一般而言，系统是某些元件或部件以某种方式构造而成的整体。系统 本身所具有的某些特性就构成了通过系统的输入信号与输出信号的某种联系。 系统的分类可有：线性系统和非线性系统，时变系统和非时变系统，集总参 数系统和分布参数系统，连续系统和离散系统。 1）线性系统：是具有均匀性和相加性的系统 2）时不变系统：满足各个参数不随时间变化。 3）空间不变系统：满足H [ f ( x a , y b )] g ( x a , y b )
H (u, v) exp j 2 (ux0 (t ) vy0 (t )dt
T 0
点扩散函数示图：
(a) 光脉冲 (b) 退化的光脉冲
5.3 无约束恢复
1.无约束和有约束复原
由退化模型得： n g Hf 最小均方误差准则： n 1) 无约束复原：
在最小二乘方意义上说，希望找到一个 f 使下式达到最小:
M 1 m 0
f e (m)he ( x m)
x 0, 1, , M 1
显然，上式也是具有周期M的序列。 如果用矩阵来表示上述离散退化模型，可写成下式之形式:
用矩阵形式表示
g e (0) he (0) g e (1) he (1) g Hf g ( M 1) h ( M 1) e e he ( 1) he (0) he ( M 2) he ( M 1) he ( M 2) he (0) f e ( 0) f e (1) f ( M 1) e
其处理方法如下：
ˆ ，使下述准则函数为最小 寻找一个 f J( f ) Q × fˆ
2
( g H fˆ
2
n
2
)
式中 为一常数，是拉格朗日系数。加上约束条件后，就可以按一般 求极小值的方法进行求解。将上式对 fˆ 微分，并使结果为零，则有:
ˆ) J ( f ˆ) 0 2H ' ( g Hf ˆ f
ˆ H 'g H ' Hf
ˆ (H ' H ) 1 H ' g f
因为 H 是一方阵，并且设 H
-1
存在，则可求得
f
:
ˆ H 1 ( H ' ) 1 H ' g H 1 g f
这种方法要求知道成象系统的表达式H。 2）有约束复原方法 在最小二乘方复原处理中，为了在数学上更容易处理，常常附加某 种约束条件。 例如，可以令 Q 为 f 的线性算子，那么，最小二乘方复原问题可看 ˆ 2 的函数，服从约束条件 2 × 成是使形式为 Q × ˆ 2 的最小化 f n g H f 问题。而这种有附加条件的极值问题可用拉格朗日乘数法来处理。
实际上，大部分系统是非线性和空间变化的，但以这样的模型处理起 来困难很大，一般都简化为线性的非时变和非空间变化的近似模型进行处 理。这样近似的优点是使线性系统理论中的许多理论可以直接用来解决图 象复原问题。 3.连续函数的退化模型 设系统H对坐标为(a,b)处的冲激函数(x-a,y-b)的冲激响应为h(x,a,y,b)， 则
(2) 均匀噪声
1 /(b a) 如果 a ≤ z ≤ b p( z ) 其他 0
其均值和方差为：
( a b) / 2
2 (b a ) 2 / 12
(3) 脉冲噪声 噪声脉冲可以是正的或负的，一 般假设a和b。都是“饱和”值双极性脉 冲噪声也称椒盐噪声。
Pa p ( z ) Pb 0
数字图像处理
(Digital Image Processing）
第五章
• • • • • 5.1 5.2 5.3 5.4 5.5
图像复原
退化及噪声 图像退化的数学模型 无约束恢复 有约束恢复 交互式恢复
5.1 退化及噪声
1. 图象复原的概念
1）图像复原的定义 图像复原也称图象恢复，是图象处理中的一大类技术。所谓图像复原， 是指去除或减轻在获取数字图像过程中发生的图像质量下降（退化）这些退 化包括由光学系统、运动等等造成图像的模糊，以及源自电路和光度学因素 的噪声。 图像复原的目标是对退化的图像进行处理，使它趋向于复原成没有退化 的理想图像。成像过程的每一个环节（透镜，感光片，数字化等等）都会引 起退化。 在进行图像复原时，既可以用连续数学，也可以用离散数学进行处理。 其次，处理既可在空间域，也可在频域进行。 2）图象恢复与图象增强的异同
由于的h(x，y)的周期性，使得[H]成为一个块循环矩阵。每一 块如下所示：
he (i,0) he (i,1) Hi h (i, N 1) e
he (i, N 1) he (i,0) he (i, N 2)
he (i,1) he (i,2) he (i,0)
H M 1 H0 H M 2
H1 f e (0) ne (0) H2 f e (1) ne (1) f ( MN 1) n ( MN 1) H0 e e
根据周期性：he(x) = he(x+M)
he ( M 1) he (0) h (1) he (0) e H he ( M 1) he ( M 2)
he (1) he (2) he (0)
由于的h(x)周期性，使得[H]成为一个循环矩阵。
聚焦不准（失焦，限制了图象锐度） 模糊（限制频谱宽度）
噪声（是一个统计过程）
抖动（机械、电子） 3. 图象复原方法分类 按采用的技术可分为：无约束和有约束 按采用的策略可分为：自动和交互 按采用的处理所在域可分为：频域和空域
图像退化举例1：
图像退化举例2：
4.噪声及其特性 噪声是最常见的退化因素之一，对信号来说，噪声是一种外部干扰。 但噪声本身也是一种信号（携带了噪声源的信息）。 1）关于噪声的度量 人们常只关心噪声的强度 ，可用信噪比(signal-to-noise ratio,SNR)、 能量比(电压平方比) 等来描述。分别表示为：
相同点：改进输入图象的视觉质量 。
不同点：图象增强目的是取得较好的视觉结果(不考虑退化原因)； 图象 恢复根据相应的退化模型和知识重建或恢复原始的图象(考虑退化原因) 。
2. 图象退化的原因 图象退化指由场景得到的图象没能完全地反映场景的真实内容，产生了 失真等问题。其原因是多方面的。如：
透镜象差/色差
C ≤ x ≤ Mwk.baidu.com1 或
则退化过程为：
g e ( x, y )
M 1 N 1 m 0 n 0
f e (m, n)he ( x m, y n)
H M 1 H0 H M 2
x 0, 1, , M 1 y 0, 1, , N 1
H0 H1 g Hf n H M 1