CN110514144B - 一种基于物体光滑表面三维形貌测量系统的测量方法 - Google Patents

一种基于物体光滑表面三维形貌测量系统的测量方法 Download PDF

Info

Publication number
CN110514144B
CN110514144B CN201910855466.2A CN201910855466A CN110514144B CN 110514144 B CN110514144 B CN 110514144B CN 201910855466 A CN201910855466 A CN 201910855466A CN 110514144 B CN110514144 B CN 110514144B
Authority
CN
China
Prior art keywords
modulation degree
dimensional
image
obtaining
smooth surface
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.)
Active
Application number
CN201910855466.2A
Other languages
English (en)
Other versions
CN110514144A (zh
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.)
Sichuan University
Original Assignee
Sichuan University
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 Sichuan University filed Critical Sichuan University
Priority to CN201910855466.2A priority Critical patent/CN110514144B/zh
Publication of CN110514144A publication Critical patent/CN110514144A/zh
Application granted granted Critical
Publication of CN110514144B publication Critical patent/CN110514144B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01B11/25Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object

Abstract

本发明公开了一种三维形貌测量系统的测量系统及基于该测量系统的测量方法。该测量系统包括一维移动载物平台,一维移动载物平台的移动方向上依次设置有显微物镜、分束镜、第一透镜和CCD相机,分束镜镜面与移动方向的夹角为45°,分束镜的一侧依次设置有第二透镜和用于按第一设定间隔时间循环显示具有3幅以上设定相位差的条纹像的投影设备,投影设备、CCD相机和一维移动载物平台均与上位机连接。

Description

一种基于物体光滑表面三维形貌测量系统的测量方法
技术领域
本发明涉及三维形貌测量技术领域,具体涉及一种基于物体光滑表面三维形貌测量系统测量方法。
背景技术
光学三维显微测量技术是光学三维测量的细分,特指利用光学三维测量的方法对微小物体(通常精度在微米或纳米级别)进行三维表面形貌测量的技术。目前比较常见的方法为基于三角测量原理的结构光显微投影技术,但是这种技术由于投影方向与观察方向之间存在一个夹角,导致了遮挡和阴影的问题,并且光学显微测量技术,一般适用于测量漫反射物体。
发明内容
针对现有技术中的上述不足,本发明旨在提供一种能够解决遮挡和阴影问题的一种物体光滑表面三维形貌测量系统及基于该测量系统的测量方法。
为了达到上述发明创造的目的,本发明采用的技术方案为:
提供一种物体光滑表面三维形貌测量系统,其包括一维移动载物平台,一维移动载物平台的移动方向上依次设置有显微物镜、分束镜、第一透镜和CCD相机,分束镜镜面与移动方向的夹角为45°,分束镜的一侧依次设置有第二透镜和用于按第一设定间隔时间循环显示具有3幅以上设定相位差的条纹像的投影设备,投影设备、CCD相机和一维移动载物平台均与上位机连接。
另一方面提供一种基于本方案设计的物体光滑表面三维形貌测量系统的测量方法,其包括:
S1、将待测物体固定在一维移动载物平台上;
S2、采用上位机启动投影设备,并采用CCD相机获取一维移动载物平台在初始位置时待测物体光滑表面对所有条纹像所成条纹虚像;
S3、控制一维移动载物平台按设定步长远离显微物镜,之后采用CCD相机获取待测物体光滑表面对所有条纹像所成条纹虚像,并将移动次数累加一次;
S4、判断移动次数是否大于预设次数,若是,进入步骤S5,否则返回步骤S3;
S5、基于CCD相机获取的图像,采用相移算法得到所有位置下图像的调制度分布,根据每个设定位置下图像的调制度分布得到每个像素点随位置变化的调制度散点图,拟合调制度散点图得到调制度拟合曲线,根据调制度拟合曲线得到每个像素点最大调制度对应的横坐标;
S6、根据横坐标和设定步长得到待测物体光滑表面三维形貌。
进一步地,第一透镜和第二透镜为远心镜头。
进一步地,分束镜呈楔形。
进一步地,基于CCD相机获取的图像,采用相移算法得到所有位置下图像的调制度分布进一步包括:
S5.1、基于待测物体位于同一位置时CCD相机获取的图像,采用相移算法得到该位置下图像的调制度分布,之后进入步骤S5.2;
S5.2、重复步骤S5.1得到所有位置下图像的调制度分布。
进一步地,相移算法的计算公式为:
Figure GDA0002538949210000021
其中,M(x,y)为像素点坐标为(x,y)的调制度,Ii(x,y)为第i张图像中像素点坐标为(x,y)的光强灰度值,L为单次循环中条纹像的副数。
进一步地,Ii的计算公式为:
Ii(x,y)=I0(x,y)+M(x,y)cos(2πfx+δn)
其中,I0(x,y)为背景光强,M(x,y)为条纹像的调制度,f为条纹像的空间频率,δn为设定相位差。
进一步地,调制度拟合曲线为采用一元二次多项式函数拟合调制度散点图得到。
进一步地,待测物体光滑表面三维形貌高度的计算公式如下:
Figure GDA0002538949210000031
其中,h(x,y)为点坐标为(x,y)的高度,△a为设定步长,z0为横坐标的值。
本发明的有益效果为:通过一维移动载物平台、显微物镜、分束镜、第一透镜和CCD相机和投影设备构成的测量系统,其成像光路中没有遮挡和阴影,能够检测物体光滑表面形貌。
附图说明
图1为本发明中测量系统的结构示意图;
图2为具体实施例中分束镜的结构示意图;
图3为具体实施例中某一待测物体的测量结果。
其中,1、一维移动载物平台;2、分束镜;3、CCD相机;4、第一透镜;5、投影设备;6、第二透镜;7、显微物镜。
具体实施方式
下面结合附图,对本发明的具体实施方式做详细说明,以便于本技术领域的技术人员理解本发明。但应该清楚,下文所描述的实施例仅仅是本发明的一部分实施例,而不是全部实施例。在不脱离所附的权利要求限定和确定的本发明的精神和范围内,本领域技术人员在没有做出任何创造性劳动所获得的所有其他实施例,都属于本发明的保护范围。
提供一种基于物体光滑表面三维形貌测量系统,如图1所示,包括一维移动载物平台1,一维移动载物平台1的移动方向上依次设置有显微物镜7、分束镜2、第一透镜4和CCD相机3,分束镜2镜面与移动方向的夹角为45°,分束镜2的一侧依次设置有第二透镜6和用于按第一设定间隔时间循环显示具有3幅以上设定相位差的条纹像的投影设备5,投影设备5、CCD相机3和一维移动载物平台1均与上位机连接。
在测量前一般需要对测量系统进行光路调节使得CCD相机3对待测物体和条纹虚像的成像均较为清晰。具体方法为使投影设备5上条纹像经过第二透镜6及分束镜2反射所成的中间条纹像与CCD相机3经过第一透镜4后的共轭相面重合,也就是该重合面为测量系统调焦位置。
在具体测量时,显微物镜7、分束镜2、第一透镜4和CCD相机3、第二透镜6和投影设备5位置不变。投影设备5上的条纹像经第二透镜6及分束镜2反射为中间条纹像b,进入分束镜2分束面的光线与分束面反射的光线垂直,待测物体光滑表面对该中间条纹像b经显微物镜7和第一透镜4反射到CCD相机3的像平面a,中间条纹像b和CCD像平面a不重合。
另外,在实施时,为避免常规平行分束镜2后表面的反射像正好与分束面的像部分重合而影响成像质量,如图2所示,分束镜2呈楔形。投影设备为DMD数字投影仪。
基于本方案提供的物体光滑表面三维形貌测量系统的测量方法包括:
S1、将待测物体固定在一维移动载物平台1上。
S2、采用上位机启动投影设备5,并采用CCD相机3获取一维移动载物平台1在初始位置时待测物体光滑表面对所有条纹像所成条纹虚像。
S3、控制一维移动载物平台1按设定步长远离显微物镜7,之后采用CCD相机获取待测物体光滑表面对所有条纹像所成条纹虚像,并将移动次数累加一次。
S4、判断移动次数是否大于预设次数,若是,进入步骤S5,否则返回步骤S3;一般初始位置与结束位置之间的距离大于待测物体沿移动方向的厚度。
S5、基于CCD相机3获取的图像,采用相移算法得到所有位置下图像的调制度分布,根据每个设定位置下图像的调制度分布得到每个像素点随位置变化的调制度散点图,拟合调制度散点图得到调制度拟合曲线,根据调制度拟合曲线得到每个像素点最大调制度对应的横坐标。
具体地,基于CCD相机3获取的图像,采用相移算法得到所有位置下图像的调制度分布进一步包括:
S5.1、基于待测物体位于同一位置时CCD相机3获取的图像,采用相移算法得到该位置下图像的调制度分布,之后进入步骤S5.2。
其中,相移算法的计算公式为:
Figure GDA0002538949210000051
其中,M(x,y)为像素点坐标为(x,y)的调制度,Ii(x,y)为第i张图像中像素点坐标为(x,y)的光强灰度值,L为单次循环中条纹像的副数。
上式中,Ii(x,y)的计算公式可以为:
Ii(x,y)=I0(x,y)+M(x,y)cos(2πfx+δn)
其中,I0(x,y)为背景光强,M(x,y)为条纹像的调制度,f为条纹像的空间频率,δn为设定相位差。
S5.2、重复步骤S5.1得到所有位置下图像的调制度分布。
其中,调制度拟合曲线为采用一元二次多项式函数拟合调制度散点图得到。一元二次多项式函数的表达式为Mz=a+bz+cz2,利用调制度散点图中离散数据zi,yi,i=1,2,3...n,n-1为设定次数,横坐标zi和z均表示循环拍摄的次数,也即初始位置时,zi和z均取1,每移动1次,zi和z均加1,纵坐标yi和Mz均表示具体调制度的值,求出拟合函数和离散数据序列的均方误差Pa,b,c,其计算公式为:
Figure GDA0002538949210000061
再通过如下极小值条件,得到系数a,b和c的值,进而求得其极大值对应的横坐标:
Figure GDA0002538949210000062
S6、根据横坐标和设定步长得到待测物体光滑表面三维形貌。
具体地,待测物体光滑表面三维形貌高度的计算公式如下:
Figure GDA0002538949210000063
其中,h(x,y)为点坐标为(x,y)的高度,△a为设定步长,z0为横坐标的值。
实施时,本方案优选第一透镜4和第二透镜6为远心镜头。显微物镜7的放大倍数是10×,数值孔径是0.25。一维载物平台为PI纳米平移台,设定步长为5μm。CCD相机3为IDS的UI-2250SE,其分辨率1600×1200pixels,像元尺寸:4.4×4.4μm;DMD数字投影设备5的分辨率:WXGA1280×800。
其中,针对某一具体待测物体光滑表面三维形貌的测量,单次循环中条纹像数量L为4帧,此时调制度计算公式可改写为
Figure GDA0002538949210000064
共计245个设定位置,每个设定位置拍摄4幅图像。通过该方法得到的其光滑表面三维形貌如图3所示。

Claims (6)

1.一种基于物体光滑表面三维形貌测量系统的测量方法,所述物体光滑表面三维形貌测量系统包括一维移动载物平台(1),所述一维移动载物平台(1)的移动方向上依次设置有显微物镜(7)、分束镜(2)、第一透镜(4)和CCD相机(3),所述分束镜(2)镜面与所述移动方向的夹角为45°,所述分束镜(2)的一侧依次设置有第二透镜(6)和用于按第一设定间隔时间循环显示具有3幅以上设定相位差的条纹像的投影设备(5),所述投影设备(5)、CCD相机(3)和一维移动载物平台(1)均与上位机连接,其特征在于,测量方法包括:
S1、将待测物体固定在一维移动载物平台(1)上;
S2、采用上位机启动投影设备(5),并采用CCD相机(3)获取一维移动载物平台(1)在初始位置时待测物体光滑表面对所有条纹像所成条纹虚像;
S3、控制一维移动载物平台(1)按设定步长远离显微物镜(7),之后采用CCD相机获取待测物体光滑表面对所有条纹像所成条纹虚像,并将移动次数累加一次;
S4、判断移动次数是否大于预设次数,若是,进入步骤S5,否则返回步骤S3;
S5、基于CCD相机(3)获取的图像,采用相移算法得到所有位置下图像的调制度分布,根据每个设定位置下图像的调制度分布得到每个像素点随位置变化的调制度散点图,拟合调制度散点图得到调制度拟合曲线,根据调制度拟合曲线得到每个像素点最大调制度对应的横坐标;
所述调制度拟合曲线为采用一元二次多项式函数拟合调制度散点图得到;
一元二次多项式函数的表达式为M( z) =a+bz+cz2,利用调制度散点图中离散数据zi,yi,i=1,2,3...n,n-1为设定次数,横坐标zi和z均表示循环拍摄的次数,也即初始位置时,zi和z均取1,每移动1次,zi和z均加1,纵坐标yi和M( z) 均表示具体调制度的值,求出拟合函数和离散数据序列的均方误差P(a,b,c),其计算公式为:
Figure FDA0002538949200000021
再通过如下极小值条件,得到系数a,b和c的值,进而求得其极大值对应的横坐标:
Figure FDA0002538949200000022
S6、根据横坐标和设定步长得到待测物体光滑表面三维形貌:
所述待测物体光滑表面三维形貌高度的计算公式如下:
Figure FDA0002538949200000023
其中,h(x,y)为点坐标为(x,y)的高度,△a为设定步长,z0为横坐标的值。
2.根据权利要求1所述的测量方法,其特征在于,所述第一透镜(4)和第二透镜(6)为远心镜头。
3.根据权利要求1所述的测量方法,其特征在于,所述分束镜(2)呈楔形。
4.根据权利要求1所述的测量方法,其特征在于,所述基于CCD相机(3)获取的图像,采用相移算法得到所有位置下图像的调制度分布进一步包括:
S5.1、基于待测物体位于同一位置时CCD相机(3)获取的图像,采用相移算法得到该位置下图像的调制度分布,之后进入步骤S5.2;
S5.2、重复步骤S5.1得到所有位置下图像的调制度分布。
5.根据权利要求4所述的测量方法,其特征在于,所述相移算法的计算公式为:
Figure FDA0002538949200000031
其中,M(x,y)为像素点坐标为(x,y)的调制度,Ii(x,y)为第i张图像中像素点坐标为(x,y)的光强灰度值,L为单次循环中条纹像的副数。
6.根据权利要求5所述的测量方法,其特征在于,所述Ii(x,y)的计算公式为:
Ii(x,y)=I0(x,y)+M(x,y)cos(2πfx+δn)
其中,I0(x,y)为背景光强,M(x,y)为条纹像的调制度,f为条纹像的空间频率,δn为设定相位差。
CN201910855466.2A 2019-09-11 2019-09-11 一种基于物体光滑表面三维形貌测量系统的测量方法 Active CN110514144B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910855466.2A CN110514144B (zh) 2019-09-11 2019-09-11 一种基于物体光滑表面三维形貌测量系统的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910855466.2A CN110514144B (zh) 2019-09-11 2019-09-11 一种基于物体光滑表面三维形貌测量系统的测量方法

Publications (2)

Publication Number Publication Date
CN110514144A CN110514144A (zh) 2019-11-29
CN110514144B true CN110514144B (zh) 2020-10-09

Family

ID=68631888

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910855466.2A Active CN110514144B (zh) 2019-09-11 2019-09-11 一种基于物体光滑表面三维形貌测量系统的测量方法

Country Status (1)

Country Link
CN (1) CN110514144B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6268923B1 (en) * 1999-10-07 2001-07-31 Integral Vision, Inc. Optical method and system for measuring three-dimensional surface topography of an object having a surface contour
CN106197322A (zh) * 2016-09-20 2016-12-07 电子科技大学 一种面结构光三维测量系统
CN108775875A (zh) * 2018-08-31 2018-11-09 苏州大学 一种基于色散光谱编码的微结构形貌测量方法及其装置
CN109341574A (zh) * 2018-09-30 2019-02-15 中国科学院光电技术研究所 一种基于结构光的微纳结构三维形貌高速检测方法
CN109596065A (zh) * 2018-12-07 2019-04-09 中国科学院光电技术研究所 一种基于时域相移算法的高精度微纳三维测量方法
CN110057841A (zh) * 2019-05-05 2019-07-26 电子科技大学 一种基于透射结构光的缺陷检测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6268923B1 (en) * 1999-10-07 2001-07-31 Integral Vision, Inc. Optical method and system for measuring three-dimensional surface topography of an object having a surface contour
CN106197322A (zh) * 2016-09-20 2016-12-07 电子科技大学 一种面结构光三维测量系统
CN108775875A (zh) * 2018-08-31 2018-11-09 苏州大学 一种基于色散光谱编码的微结构形貌测量方法及其装置
CN109341574A (zh) * 2018-09-30 2019-02-15 中国科学院光电技术研究所 一种基于结构光的微纳结构三维形貌高速检测方法
CN109596065A (zh) * 2018-12-07 2019-04-09 中国科学院光电技术研究所 一种基于时域相移算法的高精度微纳三维测量方法
CN110057841A (zh) * 2019-05-05 2019-07-26 电子科技大学 一种基于透射结构光的缺陷检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于相移与调制度比的结构光三维面形垂直测量方法;卢明腾 等;《中国激光》;20190731;第46卷(第7期);0704003-1-9 *

Also Published As

Publication number Publication date
CN110514144A (zh) 2019-11-29

Similar Documents

Publication Publication Date Title
US5307151A (en) Method and apparatus for three-dimensional optical measurement of object surfaces
Chen et al. Adaptive pixel-to-pixel projection intensity adjustment for measuring a shiny surface using orthogonal color fringe pattern projection
JP6363863B2 (ja) 情報処理装置および情報処理方法
WO2016138758A1 (zh) 一种远心成像的三维形貌测量系统标定方法
CN106989689B (zh) 大口径平面光学元件面形的子孔径拼接检测方法
CN107610183B (zh) 一种条纹投影相位高度转换映射模型的标定方法
CN110296667A (zh) 基于线结构光多角度投影的高反射表面三维测量方法
JP2010271199A (ja) 三次元形状測定装置、三次元形状測定方法及びコンピュータプログラム
CN101118155A (zh) 用于感测反光物体的表面形状的方法和系统
JP2001066223A (ja) 反射光品質を決定する方法
JP6112909B2 (ja) シャック・ハルトマンセンサーを用いた形状計測装置、形状計測方法
CN115775303B (zh) 一种基于深度学习与光照模型的高反光物体三维重建方法
Xiaoling et al. Calibration of a fringe projection profilometry system using virtual phase calibrating model planes
CN105528770A (zh) 一种投影仪镜头畸变校正方法
CN110487214A (zh) 一种基于光度立体技术和结构光技术相结合的产品合格率的检测系统及其检测方法
Yang et al. A novel projector ray-model for 3D measurement in fringe projection profilometry
JP2024029135A (ja) 対向配置チャネルを有する三次元センサ
CN110514144B (zh) 一种基于物体光滑表面三维形貌测量系统的测量方法
CN110146032A (zh) 基于光场分布的合成孔径相机标定方法
CN206905715U (zh) 一种可应用于3d曲面玻璃面形检测的测量装置
CN105841620A (zh) 基于投影云纹方法的离面位移测量系统及方法
Yu et al. A three-dimensional measurement system calibration method based on red/blue orthogonal fringe projection
US10613308B2 (en) Method and microscope for measuring and calculating heights on curved surface of microscope slide
CN112361989B (zh) 一种通过点云均匀性考量测量系统标定参数的方法
CN113298882B (zh) 一种用于显微三维形貌测量系统的摄像机标定装置和方法

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
GR01 Patent grant
GR01 Patent grant