CN104318047B - 一种激光标刻***几何畸变的校正方法 - Google Patents
一种激光标刻***几何畸变的校正方法 Download PDFInfo
- Publication number
- CN104318047B CN104318047B CN201410411814.4A CN201410411814A CN104318047B CN 104318047 B CN104318047 B CN 104318047B CN 201410411814 A CN201410411814 A CN 201410411814A CN 104318047 B CN104318047 B CN 104318047B
- Authority
- CN
- China
- Prior art keywords
- deflection angle
- valuation
- galvanometer deflection
- galvanometers
- alpha
- 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
Landscapes
- Laser Beam Processing (AREA)
Abstract
本发明提供了一种激光标刻***几何畸变的校正方法,包含下述步骤:(a)利用目标标刻坐标计算扫描振镜偏转角度的估值;(b)以扫描振镜偏转角度的估值为初值,依据扫描振镜的扫描几何公式进行迭代计算,利用迭代终止时所得的振镜偏转角度数值解与偏转角度估值的差值进行补偿运算;计算所得的补偿值与利用实际标刻坐标计算所得的扫描振镜偏转角度估值的代数和即为实际标刻坐标所对应的扫描振镜偏转角度。该方法提出的依据扫描振镜数值解建立适当规模的误差校正表,并选择合适的插值算法能够使得经校正后的几何畸变误差满足标刻加工精度的要求,而且校正精度通过离线仿真是可验证的、可控的。
Description
技术领域
本发明涉及一种几何畸变校正方法,尤其是激光标刻***的几何畸变校正方法。
背景技术
激光标刻技术是通过将激光聚焦于加工物件表面,利用相应的物理或化学变化在材料表面形成标记的一种标刻方式。与传统的化学腐蚀、机械刻划、丝印及喷墨等标刻方式相比,具有显著的优势,能够实现文字、图形、位图、二维码等信息的标记。激光标刻技术的发展经历了机械扫描、转镜扫描和振镜扫描三个发展时代。振镜扫描因其在标刻速度、标刻质量、应用范围等方面的优势业已成为激光标刻领域的主流技术。振镜扫描技术的光学属性必然带来光学设备固有的误差,如透镜常见的径向畸变与切向畸变误差、f·θ透镜特有的场曲误差、设备结构与安装精度导致的垂直度误差与零点误差、控制***硬件导致的非线性误差以及控制***控制策略导致的复合几何畸变误差等。其中,扫描振镜的复合几何畸变误差的影响最为严重,随着标刻区域的扩大,复合几何畸变也愈加突出,严重影响了激光标刻的加工精度。
针对激光扫描振镜的复合几何畸变,主要有硬件校正与软件校正两类解决方案。元器件参数的一致性、稳定性问题及参数调整的灵活性等问题限制了硬件校正的性能与适用性,故除早期有部分文献的硬件校正措施外,现多采用软件校正。最早提出的软件校正方法是增量补偿法,利用振镜扫描***扫描几何畸变公式得到的畸变量,采用计算机软件修正法,对振镜标刻坐标数据进行修正,明显改善了几何畸变但仍然存在一定的残差。利用曲面最小二乘拟合的方法及通过基于平面坐标变换几何校正算法原理的图形畸变校正方法也是常见的方法。然而,在对测量数据进行曲面拟合或坐标变换时,通常采用均方根误差为常用的误差指标,而没有考虑样本数据的概率统计特性,限制了校正的精度。
发明内容
本发明的目的是提供种一种激光标刻***几何畸变的校正方法,能够使得经校正后的几何失真误差满足标刻加工精度的要求,而且校正精度通过离线仿真是可验证的、可控的。
本发明的目的是通过下述技术方案实现的:一种激光标刻***几何畸变的校正方法,至少包含下述步骤:
(a)离线时,利用目标标刻坐标计算扫描振镜偏转角度的估值,并以该估值为初值,依据扫描振镜的扫描几何公式进行迭代计算,利用迭代终止时所得的振镜偏转角度数值解与偏转角度估值的差值建立校正表;
(b)在线时,利用目标标刻坐标计算扫描振镜偏转角度的估值,使用步骤(a)所得的校正表对所得的扫描振镜偏转角度估值进行补偿;
其中步骤(a)包括选取有限的、有代表性的目标标刻点坐标及其对应的振镜偏转角度数值解与偏转角度估值的差值分别建立x振镜与y振镜的校正表;
其中步骤(b)包括根据目标标刻坐标通过校正表的查找与插值计算获得扫描振镜偏转角度的补偿值,该补偿值与根据实际标刻坐标计算所得的扫描振镜偏转角度估值的代数和即为目标标刻坐标所对应的扫描振镜偏转角度;
所述步骤(a)、步骤(b)中的扫描振镜偏转角度的估值计算过程如下:
设目标标刻坐标为xd,yd,f·θ透镜的焦距为f,扫描振镜偏转角度的估值计算如下:
式中:分别为x振镜、y振镜偏转角度的估值
所述步骤(a)中的迭代计算过程如下:
迭代初值:
迭代公式:
式中:xd,yd为标刻点的坐标,f为f·θ透镜的焦距,αn与βn分别是第n次迭代计算所得的x振镜与y振镜偏转角度,αn+1与βn+1分别是第n+1次迭代计算所得的x振镜与y振镜偏转角度;
所述步骤(a)中与当前标刻点对应的振镜偏转角度数值解与偏转角度估值的差值计算过程如下:
式中:Δα与Δβ分别是x振镜与y振镜偏转角度数值解与偏转角度估值的差值,αd与βd分别是迭代终止时所得的x振镜与y振镜偏转角度数值解;
所述步骤(b)中的振镜偏转角度补偿计算
α=xd/2f+LUTα(xd,yd)
β=yd/2f+LUTβ(xd,yd)
式中,α与β分别是为目标标刻坐标所对应的x振镜与y振镜偏转角度,LUTα(·)与LUTβ(·)分别是根据目标标刻坐标xd,yd通过校正表的查找与插值计算获得扫描振镜偏转角度的补偿值。
所述的校正表可以只覆盖第I象限,第I象限的振镜偏转角度补偿值可以利用该校正表得出,其他3个象限的振镜偏转角度补偿值亦可以依据几何畸变的对称性根据覆盖第I象限的校正表计算得出;所述4个象限的振镜偏转角度补偿计算如下:
α=sgn(xd)×(abs(xd)/2f+LUT′α(abs(xd),abs(yd)))
β=sgn(yd)×(abs(yd)/2f+LUT′β(abs(xd),abs(yd)))
式中,α与β分别是为目标标刻坐标所对应的x振镜与y振镜的偏转角度,sgn(·)是符号函数,LUT′α(·)与LUT′β(·)分别是根据建立在第I象限的校正表经查找与插值计算所得的扫描振镜偏转角度的补偿值。
本发明提出的方法,能够使得经校正后的几何失真误差满足标刻加工精度的要求,而且校正精度通过离线仿真是可验证的、可控的。
附图说明
图1为几何畸变仿真图。
图2a为采用双线性插值计算校正后的位置误差曲面。
图2b为采用双二次插值计算校正后的位置误差曲面。
具体实施方式
下面将对本发明的具体实施方式进行描述,但是应当明白,本发明并不仅仅局限于所描述的具体实施方式,这些具体实施方式仅仅是说明性举例,而非对本发明作出限制,本领域技术人员能够根据具体情形对具体实施方式进行一些变化,而这些变化都在本发明的保护范围之内。
激光标刻***的核心部件是扫描振镜,其工作方式是典型的物镜前扫描。激光束经x振镜与y振镜的的反射后由f·θ透镜聚焦于加工工作面。振镜偏转角度与标刻坐标之间存在复杂的约束关系。
式中:x、y为标刻点的坐标,f为f·θ透镜的焦距,α与β分别是x振镜与y振镜的偏转角度。可以看出,振镜偏转角度与标刻坐标之间是复杂的非线性关系,而且标刻点的横坐标与纵坐标同时受到两只振镜偏转角度的影响。唯一的例外是,当振镜沿x轴扫描(β=0)或y轴扫描(α=0)时,这种复杂的非线性关系退化为线性关系。
x=2fα (3)
y=2fβ (4)
扫描振镜工作时,两只振镜的偏转是在控制信号的驱动下完成的。激光标刻***的控制单元需将目标标刻坐标换算为振镜偏转角度,并周期性地将偏转角度信息以模拟信号或数字信号的方式发送至扫描振镜,刷新间隔通常为10μS。由公式(1)与公式(2)的形式不难看出,在如此短的时间内完成偏转角度的解算对控制单元的处理器来说是一个极大的挑战。为减轻处理器的负荷,通常采用的方法是将公式(3)与公式(4)作为公式(1)和公式(2)的近似,以简化解算过程。这种近似必然带来相应的误差。
图1是仿真焦距为300mm的f·θ透镜在标刻正方形轨迹时的结果,正方形的尺寸为400mm×400mm。图中虚线所示的轨迹是目标标刻轨迹,实线所示的轨迹是采用公式(3)和公式(4)近似解算扫描振镜偏转角度时的实际标刻轨迹。实际标刻轨迹的失真是意料之中的,即y方向的桶形失真和x方向的枕形失真,而且这种复合几何畸变会随着振镜扫描角度的增加而加重。
本实施例以焦距为300mm的f·θ透镜和标刻范围为400mm×400mm为例,在第I象限以5mm的间隔选取多个目标标刻点,在离线的条件下,对每一个预选的标刻点做做如下处理:
首先利用目标标刻坐标计算扫描振镜偏转角度的估值,
以扫描振镜偏转角度的估值为初值,依据下式进行迭代计算,
迭代20次后,迭代终止。利用下式计算与当前标刻点对应的扫描振镜偏转角度补偿值,
根据所得的多组标刻坐标及其对应的扫描振镜偏转角度补偿值建立第I象限的校正表。
在线工作时,激光标刻***的控制单元根据实际目标标刻坐标利用下式换算振镜偏转角度,
α=sgn(xd)×(abs(xd)/2f+LUT′α(abs(xd),abs(yd)))
β=sgn(yd)×(abs(yd)/2f+LUT′β(abs(xd),abs(yd)))
式中,sgn(·)是符号函数,LUT′α(·)与LUT′β(·)分别是根据建立在第I象限的校正表查找与插值计算结果。图2给出了在本实施例中分别采用双线性插值(图2a)与双二次插值(图2b)算法校正后的位置残差曲面。
依据偏转角度补偿值建立误差校正表并结合适当的插值计算,是在实时条件下完成由目标标刻坐标到振镜偏转角度映射可行的方法。依据几何畸变的对称性,误差校正表只覆盖第I象限,其他3个象限的角度补偿值根据几何畸变的对称性由覆盖第I象限的校正表换算得出,这在一定程度上减小了误差校正表的规模。
本发明提出的依据扫描振镜数值解建立适当规模的误差校正表,并选择合适的插值算法能够使得经校正后的几何畸变误差满足标刻加工精度的要求,而且校正精度通过离线仿真是可验证的、可控的。
Claims (1)
1.一种激光标刻***几何畸变的校正方法,其特征在于:至少包含下述步骤:
(a)离线时,利用目标标刻坐标计算扫描振镜偏转角度的估值,并以该估值为初值,依据扫描振镜的扫描几何公式进行迭代计算,利用迭代终止时所得的振镜偏转角度数值解与偏转角度估值的差值建立校正表;
(b)在线时,利用目标标刻坐标计算扫描振镜偏转角度的估值,使用步骤(a)所得的校正表对所得的扫描振镜偏转角度估值进行补偿;
其中步骤(a)包括选取有限的、有代表性的目标标刻点坐标及其对应的振镜偏转角度数值解与偏转角度估值的差值分别建立x振镜与y振镜的校正表;
其中步骤(b)包括根据目标标刻坐标通过校正表的查找与插值计算获得扫描振镜偏转角度的补偿值,该补偿值与根据实际标刻坐标计算所得的扫描振镜偏转角度估值的代数和即为目标标刻坐标所对应的扫描振镜偏转角度;
所述步骤(a)、步骤(b)中的扫描振镜偏转角度的估值计算过程如下:
设目标标刻坐标为xd,yd,f·θ透镜的焦距为f,扫描振镜偏转角度的估值计算如下:
式中:分别为x振镜、y振镜偏转角度的估值
所述步骤(a)中的迭代计算过程如下:
迭代初值:
迭代公式:
式中:xd,yd为标刻点的坐标,f为f·θ透镜的焦距,αn与βn分别是第n次迭代计算所得的x振镜与y振镜偏转角度,αn+1与βn+1分别是第n+1次迭代计算所得的x振镜与y振镜偏转角度;
所述步骤(a)中与当前标刻点对应的振镜偏转角度数值解与偏转角度估值的差值计算过程如下:
式中:Δα与Δβ分别是x振镜与y振镜偏转角度数值解与偏转角度估值的差值,αd与βd分别是迭代终止时所得的x振镜与y振镜偏转角度数值解;
所述步骤(b)中的振镜偏转角度补偿计算
α=xd/2f+LUTα(xd,yd)
β=yd/2f+LUTβ(xd,yd)
式中,α与β分别是为目标标刻坐标所对应的x振镜与y振镜偏转角度,LUTα(·)与LUTβ(·)分别是根据目标标刻坐标xd,yd通过校正表的查找与插值计算获得的x振镜与y振镜偏转角度的补偿值;
所述的校正表只覆盖第I象限,第I象限的振镜偏转角度补偿值利用该校正表得出,其他3个象限的振镜偏转角度补偿值依据几何畸变的对称性根据覆盖第I象限的校正表计算得出;所述4个象限的振镜偏转角度补偿计算如下:
α=sgn(xd)×(abs(xd)/2f+LUT′α(abs(xd),abs(yd)))
β=sgn(yd)×(abs(yd)/2f+LUT′β(abs(xd),abs(yd)))
式中,α与β分别是为目标标刻坐标所对应的x振镜与y振镜的偏转角度,sgn(·)是符号函数,LUT′α(·)与LUT′β(·)分别是根据覆盖第I象限的校正表经查找与插值计算后所得的x振镜与y振镜偏转角度的补偿值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410411814.4A CN104318047B (zh) | 2014-08-18 | 2014-08-18 | 一种激光标刻***几何畸变的校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410411814.4A CN104318047B (zh) | 2014-08-18 | 2014-08-18 | 一种激光标刻***几何畸变的校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104318047A CN104318047A (zh) | 2015-01-28 |
CN104318047B true CN104318047B (zh) | 2017-04-05 |
Family
ID=52373278
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410411814.4A Active CN104318047B (zh) | 2014-08-18 | 2014-08-18 | 一种激光标刻***几何畸变的校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104318047B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105991949A (zh) * | 2015-03-06 | 2016-10-05 | 深圳市恒瀚电子科技有限公司 | 一种激光扫描图像畸变校正方法、装置及*** |
CN105252911B (zh) * | 2015-09-22 | 2017-04-12 | 深圳市创鑫激光股份有限公司 | 一种激光打标的校正方法和装置 |
CN106001933B (zh) * | 2016-05-10 | 2017-08-11 | 湖南大捷智能装备有限公司 | 激光切割修边线的优化方法 |
CN109865943A (zh) * | 2019-01-08 | 2019-06-11 | 合肥泰沃达智能装备有限公司 | 一种超高速区块链激光雕刻*** |
CN109759714B (zh) * | 2019-01-17 | 2021-04-02 | 南开大学 | 一种基于飞秒激光成丝的大幅面打标***及打标范围标定方法 |
CN111275761B (zh) * | 2020-01-17 | 2023-10-03 | 湖北三江航天红峰控制有限公司 | 一种自适应高度的视觉定位激光标刻方法 |
CN111796256B (zh) * | 2020-07-08 | 2024-03-29 | 宁波傲视智绘光电科技有限公司 | 激光雷达用扫描***的位置设置方法、设备和存储装置 |
CN115037369A (zh) * | 2022-05-31 | 2022-09-09 | 中国科学院上海光学精密机械研究所 | 基于星载协作的激光通信载荷瞄准方法及*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101718958A (zh) * | 2009-11-24 | 2010-06-02 | 中国科学院长春光学精密机械与物理研究所 | 绿光调阻的聚焦透镜非线性像差校正的控制***和方法 |
CN101804521A (zh) * | 2010-04-15 | 2010-08-18 | 中国电子科技集团公司第四十五研究所 | 振镜***校正装置及其校正方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002316276A (ja) * | 2001-04-18 | 2002-10-29 | Nec Corp | レーザマーキング装置のマーキング歪み補正装置及びマーキング補正方法 |
JP2002321072A (ja) * | 2001-04-27 | 2002-11-05 | Sunx Ltd | レーザマーキング装置 |
-
2014
- 2014-08-18 CN CN201410411814.4A patent/CN104318047B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101718958A (zh) * | 2009-11-24 | 2010-06-02 | 中国科学院长春光学精密机械与物理研究所 | 绿光调阻的聚焦透镜非线性像差校正的控制***和方法 |
CN101804521A (zh) * | 2010-04-15 | 2010-08-18 | 中国电子科技集团公司第四十五研究所 | 振镜***校正装置及其校正方法 |
Non-Patent Citations (3)
Title |
---|
Correction of field distortion of laser marking systems using surface compensation function;Ming-Fei Chen 等;《Optics & Lasers in Engineering》;20090131;第47卷(第1期);第84-89页 * |
双振镜激光扫描的误差分析及校正方法;韩万鹏 等;《光电技术应用》;20110815;第26卷(第4期);第14-18页 * |
激光打标中振镜几何畸变误差多项式拟合校正;彭留永 等;《激光技术》;20130925;第37卷(第5期);第601-604页第1小节,第2.1-2.2小节,图2 * |
Also Published As
Publication number | Publication date |
---|---|
CN104318047A (zh) | 2015-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104318047B (zh) | 一种激光标刻***几何畸变的校正方法 | |
CN109471101B (zh) | 一种弹载sar图像校正方法 | |
CN106052556B (zh) | 一种三坐标测量机空域坐标修正方法 | |
CN103901428A (zh) | 弹载sar子孔径前斜视高阶非线性调频变标成像方法 | |
CN107290735B (zh) | 一种基于自制地基激光雷达铅垂度误差的点云误差校正方法 | |
CN102538823B (zh) | 大尺寸tdiccd焦面异速成像匹配误差检测*** | |
CN204313798U (zh) | 一种激光束现场标定装置 | |
CN107768829B (zh) | 一种基于太阳跟踪的天线指向修正方法 | |
CN107314740B (zh) | 一种双振镜扫描***的位置校准***及校准方法 | |
CN105974415A (zh) | 一种机载sar方位空变运动误差的高精度补偿方法 | |
CN105068065A (zh) | 星载激光测高仪在轨检校方法及*** | |
CN101046386A (zh) | 太阳敏感器测量基准转换方法与装置 | |
CN103926548A (zh) | 一种快速测量射电望远镜反射面精度的方法 | |
CN110927965B (zh) | 用于补偿因光线偏折产生误差的补偿透镜的设计方法 | |
CN106092515A (zh) | 一种用法布里-珀罗标准具测量焦距和转角的方法 | |
CN104236857A (zh) | 基于四分之一波片法的液晶光学移相分布检测***及方法 | |
CN111679261B (zh) | 一种基于反光板的激光雷达定位方法及*** | |
CN109708559A (zh) | 一种基于隅角镜的光电自准直仪的测角方法 | |
Shimanskii et al. | Alignment of the writing beam with the diffractive structure rotation axis in synthesis of diffractive optical elements in a polar coordinate system | |
CN110108271B (zh) | 一种气动光学效应引起的星光偏折补偿方法 | |
CN102305929A (zh) | 机载合成孔径雷达杠杆臂误差补偿方法 | |
CN115951502A (zh) | 一种利用激光跟踪仪的平行光管焦面组件装调方法 | |
CN104949692A (zh) | 多波束测深***安装校正值的计算方法 | |
CN105758428B (zh) | 利用标定装置标定动态靶标动态形变测角误差的方法 | |
CN105136295B (zh) | 一种aotf同一幅图中光谱不均匀的解决方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |