CN110567395B - 一种激光光斑尺寸边缘轮廓线绘制方法 - Google Patents

一种激光光斑尺寸边缘轮廓线绘制方法 Download PDF

Info

Publication number
CN110567395B
CN110567395B CN201910875614.7A CN201910875614A CN110567395B CN 110567395 B CN110567395 B CN 110567395B CN 201910875614 A CN201910875614 A CN 201910875614A CN 110567395 B CN110567395 B CN 110567395B
Authority
CN
China
Prior art keywords
laser spot
spot size
theta
edge contour
pixel
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
CN201910875614.7A
Other languages
English (en)
Other versions
CN110567395A (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.)
Changchun University of Science and Technology
Original Assignee
Changchun University of Science and Technology
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 Changchun University of Science and Technology filed Critical Changchun University of Science and Technology
Priority to CN201910875614.7A priority Critical patent/CN110567395B/zh
Publication of CN110567395A publication Critical patent/CN110567395A/zh
Application granted granted Critical
Publication of CN110567395B publication Critical patent/CN110567395B/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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Image Analysis (AREA)
  • Electron Beam Exposure (AREA)

Abstract

本发明公开了一种激光光斑尺寸边缘轮廓线绘制方法,先求取激光光斑质心位置并以质心为坐标原点O'建立直角坐标系,通过求解θi角度对应各有效像元(x'a,y'b)到过坐标原点O'且角度为90°+θi直线的位移diab及激光光斑尺寸ωθi,从而获得θi角度方向对应的激光光斑尺寸边缘轮廓点的坐标(xθi,yθi),通过连接不同θi角度方向的激光光斑尺寸边缘轮廓点,从而绘制出一个激光光斑尺寸边缘轮廓线。本发明通过连接各个边缘轮廓点绘制出光斑尺寸边缘轮廓曲线,直观地反映出激光光斑尺寸的几何形状。

Description

