主权项 |
一种应用于数字全息重构过程的对焦方法,其特征在于,该方法包括以下步骤:步骤(1)采用数字全息方法得到全息图后,处理全息图得到物光场矩阵U<sub>0</sub>;步骤(2)用肉眼或者用常见长度测量工具粗略测量后估算出一个对焦距离区间[z<sub>1</sub>,z<sub>k</sub>],确保对焦距离在该区间内,此处对焦距离为成像对象与数字图像传感器的距离。设置区间采样数为k,计算出采样间隔距离<img file="FDA0000901022350000011.GIF" wi="349" he="111" />每隔采样间隔距离delta选取一个反向传播距离,依次得到z<sub>2</sub>=z<sub>1</sub>+delta,z<sub>3</sub>=z<sub>1</sub>+2delta,z<sub>4</sub>=z<sub>1</sub>+3delta,…,z<sub>k</sub>=z<sub>1</sub>+(k‑1)*delta;步骤(3)对反向传播距离z<sub>1</sub>,z<sub>2</sub>,z<sub>3</sub>,…,z<sub>k</sub>依次采用光学衍射算法反向衍射得到光场矩阵U<sub>1</sub>,U<sub>2</sub>,U<sub>3</sub>,…,U<sub>k</sub>;步骤(4)在每个光场矩阵U<sub>x</sub>中选相同位置需要对焦的目标区域,得到新矩阵ROI<sub>x</sub>,即对焦窗口,矩阵ROI<sub>x</sub>大小为m*n,其中,x=1,2,…,k;步骤(5)通过矩阵ROI<sub>x</sub>计算实数R(z<sub>x</sub>),公式如下;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>x</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>Σ</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mo>Σ</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>|</mo><msub><mi>ROI</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo><mo>,</mo><mrow><mo>(</mo><mi>x</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000901022350000012.GIF" wi="1046" he="134" /></maths>其中,ROI<sub>x</sub>(i,j)为矩阵ROI<sub>x</sub>第i行、第j列的元素,或者通过公式如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>x</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>Σ</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mo>Σ</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>|</mo><msub><mi>ROI</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><msup><mo>|</mo><mn>2</mn></msup><mo>,</mo><mrow><mo>(</mo><mi>x</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000901022350000013.GIF" wi="1046" he="142" /></maths>其中,ROI<sub>x</sub>(i,j)为矩阵ROI<sub>x</sub>第i行、第j列的元素;步骤(6)计算数列y(z<sub>1</sub>),y(z<sub>2</sub>),…,y(z<sub>k‑1</sub>),公式如下:y(z<sub>x</sub>)=|R(z<sub>x</sub>)‑R(z<sub>x+1</sub>)|<sup>2</sup>,(x=1,2,...,k‑1),步骤(7)对数列y(z<sub>1</sub>),y(z<sub>2</sub>),…,y(z<sub>k‑1</sub>)进行处理找到对焦面。 |