发明名称 一种基于叠层扫描的相位差波前探测和图像复原方法
摘要 本发明涉及一种基于叠层扫描的相位差波前探测和图像复原方法,可用于恢复波前畸变,并且能够有效复原受像差影响的模糊图像。本发明将一个小孔径的光阑在成像系统的光瞳平面内以叠层扫描的方式进行移动,并用图像传感器记录下相应的受不同像差影响的子图像。通过使用基于叠层扫描的相位差波前探测和图像复原算法对记录的系列子图像进行处理,可以探测出系统的波前畸变并复原图像。本发明采用小孔径的光阑进行空间扫描以产生含有相位差异的子图像,相对于目前用于相位差法的各种技术,更容易满足奈斯奎特采样频率,像差探测范围大,并且无需离焦光路和各种衍射器件,具有系统紧凑,使用方便等优点。
申请公布号 CN105716725A 申请公布日期 2016.06.29
申请号 CN201610136602.9 申请日期 2016.03.10
申请人 中国科学院光电技术研究所 发明人 谢宗良;马浩统;任戈;亓波;史建亮;崔占刚;谭玉凤;王智鹏;何小君;董理
分类号 G01J9/00(2006.01)I 主分类号 G01J9/00(2006.01)I
代理机构 代理人
主权项 一种基于叠层扫描的相位差波前探测和图像复原方法,其特征在于该方法包含以下步骤:第一步,将小孔径的光阑在成像系统的光瞳平面内以叠层扫描的方式进行平移,并用图像传感器依次记录下与各扫描位置所对应的受不同像差影响的子图像;第二步,使用基于叠层扫描的相位差波前探测和图像复原算法对序列子图像进行处理,得到待测的畸变波前和复原的图像,处理算法流程如下:1)根据极大似然估计或最小二乘优化理论建立目标函数E:<maths num="0001"><math><![CDATA[<mrow><mi>E</mi><mo>=</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub></mrow></munder><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msup><mrow><mo>|</mo><msub><mi>D</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub></mrow></munder><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msup><mrow><mo>|</mo><msub><mi>D</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><msubsup><mi>H</mi><mi>k</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msup><mrow><mo>|</mo><msubsup><mi>H</mi><mi>k</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow>]]></math><img file="FDA0000938861550000011.GIF" wi="1139" he="263" /></maths>其中,D<sub>k</sub>(f<sub>x</sub>,f<sub>y</sub>)是调制光阑处在第k个扫描位置时图像传感器采集的子图像的傅里叶变换,H<sub>k</sub>(f<sub>x</sub>,f<sub>y</sub>)是调制光阑处在第k个扫描位置时对应的光学通道的光学传递函数,<img file="FDA0000938861550000012.GIF" wi="234" he="87" />是H<sub>k</sub>(f<sub>x</sub>,f<sub>y</sub>)的共轭,K为扫描的次数,f<sub>x</sub>,f<sub>y</sub>为像平面频域的坐标;2)根据傅里叶光学,估计出第k个光学通道的点扩散函数,系统的波前相位分布<img file="FDA0000938861550000013.GIF" wi="154" he="71" />可以由泽尼克多项式来表征:<img file="FDA0000938861550000014.GIF" wi="438" he="133" />其中Z<sub>n</sub>(u,v)表示第n阶泽尼克多项式,α<sub>n</sub>表示第n阶泽尼克多项式的系数;3)目标函数E演变为关于变量α<sub>n</sub>的函数E(α),采用诸如SPGD,遗传算法,模拟退火,神经网络等最优化算法寻找使目标函数最小的解,即可求得待测的波前畸变<img file="FDA0000938861550000016.GIF" wi="63" he="50" />4)根据求得的波前畸变,利用下式可以复原图像R:<maths num="0002"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>,</mo><mi>&alpha;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>D</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><msubsup><mi>H</mi><mi>k</mi><mo>*</mo></msubsup><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>.</mo></mrow>]]></math><img file="FDA0000938861550000015.GIF" wi="836" he="263" /></maths>
地址 610209 四川省成都市双流350信箱