CN109188681B - 经F-theta聚焦透镜后引起的畸变误差补偿方法 - Google Patents

经F-theta聚焦透镜后引起的畸变误差补偿方法 Download PDF

Info

Publication number
CN109188681B
CN109188681B CN201811130506.9A CN201811130506A CN109188681B CN 109188681 B CN109188681 B CN 109188681B CN 201811130506 A CN201811130506 A CN 201811130506A CN 109188681 B CN109188681 B CN 109188681B
Authority
CN
China
Prior art keywords
focusing lens
galvanometer
theta
curved surface
sin
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
CN201811130506.9A
Other languages
English (en)
Other versions
CN109188681A (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.)
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for 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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN201811130506.9A priority Critical patent/CN109188681B/zh
Publication of CN109188681A publication Critical patent/CN109188681A/zh
Application granted granted Critical
Publication of CN109188681B publication Critical patent/CN109188681B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B26/00Optical devices or arrangements for the control of light using movable or deformable optical elements
    • G02B26/08Optical devices or arrangements for the control of light using movable or deformable optical elements for controlling the direction of light
    • G02B26/10Scanning systems
    • G02B26/101Scanning systems with both horizontal and vertical deflecting means, e.g. raster or XY scanners
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/0012Optical design, e.g. procedures, algorithms, optimisation routines

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Laser Beam Processing (AREA)
  • Mechanical Optical Scanning Systems (AREA)

Abstract

本发明涉及一种经F‑theta聚焦透镜后引起的畸变误差补偿方法,针对由添加F‑theta聚焦透镜引起的桶形畸变,搭建测试光路,研究出光束从路径BDFG映射到坐标平面后GH的距离,通过此距离和BD的方向矢量(通过振镜偏转角以及几何关系计算)计算出G点的映射坐标值,进而可得到计算机数字控制量与振镜偏转角、映射坐标值的数学关系。利用模型找到计算机数字控制量与映射坐标值的关系,实际运用时通过修改输入计算机数字控制量即可达到补偿的效果。

Description

