主权项 |
1、用于哈特曼波前探测器的基于DFT的质心偏移量检测方法,其核心是基于DFT的质心偏移量方法,该方法的流程为:(1)将标准平面波入射到子孔径单元为M×M的哈特曼传感器上获取的参考全靶面图像,取其中某一要检测质心偏移量的子孔径所对应的图像,即分辨率为N×N的参考图像IC,其中M,N均为整数;(2)将分辨率均为N×N的IC的像素点的灰度值分别在两个待检测的方向叠加,如X方向和Y方向,得到两个1×N的向量Icx和Icy:<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>Icx</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>IC</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>0,1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mi>Icy</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>IC</mi><mrow><mo>(</mo><mi>j</mi><mo>,</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>0,1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>(3)对此两向量Icx和Icy分别作DFT,得到两个1×N的向量Fcx和Fcy:<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>Fcx</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>Icx</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msubsup><mi>W</mi><mi>N</mi><mi>kn</mi></msubsup><mo>,</mo><mi>k</mi><mo>=</mo><mn>0,1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mi>Fcy</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>Icy</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msubsup><mi>W</mi><mi>N</mi><mi>kn</mi></msubsup><mo>,</mo><mi>k</mi><mo>=</mo><mn>0,1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>(4)分别计算Fcx和Fcy各分量的相位,得到两个1×N的向量Cx和Cy,即参考图像的相位向量,这样,对参考图像的处理完成;(5)将待检测波面入射到该哈特曼传感器上获取待检测全靶面图像,取其中与参考图像相同子孔径所对应的图像,即分辨率为N×N的待检测的图像ID;(6)同对参考图像的处理一样,先将待检测的图像ID的像素点的灰度值分别在两个待检测的方向叠加,如X方向和Y方向,得到两个1×N的向量Idx和Idy:<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>Idx</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>ID</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>0,1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mi>Idy</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>ID</mi><mrow><mo>(</mo><mi>j</mi><mo>,</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>0,1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>(7)对此两向量Idx和Idy分别作DFT,得到两个1×N的向量Fdx和Fdy:<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>Fdx</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>Idx</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msubsup><mi>W</mi><mi>N</mi><mi>kn</mi></msubsup><mo>,</mo><mi>k</mi><mo>=</mo><mn>0,1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mi>Fdy</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>Idy</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msubsup><mi>W</mi><mi>N</mi><mi>kn</mi></msubsup><mo>,</mo><mi>k</mi><mo>=</mo><mn>0,1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>(8)分别计算Fdx和Fdy各分量的相位,得到两个1×N的向量Dx和Dy,即待检测图像的相位向量;(9)分别计算Dx与Dy及Cx与Cy的各分量的相位差,得到两个1×N的向量DDx和DDy:<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>DDx</mi><mo>=</mo><mi>Dx</mi><mo>-</mo><mi>Cx</mi></mtd></mtr><mtr><mtd><mi>DDy</mi><mo>=</mo><mi>Dy</mi><mo>-</mo><mi>Cy</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>(10)分别对DDx和DDy解卷绕,得到DDx’和DDy’;(11)分别根据DDx’和DDy’的前T项,一般可取低频项,建议取前[N/5]至[N/4]项,[.]表取整,进行一次拟合,得到斜率kx和ky,可得待检测图像的质心与参考图像的质心在x和y方向的相对偏移量dx,dy分别为<math><mrow><mi>dx</mi><mo>=</mo><mo>-</mo><mi>kx</mi><mo>·</mo><mfrac><mi>N</mi><mn>2</mn></mfrac></mrow>和<math><mrow><mi>dy</mi><mo>=</mo><mo>-</mo><mi>ky</mi><mo>·</mo><mfrac><mi>N</mi><mn>2</mn></mfrac></mrow>个像素;(12)重复(1)到(11),对子孔径单元为M×M的哈特曼探测器的其它子孔径进行相同的处理,得到所有重构波前所需的子孔径在x和y方向的偏移量。 |