主权项 |
一种基于叠层扫描的相位差波前探测和图像复原方法,其特征在于该方法包含以下步骤:第一步,将小孔径的光阑在成像系统的光瞳平面内以叠层扫描的方式进行平移,并用图像传感器依次记录下与各扫描位置所对应的受不同像差影响的子图像;第二步,使用基于叠层扫描的相位差波前探测和图像复原算法对序列子图像进行处理,得到待测的畸变波前和复原的图像,处理算法流程如下:1)根据极大似然估计或最小二乘优化理论建立目标函数E:<maths num="0001"><math><![CDATA[<mrow><mi>E</mi><mo>=</mo><munder><mo>Σ</mo><mrow><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub></mrow></munder><munderover><mo>Σ</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>Σ</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>Σ</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>Σ</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>α</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mo>Σ</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>Σ</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> |