CN107966111A - 一种用于数字全息重构过程的自动对焦方法 - Google Patents

一种用于数字全息重构过程的自动对焦方法 Download PDF

Info

Publication number
CN107966111A
CN107966111A CN201711169584.5A CN201711169584A CN107966111A CN 107966111 A CN107966111 A CN 107966111A CN 201711169584 A CN201711169584 A CN 201711169584A CN 107966111 A CN107966111 A CN 107966111A
Authority
CN
China
Prior art keywords
focal distance
digital hologram
curve
distance
section
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
CN201711169584.5A
Other languages
English (en)
Inventor
文永富
程灏波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201711169584.5A priority Critical patent/CN107966111A/zh
Publication of CN107966111A publication Critical patent/CN107966111A/zh
Withdrawn legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Holo Graphy (AREA)

Abstract

本发明公开了一种用于数字全息重构过程的自动对焦方法。该方法包括如下步骤:得到待处理数字全息图后,估算对焦距离并设置第一对焦距离区间、第二对焦距离区间、采样数和采样间隔距离,得到一系列临时对焦距离;分别对第一对焦区间和第一对焦区间内的临时对焦距离进行数字全息重构,得到光强图像序列;分别计算光强图像序列的自聚焦评价函数,得到第一评价函数曲线和第二评价函数曲线,进而作差得到差动评价曲线;在差动评价曲线横坐标为零的点附近选取一段曲线,进行线性拟合;拟合曲线过零点的位置即为图像聚焦位置。本发明可用于数字全息重构过程中的自动对焦,具有技术简单,定位精度高、速度快的特点。

Description

一种用于数字全息重构过程的自动对焦方法
技术领域
本发明涉及一种自动对焦方法,尤其涉及一种用于数字全息重构过程的自动对焦方法,属于数字全息测量领域。
技术背景
数字全息术与传统光学全息原理一致,区别在于前者的记录介质是光电转换器件CCD或CMOS,代替了传统的全息银盐干板。与传统的光学全息相比,数字全息的显著特点是不需要显影、定影与漂白等化学处理过程,可以根据复振幅分布采用数字全息重构算法同时获得物体的光强图像以及相位信息。
在数字全息重构过程中,通过模拟光学衍射来实现对被测物体的光波场重构。在数值重构过程中,重构距离的选取至关重要。重构时只有准确知道被测物体到数字图像传感器之间的距离,才能准确重构出被测样品的光强图像以及相位信息。如果偏离了这个距离,重构出来的光强图像由于处于离焦状态而变得模糊。传统自动对焦方法,通常是通过依次改变重构距离并进行数字全息重构,得到不同重构距离下的光强图像,然后通过采用合适的自聚焦判断函数,得到一系列的函数值,进而得到一条评价函数曲线,通过找到评价函数曲线的极值点,其对应的距离即为对焦面的位置。然而,传统自动对焦方法利用评价函数得到的评价函数曲线,其曲线峰值点处往往较为平滑,不利于曲线拟合,且峰值处斜率往往比较小,线性度差,不利于准确判定对焦面的位置;此外,评价函数的曲线斜率关于曲线峰值点对称,无法根据函数值的变化快速判断调焦运动方向,因而也不利于软件快速自动搜索处理。
发明内容
本发明针对现有技术的不足,提出一种用于数字全息重构过程的自动对焦方法,其目的在于实现数字全息测量过程中的自动聚焦,便于快速准确获取被测物体表面三维形貌。
为达到上述目的,本发明通过以下技术方案实现的,一种应用于数字全息重构过程的自动对焦方法,该方法包括以下步骤:
步骤(1):在得到待处理的数字全息图后,用肉眼或者长度测量工具粗略测量后估算出一个重构距离zc,以zc为基准设置第一对焦距离区间[z1,zk];设置对焦距离区间采样数k,计算采样间隔距离dz=(zk-z1)/(k-1);针对第一对焦距离区间[z1,zk],每隔dz选取一个距离,得到不同临时对焦距离z2=z1+dz,z3=z1+2*dz,…,zk=z1+(k-1)*dz。
设置第二对焦距离区间[z1’,zk’];针对第二对焦距离区间[z1’,zk’],每隔距离dz选取一个距离,得到不同临时对焦距离z2’=z1’+dz,z3’=z1’+2*dz,…,zk’=z1’+(k-1)*dz;
步骤(2):针对第一对焦距离区间[z1,zk],采用数字全息重构算法,分别得到不同临时对焦距离下的物光场U1,U2,…,Uk。
针对第二对焦距离区间[z1’,zk’],采用数字全息数值重构算法,分别得到不同临时对焦距离下的物光场U1’,U2’,…,Uk’;
步骤(3):利用物光场U1,U2,…,Uk计算得到不同临时对焦距离下的光强图像I1,I2,…,Ik。
利用物光场U1’,U2’,…,Uk’计算得到不同临时对焦距离下的光强图像I1’,I2’,…,Ik’。
步骤(4):利用自聚焦评价函数计算出每个光强图像I1,I2,…,Ik的函数值C1,C2,…,Ck,并以函数值为纵坐标,临时对焦距离为横坐标,建立第一条评价函数曲线。
利用自聚焦评价函数计算出每个光强图像I1’,I2’,…,Ik’的函数值C1’,C2’,…,Ck’,并以函数值为纵坐标,临时对焦距离为横坐标,建立第二条评价函数曲线。
步骤(5):将步骤(4)得到的第一条评价曲线和第二条评价函数曲线进行作差处理,得到差动评价曲线。
步骤(6):在差动评价曲线横坐标为零的点附近选取一段近似于直线的曲线,并进行线性拟合。
步骤(7):搜索拟合后的直线过零点的位置,零点所在得位置准确对应图像对焦面的位置,即实现准确聚焦。
本发明的有益效果在于:由于过零点附近曲线近似满足线性,便于对测量数据进行线性拟合;同时,过零点附近曲线斜率较大,便于对焦距离的确定;另外,由于过零点前后曲线符号相反,成双极性分布,便于判断对焦搜索方向,方便快速自动对焦;此外,整个对焦过程通过数值计算实现,无需引入额外的硬件设备。
附图说明
图1为本发明一种用于数字全息重构过程的自动对焦方法的流程图;
图2为具体实施例中采用的数字全息图;
图3为具体实施例中采用说明书发明内容里步骤(4)后本发明得到的第一条评价函数曲线;
图4为具体实施例中采用说明书发明内容里步骤(4)后本发明得到的第二条评价函数曲线;
图5为具体实施例中采用说明书发明内容里步骤(5)后本发明得到的差动评价曲线;
图6为具体实施例中图2给出的数字全息图采用本发明得到的聚焦光强图像。利用本发明得到的聚焦位置进行数字全息重构得到的光强图像。
具体实施方式
下面结合附图及具体实施方式进一步详细介绍本发明。应当说明的是,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此间未构成冲突就可以相互组合。
图1是本发明一种用于数字全息重构过程的自动对焦方法的流程图。
步骤(1):得到待处理的数字全息图,如图2所示;用长度测量工具粗略测量后估算出一个重构距离zc=107mm,以zc为基准设置第一对焦距离区间[104.6mm,112.6mm];设置对焦距离区间采样数k=400,计算采样间隔距离dz=0.02mm;设置第二对焦距离区间[103.4mm,111.4mm]。
针对第一对焦距离区间[104.6mm,112.6mm],每隔距离dz=0.02mm选取一个距离,得到不同临时对焦距离z2=104.62mm,z3=104.64mm,…,zk=112.6mm;
针对第二对焦距离区间[103.4mm,111.4mm],每隔距离dz=0.02mm选取一个距离,得到不同临时对焦距离z2’=103.42mm,z3’=103.44mm,…,zk’=111.4mm;
步骤(2):针对第一对焦距离区间[104.6mm,112.6mm],采用数字全息数值重构算法,分别得到不同临时对焦距离z2=104.62mm,z3=104.64mm,…,zk=112.6mm下的物光场U1,U2,…,Uk;
针对第二对焦距离区间[103.4mm,111.4mm],采用数字全息数值重构算法,分别得到不同临时对焦距离z2’=103.42mm,z3’=103.44mm,…,zk’=111.4mm下的物光场U1’,U2’,…,Uk’。
步骤(3):利用物光场U1,U2,…,Uk计算得到不同临时对焦距离下的光强图像I1,I2,…,Ik。
利用物光场U1’,U2’,…,Uk’计算得到不同临时对焦距离下的光强图像I1’,I2’,…,Ik’。
步骤(4):利用自聚焦评价函数计算出每个光强图像I1,I2,…,Ik的函数值C1=0.0221,C2=0.0222,…,Ck=0.0001,并以函数值为纵坐标,临时对焦距离为横坐标,建立第一条评价函数曲线,如图3所示。
利用自聚焦评价函数计算出每个光强图像I1’,I2’,…,Ik’的函数值C1’=0.0002,C2’=0.002,…,Ck’=0.0447,并以函数值为纵坐标,临时对焦距离为横坐标,建立第二条评价函数曲线,如图4所示。
步骤(5):将步骤(4)得到的第一条评价曲线和第二条评价函数曲线进行作差处理,得到差动评价曲线,如图5所示。
步骤(6):在差动评价曲线横坐标为零的点附件选取一段曲线,并进行最小二乘线性拟合,得到一条拟合直线。
步骤(7):找到拟合直线过零点处对应的横坐标为108mm,即为被测样品对焦面的位置,完成自动聚焦。图6为对图2给出的数字全息图采用本发明得到的聚焦光强图像。
尽管特定参考实施例详细的描述了本发明,但在此描述的本发明实施例的意图不是详尽或者局限于所公开的具体形式。相反,所选用的用于说明问题的实施例是为了使本技术领域内的技术人员实施本发明而选择的。在不脱离本发明的权利要求所描述和限定的本发明的实质范围的情况下,存在变型例或修改例。

