专利名称：单幅图像复原的解析迭代快速频谱外推方法
技术领域：
本发明属于针对空间不变成像系统单幅退化图像复原技术，特别是针对光学遥感单 幅图像复原的解析迭代频谱外推方法。
背景技术：
全天候、全天时高分辨率遥感成像理论和技术在全世界范围内受到极大重视和飞 跃发展，成为空间对地观测发展的热点。未来遥感探测成像器将进一步朝着高分辨率、 高信噪比、轻量化方向发展。遥感成像应用效能(如目标识别)主要取决于获取的遥感 图像的视觉质量和空间分辨率，而空间分辨率是衡量遥感图像的重要指标。另外，为了 获得高分辨率图像，目前大多数高分辨率对地观测相机系统的有效载荷体积和重量庞 大，因此在对地观测和深空探测领域对遥感探测成像器相机的小型化和轻量化提出了非 常高的要求。
理论上，光学遥感图像的高分辨率，依赖于CCD阵列单元的尺寸、CCD像元个数、光 学口径等，但实际获取图像的分辨率还会受到系统光学模糊、运动模糊、系统噪声等多 种误差因素的影响。提高光学遥感图像的清晰度和分辨率通常有两种途径 一是改进和 更新硬件设备；二是通过建立数学和物理模型，利用数据处理技术提高图像的清晰度和 分辨率。
一般遥感成像系统的降质模型与系统点扩散函数/7相关。记F(/J)为/J的傅立叶变 换，并令(《，;7)e[-H为频域坐标，则F(/7)称为调制传递函数(MTF)，它一般包含如下 3部分。
1) CCD的调制传递函数。在焦平面上CCD敏感像素是一个细小的[-c/2,c/2]z的窗
口，其中C为相邻CCD像元的距离，可以用一个矩形函数来描述，这样可推导出CCD的 MTF为:
W "人_ s,c) sin(巧c) 2祸cl^U-2祸咖l
2)光学孔径系统的调制传递函数H。。可看作是一个各向同性的低通滤波器，表示
为味")
—e一2孤c^2"2
3)运动引起的调制传递函数Wm。由于CCD敏感单元会在一定的采样时间内对到
达的光子数目计数。而在采样时间假设系统运动距离为r,则将产生运动模糊，其MTF 为
W匸g 、 ■— sin(吋(^)Kd2》r)
一般而言，遥感成像光学系统可假设为圆对称系统，贝IJ总体MTF是各部分MTF的乘积。 F(/7) = 〃soAV。o〃m。
假设w(;c,力为未降质的图像，从而降质观测图像w。(x,力的成像模型为
"o(x,y) = F-1 (F一y)) 。F(/7)}(xy) +" (1)
或 u0(x，y) = u(x,y)<E>/ (x'y) + n
式中，(x，力为空间坐标，^，7)为频率域坐标，F和F—1分别表示傅立叶正反变换， n表示噪声，0表示巻积算子。
图像复原问题实质上从w。中估计合适的w。 Wiener滤波复原算法源于最小二乘估 计，通过假定图像和噪声都是属于广义平稳随机过程，得到一种简单的频域复原算法，
公
式为:
(2)
、柳o柳+ SJ"W"" 其中F(u)， F(W为t/和^的傅里叶变换，S/"^/W分别为噪声和未失真图像的功
率谱，"。"表示矩阵的元素相乘，而"*"表示复共轭。功率谱比值SJ"W";I起到了规整
化作用。但两个功率谱常常难以估计，因此利用近似Wiener滤波
u厂,f FW'。F J l (3)
式中Y是一个正数，在数值上最好为观测图像信噪比的倒数。
受限制最小二乘复原方法(RLS)，通过Laplace算子生成的循环矩阵的规整化，得 到了如下滤波公式
6(4)
其中c为Laplace算子生成的循环矩阵，WO为c的傅里叶变换。Wiener滤波算法简单、
便于实现、在噪声强度不大(高信噪比)情形能够获得较好的复原效果。但随着噪声强 度的增大，复原效果急剧下降并存在较多的"寄生波纹效应"；RLS具有Wiener滤波类 似的频域公式，实现也简单，能够减少"寄生波纹效应"，但是容易出现过度平滑现象， 不能很好的保持图像边缘。
Lucy-Richardson算法是一种基于泊松分布的最大似然图像复原算法。LR算法可表示
为
(5)
1(a "*)(, ,;)
从以上分析可看出，该算法的超分辨能力来自于每步迭代中涉及的非线性处理，而 且该算法的更新规则自动实现了有约束条件的频谱外推，可用于图像复原。 Lucy-Richardson算法是一种乘积形式的迭代恢复方法，最初主要运用于医学和天文图 像复原。由于噪声的存在，算法虽然能够外推部分高频分量，但同时也破坏了图像通带 内的低频分量，出现了一些寄生波纹，图像复原能力有限，而且收敛速度较慢。
近几年，随着图像处理中的变分偏微分方程的深入研究，出现了一系列新方法。图 像复原属于Hardmard意义下的非适定数学反问题。基于Rudin-Osher连续全变差模型 提出的图像正则化复原算法(l.T. Chan, S. Esedoglu， F. Park, A. Yip, Recent developments in total variation image restoration, in: Handbook of Mathematical Models in Computer Vision, Springer-Verlag， 2005, pp. 17-30.)是目前应用广泛的算法。该算法基于如下变分
mi+。(x,力-"0,力(8)/1(;<:，力||2 (6) 其中连续全变差模型定义为
利用变分法和梯度下降法(GD)转换为求解非线性方程
尝—"。—(7)
7在(7)中AV为勿对应的巻积算子，为避免lv"l等于零，需要引入参数ofl，将其替换为 ^7^7。全变差图像复原模型能够减少"寄生波纹效应"也能较好保持图像边缘，但是
会导致图像中出现"阶梯效应"。另外梯度下降算法存在如下缺点1)收敛速度较慢；2)
在求解过程中需要引入参数f ，而且因为参数e的引入迭代解不能严格收敛于(6)的精 确解；3)每次迭代都要计算一次/Z'H，由于线性空间不变的模糊算子的支集一般很大， 当受到模糊影响的遥感图像尺寸很大时去模糊的计算代价更高。其他变分偏微分方程的 图像复原方法往往存在类似的问题，大大限制了其在工程应用中的实用性。
频谱外推方法也是图像复原的一条有效途径，该方法主要源于解析延拓理论1) 任何空域有界函数的傅立叶变换是解析函数；2)对任何解析函数，只要能准确知道它 在有限区间内的部分信息，就可以唯一的确定整个函数。对一幅图像而言，由于其空域 有界，因此其频谱函数必然解析。根据解析延拓理论，截止频率以上的频谱信息可以通 过截止频率以下的频谱得以重建，从而能够实现图像的高分辨率复原。例如根据此思想 提出的一种针对无源毫米波成像运用的最大似然频域校正超分辨算法。该算法以 Liicy-Richardson算法作为主迭代，通过一种频域校正算法，用Wiener滤波器恢复的频谱 代替通带内的频谱，在外推图像高频分量的同时，保证图像的低频分量不被破坏，对无 源毫米波成像复原获得了应用(2.郑鑫，杨建宇，李良超.无源毫米波成像的最大似然频 域校正超分辨算法.自动化学报，2009， Vol 35,N0.1,pp 28-33)。但是由于Wiener滤波和 Lucy-Richardson算法本身存在"寄生波纹效应"明显和边缘清晰度不足等问题，直接将此 两种算法耦合外推对与复杂的遥感图像效果处理并不理想。
发明内容
本发明的目的在于提供一种单幅图像复原的解析迭代快速频谱外推方法，以具有解 析计算的迭代频谱外推方式，保证图像通带内频谱和高频频谱有效合成，克服传统方法 "寄生波纹效应"明显和边缘清晰度不足等问题，而且算法充分利用快速傅立叶变换技 术，实现图像的快速复原。
实现本发明目的的技术解决方案为 一种单幅图像复原的解析迭代快速频谱外推方 法，包括迭代系统初始化过程和解析迭代过程，即
1.1所述的迭代系统初始化过程包括输入一幅待复原的MxA/大小的降质图像W。，计算降质图像频谱/Yu";设置初始迭代解17(0)=17。，初始频谱校正图像w(w二l/。；设置
迭代系统参数初始值梯度信息阈值"">>0 ，预测贡献度控制参数"^>0，合成贡
献度参数""> 1 ，图像频谱误差预设阈值e《10—4 ;并设置空间不变成像系统的系统
调制传递函数WW;
1.2所述的解析迭代过程包括如下步骤
步骤l:首先利用第k-l次(k21)的频谱校正图像M^'、分别计算其在水平方向和 垂直方向的向前差分，并利用第k-l次的梯度信息阈值""—〃作阈值收縮判断，得到
第k次的水平方向差分估计o^ 和垂直方向差分估计o^^;
步骤2:将第k次的水平方向差分估计(v,)^和垂直方向差分估计(M「/作快速傅立
叶变换(FFT)得到水平方向高频频谱F(v(")和垂直方向高频频谱尸(v〖")，并与第 k-l次合成图像频谱F ("'O进行频谱预测得到第k次的预测频谱iTM^入
步骤3:利用图像频谱合成校正，将第k次的预测频谱F(^^和降质图像频谱F "结合光学系统调制传递函数F(W及其复共轭F(W'通过矩阵元素逐点相乘、相加和矩阵元 素逐点相除的方式进行合成校正，得到合成校正频谱fY""。；
步骤4:采用2范数的平方l.g计算第k次校正图像频谱W"f力和第K-l次校正图像
频谱w""-"的误差，如果Ifof^—'。lMi""w《"，则转至步骤5，否则转至步骤
6;
步骤5:表明频谱外推没有收敛，则进行参数自适应校正，转至步骤l, k=k+l; 步骤6:判断是否达到收敛条件|^",-/^*-"|/|/^""《"，如果到达则算法结
束，输出复原图像，对合成校正图像/f^^进行快速逆傅立叶变换得到复原图像。
本发明与现有技术相比，其显著优点(1)本发明能够实现高视觉质量的快速复 原，能够很好的克服传统Wiener滤波算法的"寄生波纹效应"和约束最小二乘算法的"过 度平滑现象"，边缘清晰度高。经大量统计实验结果表明，而本发明方法一般在1.56秒时间内，能够取得峰值信噪比(PSNR)比传统算法高1.0-2.0 dB左右，改进信噪比高
0.5-0.8dB左右，在抗噪性能、去模糊效果、实时性等三方面得到很好的平衡。(2)本
发明方法在光学遥感成像系统、合成孔径雷达成像系统、电视制导、航空探测等都有广
泛的应用前景，同时也为超分辨复原提供了新的途径。 下面结合附图对本发明作进一步详细描述。
图1是本发明的整体迭代结构的流程图。 图2是本发明解析迭代过程步骤1数据流程图。 图3是本发明解析迭代过程步骤2数据流程图。 图4是本发明解析迭代过程步骤3数据流程图。 图5是本发明的实验测试图像和MTF曲线。 图6是本发明与传统方法的一个复原实验对比效果图。 图7是本发明与传统方法的一个复原实验中的边缘清晰度对比图。
具体实施例方式
本发明以成像系统光学调制传递函数、运动模糊函数为基础，通过不断精细的估计 梯度信息、预测频谱和合成低高频频谱的方式进行频谱迭代外推，达到降质图像反演清 晰图像的目的，基本思想源于解析延拓理论1)任何空域有界函数的傅立叶变换是解 析函数；2)对任何解析函数，只要能准确知道它在有限区间内的部分信息，就可以唯 一的确定整个函数。对一幅图像而言，由于其空域有界，因此其频谱函数必然解析。根 据解析延拓理论，截止频率以上的频谱信息可以通过截止频率以下的频谱得以重建，从 而能够实现图像的高分辨率复原。
结合图1,为实现空间不变成像系统降质图像高视觉质量的快速图像复原，本发明 基于解析迭代频谱外推的方法包括迭代系统初始化过程和解析迭代过程。
1.1所述的迭代系统初始化过程包括①输入一幅待复原的MxA/大小的降质图像
W。，计算降质图像频谱/W;②初始迭代解i/W-U。，初始频谱校正图像WW-U。；③ 设置迭代系统参数初始化值梯度信息阈值"(">0,预测贡献度控制参数々^>0， 合成贡献度参数,"> 1 ，图像频谱误差预设阈值e ^ 10—4 ; 空间不变成像系统的系统调制传递函数&W ，其中FfW为各部分MTF的乘积F(/7卜〃sOW。0 。 1.2所述的解析迭代过程步骤如下
步骤l:首先利用第k-l次的频谱预测图像H^W (k21)，分别计算其在水平方向和
垂直方向的向前差分，并利用第k-l次的梯度信息阈值"^-"作阈值收縮判断，得到第 k次的水平方向和垂直方向信息估计。
如图2所示该步骤的详细过程为，首先通过具有周期边界条件的差分算子计算第k-
1次的频谱校正图像W-^在空间位置(/j')的水平和垂直方向差分当""W,"^M时， 计算第k- 1次频谱预测图像在空间位置(/J)处的水平方向的差分 (A""l =《了)-《"和垂直方向的差分(/),(")L =《 > -;在, = 7V，y = M时的图像边 界需采取周期边界条件计算水平方向差分(/>,("^-H^，和垂直方向差分
然后计算第k-l次的频谱校正图像在空间位置(/,"的梯度模 l(Dw("')),』=^((iX'l.)' + ,丫 ，
最后利用第k-l次的梯度信息阈值《("')进行梯度阈值收縮得到水平方向和垂直方向
的差分信息估计，得到第k次的水平方向差分估计o^ 和垂直方向差分估计(o;,4/。估
计方法为
水平方向
垂直方向h)
步骤2:利用髙频频谱预测得到第k次的预测频谱。如图3所示该步骤的详细过程
为将第k次的水平方向差分估计o^ 和垂直方向差分估计og^作快速傅立叶变换
(FFT)得到水平方向高频频谱F(v"0和垂直方向高频频谱F("/0 ;将水平方向差
11分算子[-l,l]巻积核/),和垂直方向差分算子
-1, 1
巻积核A作FFT;结合第k-l次复原图
像频谱F"("",通过矩阵元素逐点相加、相乘和矩阵元素逐点相除的方式作预测更
新得到第k次的预测频谱F(V^ 。计算公式为
其中F(A)， F(Z^分别为水平和垂直差分算子线性巻积核的快速傅立叶变换； 为正参数，称为预测贡献度控制参数，控制水平方向高频频谱和垂直方向高频频谱对预 测频谱的贡献度。"。"表示矩阵元素逐点相乘，"/"表示矩阵元素逐点相除。
步骤3:利用图像频谱合成校正。如图4所示，该步骤的详细过程为将第k次的
预测频谱尸6W和降质图像频谱/^M结合光学系统调制传递函数/YW及其复共轭F"/
通过通过矩阵元素逐点相乘、相加和矩阵元素逐点相除的方式进行合成校正，得到合成 校正频谱W""M。计算公式为
其中斤W为光学系统调制传递函数，其可以包括CCD的调制传递函数、光学孔 径系统的调制传递函数、运动引起的调制传递函数等多个传递函数的乘积形成的系统调 制传递函数，为FfW的复共轭。""为取正值的参数，称为合成贡献度参数。同 样"。"表示矩阵元素逐点相乘，"/"表示矩阵元素逐点相除。
步骤4:采用2范数的平方ll计算第k次校正图像频谱F(V"和第K-l次校正图像 频谱F(V"";的误差，如果||fY"，-F(V*-"|/||/=V"《>e，则转至步骤5，否则转至步骤 6。
步骤5:表明频谱外推没有收敛，则进行参数自适应校正，自适应调整参数，转至 步骤l, k=k+l。参数调整方法如下
121. 梯度信息阈值0^-O ，在迭代过程中逐渐减少，本发明方法中经实验证实"。取 特定的参数值"《"W ^25，按照特定的递减顺序W" < /进行控制即满足 濟>.... >0即可。
2. 本发明方法经研究实验测得贡献度参数;^W与需与W^满足正比关系， / W=A:'a ，A:e(^U,0.S,为实验参数。
3. 合成贡献度参数,在迭代过程中按照递增顺序进行控制 Y("/""') > 2且Y(0)e "，20/ 。
步骤6:判断是否达到收敛条件kOfY""-'《/||F(V"《"，如果到达则算法结 束，输出复原图像，对合成校正图像/f^O进行快速逆傅立叶变换得到复原图像；
步骤l、 2、 3是本发明的核心处理过程，均参与迭代更新过程。每个核心单元均具 有解析形式的数学公式控制。步骤1在图像空域中进行，而步骤2和步骤3均在傅立叶 域进行，整个频谱一次外推过程包含二次快速傅立叶正变换和一次快速傅立叶逆变换的
算法复杂度。算法的迭代执行速度快、对一幅MxM大小的图像整体复杂度为0^/'/。gi^。
下面结合图5、图6和图7所示，通过一个图像复原实施例及其效果评价来说明本
发明方法的技术效果。
在实施例中，如图S.a所示，实验图像为一幅清晰的370x370的光学遥感图像(参 考图像)。如图5.b所示为实验中遥感成像系统的MTF模拟曲线。
本发明方法中的迭代参数"^=700，々^=30,/^=^，迭代终止条件均取为 |一""0-7=7 |)/|||^""0|〖《10-4 。
实验中对比算法包括MATLAB 7.1中Wiener滤波、约束最小二乘算法(RLS)， Lucy-Richardson算法；同时我们和全变差梯度最速下降法(TVGD)进行比较。TVGD采 取中心差分的离散迭代格式7/ " — 力+ a晰-
其中参数y取为O.l，网格比r取W， Lagrange参数a取25.0。
实施例在MATLAB 7.1平台模拟仿真实现，计算环境为Dell D620笔记本电脑、Intel 1.83GHz CPU，1G内存。为了验证本发明方法的有效性。用于定量评估去模糊算法的性
能指标采取改进信噪比(A，)、峰值信噪比(PSNR)和相对误差(i e五厅)。设《。， 分别为wxw的降质图像、参考图像和去模糊图像，A，、 PSNR和i e五rr定义如下
厶咖=70/吸7。{|"—"。|;/|"—"I;}，
尸SiVR = 7。/收/。卜x W x 2552/||" - m*『},
一般而言改进信噪比A，和峰值信噪PSNR比越大越好，而相对误差ife￡>r越小越好。
图6和表1分别给出了不同算法和本发明方法的复原结果图像和各项性能指标。
如图6.a所示为根据降质模型模拟生成的降质图像，其中加入的噪声为均值为0、 方差为15高斯噪声。可看到降质图像模糊效应和噪声干扰非常严重，导致空间分辨率 非常低，严重影响人们的视觉理解。如图6.b所示为Wiener滤波结果。Wiener滤波复 原图像"寄生波纹效应"明显，复原效果较差。如图6x所示为约束最小二乘(RLS)算法的 滤波结果，如图6.d所示为Lucy-Richardson算法的滤波结果，可看到两种算法对于"寄 生波纹效应"有所改善，但"阶跃"边缘和线状边缘保持效果不好(例如图6.c，图6.d中江 面和城市的交界线)。如图6.e所示为TVGD算法结果，该算法能够极大的减少"寄生 波纹效应"，视觉质量得到较好提高，但存在视觉上的"阶梯效应"(如图6.d)。如图6.f 所示为本发明方法的复原图像，本发明方法复原图像很好的克服了"寄生波纹效应"，并 且复原图像的细节在对比算法中最清晰，视觉质量最高。
14表一给出了本文发明方法和对比算法的性能指标情况。从复原图像质量客观评价指 标来看，本发明方法复原图像质量的客观评价指数得到进一步改善。例如PSNR比TVAP
提高1dB,而a卿比TVAP提高0.8dB。
从算法所花费的CPU时间来看，Wiener滤波和RLS算法是非迭代的算法，时间最 少，分别为0.14秒和0.31秒。Lucy-Richardson算法、TVGD算法和本发明方法都属 于迭代算法，时间较长。其中TVGD算法的所用时间最长，达到84.281秒，说明算法 收敛速度较慢；Lucy-Richardson算法花费3.08秒，而本发明方法仅用了 1.56秒就达到 优于TVGD算法的图像质量。
表1:不同复原算法性能指标比较结果
MethodPSNRReErrT i me (s)
降质图像13.27160.0097841.97-
Wiener滤波20.64310.0081042.410.14
RLS算法21.18970.0064563.780.31
Lucy-Richardson算法20.24620.0073143.233.08
TVGD22.00170.0061243.8184.281
本发明方法23.41000.0051254.841.56
通过对复原图像边缘细节点集比对情况能够很好的区分复原算法的有效性。如图7 所示为通过Matlab中的Canny算子提取原始图像、降质图像和各种算法输出图像的边缘 图进行分析。图7.&为参考图像的清晰边缘图，图7.b为降质图像的边缘图。图7.c为Wiener 滤波图像边缘图，不难发现Wiener滤波图像边缘存在较多的"寄生波纹效应"引起的伪 边缘；图7.d所示为RLS算法复原图像边缘，图7.e所示为Lucy-Richardson算法复原图像 边缘，相比Wiener滤波算法，这两种算法伪边缘得到减少，但是也丧失了较多的边缘细 节。图7.f为全变差梯度最速下降法(TVGD)复原图像边缘图，该算法边缘清晰度得到很 大提高。图7.g为本发明方法的复原图像边缘图，进行综合对比，本发明方法对图像边缘 的保持效果最好，细节得到很好的复原。因此本发明方法得到了很好的图像质量-计算 时间性价比，实用性非常强。
权利要求
1、一种单幅图像复原的解析迭代快速频谱外推方法，其特征在于包括迭代系统初始化过程和解析迭代过程，即1.1所述的迭代系统初始化过程包括输入一幅待复原的M×N大小的降质图像u0，计算降质图像频谱F(u0)；设置初始迭代解u(0)＝u0，初始频谱校正图像w(0)＝u0；设置迭代系统参数初始值梯度信息阈值α(0)＞0，预测贡献度控制参数β(0)＞0，合成贡献度参数γ(0)＞1，图像频谱误差预设阈值ε≤10-4；并设置空间不变成像系统的系统调制传递函数F(h)；1.2所述的解析迭代过程包括如下步骤步骤1首先利用第k-1次(k≥1)的频谱校正图像w(k-1)，分别计算其在水平方向和垂直方向的向前差分，并利用第k-1次的梯度信息阈值α(k-1)作阈值收缩判断，得到第k次的水平方向差分估计(v1)i，j(k)和垂直方向差分估计(v2)i，j(k)；步骤2将第k次的水平方向差分估计(v1)i，j(k)和垂直方向差分估计(v2)i，j(k)作快速傅立叶变换(FFT)得到水平方向高频频谱F(v1(k))和垂直方向高频频谱F(v2(k))和第k-1次合成图像频谱F(u(k-1))进行频谱预测得到第k次的预测频谱F(w(k))；步骤3利用图像频谱合成校正，将第k次的预测频谱F(wk)和降质图像频谱F(u0)结合光学系统调制传递函数F(h)及其复共轭F(h)＊通过矩阵元素逐点相乘、相加和矩阵元素逐点相除的方式进行合成校正，得到合成校正频谱F(u(k))；步骤4采用2范数的平方||·||22计算第k次校正图像频谱F(u(k))和第K-1次校正图像频谱F(u(k-1))的误差，如果<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><msubsup> <mrow><mo>|</mo><mo>|</mo><mi>F</mi><mrow> <mo>(</mo> <msup><mi>u</mi><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </msup> <mo>)</mo></mrow><mo>-</mo><mi>F</mi><mrow> <mo>(</mo> <msup><mi>u</mi><mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo></mrow> </msup> <mo>)</mo></mrow><mo>|</mo><mo>|</mo> </mrow> <mn>2</mn> <mn>2</mn></msubsup><mo>/</mo><msubsup> <mrow><mo>|</mo><mo>|</mo><mi>F</mi><mrow> <mo>(</mo> <msup><mi>u</mi><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </msup> <mo>)</mo></mrow><mo>|</mo><mo>|</mo> </mrow> <mn>2</mn> <mn>2</mn></msubsup><mo>></mo><mi>ϵ</mi><mo>,</mo> </mrow>]]></math> id="icf0001" file="A2009100312700002C1.tif" wi="54" he="5" top= "217" left = "72" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>则转至步骤5，否则转至步骤6；步骤5表明频谱外推没有收敛，则进行参数自适应校正，调整系统参数，转至步骤1，k＝k+1；步骤6判断是否达到收敛条件<maths id="math0002" num="0002" ><math><![CDATA[ <mrow><msubsup> <mrow><mo>|</mo><mo>|</mo><mi>F</mi><mrow> <mo>(</mo> <msup><mi>u</mi><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </msup> <mo>)</mo></mrow><mo>-</mo><mi>F</mi><mrow> <mo>(</mo> <msup><mi>u</mi><mrow> <mo>(</mo> <mi>k</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo></mrow> </msup> <mo>)</mo></mrow><mo>|</mo><mo>|</mo> </mrow> <mn>2</mn> <mn>2</mn></msubsup><mo>/</mo><msubsup> <mrow><mo>|</mo><mo>|</mo><mi>F</mi><mrow> <mo>(</mo> <msup><mi>u</mi><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </msup> <mo>)</mo></mrow><mo>|</mo><mo>|</mo> </mrow> <mn>2</mn> <mn>2</mn></msubsup><mo>≤</mo><mi>ϵ</mi><mo>,</mo> </mrow>]]></math> id="icf0002" file="A2009100312700003C1.tif" wi="54" he="5" top= "31" left = "92" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>如果到达则算法结束，输出复原图像，对合成校正图像F(u(k))进行快速逆傅立叶变换得到复原图像。
2、根据权利要求1所述的单幅图像复原的解析迭代快速频谱外推方法，其特征在 于解析迭代过程步骤1中梯度信息阈值收縮满足如下关系水平方向b)〔y=附似(jl( …"),』-"("')，"垂直方向=挑《^ (Z)w"-〃). J-<formula>formula see original document page 3</formula>其中梯度算子Z)的采取向前或向后差分算子，如:A/,这里, = WJ = M时采取周期边界条件。
3、根据权利要求1所述的单幅图像复原的解析迭代快速频谱外推方法，其特征在 于解析迭代过程的步骤2中利用高频频谱预测得到第k次的预测频谱，是将第k次的水 平方向差分估计0^ 和垂直方向差分估计hf作快速傅立叶变换(FFT)得到水平方向高频频谱,(^")和垂直方向高频频谱,(^"):将水平方向差分算子[-l，l]巻积核A和垂直方向差分算子-1, 1巻积核巧作FFT;并将第k-l次复原图像频谱, ("))通过矩阵元素逐点相加、相乘和矩阵元素逐点相除的方式作预测更新得到第k次的预测频 谱F(V^)，计算关系为
4、根据权利要求1所述的单幅图像复原的解析迭代快速频谱外推方法，其特征在 于解析迭代过程的步骤3中利用图像频谱合成校正，是将第k次的预测频谱/Tv^和降 质图像频谱/^M结合光学系统调制传递函数FfW及其复共轭F^/通过矩阵元素逐点相乘、相加和矩阵元素逐点相除的方式进行合成校正，得到合成校正频谱W""。，计算关 系为
5、 根据权利要求1所述的单幅图像复原的解析迭代快速频谱外推方法，其特征在于迭代系统包含梯度信息阈值"^—"、预测贡献度参数AW—〃 、合成贡献度参数,"" 三个关键参数。
6、 根据权利要求5所述的单幅图像复原的解析迭代快速频谱外推方法，其特征是梯度信息阈值参数0^—M乍如下调整梯度信息阈值"W^，在迭代过程中逐渐减少，且 "。取特定的参数值"《"W^725，按照特定的递减顺序"W/^T"/^7进行控制即满足> >…""-">....> 0即可。
7、 根据权利要求5所述的单幅图像复原的解析迭代快速频谱外推方法，其特征是预测贡献度参数作如下调整贡献度参数々W;与需与"满足正比关系， y^"=A:'aW， A:e(^.二0.S」为实验参数。
8、根据权利要求5所述的单幅图像复原的解析迭代快速频谱外推方法，其特征是合 成贡献度参数作如下调整在迭代过程中按照递增顺序进行控制 2且Y(。)e /1，20/ 。
全文摘要
本发明公开了一种单幅图像复原的解析迭代快速频谱外推方法，包括迭代系统初始化过程和解析迭代过程。其中初始化过程计算降质图像频谱、设置初始迭代解、系统初始参数；解析迭代过程通过不断精细的梯度信息阈值收缩估计图像的水平和垂直方向梯度信息，通过解析迭代频谱和梯度信息频谱解析预测得到预测频谱，通过图像频谱合成校正将降质图像频谱和预测频谱解析合成校正得到高分辨率图像。本发明方法基于快速傅立叶变换技术，通过截止频率以下频谱，外推截止频率以上频谱，以很小的复杂度达到图像快速复原。去模糊和抑制噪声能力明显优于常规方法，信噪比得到显著提高。该方法非常便于利用快速傅立叶变换DSP芯片组成的图像处理硬件实现。
文档编号G06T5/00GK101540043SQ200910031270
公开日2009年9月23日 申请日期2009年4月30日 优先权日2009年4月30日
发明者吴泽彬, 亮 肖, 韦志辉, 黄丽丽 申请人:南京理工大学