CN110260817A - 基于虚拟标志点的复杂曲面偏折测量自定位方法 - Google Patents

基于虚拟标志点的复杂曲面偏折测量自定位方法 Download PDF

Info

Publication number
CN110260817A
CN110260817A CN201910568780.2A CN201910568780A CN110260817A CN 110260817 A CN110260817 A CN 110260817A CN 201910568780 A CN201910568780 A CN 201910568780A CN 110260817 A CN110260817 A CN 110260817A
Authority
CN
China
Prior art keywords
point
coordinate
actual measurement
measurement
index point
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
Application number
CN201910568780.2A
Other languages
English (en)
Other versions
CN110260817B (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.)
Fudan University
Original Assignee
Fudan 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 Fudan University filed Critical Fudan University
Priority to CN201910568780.2A priority Critical patent/CN110260817B/zh
Publication of CN110260817A publication Critical patent/CN110260817A/zh
Application granted granted Critical
Publication of CN110260817B publication Critical patent/CN110260817B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明属于精密测量技术领域，具体为一种基于虚拟标志点的复杂曲面偏折测量自定位方法。本发明方法步骤为：将被测元件口径中心点作为虚拟标志点；首先探测此中心点的高度坐标，计算实测标志点的世界坐标；根据偏折测量图像的相位信息计算实测标志点对应的投影屏像素坐标，以及实测标志点的法向量；利用反向追迹方法，迭代估算球的半径和对应球心坐标；最后得到测量点的坐标，以确定被测工件面形的横向坐标。本发明可有效计算被测面形最佳拟合球面参数，提供可靠的初始测量坐标，避免传统方法中繁琐困难的位置标定工作，对于提高偏折测量的效率与通用性有重要意义。

Description

