一种用于TLS多回波点云植被遮挡目标的探测和分析方法与流程

本发明涉及多回波点云探测领域，具体涉及一种用于tls多回波点云植被遮挡目标的探测和分析方法。

背景技术：

随着大场景三维激光扫描技术不断完善，地面激光扫描仪(terrestriallaserscanners,简称tls)成为继传统遥感技术之后重要的可穿透植被叶簇、且可感知与检测其遮挡目标的一种新型测距传感器。它通过向目标发射激光脉冲并接收其反射回来的回波信号，在三维空间上形成一系列点的集合，即多回波点云。每个点都具有三维坐标和强度，以此获取目标的几何位置和反射率。对于植被遮挡场景，用tls获取点云数据，可以对被树丛遮挡的隐蔽目标进行感知与提取。

tls点云数据以其高密度、高精度和高效率等优势，成为了一种精确获取空间数据的快速手段。首先，tls采用单色光且发射波束极窄，隐蔽性好，对地物和背景具有极强的抑制能力；其次，激光对隐身目标具有极高的灵敏度，抗干扰能力突出；同时，激光扫描系统波长短，与微波扫描系统相比其体积和质量都比较小；另外，它有极高的分辨率，用其做精确制导导引头具有很高的精度。

tls以一定的水平角度和垂直角度向目标发射激光脉冲，然后记录收到回波的时间，根据几何关系进行计算得出每个扫描点的位置、反射强度等信息。多个点构成一幅具有三维空间坐标的点的集合，即为点云数据。tls全回波系统数字化记录单次发出脉冲的后向散射信号，因此，用户可以决定扫描场景的范围。当激光光束的光斑打击在可穿透物体(如树叶叶片边缘)时，另一部分光斑继续穿透叶簇，落在更深处的叶片上发生反射，如此继续重复，从而产生了多次回波。

单次回波表示激光束打到该点只产生一次反射，表示该处不可穿透。一次回波主要位于目标边缘，从而可以进行下一次反射。中间回波基本上只有植被树叶冠层，这是因为树叶与树叶之间有缝隙，激光束经过树叶边缘一部分产生一次反射，另一部分继续深入缝隙中获取树冠层深层次信息，从而产生中间回波。而末次回波也表示该处不再被穿透，主要为刚性物体，其中包含了可能被植被遮挡的隐蔽目标。

目前三维点云遮挡目标探测主要存在以下问题：

1、植被遮挡目标(即隐蔽目标)因其遮挡所展现出来的形态不完整、不规则，因此激光在其表面进行散射，波形背景复杂；

2、目标的局部结构可能存在判别力弱、不具有区分性等问题，因而单一的局部特征不能完整描述目标状态；

3、由于植被遮挡目标形态的多样性以及和扫描仪角度的随机性，且回波信号受大气效应、距离范围和噪声等变化的影响，点云遮挡检测对目标特性建模提出了较高的要求。

同时，现有技术中缺乏对基于tls点云差分位置、tls点云密度和tls点云强度实现了植被遮挡目标(即隐蔽目标)的可见性分析，无法对探测模型及探测结果做出量化的评价。

技术实现要素：

本发明的目的在于提供一种用于tls多回波点云植被遮挡目标的探测方法。

为实现上述目的，本发明采用以下技术方案：

一种用于tls多回波点云植被遮挡目标的探测方法，包括：

s1、对tls多回波点云进行点云预处理；

s2、将预处理后的点云处理成数据集输入神经网络，进行神经网络的训练；

s3、利用训练好的神经网络识别中间回波，根据识别的中间回波，找出植被遮挡冠层，以此确定植被遮挡区域；

s4、利用末次回波滤除植被遮挡区域中的非植被遮挡目标，以此探测植被遮挡目标。

进一步地，步骤s1具体为，

s11、采用向上生长法进行地面滤除；

s12、采用基于体素的向上交错分析法进行植被提取。

进一步地，步骤s2具体为，将点云投影到xy平面，结合点云的强度信息，构成数据集；将数据集输入神经网络，进行神经网络的训练。

进一步地，步骤s4具体为，利用s2训练好的神经网络识别单次回波、一次回波以及末次回波，将单次回波、一次回波、中间回波以及末次回波分别标记为网络的第一层、第二层、第三层以及第四层，训练基于kd树思想的神经网络，学习植被遮挡区域末次回波的形成过程，以区分非植被遮挡区域的末次回波，从而探测末次回波中的植被遮挡目标。

本发明还提供一种用于tls多回波点云植被遮挡目标的分析方法，用于对如前述的探测方法获得的探测结果进行可见性分析，其在所述步骤s4之后进行步骤s5:

s5、对植被遮挡目标，进行可见性分析，其包括：

s51、对植被遮挡目标进行差分位置分析；

