发明名称 图像超分辨率与像质增强方法
摘要 图像超分辨率与像质增强方法,涉及空间遥感成像系统中的数字图像处理技术领域,解决了现有维纳滤波图像复原法存在的可被利用优化元素少、像质提升空间有限而无法有效提高图像质量的问题,该方法如下:利用光线分束器在两个通道分别采集在焦像和离焦像,并采集N幅相对于在焦像在水平和竖直方向上均错位1/N个像元的错位在焦像;对在焦像和离焦像进行PD处理,采用共轭梯度算法求出像差系数,求出经PD处理复原后目标图像并计算目标函数,对错位在焦像和复原后目标图像进行维纳滤波处理,对得到两幅复原图像再进行亚像元融合得到一幅具有超高分辨率的复原图像。本发明实现了图像的超分辨率和高清晰度,分辨率提高了1.4倍。
申请公布号 CN103310427B 申请公布日期 2015.11.18
申请号 CN201310253960.4 申请日期 2013.06.24
申请人 中国科学院长春光学精密机械与物理研究所 发明人 王刚;王磊;许博谦
分类号 G06T5/10(2006.01)I 主分类号 G06T5/10(2006.01)I
代理机构 长春菁华专利商标代理事务所 22210 代理人 王丹阳
主权项 图像超分辨率与像质增强方法,其特征在于,该方法的条件和步骤如下:步骤一、在光学成像系统中,扩展目标的光线经过望远物镜入射至光线分束器,分成两束分别在两个CCD接收图像,第一CCD在焦面位置上接收一幅在焦像,第二CCD在一个已知离焦量的离焦位置上接收一幅离焦像;使光线分束器相对于第一CCD的焦平面阵列方向进行扫描,采集N幅在第一CCD上形成的相对于第一CCD上的在焦像在水平和竖直方向上均错位1/N个像元的错位在焦像,N的取值范围是:2≤N≤10;步骤二、将步骤一采集的一对在焦像和离焦像传送至DSP图像处理器,分别进行傅立叶变换得到在焦像和离焦像的傅立叶变换谱G<sub>1</sub>(u,ν)和G<sub>2</sub>(u,ν),设定光瞳函数P<sub>1</sub>和P<sub>2</sub>,初始化泽尼克像差系数α<sub>j</sub>,设定精度为ε,将光瞳函数P<sub>1</sub>和P<sub>2</sub>分别进行傅立叶变换得到其所对应的振幅分布函数,取振幅分布函数模的平方,分别计算出在焦像和离焦像的点扩散函数h<sub>1</sub>和h<sub>2</sub>,将h<sub>1</sub>和h<sub>2</sub>分别进行傅立叶变换得到在焦像和离焦像的光学传递函数H<sub>1</sub>(u,v)和H<sub>2</sub>(u,v),推导得出目标函数E(a)的表达式如式(1)所示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>E</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><munder><mo>&Sigma;</mo><mrow><mi>u</mi><mo>,</mo><mi>v</mi><mo>&Element;</mo><msub><mi>&chi;</mi><mn>1</mn></msub></mrow></munder><mfrac><mrow><mo>|</mo><msub><mi>G</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><msup><msub><mi>H</mi><mn>2</mn></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>G</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><msup><msub><mi>H</mi><mn>1</mn></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><msup><mo>|</mo><mn>2</mn></msup></mrow><mrow><mo>|</mo><msub><mi>H</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><msup><mo>|</mo><mn>2</mn></msup><mo>+</mo><mo>|</mo><msub><mi>H</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><msup><mo>|</mo><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000750017600000011.GIF" wi="1434" he="205" /></maths>式(1)中,H<sub>1</sub><sup>*</sup>(u,v)为H<sub>1</sub>(u,v)的复共轭,H<sub>2</sub><sup>*</sup>(u,v)为H<sub>2</sub>(u,v)的复共轭;采用共轭梯度算法求解出泽尼克像差系数α<sub>j</sub>,利用此像差系数α<sub>j</sub>求得光学成像系统的像差分布,由此像差分布即可以得到一幅复原后的目标图像,再利用此像差系数α<sub>j</sub>求解出光学成像系统的光学传递函数,即为图像的退化函数H(u,v),对上述得到的复原后的目标图像进行维纳滤波处理,得到一幅二次复原后的目标清晰图像,同时对步骤一中得到的错位在焦像也进行维纳滤波处理,得到一幅复原后的错位在焦像;步骤三、将步骤二中得到的复原后的目标图像和错位在焦像按照实际位置迭加在一起,得到的初步合成的图像中有反映图像真实信息的实际数据点,相邻实际数据点之间有1个待插值点,采用三次B样条插值算法估计出待插值点的像素值,在初步合成的图像的空间域进行插值处理,得到超分辨率复原图像。
地址 130033 吉林省长春市东南湖大路3888号