基于虚拟标志点的复杂曲面偏折测量自定位方法
技术领域
本发明属于光学工程技术领域，具体为一种复杂曲面偏折测量自定位方法。
背景技术
高质量的光学元件越来越多地应用于各种领域，包括精密工程、航空航天等。光学元件的面形对其成像功能的实现至关重要。对于光滑的镜面，相位测量偏折术有着诸多优异的特性。其中，单目偏折术的测量系统简单，动态范围大，抗干扰能力强，可用于复杂曲面的测量，近年来得到广泛关注。其原理是在显示器上产生规则条纹，经被测表面反射后条纹发生变形，采用CCD相机拍摄变形图样，由已知被测工件面型参数估算被测点坐标，由几何关系推导可以计算出被测面形的表面梯度分布，再通过积分得到面形高度。
在实际测量过程中，待测工件的坐标系需要与工作台的世界坐标系统一，以便准确地估计被测点坐标。常用的方法是利用三坐标测量机或者激光跟踪器将工作台和待测工件进行坐标系标定，这种方法不仅标定工作繁琐困难，而且成本高昂。因此需要一种更简便的高精度工件定位的标定方法。
发明内容
本发明的目的在于提供一种复杂曲面偏折测量自定位方法，以实现对待测工件精确定位、对初始测量点的坐标的准确地估计。
由于大部分光学反射精密为凹面，承担聚焦或光束准直功能，局部形状可近似用球模型拟合；在偏折测量中，对于未知面形的回转光学元件，可以用球模型作为预设面形计算被测点坐标，以供后续面型重构。
本发明提供的复杂曲面偏折测量自定位方法，是基于虚拟标志点的技术的，具体步骤如下：
(a)利用高度计测量光学元件的回转中心(即虚拟标志点)的高度h；
(b)将光学元件水平放置在测量平台上，进行偏折测量；
(c)在获取的图像上提取测量区域的重心，即为实测标志点(Xm,Ym,Zm)，根据系统标定参数，得到该实测标志点的反射光线的空间直线表达式：
其中，(X0,Y0,Z0)为相机光心坐标，(am,bm,cm)是实测反射光线的单位方向向量；
(d)将回转中心的高度h作为实测标志点的高度Z代入公式(1)中的直线表达式中，计算实测标志点的Xm、Ym坐标：
(e)根据实际测量的图像信息，计算实测标志点对应屏幕像素坐标，也就是点光源坐标(XS,YS,ZS)；
(f)由点光源坐标(XS,YS,ZS)、实测标志点坐标(Xm,Ym,Zm)、反射光线三组参数，得到实测标志点法向量
(g)由实测标志点(Xm,Ym,Zm)和法向量得到一条空间直线：
其中，Nx、Ny、Nz分别是法向量的分量；
(h)将被测表面近似拟合为球面，则球心在实测标志点与该点法向量组成的空间直线上；球的半径利用被测点反向追迹的方法估算得到。
具体地，利用被测点反向追迹的方法估算球的半径，迭代过程为：
(i)将相机光心到标志点的距离设为初始半径R，球心坐标(Xc,Yc,Zc)满足以下条件：
(ii)由球心坐标、半径R和图像中每个测量点的像素坐标得到初始测量点坐标(Xi,Yi,Zi)：
其中，(ai,bi,ci)是每个测量点处反射光线的单位方向向量；
于是可以得到该点的法向量以及入射光线的方向向量
其中，·是内积符号；
(iii)根据每个测量点的坐标和入射光线进行反向追迹，并计算追迹的屏幕坐标，得到反向追迹误差；
(iv)判断反向追迹误差是否达到最小误差,不满足，则返回步骤(i)；满足则终止迭代过程，输出所有估计的被测点坐标，并进行后续的面型重构。
进一步地，由所拟合的球面球心及半径确定被测工件面形的位置与基本形状。
本发明可有效计算被测面形最佳拟合球面参数，提供可靠的初始测量坐标，避免传统方法中繁琐困难的位置标定工作，对于提高偏折测量的效率与通用性有重要意义。
附图说明
图1为球面标志点估算图示。
图2为球面镜测量结果图示。
图3为球面镜测量误差图示。
具体实施方式
下面通过实施例结合附图进一步说明本发明。
实施例1：利用高度计测量口径为70mm、半径为1000mm的球面镜的回转中心(虚拟标志点)的高度h，测量结果为12mm；将球面镜水平放置在测量平台的任意位置处，并进行偏折测量；根据侧轮廓的圆属性，在图像上提取测量区域的重心，即为实测标志点，根据系统标定参数，可以得到该实测标志点的反射光线的空间直线表达式；将回转中心的高度作为实测标志点的高度Z代入上述的直线表达式中，计算实测标志点的X、Y坐标分别为1.63mm和0.26mm；如图1所示，估算的轮廓中心和实际的回转中心有一定的偏移，得到反射光线为45°，球面镜矢高为1mm，球的半径为100mm，则横向偏移量为1mm，而实测标志点和虚拟标志点的纵向偏差仅为5um，所以对于大部分被测元件，实验中估算的实测标志点坐标误差不会超过误差容许范围。
根据本发明的方法，利用被测点反向追迹的方法迭代估算球的半径为1020mm，球心坐标为(1.7mm，0mm，1032.1mm)，计算所有估计的被测点坐标，并进行后续的面型重构，重构结果如图2所示。测量误差RMS为85.6nm，误差图如图3所示。

Claims (3)

