[发明专利]一种基于点云数据的树木自动分割方法

权利要求书

1.一种基于点云数据的树木自动分割方法，其特征在于，包括以下步骤：

S1、预处理：包括噪声点去除、数据抽稀、空间数据索引的建立以及地面滤波的预处理工作；

S2、空间网格化投影：将步骤S1中得到的点云数据根据离散点的空间密度分布自适应的将空间域划分为若干正方形格网，将三维数据垂直投影到二维平面；

S3、滤除不合格网格：剔除步骤S2中带有少量信息的网格；

S4、基于启发式搜索的树干位置识别：对步骤S3得到的网格进行搜索，得到树干所在位置的网格；

S5、划分单树树冠：利用步骤S4中得到的树干位置作为初始种子，通过改进的分水岭方法划分个体树冠。

2.如权利要求1所述的一种基于点云数据的树木自动分割方法，其特征在于：其还包括步骤S6、验证与矫正：对步骤S5中得到的分割结果进行验证。

3.如权利要求1所述的一种基于点云数据的树木自动分割方法，其特征在于，所述步骤S1包括以下分步骤：

S11、建立KD-tree空间数据索引；

S12、利用RANSAC自动检测和剔除地面点，设置被检测的目标函数模型为平面，同时设置距离参数阈值来决定RANSAC目标函数中的判断为内点的准则，再利用生成的KD-tree索引有效地提取出城区中的平坦地面。

4.如权利要求1所述的一种基于点云数据的树木自动分割方法，其特征在于，所述步骤S2包括以下分步骤：

S21、计算空间格网区域激光扫描点云的数量N，再根据点云数据的最小包围矩阵，矩阵范围为(x_range,y_range)，计算出该区域的面积S，其中x_range和y_range的范围分别是点云数据里的(x_min，x_max)和(y_min，y_max)。由此可以计算出激光点云数据的密度D：

D＝N/S；

S22、根据计算出的数据密度D，以及空间网格所要求的平均点云数目得到单位网格的面积Sgrid，最后获得网格大小GridSize以及空间网格所在的行列数目(Row，Col)；

S23、建立一个行乘以列大小的数组或容器，用于存储每个网格的属性信息；

S24、遍历激光扫描点云最小包围矩阵内的所有点，判断每个点所在的格网位置，对于落在格网边界上的点，规定每个格网包括落在其左边界和下边界的点，不包括其右边界和上边界的点。

5.如权利要求4所述的一种基于点云数据的树木自动分割方法，其特征在于，所述步骤S24中，所述激光点P(x，y)落在格网ImageInfo[r][c]的判别公式为：

r＝(x-x_min)/GridSize

c＝(y-y_min)/GridSize 。

6.如权利要求1所述的一种基于点云数据的树木自动分割方法，其特征在于，所述步骤S3包括以下分步骤：

S31、对于步骤S2得到的网格，如果网格中的最大高度小于设定的某一阈值，则删除该网格中的所有点；

S32、对于步骤S2得到的网格，通过计算网格中的总点数除以网格面积得到网格密度值，如果小于某一阈值，则删除该网格中的所有点。

7.如权利要求6所述的一种基于点云数据的树木自动分割方法，其特征在于，所述步骤S31中的所述网格的高度阈值根据实际树木高度进行设定，所述步骤S32中的所述网格的密度阈值根据激光点云的密度进行设定。

8.如权利要求1所述的一种基于点云数据的树木自动分割方法，其特征在于：所述步骤S4包括以下分步骤：

S41、随机的选择初始网格位置，将其设置为当前网格；

S42、对当前网格的邻域进行搜索，如果某个邻域网格的密度高于当前网格，那么当前网格将往更高密度的邻域网格移动。一次搜索过程将进行到所有的邻域网格密度都小于当前网格时停止，此时当前网格被标记为树干位置网格；

S43、重复步骤S41、S42直到该爬山搜索直到抵达设定的迭代次数。