经F-theta聚焦透镜后引起的畸变误差补偿方法
技术领域
本发明涉及一种误差处理技术,特别涉及一种经F-theta聚焦透镜后引起的畸变误差补偿方法。
背景技术
近年来双振镜激光扫描加工技术越来越普及,其加工效率高,加工精度好,因此该技术越来越流行,为了达到更高的加工精度就必须研究此项技术的畸变误差。通过分析总结,我们发现显而易见的误差有映射关系非线性引起的误差,添加F-theta聚焦透镜引起的误差,以及扫描反射镜偏置引起的误差,经过实验发现F-theta聚焦透镜引起的误差比较明显。
鉴于此有学者用f-θ数学模型近似的模拟出了映射坐标值与数字控制量的数学关系,经过实验发现此模型不能很适合的表达两者的数学关系;又有学者用OEF图形算法模拟出两者的数学关系,而此模型的不足之处在于其数学关系式与实验结果仍然有比较大的出入,而且有些待测物理量用现在的测量方法很难测得,基于以上两个数学模型,我们创建出了一种新的数学模型,优化了现有理论模型,使得误差减小并更贴近实验结果,能更清楚的表达整个光路在添加F-theta聚焦透镜后的光路变化过程。
发明内容
本发明是针对fθ模型和OEF图形算法的局限性,无法对F-theta聚焦透镜后的光路变化而引起的误差进行合理分析的问题,提出了一种经F-theta聚焦透镜后引起的畸变误差补偿方法,针对由添加F-theta聚焦透镜引起的桶形畸变,研究出光束从路径BDFG映射到坐标平面后GH的距离,通过此距离和BD的方向矢量(通过振镜偏转角以及几何关系计算)计算出G点的映射坐标值,进而可得到计算机数字控制量与振镜偏转角、映射坐标值的数学关系。
本发明的技术方案为:一种经F-theta聚焦透镜后引起的畸变误差补偿方法,光束经过两个旋转的X振镜和Y振镜反射后,出射光进入F-theta聚焦透镜,X振镜和Y振镜由两个互相垂直的伺服电机带动旋转;Q-B-V为光束进入振镜的实际路径;S、U分别是映射坐标点V在Y、X轴上的映射;B、Q分别是Y、X振镜中心点;θx、θy分别表示X、Y振镜转动的角度;光束进入F-theta聚焦透镜后分别在F-theta聚焦透镜的左曲面和右曲面发生两次折射,B-D-F-G为光束进入-theta聚焦透镜的实际路径;在F-theta聚焦透镜中光路的高度对称性,设定二维光路中各个参数:入射光BD,出射光为FG,聚焦透镜横向中心线IC为Z轴;ID是左曲面曲率半径,N是ID延长线与聚焦透镜纵向中心线ON的交点;DE是入射点D到聚焦透镜横向中心线IC的距离;FM是出射点F到光屏轴GH的距离,GH为Y轴,像在光屏面XHY上;CF是右曲面曲率半径,P是CF延长线与聚焦透镜的交点;A是PD的延长线与聚焦透镜横向中心线IC的交点;FJ是出射点F到聚焦透镜横向中心线IC的距离;O、H分别是聚焦透镜和光屏轴GH的中心点,H在聚焦透镜横向中心线IC上;K、L为聚焦透镜的左曲面和右曲面中心点,K、L均在聚焦透镜横向中心线IC上;
假设:
经过两振镜后的出射光BD与Z轴的角度∠DBC=θ;
激光束BD进入F-theta聚焦透镜的入射角∠BDI=α1
左曲面折射后的角度∠FDN=β1
右曲面入射角角度∠DFP=β2
右曲面折射后的角度∠GFC=α2
左曲面入射点中心角∠DIK=φ1
右曲面入射点中心角∠FCK=φ2
像点(x,y)与H点连线与X轴正方向成的角度φ;
DE距离h;
GH距离R;
已知物理量:
f:OH距离,即F-theta聚焦透镜焦距;
a:两面振镜中心点距离QB;
b:Y振镜到F-theta聚焦透镜的距离BK;
t:F-theta聚焦透镜厚度KL;
n1:F-theta聚焦透镜外介质折射率;
n2:F-theta聚焦透镜介质折射率;
r1:F-theta聚焦透镜左曲面曲率半径ID;
r2:F-theta聚焦透镜右曲面曲率半径CF;
kx:控制X振镜的模拟电压与振镜转角的比例系数;
ky:控制Y振镜的模拟电压与振镜转角的比例系数;
Dx:控制X振镜转动的计算机数字控制量;
Dy:控制Y振镜转动的计算机数字控制量;
Vx:控制X振镜转动的模拟控制量;
Vy:控制Y振镜转动的模拟控制量;
未知参数的数学关系式:
θx=kxVx
θy=kyVy
θx=(Dx/65535-0.5)*10*kx*π/180;
θy=(Dy/65535-0.5)*10*ky*π/180;
Figure BDA0001813426390000031
φ=tan-1(cos(2θx)*sin(2θy)/sin(2θx));θ=cos-1(cos2θx*cos2θy);
φ1=sin-1(h/r1);
α1=θ-φ1
β1=sin-1(sin(α1)*n1/n2);
β2=sin-1(sin(β11)*(h*cot(β11)-h*cotθ+b+r2+t)/r2);
α2=sin-1(sin(β2*n2/n1));
φ2=β211
Figure BDA0001813426390000041
x=R*cosφ;
y=R*sinφ;
F-theta聚焦透镜横截面是圆形的,即有效像平面是圆形;有效像平面是以H为中心R为半径的圆形,而像(x,y)坐标值通过R和φ的关系可求得,通过数学关系式,得到Dx、Dy与(x,y)相关的数学表达式,计算机上输入Dx、Dy计算出映射坐标值(x,y)的真实值,跟理论值对比找出畸变差值,通过修改输入的Dx、Dy进行差值补偿。
本发明的有益效果在于:本发明经F-theta聚焦透镜后引起的畸变误差补偿方法,利用模型找到计算机数字控制量与映射坐标值的关系,实际运用时通过修改输入计算机数字控制量即可达到补偿的效果。
附图说明
图1是光束经过F-theta聚焦透镜的二维光路示意图;
图2是光束经过振镜的光路示意图。
具体实施方式
图1所示为光束经过F-theta聚焦透镜的二维光路示意图,由于光路的高度对称性可知三维依然成立,光束在经过F-theta聚焦透镜时会因为两个曲面发生两次折射,此时的路径B-D-F-G为光束的实际路径;A是PD的延长线与聚焦透镜横向中心线IC(Z轴)的交点;ID是左曲面曲率半径,N是ID延长线与聚焦透镜纵向中心线ON的交点;DE是入射点D到聚焦透镜横向中心线IC的距离;FM是出射点F到光屏GH轴(Y轴)的距离,像在光屏面XHY上;CF是右曲面曲率半径,P是CF延长线与聚焦透镜的交点;FJ是出射点F到聚焦透镜横向中心线IC的距离;O、H分别是聚焦透镜和光屏轴GH的中心点,H在聚焦透镜横向中心线IC上;K、L为聚焦透镜的左曲面和右曲面中心点,K、L均在聚焦透镜横向中心线IC上。
图2所示为光束经过振镜的光路示意图,Q-B-V即为光束的实际路径;S、U分别是映射坐标点V在Y、X轴上的映射;B、Q分别是Y、X振镜中心点;θx、θy分别表示X、Y振镜转动的角度。
光束经过两个旋转的X振镜和Y振镜反射后,出射光进入F-theta聚焦透镜,X振镜和Y振镜由两个互相垂直的伺服电机带动旋转;推导出BD(光束经过两面振镜后)的矢量方向,如上图2所示,光线从-i(-1,0,0)方向入射,经过X振镜(X振镜偏转θx,则光线偏转2θx)反射后光线的单位方向矢量变为(-sin(2θx),cos(2θx),0),再经过Y振镜(Y振镜偏转θy,则光线偏转2θy)反射后光线的单位方向矢量变为BD(-sin(2θx),cos(2θx)×sin(2θy),cos(2θx)×cos(2θy))。而由于光束经过F-theta聚焦透镜后,BDFG近似在同一个平面内,故G点方向矢量可知。假设:
(经过振镜出射光BD与Z轴的角度)∠DBC=θ;
(激光束BD进入F-theta聚焦透镜的入射角)∠BDI=α1
(左曲面折射后的角度)∠FDN=β1
(右曲面入射角角度)∠DFP=β2
(右曲面折射后的角度)∠GFC=α2
(左曲面入射点中心角)∠DIK=φ1
(右曲面入射点中心角)∠FCK=φ2
(像点(x,y)与H点连线与X轴正方向成的角度)φ;
DE距离h;
GH距离R;
待测物理量:(待测物理量通过测量或选用参数为已知)
f:OH距离(F-theta聚焦透镜焦距);
a:两面振镜中心点距离(QB);
b:Y振镜到F-theta聚焦透镜的距离(BK);
t:F-theta聚焦透镜厚度(KL);
n1:F-theta聚焦透镜外介质折射率;
n2:F-theta聚焦透镜介质折射率;
r1:F-theta聚焦透镜左曲面曲率半径ID;
r2:F-theta聚焦透镜右曲面曲率半径CF;
kx:控制X振镜的模拟电压与振镜转角的比例系数;
ky:控制Y振镜的模拟电压与振镜转角的比例系数;
Dx:控制X振镜转动的计算机数字控制量;
Dy:控制Y振镜转动的计算机数字控制量;
Vx:控制X振镜转动的模拟控制量;
Vy:控制Y振镜转动的模拟控制量;
未知参数的数学关系式:
θx=kxVx
θy=kyVy
θx=(Dx/65535-0.5)*10*kx*π/180;
θy=(Dy/65535-0.5)*10*ky*π/180;
Figure BDA0001813426390000061
φ=tan-1(cos(2θx)*sin(2θy)/sin(2θx));θ=cos-1(cos2θx*cos2θy);
φ1=sin-1(h/r1);
α1=θ-φ1
β1=sin-1(sin(α1)*n1/n2);
β2=sin-1(sin(β11)*(h*cot(β11)-h*cotθ+b+r2+t)/r2);α2=sin-1(sin(β2*n2/n1));
φ2=β211
Figure BDA0001813426390000071
x=R*cosφ;
y=R*sinφ。
由于图1给出的是平面分析图,但是由于F-theta聚焦透镜有高度的对称性,扩展到三维坐标系的横截面是圆形的,所以只需要分析图1所示的平面图即可扩展到三维,上式中最后的R代表的是GH的距离,图1中二维坐标系GH轴只显示在Y轴上,但是由于F-theta聚焦透镜的横截面是圆形的,即有效像平面是圆形,是以H为中心GH(R)为长度的圆形,而像(x,y)坐标值通过R和φ的关系(最后两个式子)可以求得,最终通过上面的所有关系式,我们会得到Dx、Dy与(x,y)相关的数学表达式,就实现了计算机的控制,例如:计算机上输入Dx、Dy我们就可以计算出映射坐标值(x,y)的真实值,跟理论值对比找出畸变差值,通过修改输入的Dx、Dy达到补偿的效果。

Claims (1)

1.一种经F-theta聚焦透镜后引起的畸变误差补偿方法,光束经过两个旋转的X振镜和Y振镜反射后,出射光进入F-theta聚焦透镜,X振镜和Y振镜由两个互相垂直的伺服电机带动旋转;其特征在于,Q-B-V为光束进入振镜的实际路径;S、U分别是映射坐标点V在Y、X轴上的映射;B、Q分别是Y、X振镜中心点;θx、θy分别表示X、Y振镜转动的角度;光束进入F-theta聚焦透镜后分别在F-theta聚焦透镜的左曲面和右曲面发生两次折射,B-D-F-G为光束进入F-theta聚焦透镜的实际路径;在F-theta聚焦透镜中光路的高度对称性,设定二维光路中各个参数:入射光BD,出射光为FG,聚焦透镜横向中心线IC为Z轴;ID是左曲面曲率半径,N是ID延长线与聚焦透镜纵向中心线ON的交点;DE是入射点D到聚焦透镜横向中心线IC的距离;FM是出射点F到光屏轴GH的距离,GH为Y轴,像在光屏面XHY上;CF是右曲面曲率半径,P是CF延长线与聚焦透镜的交点;A是PD的延长线与聚焦透镜横向中心线IC的交点;FJ是出射点F到聚焦透镜横向中心线IC的距离;O、H分别是聚焦透镜和光屏轴GH的中心点,H在聚焦透镜横向中心线IC上;K、L为聚焦透镜的左曲面和右曲面中心点,K、L均在聚焦透镜横向中心线IC上;
假设:
经过两振镜后的出射光BD与Z轴的角度∠DBC=θ;
激光束BD进入F-theta聚焦透镜的入射角∠BDI=α1
左曲面折射后的角度∠FDN=β1
右曲面入射角角度∠DFP=β2
右曲面折射后的角度∠GFC=α2
左曲面入射点中心角∠DIK=φ1
右曲面入射点中心角∠FCK=φ2
像点(x,y)与H点连线与X轴正方向成的角度φ;
DE距离h;
GH距离R;
已知物理量:
f:OH距离,即F-theta聚焦透镜焦距;
a:两面振镜中心点距离QB;
b:Y振镜到F-theta聚焦透镜的距离BK;
t:F-theta聚焦透镜厚度KL;
n1:F-theta聚焦透镜外介质折射率;
n2:F-theta聚焦透镜介质折射率;
r1:F-theta聚焦透镜左曲面曲率半径ID;
r2:F-theta聚焦透镜右曲面曲率半径CF;
kx:控制X振镜的模拟电压与振镜转角的比例系数;
ky:控制Y振镜的模拟电压与振镜转角的比例系数;
Dx:控制X振镜转动的计算机数字控制量;
Dy:控制Y振镜转动的计算机数字控制量;
Vx:控制X振镜转动的模拟控制量;
Vy:控制Y振镜转动的模拟控制量;
未知参数的数学关系式:
θx=kxVx
θy=kyVy
θx=(Dx/65535-0.5)*10*kx*π/180;
θy=(Dy/65535-0.5)*10*ky*π/180;
Figure FDA0002554612360000021
φ=tan-1(cos(2θx)*sin(2θy)/sin(2θx));
θ=cos-1(cos2θx*cos2θy);
φ1=sin-1(h/r1);
α1=θ-φ1
β1=sin-1(sin(α1)*n1/n2);
β2=sin-1(sin(β11)*(h*cot(β11)-h*cotθ+b+r2+t)/r2);
α2=sin-1(sin(β2*n2/n1));
φ2=β211
Figure FDA0002554612360000031
x=R*cosφ;
y=R*sinφ;
F-theta聚焦透镜横截面是圆形的,即有效像平面是圆形;有效像平面是以H为中心R为半径的圆形,而像(x,y)坐标值通过R和φ的关系可求得,通过数学关系式,得到Dx、Dy与(x,y)相关的数学表达式,计算机上输入Dx、Dy计算出映射坐标值(x,y)的真实值,跟理论值对比找出畸变差值,通过修改输入的Dx、Dy进行差值补偿。
CN201811130506.9A 2018-09-27 2018-09-27 经F-theta聚焦透镜后引起的畸变误差补偿方法 Active CN109188681B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811130506.9A CN109188681B (zh) 2018-09-27 2018-09-27 经F-theta聚焦透镜后引起的畸变误差补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811130506.9A CN109188681B (zh) 2018-09-27 2018-09-27 经F-theta聚焦透镜后引起的畸变误差补偿方法

