[发明专利]基于多层级树木结构分析的树木分割方法

摘要

基于多层级树木结构分析的树木分割方法，包括如下步骤：步骤1，构建树冠高度模型；步骤2，分割上层树木；步骤3，分析分层树木及分割下层树木。本发明针对现有森林遥感技术无法分析上层树木的情况，引入树木结构的分析来改善了下层被遮挡树木提取困难的情况，同时在分割上下层树木时采用的分割方法能够获得较好的结果。

主权项

1.基于多层级树木结构分析的树木分割方法，包括如下步骤：步骤1，构建树冠高度模型；在提取树木之前，需要对LiDAR点云数据进行处理，构建出树冠高度模型并在该模型上提取出树木；具体步骤如下：11)对LiDAR点云数据进行高度归一化操作；初始的LiDAR点云数据可以被分为地面点和非地面点，进行高度归一化操作后，将会改变点云数据的高度值，使地面点的高度值为0，非地面点也会随着地面点而改变；12)构建局部CHM；使用高度归一化后的点云数据来构建不同高度的局部CHM，局部CHM是过滤掉特定高度值以下的点云数据生成的CHM，每隔5米生成相应的局部CHM；13)将多个CHM聚合为一个CHM，将所有不同高度值的局部CHM进行叠加，在每个像素位置上取最高值，最终获得CHM；步骤2，分割上层树木；为了缩小分层分析的范围，先对CHM进行预分割，提取出上层树木的轮廓信息；分割算法的主要思路是利用CHM的拓扑结构，将CHM划分为小节点并分析节点与节点之间的关系；利用权重值的关系将属于一棵树的节点划分到一起；分割算法具体步骤如下：21)将CHM的高度值划分为64位高度级；将相邻且相同高度级的像素点汇聚成块，聚集过程采用8邻域；22)利用邻域聚合将这个块视为节点构建出块有向图(P‑DAG),节点为像素聚集块，相邻且高度较高的块指向高度较低的块作为边；在块有向图中一个节点可以有多个父节点和子节点；拥有局部最大值的块视为没有父节点的块，这个顶点块也可作为一个层级块的顶点；将相邻的块进行连接，高度级较高的块指向高度级较低的；块可能属于多个层级；构建出P‑DAG后，将层级块作为一个实体节点，研究层级块与层级块之间的关系，构建出层级块有向图(H‑DAG)；层级块是以顶点块为中心，以这个块为起点所有能到达的块聚集而成的；与块有向图相同，层级块有向图也是较高的层级块指向较低的层级块；但是当高度级相同时，通过层级块的像素数量和层级块中顶点块的像素数量来判断层级块之间的指向；23)构建出层级块有向图后，需要计算层级块与层级块之间边的权重值，当层级块属于同一棵树时，层级块之间的权重值就越大；有向图中的权重值由三个部分组成：高度距离(LD)，节点距离(ND)和层级块中心距离(TD)；i.高度距离(SLD(p,q))是两个层级块中的顶点块与共有块之间最小的高度差；当高度差越小时，两个层级块属于同一棵树的可能性就越大；高度距离通过下面的公式来计算，其中，p和q是两个用于分析相互关系的层级块，rp是层级块p中的顶点块，rq是层级块q中的顶点块，Hp是层级块p中所有块的集合，Hq是层级块q中所有块的集合，dh(r,i)是块r和块i之间的高度差，r是rp和rq的其中一个，i是Hp和Hq共有的块；ii.节点距离(SND(p,q))是两个层级块之间层级块的顶点块到共有块的最小块数量；当距离值越小时，两个层级块属于同一棵树的可能性就越大；节点距离通过下面的公式来计算，其中dn(r,i)是P‑DAG中块r和块i之间块的数量；iii.层级块中心距离(STD(p,q))代表着两个层级块之间顶点块的水平距离；如果顶点块有多个像素，则用块的中心块来代替；同样，当距离越小时，两个层级块属于同一棵树的可能性就越大；节点距离通过下面的公式来计算，其中C(rp)和C(rq)是层级块中顶点块的中心坐标；iv.有了上面三个聚合条件后，可以计算出层级块之间边的权重值(W(p,q))，公式如下：W(p,q)＝WLD×SLD(p,q)+WND×SND(p,q)+WTD×STD(p,q)   (4)其中WLD，WND和WTD是权重系数，采用WLD：WND：WTD＝3:3:4的比重；v.计算出权重值后需要过滤层级块有向图中的边，权重值越大代表两个层级块之间联系得越紧密，两个层级块为一颗树木的可能性也就越大；反之亦然，因此我们将会过滤掉权重值小于0.5的边，当指向同一个层级块有多条边时，取权重值最大的边；24)通过层级块之间边的分配，可以知道哪些层级块属于同一棵树，但是在这种情况下无法分配共有块，也就没法获取每棵树木的轮廓；因此，对共有的块按距离来分配，将共有块分配给更近的层级块；分配完共有块后即可获得最终的分割结果；步骤3，分析分层树木及分割下层树木；初步分割后的结果可能代表着单棵树木或者多棵树木，获取到预分割的结果后，通过相应区域的点云数据来分析分割后的层级结构；利用点云数据生成LiDAR点云高度直方图，统计相应高度值的点数量，并且对直方图进行拟合处理；当存在着上下两层树木时，通常会在图中呈现出凹陷的现象；利用直方图中断层的特性，采用滑动窗口法求得断层高度位置；随后将上下两层树木分开处理，在获得上层树木信息的同时分析出下层树木信息；如果检测出了多个断点，则把最高的断点作为断层的位置并忽视其他断层位置；最终获得下层树木的LiDAR数据；通过上述操作，从原始数据中提取出下层树木的LiDAR数据；这时采用步骤2的分割方法对下层LiDAR数据进行分割来获得下层树木的信息。