一种激光光斑尺寸边缘轮廓线绘制方法
技术领域
本发明涉及一种激光光斑尺寸边缘轮廓线绘制方法,属于激光技术领域。
背景技术
在激光器的设计、制造和应用中,对垂轴截面内激光光斑参数是评价激光光斑质量的重要指标,激光光斑尺寸边缘轮廓线能够直观地反映出激光光斑的几何形状。目前,绘制激光光斑尺寸边缘轮廓线方法主要有两种形式:一种是以光斑光强最大值的1/e2作为阈值进行提取,以大于此阈值的光斑尺寸边缘线作为光斑尺寸边缘轮廓;另一种是采用光强二阶矩法,求取X方向或同时求取X方向和Y方向的光斑尺寸,进而进行光斑尺寸边缘轮廓的圆拟合或椭圆拟合。由于1/e2计算光斑的方法无法满足ABCD定律,于是,国际标准化组织(ISO)在ISO11146文件中明确将二阶矩法作为推荐的光斑尺寸求取方法。此外,激光器输出光束的垂轴截面内光斑并非总是圆形或关于X和Y方向对称,因此,采用圆拟合或椭圆拟合无法客观真实地反映激光光斑尺寸的边缘轮廓。
发明内容
针对现有光斑尺寸边缘轮廓线绘制方法的不足,本发明的目的在于提出一种激光光斑尺寸边缘轮廓线绘制方法,通过光强二阶矩方法求解激光光斑图像任意角度方向的光斑尺寸边缘轮廓点位置,连接各边缘轮廓点获得光斑尺寸边缘轮廓线,科学全面反映激光光斑的空间特性。
如附图1所示,本发明提供的一种激光光斑尺寸边缘轮廓线绘制方法,包括以下步骤:
步骤一、获取一幅像元数宽度W、像元数长度H的激光光斑图像;所述的图像无背景噪声,且图像中无饱和的像元;以所述的图像左上角为原点O,水平向右方向为x轴正方向,竖直向下方向为y轴正方向建立直角坐标系xOy;提取所述图像中每个像元的坐标值(xa,yb)及对应的光强值I(xa,yb);
步骤二、利用公式(1)和公式(2)计算出激光光斑质心O'的坐标位置(xc,yc);
Figure GDA0002715203020000011
Figure GDA0002715203020000012
步骤三、以所述的质心O'为坐标原点,选取水平向右方向为x'轴正方向,竖直向上方向为y'轴正方向,建立直角坐标系x'O'y';
步骤四、利用公式(3)计算直角坐标系x'O'y'下每个像元的坐标(x'a,y'b),提取每个像元对应的光强值I(x'a,y'b);
Figure GDA0002715203020000021
步骤五、在直角坐标系x'O'y'下,设一过坐标原点O'且角度为90°+θi的直线方程式为:
kix'+y'=0 (4)
步骤六、在所述直角坐标系x'O'y'中,利用公式(5)求得θi角度对应的各有效像元(x'a,y'b)到方程式(4)所述直线的位移diab,利用公式(6)所示的激光光斑尺寸与光强二阶矩间关系求取角度θi方向上激光光斑尺寸ωθi,即为θi方向上激光光斑尺寸边缘轮廓点Pθi到质心O'的距离;
所述步骤六中θi角度对应的有效像元的选取方法为:
当θi大于等于0°且小于等于90°或大于270°且小于等于360°时,如附图2所示,取diab>0的像元作为有效像元;
当θi大于90°且小于等于270°时,如附图3所示,取diab<0的像元作为有效像元;
Figure GDA0002715203020000022
Figure GDA0002715203020000023
步骤七、利用公式(7)和公式(8)求得θi角度方向对应的激光光斑尺寸边缘轮廓点的坐标(xθi,yθi);
xθi=ωθi*cosθi (7)
yθi=ωθi*sinθi (8)
步骤八、重复步骤五至步骤七,求得不同θi值对应的激光光斑尺寸边缘轮廓点的坐标(xθi,yθi);
步骤九、依次连接i=1到i=n(n优选大于3的整数)所对应的各激光光斑尺寸边缘轮廓点(xθi,yθi),并连线(xθ1,yθ1)和(xθn,yθn),从而绘制出一个激光光斑尺寸边缘轮廓线。
有益效果:本发明在求取激光光斑图像任意角度方向的光斑尺寸边缘轮廓点位置基础上,通过连接各个边缘轮廓点绘制出光斑尺寸边缘轮廓曲线,直观地反映出激光光斑的几何形状。
附图说明
图1是一种激光光斑尺寸边缘轮廓线绘制方法的流程示意图。
图2是角度θi在小于90°和大于270°范围时的有效点选取区域示意图。
图3是角度θi在90°到270°范围时的有效点选取区域示意图。
图4是i取24时激光光斑尺寸边缘轮廓线示意图。
具体实施方式
实施例1一种激光光斑尺寸边缘轮廓线绘制方法。
如附图1所示,本发明提供的一种激光光斑尺寸边缘轮廓线绘制方法,包括以下步骤:
步骤一、获取一幅像元数宽度W、像元数长度H的激光光斑图像;所述的图像无背景噪声,且图像中无饱和的像元;以所述的图像左上角为原点O,水平向右方向为x轴正方向,竖直向下方向为y轴正方向建立直角坐标系xOy;提取所述图像中每个像元的坐标值(xa,yb)及对应的光强值I(xa,yb);
步骤二、利用公式(1)和公式(2)计算出激光光斑质心O'的坐标位置(xc,yc);
Figure GDA0002715203020000031
Figure GDA0002715203020000032
步骤三、以所述的质心O'为坐标原点,选取水平向右方向为x'轴正方向,竖直向上方向为y'轴正方向,建立直角坐标系x'O'y';
步骤四、利用公式(3)计算直角坐标系x'O'y'下每个像元的坐标(x'a,y'b),提取每个像元对应的光强值I(x'a,y'b);
Figure GDA0002715203020000033
步骤五、在直角坐标系x'O'y'下,设一过坐标原点O'且角度为90°+θi的直线方程式为:
kix'+y'=0 (4)
步骤六、在所述直角坐标系x'O'y'中,利用公式(5)求得θi角度对应的各有效像元(x'a,y'b)到方程式(4)所述直线的位移diab,利用公式(6)所示的激光光斑尺寸与光强二阶矩间关系求取角度θi方向上激光光斑尺寸ωθi,即为θi方向上激光光斑尺寸边缘轮廓点Pθi到质心O'的距离;
所述步骤六中θi角度对应的有效像元的选取方法为:
当θi大于等于0°且小于等于90°或大于270°且小于等于360°时,如附图2所示,取diab>0的像元作为有效像元;
当θi大于90°且小于等于270°时,如附图3所示,取diab<0的像元作为有效像元;
Figure GDA0002715203020000034
Figure GDA0002715203020000035
步骤七、利用公式(7)和公式(8)求得θi角度方向对应的激光光斑尺寸边缘轮廓点的坐标(xθi,yθi);
xθi=ωθi*cosθi (7)
yθi=ωθi*sinθi (8)
步骤八、重复步骤五至步骤七,求得不同θi值对应的激光光斑尺寸边缘轮廓点的坐标(xθi,yθi);
步骤九、依次连接i=1到i=n(n为24)所对应的各激光光斑尺寸边缘轮廓点(xθi,yθi),并连线(xθ1,yθ1)和(xθn,yθn),从而绘制出一个激光光斑尺寸边缘轮廓线,如附图4所示。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (1)