Publications (2)

Publication Number Publication Date
CN109188681A CN109188681A (zh) 2019-01-11
CN109188681B true CN109188681B (zh) 2020-10-09

Family

ID=64907465

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811130506.9A Active CN109188681B (zh) 2018-09-27 2018-09-27 经F-theta聚焦透镜后引起的畸变误差补偿方法

Country Status (1)

Country Link
CN (1) CN109188681B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110076451B (zh) * 2019-06-04 2021-01-05 华霆(合肥)动力技术有限公司 激光加工装置及激光焦距补偿方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101487907B (zh) * 2009-03-05 2010-07-14 哈尔滨工业大学 连续浮雕结构微光学元件干法刻蚀图形传递误差补偿方法
CN102243373B (zh) * 2011-06-24 2014-06-25 南京英田光学工程有限公司 远距离特大出瞳直径透镜式检测仪光学系统
US9091843B1 (en) * 2014-03-16 2015-07-28 Hyperion Development, LLC Optical assembly for a wide field of view point action camera with low track length to focal length ratio
CN104678556A (zh) * 2015-02-06 2015-06-03 中国科学院理化技术研究所 一种透镜组条形象散光束实时整形与低阶像差补偿系统

Also Published As

Publication number Publication date
CN109188681A (zh) 2019-01-11

