专利名称：结合帧选择和盲解卷积的自适应光学图像高分辨率复原方法
技术领域：
本发明属于图像处理技术，提出了一种针对自适应光学系统部分校正图像的盲图像复 原技术。
背景技术：
图像复原就是从观测到的降质图像中得到真实图像的估计过程。其成像模型可以表示为 gO，_y) = /0,;v) /i(x，y) + "(>，:y) (1)
式中，g(x，力表示观测到的降质图像；/(x,力代表真实目标，即待复原目标；/ (;c，力表 示光学系统的点扩散函数，用来描述波前相差；"(x,力表示系统噪声；"0"代表循环巻 积。自适应光学系统由于受限于变形镜自由度不够以及波前传感器测量不准等因素，往往 不能对对降质图像完全校正。因此，需要对其部分校正的图像结果进行二次复原。波前传 感器受噪声等因素制约，由其探测的波前信息得到的点扩散函数不能完全表示实际的波前 扰动，需要对其进行修正。
自适应光学系统闭环过程中的短曝光图像序列中并非每一帧作二次校正都能取得满意 的结果，帧选择技术可以从序列中选择出最适合二次校正的帧组合，使用这样的帧进行盲 解巻积处理可以有效地提高复原质量，提高盲解巻积的收敛速度。
梁莹等曾提出过利用盲解巻积对部分校正的自适应光学图像进行事后处理，但未考虑到并 非所有记录的降质图像都可以通过解巻积得到图像质量上的提高；J. J. Green等曾提出使用帧选 择技术对短曝光序列进行事后处理，但他们只是利用波前传感器中的数据进行解巻积，而未考 虑到波前传感器本身受噪声影响并不能完全代表波前相差；他们也未考虑对自适应光学系统部 分校正后的图像进行二次校正。
发明内容
本发明的目的是克服现有技术的不足，提出一种结合帧选择和盲解巻积的自适应光 学图像高分辨率复原方法，该方法运算速度较快、收敛性有保证。使用此方法可以有效地 提高待校正图像的成像质量。
本发明的技术解决方案是利用帧选择与盲解巻积的图像复原算法，其特点在于步骤 如下
(1) 记录自适应光学闭环校正时的短曝光图像序列，计算序列中每一帧图像的香农熵, 选择熵较小的降质图像&(X，力进行盲解巻积图像复原；
(2) 使用随机相位生成点扩散函数的初始值&(X,力；
(3) 使用&(JC,力与得到的&(JC,30估计出目标六x，力，对其加入正性限制，得到估计
(4) 使用^(x，y)与/(x，力得到出点扩散函数的估计值&Oc,力，同样加入正性限制， 得到估计值《(JC,力；
(5) 检查迭代值S(;c，力、/(x,;/)是否达到迭代停止要求，若未满足则重回步骤(3); 否则停止循环，输出/(x,力和《(x,力。
本发明与现有技术相比有如下优点
(1) 帧选择从众多降质图像中选择出部分最适合二次处理的降质图像，以避免不合格
的降质图像对真实目标复原的影响；并且能够提高盲解巻积的收敛速度，减小其陷入无效 解的可能。
(2) 使用同一目标的多帧降质图像可以有效地提高复原质量，帧数越多复原效果越好。 同时，多帧降质图像可以很好地避免迭代过程中出现无效解。对较强的波前扰动，只要适 当增加短曝光降质图像的帧数就能得到较好的复原结果。
(3) 数据依赖仅仅为降质图像，越少的数据依赖，则有越强的适应性和普适性。
图l为本发明的算法流程图2是利用本发明复原受随机扰动影响的文字图像；其中，(a)为源图像，(b) (c)图 为利用帧选择技术从IOO帧加入随机扰动的含噪降质图像中遴选出的两帧，(d)图为复原图像；
图3为使用本发明对实测恒星目标的复原处理，其中，(a) (c)为经过帧选择的三 帧100像素X100像素自适应光学校正后恒星目标。(d)为经本发明方法二次校正后的结果；
图4为没有帧选择过程的实测恒星目标的复原处理，其中，(a) (c)为随机选出的三 帧降质图像；(d)为盲解巻积后的结果；
图5为图3 (a)、图3 (d)和图4 (d)的截面图。
具体实施例方式
下面对本发明做详细说明。针对现在图像复原算法的不足，本发明提出了一种利用帧
选择与盲解巻积的图像复原算法即从记录的自适应光学系统闭环时的校正图像短曝光序 列中遴选出最适合作二次处理的降质图像进行盲解巻积事后处理。使用同一目标的多帧降
质图像而非单帧可以将一个"知一求二"(已知g(x，力未知/j(x，力与/0,力)的问题转变为 "知m求m+l"(已知gmO，力未知&(x，力与/(;c,力)的问题，大大减少了问题的不确定
性；同时引入了帧选择剔除了不合适的降质图像对真实目标进行估计，增加了迭代精度， 使整个算法更容易收敛。该算法的流程图由图1所示。下面分别说明
(1) 记录W帧连续的自适应光学系统闭环校正时的短曝光降质图像序列^Oc，力。
(2) 计算^(x,力中每一单帧的香农熵，
M =-&'0，力108^(",力' (2)
其中^,(x,力代表第/帧的概率密度函数，定义为-
小力-^^V (3)
(3) 根据所需进行解巻积的帧数S，数选择M,最小的S帧形成降质图像序列gJ;c,力，
(4) 随机产生的波前残差《0c,力，使用点扩散函数的物理计算公式(4)得到其对应 的点扩散函数作为初始值-
^iT(Pe邻[见"力])12 (4) 式中，户代表瞳径函数，iT表示傅里叶变换。
(5) 使用g力,力与得到的&"力估计出目标/(x，力。定义&a,y;)、 &(/;，/;)和 户(X,/》分别为^"力、乞(x,力和的傅里叶频谱。
2X'(力，力X^o;,力)
!>m(y;,y;)|
六x,力"一c/;,/;》 (6) 其中，"*"代表复共轭；/Fr代表逆傅里叶变换。s,为一常量防止除数为0。
(6) 根据光成像的物理意义，/(x，力20。因此加入正性限制。
D,，M>0 (7) L 0，otherwise
(7) 根据维纳逆滤波原理，使用^(xj)与/(x,力得到出点扩散函数的估计值&(;c,力。
戶c/;,y;)和之o;,力)分别是i""y)和&"力的傅里叶变换；^的作用与&—样。
(8) 根据点扩散函数的物理计算公式(4)，对4,(;c，力同样加入正性限制。
<formula>complex formula see original document page 7</formula>
(9) 定义迭代停止条件。如果未达到条件要求，则重回(3)步，直至收敛为止。
<formula>complex formula see original document page 7</formula>
H代表定义在希尔柏特空间上的范数；e为一正常数。 实施例1
如图2所示，利用本发明对模糊字体进行复原处理。其中，(a)为源图像，(b) (c)图 为利用帧选择技术从100帧加入随机扰动的含噪降质图像中遴选出的两帧，(d)图为复原图像； 结果显示，本发明可以有效地复原降质图像。整个计算过程经过25次迭代，耗时20秒。 实施例2
应用本发明对经过自适应光学系统部分校正的图像进行复原，如图3所示。(a) ~ (c) 为经过帧选择的三帧100像素X100像素自适应光学校正后恒星目标。(d)为经本发明方法二 次校正后的结果。实验表明，复原结果较好。整个计算过程耗时30秒，经过了30次迭代。从 图3可以看到，在相同能量下，(d)的峰值较(a) (c)分别提高了2.56倍、2.38和2.45倍， 也就意味着斯特列尔比分别提高了2.56倍、2.38禾tl 2.45倍。图4比较了图3 (a)与(d)的截 面图。
为作比较，从降质图像短曝光序列中随机选出三帧图像进行二次处理，如图4所示(a)
(c) 为随机选出的三帧降质图像；(d)为盲解巻积后的结果。可以发现，在相同能量下，图4
(d) 的峰值只有图3 (d)的65.P/。。但经二次校正后，图4 (d)的斯特列尔比较(a)与(c) 分别提高了 1.28倍与1.14倍。
图5分析了二次校正在分辨率上的提高。可以看到，图3 (d)与图4(d)的分辨率都较图3(a) 都有明显的提高，分别为4.9个象素与5.1个象素，其中前者达到了系统衍射极限，后者也十分 拉近。
综上所述，本发明可以有效地提高成像质量，可以对自适应光学系统进行很好补充。
权利要求
1、一种结合帧选择和盲解卷积的自适应光学图像高分辨率复原方法，其特征在于(1)记录自适应光学闭环校正时的短曝光图像序列，计算序列中每一帧图像的香农熵，选择熵较小的降质图像gm(x，y)进行盲解卷积图像复原； (2)使用随机相位生成点扩散函数的初始值(3)使用gm(x，y)与得到的估计出目标，对其加入正性限制，得到估计值(4)使用gm(x，y)与得到出点扩散函数的估计值，同样加入正性限制，得到估计值(5)检查迭代值是否达到迭代停止要求，若未满足则重回步骤(3)；否则停止循环，输出和
2、 根据权利要求1所述的一种结合帧选择和盲解巻积的自适应光学图像高分辨率复原 方法，其特征在于所述步骤(I)中的短曝光图像中的曝光时间范围在5ms 40tns间。
3、 根据权利要求1所述的一种结合帧选择和盲解巻积的自适应光学图像高分辨率复原 方法，其特征在于所述步骤(1)中计算^(X,力中每一单帧的香农熵公式为<formula>see original document page 2</formula>其中y,O，力代表第Z'帧的概率密度函数，定义为<formula>see original document page 2</formula>
4、 根据权利要求1所述的一种结合帧选择和盲解巻积的自适应光学图像高分辨率复原 方法，其特征在于所述步骤(2)中点扩散函数作为初始值公式为<formula>see original document page 2</formula>式中，P代表瞳径函数，iT表示傅里叶变换。
5、 根据权利要求1所述的一种结合帧选择和盲解巻积的自适应光学图像高分辨率复原 方法，其特征在于所述步骤(3)中使用gm(x,jO与得到的&(;c,力估计出目标/Oc,力的公式为定义gm(fx,fy)、Hm(fx;,fy)和F(fx，fy)分别为gm(x,y)，力、hm(x,y)和y(x,y)力的傅里叶频 谱，<formula>see original document page 3</formula> (5)<formula>see original document page 3</formula> (6) 其中，"*"代表复共轭,IFT代表逆傅里叶变换，ε1为一常量防止除数为o。
6、根据权利要求i所述的一种结合帧选择和盲解巻积的自适应光学图像高分辨率复原方法，其特征在于所述步骤(3)中加入正性限制的公式为<formula>see original document page 3</formula>。
7、 根据权利要求1所述的一种结合帧选择和盲解巻积的自适应光学图像高分辨率复原方法，其特征在于所述步骤(4)中使用gm(x,y)与f(x,y)得到出点扩散函数的估计值hm(x,y)的公式为:<formula>see original document page 3</formula><formula>see original document page 3</formula>F( x, y)和Hm(fx,fy)分别是f(x,y)和hm(x,y)的傅里叶变换，ε2为一常量防止除数为o。
8、 根据权利要求1所述的一种结合帧选择和盲解巻积的自适应光学图像高分辨率复原方法，其特征在于所述步骤(4)中加入正性限制的公式为<formula>see original document page 3</formula>。
9、 根据权利要求1所述的一种结合帧选择和盲解巻积的自适应光学图像高分辨率复原方法，其特征在于所述步骤(5)中的停止条件为<formula>see original document page 3</formula>‖:‖表定义在希尔柏特空间上的范数，ε为一正常数。
全文摘要
一种结合帧选择和盲解卷积的自适应光学图像高分辨率复原方法(1)记录自适应光学闭环校正时的短曝光图像序列g<sub>n</sub>(x，y)，计算序列中每一帧图像的香农熵，选择熵较小的降质图像g<sub>m</sub>(x，y)进行盲解卷积图像复原；(2)使用随机相位生成点扩散函数的初始值h<sub>m</sub>(x，y)；(3)使用g<sub>m</sub>(x，y)与得到的h<sub>m</sub>(x，y)估计出目标f(x，y)，对其加入正性限制，得到估计值f(x，y)；(4)使用g<sub>m</sub>(x，y)与f(x，y)得到出点扩散函数的估计值h<sub>m</sub>(x，y)，同样加入正性限制，得到估计值h(x，y)；(5)检查迭代值h(x，y)、f(x，y)是否达到迭代停止要求，若未满足则重回步骤(3)；否则停止循环，输出f(x，y)和h(x，y)。本发明有效地提高了复原质量，加快了收敛速度，可以很好地补充自适应光学系统受硬件限制的校正能力，提升成像质量。
文档编号G06T5/00GK101206762SQ200710177560
公开日2008年6月25日 申请日期2007年11月16日 优先权日2007年11月16日
发明者雨 田, 饶长辉 申请人:中国科学院光电技术研究所