1.一种基于虚拟标志点的复杂曲面偏折测量自定位方法，其特征在于，具体步骤如下：
(a)利用高度计测量光学元件的回转中心即虚拟标志点的高度h；
(b)将光学元件水平放置在测量平台上，进行偏折测量；
(c)在获取的图像上提取测量区域的重心，即为实测标志点(Xm,Ym,Zm)，根据系统标定参数，得到该实测标志点的反射光线的空间直线表达式:
其中，(X0,Y0,Z0)为相机光心坐标，(am,bm,cm)是实测反射光线的单位方向向量；
(d)将回转中心的高度h作为实测标志点的高度Z代入公式(1)中的直线表达式中，计算实测标志点的Xm、Ym坐标：
(e)根据实际测量的图像信息，计算实测标志点对应屏幕像素坐标，也就是点光源坐标(XS,YS,ZS)；
(f)由点光源坐标(XS,YS,ZS)、实测标志点坐标(Xm,Ym,Zm)、反射光线三组参数，得到实测标志点法向量
(g)由实测标志点(Xm,Ym,Zm)和法向量得到一条空间直线:
其中，Nx、Ny、Nz分别是法向量的分量；
(h)将被测表面近似拟合为球面，球心在实测标志点与该点法向量组成的空间直线上；球的半径利用被测点反向追迹的方法估算得到。
2.根据权利要求1所述的复杂曲面偏折测量自定位方法，其特征在于，利用被测点反向追迹的方法估算球的半径，迭代过程为:
(i)将相机光心到标志点的距离设为初始半径R，球心坐标(Xc,Yc,Zc)满足以下条件：
(ii)由球心坐标、半径R和图像中每个测量点的像素坐标得到初始测量点坐标(Xi,Yi,Zi):
其中，(ai,bi,ci)是每个测量点处反射光线的单位方向向量；
于是得到该点的法向量以及入射光线的方向向量
其中，·是内积符号；
(iii)根据每个测量点的坐标和入射光线进行反向追迹，并计算追迹的屏幕坐标，得到反向追迹误差；
(iv)判断反向追迹误差是否达到最小误差,不满足，则返回步骤(i)；满足则终止迭代过程，输出所有估计的被测点坐标，并进行后续的面型重构。
3.根据权利要求2所述的复杂曲面偏折测量自定位方法，其特征在于，由所拟合的球面球心及半径确定被测工件面形的位置与基本形状。
CN201910568780.2A 2019-06-27 2019-06-27 基于虚拟标志点的复杂曲面偏折测量自定位方法 Active CN110260817B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910568780.2A CN110260817B (zh) 2019-06-27 2019-06-27 基于虚拟标志点的复杂曲面偏折测量自定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910568780.2A CN110260817B (zh) 2019-06-27 2019-06-27 基于虚拟标志点的复杂曲面偏折测量自定位方法

Publications (2)

Publication Number Publication Date
CN110260817A true CN110260817A (zh) 2019-09-20
CN110260817B CN110260817B (zh) 2020-11-20

Family

ID=67922497

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910568780.2A Active CN110260817B (zh) 2019-06-27 2019-06-27 基于虚拟标志点的复杂曲面偏折测量自定位方法

Country Status (1)

Country Link
CN (1) CN110260817B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111536873A (zh) * 2020-05-19 2020-08-14 复旦大学 一种采用辅助相机对单目偏折测量工件进行定位和测量的方法
CN114964035A (zh) * 2022-04-08 2022-08-30 中国科学院上海光学精密机械研究所 基于空变球模型自适应单目偏折的工件面形重构方法
CN118654584A (zh) * 2024-08-15 2024-09-17 恒迈光学精密机械(杭州)有限公司 一种自由曲面透明件的光学等厚测试方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07260443A (ja) * 1994-03-24 1995-10-13 Nkk Corp ３次元曲面形状の測定方法及び装置
JPH08178635A (ja) * 1994-12-21 1996-07-12 C K S Chiyuuki:Kk 背板材の形状寸法測定方法
CN102183213A (zh) * 2011-03-02 2011-09-14 中国科学院光电技术研究所 一种基于相位测量偏折术的非球面镜检测方法
CN103389072A (zh) * 2013-07-22 2013-11-13 北京信息科技大学 一种基于直线拟合的像点定位精度评估方法
CN106017863A (zh) * 2016-05-11 2016-10-12 上海应用技术学院 检测非球面的相位测量偏折方法
CN108507492A (zh) * 2018-05-14 2018-09-07 中国计量大学 一种平曲面透射元件表面的高精度大动态范围测量方法及测量系统
CN109373931A (zh) * 2018-12-14 2019-02-22 上海晶电新能源有限公司 一种太阳能热发电用光学设备反射面面形检测系统及方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07260443A (ja) * 1994-03-24 1995-10-13 Nkk Corp ３次元曲面形状の測定方法及び装置
JPH08178635A (ja) * 1994-12-21 1996-07-12 C K S Chiyuuki:Kk 背板材の形状寸法測定方法
CN102183213A (zh) * 2011-03-02 2011-09-14 中国科学院光电技术研究所 一种基于相位测量偏折术的非球面镜检测方法
CN103389072A (zh) * 2013-07-22 2013-11-13 北京信息科技大学 一种基于直线拟合的像点定位精度评估方法
CN106017863A (zh) * 2016-05-11 2016-10-12 上海应用技术学院 检测非球面的相位测量偏折方法
CN108507492A (zh) * 2018-05-14 2018-09-07 中国计量大学 一种平曲面透射元件表面的高精度大动态范围测量方法及测量系统
CN109373931A (zh) * 2018-12-14 2019-02-22 上海晶电新能源有限公司 一种太阳能热发电用光学设备反射面面形检测系统及方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111536873A (zh) * 2020-05-19 2020-08-14 复旦大学 一种采用辅助相机对单目偏折测量工件进行定位和测量的方法
CN114964035A (zh) * 2022-04-08 2022-08-30 中国科学院上海光学精密机械研究所 基于空变球模型自适应单目偏折的工件面形重构方法
CN118654584A (zh) * 2024-08-15 2024-09-17 恒迈光学精密机械(杭州)有限公司 一种自由曲面透明件的光学等厚测试方法

Also Published As

Publication number Publication date
CN110260817B (zh) 2020-11-20

Similar Documents

Publication Publication Date Title
US10690492B2 (en) Structural light parameter calibration device and method based on front-coating plane mirror
CN109163657B (zh) 一种基于双目视觉三维重建的圆形目标位姿检测方法
CN104990515B (zh) 大型物体三维形状测量系统及其测量方法
JP5270670B2 (ja) ２次元画像による３次元組立て検査
He et al. Accurate calibration method for blade 3D shape metrology system integrated by fringe projection profilometry and conoscopic holography
CN108709499A (zh) 一种结构光视觉传感器及其快速标定方法
CN108489421A (zh) 一种条纹投影检测平面元件面形方法及装置
CN109341546A (zh) 一种点激光位移传感器在任意安装位姿下的光束标定方法
CN109307480B (zh) 一种透射元件多表面面形检测方法
CN110260817A (zh) 基于虚拟标志点的复杂曲面偏折测量自定位方法
CN111366079B (zh) 一种快速确定偏折测量系统中各部件几何位置的标定方法
CN108195314B (zh) 基于多视场拼接的反射式条纹三维面形测量方法
CN111207670A (zh) 一种线结构光标定装置及方法
CN103712572A (zh) 结构光源与相机结合的物体轮廓三维坐标测量装置
CN110428471B (zh) 一种针对光学自由曲面子孔径偏折测量的精确自定位方法
Jia et al. A field measurement method for large objects based on a multi-view stereo vision system
CN110146032B (zh) 基于光场分布的合成孔径相机标定方法
CN103697811A (zh) 一种相机与结构光源结合获取物体轮廓三维坐标的方法
CN110211175B (zh) 准直激光器光束空间位姿标定方法
CN113916128A (zh) 一种基于光笔式视觉测量系统的提高精度的方法
Ye et al. A calibration trilogy of monocular-vision-based aircraft boresight system
CN105758339A (zh) 基于几何误差修正技术的光轴与物面垂直度检测方法
CN115164793A (zh) 微纳米级复合测量的坐标统一标定器及标定方法
CN105841636A (zh) 基于直线运动部件误差补偿的光轴与物面垂直度检测方法
CN115289997A (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