CN110428471B - 一种针对光学自由曲面子孔径偏折测量的精确自定位方法 - Google Patents
一种针对光学自由曲面子孔径偏折测量的精确自定位方法 Download PDFInfo
- Publication number
- CN110428471B CN110428471B CN201910606141.0A CN201910606141A CN110428471B CN 110428471 B CN110428471 B CN 110428471B CN 201910606141 A CN201910606141 A CN 201910606141A CN 110428471 B CN110428471 B CN 110428471B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- screen
- world coordinate
- measurement
- sub
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
Abstract
本发明涉及一种针对光学自由曲面子孔径偏折测量的精确自定位方法,包括以下步骤:1)计算名义面形f上的测量点P的坐标并获取测量点P对应的法向量n;2)根据法向量n逆向追迹屏幕坐标S;3)将回转台从位置A旋转一定角度至位置B,并测量两个位置下的偏折图像,构建最小二乘代价方程;4)求解最小二乘代价方程,即得到实测子孔径表面的位置。与现有技术相比,本发明具有标定简单准确等优点。
Description
技术领域
本发明涉及精密测量技术领域,尤其是涉及一种针对光学自由曲面子孔径偏折测量的精确自定位方法。
背景技术
高质量的光学元件越来越多地应用于各种领域,包括精密工程、航空航天等。光学元件的面形对其成像功能的实现至关重要。对于光滑的镜面,相位测量偏折术有着诸多优异的特性。其中,单目偏折术的测量系统简单,动态范围大,抗干扰能力强,可用于复杂曲面的测量,近年来得到广泛关注。其原理是在显示器上产生规则条纹,经被测表面反射后条纹发生变形,采用CCD相机拍摄变形图样,由已知被测工件面型参数估算被测点坐标,由几何关系推导可以计算出被测面形的表面梯度分布,再通过积分得到面形高度。
在实际测量过程中,待测工件的坐标系需要与工作台的世界坐标系统一,以便准确地估计被测点坐标。常用的方法是利用三坐标测量机或者激光跟踪器将工作台和待测工件进行坐标系标定,这种方法不仅标定工作繁琐困难,而且成本高昂。因此需要一种更简便的高精度工件定位的标定方法。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种针对光学自由曲面子孔径偏折测量的精确自定位方法。
本发明的目的可以通过以下技术方案来实现:
一种针对光学自由曲面子孔径偏折测量的精确自定位方法,包括以下步骤:
1)计算名义面形f上的测量点P的坐标并获取测量点P对应的法向量n;
2)根据法向量n逆向追迹屏幕坐标S;
3)将回转台从位置A旋转一定角度至位置B,并测量两个位置下的偏折图像,构建最小二乘代价方程;
4)求解最小二乘代价方程,即得到实测子孔径表面的位置。
所述的步骤1)中,测量点P的坐标[xp yp f(xp,yp)]T计算式为:
[xp yp f(xp,yp)]T=τr+(Rw_m·[xc yc zc]T+Tw_m)
r=Rw_m·([xc yc zc]T-Tc_w)
其中,τ为P点到相机光心的距离,r为反射光线的单位方向向量,Rw_m为从世界坐标系到模型坐标系的旋转矩阵,[xc yc zc]T为世界坐标系中相机像素C的坐标,Tw_m为从世界坐标系到模型坐标系的平移向量,Tc_w为从相机坐标系到世界坐标系的平移向量。
所述的步骤2)具体为:
根据法向量n和反射定律计算入射光线i,根据入射光线i和屏幕所在平面相交的参数ε获取屏幕S的坐标,以及对应的屏幕原点和从屏幕坐标系到世界坐标系的平移向量Ts_w。
所述的入射光线i的计算式为:
i=r-2(nT·r)n。
所述的屏幕S的坐标与从屏幕坐标系到世界坐标系的平移向量Ts_w的关系式为:
ns=Rw_mRs_w·[0 0 1]T
其中,[xs ys zs]T为屏幕S的坐标,ns为屏幕所在平面的法向量,Rs_w为从屏幕坐标系到世界坐标系的旋转矩阵。
所述的步骤3)中,最小二乘代价方程的表达式为:
x=[Rm_w,Tm_w]
其中,x为待解的模型位置参数,包括从模型坐标系到世界坐标系的旋转矩阵Rm_w和从模型坐标系到世界坐标系的平移向量Tm_w,为实际屏幕坐标,SA和SB是反向追迹的屏幕坐标,l为相机中的测量点总数,i为相机中的测量点数。
所述的步骤4)中,采用Levenberg-Marquardt法求解最小二乘代价方程。
与现有技术相比,本发明具有以下优点:
本发明可精确计算被测工件模型位置参数,省去了传统方法中繁琐困难的标定工作,对于简化偏折测量中标定工作有重要意义。
附图说明
图1为反向追迹过程图。
图2为球面镜测量结果图。
图3为球面镜测量误差图。
图4为本发明的方法流程图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
如图4所示,本发明提供一种能有效对待测工件精确定位的方法,以实现准确地估计初始测量点的坐标,具体步骤如下:
1.计算名义面形上的测量点的坐标:
相机的像素C坐标为[uc vc]T,在世界坐标系中可以表示为[xc yc zc]T。相机光心定义为Oc,则反射光线r可以定义为:
r=Rw_m·([xc yc zc]T-Tc_w) (1)
本发明所有数学表达式中,R和T分别表示旋转矩阵与平移向量,其下标表示源坐标系和目的坐标系之间的坐标转换,如Rw_c意味着从世界坐标系到相机坐标系的旋转矩阵。本专利中使用的坐标系包括相机坐标系、屏幕坐标系及模型坐标系,分别用下标c,s,m表示。
因为测量点P[xp yp f(xp,yp)]T在名义面形f上,用反射光线r所在直线方程的参数τ计算出测量点P的坐标:
[xp yp f(xp,yp)]T=τr+(Rw_m·[xc yc zc]T+Tw_m) (2)
2.逆向追迹屏幕坐标S:
被测点P的法向量n由标准方程f计算:
根据反射定律计算入射光线:
i=r-2(nT·r)n. (4)
由入射光线i和屏幕所在平面相交的参数ε求得屏幕坐标S,对应的屏幕原点为Os,也就是平移向量Ts_w。
其中,ns=Rw_mRs_w·[0 0 1]T是屏幕所在平面的法向量。
3.优化模型位置参数
将回转台旋转一定角度,位置分别记为A和B。测量两个位置下的偏折图像,假设相机中的测量点数为l,最小二乘代价方程定义为:
4.采用Levenberg-Marquardt法求解式(6),即可得到实测子孔径表面的位置。
本例中使用测量口径为70mm、半径为1000mm的球面镜作为方法验证;
将球面镜水平放置在测量平台上,并进行偏折测量;被测工件模型位置参数的初始值为:
根据系统标定参数,可以得到名义面形上的测量点的坐标,并进行反向追迹操作,得到反向追迹的屏幕坐标,最后通过优化得到实际的被测工件模型位置参数。反向追迹误差为0.07mm。计算所有估计的被测点坐标,并进行后续的面型重构,重构结果如图2所示。测量误差RMS为85.6nm,误差图如图3所示。
Claims (5)
1.一种针对光学自由曲面子孔径偏折测量的精确自定位方法,其特征在于,包括以下步骤:
1)计算名义面形f上的测量点P的坐标并获取测量点P对应的法向量n,测量点P的坐标[xp yp f(xp,yp)]T计算式为:
[xp yp f(xp,yp)]T=τr+(Rw_m·[xc yc zc]T+Tw_m)
r=Rw_m·([xc yc zc]T-Tc_w)
其中,τ为P点到相机光心的距离,r为反射光线的单位方向向量,Rw_m为从世界坐标系到模型坐标系的旋转矩阵,[xc yc zc]T为世界坐标系中相机像素C的坐标,Tw_m为从世界坐标系到模型坐标系的平移向量,Tc_w为从相机坐标系到世界坐标系的平移向量;
2)根据法向量n逆向追迹屏幕坐标S,具体为:
根据法向量n和反射定律计算入射光线i,根据入射光线i和屏幕所在平面相交的参数ε获取屏幕S的坐标,以及对应的屏幕原点和从屏幕坐标系到世界坐标系的平移向量Ts_w;
3)将回转台从位置A旋转一定角度至位置B,并测量两个位置下的偏折图像,构建最小二乘代价方程;
4)求解最小二乘代价方程,即得到实测子孔径表面的位置。
2.根据权利要求1所述的一种针对光学自由曲面子孔径偏折测量的精确自定位方法,其特征在于,所述的入射光线i的计算式为:
i=r-2(nT·r)n。
5.根据权利要求1所述的一种针对光学自由曲面子孔径偏折测量的精确自定位方法,其特征在于,所述的步骤4)中,采用Levenberg-Marquardt法求解最小二乘代价方程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910606141.0A CN110428471B (zh) | 2019-07-05 | 2019-07-05 | 一种针对光学自由曲面子孔径偏折测量的精确自定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910606141.0A CN110428471B (zh) | 2019-07-05 | 2019-07-05 | 一种针对光学自由曲面子孔径偏折测量的精确自定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110428471A CN110428471A (zh) | 2019-11-08 |
CN110428471B true CN110428471B (zh) | 2021-12-07 |
Family
ID=68410306
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910606141.0A Active CN110428471B (zh) | 2019-07-05 | 2019-07-05 | 一种针对光学自由曲面子孔径偏折测量的精确自定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110428471B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112449178B (zh) * | 2020-11-19 | 2022-07-01 | 湖北航天技术研究院总体设计所 | 一种屏幕观测设备机位标定方法及系统 |
CN115330882A (zh) * | 2021-09-16 | 2022-11-11 | 苏州大学 | 一种基于相位偏折术的标定系统及方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105783780A (zh) * | 2016-04-29 | 2016-07-20 | 浙江大学 | 一种自由曲面非常规子孔径拼接干涉检测装置与方法 |
CN108489421A (zh) * | 2018-03-12 | 2018-09-04 | 四川大学 | 一种条纹投影检测平面元件面形方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102011004376B3 (de) * | 2011-02-18 | 2012-06-21 | Carl Zeiss Smt Gmbh | Verfahren und Vorrichtung zum Bestimmen einer Form einer optischen Testfläche |
-
2019
- 2019-07-05 CN CN201910606141.0A patent/CN110428471B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105783780A (zh) * | 2016-04-29 | 2016-07-20 | 浙江大学 | 一种自由曲面非常规子孔径拼接干涉检测装置与方法 |
CN108489421A (zh) * | 2018-03-12 | 2018-09-04 | 四川大学 | 一种条纹投影检测平面元件面形方法及装置 |
Non-Patent Citations (1)
Title |
---|
Improved Subaperture Clutter Suppression Method;Zhang Yuan 等;《2015 IEEE 5th Asia-Pacific Conference on Synthetic Aperture Radar(APSAR)》;20151029;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110428471A (zh) | 2019-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110524309B (zh) | 基于四基站激光追踪系统的数控转台几何误差测量方法 | |
CN110440692B (zh) | 激光跟踪仪与结构光3d扫描仪组合式测量标定方法 | |
CN108340211B (zh) | 基于单目视觉的数控机床轮廓误差三维测量方法 | |
CN106441117B (zh) | 基于多站etalon激光跟踪系统的转台误差检测方法 | |
CN105783775A (zh) | 一种镜面及类镜面物体表面形貌测量装置与方法 | |
CN110966935B (zh) | 基于标志点的偏折测量系统一体化几何标定方法 | |
CN106767540B (zh) | 一种交会测量相机光轴与反射镜夹角误差标定方法 | |
WO2018201677A1 (zh) | 基于光束平差的远心镜头三维成像系统的标定方法及装置 | |
CN111811433B (zh) | 基于红蓝正交条纹的结构光系统标定方法及装置和应用 | |
CN108195314B (zh) | 基于多视场拼接的反射式条纹三维面形测量方法 | |
CN107957251B (zh) | 基于计算机辅助校正的反射球面通用化检测方法 | |
CN111207670A (zh) | 一种线结构光标定装置及方法 | |
CN110428471B (zh) | 一种针对光学自由曲面子孔径偏折测量的精确自定位方法 | |
CN111366079B (zh) | 一种快速确定偏折测量系统中各部件几何位置的标定方法 | |
CN106546193B (zh) | 一种高反射物体表面三维测量方法和系统 | |
CN103697811B (zh) | 一种相机与结构光源结合获取物体轮廓三维坐标的方法 | |
CN110146032B (zh) | 基于光场分布的合成孔径相机标定方法 | |
CN114943774B (zh) | 一种红外双目相机温度补偿计算方法 | |
CN109773589B (zh) | 对工件表面进行在线测量和加工导引的方法及装置、设备 | |
CN110211175B (zh) | 准直激光器光束空间位姿标定方法 | |
CN109506629B (zh) | 一种水下核燃料组件检测装置旋转中心标定的方法 | |
Liu et al. | Portable light pen 3D vision coordinate measuring system-probe tip center calibration | |
CN113916128A (zh) | 一种基于光笔式视觉测量系统的提高精度的方法 | |
CN110260817B (zh) | 基于虚拟标志点的复杂曲面偏折测量自定位方法 | |
CN105758339A (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 |