CN110836634A - 可适应多种光束的四象限探测器标定方法 - Google Patents
可适应多种光束的四象限探测器标定方法 Download PDFInfo
- Publication number
- CN110836634A CN110836634A CN201910871989.6A CN201910871989A CN110836634A CN 110836634 A CN110836634 A CN 110836634A CN 201910871989 A CN201910871989 A CN 201910871989A CN 110836634 A CN110836634 A CN 110836634A
- Authority
- CN
- China
- Prior art keywords
- light
- quadrant
- detector
- noise
- spot
- 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.)
- Granted
Links
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
- G01B21/04—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness by measuring coordinates of points
- G01B21/042—Calibration or calibration artifacts
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明公开了一种可适应多种光束的四象限探测器标定方法,具体步骤为:根据入射光的能量分布规律,确定能量分布模型,根据能量分布模型确定白噪声背景光;根据四象限探测器的象限尺寸、盲区宽度、光敏响应度以及白噪声背景光,确定各象限的响应值;利用光斑位置测量算法计算由高斯光斑引起的各象限的输出电压值;根据由高斯光斑引起的各象限的输出电压值,利用四象限定位算法计算光斑质心的实际偏移量,根据光斑质心的实际偏移量确定探测器的测量误差。本发明将探测器盲区宽度加入模型进行计算,使计算精度更高。
Description
技术领域
本发明属于光电测量领域,具体为一种可适应多种光束的四象限探测器标定方法。
背景技术
四象限探测器(4-QD)是一种位置敏感器件,具有高灵敏度、高精度、计算简便等优点,被广泛应用于激光制导武器、激光雷达、空间光通信等激光跟踪领域。在对点源目标的跟踪与探测系统中,通常会在成像透镜的焦平面处放置四象限探测器,以测量焦平面处光斑质心的移动从而动态地解算出入射波前的倾斜量。信号光的光子噪声、背景光和探测器本身噪声使得四象限探测器存在跟踪误差;而四象限探测器在制造过程中,需要对探测器的光敏面进行象限分割,因此象限间不可避免地存在盲区(死区),当跟踪系统用于天文观测或星-地通信等领域时,由于接收到的光能量大部分需用于成像系统,因而探测器大都工作在光子技术模式下,此时盲区对焦平面光斑形状的改变和能量遮挡都会影响到四象限探测器的灵敏度和探测精度,同时在长波红外波段,由碲镉汞材料制作的四象限探测器由于材料面小,常需要用拼接的方式组成足够大的光敏面,此时的盲区面积则更加不可忽略;传统的四象限探测器解算方法的原理是建立在光斑的形状和半径是一个较理想的能量均匀分布的圆形光斑的基础上,其优势在于计算简便,便于应用在动态测量系统中。
实际情况中光斑的能量分布通常呈高斯分布,而且形状并不都是圆形,而在实际的测量情况中,光斑的形状和光强分布都会影响探测器的测量精度。传统的光斑能量分布模型仅能描述理想的圆形高斯光斑,而对于椭圆形的高斯光斑,传统的各象限的响应情况和测量算法的解算误差则不能够精准的描述。
发明内容
本发明的目的在于提出了一种可适应多种光束的四象限探测器标定的方法。
实现本发明的技术解决方案为:一种可适应多种光束的四象限探测器标定的方法,具体步骤为:一种可适应多种光束的四象限探测器标定方法,其特征在于,具体步骤为:
步骤1、根据入射光的能量分布规律,确定能量分布模型,根据能量分布模型确定白噪声背景光;
步骤2、根据四象限探测器的象限尺寸、盲区宽度、光敏响应度以及白噪声背景光,确定各象限的响应值;
步骤3、利用光斑位置测量算法计算由高斯光斑引起的各象限的输出电压值;
步骤4、根据由高斯光斑引起的各象限的输出电压值,利用四象限定位算法计算光斑质心的实际偏移量,将光斑质心的实际偏移量与光斑质心移动量实际值作差确定探测器的测量误差。
优选地,步骤1确定的能量分布模型为:
式中,EA(xc,yc)、EB(xc,yc)、EC(xc,yc)、ED(xc,yc)分别表示各个象限受到的光照强度,EElliptical(x,y)表示椭圆光斑光强,ENoise为背景光造成的白噪声光强,R为4-QD光敏面各象限的边长,G表示象限间盲区的半宽度。
优选地,椭圆光斑光强EElliptical(x,y)具体为:
式中A0为光束的中心光强,xc和yc分别为光斑中心相对探测器中心在x方向上的偏移和y方向上的偏移,rx和ry分别为椭圆光斑长轴和短轴的半径,θ为光斑绕质心旋转的角度,ENoise为背景光造成的白噪声光强。
优选地,步骤2确定的各象限的响应值具体为:
式中,UNoisei为各象限的背景光噪声和探测器本身噪声引起电压信号的强度,UNoise为背景光噪声和探测器本身噪声引起电压信号的强度,ρ为探测器光敏响应度。
优选地,步骤3利用光斑位置测量算法计算的由高斯光斑引起的各象限的输出电压值具体为:
式中,UA、UB、UC、UD分别表示由高斯光斑引起的各象限的输出电压值, A0为光束的中心光强,R为4-QD光敏面各象限的边长,G表示象限间盲区的半宽度,UNoise为背景光噪声和探测器本身噪声引起电压信号的强度。
优选地,步骤4中光斑质心的实际偏移量为:
式中,k为包涵了4-QD光响应度、电流-电压转换系数的比例系数
本发明与现有技术相比,其显著优点为:本发明将探测器盲区宽度加入模型进行计算,使计算精度更高;本发明的误差模型在实际应用中能够更好地对探测系统中的入射光进行拟合和描述,且对探测器的响应情况计算的更加精确,能够在多种场合和情况下对探测系统的测量误差进行计算。
下面结合附图对本发明做进一步详细的描述。
附图说明
图1为本发明的流程图。
图2为四象限探测器的感光表面的结构和光强分布示意图。
图3为光斑能量分布以及残差示意图。
图4为探测器各象限的输出电压值与各象限的实际输出电压值示意图。
具体实施方式
本发明对光斑能量分布函数进行优化,使其推广到能够适应带有倾斜的椭圆高斯分布的光斑,并将探测器盲区宽度加入模型中,从而能够对四象限探测器各象限的输出值进行计算,并分析定位算法的计算误差。
如图1、2所示,一种可适应多种光束的四象限探测器标定方法,具体步骤为:
步骤1、根据入射光的能量分布规律,确定能量分布模型,根据能量分布模型确定白噪声背景光,具体方法为:
在分析探测器光敏面的光能分布时,利用公式(1)对入射光的光强分布进行标定:
其中:xθ=(x-xc)cos θ+(y-yc)sin θ,yθ=-(x-xc)sin θ+ (y-yc)cos θ。EElliptical(x,y)表示椭圆光斑光强,A0为光束的中心光强,xc和yc分别为光斑中心相对探测器中心在水平(x)和垂直(y)方向上的偏移,R为4-QD 光敏面各象限的边长,rx和ry分别为椭圆光斑长轴和短轴的半径,θ为光斑绕质心旋转的角度,ENoise为背景光造成的白噪声光强。对于圆形高斯光斑,仍可用以上模型表示,此时长轴rx与短轴ry相等。
根据标定所得的光斑能量分布参数,结合四象限探测器的象限宽度,并且将探测器光敏面的盲区宽度考虑在内,则能够推导出四象限探测器的各象限受到的光照强度模型为:
其中,R为4-QD光敏面各象限的边长,G为象限间沟道的半宽度。
步骤2、根据四象限探测器的象限尺寸、盲区宽度、光敏响应度以及白噪声背景光,确定各象限的响应值,具体方法为:
在某些实施例中,应用场景中的背景光为白噪声背景光,则此时探测器各象限的噪声对探测器输出的响应情况为:
其中,UNoiseA,B,C,D分别为各象限的背景光噪声和探测器本身噪声引起电压信号的强度,UNoise为背景光噪声和探测器本身噪声引起电压信号的强度,ρ为探测器光敏响应度。
步骤3、利用光斑位置测量算法计算由高斯光斑引起的各象限的输出电压值,具体为:
根据步骤2得到的各象限的响应值,确定白噪声背景光造成的各象限噪声的电压在数值上基本相同。
此时由高斯光斑引起的各象限的输出电压值,即由高斯光斑引起的各象限的输出电压值为:
式中,UA、UB、UC、UD分别表示由高斯光斑引起的各象限的输出电压值, A0为光束的中心光强,R为4-QD光敏面各象限的边长,G表示象限间盲区的半宽度,UNoise为背景光噪声和探测器本身噪声引起电压信号的强度。
步骤4、根据由高斯光斑引起的各象限的输出电压值,利用四象限定位算法计算光斑质心的实际偏移量,将光斑质心的实际偏移量与光斑质心移动量实际值作差确定探测器的测量误差。具体为:
将由高斯光斑引起的各象限的输出电压值代入四象限定位算法,计算光斑质心的实际偏移量(xc,yc),根据光斑质心的实际偏移量(xc,yc)确定由高斯光斑引起的各象限的输出电压值(Δx,Δy)的误差,其中,光斑质心的实际偏移量(xc,yc)为:
式中,k为包涵了4-QD光响应度、电流-电压转换系数的比例系数。
本发明能够对造成四象限探测器出现测量误差的多种因素进行分析。因此,在实际应用中能够更好地对探测系统中的入射光进行拟合和描述,且对探测器的响应情况计算的更加精确,能够在多种场合和情况下对探测系统的测量误差进行计算。
实施例1
(1)光强分布模型的验证与对比分析:本实施例首先验证光强分布模型的正确性。采用两种激光器,其光源形状分别为椭圆形和圆形,各自利用式(1),结合光束质量分析仪所获得的光斑能量分布数据,如图3(a)、(d)、(g),拟合出光能分布模型的各项参数,将模型还原的光斑分布图像,如图3(b)、(e)、(i),与光斑图像对比,得出残差,如图3(c)、(f)、(j)。a)d)表示光源实际形成的椭圆形光斑,g)表示光源实际形成的圆形光斑,b)表示通过常用的圆形模型还原的光斑能量分布,e)i)表示通过本文的椭圆形模型还原的光斑能量分布,c)f)j) 分别为实际光斑与还原光斑图像的残差。由此可见,对于能量呈高斯分布的光斑,无论是椭圆或圆形高斯分布,都能够利用模型(式1)对其进行较准确地描述。由此说明,对于实际情况中常出现的椭圆高斯光斑,可由本发明得出更好的拟合结果。
(2)象限响应模型的验证与分析:在对光斑能量分布的各项参数进行标定之后,在实验的过程中,四象限探测器垂直于光轴沿水平或垂直方向移动,每次沿x轴移动0.2*R(R表示探测器半径)的距离,沿y轴移动0.1*R的距离,并记录各象限的输出电压,然后继续移动探测器,直到光斑质心基本移出探测器的光敏面为止。试验共分三组,每组利用变焦镜头调整光斑的能量分布参数,使红、绿、蓝三组所测量的光斑半径逐渐增大。
然后,将获取的光斑参数,结合探测器的固有参数(盲区宽度、象限宽度、光敏响应度和电路的电流-电压转换系数),并带入公式(4),计算出该探测器各象限的输出电压值,并将其与各象限的实际输出电压值进行对比,其结果分别如图4(a)、(b)、(c)、(d)所示,其中(a)、(b)、(c)、(d)分别表示各象限输出电压的测量值(即实际值)和当前模型及常用模型计算值,(e)、(f)分别表示x轴和y轴方向质心移动量的实际值和在本模型中的计算值,(g)、(h)分别表示x轴和y轴方向上,本文模型及常用模型的计算结果与质心移动量实际值的误差。
可以看出,本发明的计算值与实际测量值几乎完全吻合,相比常用模型误差更精确。
再将由模型得出的各象限输出值带入式(5)计算光斑质心沿x轴与y轴的移动量,并且也与实际质心的移动量进行对比,其结果分别如图4(e)、(f)。在x轴和y轴方向上,将本发明及常用模型的计算结果与质心移动量实际值作差,所得的误差大小如图4(g)、(h)所示。显然,本发明误差相对常用模型误差更小。
Claims (6)
1.一种可适应多种光束的四象限探测器标定方法,其特征在于,具体步骤为:
步骤1、根据入射光的能量分布规律,确定能量分布模型,根据能量分布模型确定白噪声背景光;
步骤2、根据四象限探测器的象限尺寸、盲区宽度、光敏响应度以及白噪声背景光,确定各象限的响应值;
步骤3、利用光斑位置测量算法计算由高斯光斑引起的各象限的输出电压值;
步骤4、根据由高斯光斑引起的各象限的输出电压值,利用四象限定位算法计算光斑质心的实际偏移量,根据光斑质心的实际偏移量确定探测器的测量误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910871989.6A CN110836634B (zh) | 2019-09-16 | 2019-09-16 | 可适应多种光束的四象限探测器标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910871989.6A CN110836634B (zh) | 2019-09-16 | 2019-09-16 | 可适应多种光束的四象限探测器标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110836634A true CN110836634A (zh) | 2020-02-25 |
CN110836634B CN110836634B (zh) | 2021-09-03 |
Family
ID=69574575
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910871989.6A Active CN110836634B (zh) | 2019-09-16 | 2019-09-16 | 可适应多种光束的四象限探测器标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110836634B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111272084A (zh) * | 2020-03-03 | 2020-06-12 | 丹阳钒曦光电科技有限公司 | 一种四象限探测器的标定方法 |
CN115218788A (zh) * | 2022-07-22 | 2022-10-21 | 天津大学 | 工作在波长大于1.1微米的硅四象限探测器及定位系统 |
CN117288246A (zh) * | 2023-11-27 | 2023-12-26 | 南昌大学 | 一种基于热电子效应的多象限探测器校正检测方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6438805A (en) * | 1987-08-05 | 1989-02-09 | Mitsubishi Electric Corp | Control method for working device |
JP2006172210A (ja) * | 2004-12-16 | 2006-06-29 | Matsushita Electric Works Ltd | 車両用距離画像センサおよびそれを用いる障害物監視装置 |
US20070080285A1 (en) * | 2005-09-08 | 2007-04-12 | Rafael-Armament Development Authority Ltd. | Arrangement for increasing the fill factor in a four-quadrant-type detector |
CN103630250A (zh) * | 2013-11-04 | 2014-03-12 | 中国科学院光电技术研究所 | 一种能自适应校正装调误差的四象限倾斜跟踪传感器 |
CN106100393A (zh) * | 2016-07-14 | 2016-11-09 | 中车大连电力牵引研发中心有限公司 | 单相四象限整流器死区补偿方法和装置 |
CN107356197A (zh) * | 2016-05-09 | 2017-11-17 | 南京理工大学 | 一种四象限光电探测器基于高斯分布的光斑定位方法 |
CN108240791A (zh) * | 2016-12-26 | 2018-07-03 | 南京理工大学 | 一种基于四象限探测器获取高斯光束高精度定位的方法 |
CN108681992A (zh) * | 2018-04-23 | 2018-10-19 | 南京理工大学 | 用于探测器阵列法测量激光光斑的图像插值算法 |
-
2019
- 2019-09-16 CN CN201910871989.6A patent/CN110836634B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6438805A (en) * | 1987-08-05 | 1989-02-09 | Mitsubishi Electric Corp | Control method for working device |
JP2006172210A (ja) * | 2004-12-16 | 2006-06-29 | Matsushita Electric Works Ltd | 車両用距離画像センサおよびそれを用いる障害物監視装置 |
US20070080285A1 (en) * | 2005-09-08 | 2007-04-12 | Rafael-Armament Development Authority Ltd. | Arrangement for increasing the fill factor in a four-quadrant-type detector |
CN103630250A (zh) * | 2013-11-04 | 2014-03-12 | 中国科学院光电技术研究所 | 一种能自适应校正装调误差的四象限倾斜跟踪传感器 |
CN107356197A (zh) * | 2016-05-09 | 2017-11-17 | 南京理工大学 | 一种四象限光电探测器基于高斯分布的光斑定位方法 |
CN106100393A (zh) * | 2016-07-14 | 2016-11-09 | 中车大连电力牵引研发中心有限公司 | 单相四象限整流器死区补偿方法和装置 |
CN108240791A (zh) * | 2016-12-26 | 2018-07-03 | 南京理工大学 | 一种基于四象限探测器获取高斯光束高精度定位的方法 |
CN108681992A (zh) * | 2018-04-23 | 2018-10-19 | 南京理工大学 | 用于探测器阵列法测量激光光斑的图像插值算法 |
Non-Patent Citations (4)
Title |
---|
QINGSONG GAO等: "Research on correction algorithm of laser positioning system based on four quadrant detector", 《PROCEEDINGS OF SPIE》 * |
唐彦琴等: "四象限探测器基于高斯分布的激光光斑中心定位算法", 《红外与激光工程》 * |
李生民等: "空间激光通信中四象限光电探测器环形光斑检测及误差补偿", 《中国激光》 * |
杨桂栓等: "死区对四象限探测器探测范围和灵敏度影响的研究", 《激光与光电子学进展》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111272084A (zh) * | 2020-03-03 | 2020-06-12 | 丹阳钒曦光电科技有限公司 | 一种四象限探测器的标定方法 |
CN115218788A (zh) * | 2022-07-22 | 2022-10-21 | 天津大学 | 工作在波长大于1.1微米的硅四象限探测器及定位系统 |
CN117288246A (zh) * | 2023-11-27 | 2023-12-26 | 南昌大学 | 一种基于热电子效应的多象限探测器校正检测方法及系统 |
CN117288246B (zh) * | 2023-11-27 | 2024-03-29 | 南昌大学 | 一种基于热电子效应的多象限探测器校正检测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110836634B (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110836634B (zh) | 可适应多种光束的四象限探测器标定方法 | |
US9574967B2 (en) | Wavefront measurement method, shape measurement method, optical element manufacturing method, optical apparatus manufacturing method, program, and wavefront measurement apparatus | |
CN108362469B (zh) | 基于压敏漆与光场相机的尺寸与表面压力测量方法与装置 | |
CN108121049B (zh) | 一种多谱段多通道遥感相机镜头的装调测试方法 | |
CN106404352B (zh) | 一种大视场望远镜光学系统畸变与场曲的测量方法 | |
CN105486489B (zh) | 电视成像系统调制传递函数测试装置及方法 | |
CN106644105B (zh) | 基于双螺旋点扩散函数的波前传感器、探测方法及系统 | |
CN106768882B (zh) | 一种基于夏克-哈特曼波前传感器的光学系统畸变测量方法 | |
CN110662020B (zh) | 一种基于自准直原理的传函测试系统及方法 | |
CN109323851B (zh) | 一种太赫兹焦平面响应率及响应率不均匀性测试系统及方法 | |
CN101672726B (zh) | 空间光通信终端通信探测器定位测试装置及方法 | |
CN105913414B (zh) | 一种红外摄像头视觉系统的标定装置及标定方法 | |
CN103226243A (zh) | 一种校正自适应光学系统非共光路误差的方法和系统 | |
CN103676487A (zh) | 一种工件高度测量装置及其校正方法 | |
CN107764518B (zh) | 一种光学镜头焦距测量装置及方法 | |
CN106705888B (zh) | 干涉检测中的ccd坐标系与镜面坐标系非线性关系标定方法 | |
CN103134443B (zh) | 一种大口径大径厚比反射镜面形自准直检测装置及方法 | |
CN111457911A (zh) | 一种基于偏振二维残差信息的仿生偏振罗盘标定方法 | |
CN114235352B (zh) | 一种四相机实时偏振成像系统的相对偏振角度测试方法 | |
CN111220088B (zh) | 测量系统和方法 | |
US9442006B2 (en) | Method and apparatus for measuring the shape of a wave-front of an optical radiation field | |
US10378995B2 (en) | Wavefront testing system | |
CN111272084A (zh) | 一种四象限探测器的标定方法 | |
CN104764588A (zh) | 单脉冲激光动态焦斑位置测量装置及测量方法 | |
CN108731808B (zh) | Ims型快照式成像光谱仪子孔径中心位置定标方法及装置 |
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 |