CN106502074A - 一种用于像面数字全息显微测量的自动聚焦方法 - Google Patents
一种用于像面数字全息显微测量的自动聚焦方法 Download PDFInfo
- Publication number
- CN106502074A CN106502074A CN201610983732.6A CN201610983732A CN106502074A CN 106502074 A CN106502074 A CN 106502074A CN 201610983732 A CN201610983732 A CN 201610983732A CN 106502074 A CN106502074 A CN 106502074A
- Authority
- CN
- China
- Prior art keywords
- image
- auto
- hologram
- measuring
- image planes
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000012360 testing method Methods 0.000 claims abstract description 17
- 238000001914 filtration Methods 0.000 claims abstract description 4
- 238000001228 spectrum Methods 0.000 claims abstract description 4
- 230000009466 transformation Effects 0.000 claims abstract description 4
- 230000003287 optical effect Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims 1
- 238000012876 topography Methods 0.000 abstract description 7
- 238000001514 detection method Methods 0.000 abstract description 5
- 238000005259 measurement Methods 0.000 description 18
- 239000000523 sample Substances 0.000 description 14
- 230000008569 process Effects 0.000 description 12
- 238000003384 imaging method Methods 0.000 description 6
- 238000009647 digital holographic microscopy Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000012472 biological sample Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000001093 holography Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/0005—Adaptation of holography to specific applications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/0005—Adaptation of holography to specific applications
- G03H2001/0033—Adaptation of holography to specific applications in hologrammetry for measuring or analysing
Abstract
本发明公开了一种用于像面数字全息显微测量的自动聚焦方法,该方法包括如下步骤:沿物镜主光轴方向移动待测样品并记录一系列等间隔的全息图像;通过傅里叶变换、角谱滤波、逆傅里叶变换等操作获得每张全息图像所对应的物光强度图像;计算每张强度图像的自相关算子;使用多项式拟合自相关算子相对于全息图记录位置的数据变化趋势,拟合曲线的最高点即自相关算子最大值,其对应的图像位置即聚焦平面位置;通过程序控制电动升降台,将待测样品表面移动到聚焦平面的位置,获得聚焦的像面全息图像。本发明一种用于像面数字全息显微测量的自动聚焦方法,实现了自动像面全息显微检测,获得精度更高,实时性更好的表面形貌测量方法。
Description
技术领域
本发明属于数字全息显微测量领域,更具体地,涉及一种用于像面数字全息显微测量的自动聚焦方法。
背景技术
伴随着先进制造技术的发展,以微机电系统(MEMS)、微光机电系统(MOEMS)为代表的元器件总体尺寸越来越小,器件的表面形貌和微观缺陷对其性能的影响也越来越明显,变得不可忽视。表面形貌测量最常用的方法包括机械探针法、光学测量法等。其中,数字全息显微测量技术(DHM),凭借其可实时测量、非接触测量、三维测量等优势,十分适合于微小元器件表面形貌实时检测。
目前常用的数字全息显微测量原理中,多数使用的是菲涅尔全息图,其采用了图像传感器记录平面与显微光路成像平面不重合的记录方式,并在后续的三维图像重建过程中通过光场传播的方法计算成像平面所成的光场数据。表面上这种检测方法具有检测过程中无需聚焦的优势,但实际上测量结果的三维数据计算过程中仍然需要知道图像传感器记录平面与显微光路成像平面之间的距离。一方面,由于该距离无法精确测量,目前多采用迭代计算成像的方法寻找成像距离,其计算量过大不利于实时测量;另一方面,光场传播的计算通常采用麦克斯韦方程组的近似解,对重建图像的精度损失相对更大。现有理论证明,采用像面全息记录方式能获得最高的精度,且由于图像传感器记录平面与显微光路成像平面之间距离为0,无需后续的光场传播计算,有效减少了三维成像时间,更好的保证了实时测量效果。故将图像聚焦过程从后续图像计算转移到测量前的聚焦准备过程中,更利于得到精度高、速度快的测量结果。然而,数字全息测量中,图像传感器获得的是带干涉条纹的全息图像,与普通成像系统获得的自然图像有本质上的区别,传统的聚焦判断方法无法直接使用到像面数字全息显微测量技术中。
专利文献CN105446111A公开了一种应用于数字全息重构过程的对焦方法,其是在数字全息计算处理得到物光场后,物光场反向轴传播,对两相邻截面光场强度的同一对焦区域矩阵元素作差,然后对所得到的矩阵的元素平方,再对平方后的矩阵所有元素求和得到一个函数值,依次对相邻等间距的轴向截面图像上述处理即可得到一组函数值,找到这组函数值的最小值,其对应的距离处即为对焦面的位置。该方法计算简单、适用方位广,对焦准确,但也存在如下缺陷或不足:
(1)对焦过程在拍摄完全息图像之后进行,使用的是菲涅尔全息图,理论上分辨率不如像面全息图;
(2)对焦过程结束之后的测量中,每张图像都要进行光场传播的模拟计算,增加了计算量,测量速度不如无需进行光场传播计算的像面全息技术。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种用于像面数字全息显微测量的自动聚焦方法,其目的在于提供一种可用于全息显微测量前自动聚焦的方法,由此实现自动像面全息显微检测,获得精度更高,实时性更好的表面形貌测量方法。
为了实现上述目的,本发明提供了一种用于像面数字全息显微测量的自动聚焦方法,该方法包括如下步骤:
S1:将待测样品放于载物台上,调节载物台,使待测样品表面处于焦平面下方位置;
S2:控制电动载物台沿物镜光轴方向向上等间隔移动待测样品n次,每移动一次,拍摄一张全息图,获得一系列聚焦位置不同的全息图;
S3:获得每张全息图所对应的物光强度图像Ii,按照下列公式计算每张强度图像的自相关算子Ci:
其中,M,N表示物光强度图像Ii的大小,p,q为整数,f(p,q)表示像素点(p,q)处的图像灰度值,k为一个描述自相关程度的系数;
S4:判断Ci是否单调递增,如果Ci单调递增,则表明扫描还未经过焦平面,则继续增大n值,重复步骤S2和S3;如果Ci先增后减,则说明扫描经过了焦平面,继续步骤S5;
S5:以该初始位置为原点,各个全息图记录的位置为横坐标,各全息图所对应的自相关算子Ci为纵坐标值,对数据进行多项式拟合,拟合曲线的最高点所对应的横坐标即为待测样品的聚焦位置,通过控制电动升降台,直接将待测样品移动到此位置,即可获得聚焦的像面全息图。
进一步地,步骤S2中所述聚焦为测量前聚焦。
进一步地,步骤S3中所述物光强度图像Ii通过对全息图像进行傅里叶变换、角谱滤波、逆傅里叶变换数据处理获得。
进一步地,步骤S3中所述自相关算子Ci也可为方差灰度算子、Roberts梯度算子、Sobel梯度算子或自然图像聚焦判断算子。
进一步地,步骤S3中所述自相关程度的系数k取值范围为5~30,k值越大,图像自相关的程度越小。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
(1)本发明的方法,提供一种可用于全息显微测量前自动聚焦的方法,由此实现自动像面全息显微检测,获得精度更高,实时性更好的表面形貌测量方法。
(2)通过提前聚焦,将通常采用的菲涅尔全息图转换为像面全息图,避免了光场传播的计算,提高了测量精度与三维数据重建速度。
(3)本发明在实施的过程是在现有的设备——数字全息显微镜下实施的,没有使用额外的光路设备,同时降低了数据处理过程中的运算量,节省了硬件成本。
(4)本发明通过采用像面数字全息显微测量自动聚焦方法,能够有效的通过目前常用的数字全息显微装置对被测样品表面形貌及微小表面缺陷进行高精度、快速测量与识别,对于表面形貌参数、微观表面缺陷的检测有着重要的意义,有效地推动微器件可靠性的进一步发展。
附图说明
图1为本发明实施例的一种用于像面数字全息显微测量的自动聚焦方法的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
图1为本发明实施例的一种用于像面数字全息显微测量的自动聚焦方法的流程图。如图1所示,该方法的具体步骤:
步骤1:将待测样品放于由计算机控制的电动载物台上,通过肉眼观察结合电动调节载物台,使待测区域位于图像传感器测量范围内,并连续、快速上下调节载物台,观察聚焦平面大致位置,最终使待测样品表面处于焦平面下方1mm附近某位置,记该位置为X0;
步骤2:让程序自动控制电动载物台沿物镜光轴方向向上等间隔0.1mm移动待测样品10次,每移动一次,拍摄一张全息图,获得一系列聚焦位置不同的全息图Hi(i=1,2...10);
步骤3:通过傅里叶变换、角谱滤波、逆傅里叶变换等操作获得每张全息图所对应的物光强度图像Ii;按照下列公式计算每张强度图像的自相关算子Ci,
其中,M,N表示强度图Ii的大小(横纵坐标像素值,本实施例中M,N均取512),p,q为整数,f(p,q)表示像素点(p,q)处的图像灰度值,k为一个描述自相关程度的系数,本实施例中取5,Ci的最大值表示最佳的清晰度,即图像聚焦的判断依据;
步骤4:判断Ci是否单调递增,如果是,则说明扫描还未经过焦平面,则继续增大n值,重复步骤2和步骤3;如果Ci是先增后减,则说明扫描经过了焦平面,继续步骤5;
步骤5:以位置X0为原点,各个全息图记录的位置为横坐标,各全息图所对应的强度图计算得到的自相关算子Ci为纵坐标值,对数据进行多项式拟合。拟合曲线的最高点所对应的横坐标Xfoc即为待测样品的聚焦位置,通过控制电动升降台,直接将待测样品移动到此位置,即可获得聚焦的像面全息图。
本发明提供一种可用于全息显微测量前自动聚焦的方法,由此实现自动像面全息显微检测,获得精度更高,实时性更好的表面形貌测量方法。数字全息显微测量装置可为反射式或透射式。一次聚焦后,样品的动态变化过程中(如生物样品的水平运动,MEMS样品的振动测试,流水生产线上样品的水平移动)无需重新聚焦,故三维测量数据重建过程中无需光场传播计算,提高了测量效率。以上发明构思中无需使用额外的参考光遮挡设备,无需高性能的数据处理设备,成本低。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种用于像面数字全息显微测量的自动聚焦方法,其特征在于:该方法包括下列步骤:
S1:将待测样品放于载物台上,调节载物台,使待测样品表面处于焦平面下方位置;
S2:控制电动载物台沿物镜光轴方向向上等间隔移动待测样品n次,每移动一次,拍摄一张全息图,获得一系列聚焦位置不同的全息图;
S3:获得每张全息图所对应的物光强度图像Ii,按照下列公式计算每张强度图像的自相关算子Ci:
其中,M,N表示物光强度图像Ii的大小,p,q为整数,f(p,q)表示像素点(p,q)处的图像灰度值,k为一个描述自相关程度的系数;
S4:判断Ci是否单调递增,如果Ci单调递增,则表明扫描还未经过焦平面,则继续增大n值,重复步骤S2和S3;如果Ci先增后减,则说明扫描经过了焦平面,继续步骤S5;
S5:以该初始位置为原点,各个全息图记录的位置为横坐标,各全息图所对应的自相关算子Ci为纵坐标值,对数据进行多项式拟合,拟合曲线的最高点所对应的横坐标即为待测样品的聚焦位置,通过控制电动升降台,直接将待测样品移动到此位置,即可获得聚焦的像面全息图。
2.根据权利要求1所述的一种用于像面数字全息显微测量的自动聚焦方法,其特征在于:步骤S2中所述聚焦为测量前聚焦。
3.根据权利要求1或2所述的一种用于像面数字全息显微测量的自动聚焦方法,其特征在于:步骤S3中所述物光强度图像Ii通过对全息图像进行傅里叶变换、角谱滤波、逆傅里叶变换数据处理获得。
4.根据权利要求1-3中任一项所述的一种用于像面数字全息显微测量的自动聚焦方法,其特征在于:步骤S3中所述自相关算子Ci也可为方差灰度算子、Roberts梯度算子、Sobel梯度算子或自然图像聚焦判断算子。
5.根据权利要求1-4中任一项所述的一种用于像面数字全息显微测量的自动聚焦方法,其特征在于:步骤S3中所述自相关程度的系数k取值范围为5~30,k值越大,图像自相关的程度越小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610983732.6A CN106502074B (zh) | 2016-11-08 | 2016-11-08 | 一种用于像面数字全息显微测量的自动聚焦方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610983732.6A CN106502074B (zh) | 2016-11-08 | 2016-11-08 | 一种用于像面数字全息显微测量的自动聚焦方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106502074A true CN106502074A (zh) | 2017-03-15 |
CN106502074B CN106502074B (zh) | 2019-04-12 |
Family
ID=58324130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610983732.6A Active CN106502074B (zh) | 2016-11-08 | 2016-11-08 | 一种用于像面数字全息显微测量的自动聚焦方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106502074B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109828444A (zh) * | 2019-03-28 | 2019-05-31 | 哈尔滨工业大学 | 一种用于多目标物数字全息测量的自动聚焦方法 |
CN111044455A (zh) * | 2019-12-27 | 2020-04-21 | 河北工程大学 | 一种数字全息显微成像设备光路共焦装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130342673A1 (en) * | 2011-11-09 | 2013-12-26 | Carl Zeiss Ag | Method and device for focussing a microscope automatically |
-
2016
- 2016-11-08 CN CN201610983732.6A patent/CN106502074B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130342673A1 (en) * | 2011-11-09 | 2013-12-26 | Carl Zeiss Ag | Method and device for focussing a microscope automatically |
Non-Patent Citations (5)
Title |
---|
刘义鹏: "《显微镜自动对焦系统设计与应用》", 《中国优秀博硕士论文全文期刊库(硕士)工程科技Ⅱ辑》 * |
易秋实,张红民,吴萍,等: "《一种光学显微镜自动聚焦系统设计与实现》", 《计算机工程与应用》 * |
王华英: "《像面数字全息显微系统的成像分辨率及成像特点分析》", 《2012中国光学学会全息与光信息处理专业委员会学术年会论文集,》 * |
王大勇: "《数字全息显微成像的理论和实验研究》", 《中国博士学位论文全文数据库信息科技辑》 * |
黄德天: "《基于图像技术的自动调焦方法研究》", 《中国博士学士论文全文数据库信息科技辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109828444A (zh) * | 2019-03-28 | 2019-05-31 | 哈尔滨工业大学 | 一种用于多目标物数字全息测量的自动聚焦方法 |
CN111044455A (zh) * | 2019-12-27 | 2020-04-21 | 河北工程大学 | 一种数字全息显微成像设备光路共焦装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106502074B (zh) | 2019-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tang et al. | Real-time detection of surface deformation and strain in recycled aggregate concrete-filled steel tubular columns via four-ocular vision | |
US7929044B2 (en) | Autofocus searching method | |
CN102840829B (zh) | 基于人工标记的高温物体面内位移场的测量方法 | |
Dai et al. | Strain field estimation based on digital image correlation and radial basis function | |
CN109751964B (zh) | 一种高精度非接触式管径测量方法及装置 | |
Zhang et al. | Accuracy improvement in laser stripe extraction for large-scale triangulation scanning measurement system | |
Thomason et al. | Calibration of a microlens array for a plenoptic camera | |
CN112019751B (zh) | 基于标定信息的自动聚焦方法 | |
CN109141823A (zh) | 一种基于清晰度评价的显微镜系统景深测量装置和方法 | |
CN106502074A (zh) | 一种用于像面数字全息显微测量的自动聚焦方法 | |
Yuan et al. | Dual-biprism-based digital image correlation for defect detection of pipelines | |
JP2007240465A (ja) | 3次元変位ひずみ計測方法及び装置 | |
TWI521295B (zh) | 斜軸式顯微自動對焦系統及方法 | |
Lee et al. | Real-time application of critical dimension measurement of TFT-LCD pattern using a newly proposed 2D image-processing algorithm | |
Solav et al. | Duodic: 3d digital image correlation in Matlab | |
JP6533914B2 (ja) | 計測方法、計測装置、計測プログラム及び計測プログラムを記録した、コンピュータ読み取り可能な記録媒体 | |
Angel-Toro et al. | Speckle decorrelation influence on measurements quality in vortex metrology | |
CN107588741B (zh) | 一种基于莫尔条纹的相机深度变化的测量方法及系统 | |
CN112903121B (zh) | 一种基于交叉迭代自动对焦的波前检测方法 | |
CN106556350A (zh) | 显微镜载玻片曲面高度值的测算方法以及一种显微镜 | |
Chen et al. | In-situ volumetric topography of IC chips for defect detection using infrared confocal measurement with active structured light | |
Nishimura | Characterisation of soil deformation over wide strain ranges in triaxial test with high-precision stereophotogrammetry | |
Zhan | The research of three-dimensional morphology recovery of image sequence based on focusing method | |
CN112857750B (zh) | 一种基于边缘增强的扩展目标波前探测方法 | |
Wu et al. | Topological similarity measurement: a methodology for fast full-field deformation estimation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |