CN103727897B - 类镜面表面测试方法 - Google Patents
类镜面表面测试方法 Download PDFInfo
- Publication number
- CN103727897B CN103727897B CN201410026015.5A CN201410026015A CN103727897B CN 103727897 B CN103727897 B CN 103727897B CN 201410026015 A CN201410026015 A CN 201410026015A CN 103727897 B CN103727897 B CN 103727897B
- Authority
- CN
- China
- Prior art keywords
- testee
- class
- ccd
- phase change
- change value
- 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
Links
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开本发明提供了一种类镜面表面测试方法,属于测量技术领域。它解决了现有的类镜面测试方法测试成本高的问题。本发明类镜面表面测试方法其包括:在显示屏上显示余弦条纹,将CCD分别沿被测物体表面X轴和Y轴方向移动移动以获得每一移动位置移动前后CCD摄取的余弦条纹图,比较每一移动位置移动前后的余弦条纹图以获得相位变化值,通过相位变化值得到被测物体类镜面表面X方向和Y方向的表面斜率,通过对斜率求解积分得到被测物体类镜面表面X方向和Y方向的初始高度分布值。本发明不需要设置辅助反射镜与干涉仪,在测试时仅需要一个CCD摄像机即可而且不需要对CCD摄像机进行标定,有效地降低了测量成本,简化了***结构。本发明还提供一种类镜面表面测试装置。
Description
技术领域
本发明涉及测量技术领域,尤其涉及一种类镜面的三维曲面测试方法及装置。
背景技术
光学三维曲面测量技术以其高精度、高效率和非接触性等优点已被广泛应用于高速检测、产品开发、质量控制、反求工程以及CAD/CAM等领域。在工程中,特别是在现代制造业,存在大量表面为类镜面的物体需要测量,例如,汽车工业中的非球面镜、喷漆打蜡后的车身表面、抛光模具表面,建筑陶瓷行业中的瓷砖表面等,该些物体表面均具有类镜面的反射性质。
目前,对于类镜面表面的测量,尤其是对深度类镜面表面的测量非常困难,其主要原因是:一、类镜面表面不具有漫反射特性,对于类镜面表面测量不能采用已成熟的条纹投影轮廓技术;二、现有的类镜面表面测量需要复杂的光学***,比如采用计算全息光学元件对波面进行补偿,一个具体地例子是,在对抛物面镜表面进行测量时,一般采用相移式干涉仪辅助反射镜完成。这种方法具有很高的精度,但在测量过程中需要一块与被检测元件尺寸相同或更大的辅助反射镜,而这样一个高精度反射镜的加工难度大,制造成本高。
发明内容
有鉴于此,本发明的主要目的就是针对现有技术的不足,提出一种成本较低、测试精度高的类镜面表面测试方法及类镜面表面测试。
本发明解决其技术问题采用的技术方案是,提出一种类镜面表面的测试方法,其包括步骤:
S1:在显示屏上显示余弦条纹;将CCD分别沿被测物体类镜面表面X轴方向和Y轴方向移动多次并记录每一移动位置移动前后CCD摄取的余弦条纹图;
S2:比较每一移动位置移动前后的余弦条纹图以获得相位变化值
S3:将相位变化值代入公式:
中,得到被测物体类镜面表面X方向和Y方向的表面斜率;
其中,α表示CCD中轴线与被测物体类镜面表面法线方向的夹角;P表示一副余弦条纹图中相邻条纹间距;L表示被测物体与显示屏的垂直距离;
S4:对X方向和Y方向的斜率分别求解积分,得到被测物体类镜面表面X方向和Y方向的初始高度分布值,X和Y方向初始高度分布值即可表征被测物体类镜面表面的形态。
进一步地,步骤S1中,显示屏与CCD处于同一水平位置。
进一步地,步骤S2中,通过比较每一移动位置移动前后的余弦条纹图的光强分布得到相位变化值其中光强分布I(x,y)与相位变化值的关系式为:
a为零频背景光强度,b为余弦条纹的振幅,P为相邻条纹间距。
进一步地,步骤S1中,CCD分别沿被测物体类镜面表面X轴方向和Y轴方向各移动四次,记录X方向四个移动位置移动前后的余弦条纹图和Y方向四个移动位置移动前后的余弦条纹图。
进一步地,对步骤S4中获得的被测物体类镜面表面X和Y方向的初始高度分布值乘以常数K后得到被测物体类镜面表面的实际高度分布值,其中K=R*δ,R表示CCD影像解析度,δ表示影像像素。
本发明还提供一种类镜面表面测试装置,其包括:
显示装置,用于显示余弦条纹;
CCD,与所述显示装置相同水平高度放置,所述CCD用于沿被测物体类镜面表面X轴方向和Y轴方向移动多次并记录每一移动位置移动前后CCD摄取的余弦条纹图;
处理***,用于比较每一移动位置移动前后的余弦条纹图以获得相位变化值将相位变化值通过公式计算得出被测物体类镜面表面X方向和Y方向的表面斜率,对X方向和Y方向的表面斜率分别求解积分后得到被测物体类镜面表面的X和Y方向的初始高度分布值;其中α表示CCD中轴线与被测物体类镜面表面法线方向的夹角;P表示相邻条纹间距;L表示被测物体与显示屏的垂直距离。
进一步地,所述CCD还根据每一移动位置移动前后的余弦条纹图获得每一移动位置的光强分布I(x,y),所述处理***根据该光强分布I(x,y)获得相位变化值所述光强分布I(x,y)与相位变化值的关系式为:
a为背景光强度,b为振幅,P为相邻条纹间距。
进一步地,所述处理***还将被测物体类镜面表面X和Y方向的初始高度分布值乘以常数K后得到被测物体类镜面表面X和Y方向的实际高度分布值,其中K=R*δ,R表示CCD影像解析度,δ表示影像像素。
本发明通过利用CCD在被测物体类镜面表面多个位置移动以获得移动前后的余弦条纹图,通过该余弦条纹图结合相移法得到被测物体表面各点的斜率信息并通过积分同时得到被测物体表面各点的坐标信息从而能够精确地得出被测物体类镜面表面的参数值。本发明不需要设置辅助反射镜与干涉仪,在测试时仅需要一个CCD摄像机即可而且不需要对CCD摄像机进行标定,有效地降低了测量成本,简化了***结构,能更加方便灵活地测量各类类镜面表面。并且本发明有较好的抗噪性能,在噪声较大的情况下,仍能对类镜面表面进行测量。
附图说明
图1为本发明一较佳实施例中类镜面表面测试方法的流程图;
图2为由CCD获取被测物体类镜面表面条纹偏移量的计算方法图;
图3为本发明一较佳实施例中类镜面表面测试装置的示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
请参照图1,图1为本发明一较佳实施例中类镜面表面测试方法的流程图。图1中,类镜面表面测试方法包括步骤:
S1:在显示屏上显示余弦条纹;将CCD分别沿被测物体类镜面表面X轴方向和Y轴方向移动多次并记录每一移动位置移动前后CCD摄取的余弦条纹图。
由于被测物体表面为类镜面,也就是为非平面,因此被测物体表面各处的斜率和高度存在偏差。将CCD分别沿被测物体类镜面表面X轴方向和Y轴方向移动并每一移动位置移动前后CCD摄取的余弦条纹图就可得知被测物体表面各处因为斜率变化和高度变化带来的条纹偏移量。通过分析条纹偏移量反推出被测物体类镜面表面的各项参数。
下面通过分析X轴方向被测物体表面斜率和高度的变化带来的条纹偏移量的计算。同理可以得出Y轴方向被测物体表面斜率和高度的变化带来的条纹偏移量。
首先描述由于被测物体类镜面表面高度变化带来的条纹偏移量的计算。
由CCD获取被测物体类镜面表面条纹偏移量的具体方法如图2所示。图2中,CCD10以α角度获取被测物体表面的影像,α即CCD中轴线与被测物体类镜面表面法线方向的夹角;CCD移动前摄取的余弦条纹图对应的被测物体表面设为第一镜面11,CCD移动后摄取的余弦条纹图对应的被测物体表面设为第二镜面12。第一镜面11和第二镜面12之间的高度差设为h,显示屏与被测物体表面的垂直距离设为L。
当CCD以固定角度α拍摄第一镜面11,第一镜面11对应于显示屏余弦条纹图中的A点,当CCD移动后以仍以固定角度α拍摄第二镜面12时,第二镜面12对应于显示屏余弦条纹图中的B点,即由于第一镜面11和第二镜面12之间的高度差h使得余弦条纹具有条纹偏离量dAB,条纹偏离量dAB与高度差h的几何关系为:
dAB=2h·tanα (1)
接下来描述由于被测物体类镜面表面斜率变化带来的条纹偏移量的计算。
请继续参照图2,当CCD移动后以仍以固定角度α拍摄第二镜面12时,此时由于第二镜面12还存在斜率变化θ,此斜率变化在反射时会造成2θ的角度偏移,此角度偏移会使第二镜面12对应于余弦条纹图中的B点偏移至C点,其偏移距离为dBC,dAB可由如下关系式表示:
dBC=(L-h)[tan(α+2θ)-tan(α)] (2)
由式(1)和式(2)可知,由于高度差h以及斜率变化θ在X方向所造成的的总偏移量dX为:
dx=(L-h)[tan(α+2θ)-tan(α)]-2htan(α) (3)
若显示屏与被测物体表面的垂直距离L远大于被测物体表面的高度差,则可将式(3)简化为:
dx=L[tan(α+2θ)-tan(α)] (4)
利用合角公式将式(4)整理后可以得到:
若θ在2.5度以内,且tan(α)·tan(2θ)<<1,因此可以将式(5)简化为
dx=2Lθ·sec2(α) (6)
S2:比较每一移动位置移动前后的余弦条纹图以获得相位变化值
相位变化值可通过计算机程序直接获取。优选地,相位变化值通过比较每一移动位置移动前后的余弦条纹图的光强分布得到。
假如物体表面为完美平面,则条纹影像的光强分布为一个余弦函数,因此由CCD摄取到的条纹影像的光强分布为
其中,a为背景光强度,b为振幅,P为相邻条纹间距。
若物体表面有斜率变化时,此斜率变化会造成条纹的变形,变形后的条纹的光强分布可以表示为:
即条纹变形所产生的相位变化值。
S3:将相位变化值代入公式:
中,得到被测物体类镜面表面X方向和Y方向的表面斜率;
其中,α表示CCD中轴线与被测物体类镜面表面法线方向的夹角;P表示相邻条纹间距;L表示被测物体与显示屏的垂直距离。
通过如下方式得出:
相位变化值的单位为为弧度(rad),其相位变化值与长度单位之间的关系为:
将式(6)代入式(9)可得到:
将式(10)移项整理后即可得到
S4:对X方向和Y方向的斜率分别求解积分,得到被测物体类镜面表面X方向和Y方向的初始高度分布值,X和Y方向的初始高度分布值即可表征被测物体类镜面表面的形态。
具体地,对X方向斜率求解积分的公式如下:
h(x,y)=∫θ(x,y)dx
同理可得出对Y方向斜率求解积分的公式。
另外,由于CCD存在像素和解析度,由计算机得出的X和Y方向的初始高度分布值虽然已可以表征被测物体类镜面表面的形态,但是当需要获得被测物体类镜面表面的实际形态的参数值时还需要考虑CCD的像素和解析度。将步骤S4中获得的初始高度分布值乘以常数K即可获得被测物体类镜面表面的实际形态的参数值,其中K=R*δ,R表示CCD影像解析度,δ表示影像像素。
本发明还提供一种类镜面表面的测试装置,装置示意图如图3所示。图3中,类镜面表面的测试装置包括:
显示装置100,用于显示余弦条纹;
CCD200,与所述显示装置100相同水平高度放置,所述CCD用于沿被测物体300类镜面表面X轴方向和Y轴方向移动多次并记录每一移动位置移动前后CCD摄取的余弦条纹图;
处理***400,用于比较每一移动位置移动前后的余弦条纹图以获得相位变化值将相位变化值通过公式计算得出被测物体类镜面表面X方向和Y方向的表面斜率,对X方向和Y方向的表面斜率分别求解积分后得到被测物体类镜面表面的X和Y方向的初始高度分布值;其中α表示CCD中轴线与被测物体类镜面表面法线方向的夹角;P表示相邻条纹间距;L表示被测物体与显示屏的垂直距离。
优选地,所述CCD200还根据每一移动位置移动前后的余弦条纹图获得每一移动位置的光强分布I(x,y),所述处理***400根据该光强分布I(x,y)获得相位变化值所述光强分布I(x,y)与相位变化值的关系式为:
a为背景光强度,b为振幅,P为相邻条纹间距。
优选地,所述处理***还将被测物体类镜面表面X和Y方向的初始高度分布值乘以常数K后得到被测物体类镜面表面X和Y方向的实际高度分布值,其中K=R*δ,R表示CCD影像解析度,δ表示影像像素。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (8)
1.一种类镜面表面的测试方法,其特征在于,包括步骤:
S1:在显示屏上显示余弦条纹;将CCD分别沿被测物体类镜面表面X轴方向和Y轴方向移动多次并记录每一移动位置移动前后CCD摄取的余弦条纹图;
S2:比较每一移动位置移动前后的余弦条纹图以获得相位变化值
S3:将相位变化值代入公式:中,得到被测物体类镜面表面X方向和Y方向的表面斜率;
其中,α表示CCD中轴线与被测物体类镜面表面法线方向的夹角;P表示相邻条纹间距;L表示被测物体与显示屏的垂直距离;
S4:对X方向和Y方向的斜率分别求解积分,得到被测物体类镜面表面X方向和Y方向的初始高度分布值,X和Y方向初始高度分布值即可表征被测物体类镜面表面的形态。
2.如权利要求1所述的测试方法,其特征在于,步骤S1中,显示屏与CCD处于同一水平位置。
3.如权利要求1所述的测试方法,其特征在于,步骤S2中,通过比较每一移动位置移动前后的余弦条纹图的光强分布得到相位变化值其中光强分布I(x,y)与相位变化值的关系式为:
a为背景光强度,b为振幅,P为相邻条纹间距。
4.如权利要求1所述的测试方法,其特征在于,步骤S1中,CCD分别沿被测物体类镜面表面X轴方向和Y轴方向各移动四次,记录X方向四个移动位置移动前后的余弦条纹图和Y方向四个移动位置移动前后的余弦条纹图。
5.如权利要求1至4任一项所述的测试方法,其特征在于,对步骤S4中获得的被测物体类镜面表面X和Y方向的初始高度分布值乘以常数K后得到被测物体类镜面表面的实际高度分布值,其中K=R*δ,R表示CCD影像解析度,δ表示影像像素。
6.一种类镜面表面测试装置,其特征在于,包括:
显示装置,用于显示余弦条纹;
CCD,与所述显示装置相同水平高度放置,所述CCD用于沿被测物体类镜面表面X轴方向和Y轴方向移动多次并记录每一移动位置移动前后CCD摄取的余弦条纹图;
处理***,用于比较每一移动位置移动前后的余弦条纹图以获得相位变化值将相位变化值通过公式计算得出被测物体类镜面表面X方向和Y方向的表面斜率,对X方向和Y方向的表面斜率分别求解积分后得到被测物体类镜面表面的X和Y方向的初始高度分布值;其中α表示CCD中轴线与被测物体类镜面表面法线方向的夹角;P表示相邻条纹间距;L表示被测物体与显示屏的垂直距离。
7.如权利要求6所述的测试装置,其特征在于,所述CCD还根据每一移动位置移动前后的余弦条纹图获得每一移动位置的光强分布I(x,y),所述处理***根据该光强分布I(x,y)获得相位变化值所述光强分布I(x,y)与相位变化值的关系式为:a为背景光强度,b为振幅,P为相邻条纹间距。
8.如权利要求6或7所述的测试装置,其特征在于,所述处理***还将被测物体类镜面表面X和Y方向的初始高度分布值乘以常数K后得到被测物体类镜面表面X和Y方向的实际高度分布值,其中K=R*δ,R表示CCD影像解析度,δ表示影像像素。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410026015.5A CN103727897B (zh) | 2014-01-21 | 2014-01-21 | 类镜面表面测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410026015.5A CN103727897B (zh) | 2014-01-21 | 2014-01-21 | 类镜面表面测试方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103727897A CN103727897A (zh) | 2014-04-16 |
CN103727897B true CN103727897B (zh) | 2016-11-16 |
Family
ID=50452077
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410026015.5A Active CN103727897B (zh) | 2014-01-21 | 2014-01-21 | 类镜面表面测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103727897B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109405770B (zh) * | 2018-11-02 | 2020-12-25 | 上海华力微电子有限公司 | 一种通过光阻膜厚监控涂胶显影机台涂胶槽水平程度的方法 |
CN110910506B (zh) * | 2019-12-03 | 2021-02-09 | 江苏集萃华科智能装备科技有限公司 | 基于法线检测的三维重建方法、装置、检测装置及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002039726A (ja) * | 2000-07-28 | 2002-02-06 | Yamatake Corp | 表面形状測定装置 |
CN101893766A (zh) * | 2010-06-23 | 2010-11-24 | 北京航空航天大学 | 一种变入射光角度实现相移的方法及其装置 |
CN102183213A (zh) * | 2011-03-02 | 2011-09-14 | 中国科学院光电技术研究所 | 一种基于相位测量偏折术的非球面镜检测方法 |
-
2014
- 2014-01-21 CN CN201410026015.5A patent/CN103727897B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002039726A (ja) * | 2000-07-28 | 2002-02-06 | Yamatake Corp | 表面形状測定装置 |
CN101893766A (zh) * | 2010-06-23 | 2010-11-24 | 北京航空航天大学 | 一种变入射光角度实现相移的方法及其装置 |
CN102183213A (zh) * | 2011-03-02 | 2011-09-14 | 中国科学院光电技术研究所 | 一种基于相位测量偏折术的非球面镜检测方法 |
Non-Patent Citations (1)
Title |
---|
基于条纹反射的类镜面三维面形测量方法;刘元坤,苏显渝;《光学学报》;20061130;第26卷(第11期);1636-1640 * |
Also Published As
Publication number | Publication date |
---|---|
CN103727897A (zh) | 2014-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Three-dimensional shape measurements of specular objects using phase-measuring deflectometry | |
CN102183214B (zh) | 一种大口径非球面镜结构光检测方法 | |
CN102183213B (zh) | 一种基于相位测量偏折术的非球面镜检测方法 | |
CN105783775A (zh) | 一种镜面及类镜面物体表面形貌测量装置与方法 | |
CN101865675B (zh) | 一种三维视觉检测圆结构光标定点获取方法 | |
CN101655360A (zh) | 一种可绝对位相展开的十六步双频光栅相移轮廓术 | |
CN109059806A (zh) | 一种基于红外条纹的镜面物体三维面形测量装置及方法 | |
CN104111039A (zh) | 用于任意摆放条纹投影三维测量***的标定方法 | |
CN102494634A (zh) | 一种基于条纹反射的离轴非球面镜检测方法 | |
CN106949851A (zh) | 一种基于支持向量机的线结构光视觉传感器标定方法 | |
TWI471522B (zh) | 使用陰影雲紋法線上即時檢測表面形貌與平面外變形之方法與系統 | |
CN109631798A (zh) | 一种基于π相移方法的三维面形垂直测量方法 | |
CN104614342B (zh) | 一种高温气流扰动下空气折射率三维重构测量方法 | |
CN107063122B (zh) | 光学非球面面形的检测方法及其装置 | |
CN104833319A (zh) | 透镜中心偏差测量***及测量方法 | |
KR101097716B1 (ko) | 3차원 형상 측정방법 | |
Zou et al. | High-accuracy calibration of line-structured light vision sensors using a plane mirror | |
CN105953749A (zh) | 一种光学三维形貌测量方法 | |
CN102410819A (zh) | 一种测量膜基反射镜三维面形的方法 | |
CN106546193B (zh) | 一种高反射物体表面三维测量方法和*** | |
CN103925890B (zh) | 一种基于光束畸变的三维角度测量*** | |
CN105491315B (zh) | 一种投影仪伽马校正方法 | |
CN103727897B (zh) | 类镜面表面测试方法 | |
CN110146032B (zh) | 基于光场分布的合成孔径相机标定方法 | |
CN109855559A (zh) | 一种全空间标定***及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address |
Address after: 311258 No. 1398, Wenbin Road, Wenyan street, Xiaoshan District, Hangzhou, Zhejiang. Patentee after: First facing the three dimensional Polytron Technologies Inc Address before: 311215 Zhejiang Hangzhou Xiaoshan economic and Technological Development Zone Construction No. 66 Huarui Center 1 building 18 floor. Patentee before: Hangzhou Shining 3D Technology Co., Ltd. |
|
CP03 | Change of name, title or address |