CN108333592B - 基于数字恒虚警的星载激光测距仪探测及反馈方法和系统 - Google Patents
基于数字恒虚警的星载激光测距仪探测及反馈方法和系统 Download PDFInfo
- Publication number
- CN108333592B CN108333592B CN201810072848.3A CN201810072848A CN108333592B CN 108333592 B CN108333592 B CN 108333592B CN 201810072848 A CN201810072848 A CN 201810072848A CN 108333592 B CN108333592 B CN 108333592B
- Authority
- CN
- China
- Prior art keywords
- noise
- value
- apd
- mean square
- root mean
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/4802—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明一种基于数字恒虚警的星载激光测距仪探测及反馈方法和系统,对特定区域内接收机噪声的多次测量统计获得均值和均方根值,根据获取的噪声统计结果,对雪崩二极管的偏压进行调整,使其随环境、温度等的变化而变化,实现虚警率的恒定控制;同时对地形参数快速变化时,采用调整回波预鉴别阈值来控制虚警恒定。该反馈算法克服了星载激光测距仪由于太阳高度角的变化引起输出噪声剧烈变化的影响;同时对噪声多次累加统计,消除了地形起伏较大而引入奇异点,避免奇异点对后续测量造成影响;识别云层散射虚假回波并滤除,以上均通过数字控制方法,满足系统探测概率要求情况下,降低了在轨回波的虚警率,快速并精准获取测距值。
Description
技术领域
本发明属于星载激光雷达光学遥感信号探测领域,涉及一种基于数字恒虚警的探测及在轨反馈方法和系统,提高卫星平台高速运行产生太阳高度角的变化适应性;提高激光测距仪对粗糙度起伏大、反射率较大的复杂地形的适应性,并且提高识别云层散射虚假回波的能力,满足虚警率和探测概率的要求。
背景技术
星载激光测距仪通常选用雪崩二极管(APD)作为接收机焦面探测器。APD对背景噪声极为敏感,输出的信号和噪声随APD倍增因子的变化而变化,倍增因子存在着一个最佳值。在此最佳倍增下,系统获得最佳信噪比。由此控制APD偏压进而影响输出的噪声大小,使其工作在最佳倍增因子状态。
星载激光测距仪系统虚警概率随APD输出噪声增大而降低,为降低虚警率,同时保证较高的探测概率,传统方法是星载激光测距仪采用模拟探测体制,一方面是运用回波信号送入比较器,只针对回波阈值进行动态调整,根据噪声事件计数来动态调整回波鉴别阈值,但是随着鉴别阈值的调整,会滤除激光回波携带的地表信息;另一方面是针对该种缺陷,在调整鉴别阈值的同时,还可通过动态调整APD的偏置高压,探测器的倍增因子随之降低,这样会对探测器灵敏度产生影响,在降低背景噪声灵敏度的同时,对实际目标的探测能力也随之下降。模拟电路调整方式,根据统计APD探测器噪声事件个数进行动态实时调整APD偏压,但是在实际工作中,星载平台由于速度相对较大,导致地面激光足印间的间距较大,这样出现一个噪声奇异点就会影响到下一次的APD偏压设定,进而影响测距性能。
星载激光测距仪随卫星平台运行过程中,太阳高度角的变化会引起APD探测器输出散粒噪声电流的变化,如果噪声过大将大量的地形信息湮没,地面也无法提取。另一方面,云层散射太阳光造成虚假回波的出现,需有效识别噪声并进行滤除。此外,对于陆地、森林等粗糙度较大的地形,会出现多个回波现象;对于地表反射率不同的地形,根据当时的信噪比精准的对虚警率进行调整,亟待提出一种方法以满足在轨从噪声中提取出真实回波信号,降低在轨虚警率的需求。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供一种基于数字恒虚警的探测及在轨反馈方法及系统,解决了星载测距仪易受卫星平台速度大、抖动大的影响,造成回波信号不易提取,云层散射造成虚假回波,对粗糙度起伏大、反射率较大的复杂地形虚警率较高等问题。
本发明的技术方案是:一种基于数字恒虚警的星载激光测距仪探测及反馈方法,步骤如下:
(1)读取并执行控制算法函数;
(2)根据星载激光测距仪系统要求的虚警率与探测概率计算地面回波信号最小阈噪比和最小信噪比;
(3)读取该次测量的噪声统计区域偏移量与噪声采样宽度,对噪声统计区域进行采样;
(4)计算得到步骤(3)中统计后的噪声平均值与均方根值;
(5)对噪声区域采样波形数据与云判区间获得的数据进行统计对比,若两者在均值、均方根值相近,且均值均较大则可认为存在云层反射,忽略该次噪声统计;否则,进入步骤(6);
(6)判断步骤(4)中的获得的噪声平均值是否超过存储器中对应地形的噪声均值初始值,噪声均方根值是否超过存储器中噪声均值初始值;若噪声平均值或均方根值不超过初始值则进入步骤(7),若噪声平均值或均方根值超过初始值则进入步骤(8);
(7)记录并保存此次测量的噪声值平均值与均方根值,将此次测量噪声结果分别与前若干次累加计算,获取的数值作为下次测量的噪声均值初始值与均方根初始值;
(8)读取设置APD偏压与调整ADC采集回波的波形识别阈值子程序;
(9)读取当前回波幅度峰值,根据系统最小信噪比与相关函数关系,计算出满足信噪比要求所需APD高压值;
(10)根据系统最小阈噪比与相关函数关系,计算得到满足阈噪比要求所需波形识别阈值。
所述步骤(1)的具体过程为:
步骤11:处理器读取该次测量恒虚警反馈控制算法函数,当前处于测量模式下,加载存储器中的初始设定参数,包括星载激光测距仪系统要求的探测虚警率、探测概率,当前所测量地形类型编号,当前地形相应的噪声均值初始值、噪声均方根初始值;
步骤12:加载存储器中该次测量的设定参数,包括噪声统计区域偏移量、噪声采样宽度、噪声采样系数,云判区间波形数据均值、均方根值,ADC采集的全波形数据回波峰值幅度、当前太阳高度角、地形反射率,APD温度值及高压调整系数、阈值调整系数。
所述步骤(2)的具体过程为:根据系统要求的虚警率与探测概率,计算信道的阈噪比最小值、信噪比最小值,计算方法为:
其中,Pfa0为虚警率,Pd0为探测概率,TNRmin为阈噪比,SNRmin为信噪比。
所述步骤(3)的具体过程为:
步骤31:读取该次测量的噪声统计区域偏移量即ADC采集噪声信号数据的起始位置,噪声统计起始位置位于回波门关闭之后,然后读取噪声采集宽度,计算出ADC采样噪声信号点数量nk,计算方法如下:
其中,WN为噪声采样宽度,c为光速,fs为ADC采样频率;
步骤32:ADC对噪声统计区域进行全波形采样并存于存储器中。
所述步骤(4)的具体过程为:
步骤41:计算噪声统计的平均值,计算方法为:
其中,yi,k是指第k次测距高速ADC采集的第i个波形幅度;
步骤42:计算噪声统计的均方根值,计算方法为:
所述步骤(5)的具体过程为:
步骤51:读取云判区间获得的数据统计值包括均值、均方根值;
步骤52:读取该次测量ADC采集的回波峰值幅度;
步骤53:将该次噪声统计的均值与均方根值分别与云判区间获得的均值、均方根进行比较,若满足如下准则,则认为存在云层反射,且不保存该次噪声统计;否则,进入步骤(6);判断云层反射存在准则如下:
|μk-μc|≤0.1μk,
|σk-σc|≤0.1σk,
μk≥0.5Ar,
其中,Ar为ADC采集的全波形数据中回波峰值幅度。
所述步骤(6)的具体过程为:
步骤61:根据当前地形类型编号读取对应地形的噪声均值初始值,噪声均值初始值;
步骤62:比较步骤(4)中的此次测量的噪声平均值与存储器中对应地形的噪声均值初始值的大小,比较噪声均方根值与存储器中噪声均值初始值的大小,若噪声平均值与均方根值均不超过初始值则进入步骤(7),若噪声平均值或均方根值超过初始值则进入步骤(8)。
所述步骤(7)的具体过程为:
步骤71:记录并保存此次测量的噪声值平均值与均方根值;
步骤72:将此次测量噪声均值、均方根值分别与前若干次累加计算,所得数值作为下次测量的噪声均值初始值与均方根初始值。
所述步骤(8)的具体过程为:
步骤81:根据系统最小信噪比、当前噪声统计均值与步骤(5)中的回波峰值幅度计算出最小噪声均方根值:
其中,SNRmin为系统最小信噪比,Ar为ADC采集的回波全波形数据中幅度值,μk为该次测量噪声统计均值,σn为最小噪声均方根值;
步骤82:读取当前测量太阳高度角,地表反射率,计算出APD所接收的背景光功率,计算方法为:
其中,已知常数包括:Es为太阳在1064nm处光谱辐照度,Δλ为窄带滤光片带宽,τa为大气透过率,ηr为接收光学效率,θr为接收视场;设置参数包括:θi为太阳高度角,ρT为对应地形的反射率;
步骤83:根据步骤81中得出的最小噪声均方根值、步骤82得出的APD接收到背景光功率与当前APD温度值的函数关系,计算出APD的倍增系数:
σn≈α0In 2
其中,M为APD倍增系数,In为APD输出噪声电流,F为APD噪声系数,已知常数包括:α0为噪声采样系数,Ids为APD的表面暗电流,Idb为APD的体暗电流,Ri为APD电流响应度,B为噪声带宽,e为电子电荷,k为玻尔兹曼常数,Rl为APD负载,keff为电子-空穴电离率;设置参数:T为APD当前温度;
步骤84:读取APD当前温度遥测值,根据步骤83计算获得的APD倍增系数,通过温度与APD偏压函数关系,计算出适应当前温度的偏压值:
VB=VB0+C(T-T0)
其中,VH为APD适应当前温度的偏压,VB为APD工作电压,已知常数包括:VB0是环境温度为T0时的工作电压,C为APD工作电压温度系数,n取值为1~3整数;
步骤85:读取高压调整系数根据步骤84得出的APD偏压值对其线性修正,计算出满足信噪比要求所需APD高压最佳值,修正原则:
VH0=α1VH+α2
其中,α1、α2是高压调整系数。
所述步骤(9)的具体过程为:
步骤91:根据系统最小阈噪比、此次测量噪声统计值,计算出满足阈噪比要求所需波形识别阈值:
因此需根据噪声统计值与阈噪比合理设置回波的预鉴别阈值。因需满足系统最小阈噪比TNRmin的要求,则阈值函数如下:
Vthmin=TNRmin*σ0
Vth0=max(αk*σk,Vthmin)+μk
其中,Vth0回波的预鉴别阈值,TNRmin为系统最小信噪比,σ0为存储器中多次统计后的噪声均方根值,若首次测量时,σ0为当前地形噪声均方根初始值,αk是阈值调整系数,为常数,σ0、μk为此次测量噪声值均方根值与平均值。
一种基于数字恒虚警的星载激光测距仪探测及反馈系统,包括APD探测器、放大及信号调理电路、高速数模转换电路、FPGA采集控制电路、模数转换电路、高压模块;
所述APD探测器用于接收机焦面回波探测,将接收到的回波光信号转换为电流信号;
所述放大及信号调理电路用于对APD输出信号进行放大,并将电平转换为模拟差分信号送入高速ADC;
所述高速数模转换电路用于对回波及噪声模拟差分信号进行模数转换,进行全波形采集;
所述FPGA采集控制电路用于对采集的噪声值统计,设置APD偏压与回波波形识别算法中预鉴别阈值;
所述模数转换电路用于根据FPGA设置的APD偏压进行数模转换,获得APD偏压控制电压;
所述高压模块用于APD偏压输入,受DAC控制输出电压大小。
本发明与现有技术相比的优点在于:采用高速ADC芯片对激光回波区间的信号进行高速采集,并对回波距离门外区域噪声多次统计,根据统计结果利用DAC动态调整APD偏压。若噪声均方根值过高,在轨调节APD高压,APD的增益因子的变化,使得噪声统计值处于合理范围内,同时,结合用高速ADC对回波脉冲进行全波形采集,调整波形预鉴别阈值,有效识别真实回波波形,滤除虚假回波。
本发明采用根据高速ADC采集的特定区域噪声统计值来调整APD的偏压,并调整回波波形预鉴别阈值,避免了卫星平台速度较大,地形变化造成回波出现奇异点造成下一次测距失效的情况;其次,提高了由于卫星平台抖动过大造成探测器输出信号的信噪比,便于测距信号提取,提高卫星平台高速运行产生太阳高度角的变化适应性;对于云层散射造成的虚假回波,通过云判区间的识别,提高了识别云层散射虚假回波的能力;实现慢虚警恒定控制,使其随环境、温度等的变化而变化;提高激光测距仪对粗糙度起伏大、反射率较大的复杂地形的适应性,满足低虚警率和高探测概率要求的同时快速精确的获取测距值。
附图说明
图1为数字恒虚警控制APD偏压电路;
图2为控制方法流程图;
图3为全波形采样时序图。
具体实施方式
下面结合附图对本发明的恒虚警探测及反馈方法作进一步详细说明。
如图1所示,从探测器输出的信号经放大后,由高速ADC数字化全波形采集,进入FPGA并运行反馈控制算法,具体步骤如下:
处理器读取该次测量恒虚警反馈控制算法函数,当前处于测量模式下,加载存储器中的初始设定参数。初始设定参数为:星载激光测距仪系统要求的探测虚警率、探测概率,当前所测量地形类型编号,当前地形相应的噪声均值初始值、噪声均方根初始值。加载存储器中该次测量的设定参数如:噪声统计区域偏移量、噪声采集宽度,噪声采样系数、云判区间波形数据均值、均方根值,ADC采集的全波形数据中回波幅度峰值、当前太阳高度角、地形反射率,APD温度值及高压调整系数、阈值调整系数。
一种基于数字恒虚警的星载激光测距仪探测及反馈方法,其特征在于步骤如下:
(1)读取并执行控制算法函数;
(2)根据星载激光测距仪系统要求的虚警率与探测概率计算地面回波信号最小阈噪比和最小信噪比;
(3)读取该次测量的噪声统计区域偏移量与噪声采样宽度,对噪声统计区域进行采样;
(4)计算得到步骤(3)中统计后的噪声平均值与均方根值;
(5)对噪声区域采样波形数据与云判区间获得的数据进行统计对比,若两者在均值、均方根值相近,且均值均较大则可认为存在云层反射,忽略该次噪声统计;否则,进入步骤(6);
(6)判断步骤(4)中的获得的噪声平均值是否超过存储器中对应地形的噪声均值初始值,噪声均方根值是否超过存储器中噪声均值初始值;若噪声平均值或均方根值不超过初始值则进入步骤(7),若噪声平均值或均方根值超过初始值则进入步骤(8);
(7)记录并保存此次测量的噪声值平均值与均方根值,将此次测量噪声结果分别与前若干次累加计算,获取的数值作为下次测量的噪声均值初始值与均方根初始值;
(8)读取设置APD偏压与调整ADC采集回波的波形识别阈值子程序;
(9)读取当前回波幅度峰值,根据系统最小信噪比与相关函数关系,计算出满足信噪比要求所需APD高压值;
(10)根据系统最小阈噪比与相关函数关系,计算得到满足阈噪比要求所需波形识别阈值。
如图2所示,所述步骤(1)的具体过程为:
步骤11:处理器读取该次测量恒虚警反馈控制算法函数,当前处于测量模式下,加载存储器中的初始设定参数,包括星载激光测距仪系统要求的探测虚警率、探测概率,当前所测量地形类型编号,当前地形相应的噪声均值初始值、噪声均方根初始值;
步骤12:加载存储器中该次测量的设定参数,包括噪声统计区域偏移量、噪声采样宽度、噪声采样系数,云判区间波形数据均值、均方根值,ADC采集的全波形数据回波峰值幅度、当前太阳高度角、地形反射率,APD温度值及高压调整系数、阈值调整系数。
所述步骤(2)的具体过程为:根据系统要求的虚警率与探测概率,计算信道的阈噪比最小值、信噪比最小值,计算方法为:
其中,Pfa0为虚警率,Pd0为探测概率,TNRmin为阈噪比,SNRmin为信噪比。所述步骤(3)的具体过程为:
步骤31:读取该次测量的噪声统计区域偏移量即ADC采集噪声信号数据的起始位置,噪声统计起始位置位于回波门关闭之后,然后读取噪声采集宽度,计算出ADC采样噪声信号点数量nk,计算方法如下:
其中,WN为噪声采样宽度,c为光速,fs为ADC采样频率;
步骤32:ADC对噪声统计区域进行全波形采样并存于存储器中。
所述步骤(4)的具体过程为:
步骤41:计算噪声统计的平均值,计算方法为:
其中,yi,k是指第k次测距高速ADC采集的第i个波形幅度;
步骤42:计算噪声统计的均方根值,计算方法为:
所述步骤(5)的具体过程为:
步骤51:读取云判区间获得的数据统计值包括均值、均方根值;
步骤52:读取该次测量ADC采集的回波峰值幅度;
步骤53:将该次噪声统计的均值与均方根值分别与云判区间获得的均值、均方根进行比较,若满足如下准则,则认为存在云层反射,且不保存该次噪声统计;否则,进入步骤(6);判断云层反射存在准则如下:
|μk-μc|≤0.1μk,
|σk-σc|≤0.1σk,
μk≥0.5Ar,
其中,Ar为ADC采集的全波形数据中回波峰值幅度。
所述步骤(6)的具体过程为:
步骤61:根据当前地形类型编号读取对应地形的噪声均值初始值,噪声均值初始值;
步骤62:比较步骤(4)中的此次测量的噪声平均值与存储器中对应地形的噪声均值初始值的大小,比较噪声均方根值与存储器中噪声均值初始值的大小,若噪声平均值与均方根值均不超过初始值则进入步骤(7),若噪声平均值或均方根值超过初始值则进入步骤(8)。
所述步骤(7)的具体过程为:
步骤71:记录并保存此次测量的噪声值平均值与均方根值;
步骤72:将此次测量噪声均值、均方根值分别与前若干次累加计算,所得数值作为下次测量的噪声均值初始值与均方根初始值。
所述步骤(8)的具体过程为:
步骤81:根据系统最小信噪比、当前噪声统计均值与步骤(5)中的回波峰值幅度计算出最小噪声均方根值:
其中,SNRmin为系统最小信噪比,Ar为ADC采集的回波全波形数据中幅度值,μk为该次测量噪声统计均值,σn为最小噪声均方根值;
步骤82:读取当前测量太阳高度角,地表反射率,计算出APD所接收的背景光功率,计算方法为:
其中,已知常数包括:Es为太阳在1064nm处光谱辐照度,Δλ为窄带滤光片带宽,τa为大气透过率,ηr为接收光学效率,θr为接收视场;设置参数包括:θi为太阳高度角,ρT为对应地形的反射率;
步骤83:根据步骤81中得出的最小噪声均方根值、步骤82得出的APD接收到背景光功率与当前APD温度值的函数关系,计算出APD的倍增系数:
σn≈α0In 2
其中,M为APD倍增系数,In为APD输出噪声电流,F为APD噪声系数,已知常数包括:α0为噪声采样系数,Ids为APD的表面暗电流,Idb为APD的体暗电流,Ri为APD电流响应度,B为噪声带宽,e为电子电荷,k为玻尔兹曼常数,Rl为APD负载,keff为电子-空穴电离率;设置参数:T为APD当前温度;
步骤84:读取APD当前温度遥测值,根据步骤83计算获得的APD倍增系数,通过温度与APD偏压函数关系,计算出适应当前温度的偏压值:
VB=VB0+C(T-T0)
其中,VH为APD适应当前温度的偏压,VB为APD工作电压,已知常数包括:VB0是环境温度为T0时的工作电压,C为APD工作电压温度系数,n取值为1~3整数;
步骤85:读取高压调整系数根据步骤84得出的APD偏压值对其线性修正,计算出满足信噪比要求所需APD高压最佳值,修正原则:
VH0=α1VH+α2
其中,α1、α2是高压调整系数。
所述步骤(9)的具体过程为:
步骤91:根据系统最小阈噪比、此次测量噪声统计值,计算出满足阈噪比要求所需波形识别阈值:
因此需根据噪声统计值与阈噪比合理设置回波的预鉴别阈值。因需满足系统最小阈噪比TNRmin的要求,则阈值函数如下:
Vthmin=TNRmin*σ0
Vth0=max(αk*σk,Vthmin)+μk
其中,Vth0回波的预鉴别阈值,TNRmin为系统最小信噪比,σ0为存储器中多次统计后的噪声均方根值,若首次测量时,σ0为当前地形噪声均方根初始值,αk是阈值调整系数,为常数,σ0、μk为此次测量噪声值均方根值与平均值。
步骤87后APD偏压值设定完成,作为下次测量的APD偏压加载值。步骤91后回波波形识别算法中预鉴别阈值的设定实现下次测量的回波波形识别。
一种基于数字恒虚警的探测及在轨反馈系统,包括APD探测器、放大及信号调理电路、高速数模转换电路、FPGA采集控制电路、模数转换电路、高压模块;
所述APD探测器用于接收机焦面回波探测,将接收到的回波光信号转换为电流信号;
所述放大及信号调理电路用于对APD输出信号进行放大,并将电平转换为模拟差分信号送入高速ADC;
所述高速数模转换电路用于对回波及噪声模拟差分信号进行模数转换,进行全波形采集;
所述FPGA采集控制电路用于对采集的噪声值统计,进行设置APD偏压与回波波形识别算法中预鉴别阈值;
所述模数转换电路用于根据FPGA设置的APD偏压进行数模转换,获得APD偏压控制电压;
所述高压模块用于APD偏压输入,受DAC控制输出电压大小。
本发明说明书中未作详细描述的内容属本领于技术人员的公知技术。
Claims (7)
1.一种基于数字恒虚警的星载激光测距仪探测及反馈方法,其特征在于步骤如下:
(1)读取并执行控制算法函数;
(2)根据星载激光测距仪系统要求的虚警率与探测概率计算地面回波信号最小阈噪比和最小信噪比;
(3)读取该次测量的噪声统计区域偏移量与噪声采样宽度,对噪声统计区域进行采样;
(4)计算得到步骤(3)中统计后的噪声平均值与均方根值;
(5)对噪声区域采样波形数据与云判区间获得的数据进行统计对比,若两者在均值、均方根值相近,且均值均较大则可认为存在云层反射,忽略该次噪声统计;否则,进入步骤(6);
(6)判断步骤(4)中的获得的噪声平均值是否超过存储器中对应地形的噪声均值初始值,噪声均方根值是否超过存储器中噪声均值初始值;若噪声平均值或均方根值不超过初始值则进入步骤(7),若噪声平均值或均方根值超过初始值则进入步骤(8),具体过程为:
步骤61:根据当前地形类型编号读取对应地形的噪声均值初始值;
步骤62:比较步骤(4)中的此次测量的噪声平均值与存储器中对应地形的噪声均值初始值的大小,比较噪声均方根值与存储器中噪声均值初始值的大小,若噪声平均值与均方根值均不超过初始值则进入步骤(7),若噪声平均值或均方根值超过初始值则进入步骤(8);
(7)记录并保存此次测量的噪声值平均值与均方根值,将此次测量噪声结果分别与前若干次累加计算,获取的数值作为下次测量的噪声均值初始值与均方根初始值;
(8)读取设置APD偏压与调整ADC采集回波的波形识别阈值子程序;
(9)读取当前回波幅度峰值,根据系统最小信噪比与相关函数关系,计算出满足信噪比要求所需APD高压值;具体过程为:
具体过程为:
步骤91:根据系统最小信噪比、当前噪声统计均值与步骤(5)中的回波峰值幅度计算出最小噪声均方根值:
其中,SNRmin为系统最小信噪比,μk为该次测量噪声统计均值,σn为最小噪声均方根值;
步骤92:读取当前测量太阳高度角,地表反射率,计算出APD所接收的背景光功率,计算方法为:
其中,已知常数包括:Es为太阳在1064nm处光谱辐照度,Δλ为窄带滤光片带宽,τa为大气透过率,ηr为接收光学效率,θr为接收视场;设置参数包括:θi为太阳高度角,ρT为对应地形的反射率;
步骤93:根据步骤91中得出的最小噪声均方根值、步骤92得出的APD接收到背景光功率与当前APD温度值的函数关系,计算出APD的倍增系数:
σn≈α0In 2
其中,M为APD倍增系数,In为APD输出噪声电流,F为APD噪声系数,已知常数包括:α0为噪声采样系数,Ids为APD的表面暗电流,Idb为APD的体暗电流,Ri为APD电流响应度,B为噪声带宽,e为电子电荷,k为玻尔兹曼常数,Rl为APD负载,keff为电子-空穴电离率;设置参数:T为APD当前温度;
步骤94:读取APD当前温度遥测值,根据步骤93计算获得的APD倍增系数,通过温度与APD偏压函数关系,计算出适应当前温度的偏压值:
VB=VB0+C(T-T0)
其中,VH为APD适应当前温度的偏压,VB为APD工作电压,已知常数包括:VB0是环境温度为T0时的工作电压,C为APD工作电压温度系数,n取值为1~3整数;
步骤95:读取高压调整系数根据步骤94得出的APD偏压值对其线性修正,计算出满足信噪比要求所需APD高压最佳值,修正原则:
VH0=α1VH+α2
其中,α1、α2是高压调整系数;
步骤96:根据系统最小阈噪比、此次测量噪声统计值,计算出满足阈噪比要求所需波形识别阈值,具体为:
根据噪声统计值与阈噪比合理设置回波的预鉴别阈值;因需满足系统最小阈噪比TNRmin的要求,则阈值函数如下:
Vthmin=TNRmin*σ0
Vth0=max(αk*σk,Vthmin)+μk
其中,Vth0回波的预鉴别阈值,TNRmin为系统最小信噪比,σ0为存储器中多次统计后的噪声均方根值,若首次测量时,σ0为当前地形噪声均方根初始值,αk是阈值调整系数,为常数,σ0、μk为此次测量噪声值均方根值与平均值;
(10)根据系统最小阈噪比与相关函数关系,计算得到满足阈噪比要求所需波形识别阈值。
2.根据权利要求1所述的一种基于数字恒虚警的星载激光测距仪探测及反馈方法,其特征在于:所述步骤(1)的具体过程为:
步骤11:处理器读取该次测量恒虚警反馈控制算法函数,当前处于测量模式下,加载存储器中的初始设定参数,包括星载激光测距仪系统要求的探测虚警率、探测概率,当前所测量地形类型编号,当前地形相应的噪声均值初始值、噪声均方根初始值;
步骤12:加载存储器中该次测量的设定参数,包括噪声统计区域偏移量、噪声采样宽度、噪声采样系数,云判区间波形数据均值、均方根值,ADC采集的全波形数据回波峰值幅度、当前太阳高度角、地形反射率,APD温度值及高压调整系数、阈值调整系数。
6.根据权利要求5所述的一种基于数字恒虚警的星载激光测距仪探测及反馈方法,其特征在于:所述步骤(5)的具体过程为:
步骤51:读取云判区间获得的数据统计值包括均值、均方根值;
步骤52:读取该次测量ADC采集的回波峰值幅度;
步骤53:将该次噪声统计的均值与均方根值分别与云判区间获得的均值、均方根进行比较,若满足如下准则,则认为存在云层反射,且不保存该次噪声统计;否则,进入步骤(6);判断云层反射存在准则如下:
|μk-μc|≤0.1μk,
|σk-σc|≤0.1σk,
μk≥0.5Ar,
其中,Ar为ADC采集的全波形数据中回波峰值幅度。
7.根据权利要求1-6任一所述的一种基于数字恒虚警的星载激光测距仪探测及反馈方法,其特征在于:所述步骤(7)的具体过程为:
步骤71:记录并保存此次测量的噪声值平均值与均方根值;
步骤72:将此次测量噪声均值、均方根值分别与前若干次累加计算,所得数值作为下次测量的噪声均值初始值与均方根初始值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810072848.3A CN108333592B (zh) | 2018-01-25 | 2018-01-25 | 基于数字恒虚警的星载激光测距仪探测及反馈方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810072848.3A CN108333592B (zh) | 2018-01-25 | 2018-01-25 | 基于数字恒虚警的星载激光测距仪探测及反馈方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108333592A CN108333592A (zh) | 2018-07-27 |
CN108333592B true CN108333592B (zh) | 2020-05-12 |
Family
ID=62926553
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810072848.3A Active CN108333592B (zh) | 2018-01-25 | 2018-01-25 | 基于数字恒虚警的星载激光测距仪探测及反馈方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108333592B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108919870A (zh) * | 2018-09-19 | 2018-11-30 | 明峰医疗系统股份有限公司 | 一种通过常规芯片器件实时调节spm负高压偏压的方法 |
US11041944B2 (en) * | 2019-03-01 | 2021-06-22 | Beijing Voyager Technology Co., Ltd. | Constant false alarm rate detection in pulsed LiDAR systems |
CN113614566A (zh) * | 2019-03-27 | 2021-11-05 | 松下知识产权经营株式会社 | 测距方法、测距装置以及程序 |
CN110954913A (zh) * | 2019-11-22 | 2020-04-03 | 同济大学 | 一种星载激光测高地面数据处理系统 |
CN111830488B (zh) * | 2020-06-28 | 2023-06-27 | 中山大学 | 一种基于gm-apd激光测距系统的回波信号数据处理方法及系统 |
CN112630746B (zh) * | 2020-12-01 | 2023-09-15 | 北京遥感设备研究所 | 一种用于远距目标测量的脉冲多普勒激光雷达 |
CN114047522B (zh) * | 2022-01-12 | 2022-03-22 | 中大检测(湖南)股份有限公司 | 一种智能激光测距仪控制系统 |
CN114578384B (zh) * | 2022-05-07 | 2022-09-09 | 成都凯天电子股份有限公司 | 一种用于激光大气系统的自适应恒虚警检测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104252194A (zh) * | 2014-08-29 | 2014-12-31 | 北京航天控制仪器研究所 | 一种apd偏压自动调节装置和方法 |
CN106526570A (zh) * | 2016-10-19 | 2017-03-22 | 北京空间机电研究所 | 星载激光测距仪回波探测电路 |
CN106569249A (zh) * | 2016-10-14 | 2017-04-19 | 北京空间机电研究所 | 一种星载Si‑APD探测器反向偏压自动调节方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006041628A (ja) * | 2004-07-22 | 2006-02-09 | Sumitomo Electric Ind Ltd | 光受信回路 |
-
2018
- 2018-01-25 CN CN201810072848.3A patent/CN108333592B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104252194A (zh) * | 2014-08-29 | 2014-12-31 | 北京航天控制仪器研究所 | 一种apd偏压自动调节装置和方法 |
CN106569249A (zh) * | 2016-10-14 | 2017-04-19 | 北京空间机电研究所 | 一种星载Si‑APD探测器反向偏压自动调节方法 |
CN106526570A (zh) * | 2016-10-19 | 2017-03-22 | 北京空间机电研究所 | 星载激光测距仪回波探测电路 |
Non-Patent Citations (4)
Title |
---|
一种应用于高精度脉冲激光测距的自动增益控制方法;郭颖 等;《激光与红外》;20100415;345-348 * |
基于APD C30927E的激光回波探测电路噪声抑制技术研究;刘志豪;《中国优秀硕士学位论文全文数据库 信息科技辑》;20140115;I136-1071 * |
基于FPGA的激光雷达恒虚警率控制技术研究;欧阳俊华;《红外与毫米波学报》;20090215;50-53 * |
星载激光测距仪APD最佳雪崩增益控制技术研究;李旭 等;《红外与激光工程》;20160515;1-6 * |
Also Published As
Publication number | Publication date |
---|---|
CN108333592A (zh) | 2018-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108333592B (zh) | 基于数字恒虚警的星载激光测距仪探测及反馈方法和系统 | |
Turman | Analysis of lightning data from the DMSP satellite | |
Christian et al. | The detection of lightning from geostationary orbit | |
McManamon et al. | Comparison of flash lidar detector options | |
EP2446301B1 (en) | Pulsed light optical rangefinder | |
CN112859098B (zh) | 一种光子数分辨测量增强单光子激光雷达系统及测距方法 | |
CN109375194A (zh) | 用于激光雷达的模拟前端读出电路 | |
CN108919282A (zh) | 一种激光雷达信号时刻鉴别系统 | |
WO2022160610A1 (zh) | 一种飞行时间测距方法、系统和设备 | |
KR102317072B1 (ko) | 라이다 시스템에서의 시간-디지털 변환 방법 및 장치 | |
US10690448B2 (en) | Method and apparatus for variable time pulse sampling | |
CN110501716B (zh) | 基于单光子激光雷达背景噪声率的地表分类方法 | |
JP2014081254A (ja) | 光学的測距装置 | |
CN113196105A (zh) | 数字像素 | |
WO2022161481A1 (zh) | 一种飞行时间测距方法、系统和设备 | |
Zhang et al. | Evaluation of the performance characteristics of the lightning imaging sensor | |
Suszcynsky et al. | Coordinated observations of optical lightning from space using the FORTE photodiode detector and CCD imager | |
Kirkland et al. | Optical observations of terrestrial lightning by the FORTE satellite photodiode detector | |
Jacobson et al. | Revisiting" Narrow Bipolar Event" intracloud lightning using the FORTE satellite | |
Jacobson et al. | Joint radio and optical observations of the most radio-powerful intracloud lightning discharges | |
US20220373690A1 (en) | Atmospheric sensor using programmable time-gated detection aperture | |
JP2820095B2 (ja) | レーザ計測装置 | |
Yan et al. | Pulse-based machine learning: Adaptive waveform centroid discrimination for LIDAR system | |
Kaifler et al. | The polar mesospheric cloud dataset of the Balloon Lidar Experiment BOLIDE | |
Huntington et al. | 512-element linear InGaAs APD array sensor for scanned time-of-flight lidar at 1550 nm |
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 |