1 武汉大学 电子信息学院,湖北 武汉 430072
2 中国空间技术研究院,北京 100098
3 武汉量子技术研究院,湖北 武汉 430010
硅光电倍增管(Silicon Photomultiplier,SiPM)具有极高的探测灵敏度和响应速度,且在多光子条件下具有较高的动态范围以及线性响应的特性,在光子计数激光雷达应用中有独特的优势。然而,由于SiPM多像元、单时间通道的工作模式,其输出电压信号相较于其他单光子探测器更大概率出现脉冲堆叠现象,不同鉴别阈值条件下的SiPM探测过程更为复杂。针对该问题,本文建立了SiPM光子事件响应模型,以此为基础分析了由脉冲堆叠引发的屏蔽效应和触发效应两种特殊情况的时域分布,最终建立了SiPM半解析的探测概率与虚警概率模型。同时,搭建了基于SiPM探测器的光子计数雷达系统,通过观察实测输出电压波形以及光子点云分布与理论模型相符(R2>0.95)。通过查全率与查准率对不同鉴别阈值的SiPM光子点云分布进行定量评价,并给出最优鉴别阈值区间,这对基于SiPM的光子计数激光雷达系统硬件参数的优化设计以及探测性能的定量分析具有重要的指导意义。
硅光电倍增管 脉冲堆叠 鉴别阈值 探测概率模型 photomultiplier tubes pulse pile up discrimination threshold detection probability model
美国NASA于2018年发射的ICESat-2 (The Ice, Cloud, and land Elevation Satellite-2)卫星上搭载的ATLAS (Advanced Topographic Laser Altimeter System)是目前为止全球唯一一个对地观测的星载光子计数激光雷达,具有较高的轨向空间采样率,为用遥感的方法探测海浪要素提供了可能。光子计数激光雷达用于海浪探测的前提是能够准确地提取来自海面的信号光子,并确定瞬时的海面廓线。迄今为止,用星载光子计数激光雷达探测海面形态和海浪要素的研究鲜见报道,也缺少专门针对海面信号光子的提取方法。基于海面信号光子的分布特点,文中提出了一种新的信号提取算法:首先通过直方图统计及自适应的阈值选取完成对海面回波光子的粗去噪;然后基于激光雷达光斑尺寸和海面波动特点,选取合适的搜索邻域计算信号点和噪声点密度,根据两者点密度差异对信号光子和噪声光子分类;最后用高斯函数拟合的方法进一步去除密度较大的后向散射噪声光子,最终得到来自海面反射的信号光子。利用上述算法提取了太平洋7个不同海况区域的海面信号光子和瞬时海面廓线并进一步计算出当地海浪的峰值波长和周期。将计算结果与同期欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasts, ECMWF)的全球大气再分析ERA5(ECMWF Re-Analysis5)数据作对比,在不同风速、水深的海域都获得了基本一致的结果,超过半数区域的海浪周期误差在5%以内,初步证明了星载光子计数激光雷达观测成果用于海浪要素计算的可行性。
激光雷达遥感 ICESat-2/ATLAS 海面信号去噪 海浪要素计算 lidar remote sensing ICESat-2/ATLAS sea surface signal denoising calculation of ocean wave elements 红外与激光工程
2023, 52(2): 20220366
1 武汉大学 电子信息学院,湖北 武汉 430072
2 中国空间技术研究院,北京 100098
3 武汉量子技术研究院,湖北 武汉 430010
光电倍增管(Photomultiplier Tubes,PMT)具有光子级别的灵敏度,而且没有物理意义上的死区时间,在光子计数雷达应用中有独有的优势,但PMT响应单光子输出脉冲高度呈高斯随机分布,且不同脉冲之间还有可能产生堆叠,传统的单光子模型不能准确描述PMT的光子探测过程。针对PMT输出脉冲高度随机分布,脉冲堆叠现象和光子事件鉴别阈值的设置对光子事件探测概率的影响提出了一个新的PMT光子探测模型,并根据实际应用场景进行了简化,通过蒙特卡洛仿真验证了简化模型的适用性。搭建了光子计数雷达系统,通过和盖革模式APD进行实验对比,发现PMT光子探测模型虽然有略微的探测概率损失,但在测距应用中具有更小行走误差和更高测距精确度。实验也证明了新模型相比于传统单光子探测模型更符合PMT的光子事件探测概率分布,新模型对基于PMT光子计数雷达的系统设计和理论分析具有重要的指导意义。
光电倍增管 脉冲堆叠 脉冲高度分布 探测概率模型 photomultiplier tubes pulse pile up pulse height distribution detection probability model
红外与毫米波学报
2022, 41(6): 1051
红外与激光工程
2022, 51(10): 20220003
强激光与粒子束
2021, 33(7): 074001
红外与激光工程
2020, 49(11): 20200214
单光子激光雷达所获取的光子事件存在随机分布的特点,使得其激光测距值出现不确定性,从而降低了单光子激光雷达的测距精度。在不减小采样分辨率的情况下,采用累积邻近多光斑的光子事件来构建光子累积直方图,并基于反算的目标响应函数的时间重心来确定激光测距值。针对ICESat-2单光子激光雷达,以测距均方根误差(RMSE)和平均绝对误差(MAE)为性能指标,构建出一种综合考虑陆地地形、累积光斑数目和回波光子数等多因素对其测距误差影响的评估方法。同时,选取ICESat-2过境美国犹他州西瓦利城的某观测条带的随机地形数据进行验证分析。结果表明,所提出的激光测距值解算方法能够使该条带的测距RMSE值由114.25 cm降低到63.84 cm,MAE值由70.97 cm降低到48.52 cm,均优于ICESat-2数据产品提供的137.96 cm RMSE值和97.24 cm MAE值,这对提升单光子激光雷达在陆地区域的测距精度具有一定的借鉴作用。
单光子激光雷达 累积直方图 目标响应函数 测距误差 single-photon lidar cumulative histogram target response function ranging error 红外与激光工程
2020, 49(11): 20200247
星载光子计数激光雷达作为一种新的探测体制激光雷达,已开始应用于海面测量。然而受海风等多种因素的影响,海面存在一定的粗糙度和较大的起伏变化,因此光子计数激光雷达返回的信号点云在返回能量和信号光子分布上存在较大的变化,潜在的影响到了海面高程测量精度。本文基于JONSWAP海浪谱和微面元模型理论,结合蒙特卡洛方法建立了光子计数激光雷达海洋目标的仿真模型。以ICESat-2星载光子计数激光雷达的系统参数作为输入,仿真了不同风速条件下海面的信号光子分布,通过与ICESat-2实测结果对比证明了仿真方法的正确性。基于仿真模型,分析了不同风速条件下,光子计数激光雷达的测距误差分布。结果表明,光子计数激光雷达测得的海面高程小于实际参考海面,且测量偏差和标准差随风速增加而增大,当风速为10m/s,累计脉冲次数为100次时,测量偏差约为-2.5cm,标准差为3.6cm。所建立的仿真模型和分析结果对优化针对海面观测的星载光子计数激光雷达的系统参数设计和平均海面观测结果修正具有重要的参考意义。
光子计数激光雷达 海洋仿真 微面元理论 误差分析 photon-counting sea surface simulation micro plane element theory error analysis
激光雷达脉冲回波的波形分解方法是提取其波形参数的重要手段,也为反演目标高度、倾斜度和粗糙度、反射率提供直接的参数来源。针对部分信噪比较差且具有一定混叠程度的脉冲回波,提出一种基于可变分量的参数随机抽样方法的波形分解算法(WDVCM)。该算法以高斯混合函数为优化模型,通过随机产生高斯分量的特征参数以及删减或生成高斯分量等操作,并分别基于能量函数和拟合标准差作为参数优化的判据,从而实现波形的分解及其参数提取。利用该算法对美国国家航空航天局(NASA)的对地观测星载激光雷达(GLAS)一个条带中的4584个原始波形进行了处理分析。结果发现,约99%的WDVCM和97%的NASA拟合波形结果的相关系数均超过0.95,其中两者相关系数差异不超过0.05占98%。同时,WDVCM和NASA拟合波形的标准差系数均值分别为2.21和3.28,约89%的WDVCM拟合波形的标准差系数均小于NASA拟合波形的标准差系数。所得结果表明,WDVCM对混叠高斯波形的拟合效果更好,适用性更强。
激光雷达 波形分解 特征参数 相关系数 拟合标准差 Lidar waveform decomposition characteristic parameter correlation coefficient standard deviation of fitting 红外与激光工程
2019, 48(10): 1005009