收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

图像反卷积算法研究

杨航  
【摘要】:近半个世纪以来,图像反卷积算法始终是数字图像处理领域研究的热点之一.图像反卷积算法的研究是一个重要但又困难的问题,也是具有理论挑战性的分支以及具有重要的实际应用价值.图像反卷积作为一种改善图像质量的技术,从失真和降质的观测量来估计原来不失真图像,尽量降低和消除图像的失真和噪声.在图像复原的问题中,造成图像模糊的点扩散函数总是一个低通滤波器,这就使输入图像的高频成分受到抑制甚至丧失.反卷积是一个逆过程,目的是要”找回”丢失了的高频成分.不难理解,反卷积过程包含在观测中的噪声会被放大,这意味着反卷积的结果可能偏离真实的解.为了获得尽可能真实的解,在反卷积方法上要加进抑制噪声的考虑,这就需要在复原图像和噪声放大之间做出适当的折衷. 一.正六边形点扩散函数图像反卷积 现有图像复原方法的研究中,点扩散函数有有限的几种,这可以简化‘大类点扩散函数的问题.关于图像退化点扩散函数的三种常用模型:线性运动、圆盘离焦和高斯模型.除了线性运动模型,另外两种模型也都是基于圆形光圈建立起来的.有些图像复原的研究工作虽然超出了这三种模型,但仍然要求点扩散函数具有圆对称性质. 我们注意到,随着科技的不断进步,具有可变光圈(一般是正多边形光圈)的光学系统越来越广泛地应用于各行各业.为了适应客观实际需求,有必要在图像反卷积方法的研究中,根据这一类降质图像的成像机理,研究相应的反卷积方法,以提高图像复原的精度和效率.而在现有文献中,关于可变(正多边形)光圈降质图像的反卷积方法的研究至今未见相关报导.本论文的主要研究之一是对正六边形光圈离焦图像反卷积方法的研究. 图像的降质模型可以写成如下形式:y(l,s)=(h*x)(t,s)+γ(t,s)(1)其中,y,h,x和r分别表示观测图像,真实图像,点扩散函数和噪声.图像复原的目的就是通过观测图像y(t,s)来得到真实图像x(t,s)的一个精确估计.根据Fourier变换的卷积定理,上述模型在频域的表达为:Y(υ,υ)=H(υ,υ)X(υ,υ)+Γ(υ,υ)(2)其中Y,H,X,Γ分别是y,h,x,r的Fourier变换. 对于处理正六边形可变光圈造成的图像离焦,本文的贡献有如下几个方面: 1)基于正六边形光圈光学系统的成像机理,提出了正六边形离焦模型(3.2节):其中a为正六边形的边长,Rθ=(?). 该点扩散函数模型刻画了正六边形光圈离焦的本质,对于精确估计正六边形光圈离焦点扩散函数是必要的. 2)本文通过分析和推导,我们从理论上提出并证明了正六边形离焦点扩散函数Fourier变换H(u,υ)的零点分布定理(3.3节). 第一圈零点分布定理(定理3.7):当0x(?)时,f(x)有且仅有一个零点. 第二圈零点分布定理(定理3.14):f(x)的第二个零点在区间(?)≤x≤(?)上.在这个区间上f(x)有且仅有一个零点. 其中f(x)是H(u,υ)的一种数学变形. 这两个定理从理论上证明了正六边形光圈离焦点扩散函数在频域存在零点,并可以对零点具体位置的搜索提供指导.根据这两个定理,用数值方法先一次性计算出H(u,υ)在频域中心附近的第一圈零点集合L1和第二圈零点集合L2. 3)在上面工作的基础上,本文提出了正六边形光圈离焦图像快速复原方法(3.4节),方法的主要步骤如下: 为了得到点扩散函数,利用观测图像y(t,s)的Fourier变换和H(u,υ)的第一圈、第二圈零点集合L,可以辨识出正六边形点扩散函数发生离焦现象时正六边形的边长a和偏转角度θ.接下来的问题就是一个已知点扩散函数的非盲反卷积问题了.为了算法的执行效率,加之图像的信噪比(SNR)比较高,可以采用Wiener滤波的方法.具体步骤如下: a).计算观测图像y(t,s)的Fourier变换Y(u,υ). b).求解问题得到点扩散函数的两个参数的估计:正六边形的边长a和偏转角度θ.并由此得到点扩散函数的估计Ha,θ(υ,υ)=H(n(υ,υ)Rθ). c).利用Wiener滤波(或者其他有效的非盲图像反卷积算法),通过Ha,θ来估计真实图像. 应用此方法对于有较高信噪比(SNR)的图像进行了图像复原实验,得到比较满意的效果. 二.保持纹理结构的图像反卷积算法 在自然图像中,纹理特征广泛的存在.例如指纹图像,地震剖面图像等等,这些图像中包含着丰富的弯曲的振荡纹理和定向纹理信息.由于纹理特征是一种高频信息,而图像模糊过程是一个低通滤波的过程,因此在图像反卷积时,就需要将这种重要的图像特征恢复出来.现今的绝大多数的反卷积算法,采用的先验约束主要是小波、曲波、全变差等,由于它们对纹理信息表达的局限性,因此这些算法难以将图像的纹理信息很好的复原.这些图像反卷积算法往往会将图像中的纹理丢失,造成图像失真和视觉上的不自然.为了很好的解决纹理图像复原这一类特殊而且重要的问题,就需要采用适合纹理信息表达的数学工具.波原子变换是新进提出的一种半尺度半方向变换,它的重要特性就是能够较其他多尺度工具更能稀疏的表达振荡纹理和定向纹理特征.由于这一重要特点,很自然的就会想到将波原子这一工具应用到纹理保持的图像处理中. 在一些基于多尺度变换的图像反卷积算法中,我们发现得到的复原图像包含一些振铃,产生这一现象的原因是阈值收缩时导致的振铃效应.同样,振铃效应也出现在了基于波原子的反卷积算法中.因此为了即能够保证图像的纹理信息在反卷积过程中得到很好的保持,又能够减少处理后的振铃效应,我们提出: 1).提出了一种基于波原子变换的图像反卷积算法.这个算法可以有效的处理纹理信息丰富的图像,而且能够很好的保持图像的纹理结构. 对于图像反卷积而言,Fourier域是一个基本的选择,因为根据卷积定理有:Y(κ1, κ2)=H(κ1,κ2)X(κ1,κ2)+Γ(κ1,κ2)(6)Y,H,X,Γ分别是y,h,x,r的离散Fourier变换. a). Fourier正则化反卷积(FoRD).图像FoRD估计x。的Fourier形式为其中,其中H表示H的复共轭,(?)(k1,k2)≥0是正则项,用来控制收缩的幅度. b).波原子域Wiener滤波. 对于含有有色噪声的图像x。,波原子域的Wiener收缩估计为:其中cα,μ是xα在下标为μ的波原子变换系数,cα,μh是波原子域硬阈值估计xαh的波原子变换系数,σα,μ2是波原子域下的噪声方差,β是正则化参数.计算cα,μw的逆波原子变换,这样我们就得到了基于波原子域Wiener滤波估计xαw. 2).提出联合非局部均值滤波来解决振铃效应,并且依然能够很好的保持图像的纹理信息.使用联合非局部均值滤波去处理由波原子反卷积得到的图像,这种方法有效的解决了由波原子收缩导致的振铃现象. a).使用估计图像xαw作为联合非局部均值滤波器中的参考图像.估计图像xαw保留了图像大部分的纹理特征和细节,所以联合非局部均值滤波能够有效的改进图像质量,得到的估计图像xJ按如下公式计算: b).加权估计.如果直接使用联合非局部均值滤波,那么在输出的反卷积图像中会出现一些噪声圆斑,尤其是当噪声水平很高时(噪声方差很大时).为了抑制这些圆斑和保持处理上的简单有效,我们在联合非局部均值滤波估计的基础上加上基于波原子的反卷积估计.也就是说,最终的反卷积结果x*由下式得到:χ*=βχαω+(1-β)χJ,β∈[0,1)(11)这个加权方法同时也可以平衡波原子估计和联合非局部均值估计,改进图像质量.


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 张凤生;王海峰;时洪光;;基于CCD的乳腺X线摄影图像预处理方法的研究[J];应用光学;2010年03期
2 桑梓勤,丁明跃,张天序;雨雾天气下的户外场景成像[J];电子学报;2000年03期
3 刘正军,王长耀,骆成凤;CBERS-1 PSF估计与图像复原[J];遥感学报;2004年03期
4 李奇;冯华君;徐之海;;用于全数字对焦的点扩散函数性能分析与评价[J];浙江大学学报(工学版);2006年06期
5 孙凤国;李月卿;王昌元;刘伟;孙克英;;基于小波变换的提高螺旋CT图像质量的算法[J];泰山医学院学报;2006年03期
6 叶嘉鲁;李福德;肖伟;于枝伶;;含噪声的运动模糊图像点扩散函数尺度鉴别[J];电子器件;2007年02期
7 孙天一;杨盈昀;李楠;;基于POCS算法的块效应消除方法研究[J];中国传媒大学学报(自然科学版);2007年01期
8 邹荣;赵乐军;;基于图像模糊度的雾含水量测量方法[J];计算机工程;2007年14期
9 梅益君;王元庆;;基于点扩散函数的多聚焦图像融合方法[J];计算机工程;2007年19期
10 吕清花;翟中生;;轴锥镜对无衍射成像系统景深的影响[J];光电子技术;2008年01期
11 杨文婧;易鹏;;运动模糊图像复原技术研究[J];大众科技;2009年04期
12 林芳;潘哲朗;;维纳滤波在图像去卷积上的应用[J];科技资讯;2010年03期
13 谢云峰;;运动退化图象的一种复原方法[J];遥测遥控;1985年02期
14 梁路江,付全喜,陶建峰,向正义;一种基于冲击响应的目标识别反卷积算法[J];探测与控制学报;1999年02期
15 朱雯兰,隋成华,唐轶峻,申小阳,胡景耀;退化图像复原技术在天文图片处理中的应用[J];红外;2005年03期
16 刘义鹏,裴锡宇,冯华君,李奇;一种基于DFD的自动对焦算法[J];光学仪器;2005年04期
17 吴友朋;王红霞;周战荣;;无透镜傅里叶变换数字全息图再现像质的影响因素分析[J];红外与激光工程;2006年S4期
18 马宗峰;张春熹;王夏霄;王纪强;;基于1.55μm的全光纤测速激光雷达[J];北京航空航天大学学报;2008年05期
19 翟中生;陈慧;赵斌;;基于轴锥镜大景深成像系统的研究[J];半导体光电;2008年03期
20 周祚峰;水鹏朗;;交替使用小波去噪和全变差正则化的盲图像恢复算法[J];电子与信息学报;2008年12期
中国重要会议论文全文数据库 前10条
1 姜珺;吕帆;金成鹏;毛欣杰;;人眼角膜散光矫正的点扩散函数(PSF)分析表达[A];2006年浙江省眼科学术会议论文集[C];2006年
2 李锐钢;邓伟杰;;采用局部反卷积算法提高数控光学加工的效率[A];中国光学学会2011年学术大会摘要集[C];2011年
3 李勇;李平;文玉梅;;光电成像系统获取图像的多项式能量约束复原[A];中国仪器仪表学会第六届青年学术会议论文集[C];2004年
4 李勇;范承玉;时东锋;冯晓星;乔春红;徐博;;基于加速阻尼Richardson-Lucy算法的湍流退化图像盲复原算法[A];中国光学学会2010年光学大会论文集[C];2010年
5 谢敬辉;孙萍;莫晓丽;;光学扫描全息术系统点扩散函数的优化设计[A];中国光学学会2006年学术大会论文摘要集[C];2006年
6 张旭明;熊有伦;尹周平;;噪声污染运动模糊图像点扩散函数估计[A];图像图形技术研究与应用(2010)[C];2010年
7 张琳琳;姜志国;孟如松;;显微图像复原中的点扩散函数估计[A];第十一届中国体视学与图像分析学术会议论文集[C];2006年
8 刘力;陈文霞;王桂英;;受激发射损耗荧光显微术的系统点扩散函数研究[A];中国光学学会2006年学术大会论文摘要集[C];2006年
9 刘勇;陈家璧;;多模光纤共聚焦系统层析能力的分析[A];第六届全国信息获取与处理学术会议论文集(2)[C];2008年
10 朱宏权;王奎禄;宋顾周;马继明;;LSO闪烁体荧光弥散效应的计算和校正[A];全国计算物理学会第六届年会和学术交流会论文摘要集[C];2007年
中国博士学位论文全文数据库 前10条
1 杨航;图像反卷积算法研究[D];吉林大学;2012年
2 乌大琨;空间无源对抗中的沾染研究[D];哈尔滨工业大学;2009年
3 乌大琨;空间无源对抗中的沾染研究[D];哈尔滨工业大学;2009年
4 孙韶杰;模糊图像中感兴趣信息的盲复原方法研究[D];国防科学技术大学;2010年
5 王忠生;复杂光瞳成像计算机仿真与实验验证技术研究[D];中国科学院研究生院(长春光学精密机械与物理研究所);2011年
6 刘钦晓;波前编码成像系统解码算法理论研究及其应用[D];浙江大学;2010年
7 朱宏权;γ射线针孔照相系统的图像复原研究[D];清华大学;2008年
8 龚绍润;X光同轴相衬成像原理的模拟与实验研究[D];天津大学;2010年
9 马卫红;基于图像分析的光学传递函数测试技术研究[D];中国科学院研究生院(西安光学精密机械研究所);2005年
10 柯才军;微透镜阵列的设计、制作及与CCD的集成技术[D];华中科技大学;2005年
中国硕士学位论文全文数据库 前10条
1 洪鹏飞;基于盲反卷积算法的核能谱获取系统[D];兰州大学;2013年
2 牛延棚;高分辨率扩散加权磁共振成像球面反卷积算法研究[D];浙江工业大学;2013年
3 单晓明;数字图像的盲复原研究[D];南京理工大学;2005年
4 安岗;CCD光学成像系统的点扩散函数及其在亚像素边缘定位中的应用[D];吉林大学;2008年
5 陈涛;卫星遥感图像去薄云薄雾计算机算法研究[D];南京理工大学;2004年
6 孙艳丽;运动模糊图像的恢复与处理[D];大连海事大学;2008年
7 贾昉;数字全息原理及应用的研究[D];西北大学;2008年
8 郝林斌;复色光照明条件下成像系统的点扩散函数[D];山西师范大学;2010年
9 王婷婷;基于多幅图像的有效点扩散函数的高精度测量[D];广东工业大学;2011年
10 刘世豪;基于LO范数和全变差的正则化反卷积算法[D];吉林大学;2013年
中国重要报纸全文数据库 前10条
1 ;红网3D眼镜[N];中国电脑教育报;2001年
2 张京;专家详解数字电视[N];中国消费者报;2000年
3 ;小巧玲珑 品质卓越[N];中国电脑教育报;2002年
4 本报记者 吴洣麓;再见,显像管电视机朋友[N];北京科技报;2008年
5 宣庆庆;请关爱您的电脑[N];中国石化报;2003年
6 新文;DVD、彩电如何搭配?[N];中国文化报;2001年
7 湖北 曾安君;碟机故障检修实例集萃(八)[N];电子报;2006年
8 云龙;让失真、抖动远离IPTV[N];通信产业报;2007年
9 冯华中;显示器:最活跃的桌面舞者[N];电脑商报;2008年
10 黄绍平;PDP真正的平面电视[N];中国电子报;2000年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978