CN110567395B - 一种激光光斑尺寸边缘轮廓线绘制方法 - Google Patents
一种激光光斑尺寸边缘轮廓线绘制方法 Download PDFInfo
- 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
- edge contour
- theta
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000006073 displacement reaction Methods 0.000 claims abstract description 4
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 238000010187 selection method Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 3
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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)
- Electron Beam Exposure (AREA)
- Image Analysis (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)；
步骤三、以所述的质心O'为坐标原点，选取水平向右方向为x'轴正方向，竖直向上方向为y'轴正方向，建立直角坐标系x'O'y'；
步骤四、利用公式(3)计算直角坐标系x'O'y'下每个像元的坐标(x'a，y'b)，提取每个像元对应的光强值I(x'a，y'b)；
步骤五、在直角坐标系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的像元作为有效像元；
步骤七、利用公式(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)；
步骤三、以所述的质心O'为坐标原点，选取水平向右方向为x'轴正方向，竖直向上方向为y'轴正方向，建立直角坐标系x'O'y'；
步骤四、利用公式(3)计算直角坐标系x'O'y'下每个像元的坐标(x'a，y'b)，提取每个像元对应的光强值I(x'a，y'b)；
步骤五、在直角坐标系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的像元作为有效像元；
步骤七、利用公式(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)；
步骤三、以所述的质心O'为坐标原点，选取水平向右方向为x'轴正方向，竖直向上方向为y'轴正方向，建立直角坐标系x'O'y'；
步骤四、利用公式(3)计算直角坐标系x'O'y'下每个像元的坐标(x'a，y'b)，提取每个像元对应的光强值I(x'a，y'b)；
步骤五、在直角坐标系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的像元作为有效像元；
步骤七、利用公式(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)，从而绘制出一个激光光斑尺寸边缘轮廓线。
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)
| 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)
| 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 | 王菲 | 一种激光光斑尺寸全角度评价与表征方法 |
-
2019
- 2019-09-18 CN CN201910875614.7A patent/CN110567395B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110567395A (zh) | 2019-12-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| TWI539184B (zh) | 自由曲面成像系統的設計方法 | |
| CN110567395B (zh) | 一种激光光斑尺寸边缘轮廓线绘制方法 | |
| CN103996173B (zh) | 一种基于变长轴椭圆拟合的鱼眼图像校正方法 | |
| CN103639599B (zh) | 一种激光去毛刺系统及方法 | |
| CN104299249B (zh) | 高鲁棒性的标志点解码方法及系统 | |
| CN102411778A (zh) | 一种机载激光点云与航空影像的自动配准方法 | |
| CN102927921B (zh) | 基于光学相似三角形法的立木胸径测量方法 | |
| CN107590825A (zh) | 一种基于sfm的点云孔洞修补方法 | |
| CN103366006B (zh) | 一种基于三维激光点云建立地理空间索引的方法 | |
| CN104898275A (zh) | 自由曲面离轴三反光学系统的设计方法 | |
| CN105427304B (zh) | 基于多种特征联合的目标sar图像和光学图像配准方法 | |
| CN109579701A (zh) | 基于结构光视觉测量系统的椭圆中心投影畸变消除方法 | |
| CN106570905A (zh) | 一种非合作目标点云初始姿态验证方法 | |
| CN108152939B (zh) | 离轴非球面三反光学系统 | |
| CN104182933A (zh) | 一种基于逆向除法模型的广角镜头图像畸变校正方法 | |
| CN107689033A (zh) | 一种基于椭圆分割的鱼眼图像畸变校正方法 | |
| CN105488807A (zh) | 一种远心镜头的标定和矫正方法 | |
| CN107369172B (zh) | 一种智能设备和输出深度图像的方法 | |
| CN104853180A (zh) | 图像获取系统的校正方法 | |
| CN107492080A (zh) | 免定标便捷的单目镜头图像径向畸变矫正方法 | |
| CN103884278B (zh) | 一种激光跟踪仪几何角度误差综合修正方法 | |
| CN102903092A (zh) | 一种基于四点变换的图像自适应校正方法 | |
| CN117522803A (zh) | 基于双目视觉和目标检测的桥梁部件精确定位方法 | |
| WO2010066123A1 (zh) | 障碍物分割方法及装置 | |
| 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 |