s52、对植被遮挡目标进行回波密度分析；

s53、对植被遮挡目标进行回波强度分析。

进一步地，步骤s51中，点云差分位置分析具体为，

s511、计算植被遮挡区域和植被遮挡目标相邻点之间的三维坐标差值，获得差分位置矩阵；

s512、将得到的差分位置矩阵求特征值，得到复数矩阵，分别用quiver，feather，compass，angle函数，得到植被遮挡区域和植被遮挡目标的差分位置矩阵特征值可视化图形；

s513、对三维坐标差值做二次幂变换和log函数变换，得到植被遮挡区域和植被遮挡目标的二次幂矩阵特征值和特征向量的可视化图形。

进一步地，步骤s52具体为，

s521、获取植被遮挡区域与植被遮挡目标在x、y、z上的点云密度值，并进行归一化；

s522、计算多回波ndi指数，则有，

其中，n一次回波代表一次回波的距离范围均值，n末次回波代表末次回波的距离范围均值；

s523、以归一化的点云密度值为横坐标，以ndi指数为纵坐标，获得ndi指数随点云密度变化的示意图。

进一步地，步骤s53具体为，

s531、以点云强度值为横坐标，以点个数为纵坐标，同时，以不同颜色区分不同强度值，获得植被遮挡区域与植被遮挡目标的强度特性图；

s532、以xy平面植被遮挡区域的差分位置距离为横坐标，以植被遮挡目标的点云强度为纵坐标，获得二者的相关性二维特性图；

s533、分别以植被遮挡区域在x轴与y轴方向上的差分位置距离dx和dy为x轴与y轴，以植被遮挡目标的点云强度为纵坐标，建立xy平面上差分位置与点云强度的相关性三维特性图。

采用上述技术方案后，本发明与背景技术相比，具有如下优点：

1、本发明提出的全波形tls激光束为小光斑，小光斑具有对形状特征刻画精细的特点。小光斑和激光点云的特征层融合提供了新的理论框架，将提高了三维对象类检测识别的类别数量和识别性能；同时，本发明把全波段分解成多个独立回波，然后基于神经网络对自动对回波进行识别，从而实现相关隐蔽目标的探测。

2、本发明通过建立分解拟合，多回波数据和点云局部提取特征，克服配准、遮挡等因素带来的影响，提高对三维空间对象的局部描述能力，实现不确定性联合关系的匹配。

3、本发明把全波段分解成多个独立回波，然后基于tls多回波点云的位置、强度和能量等信息获取回波特性，可为边缘检测及树木等大场景目标提供重要研究基础。

4、本发明能在最小人工干预条件下实现tls点云数据的智能化场景分析和解译,研究成果将在森林调查、隐蔽目标分析、军事作战以及救援救灾等方面发挥巨大作用，具有重要的学术价值和现实意义。

附图说明

图1为本发明流程图。

图2为点云预处理的示意图，图2(a)为原始点云，图2(b)为地面滤除与植被提取后的示意图。

图3为神经网络的训练表现示意图，浅色为损失率，深色为精确率。

图4植被遮挡区域及其植被遮挡目标的探测结果示意图，其中深色部分为隐蔽目标的探测结果。

图5点云目标位置差分特性图，纵坐标对应位置差分幅度，横坐标对应点云目标，图5(a)对应植被遮挡区域，图5b对应植被遮挡目标。

图6植被遮挡区域的差分位置矩阵特征值可视化图，图6(a)-图6(b)分别对应quiver、feather、compass、angle函数。

图7为植被遮挡目标的差分位置矩阵特征值可视化图，图7(a)-图7(b)分别对应quiver、feather、compass、angle函数。

图8为目标位置差值二次幂矩阵特征值和特征向量可视化图，图8(a)对应植被遮挡区域，图8(b)对应植被遮挡目标。

图9为植被遮挡区域的密度特性图，图9(a)-图9(c)分别对应x、y及z轴方向。

图10为植被遮挡目标的密度特性图，图10(a)-图10(c)分别对应x、y及z轴方向。

图11为植被ndi指数分析图。

图12为植被遮挡目标的可见性分析图，图12(a)为正面，图12(b)为斜侧面。

图13为强度特性图，图13(a)对应植被遮挡区域，图13(b)对应植被遮挡目标。

图14为xy平面上差分位置与能量强度的相关图，图14(a)为隐蔽目标点云强度与xy平面遮挡目标差分位置距离相关性特性图，图14(b)为植被遮挡区域在x轴和y轴方向上差分位置距离与点云强度的相关性特性图；图14(c)对图14(b)根据点云强度进行划分；图14(d)为建立xy平面上差分位置与点云强度的相关性的三维特性图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

请参考图1所示，本发明包括两块核心内容，分别为植被遮挡目标探测与植被遮挡目标分析。下面结合具体实施例对二者进行论述。

实施例1

一种用于tls多回波点云植被遮挡目标的探测方法主要包括以下步骤：

s1、对tls多回波点云进行点云预处理。

根据全波形tls获取的多回波点云，为了消除冗余信息，用向上生长法进行地面滤除，然后用基于体素的向上交错分析法进行植被提取。原始点云如图2(a)所示,滤除冗余地面以及植被提取结果如图2(b)所示。

s2、将预处理后的点云处理成数据集输入神经网络，进行神经网络的训练。

将点云分别投影到xy、yz及xz平面，结合点云的强度信息，构成点云的不同三维属性(x，y，i)、(y，z，i)和(x，z，i)并进行训练，其后输入神经网络进行训练。经实验论证，(x，y，i)的训练效果最好，损失率和精确率结果如图3所示，因此，本实施例采用(x，y，i)数据集进行训练。

s3、利用训练好的神经网络识别中间回波，根据识别的中间回波，找出植被遮挡冠层，以此确定植被遮挡区域。

网络训练完成以后，就可以对tls多回波中的单次回波、一次回波、中间回波和最后回波进行识别了。由于本发明需要提取植被遮挡区域(即植被遮挡冠层)，因此主要针对中间回波进行识别，中间回波识别结果合实际中间回波进行对比，发现精确率达到90％以上。

其后，利用deepkd-networks[10]的生成比拟树木遮挡区域tls多回波点云的生成。将单次回波、一次回波、中间回波和末次回波分别标记为网络的第一层(叶节点)、第二层、第三层和最四层。训练类似于kd树的神经网络，让其学习植被遮挡区域末次回波的形成过程，以区别非植被遮挡区域的末次回波，从而滤除非植被遮挡目标(即显性目标)，探测末次回波里面的植被遮挡目标(即隐蔽目标)。植被遮挡区域与植被遮挡目标的探测结果如图4所示，其中，深色部分为植被遮挡目标。如此，即可完成植被遮挡目标的探测。

实施例2

一种用于tls多回波点云植被遮挡目标的分析方法在实施例1的基础上进行，其包括后续步骤s5。

s5、对植被遮挡目标，进行可见性分析。

s51、对植被遮挡目标进行差分位置分析。

s511、计算植被遮挡区域和植被遮挡目标相邻点之间的三维坐标差值，获得差分位置矩阵。

为了能够探寻激光束在穿透树叶间隙时，所贯穿的树木冠层和隐蔽目标三维坐标的变化并挖掘其位置特性。本发明分别计算树冠遮挡区域和隐蔽目标相邻点之间的三维坐标差值：

dx＝xi-xi-1；

dy＝yi-yi-1；

dz＝zi-zi-1；

式中，其中i为样本中第i个点，xi、yi及zi分别为第i个点的三维坐标，dx、dy及dz分别为三维坐标的的差值，a为差值幅度。随之，本发明利用目标点云的差分位置矩阵特征值大小和特征向量方向的变化，分析遮挡目标对隐蔽目标的影响。

s512、将得到的差分位置矩阵求特征值，得到复数矩阵，分别用quiver，feather，compass，angle函数，得到植被遮挡区域和植被遮挡目标的差分位置矩阵特征值可视化图形。

图5为截取的部分植被遮挡区域的位置差分特性和植被遮挡目标的位置差分特性。横坐标代表点的三维坐标差分个数，纵坐标代表点云的差分幅度a。当植被遮挡区域点与点之间距离变大时(图中间区域)，植被遮挡目标的点间距很小，说明捕获的植被遮挡目标的点的密度更大，对于其特性的分析更加具有优势。可以看出植被遮挡区域在位置分布上更加均匀分布，而植被遮挡目标在分布上具有很大的不确定性，主要是受植被郁闭度的影响。将得到的差值矩阵求特征值，得到复数矩阵，分别用quiver，feather，compass，angle函数，得到植被遮挡区域和植被遮挡目标的特征值特性，如图6和图7。

s513、对三维坐标差值做二次幂变换和log函数变换，得到植被遮挡区域和植被遮挡目标的二次幂矩阵特征值和特征向量的可视化图形。

为了进一步可视化分析遮挡目标和隐蔽目标差分位置的相对变化。本发明对坐标差值做二次幂变换和log函数变换，得到图8。其中eval代表特征值样本，evec代表特征向量样本。用颜色标记差分位置特征值的大小和特征向量的变化的方向幅度。

