- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第五章பைடு நூலகம்图 像 复 原 与 重 建
1.退化模型
2.代数恢复方法 3.频率域恢复方法 4.几何校正 5.图像重建
数字图像处理
电子信息与自动化学院
1
第五章
图像复原与重建
什么是图像复原? 什么是图像重建? 数字图像如何进行几何变换(缩放、旋转等)
数字图像处理
电子信息与自动化学院
2
5.1 退化模型
g Hf
g、f都是M维列向量,H是M×M阶矩阵,矩阵中的每一行 元素均相同,只是每行以循环方式右移一位,因此矩阵H 是循环矩阵。循环矩阵相加或相乘得到的还是循环矩阵。
数字图像处理 电子信息与自动化学院
18
5.1.2 退化的数学模型
二维离散模型 设输入的数字图像f(x, y)大小为A×B,点扩展函数h(x, y)被均 匀采样为C×D大小。为避免交叠误差,仍用添零扩展的方法, 将它们扩展成M=A+C-1和N=B+D-1个元素的周期函数。
数字图像处理
电子信息与自动化学院
12
5.1.2 退化的数学模型
退化的数学模型
f (x, y) n (x, y)
h(x,y)
g (x, y)
在时域
g ( x, y) f ( x, y) * h( x, y) n( x, y)
f ( , )h( x , y )dd n( x, y)
二维离散退化模型同样可以表示为:
g Hf
式中,g、 f是MN×1维列向量,H是MN×MN维矩阵。其方法 是将g(x, y)和f(x, y)中的元素排成列向量。
数字图像处理 电子信息与自动化学院
20
5.1.2 退化的数学模型
H 0 H M 1 H M 2 H M 1 H1 H 0 H H M 1 H M 2 H M 3
f ( x) f e ( x) 0 h( x ) he ( x) 0
0 x A 1 A x M 1
0 x B 1 B x M 1
M 1 m 0
ge ( x ) f e ( x ) he ( x ) f e (m)he ( x m)
若把噪声考虑进去, 则离散图像退化模型为
g e ( x, y )
m 0
数字图像处理
M 1 N 1 n 0
f (m, n)h ( x m, y n) n ( x, y )
e e e
电子信息与自动化学院
22
5.1.2 退化的数学模型
写成矩阵形式为
g Hf n
上述线性空间不变退化模型表明,在给定了g(x, y),并且知 道退化系统的点扩展函数h(x, y)和噪声分布n(x, y)的情况下,可
19
5.1.2 退化的数学模型
则输出的降质数字图像为
g e ( x, y)
m 0
M 1 N 1 n 0
f (m, n)h ( x m, y n) f ( x, y) * h( x, y)
e e
式中:x=0, 1, 2, …, M-1; y=0, 1, 2, …, N-1。
H1 H2 H0
Hi(i=0, 1, 2,…, M-1)为子矩阵,大小为N×N,即H矩阵由M×M 个大小为N×N的子矩阵组成, 称为分块循环矩阵。分块矩阵是
由延拓函数he(x, y)的第j行构成的,构成方法如下:
数字图像处理 电子信息与自动化学院
21
5.1.2 退化的数学模型
f ( x, y ) 0 x A 1且0 y B 1 f e ( x, y ) 其他 0 h ( x , y ) 0 x C 1 且 0 y D 1 he ( x, y ) 其他 0
数字图像处理 电子信息与自动化学院
分析退化原因
数字图像处理
建立退化模型
反向推演
恢复图像
10
电子信息与自动化学院
5.1.1 退化 图像复原与图像增强
目的: 都是为了改善图像的质量。 不同之处:
图像复原是试图利用退化过程的先验知识使已退化的图像 恢复本来面目,即根据退化的原因, 分析引起退化的环境因 素,建立相应的数学模型, 并沿着使图像降质的逆过程恢复 图像。从图像质量评价的角度来看, 图像复原就是提高图像 的可理解性。
数字图像处理
电子信息与自动化学院
17
5.1.2 退化的数学模型
ge ( x ) f e ( x ) he ( x ) f e (m)he ( x m)
m 0 M 1
g (0) he (0) g (1) he (1) g ( 2) h ( 2) e g ( M 1) he ( M 1)
而图像增强不考虑图像如何退化,只通过试探各种技术来 提高视觉效果,图像增强的过程基本上是一个探索的过程, 它利用人的心理状态和视觉系统去控制图像质量, 直到人们 的视觉系统满意为止。
数字图像处理 电子信息与自动化学院
11
5.1.1 退化
评价准则(对图像复原结果的评价) 最小均方准则 加权均方准则 最大熵准则
结论
f (x, y) n (x, y)
h(x,y)
g (x, y)
离散图像的时域数学模型
g ( x, y) f ( x, y) * h( x, y) n( x, y)
可写为
g Hf n
g、 f、n是MN×1维列向量,H是MN×MN维矩阵
数字图像处理 电子信息与自动化学院
24
5.2 代数恢复法
16
电子信息与自动化学院
5.1.2 退化的数学模型
因为he(x)的周期为M,所以he(x)=he(x+M),即
he ( 1) he ( M 1) he ( 2) he ( M 2) he ( 3) he ( M 3) he ( M 1) he (1)
he (0) he (1) H he (2) h ( M 1) e he (1) he (0) he (2) he (1) he (3) he ( M 2) he (0) he ( M 1)
电子信息与自动化学院
7
5.1.1 退化
造成退化的常见因素
…
数字图像处理
电子信息与自动化学院
8
5.1.1 退化
产生退化的具体原因 光学系统的像差 光学成像衍射 成像系统的非线性畸变 摄影胶片的感光的非线性 成像过程的相对运动 大气的湍流效应 环境随机噪声
数字图像处理 电子信息与自动化学院
式中,x=0, 1, 2, …, M-1。
数字图像处理 电子信息与自动化学院
15
5.1.2 退化的数学模型
因为fe(x)和he(x)已扩展成周期函数,故ge(x)也是周期性函数, 用矩阵表示为
g (0) he (0) g (1) he (1) g ( 2) h ( 2) e g ( M 1) he ( M 1)
1 称为逆滤波器 H u, v
数字图像处理 电子信息与自动化学院
27
5.3.1 逆滤波恢复法
逆滤波恢复法的基本步骤
(1)对退化图像g(x,y)作二维离散傅立叶变换,得到G(u,v); (2)计算系统点扩散函数h(x,y)的二维傅立叶变换,得到
9
5.1.1 退化
克服退化的措施——图像复原
采取一定的方法尽可能地减少或消除图像质 量的下降,恢复图像的本来面目,这就是图像复 原,也称为图像恢复。 典型图像复原是利用退化现象的某种先验知 识,建立退化现象的数学模型,再根据模型进行 反向的推演运算,以恢复原来的景物图像。
图像复原的一般过程
he ( 1) he (0) he (1) he ( M
he ( M 1) f e (0) he ( M 2) f e (1) he ( M 3) f e (2) f ( M 1 ) 2) he (0) e
g x, y f x, y hx, y nx, y
傅立叶变换
Gu, v H u, vF u, v N u, v
G u, v N u, v F u, v H u, v H u, v
恢复原图像
F u, v Gu, v H u, v 不考虑噪声
估计出原始图像f(x, y)。
假 设 图 像 大 小 M=N=512 , 相 应 矩 阵 H 的 大 小 为
MN×MN=262 144×262 144,这意味着要解出f (x, y)需要解262
144个联立方程组,其计算量十分惊人。
数字图像处理 电子信息与自动化学院
23
5.1.2 退化的数学模型
he ( j,0) he ( j,1) Hj he ( j, N 1)
he ( j, N 1) he ( j, N 2) he ( j,1) he ( j,0) he ( j, N 1) he ( j,0) he ( j , N 2) he ( j, N 3) he ( j,0)
5.1.1退化
退化的概念 图像在形成、传输和记录过程中,由于成 像系统、传输介质和设备的不完善,使图像的 质量下降,这一过程称为图像的退化。
核心:图像质量下降 原因:成像系统、传输介质和设备不完善 产生环节:形成、传输和记录