1.一种激光光斑尺寸边缘线绘制方法,其特征在于,包括以下步骤:
步骤一、获取一幅像元数宽度W、像元数长度H的激光光斑图像;所述的图像无背景噪声,且图像中无饱和的像元;以所述的图像左上角为原点O,水平向右方向为x轴正方向,竖直向下方向为y轴正方向建立直角坐标系xOy;提取所述图像中每个像元的坐标值(xa,yb)及对应的光强值I(xa,yb);
步骤二、利用公式(1)和公式(2)计算出激光光斑质心O'的坐标位置(xc,yc);
Figure FDA0002715203010000011
Figure FDA0002715203010000012
步骤三、以所述的质心O'为坐标原点,选取水平向右方向为x'轴正方向,竖直向上方向为y'轴正方向,建立直角坐标系x'O'y';
步骤四、利用公式(3)计算直角坐标系x'O'y'下每个像元的坐标(x'a,y'b),提取每个像元对应的光强值I(x'a,y'b);
Figure FDA0002715203010000013
步骤五、在直角坐标系x'O'y'下,设一过坐标原点O'且角度为90°+θi的直线方程式为:
kix'+y'=0 (4)
步骤六、在所述直角坐标系x'O'y'中,利用公式(5)求得θi角度对应的各有效像元(x'a,y'b)到方程式(4)所述直线的位移diab,利用公式(6)所示的激光光斑尺寸与光强二阶矩间关系求取角度θi方向上激光光斑尺寸ωθi,即为θi方向上激光光斑尺寸边缘轮廓点Pθi到质心O'的距离;
所述步骤六中θi角度对应的有效像元的选取方法为:
当θi大于等于0°且小于等于90°或大于270°且小于等于360°时,取diab>0的像元作为有效像元;
当θi大于90°且小于等于270°时,取diab<0的像元作为有效像元;
Figure FDA0002715203010000014
Figure FDA0002715203010000015
步骤七、利用公式(7)和公式(8)求得θi角度方向对应的激光光斑尺寸边缘轮廓点的坐标(xθi,yθi);
xθi=ωθi*cosθi (7)
yθi=ωθi*sinθi (8)
步骤八、重复步骤五至步骤七,求得不同θi值对应的激光光斑尺寸边缘轮廓点的坐标(xθi,yθi);
步骤九、依次连接i=1到i=n(n为大于3的整数)所对应的各激光光斑尺寸边缘轮廓点(xθi,yθi),并连线(xθ1,yθ1)和(xθn,yθn),从而绘制出一个激光光斑尺寸边缘轮廓线。
CN201910875614.7A 2019-09-18 2019-09-18 一种激光光斑尺寸边缘轮廓线绘制方法 Active CN110567395B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910875614.7A CN110567395B (zh) 2019-09-18 2019-09-18 一种激光光斑尺寸边缘轮廓线绘制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910875614.7A CN110567395B (zh) 2019-09-18 2019-09-18 一种激光光斑尺寸边缘轮廓线绘制方法

Publications (2)