Claims (5)

1.一种用于数字全息重构过程的自动对焦方法,其特征在于,该方法包括如下步骤:
步骤(1):在得到待处理的数字全息图后,用肉眼或者长度测量工具粗略测量后估算出一个重构距离zc,以zc为基准设置第一对焦距离区间[z1,zk];设置对焦距离区间采样数k,计算采样间隔距离dz=(zk-z1)/(k-1);针对第一对焦距离区间[z1,zk],每隔dz选取一个距离,得到不同临时对焦距离z2=z1+dz,z3=z1+2*dz,…,zk=z1+(k-1)*dz。
设置第二对焦距离区间[z1’,zk’];针对第二对焦距离区间[z1’,zk’],每隔距离dz选取一个距离,得到不同临时对焦距离z2’=z1’+dz,z3’=z1’+2*dz,…,zk’=z1’+(k-1)*dz;
步骤(2):针对第一对焦距离区间[z1,zk],采用数字全息重构算法,分别得到不同临时对焦距离下的物光场U1,U2,…,Uk。
针对第二对焦距离区间[z1’,zk’],采用数字全息数值重构算法,分别得到不同临时对焦距离下的物光场U1’,U2’,…,Uk’;
步骤(3):利用物光场U1,U2,…,Uk计算得到不同临时对焦距离下的光强图像I1,I2,…,Ik。
利用物光场U1’,U2’,…,Uk’计算得到不同临时对焦距离下的光强图像I1’,I2’,…,Ik’。
步骤(4):利用自聚焦评价函数计算出每个光强图像I1,I2,…,Ik的函数值C1,C2,…,Ck,并以函数值为纵坐标,临时对焦距离为横坐标,建立第一条评价函数曲线。
利用自聚焦评价函数计算出每个光强图像I1’,I2’,…,Ik’的函数值C1’,C2’,…,Ck’,并以函数值为纵坐标,临时对焦距离为横坐标,建立第二条评价函数曲线。
步骤(5):将步骤(4)得到的第一条评价曲线和第二条评价函数曲线进行作差处理,得到差动评价曲线。
步骤(6):在差动评价曲线横坐标为零的点附近选取一段近似于直线的曲线,并进行线性拟合。
步骤(7):搜索拟合后的直线过零点的位置,零点所在得位置准确对应图像对焦面的位置,即实现准确聚焦。
2.根据权利要求1所述的一种用于数字全息重构过程的自动对焦方法,其特征在于,所述重构距离zc是指被测物体和数字图像传感器的距离。
3.根据权利要求1所述的一种用于数字全息重构过程的自动对焦方法,其特征在于,所述的第一对焦距离区间[z1,zk],满足关系式z1<zc<zk。
4.根据权利要求1所述的一种用于数字全息重构过程的自动对焦方法,其特征在于,所述的第二对焦距离区间[z1’,zk’],满足关系式z1’<zc<zk’且同时还满足z1’<z1,zk’<zk。
5.根据权利要求1所述的一种用于数字全息重构过程的自动对焦方法,其特征在于,所述的自聚焦评价函数可以是方差灰度算子、Roberts梯度算子、Sobel梯度算子、梯度向量平方算子、二阶拉普拉斯算子或者傅里叶变换频谱加权对数算子。
CN201711169584.5A 2017-11-21 2017-11-21 一种用于数字全息重构过程的自动对焦方法 Withdrawn CN107966111A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711169584.5A CN107966111A (zh) 2017-11-21 2017-11-21 一种用于数字全息重构过程的自动对焦方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711169584.5A CN107966111A (zh) 2017-11-21 2017-11-21 一种用于数字全息重构过程的自动对焦方法

