发明名称 一种应用于数字全息重构过程的对焦方法
摘要 本发明公开了一种应用于数字全息重构过程的对焦方法。该方法是在数字全息计算处理得到物光场后,物光场反向轴向传播,对两相邻截面光场强度(或振幅)的同一对焦区域矩阵元素作差,然后对所得到的矩阵的元素平方,再对平方后的矩阵所有元素求和得到一个函数值,依次对相邻等间距的轴向截面图像作上述处理即可以得到一组函数值,找到这组函数值的最小值,其对应的距离处即为对焦面的位置。该方法可应用于数字全息仪器装置的重构过程中的对焦,具有计算简单,适应范围广,对焦准确,不需要额外增加硬件装置等优点。
申请公布号 CN105446111A 申请公布日期 2016.03.30
申请号 CN201610005898.0 申请日期 2016.01.06
申请人 中国科学院上海光学精密机械研究所 发明人 吕萌;司徒国海
分类号 G03H1/22(2006.01)I 主分类号 G03H1/22(2006.01)I
代理机构 上海新天专利代理有限公司 31213 代理人 张泽纯;张宁展
主权项 一种应用于数字全息重构过程的对焦方法,其特征在于,该方法包括以下步骤:步骤(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>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mo>&Sigma;</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>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mo>&Sigma;</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>)进行处理找到对焦面。
地址 201800 上海市嘉定区上海市800-211邮政信箱
您可能感兴趣的专利