Similar Documents

Publication Publication Date Title
Clarke et al. The development of camera calibration methods and models
CN105629359B (zh) 一种高精度五棱镜的制作方法
CN103093676B (zh) 分光计数字化读数装置
CN109188681B (zh) 经F-theta聚焦透镜后引起的畸变误差补偿方法
CN106092515B (zh) 一种用法布里-珀罗标准具测量焦距和转角的方法
CN107991770B (zh) 用于确定圆柱形光学对象的折射率分布图的方法
CN106767675B (zh) 基于细光束的f-theta测量系统的优化方法
CN111044077A (zh) 星敏感器测量坐标系与基准立方镜坐标系间的标定方法
US3910704A (en) Compensating device for sighting instruments
CN105976374B (zh) 摄影测量相机的场相关畸变标定方法
Beutler Comparison of 2D and 3D measurements of aspheres with a tactile and optical sensor on one measuring instrument
Zuercher et al. Improved reflector for interferometric tracking in three dimensions
CN109708559A (zh) 一种基于隅角镜的光电自准直仪的测角方法
JPH0760221B2 (ja) 光走査装置
CN110361706B (zh) 一种雷达天线阵面法线方位角测量装置及方法
Chen et al. Fringe Projection Profilometry for Three-Dimensional Measurement of Aerospace Blades
RU2569072C2 (ru) Датчик угла поворота
CN112598583B (zh) 一种折射率模型畸变修正方法
JP2006003184A (ja) 面法線計測方法及びその装置
Bossemeyer et al. Calibration of focus-adjustable cameras for triangulation sensors by model-based camera parameter estimation: Polynomial estimation of the calibration parameters over the focal range of an objective lens
CN107830836B (zh) 一种双冗余psd自准直仪吊索摆角的测量方法
RU2773806C1 (ru) Способ измерения дисторсии в интерферограмме оптической асферической поверхности
CN111595271B (zh) 一种角度测量设备的标校方法及装置
Wright The formation of interference figures. a study of the phenomena exhibited by transparent inactive crystal plates in convergent polarized light
Sha et al. Study of Modified Offset Trajectory for Bonnet Polishing Based on Lifting Bonnet Method

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