图像复原方法综述
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
图像复原方法综述
1、摘要
图像是人类视觉的基础,给人具体而直观的作用。图像的数字化包括取样和量化两个步骤。数字图像处理就是将图像信号转换成数字格式,并利用计算机进行加工和处理的过程。
图像复原是图像处理中的一个重要问题,对于改善图像质量具有重要的意义。解决该问题的关键是对图像的退化过程建立相应的数学模型,然后通过求解该逆问题获得图像的复原模型并对原始图像进行合理估计。
本文主要介绍了图像退化的原因、图像复原技术的分类和目前常用的几种图像复原方法,详细的介绍了维纳滤波、正则滤波、LR算法和盲区卷积,并通过实验证明了该方法的可行性和有效性。
关键词:图像退化、图像复原、维纳滤波、正则滤波、LR算法、盲区卷积、
2、图像复原概述
在图像的获取、传输以及保存过程中,由于各种因素,如大气的湍流效应、摄像设备中光学系统的衍射、传感器特性的非线性、光学系统的像差、成像设备与物体之间的相对运动、感光胶卷的非线性及胶片颗粒噪声以及电视摄像扫描的非线性等所引起的几何失真,都难免会造成图像的畸变和失真。通常,称由于这些因素引起的质量下降为图像退化。
图像退化的典型表现是图像出现模糊、失真,出现附加噪声等。由于图像的退化,在图像接受端显示的图像已不再是传输的原始图像,图像效果明显变差。为此,必须对退化的图像进行处理,才能恢复出真实的原始图像,这一过程就称为图像复原[1]。
图像复原技术是图像处理领域中一类非常重要的处理技术,与图像增强等其他基本图像处理技术类似,也是以获取视觉质量某种程度的改善为目的,所不同的是图像复原过程实际上是一个估计过程,需要根据某些特定的图像退化模型,对退化图像进行复原。简言之,图像复原的处理过程就是对退化图像品质的提升,并通过图像品质的提升来达到图像在视觉上的改善。
由于引起图像退化的因素众多,且性质各不相同,目前没有统一的复原方法,众多研究人员根据不同的应用物理环境,采用了不同的退化模型、处理技巧和估计准则,从而得到了不同的复原方法。
图像复原算法是整个技术的核心部分。目前,国内在这方面的研究才刚刚起步,而国外却已经取得了较好的成果。早期的图像复原是利用光学的方法对失真的观测图像进行校正,
而数字图像复原技术最早则是从对天文观测图像的后期处理中逐步发展起来的。其中一个成功例子是NAS A的喷气推进实验室在1964年用计算机处理有关月球的照片。照片是在空间飞行器上用电视摄像机拍摄的,图像的复原包括消除干扰和噪声,校正几何失真和对比度损失以及反卷积。另一个典型的例子是对肯尼迪遇刺事件现场照片的处理。由于事发突然,照片是在相机移动过程中拍摄的,图像复原的主要目的就是消除移动造成的失真[2]。
早期的复原方法有:非邻域滤波法,最近邻域滤波法以及效果较好的维纳滤波和最小二乘滤波等。随着数字信号处理和图像处理的发展,新的复原算法不断出现,在应用中可以根据具体情况加以选择。
目前国内外图像复原技术的研究和应用主要集中于诸如空间探索、天文观测、物质研究、遥感遥测、军事科学、生物科学、医学影象、交通监控、刑事侦察等领域。如生物方面,主要是用于生物活体细胞内部组织的三维再现和重构,通过复原荧光显微镜所采集的细胞内部逐层切片图,来重现细胞内部构成;医学方面,如对肿瘤周围组织进行显微观察,以获取肿瘤安全切缘与癌肿原发部位之间关系的定量数据;天文方面,如采用迭代盲反卷积进行气动光学效应图像复原研究等。
3、图像退化模型
图像复原问题的有效性关键之一取决于描述图像退化过程模型的精确性。要建立图像的退化模型,则首先必须了解、分析图像退化的机理并用数学模型表现出来。在实际的图像处理过程中,图像均需以数字离散函数表示,所以必须将退化模型离散化
[3]
。
对于退化图像),(y x g :
⎰⎰+∞∞-+∞
∞
-+--=
)
,(),(),(),(y x n d d y x h f y x g βαβαβα
(1)
如果上式中f ,h ,n ,g 按相同间隔采样,产生相应的阵列[]AB j i f ),(、[]CD j i h ),(、
[]AB j i n ),(、[]AB j i g ),(,然后将这些阵列补零增广得到大小为N M ⨯的周期延拓阵列,为了
避免重叠误差,这里1-+≥C A M ,1-+≥D B N 。由此,当k=0,1,L,M-1;l =0,1,L,N-1时,即可得到二维离散退化模型形式:
∑∑-=-=+--=
101
),(),(),(),(M i N j e
e
e
e l k n j l i k h j i
f l k g
(2)
如果用矩阵表示上式,则可写为:
n
Hf g +=
(3)
其中,f ,g ,n 为一个行堆叠形成的1⨯MN 列向量,H 为MN MN ⨯阶的块循环矩阵。 现实中造成图像降质的种类很多,常见的图像退化模型及点扩展函数有如下情景[15]: (1) 线性移动降质
在拍照时,成像系统与目标之间有相对直线移动会造成图像的降质。水平方向线性移动可以用以下降质函数来描述:
⎪⎩⎪
⎨⎧=≤≤=其他
若0
01),(n and d m d
n m h
(4)
式中,d 是降质函数的长度。在应用中如果线性移动降质函数不在水平方向,则可类似地定义移动降质函数。
(2) 散焦降质
当镜头散焦时,光学系统造成的图像降质相应的点扩展函数是一个均匀分布的圆形光斑。此时,降质函数可表示为:
⎪⎩⎪
⎨⎧=+=其他
若0
R 1),(2
222
n m R n m h π
(5)
式中,R 是散焦半径。
(3) 高斯(Gauss)降质
Gau ss降质函数是许多光学测量系统和成像系统最常见的降质函数。对于这些系统,决定系统点扩展函数的因素比较多。众多因素综合的结果总是使点扩展函数趋于Gau ss 型。典型的系统可以举出光学相机和CCD 摄像机、γ相机、CT 相机、成像雷达、显微光学系统等。G aus s降质函数可以表达为:
⎩
⎨
⎧∈+-=其他若0),()](exp[),(22C n m n m K n m h α
(6)
式中,K 是归一化常数,α是一个正常数,C 是),(n m h 的圆形支持域。 4、几种较经典的复原方法介绍