图像复原方法综述Word文档格式.doc
《图像复原方法综述Word文档格式.doc》由会员分享,可在线阅读,更多相关《图像复原方法综述Word文档格式.doc(11页珍藏版)》请在冰豆网上搜索。
简言之,图像复原的处理过程就是对退化图像品质的提升,并通过图像品质的提升来达到图像在视觉上的改善。
由于引起图像退化的因素众多,且性质各不相同,目前没有统一的复原方法,众多研究人员根据不同的应用物理环境,采用了不同的退化模型、处理技巧和估计准则,从而得到了不同的复原方法。
图像复原算法是整个技术的核心部分。
目前,国内在这方面的研究才刚刚起步,而国外却已经取得了较好的成果。
早期的图像复原是利用光学的方法对失真的观测图像进行校正,而数字图像复原技术最早则是从对天文观测图像的后期处理中逐步发展起来的。
其中一个成功例子是NASA的喷气推进实验室在1964年用计算机处理有关月球的照片。
照片是在空间飞行器上用电视摄像机拍摄的,图像的复原包括消除干扰和噪声,校正几何失真和对比度损失以及反卷积。
另一个典型的例子是对肯尼迪遇刺事件现场照片的处理。
由于事发突然,照片是在相机移动过程中拍摄的,图像复原的主要目的就是消除移动造成的失真[2]。
早期的复原方法有:
非邻域滤波法,最近邻域滤波法以及效果较好的维纳滤波和最小二乘滤波等。
随着数字信号处理和图像处理的发展,新的复原算法不断出现,在应用中可以根据具体情况加以选择。
目前国内外图像复原技术的研究和应用主要集中于诸如空间探索、天文观测、物质研究、遥感遥测、军事科学、生物科学、医学影象、交通监控、刑事侦察等领域。
如生物方面,主要是用于生物活体细胞内部组织的三维再现和重构,通过复原荧光显微镜所采集的细胞内部逐层切片图,来重现细胞内部构成;
医学方面,如对肿瘤周围组织进行显微观察,以获取肿瘤安全切缘与癌肿原发部位之间关系的定量数据;
天文方面,如采用迭代盲反卷积进行气动光学效应图像复原研究等。
3、图像退化模型
图像复原问题的有效性关键之一取决于描述图像退化过程模型的精确性。
要建立图像的退化模型,则首先必须了解、分析图像退化的机理并用数学模型表现出来。
在实际的图像处理过程中,图像均需以数字离散函数表示,所以必须将退化模型离散化[3]。
对于退化图像:
(1)
如果上式中,,,按相同间隔采样,产生相应的阵列、、、,然后将这些阵列补零增广得到大小为的周期延拓阵列,为了避免重叠误差,这里,。
由此,当k=0,1,L,M-1;
l=0,1,L,N-1时,即可得到二维离散退化模型形式:
(2)
如果用矩阵表示上式,则可写为:
(3)
其中,,,为一个行堆叠形成的列向量,H为阶的块循环矩阵。
现实中造成图像降质的种类很多,常见的图像退化模型及点扩展函数有如下情景[15]:
(1)线性移动降质
在拍照时,成像系统与目标之间有相对直线移动会造成图像的降质。
水平方向线性移动可以用以下降质函数来描述:
(4)
式中,d是降质函数的长度。
在应用中如果线性移动降质函数不在水平方向,则可类似地定义移动降质函数。
(2)散焦降质
当镜头散焦时,光学系统造成的图像降质相应的点扩展函数是一个均匀分布的圆形光斑。
此时,降质函数可表示为:
(5)
式中,R是散焦半径。
(3)高斯(Gauss)降质
Gauss降质函数是许多光学测量系统和成像系统最常见的降质函数。
对于这些系统,决定系统点扩展函数的因素比较多。
众多因素综合的结果总是使点扩展函数趋于Gauss型。
典型的系统可以举出光学相机和CCD摄像机、相机、CT相机、成像雷达、显微光学系统等。
Gauss降质函数可以表达为:
(6)
式中,K是归一化常数,是一个正常数,C是的圆形支持域。
4、几种较经典的复原方法介绍
图像复原算法有线性和非线性两类。
线性算法通过对图像进行逆滤波来实现反卷积,这类方法方便快捷,无需循环或迭代,直接可以得到反卷积结果,然而,它有一些局限性,比如无法保证图像的非负性。
而非线性方法通过连续的迭代过程不断提高复原质量,直到满足预先设定的终止条件,结果往往令人满意。
但是迭代程序导致计算量很大,图像复原时耗较长,有时甚至需要几个小时。
所以实际应用中还需要对两种处理方法综合考虑,进行选择[4]。
(1)维纳滤波法
维纳滤波法是由Wiener首先提出的,应用于一维信号处理,取得了很好的效果。
之后,维纳滤波法被用于二维信号处理,也取得了不错的效果,尤其在图像复原领域,由于维纳滤波计算量小,复原效果好,从而得到了广泛的应用和发展。
维纳滤波器寻找一个使统计误差函数
(7)
最小的估计。
E是期望值操作符,是未退化的图像。
该表达式在频域可表示为
(8)
其中,
表示退化函数
表示的复共轭
表示噪声的功率谱
表示未退化图像的功率谱
比率称为信噪功率比。
在IPT中维纳滤波使用函数deconvwnr来实现的。
模拟实验结果如下:
(2)正则滤波法
另一个容易实现线性复原的方法称为约束的最小二乘方滤波,在IPT中称为正则滤波,并且通过函数deconvreg来实现。
在最小二乘复原处理中,常常需要附加某种约束条件。
例如令Q为f的线性算子,那么最小二乘方复原的问题可以看成使形式为的函数,服从约束条件的最小化问题,这种有附加条件的极值问题可以用拉格朗日乘数法来处理。
寻找一个,使下述准则函数为最小:
(9)
式中叫拉格朗日系数。
通过指定不同的Q,可以得到不同的复原目标。
模拟实验结果如下:
(3)Lucy-Richardson算法
L-R算法是一种迭代非线性复原算法,它是从最大似然公式印出来的,图像用泊松分布加以模型化的。
当下面这个迭代收敛时模型的最大似然函数就可以得到一个令人满意的方程:
(10)
*代表卷积,代表未退化图像的估计,g和h和以前定义一样。
在IPT中,L-R算法由名为deconvlucy的函数完成的。
(4)盲去卷积
在图像复原过程中,最困难的问题之一是,如何获得PSF的恰当估计。
那些不以PSF为基础的图像复原方法统称为盲区卷积。
它以MLE为基础的,即一种用被随机噪声所干扰的量进行估计的最优化策略。
工具箱通过函数deconvblind来执行盲区卷积。
图像分别迭代5次,10次,20次的结果。
5、参考文献
[1]罗军辉冯平等。
MATLAB7.0在图像处理中的应用[M],北京:
机械工业出版社
2006,202~204,257~278。
[2]刘维一,于德月,王肇圻等。
用迭代法消除数字图像放大后的模糊[J]。
光电子.
激光,2002,13(4):
398~400。
[3]陶洪。
数字共焦显微技术及其图像复原算法研究四川大学硕士学位论文
CNKI:
:
CDMD:
10610.2.2003.6632。
[4]冈萨雷斯等著。
数字图像处理(MATLAB版)[M],电子工业出版社2009。
6、程序附录清单
(1)维纳滤波
clc
clearall
closeall
I=imread('
C:
\DocumentsandSettings\Administrator\桌面\原始图.jpg'
);
ubplot(231);
Imshow(I);
I=rgb2gray(I);
subplot(232);
imshow(I);
noise=0.1*randn(size(I));
PSF=fspecial('
motion'
21,11);
Blurred=imfilter(I,PSF,'
circular'
BlurredNoisy=im2uint8(Blurred);
NP=abs(fftn(noise)).^2;
NPOW=sum(NP(:
)/prod(size(noise)));
NCORR=fftshift(real(ifftn(NP)));
IP=abs(fftn(I)).^2;
IPOW=sum(IP(:
ICORR=fftshift(real(ifftn(IP)));
ICORR1=ICORR(:
ceil(size(I,1)/2));
NSR=NPOW/IPOW;
subplot(233);
imshow(BlurredNoisy,[]);
title('
A=BlurredandNoisy'
subplot(234);
imshow(deconvwnr(BlurredNoisy,PSF),[]);
deconbwnr(A,PSF,NSR)'
subplot(235);
imshow(deconvwnr(BlurredNoisy,PSF,NCORR,ICORR),[]);
deconbwnr(A,PSF,NCORR,ICORR)'
subplot(236);
imshow(deconvwnr(BlurredNoisy,PSF,NPOW,ICORR1),[]);
deconbwnr(A,PSF,NPOW,ICORR_1_D)'
(2)正则滤波
subplot(231);
gaussian'