从图8可以看出，当遮挡目标特征值和特征向量变化幅度很大时(图8(a)样本1和样本2)时，对应的隐蔽目标特征值变化幅度很小，特征向量的变化更加连续(图8(b):evec样本1和evec样本2整体连贯性更强)。这说明当树木冠层的叶簇稀疏的频繁变化有利于隐蔽目标点云的生成，被叶簇遮挡的隐蔽目标信息可能更加完善。当遮挡目标特征向量变化幅度较小时(图8(a)样本3)时，对应的隐蔽目标特征向量方向更加具有不连贯性(图8(b):evec样本3中间部分明显缺失)。

s52、对植被遮挡目标进行回波密度分析。

s521、获取植被遮挡区域与植被遮挡目标在x、y、z上的点云密度值，并进行归一化。

实验统计了多个树冠层中植被遮挡区域和植被遮挡目标在x,y,z上的密度特性如图9和10。其中颜色的深浅和线条密集程度反映了密度的大小。线条越密集、颜色越深，代表密度越大。从中可以看出不管在x,y还是z轴，如果植被遮挡区域的密度变大，相应的隐蔽目标会变得稀疏。相反，如果植被遮挡区域变稀疏，那么隐蔽目标的密度会有所增加。

s522、为对植被遮挡目标进行可见性分析，本发明为植被遮挡区域引入ndi指数。

其中，n一次回波代表一次回波的距离范围均值，n末次回波代表末次回波的距离范围均值。

s523、以归一化的点云密度值为横坐标，以ndi指数为纵坐标，获得ndi指数随点云密度变化的示意图，以此达到可见性分析的需要，其结果如图11所示。

由图11中曲线可知，ndi指数随着植被点云密度的增大而减少，特别是当植被冠层点云密度达到一定程度时，ndi指数骤降。因此，本发明中遮挡可见性分析方法适用于高层次少间隙的叶簇穿透模型。由图12选取的三维点云来看，当树叶茂盛时，对于后面的建筑物，该树叶冠层将变得几乎不可穿透。

s53、对植被遮挡目标进行回波强度分析。

点云场景中，位于不同位置的点，会因为tls的相对位置与距离范围、点的质地属性等人为因素和自然因素而后向散射出不同强度值的激光束。激光的能量后向散射传输公式为：

其中pt是发射功率，ρ是反射率，r是有效散射半径。这个方程式描述相关功率密度的高低。随着激光束的传播。获得散射体截获的功率密度必须乘以有效散面积as。因为激光波长总是小于散射元素的大小,有效面积只是激光束与目标接触的投影面积。一些能量被散射体接收而被吸收，其余的是向不同方向重新辐射。

s531、以点云强度值为横坐标，以点个数为纵坐标，同时，以不同颜色区分不同强度值，获得植被遮挡区域与植被遮挡目标的强度特性图。

树冠层中遮挡目标和隐蔽目标的强度特性如图13，其中横坐标代表点云强度值，纵坐标为该强度值下点云个数。颜色变化的深浅也代表强度值变化的趋势，如红黄色代表强度高，绿色代表强度适中，蓝色代表强度很低。从图13中可以看出遮挡目标强度均值明显高于隐蔽目标。遮挡目标点的众数落在[3200,3300]区间，而隐蔽目标点的众数落在[2400,2500]区间。确定强度阈值，可以划分具有不同质地属性的目标。

s532、以xy平面植被遮挡区域的差分位置距离为横坐标，以植被遮挡目标的点云强度为纵坐标，获得二者的相关性二维特性图。

s533、分别以植被遮挡区域在x轴与y轴方向上的差分位置距离dx和dy为x轴与y轴，以植被遮挡目标的点云强度为纵坐标，建立xy平面上差分位置与点云强度的相关性三维特性图。

由基于pointnet的中间回波植被遮挡区域提取的实验发现，单靠目标的三维坐标区分遮挡区域和非遮挡区域效果并不理想。而且实验发现对于中间回波，其点云坐标x,y和回波强度i相关性较好。为继续探索和验证该结论，本发明选取了一段植被遮挡区域实验样本，如图14。图14(a)描述了遮挡目标xy平面上点云差分位置距离和隐蔽目标点云回波强度的关系。实现发现隐蔽目标回波强度与植被点云差分位置距离成正相关。如图14(b)(c)，横坐标和纵坐标分别代表植被等遮挡目标在x轴和y轴方向上的差分位置距离dx和dy，黄色区域代表隐蔽目标点云强度最强的点，蓝色区域代表隐蔽目标点云强度最弱的点。由图14(d)可知，如果植被点云中点与点之间在xy平面上的距离差值越大，那么激光束穿透的几率越高，由此隐蔽目标产生的回波强度也不会因为遮挡效果而发生巨大衰减。而由图像蓝色区域可知，点与点之间间隔很小，即在树丛密集的情况下，树冠层会很大程度地衰减激光束强度，从而大大影响了隐蔽目标的探测。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。