CN113806678B - 一种基于相关系数匹配的曲率半径测量方法 - Google Patents
一种基于相关系数匹配的曲率半径测量方法 Download PDFInfo
- Publication number
- CN113806678B CN113806678B CN202111042933.3A CN202111042933A CN113806678B CN 113806678 B CN113806678 B CN 113806678B CN 202111042933 A CN202111042933 A CN 202111042933A CN 113806678 B CN113806678 B CN 113806678B
- Authority
- CN
- China
- Prior art keywords
- curvature radius
- ring
- radius
- correlation coefficient
- interferogram
- 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
-
- 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
- G01B11/255—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures for measuring radius of curvature
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于相关系数匹配的曲率半径测量方法,包括(1)干涉图预处理;(2)初始曲率半径求解;(3)相关系数求解,以R∈{R初±ΔR}为索引范围,仿真出索引范围内每一个曲率半径R对应的牛顿环干涉图IR(x,y),分别计算IR(x,y)与I1(x,y)的相关系数,找出相关系数最大值对应的曲率半径R,即为待测球面镜的曲率半径。本发明采用单幅牛顿环干涉图,根据暗环半径求得一个初始曲率半径,由初始曲率半径设置索引范围,再进行相关系数匹配出高精度的待测曲率半径,其测量过程简单高效,避免了因繁琐的调试步骤而引入测量误差,使得测量方法简便高效的同时又确保了测量的精度,其测量精度优于0.5%。
Description
技术领域
本发明涉及一种曲率半径测量方法,尤其涉及一种基于相关系数匹配的曲率半径测量方法。
背景技术
球面镜在光刻物镜工程研制领域、高功率激光驱动系统、大型望远系统、高精密光学系统等方面得到广泛应用,球面镜的加工质量对整个系统的成像效果有着关键影响。曲率半径是表征球面镜的重要参数,曲率半径在加工过程中,准确测得曲率半径是光学车间检验的重要工作。因而,简便、快捷、稳定、高精度的曲率半径测量方法是研究人员们不断探索的研究方向。
中国专利公开号CN107966114A公开了一种非接触测量曲率半径的装置。包括:计算机、第一数据传输线、第二数据传输线、光栅尺定尺、光栅尺动尺、激光位移传感器、处理器、光栅数据采集卡、被测器件。通过调节激光位移传感器测出被测件距激光传感器的距离,移动光栅尺位置读取被测件表面三点的位置坐标并代入圆的方程从而求出曲率半径。该方法虽然操作简单,但是光栅尺位置的确定易产生误差,影响曲率半径测量结果。
中国专利公开号CN110455221A公开了一种快速测量光学镜片曲率半径的光路结构和设备,包括:发光LED、切换式透镜单元、棱镜、半反半透镜﹑保护玻璃、光阑和CCD摄像机。通过采用CCD图像采集和亚像素拟合圆的算法结合,并结合工控机得到被测样本的曲率半径。首先,该结构中存在多次反射改变光迹,使得适用范围变广,但是存在传输损耗,对成像质量存在影响;其次,专利中测量曲率半径的方法是通过工控机对图像进一步处理后再进行亚像素拟合圆得到的曲率半径精度较高,精度受图像质量影响大,并且稳定性降低。
中国专利公开号CN102128600B公开了一种利用激光测量透镜曲率半径的方法及其装置。包括:安装座、上盖、圆口、凹槽、凸透镜、反射镜、读数显微镜。通过度数显微镜测量出干涉条纹安环级数并用逐差法算出暗环直径,最后再由公式计算出曲率半径。该方法虽然解决了测量时透镜产生应力变形的问题,但该装置操作要求严谨,无辅助器件,调试阶段较为繁琐,不便于高效地测量透镜的曲率半径。
发明内容
发明目的:本发明目的在于提供一种基于相关系数匹配的曲率半径测量方法,在保证测量精度的情况下简化测量方法、提高测量效率和扩大测量范围。
技术方案:本发明包括以下步骤:
(1)干涉图预处理;
(2)初始曲率半径求解;
(3)相关系数求解,以R∈{R初±ΔR}为索引范围,仿真出索引范围内每一个曲率半径R对应的牛顿环干涉图IR(x,y),分别计算IR(x,y)与I1(x,y)的相关系数,找出相关系数最大值对应的曲率半径R,即为待测球面镜的曲率半径。
所述步骤(1)中的干涉图预处理为:提取采集的牛顿环干涉图有效区域I0(x,y),去除I0(x,y)的常数项得到I1(x,y),由I1(x,y)经Hough变换来标定干涉图的中心坐标。
所述的干涉图预处理具体包括:
步骤1:提取采集的牛顿环干涉图中的有效部分I0(x,y),去除牛顿环干涉图I0(x,y)的常数项,得到I1(x,y),其表达式为:
步骤2:通过边缘识别后提取圆环,经Hough变换后确定干涉图中心坐标(X1,Y1)。
所述步骤2具体包括:
1)提取圆环(x-a)2+(y-b)2=r2,其中(a,b)为提取圆环的圆心坐标,r为提取圆环的半径,(x,y)为xoy坐标系下提取圆环上任意一点的坐标;
2)提取圆环上所有的点(xi,yi)(i=1,2,3,…,n),其任意一点经Hough变换到参数空间坐标系后变成圆锥面,参数空间下的圆锥面表达式为:
(a-xi)2+(b-yi)2=r2
其中,(xi,yi)为xoy坐标系中所提取的圆环上第i个点的坐标,a、b、r为参数空间坐标系中满足圆锥表达式的变量;
对于提取的圆环上的n个点,在参数空间会形成对应的n个圆锥,且n个圆锥都会相交于一点,通过寻找参数空间上的累积最大值,得到坐标(a,b,r),其中(a,b)为干涉图的中心坐标(X1,Y1),r为提取圆环的半径。
所述步骤(2)中初始曲率半径求解为:取通过干涉图I1(x,y)中心坐标的轮廓曲线,并确定暗环半径R暗的大小,根据暗环半径求得初始曲率半径R初。
所述的初始曲率半径求解具体为:
1)取通过I1(x,y)中心坐标(X1,Y1)的轮廓曲线,提取远离中心的6~10组暗环半径R暗,对R暗进行像素标定后,根据暗环半径得到3~5个曲率半径值R0,其计算公式为:
其中,Rm和Rn分别为I1(x,y)的第m级和第n级暗环半径,λ为工作波长;
2)取R0的平均值作为初始曲率半径R初。
所述步骤(3)中相关系数的求解具体包括:
步骤1:根据初始曲率半径R初确定索引范围R∈{R初±ΔR},其中,ΔR=c·R初,0≤c≤0.2;
步骤2:仿真出每一个索引曲率半径R对应的牛顿环干涉图IR(x,y),计算IR(x,y)和I1(x,y)的相关系数,其相关公式为:
其中R为索引曲率半径,(X1,Y1)为干涉图I1(x,y)的中心坐标,(X2,Y2)为干涉图IR(x,y)的中心坐标,Cov(I1,IR)为I1和IR的协方差,δI1和δIR分别为I1(x,y)和IR(x,y)的方差;
步骤3:找出相关系数最大值ρmax对应的索引曲率半径R,即待测球面镜的曲率半径。
有益效果:本发明采用单幅牛顿环干涉图,根据暗环半径求得一个初始曲率半径,由初始曲率半径设置索引范围,再进行相关系数匹配出高精度的待测曲率半径,其测量过程简单高效,避免了因繁琐的调试步骤而引入测量误差,使得测量方法简便高效的同时又确保了测量的精度,其测量精度优于0.5%。
附图说明
图1为本发明的流程图;
图2中(a)为仿真的牛顿环干涉图,
(b)为通过仿真的牛顿环干涉图中心坐标的轮廓曲线,
(c)为仿真干涉图的暗环半径数据和求得的R初,
(d)为仿真的牛顿环干涉图在索引范围内相关系数匹配结果;
图3中(a)为实验采集的牛顿环干涉图,
(b)为通过采集的牛顿环干涉图中心坐标的轮廓曲线,
(c)为采集干涉图的暗环半径数据和求得的R初,
(d)为采集的牛顿环干涉图在索引范围内相关系数匹配结果。
具体实施方式
下面结合附图对本发明作进一步说明。
如图1所示,本发明的测量方法包括以下步骤:(1)干涉图预处理:提取采集的单幅牛顿环干涉图有效区域I0(x,y),去除I0(x,y)的常数项得到I1(x,y),由I1(x,y)经Hough变换来标定干涉图的中心坐标;(2)初始曲率半径求解:取通过干涉图I1(x,y)中心坐标的轮廓曲线,并确定暗环半径R暗的大小,根据暗环半径求得初始曲率半径R初;(3)相关系数求解:以R∈{R初±ΔR}为索引范围,仿真出索引范围内每一个曲率半径对应的牛顿环干涉图IR(x,y),分别计算IR(x,y)与I1(x,y)的相关系数,找出相关系数最大值ρmax对应的曲率半径R,即为待测球面镜的高精度曲率半径。
具体过程如下:
(1)干涉图预处理
步骤1:提取采集的牛顿环干涉图中的有效部分I0(x,y),去除牛顿环干涉图I0(x,y)的常数项,得到I1(x,y),其表达式为:
步骤2:通过边缘识别后提取一个合适的圆环,经Hough变换确定干涉图中心坐标(X1,Y1),具体为:
提取圆环的表达式:(x-a)2+(y-b)2=r2;
其中(a,b)为提取圆环的圆心坐标,r为提取圆环的半径,(x,y)为xoy坐标系下提取圆环上任意一点的坐标;
而提取圆环上所有的点(xi,yi)(i=1,2,3,…,n),其任意一点经Hough变换到参数空间坐标系下就会变成圆锥面:
参数空间下的圆锥面表达式:(a-xi)2+(b-yi)2=r2;
其中(xi,yi)为xoy坐标系中所提取的圆环上第i个点的坐标,a、b、r为参数空间坐标系中满足圆锥表达式的变量;
对于提取的圆环上的n个点,在参数空间会形成对应的n个圆锥,且n个圆锥都会相交于一点,通过寻找参数空间上的累积最大值,就可以得到坐标(a,b,r),其中(a,b)为干涉图的中心坐标(X1,Y1),r为提取圆环的半径。
(2)初始曲率半径求解
步骤1:取通过牛顿环干涉图中心坐标(X1,Y1)的轮廓曲线,选择离中心较远的6~10组暗环半径R暗,对R暗进行像素标定后,根据暗环半径得到3~5个曲率半径值R0,取R0的平均值作为初始曲率半径R初,具体为:
其中Rm和Rn分别为I1(x,y)的第m级和第n级暗环半径,λ为工作波长。
(3)相关系数求解
步骤1:根据初始曲率半径R初确定一个索引范围R∈{R初±ΔR},ΔR=c·R初,其中0≤c≤0.2;
步骤2:仿真出每一个索引曲率半径R对应的牛顿环干涉图IR(x,y),再计算IR(x,y)和I1(x,y)的相关系数,其相关公式为:
其中R为索引曲率半径,(X1,Y1)为干涉图I1(x,y)的中心坐标,(X2,Y2)为干涉图IR(x,y)的中心坐标,Cov(I1,IR)为I1和IR的协方差,δI1和δIR分别为I1(x,y)和IR(x,y)的方差;
步骤3:找出相关系数最大值ρmax对应的索引曲率半径R,ρmax对应的索引曲率半径就是待测球面镜的高精度曲率半径。
实施例1
本实施例为仿真例,球面镜的有效口径D=50mm,波长λ=632.8nm,干涉图I1(x,y)的像素大小为512×512,索引范围R∈{R初±ΔR},其ΔR=0.15·R初,标准镜的曲率半径为50000mm。
曲率半径为50000mm的仿真牛顿环干涉图I1(x,y)如图2(a)所示;提取采集牛顿环干涉图,提取通过牛顿环干涉图中心坐标的轮廓曲线如图2(b)所示;根据轮廓曲线提取离中心较远的10组暗环半径R暗,并求得初始曲率半径R初,如图2(c)所示;由初始曲率半径R初,设定索引范围内的相关系数匹配结果,如图2(d)所示。结果显示相关系数最大值ρmax对应的索引曲率半径为50001.3955mm。
实施例2
本实施例为实验例,球面镜的有效口径D=50mm,波长λ=632.8nm,干涉图I1(x,y)的像素大小为512×512,索引范围R∈{R初±ΔR},其ΔR=0.15·R初,标准镜的曲率半径为41400mm。
实验采集的牛顿环干涉图I1(x,y)如图3(a)所示;提取通过牛顿环干涉图中心坐标的轮廓曲线如图3(b)所示;根据轮廓曲线提取离中心较远的10组暗环半径R暗,并求得初始曲率半径R初,如图3(c)所示;在索引范围内的相关系数匹配结果如图3(d)所示。结果显示相关系数最大值ρmax对应的索引曲率半径为41219.8410mm,相较于标准球面镜曲率半径41400mm,其精确度优于0.44%。
Claims (4)
1.一种基于相关系数匹配的曲率半径测量方法,其特征在于,包含以下步骤:
(1)干涉图预处理:提取采集的牛顿环干涉图有效区域I0(x,y),去除I0(x,y)的常数项得到I1(x,y),由I1(x,y)经Hough变换来标定干涉图的中心坐标;
(2)初始曲率半径求解:取通过干涉图I1(x,y)中心坐标的轮廓曲线,并确定暗环半径R暗的大小,根据暗环半径求得初始曲率半径R初;
(3)相关系数求解,以R∈{R初±ΔR}为索引范围,仿真出索引范围内每一个曲率半径R对应的牛顿环干涉图IR(x,y),其中,R初为初始求解的曲率半径,ΔR是曲率半径搜索的界限,分别计算IR(x,y)与I1(x,y)的相关系数,找出相关系数最大值对应的曲率半径R,即为待测球面镜的曲率半径,具体包括:
步骤3.1:根据初始曲率半径R初确定索引范围R∈{R初±ΔR},其中,ΔR=c·R初,0≤c≤0.2;
步骤3.2:仿真出每一个索引曲率半径R对应的牛顿环干涉图IR(x,y),计算IR(x,y)和I1(x,y)的相关系数,其相关公式为:
其中R为索引曲率半径,(X1,Y1)为干涉图I1(x,y)的中心坐标,(X2,Y2)为干涉图IR(x,y)的中心坐标,Cov(I1,IR)为I1和IR的协方差,和分别为I1(x,y)和IR(x,y)的方差,λ为工作波长;
步骤3.3:找出相关系数最大值ρmax对应的索引曲率半径R,即待测球面镜的曲率半径。
3.根据权利要求2所述的一种基于相关系数匹配的曲率半径测量方法,其特征在于,所述步骤1.2具体包括:
1)提取圆环(x-a)2+(y-b)2=r2,其中(a,b)为提取圆环的圆心坐标,r为提取圆环的半径,(x,y)为xoy坐标系下提取圆环上任意一点的坐标;
2)提取圆环上所有的点(xi,yi)(i=1,2,3,…,n),其任意一点经Hough变换到参数空间坐标系后变成圆锥面,参数空间下的圆锥面表达式为:
(a-xi)2+(b-yi)2=r2
其中,(xi,yi)为xoy坐标系中所提取的圆环上第i个点的坐标,a、b、r为参数空间坐标系中满足圆锥表达式的变量;
对于提取的圆环上的n个点,在参数空间会形成对应的n个圆锥,且n个圆锥都会相交于一点,通过寻找参数空间上的累积最大值,得到坐标(a,b,r),其中(a,b)为干涉图的中心坐标(X1,Y1),r为提取圆环的半径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111042933.3A CN113806678B (zh) | 2021-09-07 | 2021-09-07 | 一种基于相关系数匹配的曲率半径测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111042933.3A CN113806678B (zh) | 2021-09-07 | 2021-09-07 | 一种基于相关系数匹配的曲率半径测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113806678A CN113806678A (zh) | 2021-12-17 |
CN113806678B true CN113806678B (zh) | 2022-08-26 |
Family
ID=78940555
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111042933.3A Active CN113806678B (zh) | 2021-09-07 | 2021-09-07 | 一种基于相关系数匹配的曲率半径测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113806678B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106595529A (zh) * | 2016-03-15 | 2017-04-26 | 南京理工大学 | 基于虚拟牛顿环的大曲率半径非零位干涉测量方法及装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4183220B2 (ja) * | 2000-02-22 | 2008-11-19 | フジノン株式会社 | 光学球面曲率半径測定装置 |
US6732046B1 (en) * | 2001-10-03 | 2004-05-04 | Navigation Technologies Corp. | Application of the hough transform to modeling the horizontal component of road geometry and computing heading and curvature |
DE102008049159B4 (de) * | 2008-09-25 | 2010-09-09 | Jenoptik Laser, Optik, Systeme Gmbh | Verfahren und Vorrichtung zur optischen Messung langer Krümmungsradien optischer Prüfflächen |
CN102128600B (zh) * | 2010-12-10 | 2012-02-01 | 西安科技大学 | 一种利用激光测量透镜曲率半径的方法及其装置 |
CN108895987A (zh) * | 2018-07-17 | 2018-11-27 | 苏州大学 | 基于复合涡旋光干涉的透镜曲率半径测量方法 |
CN110455221B (zh) * | 2019-09-11 | 2024-04-09 | 大连鉴影光学科技有限公司 | 一种快速测量光学镜片曲率半径的光路结构和设备 |
CN112066913B (zh) * | 2020-08-28 | 2021-09-21 | 四川大学 | 一种测量球面曲率半径的方法 |
-
2021
- 2021-09-07 CN CN202111042933.3A patent/CN113806678B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106595529A (zh) * | 2016-03-15 | 2017-04-26 | 南京理工大学 | 基于虚拟牛顿环的大曲率半径非零位干涉测量方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN113806678A (zh) | 2021-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103499297B (zh) | 一种基于ccd的高精度测量方法 | |
Liu et al. | An improved online dimensional measurement method of large hot cylindrical forging | |
WO2006026214A1 (en) | Measuring sub-wavelength displacements | |
US20210358160A1 (en) | Method and system for determining plant leaf surface roughness | |
AU2021104874A4 (en) | A surface roughness detection method based on machine vision and machine learning | |
CN110727093A (zh) | 多光谱显微自动聚焦装置及方法 | |
TWI407078B (zh) | Micro - lens array surface profile detection system and its detection method | |
CN113806678B (zh) | 一种基于相关系数匹配的曲率半径测量方法 | |
JP2006343234A (ja) | 非球面形状誤差の評価方法および形状評価装置 | |
CN108225187B (zh) | 一种基于波前传感的非球面透镜误差检测方法 | |
CN111986315B (zh) | 一种基于深度学习的超分辨白光干涉图三维重建方法 | |
Zhang et al. | Accurate profile measurement method for industrial stereo-vision systems | |
CN116593137B (zh) | 一种基于干涉仪的光学镜片质量测试方法及系统 | |
WO2021179400A1 (zh) | 一种基于计算机视觉的装配过程几何参数自适应测量系统及方法 | |
CN112990190A (zh) | 指针式仪表读数的识别方法、装置 | |
CN109781153B (zh) | 物理参数估计方法、装置和电子设备 | |
Zhou et al. | Multitarget fiber high-precision position detection method based on a front light source | |
Zhao et al. | Sub-microscale precision repeatability position measurement using integrated polar microstructure and feature extraction method | |
CN114509018A (zh) | 一种全场实时桥梁挠度测量方法 | |
CN210719032U (zh) | 位移测试系统 | |
Li et al. | An on-line measurement method of the medium thickness steel plate based on structured light vision sensor | |
CN109489546B (zh) | 基于弹性特征图像尺的定位量测系统 | |
CN112949197B (zh) | 一种基于深度学习的平凸透镜曲率半径在线测量方法 | |
CN110686602A (zh) | 位移测试系统及位移测试方法 | |
CN114791266B (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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230920 Address after: 212399 Danyang hi tech Innovation Park, South Third Ring Road, Yunyang street, Danyang City, Zhenjiang City, Jiangsu Province Patentee after: Danyang Juxiang Vision Technology Co.,Ltd. Address before: No.2, Mengxi Road, Jingkou District, Zhenjiang City, Jiangsu Province, 212008 Patentee before: JIANGSU University OF SCIENCE AND TECHNOLOGY |
|
TR01 | Transfer of patent right |