基于微通道板的单光子激光测高技术研究 下载: 1055次封面文章
1 引言
激光测高技术是一种先进的测量手段,在激光器出现之初,其作为传统雷达技术与激光技术的结合产物就已被提出来[1]。卫星激光测高技术以卫星为平台,搭载激光器、探测器和计时器,通过测量激光从激光器发射到地面反射之后再被探测器接收的时间,就可以知道激光传输路径的精确距离。再结合精确的卫星轨道和姿态信息,就可以获取地面三维模型。激光测高技术可以被应用于气候研究、调查森林覆盖率、城市规划等。1994年,美国第一次在探月卫星Clementine号上搭载激光测高仪[2-4]来获取高精度的月球表面特征信息。目前,激光测高技术已得到快速发展和广泛应用[5]。
美国的Clementine、ICESat-1、LOLA等[3-4,6-7],以及我国的嫦娥一号、嫦娥二号、资源三号02星[8-10]都搭载有激光测高仪,日本、印度等也有成功发射星载激光测高仪的记录[11-13],这表明激光测高仪在遥感观测中具有重要意义。目前发射的激光测高仪都采用工作在线性模式的雪崩光电二极管(APD),可探测的峰值功率在10-9~10-8 W量级,需要激光器的单脉冲能量在10-2~10-1 J量级[14],大脉冲能量制约了激光器脉冲的重复频率,同时也会影响激光器的寿命。激光测高仪的发展趋势是采用单光子探测器,单光子探测器的灵敏度可达10-15~10-14 W量级[14],可以提高脉冲的重复频率,减小足迹间隔,获得更高的空间分辨率。目前,基于卫星平台的单光子激光测高技术还在研究阶段,NASA计划中的ICESat-2和LIST卫星都将搭载单光子激光测高仪[15-16]。
单光子激光测高仪是激光测高仪的发展趋势,单光子激光测高技术与以前的技术有很大区别,点云数据的处理是一个难题,目前国内对这方面的研究还很少。本文对单光子激光测高技术进行了数学建模,对误差进行了理论分析和数值计算,对单光子器件死时间带来的系统误差进行了理论分析和数值计算,并建立了地物模型进行仿真计算,提出一种优化高程信息的方法。
2 单光子激光测高技术
2.1 单光子激光测高原理
星载激光测高仪类似于激光测距仪,用飞行时间法(TOF)对地面成像,可以获得地面目标的高程信息[17]。
单光子激光测高的原理是用窄脉冲、高重复频率的激光照射地面,然后用探测器读出返回的光子,每个光子可以精确定位三维坐标。
2.2 探测器
单光子探测器主要有两种:微通道板(MCP)和APD[16,18-19]。APD根据其淬火电路的不同又可分为被动淬火APD(PQ-APD)和主动淬火APD(AQ-APD)。
APD是一种高性能的光电传感器,普遍应用在各种场景中。当APD偏置电压超过其雪崩击穿电压时,称为盖革模式APD,可以响应单光子[20]。APD由于特殊的工作原理,在雪崩发生之后不能自行恢复到探测状态,需要配合淬火电路来重置[21]。重置过程中,探测器不能响应到来的光子,导致其死时间较长,由于每个APD单元都有独立的淬火电路,导致APD阵列的集成度不够,APD阵列像素尺寸在0.5~2 mm之间,像素间隙为0.3~0.4 mm,典型的分辨率有4 pixel×4 pixel、8 pixel×8 pixel、16 pixel×16 pixel,像素总数只有几百[22-23]。集成度不够高导致APD阵列体积大且很重,极大地增加了测高系统的设计难度。
MCP是一种二维平面型电子倍增器,由大量平行排布的单通道电子倍增器组成。通道内壁有较高的倍增系数,光子从光阴极入射转化为电子,电子在电场作用下连续碰撞内壁实现多次倍增,从而增强光子图像[24]。通道的孔径通常为6~12 μm,器件像素总数可以达到百万量级[25]。
探测器的死时间会影响计数率。在强烈的阳光背景下,背景光子可能会触发探测器,导致探测器进入死时间,无法被信号光子触发,所以信号光子的计数率会下降,用信号光子计数率下降因子
式中:
从
与APD相比,MCP有一个缺点,即提取过量的电荷会降低其寿命。此外,MCP的电流会饱和,在强光刺激下,MCP输出电流接近饱和时,增益会降低,这种饱和机制能在强光环境下保护探测器。
2.3 信号预滤波
对于单光子激光测高仪,有三种预滤波技术可以使用[26],它们分别是光谱滤波、空间滤波和时间滤波。光谱滤波是指选用与激光器波长匹配的窄带滤光片,在保证信号光子透过率的情况下,尽量缩小带宽,从而有效滤除背景噪声。空间滤波是指根据激光器投射在地面的光斑尺寸来设计探测器的光学系统,使视场与光斑匹配,获取尽可能多的信号光子。时间滤波是指探测器受选通门控制,只关心地面高程范围内返回的光子,而不响应云层、大气散射的光子。
空间滤波能增大信号的能量,增加探测到的事件数;时间滤波能有效避免云层和大气散射的干扰;空间滤波和时间滤波能提高信噪比,但不能提高探测器的
2.4 探测器的选择
与APD相比,MCP有诸多优势:MCP是二维元件,很容易达到百万量级的像素数(APD想要获取二维信息就必须拼接成APD阵列,其像素数也仅有几百量级);MCP具有高增益、增益可变、电流自饱和、防强光等优点;MCP搭配交叉延迟线(CDL)可以获取光子的(
3 激光测距的性能
3.1 随机误差
3.1.1 随机误差模型
若要获得三维的地面高程信息,就需要建立地面坐标和高程信息之间的映射。垂直于测距方向的误差,即被测距点在水平坐标上的定位误差,会影响测距精度。在测距过程中,对卫星精确定位、定轨以及姿态控制是星载激光测高仪精确工作的基本要求。将一组特定的激光测高信息考虑为成千上万激光测高点云中的一个点,则每一组测高信息应该包括水平方向上的坐标(
式中:
取光斑的中心为理想(
式中:
大气湍流引起的横向误差为
式中:
卫星指向不准确引起的水平方向误差可以表示为
式中:
沿测距方向的误差
式中:
飞行时间误差引起的测距误差可以表示为
式中:
地表斜面引起的测距误差可以表示为
式中:tan
随机误差的模型为加性高斯噪声,概率密度分布符合均值为0,半峰全宽为
3.1.2 随机误差数值计算
为了计算星载激光测高仪的随机误差,以ICESat-2上搭载的激光测高仪ATLAS参数作为基础[27]设置仿真参数。ATLAS的性能参数如
表 1. ATLAS的性能参数
Table 1. Performance parameters of ATLAS
|
光斑直径即(2)和(3)式中的
在(8)式中,
3.2 回波信号分析
3.2.1 回波光子数
太阳光的成分比较复杂,是激光探测的主要干扰光源,太阳光的一部分通过大气散射进入探测器的镜头,一部分被地表反射进入镜头,通过探测器的窄带滤波器后被单光子器件响应。由于太阳发光在时间上是均匀的,所以这部分光子均匀地散布在整个空间,不带有特定的时间特征,这部分噪声就是背景噪声。
对于单光子计数,入瞳处接收到的背景噪声光子数可以表示为
式中:
激光器发射脉冲后经地面反射回到望远镜入瞳处的信号光子数可以表示为
式中:
激光的带宽可以表示为
式中:Δ
3.2.2 信噪比的数值计算
对于单光子探测的激光脉冲,其单脉冲的能量低于传统线性探测巨脉冲的能量,但脉冲的重复频率高,ATLAS脉冲的重复频率为10 kHz[27],两个相邻脉冲的时间间隔为0.1 ms。为了防止混淆,用不同的时间窗口来区分不同的出射脉冲,ICESat-2的轨道高度
ATLAS激光器的中心波长
使用ModTran来计算探测器接收到的背景光强度,正午太阳直射地表时对探测器的影响最大,本研究比较关心太阳光中532 nm的分量,所以选用6种不同的大气模型。在太阳直射地表和探测器天顶角为180°的条件下,计算经大气散射和地表反射后探测器观察到的地表方向的辐亮度,以及大气的透过率,结果如
保守估计下,选择最坏的参数,滤波器Δ
表 2. 典型大气模型的影响
Table 2. Influence of typical atmosphere model
|
将同样的参数代入(11)式可得:当地面反射率为1时,探测器接收到的信号的单脉冲能量为9.641×10-15 J,光子数为25803.1。
不同地表反射率下信噪比的变化如
由
3.2.3 回波光子分布
在激光测距中,时间
式中:
假设被测地面的高度为100 m,测高误差为
背景噪声的光子总数为729.2,背景噪声分布如
3.3 探测概率分析
3.3.1 探测概率模型
MCP的死时间约为2 ns,CDL也会引入纳秒级别的死时间,根据以上的计算可知回波脉冲宽度为2.6 ns,MCP/CDL的系统死时间大于回波脉冲宽度,整个回波脉冲最多只能触发一次计数,所以设定不合理的探测率可能会引入系统误差。
每个独立的光子触发探测器是独立事件,探测器探测到的光子的概率为
式中:
探测器被光子触发的概率分布满足几何分布,可以表示为
式中:
式中:
3.3.2 探测概率的数值计算
在整个回波脉冲中,若干个光子进入望远镜,探测率过低会导致探测的光子数不足,信息量不足,很难重建地面的高程模型。但探测率也不是越高越好,若探测率过高,探测器容易被回波脉冲的前沿触发,引入系统误差,甚至可能被噪声光子触发,导致探测器一直处于死时间,回波脉冲就不能以正确的时间触发探测器。
考虑探测器第一次被触发的时间以及引入的误差的分布,以几何分布的模型为基础,计算得到测距误差的概率分布,结果如
图 5. 不同探测率下的误差概率分布
Fig. 5. Probability distribution of errors at different detection rates
不同探测率下峰值探测率和半峰全宽如
表 3. 不同探测率下的峰值探测率和半峰全宽
Table 3. Peak detection rate and FWHM at different detection rates
|
测距误差的分布不满足高斯分布,但呈类似于高斯分布的形状,过高的探测率会导致探测器响应回波脉冲的前沿,使测量得到的飞行时间偏小,测量的高程偏高。在
由于探测器响应回波脉冲的前沿,所以会导致测量距离偏小,其带来的误差是正值,而且随着探测率增加,误差增大。在一般的激光测距系统中,单光子探测器的探测率选用80%就能兼顾误差和测距效率,如
随着探测率增加,峰值高度增高,半峰全宽降低。峰值半峰全宽积定义为峰值高度与半峰全宽的乘积,是一个无量纲的量,从
4 仿真
4.1 提取高程信息
为了验证系统的性能,以及验证点云信息处理算法,对激光测高系统进行建模仿真。首先建立地物模型,地物模型如
图 7. 不同参数与探测率的关系。(a)半峰全宽与探测率的关系;(b)峰值探测率与探测率的关系;(c)峰值半峰全宽积与探测率的关系
Fig. 7. Relationship between different parameters and detection rate. (a) Relationship between FWHM and detection rate; (b) relationship between peak detection rate and detection rate; (c) relationship between detection rate and the product of peak height and FWHM
单地元尺寸为10 m,整个地物模型包含64×64个地元,分为地面背景和前景,地面背景高度为0,前景有4个部分:B1是横截面尺寸为300 m×400 m的方盒,高度为1 m;B2是横截面尺寸为150 m×150 m的方盒,高度为1.75 m;B3是横截面尺寸为300 m×20 m的窄条形方盒,高度为1.4 m;C1是直径为80 m的圆柱体,高度为1 m。
使用蒙特卡罗模型进行仿真,光斑直径为10 m,则匹配的地元尺寸为10 m。根据600 km的近地轨道以及10 kHz的重复频率,可以计算得到两个光斑之间的步长为0.688 m。对每个像元都进行14次脉冲发射,因为单脉冲触发探测器的概率为80%,所以平均有11个事件被记录。
点云图是一系列带有(
显然,误差越大,从点云中提取高程信息就会越困难。为了滤除背景噪声,以及抵消横向误差
图 9. 事件点云图。(a) Δz=0.1 m; (b)Δz=0.2 m; (c) Δz=0.3 m; (d) Δz=0.39 m
Fig. 9. Point clouds of event. (a) Δz=0.1 m; (b) Δz=0.2 m; (c) Δz=0.3 m; (d) Δz=0.39 m
滤波之后基本能还原地物模型的形状,但因为误差和噪声的影响,平台处不平整。高程信息的直方图如
从
评价测高精度一般用均方根误差,其定义为
式中:
4.2 优化的高程信息
可以注意到,激光测高与图像分割有相似的形式,高程信息的突变处往往伴随光影的明显变化,比如森林的边界、楼房的阴影、海洋沿岸。对于大部分卫星,激光测高仪是卫星上众多载荷中的一个,一般具有遥感观测能力的卫星,都配有遥感相机。典型的遥感相机有R、G、B三个通道,或加一个近红外NI通道共4个通道。可以使用多光谱相机所成的像来提取图像的边界,并对带有误差的高程信息进行优化。优化的目标函数表示为[29]
式中:
式中:
使目标函数极小的
式中:
把地物模型转化为RGB三通道的遥感图像,并给地面、B1、B2、B3、C1涂上不同的颜色,如
取
优化后的高程直方图如
这样的修正可以消除大部分噪声,并可以还原地物模型的真实形状。优化后,高程的均方根误差为2.6 cm,说明用遥感图像优化高程信息的方法是有效的。
5 结论
目前,基于单光子探测器的激光测高仪尚未应用在卫星平台上,理论研究也不完善,本研究对单光子激光测高技术进行梳理,从激光的数学模型出发,到大气模型、地物模型,对激光测高的各个环节都进行了讨论。代入典型值进行数值计算,得到激光回波的半峰全宽为0.39 m,标准差为0.17 m。为研究激光测高点云的处理方法,搭建仿真地物模型,用蒙特卡罗方法进行仿真计算,采用中值滤波处理点云,计算得到高程的均方根误差为6.1 cm,低于激光的标准差,这是因为高程信息是由多个光子事件联合计算得到的,相当于多次采样取平均值。对于单光子激光测高仪,如果能提高激光脉冲的重复频率,就可以进一步提高测高精度。最后提出一种利用遥感的多光谱图像来优化高程信息的算法,先验知识表明高程信息变化的地方往往是遥感图像RGB通道发生变化的地方,所以用遥感图像的边界来约束高程信息进行优化。仿真结果表明优化后的高程的均方根误差为2.6 cm,优于中值滤波的结果。
[1] 李然, 王成, 苏国中, 等. 星载激光雷达的发展与应用[J]. 科技导报, 2007, 25(14): 58-63.
[4] 唐新明, 李国元. 激光测高卫星的发展与展望[J]. 国际太空, 2017( 11): 13- 18.
Tang XM, Li GY. Development and prospect of laser altimetry satellite[J]. Space International, 2017( 11): 13- 18.
[5] 杨馥, 贺岩, 周田华, 等. 基于伪随机码调制和单光子计数的星载测高计仿真[J]. 光学学报, 2009, 29(1): 21-26.
[6] Schutz BE. Laser altimetry and lidar from ICESat/GLAS[C]. 2001 International Geoscience and Remote Sensing Symposium, Sydney, 2001: 1016- 1019.
[8] 王建宇, 舒嵘, 陈卫标, 等. 嫦娥一号卫星载激光高度计[J]. 中国科学: 物理学力学天文学, 2010, 40(8): 1063-1070.
[9] . 嫦娥一号卫星的初步科学成果与嫦娥二号卫星的使命[J]. 航天器工程, 2010, 19(5): 1-6.
[10] 宋博, 李旭, 郑伟, 等. 02)星载高精度激光测距技术的实现[J]. 光电子技术, 2017, 37(1): 61-65.
[14] 季云飞, 耿林, 冯国旭, 等. 激光测高技术的发展趋势[J]. 激光与红外, 2011, 41(8): 830-833.
[16] Yu AW, Harding DJ, KrainakM, et al. Development of an airborne lidar surface topography simulator[C]. 2011 Laser Applications to Photonic Applications, CLEO: Applications and Technology, Baltimore, 2011: 1- 2.
[17] 蔡厚智, 刘进元, 付文勇, 等. 基于微通道板选通的飞行时间测量技术[J]. 光学学报, 2018, 38(2): 0204002.
[18] Priedhorsky WC, Smith RC, ChengH. Laser ranging and mapping with a photon-counting detector[C].Proceedings of SPIE, 1995: 441- 52.
[19] Baron MH, Priedhorsky WC. Crossed-delay line detector for ground- and space-based applications[C]. Proceedings of SPIE, 1993: 188- 198.
[20] Aull BF, Marino RM. Three-dimensional imaging with arrays of Geiger-mode avalanche photodiodes[C]. Proceedings of SPIE, 2005, 6014: 60140D.
[21] 张国青, 张英堂, 翟学军, 等. 多像素光子计数器的信噪比特性[J]. 光学学报, 2013, 33(3): 0304001.
[24] 白廷柱, 金伟其. 光电成像原理[M]. 北京: 北京理工大学出版社, 2013: 127- 128.
Bai TZ, Jin WQ. Photoelectric imaging principle[M]. Beijing: Beijing Institute of Technology Press, 2013: 127- 128.
[25] Vallerga J V. Siegmund O H W. 2 K×2 K resolution element photon counting MCP sensor with >200 kHz event rate capability[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2000, 442(1/2/3): 159-163.
[26] DegnanJ. Impact of receiver deadtime on photon-counting SLR and altimetry during daylight operations[C]. 16th International Workshop on Laser Ranging, 2008: 339- 346.
[27] 许艺腾. 单光子激光测高数据处理技术研究[D]. 西安: 西安科技大学, 2017: 5- 6.
Xu YT. Research on data processing technology of single photon laser altimetry[D]. Xi'an: Xi'an University of Science andTechnology, 2017: 5- 6.
[28] 陈聪, 陈海燕. 洛伦兹与超高斯矩形超短脉冲的时间带宽积[J]. 长江大学学报(自然科学版), 2013, 10(1): 10-11.
[29] He KM, SunJ, Tang XO. Single image haze removal using dark channel prior[C]. Miami: 2009 IEEE Conference on Computer Vision and Pattern Recognition, 2009: 1956- 1963.
[30] Levin A, Lischinski D, Weiss Y. A closed-form solution to natural image matting[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2008, 30(2): 228-242.
Article Outline
米小什, 赵惠, 樊学武, 盛立志. 基于微通道板的单光子激光测高技术研究[J]. 光学学报, 2018, 38(12): 1228001. Xiaoshi Mi, Hui Zhao, Xuewu Fan, Lizhi Sheng. Photon-Counting Laser Altimetry Based on Microchannel Plate[J]. Acta Optica Sinica, 2018, 38(12): 1228001.