Publication Number Publication Date
CN110567395A CN110567395A (zh) 2019-12-13
CN110567395B true CN110567395B (zh) 2021-04-06

Family

ID=68780606

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910875614.7A Active CN110567395B (zh) 2019-09-18 2019-09-18 一种激光光斑尺寸边缘轮廓线绘制方法

Country Status (1)

Country Link
CN (1) CN110567395B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114034471B (zh) * 2021-11-27 2024-04-12 北京理工大学 一种激光光路轮廓的测量方法
CN117123938B (zh) * 2023-10-12 2024-01-12 武汉金运激光股份有限公司 一种用于激光切割薄膜的功率补偿方法及其补偿系统

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2838513B1 (fr) * 2002-04-12 2004-09-10 Essilor Int Procede pour relever la forme d'un contour d'une lentille ophtalmique prealablement usinee
CN101871772B (zh) * 2010-07-14 2011-11-16 哈尔滨工业大学 基于轨迹重构的光斑轮廓测量方法
CN104316049B (zh) * 2014-10-28 2017-06-23 中国科学院长春光学精密机械与物理研究所 高精度低信噪比椭圆化星点光斑细分定位方法
CN106595482B (zh) * 2016-12-14 2019-06-04 海信集团有限公司 激光投影系统的光斑测量方法及装置
CN108844463A (zh) * 2018-06-12 2018-11-20 南京理工大学 基于激光光斑位置的地下管道空间方位测量系统及方法
CN109741266B (zh) * 2018-12-03 2021-04-02 西北核技术研究所 一种阵列探测法激光光斑图像的复原显示方法
CN109785245B (zh) * 2018-12-06 2023-03-24 江苏大学 一种光斑图像修整方法
CN110260787B (zh) * 2019-06-26 2020-12-01 王菲 一种激光光斑尺寸全角度评价与表征方法

Also Published As

Publication number Publication date
CN110567395A (zh) 2019-12-13

Similar Documents

Publication Publication Date Title
CN110567395B (zh) 一种激光光斑尺寸边缘轮廓线绘制方法
CN103639599B (zh) 一种激光去毛刺系统及方法
CN109596059B (zh) 一种基于平行线结构光的飞机蒙皮间隙与阶差测量方法
CN102411778A (zh) 一种机载激光点云与航空影像的自动配准方法
CN114119553B (zh) 一种以十字激光为基准的双目视觉异面圆孔检测方法
CN104570340B (zh) 自由曲面成像系统的设计方法
CN110487216A (zh) 一种基于卷积神经网络的条纹投影三维扫描方法
KR101053506B1 (ko) 어안 렌즈를 구비하는 카메라 모듈에 의해 획득된 영상의 왜곡 보정을 위한 파라미터 조정 방법
CN110097593A (zh) 一种从多线激光雷达点云数据中识别圆柱面的方法
CN102208025A (zh) 一种文本图像几何畸变的矫正方法
CN104748683A (zh) 一种数控机床工件在线自动测量装置及测量方法
US20110010098A1 (en) Method of pre-stack two-dimension-like transformation of three-dimensional seismic record
CN105701776B (zh) 一种用于自动光学检测的镜头畸变矫正方法及系统
CN109583365A (zh) 基于成像模型约束非均匀b样条曲线拟合车道线检测方法
CN101038163A (zh) 可变焦摄像机空间三维位姿的单视测量方法
KR102132242B1 (ko) 스테레오 카메라 보정 시스템
CN110763306A (zh) 一种基于单目视觉的液位测量系统及方法
CN107689033A (zh) 一种基于椭圆分割的鱼眼图像畸变校正方法
CN111105497A (zh) 一种三维建图智能地点识别方法
CN104537627A (zh) 一种深度图像的后处理方法
CN110889874A (zh) 一种双目相机标定结果的误差评估方法
CN109372545B (zh) 圆形盾构通缝隧道管片底部接缝张开程度评价方法
CN103884278A (zh) 一种激光跟踪仪几何角度误差综合修正方法
CN106097423A (zh) 基于k近邻的LiDAR点云强度校正方法
US8666531B2 (en) Slicing method for rapid prototyping apparatus

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