-
文中方法具体工作流程如图1所示,可划分为以下6个步骤:
1)读取波形数据。对ALB获取的LAS 1.3格式文件进行数据解析,得到波形数据和点云数据。
2)划分波形类别。结合实际影像和波形的形态,明确要分类的波形类型:异常波形、过拟合波形、陆地波形、海面波形和测深波形5种,波形的具体形态见图2。
3)提取波形形态特征。分析不同波形的形态特点,对波形进行特征提取。
4)评估波形特征。利用RFFS模型对提取的特征进行特征重要性排序。
5)自适应确定最佳波形特征组合。将得到的特征重要性排名以不同的组合输入分类器中,确定最佳的特征组合。
6)波形分类实验及评估。将最佳特征组合输入分类器进行分类实验和精度评估检验。
总体方法流程如图1所示。
-
目前全波形分类集中在海陆波形的划分[10,17]。纪雪等[10]尝试对除水、陆地波形以外的过拟合波形和孤立波形进行分类。但是,分类过程复杂,特征也较为常规,未针对不同类别的波形空间形态差异进行针对性提取。因此,文中在水陆分类过程中,尝试对干扰波形进行识别和剔除,以减少后续对信号的干扰和坐标解算的计算量。结合实际影像、波形形态,明确了要分类的5种波形类型:异常波形、过拟合波形、陆地波形、海面波形和测深波形,下面将对其分别进行介绍。
1)异常波形,如图2(a)所示,其形态整体呈现了先下降后上升再下降的趋势;水域、陆地范围都有分布,较为随机,没有规律。
2)过拟合波形,如图2(b)所示,分布在激光透射性较差或亮反射区域,如建筑物、柏油马路和光滑的海面(镜反射)等。从形态上可以直观地看出:相邻两点间振幅发生了剧烈的变化,这是由于目标物的回波信号超过了传感器的接收阈值。
3)陆地波形,如图2(c)所示,指激光最后发射到陆地再将信息反射回接收器接收到的完整波形。对于植被等透射性较好的地物,陆地回波可能出现多个波峰。
4)海面波形,如图2(d)所示,这里定义为未穿透水体,直接从海面返回接收器的波形信号,因此,海面波形只有一个波峰。
5)测深波形,如图2(e)所示,是激光到达海底所接收的波形信息。测深波形一般会有两个波峰,第一个波峰是海面反射的回波信号,一般峰值较高,第二个回波峰值是由海底反射的,一般较小。
-
波形特征提取是海陆波形分类的重要步骤。文中在常规特征参数的基础上总结了不同波形的形态和时域特征,并依据这些特征进行了创新和改进,共选取了24维特征参数。
1)最大最小振幅时间差
其中,异常波形在形态上与其他波形具有显著的差异。强度的最大值和最小值一般分布在首、尾位置,因此,最大最小振幅时间差${{T}}_{\mathrm{m}\mathrm{a}\mathrm{x}-\mathrm{m}\mathrm{i}\mathrm{n}}$可以较好地区分异常波形。
2)相邻两点间的振幅偏差
过拟合波形在形态上出现“截断”现象,有两个尖锐的波峰,如图3所示,基于此,提出了相邻两点间振幅偏差为$ \Delta $A。
3)最大强度
通常情况下,首个回波的信号强度最大,选取最大强度Wf-I特征,可以反映不同地物的属性。陆域回波信号通常具有更高的回波强度[17],此特征有助于区分陆地和海洋波形。
4)最大强度位置
激光脉冲回波信号,除了信号幅值外,脉冲的时间也是反映回波信号的主要特征量[18]。最大强度位置It表示回波信号最大强度所出现的时间。
5)峰值数量
不同目标物具有的峰值数量[10]Num不尽相同。例如,陆地波形多为一个,当遇到植被等由于光束透射会有多个回波;测深波形为两个峰值;海面回波波形有一个峰值。
6)首个回波峰值
文献[10]中的10维最重要特征主要分布在首个回波之中。因此,文中提取首个回波峰值Fpk表示多个回波峰值中的第一个回波的峰值数值。
7)首个回波峰值的位置
首个回波峰值的位置Ft主要描述了首个回波峰值所出现的时间,以进一步了解激光首次接触地物的时刻。
8)首个回波偏度
首个回波偏度Sfirst用于描述波形的偏移程度。观察测深波形形态多呈右偏态[19],因此提取此特征可用于区别其他波形。
9)半峰宽
此特征的方法是先找到最大的峰值计算出峰值高度一半的数值,然后根据这个数值计算出其出现的第一个和最后一个时间,作差即可得到半峰宽FW-HM[20]。
10)回波波宽
回波波宽Wfw是以数据形式表达不同地物波形的展宽程度。相比于陆地回波而言,海洋浅水和海洋深水的回波展宽程度较为明显[8]。其表达式为:
$$ {W f}_{w}=\frac{F W H M}{2\sqrt{2{\ln}2}} $$ (1) 11)衰减常数
衰减常数a[21],即首个峰值A1与第二个峰值A2之差,再将差值与首个峰值A1作商而得到的。如图4所示。此特征表征激光在水体中传输时的衰减[22],考虑到首个振幅、第二个振幅对区分波形的影响,其公式为:
$$ a=\frac{{A}_{1}-{A}_{2}}{{A}_{1}} $$ (2) 12)峰值之比
峰值之比Rp,即将最小峰值与最大峰值进行作商,一般首个峰值即为最大波峰,随着信号的逐渐减弱,最小峰值一般会出现在最后一个回波,因此,为探究不同波形的吸收、衰减作用提出了此特征。
13)震荡主频
震荡主频f,如图4所示,采用回波波形中的首个峰值和第二个峰值所出现的时间差取倒数而得到的[21]。其公式可表示为:
$$ f=\frac{1}{{t}_{2}-{t}_{1}} $$ (3) 式中:t1为首个峰值出现的时间;t2为第二个峰值出现的时间。
14)波形总能量
一条波形是由激光接触到地物所产生的回波和激光传输过程中产生的背景噪声所组成的,不同地物产生的能量各不相同。因此,波形总能量Sall表征回波中所有能量的累积。Sall的计算公式为:
$$ {S}{ }_{\rm all}=\frac{1}{2}{\sum }_{t=1}^{t=288}\left(W f\left(t\right)+W f\left(t+1\right)\right) $$ (4) 式中:Wf(t)为采用的原始波形;t=1表示波形包含噪声的开始时间为1 ns;t=288表示波形包含背景噪声的结束时间是288 ns。
15)有效回波能量
有效回波能量Sw,即除背景噪声外组成目标物的信号总能量[19]。此特征主要用来区分陆地、海面和测深回波所产生的有效能量。其公式可以表示为:
$$ {S}_{w}=\frac{1}{2}{\sum }_{t={t}_{\rm s}}^{t={t}_{\rm o}}\left(W f\left(t\right)+W f\left(t+1\right)\right) $$ (5) 式中:ts表示波形有效范围的开始时间;to表示波形有效范围的结束时间。
16)波形有效范围
波形的有效范围L是描述回波信号的持续时间,在一定程度上反映信号的重叠情况[17],如图5所示。L的计算公式为:
$$ L={t}_{\rm o}-{t}_{\rm s} $$ (6) 17)峰值20%以上的面积
选取最大峰值的20%这部分回波的能量[17],如图5所示,其公式可以表示为:
$$ {S}_{20 \text{%}}={\int }_{t={t}_{\rm 2s}}^{t={t}_{\rm 2o}}W f\left(t\right){\rm{d}}t $$ (7) 式中:t2s表示波形峰值20%的开始时间;t2o表示波形峰值20%的结束时间。
18)最大振幅的10%所持续的时间
最大振幅的10%所持续的时间[23]T10%,即利用最大振幅的10%作为阈值,确定阈值范围的时间差即为T10%,如图5所示。
19)最大振幅的50%所持续的时间
最大振幅的50%所持续的时间[23]T50%,与最大振幅的10%所持续时间T10%相似,如图5所示。
20)脉冲上升沿10%~90%的时间
选取波形数据中第一个大于最大峰值10%和第一个大于最大峰值90%的数据所对应的时间差值作为脉冲上升沿10%~90%的时间[24]Tup1,如图6所示。
21)脉冲上升沿10%~90%的陡度
脉冲上升沿10%~90%的陡度[24]ST1选用了80%的回波信号最大强度与脉冲上升部分的10%~90%的时间差作商。
$$ {S T}_{1}=\frac{0.8\times W f{\text{-}}I}{{T}_{\rm up1}} $$ (8) 式中:Wf-I表示回波信号最大强度;Tup1表示脉冲上升沿10%~90%的持续时间。
22)脉冲下降沿10%~90%的时间
脉冲下降沿10%~90%的时间[24]Tdown1可理解为波形下降的部分选取最大幅值的10%~90%对应的时间列的时间差即为脉冲下降沿10%~90%的时间Tdown1,如图6所示。
23)脉冲上升沿30%~90%的时间
此特征与Tup1特征相似,即将10%~90%换为30%~90%即为脉冲上升沿30%~90%的时间[24]Tup3,如图6所示。
24)脉冲上升沿30%~90%的陡度
选取脉冲上升波形部分的30%~90%的数据计算其数值变化情况[24],ST3的公式为:
$$ {S T}_{3}=\frac{0.6 \times W f {\text{-}} I}{{T}_{\rm up3}} $$ (9) 式中:Tup3表示脉冲上升沿30%~90%的持续时间。
-
随机森林(Random Forest,RF)是Breiman于2001年提出的一种统计学习理论方法[25-27],属于基于决策树的分类器[25]。目前,随机森林在分类、回归和降维等问题方面得到了广泛应用。
RFFS算法[26]对于重要测度指标的计算主要有两种:基于Gini指数的平均不纯度减少指标(Mean De-crease Impurity,MDI)和基于OOB数据置换的平均精度准确率减少指标(Mean Decrease Accuracy,MDA)[27]。利用第二种方式计算模型中各特征的重要性,首先在原始训练样本中通过Bootstrap的方式有放回地进行样本的抽取。假设训练数据中共有N个对象的数据,从这些数据样本中每次有放回地取出一些样本数据,每次取出的并不完全相同,这些样本将作为决策树的训练数据组,这些训练数据组被分为n个样本集,每个样本集生成一棵决策树;另外1/3的原始样本并未被抽取,这部分样本称为袋外(Out of Bag, OOB)数据,可用于预测错误率E1;然后随机置换变量$ i $的观测值构建决策树,计算袋外数据的预测错误率E2i,最后计算OOB误差($ {OOB}_{\mathrm{e}\mathrm{r}\mathrm{r}\mathrm{o}\mathrm{r}} $),其公式为:
$$ {OOB}_{\mathrm{e}\mathrm{r}\mathrm{r}\mathrm{o}\mathrm{r}}={E2}_{i}-E1 $$ (10) 使用OOB误差进行模型评估以确定模型最优特征数量。通过计算不同特征组合模型的OOB误差,综合考虑模型精度和计算复杂度来确定最优特征个数[28]并得出变量重要性排序。特征$ i $的置换重要性importance(i)公式为:
$$ importance\left(i\right)=\sum {OOB}_{\mathrm{e}\mathrm{r}\mathrm{r}\mathrm{o}\mathrm{r}}/T $$ (11) 式中:T表示分类树的数量。
-
随机森林是一种非参数分类与回归方法,不需要先验知识,易于使用,能在保证良好精度的前提下并行化处理高维海量数据,是一种高效的机器学习算法[28]。
随机森林决策树构建中的一个关键步骤——节点分裂时的特征选择:假设每个数据点都有M个特征,则可以在每个节点处随机选择m(m<M)个特征,选用信息熵、信息增益或者基尼指数来选择最好的特征进行节点分裂。然后将m分成两部分,用最佳分裂方案对节点进行分裂,决策树生长期内m大小始终不变,每棵树均进行完全的生长且不被修剪。最后,将上述步骤进行重复,通过k次样本抽取和样本训练得到k个决策树模型,利用集成学习理论将这k个决策树进行线性组合,其中每个决策树占相等的权重。组合多棵决策树,从每个决策树中获取预测结果,对每个预测结果进行投票,通过投票得出最终的预测结果。
分类问题中CART使用$ Gini $值作为评判标准。定义为:
$$ Gini=1-\sum \left(P(i)\times P\left(i\right)\right) $$ (12) 式中:$ P(i) $为当前节点上数据集中第$ i $类样本的比例。
-
文中的研究区域位于中国海南省三亚市崖州区的东锣岛,其中扫描条带覆盖东锣岛范围,如图7所示。东锣岛陆域面积约0.1180 km2,地理坐标为北纬18°19′33.20″,东经108°59′24.13″,岸线长约1.432 km,其附近的水深约为4.3~10.2 m。研究区数据由Optech公司设计的Aquarius型号机载激光雷达进行采集。此型号的机载激光雷达专门为沿海和内河航道测绘而设计。
-
文中在总数据集中随机选取了10万个样本数据,每个类别各随机选取了1万个训练集,1万个测试集。对每个波形都提取了24个特征参数,生成了1×24的特征矩阵,输入到随机森林特征重要性评估模型中进行特征重要性排序。为避免特征冗余,增强所选特征的可靠性,随机选取每类波形的样本数据各1万条进行了10次随机森林特征重要性实验,根据10次实验的权重系数取平均值得到最终的特征重要性排名,如图8所示。
-
在大量特征中选取最可靠的特征,是提高分类精度和表达地物类型空间分布的关键。根据特征重要性排序结果选择不同特征进行组合、分类和精度评估以确定最优的特征组合,如图9所示。
从图9可以看出,Top15由于特征的冗余,精度最低;Top10较Top15精度有所上升,说明11~15维特征对于5种波形的区分效果不明显;Top9与Top10相较,精度有0.32%的下降,这说明第10维特征对此随机数据的精度有贡献;Top9~Top5随着特征数量的降低,精度在逐渐提高,最大可达到98.10%。特征重要性前5维特征和第10维特征对分类的贡献最高,重新进行随机实验,精度达到了98.55%。震荡主频f、峰值之比Rp、相邻两点间振幅偏差$ \Delta{A} $、最大强度Wf-I、衰减常数$ a $、首个回波峰值Fpk作为最终的特征优选结果。其在5种波形上的箱型分布如图10所示。
-
为了进一步检验随机森林分类方法的性能,选取了K最邻近值(KNN)和支持向量机(SVM)分类算法进行对比,分类结果如表1所示。
表 1 不同分类方法结果对比
Table 1. Comparison of the results of different classi-fication methods
Waveform types RF SVM KNN Anomalies waveform 100% 100% 100% Over-fitted waveform 100% 100% 100% Land waveform 99.46% 95.9% 93.42% Sea surface waveform 94.29% 84.5% 72.86% Bathymetry waveform 99.02% 90% 88.77% OA 98.55% 94% 91.01% 实验结果表明:RF的分类效果最好,总体分类准确率(OA)为98.55%。其中,异常波形和过拟合波形的分类精度均为100%,陆地波形的分类精度为99.46%,海面波形的分类精度为94.29%,测深波形的分类精度为99.02%,Kappa系数为0.9820。OA和Kappa系数的计算公式分别为:
$$ OA=\frac{\displaystyle\sum _{i=1}^{n}{n}_{ii}}{N} $$ (13) $$ Kappa=\frac{N\displaystyle\sum _{i=1}^{n}{n}_{ii}-\displaystyle\sum _{i=1}^{n}{n}_{i.}{n}_{.i}}{N\times N-\displaystyle\sum _{i=1}^{n}{n}_{i.}{n}_{.i}} $$ (14) 式中:$ N $是波形总数;$ {n}_{ii} $代表对角线上的总波形数;$ {{n}}_{i.} $则是第$ i $行的总波形数;$ {n}_{.i} $则代表第$ i $列的总波形数。
测试数据的分类结果相对应的混淆矩阵如图11所示。研究区的波形分类结果如图12所示。其中,过拟合波形分布呈现出一定的航带中央集中性,这主要是由于水域范围内,过拟合波形的产生主要是由于发生海面镜面反射造成的高返回信号,这与激光发射角度有一定关系,随着发射角增大,镜面反射现象逐渐弱化,因此记录的过拟合波形减少;而在陆地附近的过拟合波形分布主要是由于地物类别引起的,对于沙滩、裸露礁石、建筑物、柏油马路等地物得到高返回信息,因此过拟合波形较多,这使得在陆地区域不存在过拟合波形航带中央集中特性。
为了验证文中分类的准确性、模型的迁移性和普适性,另外选择了一个实验区域的条带数据进行波形分类,结果如图13所示。从图中可以看出:过拟合波形主要分布在裸露地面和建筑物上,异常波形的分布比较均匀。任意选取10块区域进行海陆分类精度验证,海洋的分类精度为97.43%,陆地的分类精度为96.2%,水陆总体分类精度为96.81%,满足生产需求。但是,精度发生了大约2%的降低,主要原因是研究区和实验区的数据并不是同一时间、同一区域采集的,而且飞行条件和飞行高度等参数设置不一致、区域的地物类别引起的地物振幅的变化等造成了一定的偏差。
Classification of full waveform data for monochromatic airborne LiDAR bathymetry based on waveform morphological features
-
摘要: 单频机载激光雷达测深系统凭借低成本、低负载、高采样率等优势成为大范围海岸带地形地貌探测的理想选择。然而如何解决单频局限,在不依赖辅助传感器情况下实现全波形数据的准确划分成为精确点位坐标解算的关键环节。目前基于全波形形态特征进行波形分类研究缺乏系统性评估分析和普遍性结论。该研究尝试从全波形空间形态入手,细化了波形类别(异常波形、过拟合波形、陆地波形、海面波形和测深波形),在已有波形特征基础上,系统分析了不同类别波形的形态特征差异,有针对性地提取了24维波形特征并基于随机森林特征选择和分类模型完成了各特征分类性能及最佳特征组合评估与定量分析。研究证明,包括相邻两点间振幅偏差、震荡主频等在内的6维特征组合对5种波形的分类效果最好,总体分类精度可达98.55%,Kappa系数为0.9820。为了验证特征的普适性,另外选取了一块实验区域进行验证,得到水陆分类的总体精度为96.81%。Abstract:
Objective Monochromatic airborne LiDAR bathymetry becomes considerably favorable for topography and geomorphology detection over coastal area by means of its low cost, low load and high sampling rate. However, addressing the limitation of single wavelength to realize the accurate division of full waveform data independently from auxiliary sensor becomes the critical part of coordinate calculation. Given the existing literatures, there is a lack of systematic evaluation analysis and general conclusions for waveform classification contraposing to full waveform morphological features. Methods In view of the latest development of waveform features extraction, refined waveform categories (anomalies, over-fitted, land, sea surface and bathymetry waveforms), 24-dimensional waveform features are designed and calculated upon systematic analysis on morphological characteristics of different waveforms, and then their classification performance and optimal feature combination are evaluated and quantitatively analyzed utilizing random forest feature selection and classification model. Results and Discussions The results proved that the combination of 6-dimensional features (Fig.8-11), including deviation of amplitude between two adjacent points and oscillating main frequency, is the most effective in classifying five waveforms, with an overall classification accuracy of 98.55% and a Kappa coefficient of 0.982 0 (Fig.9-10, Fig.12, Tab.1). To verify the universality of the features, an additional experimental area was selected for validation and the overall accuracy of water and land classification was 96.81% (Fig.13). Conclusions To accurately identify waveforms, a systematic analysis was conducted to determine the morphological differences between different types of waveforms, and 24-dimensional feature parameters were extracted. After the optimal feature combination and classification performance evaluation, it was found that the 6-dimensional features of oscillating main frequency f, ratio of peak Rp, deviation of amplitude between two adjacent points ∆A, maximum intensity Wf-I, decay constant a, and first echo peak Fpk were highly effective in distinguishing these five types of waveforms, where 100% of the anomalies and over-fitted waveforms were extracted, strongly confirming the relevance and validity of the morphological features. After replacing the experimental area, the accuracy of the water and land classification reached 96.81%, proving that the features and methods used were adaptable and generalizable, and could meet the production requirements. The 2% decrease in waveform classification accuracy after changing the study area is mainly due to the varying equipment parameter settings and coverage feature categories in different experimental areas. Limited sample selection further compounds this issue. To maintain accuracy, the sample data can be appropriately supplemented according to the actual situation in the experimental area. Although the waveform morphology has been studied thoroughly, additional experimental evidence is necessary to ascertain the impact of intrinsic factors such as signal reception systems on the echoes. To this end, future research will focus on expanding the study's data and signal detection. -
表 1 不同分类方法结果对比
Table 1. Comparison of the results of different classi-fication methods
Waveform types RF SVM KNN Anomalies waveform 100% 100% 100% Over-fitted waveform 100% 100% 100% Land waveform 99.46% 95.9% 93.42% Sea surface waveform 94.29% 84.5% 72.86% Bathymetry waveform 99.02% 90% 88.77% OA 98.55% 94% 91.01% -
[1] 刘永明, 邓孺孺, 秦雁, 等. 机载激光雷达测深数据处理与应用[J]. 遥感学报, 2017, 21(06): 982-995. Liu Yongming, Deng Ruru, Qin Yan, et al. Data processing methods and applications of airborne LiDAR bathymetry [J]. Journal of Remote Sensing, 2017, 21(6): 982-995. (in Chinese) [2] 陈坚, 金翔龙. 机载激光测深技术进展及应用[J]. 海洋通报, 2002(06): 75-82+90. doi: 10.3969/j.issn.1001-6392.2002.06.011 Chen Jian, Jin Xianglong. Advances in and application of airborne laser bathymetry technology [J]. Marine Science Bulletin, 2002(6): 75-82, 90. (in Chinese) doi: 10.3969/j.issn.1001-6392.2002.06.011 [3] 翟国君, 王克平, 刘玉红. 机载激光测深技术[J]. 海洋测绘, 2014, 34(02): 72-75. doi: 10.3969/j.issn.1671-3044.2014.02.021 Zhai Guojun, Wang Keping, Liu Yuhong. Technology of airborne laser bathymetry [J]. Hydrographic Surveying and Charting, 2014, 34(2): 72-75. (in Chinese) doi: 10.3969/j.issn.1671-3044.2014.02.021 [4] 李洪鹏, 李国元, 蔡志坚, 等. 全波形激光雷达回波分解方法[J]. 遥感学报, 2019, 23(01): 89-98. Li Hongpeng, Li Guoyuan, Cai Zhijian, et al. Full-waveform LiDAR echo decomposition method [J]. Journal of Remote Sensing, 2019, 23(1): 89-98. (in Chinese) [5] 刘智敏, 杨安秀, 阳凡林, 等. 机载LiDAR测深在海洋测绘中应用的可行性分析[J]. 海洋测绘, 2018, 38(04): 43-47. doi: 10.3969/j.issn.1671-3044.2018.04.009 Liu Zhimin, Yang Anxiu, Yang Fanlin, et al. The feasibility analysis for the airborne LiDAR bathymetry application in marine surveying and charting [J]. Hydrographic Surveying and Charting, 2018, 38(4): 43-47. (in Chinese) doi: 10.3969/j.issn.1671-3044.2018.04.009 [6] 赵建虎, 欧阳永忠, 王爱学. 海底地形测量技术现状及发展趋势[J]. 测绘学报, 2017, 46(10): 1786-1794. doi: 10.11947/j.AGCS.2017.20170276 Zhao Jianhu, Ouyang Yongzhong, Wang Aixue. Status and development tendency for seafloor terrain measurement technology [J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(10): 1786-1794. (in Chinese) doi: 10.11947/j.AGCS.2017.20170276 [7] 胡海瑛, 惠振阳, 李娜. 基于多基元特征向量融合的机载LiDAR点云分类[J]. 中国激光, 2020, 47(08): 237-247. doi: 10.3788/CJL202047.0810002 Hu Haiying, Hui Zhenyang, Li Na. Airborne LiDAR point cloud classification based on Multiple-entity eigenvector fusion [J]. Chinese Journal of Lasers, 2020, 47(8): 0810002. (in Chinese) doi: 10.3788/CJL202047.0810002 [8] 黄田程, 陶邦一, 毛志华, 等. 基于多通道海洋激光雷达的海陆波形分类[J]. 中国激光, 2017, 44(06): 294-303. doi: 10.3788/CJL201744.0610002 Hu Tiancheng, Tao Bangyi, Mao Zhihua, et al. Classification of sea and land waveform based on multi-channel ocean lidar [J]. Chinese Journal of Lasers, 2017, 44(6): 0610002. (in Chinese) doi: 10.3788/CJL201744.0610002 [9] 邓迁, 吴德成, 况志强, 等. 用于水汽混合比自标定的532 nm/660 nm双波长激光雷达[J]. 红外与激光工程, 2018, 47(12): 102-106. doi: 10.3788/IRLA201847.1230004 Deng Qian, Wu Decheng, Kuang Zhiqiang, et al. 532 nm/660 nm dual wavelength lidar for self-calibration of water vapor mixing ratio [J]. Infrared and Laser Engineering, 2018, 47(12): 1230004. (in Chinese) doi: 10.3788/IRLA201847.1230004 [10] Ji X, Tang Q, Xu W, et al. Island feature classification for single-wavelength airborne lidar bathymetry based on full-waveform parameters [J]. Applied Optics, 2021, 60(11): 3055-3061. [11] 张扬, 黄卫东, 董长哲, 等. 海洋激光雷达探测卫星技术发展研究[J]. 红外与激光工程, 2020, 49(11): 28-39. doi: 10.3788/IRLA20201045 Zhang Yang, Huang Weidong, Dong Changzhe, et al. Research on the development of the detection satellite technology in oceanographic lidar [J]. Infrared and Laser Engineering, 2020, 49(11): 20201045. (in Chinese) doi: 10.3788/IRLA20201045 [12] Fernandez-Diaz J C, Glennie C L, Carter W E, et al. Early results of simultaneous terrain and shallow water bathymetry mapping using a single-wavelength airborne LiDAR sensor [J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2014, 7(2): 623-635. doi: 10.1109/JSTARS.2013.2265255 [13] 曹彬才, 朱述龙, 邱振戈, 等. 利用回波波形特征实现激光雷达的水陆分离[J]. 海洋测绘, 2018, 38(03): 12-16. Cao Bincai, Zhu Shulong, Qiu Zhenge, et al. Water-land classification for LiDAR bathymetric data based on echo waveform characteristics [J]. Hydrographic Surveying and Charting, 2018, 38(3): 12-16. (in Chinese) [14] Guenther G C, LaRocque P E, Lillycrop W J. Multiple surface channels in Scanning Hydrographic Operational Airborne Lidar Survey (SHOALS) airbornZe lidar[C]//Proceedings of SPIE, Ocean Optics XII, 1994, 2258: 1-9. [15] 王金虎, 李传荣, 周梅. 全波形激光雷达数据在点云分类中的应用研究[J]. 遥感信息, 2013, 28(05): 21-27. Wang Jinhu, Li Chuanrong, Zhou Mei. Analysis of airborne full-waveform LiDAR data for supervised point cloud classification [J]. Remote Sensing Information, 2013, 28(5): 21-27. (in Chinese) [16] 范士俊. 基于随机森林的全波形点云数据分类研究[J]. 首都师范大学学报(自然科学版), 2013, 34(05): 71-73+78. Fan Shijun. Research on classification for airborne full waveform Lidar data based on random forest [J]. Journal of Capital Normal University (Natural Science Edition), 2013, 34(5): 71-73, 78. (in Chinese) [17] 王丹菂, 邢帅, 徐青, 等. 单频机载激光测深海陆回波自动分类方法[J]. 测绘学报, 2022, 51(05): 750-761. doi: 10.11947/j.issn.1001-1595.2022.5.chxb202205013 Wang Dandi, Xing Shuai, Xu Qing, et al. Automatic sea-land waveform classification method for single-wavelength airborne LiDAR bathymetry [J]. Acta Geodaetica et Cartographica Sinica, 2022, 51(5): 750-761. (in Chinese) doi: 10.11947/j.issn.1001-1595.2022.5.chxb202205013 [18] 杨雨川, 谭碧涛, 龙超, 等. 主动探测回波激光脉冲时域特性[J]. 红外与激光工程, 2013, 42(12): 3228-3233. doi: 10.3969/j.issn.1007-2276.2013.12.013 Yang Yuchuan, Tan Bitao, Long Chao, et al. Time-domain characteristics of active detecting echo laser pulse [J]. Infrared and Laser Engineering, 2013, 42(12): 3228-3233. (in Chinese) doi: 10.3969/j.issn.1007-2276.2013.12.013 [19] 左志强, 唐新明, 李国元, 等. GF-7星载激光测高仪全波形自适应高斯滤波[J]. 红外与激光工程, 2020, 49(11): 132-142. doi: 10.3788/IRLA20200251 Zuo Zhiqiang, Tang Xinming, Li Guoyuan, et al. Adaptive Gaussian filtering of the full waveform of GF-7 satellite laser altimeter [J]. Infrared and Laser Engineering, 2020, 49(11): 20200251. (in Chinese) doi: 10.3788/IRLA20200251 [20] 徐志钮, 胡志伟, 赵丽娟, 等. 采用Voigt模型的布里渊散射谱关键特征高精度提取方法(英文)[J]. 红外与激光工程, 2018, 47(S1): 74-81. doi: 10.3788/IRLA201746.S122004 Xu Zhiniu, Hu Zhiwei, Zhao Lijuan, et al. Highly accurate key parameters extraction algorithm for Brillouin scattering spectrum using Voigt profile [J]. Infrared and Laser Engineering, 2018, 47(S1): S122004. (in Chinese) doi: 10.3788/IRLA201746.S122004 [21] 吴田, 付道睿, 彭勇, 等. 1000 kV输电线路耐张塔带电作业电位转移电流脉冲波形特征参数研究[J]. 电网技术, 2022, 46(05): 2009-2016. Wu Tian, Fu Daorui, Peng Yong, et al. Characteristic parameters of pulse waveform of potential transfer current for live working on tension tower of 1000 kV transmission line [J]. Power System Technology, 2022, 46(5): 2009-2016. (in Chinese) [22] 徐沛拓, 陶雨婷, 刘志鹏, 等. 海洋激光雷达实验与仿真结果的对比[J]. 红外与激光工程, 2020, 49(02): 71-77. doi: 10.3788/IRLA202049.0203007 Xu Peituo, Tao Yuting, Liu Zhipeng, et al. Comparison of oceanic lidar experiments and simulation results [J]. Infrared and Laser Engineering, 2020, 49(2): 0203007. (in Chinese) doi: 10.3788/IRLA202049.0203007 [23] 王国利, 郝艳捧, 李彦明. 油中典型局部放电模型放电波形特征参数的提取[J]. 电网技术, 2001(10): 34-37+46. doi: 10.3321/j.issn:1000-3673.2001.10.008 Wang Guoli, Hao Yanpeng, Li Yanming. Pulse current parameter evaluation of typical partial discharge models in oil [J]. Power System Technology, 2001(10): 34-37, 46. (in Chinese) doi: 10.3321/j.issn:1000-3673.2001.10.008 [24] 赵阳, 郄秀书, 孔祥贞, 等. 人工触发闪电电流波形特征参数分析[J]. 物理学报, 2009, 58(09): 6616-6626. doi: 10.7498/aps.58.6616 Zhao Yang, Qie Xiushu, Kong Xiangzhen, et al. Analysis on the parameters of the current waveforms of triggered lightning [J]. Acta Physica Sinica, 2009, 58(9): 6616-6626. (in Chinese) doi: 10.7498/aps.58.6616 [25] 孙杰, 赖祖龙. 利用随机森林的城区机载LiDAR数据特征选择与分类[J]. 武汉大学学报(信息科学版), 2014, 39(11): 1310-1313. Sun Jie, Lai Zulong. Airborne LiDAR feature selection for urban classification using random forests [J]. Geomatics and Information Science of Wuhan University, 2014, 39(11): 1310-1313. (in Chinese) [26] 熊艳, 高仁强, 徐战亚. 机载LiDAR点云数据降维与分类的随机森林方法[J]. 测绘学报, 2018, 47(04): 508-518. Xiong Yan, Gao Renqiang, Xu Zhanya. Random forest method for dimension reduction and point cloud classification based on airborne LiDAR [J]. Acta Geodaetica et Cartographica Sinica, 2018, 47(4): 508-518. (in Chinese) [27] 宋述芳, 何入洋. 基于随机森林的重要性测度指标体系[J]. 国防科技大学学报, 2021, 43(02): 25-32. doi: 10.11887/j.cn.202102004 Song Shufang, He Ruyang. Importance measure index system based on random forest [J]. Journal of National University of Defense Technology, 2021, 43(2): 25-32. (in Chinese) doi: 10.11887/j.cn.202102004 [28] 何云, 黄翀, 李贺, 等. 基于Sentinel-2 A影像特征优选的随机森林土地覆盖分类[J]. 资源科学, 2019, 41(05): 992-1001. He Yun, Huang Chong, Li He, et al. Land-cover classification of random forest based on Sentinel-2A image feature optimization [J]. Resources Science, 2019, 41(5): 992-1001. (in Chinese)