激光雷达的测距特性及其测距精度研究_姜海娇
第38卷 第5期2011年5月
中 国 激 光
Vol.38,No.5
,Ma2011激光雷达的测距特性及其测距精度研究
姜海娇 来建成* 王春勇 李振华
()南京理工大学理学院,江苏南京210094
摘要 通过实验研究了阈噪比(不变时信噪比(对激光雷达测距特性的影响,结果表明,信噪比远高于TNR)SNR)阈噪比时,激光雷达测距数据的离散分布服从正态分布,且随着信噪比的降低,测距数据离散分布逐渐偏离正态分呈现上升沿陡下降沿缓的形式;当信噪比接近阈噪比时,测距数据的离散分布又逐渐趋向于正态分布,但此时布,
采用置信概率6激光雷达的漏测率迅速增加。在此基础上,8.2%对应置信区间的一半来衡量激光雷达的测距精度,并给出了测距精度计算方法,它能更准确地表征激光雷达测距数据的密集度,具有更好的普适性。关键词 激光技术;激光雷达;测距精度;置信概率;置信区间;信噪比
:/中图分类号 TN958.98 文献标识码 A doi10.3788CJL201138.0514001
ResearchonRaninProertofLaserRadarandItsRaneAccurac ggpygy
JianHaiiaoaiJianchenanChunoniZhenhua L W L gjggyg
(SchooloScience,NaninUniversitoScienceand TechnoloNaninJiansu210094,China) fjgyfgy,jg,g Abstractheinfluencesofsinaltonoiseratio(SNR)ontheraninroertoflaserradarwithconstant T gggppy
thresholdtonoiseratio(TNR)arestudiedcomrehensivelbexeriments.TheresultsindicatethatwhenSNRis pyyp hiherthanTNR,thedistributionoflaserradar′sranindataobestheGaussiandistribution.Andwiththemuch gggy
,decreaseoftheSNR,thedistributionradualldeviatesfromtheGaussiandistributionanditsrisinedebecomes gygg steeerandthefallinedebecomesmoreslowl.WhiletheSNRisaroachintheTNR,thedistributionreturnsto pggyppg
,,theGaussiandistributionbutthemissinrateincreasesraidl.Basedontheexerimentsthehalfoftheconfidence gpyp intervalrobabilitwithconfidenceof68.2%isalsointroducedtoevaluatetheraneaccuracoflaserradar.The pygy resentedcorresondinmethodforcalculatintheraneaccuracisandusedtoevaluatetheraneaccuracoflaser ppgggygy underdifferentSNRs.radar
;;;;;Kewords lasertechniuelaserradarraneaccuracconfidencerobabilitconfidenceintervalsinalto qgypygy noiserate
OCIScodes80.3640;280.3400;280.5600 2
1 引 言
相对于常规的固体激光雷达,激光二极管(LD)雷达在高重复频率、低功耗和小型化等方面有着不可替代的优势,是集成距离传感、立体成像或者视觉在文物复原、精确制导、机器视觉、传感的重要途径,
靶场测试、车辆防撞、环境监测等领域有着极为重要
1~7]
。测距型固体激光雷达以调Q的脉的应用价值[
阈值鉴别和高精度的距离计时实现目标距离的测量过程。调Q激光脉冲具有兆瓦量级的峰值功率和纳秒量级的脉冲宽度,激光雷达测距数据的离散分布主要受噪声分布形式的影响,呈现近似正态分布形式。为此,通常用方差来衡量激光雷达的测距精
[]
度。Johnson等8研究了散粒噪声限的激光雷达系
统测距精度边界问题,给出了泊松噪声分布假设下
[]测距精度的Cramerao下界;Steinvall9运用数值-R
冲激光为光源,以飞行时间测距原理为基础,通过恒
;收到修改稿日期:2010113020110214 收稿日期:----
)基金项目:江苏省普通高校研究生科研创新计划(资助课题。AD20320
,:_作者简介:姜海娇(男,博士研究生,主要从事激光雷达方面的研究。E-m1985—)ailxiaofeixia900@163.com,导师简介:李振华(男,博士,教授,主要从事激光探测与目标定位等方面的研究。1964—):E-maillizhenhuaail.nust.edu.cn@mj
:aillaiianchenail.nust.edu.cn *通信联系人。E-m@mjgj
40
中国激光
仿真的方法给出了目标形状、探测器噪声以及湍流
10]
姚金良等[提出了运对激光雷达测距精度的影响;
的需要设定。因此,本文选用激光二极管雷达作为研究不同信噪比条件下激光雷达的测距实验对象,
特性及其测距精度指标。
图1给出实验所用激光二极管雷达的原理框图,它由激光发射系统、激光回波接收放大系统、脉冲鉴别与距离计时电路等部分组成。其中发射天线焦距2口径1接收天线焦距6口0mm、6mm,0mm、径2鉴别电平1探测器为C4mm,650mV,30950光电二极管。其测距过程的工作原理为:在时序电路脉冲激光二极管发射一个高峰值功率的的控制下,
并启动计时电路开始计时;激光脉冲经光激光脉冲,
学天线准直后,射向待测目标;目标散射照在其上的激光脉冲,形成后向散射回波,并被激光雷达的接收天线接收,聚焦到光电二极管的光敏面上,将光信号转换成电信号,形成激光回波的电脉冲信号。电信号经前置放大电路放大后,进行恒阈值的回波鉴别;当鉴别到回波信号时,立即停止计时电路,并根据计时结果计算目标距离。最后,通过串口实时将测距进行统计和分析
。数据传给计算机,
用图像配准叠加提高成像激光雷达测距精度的方法。然而,激光二极管雷达由于峰值功率较低,导致信噪比(较低。此时,回波波形对激光雷达测SNR)距特性的影响不可忽略,测距数据的离散分布将明显偏离正态分布,以方差为基础定义的测距精度概念不再适用,至今未见有对上述现象进行深入研究的报道。为此,本文在实验上研究了阈噪比(TNR)相同信噪比不同时激光雷达测距的统计特性,并在引入置信概率6此基础上,8.2%对应的置信区间一半来衡量激光雷达的测距精度的新思路,有效解决了特殊条件下激光雷达测距精度的评价问题,给出了离散分布距离数据测距精度的计算方法,并将其用于激光雷达测距精度随信噪比变化关系的研究中。
2 激光雷达测距的统计特性
脉冲激光二极管很容易通过驱动电路的参数设置改变输出脉冲功率,能够在较小的空间范围内(几十到几百米)使激光雷达的信噪比可根据实验测量
图1激光二极管雷达的实验原理框图Fi.1SchematicdiaramofLDradarexeriment ggp
0m远的一堵白墙为测试目标, 实验选择约7
通过标准白板的比照测试,白墙的半球反射率约为大气的能见度约为1白墙表面90%。实验时,0km,
2
/(·m·μ。脉冲的背景光辐照度约为800Wsrm)激光二极管的额定输出峰值功率为7通过改变0W,激光二极管的驱动电压,可以直接控制它的输出功率,达到改变激光雷达信噪比的目的。
当激光二极管以额定功率20W输出激光脉冲时,激光二极管雷达工作时的信噪比为3激光0.4,回波波形如图2所示。为了尽可能减小电路噪声和在大量实验的基础上,将阈值电平设置干扰的影响,
为1此时激光二极管雷达的工作状态较为650mV,
理想。此时的阈噪比为1错测率接近于0,虚6.72,对同一目标进行1警率极低。在此条件下,0万次定点测距,统计测距数据的离散分布直方图,并做归一化处理。图3给出了对应测距数据的统计分布曲线,纵坐标为10万次测量落在该距离上的百分数。从图中数据可以看到,测距数据的离散分布与正态分布吻合较好,吻合程度达到了99.99%。
改变激光二极管的驱动电路参数,调节激光二极管的输出功率来改变激光雷达测距的信噪比。图4给出了8种不同信噪比条件下激光雷达测距数并据的离散分布曲线。统计样本量同样为10万个,将曲线做归一化处理。
40
姜海娇等:
激光雷达的测距特性及其测距精度研究
图2发射功率为20W时的激光回波和噪声信号Fi.2Laserechowaveformandnoisesinalwhen gg
oweremittinis20
W pg
图3发射功率为20W时的测距数据统计分布Fi.3Ranedistributionwhenemittinower gggp
20
Wis
图4不同信噪比时激光雷达测距数据的离散分布曲线Fi.4LaserradarranindistributioncurveswithdifferentSNRs ggg
40
中国激光
对比图4中不同信噪比条件下的激光雷达测距
可以看到,信噪比远高于阈噪比数据离散分布曲线,
时,测距数据的离散分布较好地服从正态分布;随着))信噪比的逐渐降低,如图4(所示,测距数据ad~(的离散分布范围明显增大,测距数据的密集度逐渐变差,且分布形式逐渐偏离正态分布。具体表现为测距统计分布的不对称性,离散分布的上升沿较陡,且与正态分布吻合较好,下降沿却比上升沿缓;当信噪比接近阈噪比值时,激光雷达的漏测率显著增加,此时测距数据的离散分布又逐渐趋向于正态分布,且测距数据的离散分布范围又逐渐变小,在剔除漏测数据的情况下,测距数据的密集度有所改善。
激光雷达测距数据分布的上述统计特性直接由激光雷达接收机的噪声特性和激光回波的波形特征所决定。在信噪比远高于阈噪比时,噪声是引起测距数据离散分布的主导因素。信噪比远高于阈噪比时,激光脉冲回波波形的上升沿较陡,噪声引起阈值鉴别点处回波信号的时间抖动较小,这使得测距数在信噪比远高于阈据的离散分布范围较小。此外,
阈值点附近信号上升沿可由直线近似描述。噪比时,
在这种情况下,测距统计分布的结果主要由噪声的分布形式决定,即满足正态分布形式。这是以往激光测距研究所遇到的一般情况,在这种情况下测距数据是服从正态分布的,这是用标准差来衡量激光雷达测距精度的理论基础。
随着信噪比的逐渐降低,阈值点附近信号上升沿的波形也将显著影响激光雷达测距数据的离散分这将导致测距数据的分布形式偏离正态分布形式,
布。由于激光回波上升沿较缓,噪声能够引起阈值点的时间抖动明显增大,这使得测距数据的分布范围变宽。测距数据分布的不对称性(偏离正态分布)是因为阈值点上下斜率的绝对值明显不同,而且阈值点以上的上升沿形状决定了测距数据离散分布的后半截,阈值点以下上升沿形状决定了测距数据离散分布前半截。对于洛伦兹信号的上升沿而言,阈而值鉴别点以上位置的斜率明显小于阈值点斜率,阈值鉴别点以下位置的斜率与阈值点的斜率较为接近,这使得测距数据分布的上升沿接近于正态分布,而下降沿明显缓于上升沿。以往研究没有足够重视这种情况的存在,主要是因为传统的固体激光雷达采用的是单脉冲测距方式,这种方式要求激光雷达才能保证高的有效数据率;的信噪比远高于阈噪比,
但是激光二极管雷达(测距机)采用了多脉冲相关测距,这极大放宽了对信噪比的要求,使得激光二极管
雷达工作的信噪比有可能接近阈噪比,且具有较高因此讨论信噪比接近阈噪比时激光的有效数据率,
雷达的测距特性同样具有重要的意义。
当信噪比接近或低于阈噪比时,只有激光回波峰值附近很小区域的信号才有可能达到阈值电平位置,这时激光二极管雷达漏测率将明显上升。此外,由于激光回波峰值附近上升沿与下降沿具有绝对值几乎相等的斜率,这使得测距数据的离散分布又趋向于对称分布,同时在峰值点附近小区域内斜率几且变化非常缓慢,这使得测距数据的离散分乎为0,
且在剔除漏测数据后,测距数布又趋近于正态分布,
据的分布范围又有所减小。在这种情况下,漏测率很高,单脉冲测距方式是无法进行的,但对于能高重复频率工作的激光二极管雷达(几十或几百千赫兹)而言,通常可以通过多脉冲相关探测的方式将漏测率降到可实用的要求之下,在激光二极管雷达(测距机)的工作中通常会遇到这种情况。
此外,从图4还可以看到,随着信噪比的下降,激光雷达的测距值呈现单调上升的趋势。这是恒阈信号A的幅值鉴别方式的固有特征。如图5所示,
值明显高于信号B,恒阈值鉴别得到的时间点ta和这正好与图4给出的结果一致
。ttb的关系为ta<b,
图5激光雷达恒阈值鉴别方式原理图Fi.5Schematicdiaramoflaserradarin gg
thresholddiscrimination
3 不同信噪比条件下激光雷达的测距
精度
“在现有激光雷达测距性能评价体系中,测距精度”定义为测量的目标距离与测距均值的相对误差
11~13]
。此定义要求测距数分布,一般按标准差计算[
据的离散分布必须服从正态分布。从2节的实验结果可以看到,在信噪比远高于阈噪比时,测距数据离散分布服从正态分布是比较容易满足的,但在信噪比较低时,测距数据分布将明显偏离正态分布,上述定义不再适用。要想在不同信噪比条件下,全面评
40
姜海娇等:激光雷达的测距特性及其测距精度研究
定激光雷达的测距精度,必须拓展激光雷达测距精度的内涵。
事实上,对于测距数据统计偏离正态分布的情况,可以将“方差”的概念推广为置信概率68.2%所对应的置信区间一半。置信概率68.2%所对应的置信区间[定义为a,b]
b
应的数据分布范围正好是峰值附近与置信概率)如图6(所示。对68.2%所对应的置信区间一半,a
于测距数据分布明显偏离正态分布的情况,同样可以用置信概率68.2%对应的置信区间的一半来衡量测距数据的离散分布情况。
基于上述考虑,选用置信概率68.2%所对应的置信区间的一半代替方差来衡量激光雷达的“测距。当测距数据服从正态分布时,精度”该置信区间的一半正好与方差相一致,与现行通用的“测距精度”概念相兼容
。
F(x)=f(x)dx=0.682
,()1a
a)=f(b)(式中f(为测距数据离散分布的概率密度函数。x)
“在测距数据服从正态分布时,所对1倍方差”
∫
))图6正态(和非正态(分布情况下置信概率6ab8.2%所对应的置信区间
))Fi.6Confidenceintervalwithconfidenceof68.2%ofnormal(aandasmmetrical(bdistributionsrobabilit gypy
激光雷达测距数据是以测距分辨率为 实际上,
间隔的离散数据点,置信概率68.2%对应的置信区间的边界正好落在测距数据点上的概率极低,因此如何确定置信概率68.2%对应的置信区间的上下边界是准确计算激光雷达测距精度的关键所在。为此,假设每个距离值上是以测距分辨率为宽度的对分布区间的积分值为该点的归一化测称均匀分布,
)距数目,分别如图7(和(ab)所示。以测距统计峰值为中心,逐步扩大置信区间的上下界,直至累加概率为6假如左右两侧同时遇到测距数据8.2%为止;点,则以该数据点为中心均匀增大或减小置信区间的边界,直至测距数据的概率和为6此时对8.2%,应的距离值恰好为置信区间的上下边界
。
))图7(对称分布与(非对称分布的离散测距置信区间ab
))Fi.7Confidenceintervalofdiscreteranedateof(anormaland(basmmetricaldistributions ggy
置信概率 图8给出了在阈噪比不变的情况下,68.2%所对应的置信区间随信噪比的变化曲线。从曲线的分布可以看到,在信噪比远高于阈噪比的情况下,置信区间相对于测距数据的峰值位置对称分布,且置信区间较小,测距精度较高;当信噪比逐渐
减小时,置信区间相对于测距数据的峰值位置分布不再对称,正的偏差明显大于负的偏差,同时置信区间较大,测距数据的密集度明显变差。这恰好与低信噪比条件下,激光雷达测距数据统计分布上升沿陡而下降沿缓的特征相吻合;当信噪比进一步接近
40
中国激光
阈噪比时,测距数据的置信区间达到极大值后,又开这与2节的讨论相一致
。始逐渐减小,
4 结 论
建立了输出功率可调的激光二极管雷达测距特性研究的实验装置,并在不同信噪比条件下进行了大量的定点测距实验和分析比较。结果表明,当信噪比远高于阈噪比时,测距数据的离散分布服从正态分布;随着信噪比的降低,测距数据离散分布逐渐呈现上升沿陡下降沿缓的形式;当信偏离正态分布,
噪比接近阈噪比时,测距数据的离散分布又趋向于但漏测率迅速增加。在此基础上,提出了正态分布,
以68.2%置信概率对应置信区间的一半衡量激光雷达测距精度的新方法,给出了离散数据点测距精度的计算方法,并将其应用于实际激光雷达系统的测距精度评价。结果表明,当信噪比逐渐减小至阈测距精度先变差到极值然后又会出现一定噪比时,
程度的改善,以置信区间表征的测距精度能很好地具有很强的普适反映测距数据离散分布的密集度,性。
参
考
文
献
图88.2%置信概率对应的置信区间随信噪比的变化 6Fi.8Relationshibetweenconfidenceintervalwith gp
confidenceof68.2%andSNRrobabilitie p
在此基础上,根据“测距精度”的新定义,计算了不同信噪比条件下激光雷达的测距精度,如表1所示,它们分别与图4中不同信噪比条件下激光雷达测距数据的离散分布相对应。从表中数据可以看到,当信噪比为5激光雷达的测距精度(置0.66时,为0.信概率68.2%对应的置信区间的一半)0607,此时如果用方差衡量激光雷达的测距精度,则为两者吻合较好。由此可见,激光雷达测距精0.0658,
度在信噪比远高于阈噪比时具有很好的一致性,这随也再次验证了激光雷达测距精度定义的正确性;着信噪比的降低,激光雷达的测距精度逐渐变差,当激光雷达的测距精度为信噪比下降到19.76时,测距数据的密集度最差;随着信噪比进一步0.2170,
接近阈噪比,测距数据的密集性又有一定程度的改善。
表1不同信噪比条件下激光雷达的测距精度Table1LaserradarraneaccuracofdifferentSNRs gy
SNR 50.66 40.53 30.40 28.37 26.34 24.32 21.28 22.29 20.77 19.76 18.74 17.73 16.72
Raneaccuracm gy/
0.06070.08670.11240.11260.13440.16300.13130.17170.20380.21700.15900.15970.1553
,F.B,P.B1J.A.Beraldinlaisoulanerealworldetal..R g
modellinthrouhhihresolutiondiital3Dimainofobects ggggggj
:,andstructures[J].2000,55(4)ISPRS J.Photoramm. g230~250
,J,C.A.A.Beraldintzeni.Hihaccurac32G.Guidi--Dgy
:tmodelinofulturaleritaeheiitizinfonatello′s c h d dgggg o
"m"[addalenaJ].IEEE TransactionsonImaeProcessin gg,
():2004,133370~380
,,L.M.KlasenC.A.Gronwall3Dlaser3O.K.Steinvalletal.. :sensinatFOIoverviewandasstemC].ersective[SPIE, gypp
2004,5412:294~309
,K.H,K.Kanawaisemae4K.Saneoshietal..3 -D Igy
[reconitionsstemfordriveassistC].IntellientVehicles′93 gyg
,,Smosium199360~65yp
,eiselber.Gatedviewinandhihaccurac5J.BusckH.H -gggy
[],(),threedimensionallaserradarJ.0044324:Al.Ot.2- ppp4705~4710
,W,KTian′ananYuncaionLininChaoticlidar6Gonetal.. gggqg
[]forautomotivecollisionwarninsstemJ.ChineseJ.Lasers, gy
():2009,3692426~2430,,]龚天安王云才孔令琴等.面向汽车防撞的混沌激光雷达[J.
,():中国激光,20093692426~2430,,7DonYunsheniu WeniniuJianuolicationetal..A ggLqgLgpp
[]studoflidarinurbantrafficollutionJ.Acta OticaSinica, yp p
():2010,302315~320,,董云升刘文清刘建国等.激光雷达在城市交通污染中应用研
]():究[J.光学学报,2010,302315~320
,,L.NicholsP.GattRane8S.E.JohnsonT.recisionetal.. gp
[],,ofdirectdetectionlaserradarsstemsC.0045412:SPIE2 - y
72~869O.Steinvall.Effectsoftaretshaeandreflectiononlaserradar gp
[],():,crosssectionsJ.00039244381~4391Al.Ot.2 ppp,,10YaoinliananuiminZhaniudamaeetal..I J H gYgg X
reistrationandsuerositionforimrovinraninaccuracof gpppgggy
[],(),imainlaserradarJ.C010376:hineseJ.Lasers2 gg
1613~1617
40
姜海娇等:激光雷达的测距特性及其测距精度研究
姚金良,严惠民,张秀达等.一种应用图像配准叠加提高成像激
:光雷达测距精度的方法[J].中国激光,2010,37(6)
1613~1617
,T11CommissionofScienceechnoloandIndustrforNational gyy
Defense.Generalecificationsorolidulsedaser S f S P Lp
[]:,RanefindersGJB132491S.BeiinhinaStandardPress - gjgC
1991
,T12CommissionofScienceechnoloandIndustrforNational gyy
eriesorulsedaseranefindersDefense.Parameter S f P L Rg
[]:,GJB274096S.BeiinhinaStandardPress1996- jgC
,13CommissionofScienceTechnoloandIndustrforNational gyy
Defense.PerformanceTestMethodforPulseLaserRanefinder g
[]:,GJB2241.22006S.BeiinhinaStandardPress2006- jgC
檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱檱
《光电产品与资讯》第1中国的大光学产业”专刊征稿启事1期“
光电产品与资讯》杂志是中国激光杂志社2010年隆重推出的一本以产品和市场资讯为主的信息类期刊。该刊立足光电 《
领域,全年12期分别以不同的主题全面介绍技术与市场的最新进展。凭借中国激光杂志社的深厚学术底蕴及强大的采编力量,深入浅出地剖析光电行业事件,详细介绍光电产品的实际应用,在科研与产业之间架起一座互通的桥梁。并致力于成为选购产品的重要参考,企业推广产品的首选平台。读者了解市场、
。内容涵盖国内各地区的光学产业发展规划,我们2中国的大光学产业”专刊(英文版)全面展示最新011年第11期策划“光电技术的发展趋势、优秀产品的应用案例及实用指南等,特向广大读者征稿。该期杂志将于2012年1月份参展美国西部光电展。
征稿范围:光学、光电子学领域
)稿件类型:国内外知名专家、企业高管撰写的最新光电技术或市场发展趋势;1
)光电产品及技术的行业应用案例分析;2
)产品使用心得、经验技巧、故障排除等实用指南类型的文章。3
字 数:2500~4000字专刊截稿日期:2011年8月31日
稿件格式要求:投稿文体为英文,请使用w文章最后列出所有作者姓名、单位名称、职务或者职称、联系电话及ord格式,、邮寄地址和邮编。Eail-m
,请将稿件发至o邮件主题请标明“中国的大光学产业”专刊投稿。有任何问题请咨询李洪丹编辑,电话:eniom.ac.cn@sp。02169918166-
SPECIALISSUECALLFORPAPER
”2011ISSUE11“China′sBiOticalIndustr gpy
,EProduct&Newsisanew maazineofChineseLaserPress.It′sfeaturedbotoelectronicsnewsandmarketroducts O gypp
,,reortsandthecirculationcoversOErelatedcomaniesuniversitiesandinstitutes.ItaimstoaccelerateChinaOEresearch - pp
,v,,t,pandindustrialdeveloment.Themaincolumnsincludenewsiewointinterview,focusechnoloexchaneroduct ppgyg
,,etc.OEPNrovidesachanneltoknowtheOEmarketandbudevicesanditisalsoafavorablelatformofinformation pyp
romotinroductsforcomanies. pgpp
,ThetoicofNO.11isChina′sBiOticalIndustr.ItwillreortthedevelomentofoticalindustrinChinaalication pgpypppypp
analsisofoticalroductsandsoon.ThemaazinewillbedistributedinPhotonicsWest2012.Wehoeoucanwritesome ppgpyyarticlestous.
:,Toicsincludeoticsotoelectronics ppp
:)Te1latestotoelectronicstechnoloormarketdevelomenttrendinChinawrittenbexertsorenterriseeneral pgypyppgyp
;manaersg
);2industralicationanalsisofotoelectronicsroductsandtechnolo yppyppgy ),,,3racticaltesofarticlese..ourexeriencecommentsabouttheroducttroubleremoval. pypgypp
Articlelenth:2500~4000words g
:Deadline2011.8.31
:,,,Form:leasewritetheaerinEnlishandsendittooeniom.ac.cn.Theauthor′snamecomanositionortitle @s ppppgppy
,,,:telehoneEaildeliveraddressandostcodeshouldbeiven.Ifouhaveanuestionleasecall02169918166.-m -pypgyyqp
40