Publications (1)

Publication Number Publication Date
CN107966111A true CN107966111A (zh) 2018-04-27

Family

ID=61999822

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711169584.5A Withdrawn CN107966111A (zh) 2017-11-21 2017-11-21 一种用于数字全息重构过程的自动对焦方法

Country Status (1)

Country Link
CN (1) CN107966111A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108710205A (zh) * 2018-05-16 2018-10-26 电子科技大学 一种基于边缘灰度差异函数的光学扫描全息自聚焦方法
CN109581373A (zh) * 2018-12-17 2019-04-05 中国人民解放军国防科技大学 一种全息穿透成像雷达的自聚焦成像方法
CN109828444A (zh) * 2019-03-28 2019-05-31 哈尔滨工业大学 一种用于多目标物数字全息测量的自动聚焦方法
CN110675383A (zh) * 2019-09-24 2020-01-10 哈尔滨工业大学 用于数字全息中存在干扰物的目标物快速识别聚焦方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001050201A1 (en) * 2000-01-04 2001-07-12 Ut-Battelle L.L.C. Improvements to acquisition and replay systems for direct-to-digital holography and holovision
CN102445812A (zh) * 2010-10-13 2012-05-09 致茂电子(苏州)有限公司 高速自动对焦方法
CN103038692A (zh) * 2010-06-24 2013-04-10 皇家飞利浦电子股份有限公司 基于差分测量的自动聚焦
EP2615484A1 (en) * 2012-01-13 2013-07-17 Sony Corporation Automatic focusing apparatus and method with calibration and slope correction
CN103914870A (zh) * 2014-02-28 2014-07-09 天津工业大学 一种基于新聚焦评价函数的通用性强的全息图自动重建方法
CN105446111A (zh) * 2016-01-06 2016-03-30 中国科学院上海光学精密机械研究所 一种应用于数字全息重构过程的对焦方法
CN105578029A (zh) * 2015-09-01 2016-05-11 闽南师范大学 一种多尺度变步长的自动对焦搜索算法据传输装置和方法
CN105635590A (zh) * 2016-02-29 2016-06-01 中国工程物理研究院流体物理研究所 一种基于数字全息重构算法的对焦方法及装置
CN107179645A (zh) * 2017-06-30 2017-09-19 天津市亚安科技有限公司 一种用于相机的自动聚焦方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001050201A1 (en) * 2000-01-04 2001-07-12 Ut-Battelle L.L.C. Improvements to acquisition and replay systems for direct-to-digital holography and holovision
CN103038692A (zh) * 2010-06-24 2013-04-10 皇家飞利浦电子股份有限公司 基于差分测量的自动聚焦
CN102445812A (zh) * 2010-10-13 2012-05-09 致茂电子(苏州)有限公司 高速自动对焦方法
EP2615484A1 (en) * 2012-01-13 2013-07-17 Sony Corporation Automatic focusing apparatus and method with calibration and slope correction
CN103914870A (zh) * 2014-02-28 2014-07-09 天津工业大学 一种基于新聚焦评价函数的通用性强的全息图自动重建方法
CN105578029A (zh) * 2015-09-01 2016-05-11 闽南师范大学 一种多尺度变步长的自动对焦搜索算法据传输装置和方法
CN105446111A (zh) * 2016-01-06 2016-03-30 中国科学院上海光学精密机械研究所 一种应用于数字全息重构过程的对焦方法
CN105635590A (zh) * 2016-02-29 2016-06-01 中国工程物理研究院流体物理研究所 一种基于数字全息重构算法的对焦方法及装置
CN107179645A (zh) * 2017-06-30 2017-09-19 天津市亚安科技有限公司 一种用于相机的自动聚焦方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
BING-ZHAO ET AL: "Spectral analysis of sampled signals in the linear cononical transform doman", 《MATHEMATICAL PROBLEMS IN ENGINNERING》 *
刘长庚等: "数字全息成像中基于导数的自动对焦算法", 《中国激光》 *
高帆等: "基于图像识别的形态学检测显微镜自动对焦系统", 《企业技术开发》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108710205A (zh) * 2018-05-16 2018-10-26 电子科技大学 一种基于边缘灰度差异函数的光学扫描全息自聚焦方法
CN108710205B (zh) * 2018-05-16 2019-06-21 电子科技大学 一种基于边缘灰度差异函数的光学扫描全息自聚焦方法
CN109581373A (zh) * 2018-12-17 2019-04-05 中国人民解放军国防科技大学 一种全息穿透成像雷达的自聚焦成像方法
CN109828444A (zh) * 2019-03-28 2019-05-31 哈尔滨工业大学 一种用于多目标物数字全息测量的自动聚焦方法
CN110675383A (zh) * 2019-09-24 2020-01-10 哈尔滨工业大学 用于数字全息中存在干扰物的目标物快速识别聚焦方法
CN110675383B (zh) * 2019-09-24 2023-04-14 哈尔滨工业大学 用于数字全息中存在干扰物的目标物快速识别聚焦方法

Similar Documents

Publication Publication Date Title
CN107966111A (zh) 一种用于数字全息重构过程的自动对焦方法
CN106767489B (zh) 数字散斑干涉面内微小动态形变测量系统及测量方法
CN107991242B (zh) 一种基于偏振分光棱镜的测量样品偏振态的方法与系统
CN111238403A (zh) 一种基于光场子孔径条纹图像的三维重建方法及装置
CN112097678B (zh) 基于频率盲估计的多表面面形测量方法
CN109035345A (zh) 基于高斯过程回归的tof相机距离校正方法
CN109242774A (zh) 一种基于多维空间不变特征的平板类零件点云拼接方法
Liu et al. Coaxial projection profilometry based on speckle and fringe projection
Wang et al. Quantitative interferometric microscopy with two dimensional Hilbert transform based phase retrieval method
Dong et al. Automatic filtering for zero-order and twin-image elimination in off-axis digital holography
CN105403508B (zh) 基于合成相位传递函数的非干涉相位成像方法
CN113298700A (zh) 一种在散射场景中的高分辨图像重构方法
CN105446111B (zh) 一种应用于数字全息重构过程的对焦方法
CN104034435B (zh) 一种测量高斯光束参数的方法
Yuan et al. Phase-shifting coding method for absolute phase retrieval
CN107730565B (zh) 一种基于oct图像的材料内在光谱特征提取方法
CN106502074B (zh) 一种用于像面数字全息显微测量的自动聚焦方法
Lertrusdachakul et al. Dynamic (de) focused projection for three-dimensional reconstruction
Hong et al. Phase retrieval based on transport of intensity equation and image interpolation
Singh et al. Nanoscale measurement of in-plane and out-of-plane displacements of microscopic object by sensor fusion
US11385164B2 (en) Method for calibrating an analysis device, and associated device
CN108051089B (zh) 一种基于渥拉斯顿棱镜的测量样品偏振态的方法与系统
Oh et al. Quantitative phase imaging via the holomorphic property of complex optical fields
Wang et al. Unwrapping free Hilbert transform-based phase retrieval method in quantitative interferometric microscopy
Li et al. Flexible and universal autofocus based on amplitude difference of fractional Fourier transform

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20180427