CN108198244A - 一种苹果叶片点云精简方法及装置 - Google Patents

Abstract

本发明提供一种苹果叶片点云精简方法及装置，其中所述方法包括：利用包围盒法进行快速K‑近邻搜索，建立点云的kd‑tree空间存储结构，设定不同的阈值识别叶片点云边界并提取；通过计算点的特征参数法向量、曲率等，并进行邻域点位置信息计算，区别特征点及非特征点，并进一步对非特征点进行精简处理。本发明能够快速方便地得到边界点云及非边界点云，并进一步得到精简结果，过程中可以根据需要设定不同的K值及多个阈值，所得点云精简结果精确度较高，且计算过程便捷、计算方法合理，并适宜自动化编程实现，有效降低了计算机资源的浪费，可在一定程度上提高作业效率。

Description

一种苹果叶片点云精简方法及装置

技术领域

本发明涉及计算机图形学中植物对象的三维建模和虚拟农业技术领域，更具体地，涉及一种苹果叶片点云精简方法及装置。

背景技术

果树叶片的形态结构以及空间分布直接影响果实的质量和产量，对果树叶片的形态结构进行分析和三维重建能为冠层光照分布计算以及果树整形修剪等提供理论基础。国内外学者针对植物叶片三维重建及可视化开展了大量的工作，随着三维激光扫描设备的生产数量增多和价格的降低，其能够被逐渐应用到植物建模领域中，对植物空间信息的表达精准且快速。

通过三维激光扫描设备可获取到植物叶片大规模的点云数据，点云是在同一空间参考系下表达目标空间分布和目标表面特性的海量点的集合。若将大规模的点云数据直接用于坐标变换、参数提取、颜色映射、曲面重建等植物叶片三维重建的几何操作中将占用大量的计算机资源并且运算效率低下，因此，有必要对获取到的点云数据进行精简。目前，针对苹果叶片的点云精简方法相对较少，较多数的点云精简方法是针对刚性器械元件，少量针对植物叶片的点云精简方法也多倾向于叶片较大、特征简单的植物，如玉米、黄瓜等，通常采用随机采样、统一采样、等距采样、曲率采样和格栅法等方法实现点云数据的精简。

由于苹果叶片形态特征复杂且对细节描述要求较高，上述各单一采样方法可能会存在细节特征丢失严重、空洞问题，不适用于苹果叶片点云数据的精简。

发明内容

为了克服现有点云精简方法存在的细节特征丢失严重、空洞，不适用于苹果叶片的问题，本发明提供一种苹果叶片点云精简方法及装置。

根据本发明的一个方面，提供一种苹果叶片点云精简方法，包括：

S1，使用三维激光扫描设备获取苹果叶片的空间三维信息，并对所述空间三维信息进行去噪处理，获得原始点云；

S2，利用基于包围盒的K-邻域搜索方法计算所述原始点云中的任一点的K-邻域，并通过计算每个点与其K-邻域内各邻近点之间的平均距离确定是否为离群点，将所有离群点从所述原始点云中删除，获得当前点云；

S3，基于最小二乘法计算所述当前点云中的任一点的法向量、K-邻域内的局部平均曲率以及所述当前点云的全局平均曲率和曲率方差；

S4，通过建立k-d树空间拓扑结构并利用所述当前点云中各个点与其K-邻域内各邻近点之间的位置关系，确定当前点云中的边界点；

S5，对于当前点云中的任一非边界点，根据所述全局平均曲率、曲率方差以及预设的点云精简率确定曲率阈值，并根据所述非边界点的K-邻域内的局部平均曲率与所述曲率阈值间的大小关系以及所述非边界点与其K-邻域内所有邻近点的法向量的点积和的正负情况，将所述非边界点划分为特征点及非特征点；

S6，根据所述预设的点云精简率，将所有非特征点所组成的点云划分为多个边长等长的子立方体，计算每个所述子立方体的重心点，并将所述边界点、特征点及重心点存储为最终的精简点云结果。

其中，所述步骤S2进一步包括：

S21，建立所述原始点云的最大包围盒；

S22，将所述最大包围盒划分为边长等长的多个子立方体，从所述原始点云中任意选取一个点作为当前点，根据当前点的坐标值，计算当前点所在子立方体的索引并进行存储；

S23，在当前点所在子立方体及相邻的26个子立方体中搜索距离当前点最近的前K个点作为当前点的K-邻域，并在搜索过程中，通过创建两个向量分别存储当前点的K-邻域内各邻近点的索引以及当前点与其K-邻域内各邻近点间的欧氏距离；

S24，根据当前点与其K-邻域内各邻近点间的欧氏距离计算当前点与其K-邻域内各邻近点间的平均距离；

S25，重复步骤S21至S24直至遍历完所述原始点云，根据所述原始点云中的各个点与其K-邻域内各邻近点间的平均距离计算全局平均距离和距离方差；

S26，根据所述全局平均距离和距离方差确定距离阈值，将所获得的各个点与其K-邻域内各邻近点间的平均距离与所述距离阈值逐一进行比较，若所述平均距离大于距离阈值，则将所述平均阈值对应的点标记为离群点；

S27，将所有离群点从所述原始点云中删除，获得当前点云。

其中，所述步骤S3进一步包括：

S31，确定曲面拟和方程，并根据最小二乘拟合原理，使所述曲面拟合方程取最小值，获得下式：

S32，分别对上式中的常系数a,b,c求偏导并取0，获得常系数a，b，c的值；

S33，记

分别记偏导为f_x,f_y,f_xy,f_xx,f_yy，则所述当前点云中的任一点P_i的法向量

S34，根据上述各偏导f_x，f_y，f_xy，f_xx及f_yy的值，分别计算所述点P_i的高斯曲率K_Gauss、K-邻域内的局部平均曲率H_i；

S35，重复步骤S31至S34，直至遍历完所述当前点云，并根据当前点云各点的K-邻域内的局部平均曲率计算当前点云的全局平均曲率H_global以及曲率方差。

其中，所述步骤S4进一步包括：

S41，以当前点云中的任一点P_i(x_i,y_i,z_i)作为中心点，搜索所述中心点的K-邻域N_b(P_i)＝{P₁,P₂,…,P_k}，确定该中心点的k-邻域中距离该中心点最近的邻近点P_m(x_m,y_m,z_m)；

S42，以所述最近的邻近点P_m(x_m,y_m,z_m)与中心点P_i(x_i,y_i,z_i)组成的向量为基准向量，其他邻近点P_j(x_j,y_j,z_j)与中心点的向量为计算向量与基准向量间的夹角的余弦cosθ，并设置标识f⁺用于记录cosθ＞0的邻近点个数；

S43，设定标识阈值f_σ，当f⁺＞f_σ时，将该中心点标记为边界点进行提取，反之标记为非边界点。

其中，所述步骤S5进一步包括：

S51，根据所述全局平均曲率H_global和曲率方差确定曲率阈值的范围，在所述范围内基于预设的点云精简率选取曲率阈值的具体值H_α；

S52，对所述当前点云中的任一非边界点，若所述非边界点的K-邻域内的局部平均曲率H_i大于所述曲率阈值H_α，则将所述非边界点标记为保留点；

S53，计算所述非边界点与其K-邻域内所有邻近点的法向量的点积和，若所述点积和为非正数，则将所述非边界点标记为保留点；

S54，取步骤S52和步骤S53所获得的保留点的并集，标记为特征点，所述当前点云的非边界点中剩余的点则作为非特征点。

其中，在所述步骤S4之后并在所述步骤S5之前，还包括：

通过设置不同的K值获取多种边界提取效果，并通过综合分析各K值下的相对提取点效率，确定可获得最佳边界提取效果的K值，并将该K值下所获得的当前点云的边界点进行提取。

其中，所述点P_i的高斯曲率K_Gauss计算公式为：

所述点P_i的K-邻域内的局部平均曲率H_i的计算公式为：

上述各式中，E＝f_x·f_x，F＝f_x·f_y，G＝f_y·f_y，L＝f_xx·n，M＝f_yx·n，N＝f_yy·n，s为当前点云中点的个数。

其中，所述步骤S23还包括：

若在所述当前点所在子立方体及相邻的26个子立方体中搜索到的近邻点的个数未达到K个，则将当前点作为离散点从所述原始点云中删除，并开始遍历下一个点。

根据本发明的另一个方面，提供一种苹果叶片点云精简装置，包括：

去噪模块，用于使用三维激光扫描设备获取苹果叶片的空间三维信息，并对所述空间三维信息进行去噪处理，获得原始点云；

去离群点模块，用于利用基于包围盒的K-邻域搜索方法计算所述原始点云中的任一点的K-邻域，并通过计算每个点与其K-邻域内各邻近点之间的平均距离确定是否为离群点，将所有离群点从所述原始点云中删除，获得当前点云；

曲率计算模块，用于基于最小二乘法计算所述当前点云中的任一点的法向量、K-邻域内的局部平均曲率以及所述当前点云的全局平均曲率和曲率方差；

边界提取模块，用于通过建立k-d树空间拓扑结构并利用所述当前点云中各个点与其K-邻域内各邻近点之间的位置关系，确定当前点云中的边界点；

特征点提取模块，用于对于当前点云中的任一非边界点，根据所述全局平均曲率、曲率方差以及预设的点云精简率确定曲率阈值，并根据所述非边界点的K-邻域内的局部平均曲率与所述曲率阈值间的大小关系以及所述非边界点与其K-邻域内所有邻近点的法向量的点积和的正负情况，将所述非边界点划分为特征点及非特征点；

非特征点精简模块，用于根据所述预设的点云精简率，将所有非特征点所组成的点云划分为多个边长等长的子立方体，计算每个所述子立方体的重心点，并将所述边界点、特征点及重心点存储为最终的精简点云结果。

其中，还包括：

边界提取优化模块，通过设置不同的K值获取多种边界提取效果，并通过综合分析各K值下的相对提取点效率，确定可获得最佳边界提取效果的K值，并将该K值下所获得的当前点云的边界点进行提取。

本发明提出的一种苹果叶片点云精简方法及装置，首先基于K邻域搜索算法对经过去噪处理后的原始点云进行去离群点的处理，然后再进行边界提取，并将非边界点再次划分为特征点与非特征点，能够避免特征信息的流失，减少了空洞现象的发生，可以有效的保留局部曲面变化明显的特征区域，并对非特征点进行基于子立方体重心点保留的精简方法，也进一步地避免了点云的空洞现象的发生，能够为进一步的苹果叶片三维重建研究减少数据优化处理的步骤，可有效提升处理效率。

附图说明

图1为根据本发明一实施例提供的一种苹果叶片点云精简方法的流程示意图；

图2为根据本发明一实施例提供的苹果样本叶片及三维扫描仪获取的样本叶片的三维点云数据；

图3为根据本发明一实施例提供的苹果叶片手动剔除噪点过程图；

图4为根据本发明一实施例提供的去除离散点前后邻域点集距中心点平均距离折线图；

图5为根据本发明一实施例提供的不同K值边界提取效果图；

图6为根据本发明一实施例提供的苹果叶片点云精简方法与其他点云精简方法的效果对比图；

图7为根据本发明一实施例提供的苹果叶片点云精简装置的结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他的实施例，都属于本发明保护的范围。

如图1所示，为本发明一实施例提供的一种苹果叶片点云精简方法的流程示意图，包括：

S1，使用三维激光扫描设备获取苹果叶片的空间三维信息，并对所述空间三维信息进行去噪处理，获得原始点云；

S2，利用基于包围盒的K-邻域搜索方法计算所述原始点云中的任一点的K-邻域，并通过计算每个点与其K-邻域内各邻近点之间的平均距离确定是否为离群点，将所有离群点从所述原始点云中删除，获得当前点云；

S3，基于最小二乘法计算所述当前点云中的任一点的法向量、K-邻域内的局部平均曲率以及所述当前点云的全局平均曲率和曲率方差；

S4，通过建立k-d树空间拓扑结构并利用所述当前点云中各个点与其K-邻域内各邻近点之间的位置关系，确定当前点云中的边界点；

S5，对于当前点云中的任一非边界点，根据所述全局平均曲率、曲率方差以及预设的点云精简率确定曲率阈值，并根据所述非边界点的K-邻域内的局部平均曲率与所述曲率阈值间的大小关系以及所述非边界点与其K-邻域内所有邻近点的法向量的点积和的正负情况，将所述非边界点划分为特征点及非特征点；

S6，根据所述预设的点云精简率，将所有非特征点所组成的点云划分为多个边长等长的子立方体，计算每个所述子立方体的重心点，并将所述边界点、特征点及重心点存储为最终的精简点云结果。

具体地，S1中，首先使用三维激光扫描设备，如ArtecSpider，获取苹果叶片的空间三维信息即点云，如图2中所示，图2左边的图为苹果叶片的原始图像，图2右边的图为苹果叶片的点云图。由于扫描设备的限制和随机误差的存在，获取的叶片点云数据在叶片边缘和叶片狭缝处不可避免的会有噪点产生，为了最大程度地还原真实苹果叶片的三维形态，需对叶片三维点云进行去噪处理，可以应用平滑滤波算法如高斯滤波对点云进行去噪处理，在进行滤波处理后叶片边缘仍有噪点，再采用手动交互方法剔除边缘噪点，最后应用孤立点排异法去除其余部分。如图3所示，为手动剔除噪点的过程图，其中，图3的左图为苹果叶片的点云图，中图为噪声点，右图为去除噪声点后的点云图，为方便描述，将去除噪声点后的点云图作为原始点云图。记原始点云P＝{P_i|i＝1,2,3,…,N}，其中，P_i(x_i,y_i,z_i)为苹果叶片原始点云中的任一点，N为所述原始点云中点的总数。

S2，采用基于K-邻域搜索的平均距离计算的方法，确定离散点并去除。对经过去噪处理后的所述原始点云中的任一点均执行以下步骤：利用K-邻域搜索方法计算当前点P_i(x_i,y_i,z_i)的K-邻域，记为N_b(P_i)＝{P₁,P₂,…,P_k},计算当前点P_i与其K-邻域内各邻近点之间的距离，并取平均值D(P_i)作为平均距离，预先设定的距离阈值为D_σ，若D(P_i)＞D_σ，则标记P_i为离群点。然后，将所有离群点从所述经过去噪处理后的所述原始点云中删除。

S3，基于最小二乘法计算所述当前点云中的任一点的法向量、K-邻域内的局部平均曲率以及所述当前点云的全局平均曲率和曲率方差；所述的最佳曲率计算方法为最小二乘拟合法，具有很好的稳定性且计算简单。

S4中通过建立散乱点云的kd-tree空间拓扑结构，利用点与其邻域点位置关系信息，确定边界点并提取。主要步骤为：基于步骤S2的快速近邻搜索方法，重新存储K-邻域信息，以距离中心点最近的点与其形成的向量为基准向量，，计算其他邻域点与中心点的向量与基准向量之间的夹角θ，设置标识f⁺记录cosθ＞0的邻域点个数，即以向量角作为角的临界；设定标识阈值f_σ，当f⁺＞f_σ，说明向量角中的点个数小于标识阈值，可以认为中心点P_i(x_i,y_i,z_i)的邻域点分布相对集中于中心点的单侧，故标记为边界点提取，反之标记为非边界点保留。本步骤中通过设定的角度阈值以及标识阈值判断是否为边界点，提高了算法的包容性，减小了小部分点的干扰作用导致的结果偏差。

S5是对非边界点进行精简处理，为了更好的保留叶片的几何形状和特征，在曲率变化较明显的部分保留较多的点，此部分为特征区域；在曲率变化较小的部分，即平坦区域应删除较多的点。特征点的识别主要利用了两次平行判断，其一为根据所述全局平均曲率H_global和曲率方差σ_H确定曲率阈值的范围为(H_global±σ_H)，在所述范围内根据预设的点云精简率选取曲率阈值的具体值H_α，并将其与当前非边界点的K邻域内的局部平均曲率H_i作比较，若H_α<H_i，说明此处为曲面变化相对明显区域，为特征点，反之为非特征点。其二为计算所述非边界点与其K-邻域内所有邻近点的法向量的点积和若满足点积关系说明局部曲面变化明显，则当前非边界为特征点；取所述两次判断所获得的特征点的并集，标记为强制保留，则可将所述当前点云的非边界点中剩余的点则作为非特征点，标记为待删除。

S6，对非特征点进行精简处理，首先提取所有非特征点，根据精简率及特征点的数量，确定子立方体的边长，将非边界点云划分为等大的子立方体，计算每个子立方体的重心点，作为每簇立方体内点集的代表点进行保留，子立方体中的其他点则删除。

最后将边界点、特征点及每簇点集的重心点存储为最终的精简点云结果。

本发明提出的一种苹果叶片点云精简方法，首先基于K邻域搜索算法对经过去噪处理后的原始点云进行去离群点的处理，然后再进行边界提取，并将非边界点再次划分为特征点与非特征点，能够避免特征信息的流失，减少了空洞现象的发生，可以有效的保留局部曲面变化明显的特征区域，并对非特征点进行基于子立方体重心点保留的精简方法，也进一步地避免了点云的空洞现象的发生，能够为进一步的苹果叶片三维重建研究减少数据优化处理的步骤，可有效提升处理效率。

基于上述实施例，所述步骤S2进一步包括：

S21，建立所述原始点云的最大包围盒：设所述原始点云为P＝{P_i|i＝1,2,3,…,N}，其中，P_i(x_i,y_i,z_i)为所述原始点云中的任一点，N为所述原始点云中点的总数，其在X,Y,Z坐标轴方向上的最大值和最小值分别记为x_max,x_min,y_max,y_min,z_max,z_min，则所述原始点云的最大包围盒可表示为：

S22，将所述最大包围盒划分为边长等长的多个子立方体，其中，所述子立方体的边长为l，从所述原始点云中任意选取一个点作为当前点，根据当前点P_i的坐标值，计算当前点P_i所在子空间的索引index(P_i)＝(x,y,z)，其中，索引计算公式为：

S23，在当前点所在子立方体及相邻的26个子立方体中搜索距离当前点最近的前K个点作为当前点的K-邻域，记为N_b(P_i)＝{P₁,P₂,…,P_k}并在搜索过程中，通过创建两个向量pointIdxNKNSearch和pointNKNdist分别存储当前点的K-邻域内各邻近点的索引以及当前点与其K-邻域内各邻近点间的欧氏距离，为后续的邻域点信息调用等降低了计算机开销。

S24，根据当前点与其K-邻域内各邻近点间的欧氏距离计算当前点与其K-邻域内各邻近点间的平均距离；

具体地，首先计算当前点P_i与其K-邻域内各邻近点间的欧式距离dist(P_i)，并取平均值D(P_i)作为平均距离，其中，

S25，重复步骤S21至S24直至遍历完所述原始点云，根据所述原始点云中的各个点与其K-邻域内各邻近点间的平均距离D(P_i)计算全局平均距离D_global和距离方差；

S26，根据所述全局平均距离D_global和距离方差σ_D确定距离阈值D_σ，即距离阈值的取值范围为(D_global±σ_D)，可以根据需要在所述取值范围(D_global±σ_D)内进行设定；将所获得的各个点与其K-邻域内各邻近点间的平均距离D(P_i)与所述距离阈值D_σ逐一进行比较，若所述平均距离大于距离阈值，则将所述平均阈值对应的点P_i标记为离群点；

S27，将所有离群点从所述原始点云中删除，获得当前点云。

其中，在步骤S23中，若在所述当前点所在子立方体及相邻的26个子立方体中搜索到的近邻点的个数未达到K个，则将当前点作为离散点从所述原始点云中删除，并开始遍历下一个点。

可以根据需要确定子包围盒的边长l、K-近邻搜索的K值、距离阈值等多种参数，以实现对苹果叶片进行不同程度的精简。离散点去除结果分析的参数为邻域点距中心点的平均距离，分析结果如图4所示，图4(a)为去除离散点前邻域点集距中心点平均距离折线图，图4(b)为去除离散点前邻域点集距中心点平均距离折线图，从图4(b)中可以看出，折线图相对平稳，说明大量的离散点被去除。

基于上述实施例，所述步骤S3进一步包括：

S31，确定曲面拟和方程，并根据最小二乘拟合原理，使所述曲面拟合方程取最小值，分别如下所示：

S32，分别对(6)式中的常系数a,b,c求偏导并取0得到方程组(7)，可计算出常系数a，b，c的值：

S33，将公式(5)转换为

分别记偏导为f_x,f_y,f_xy,f_xx,f_yy，则所述当前点云中的任一点P_i的单位法向量可近似为

S34，根据上述各偏导f_x，f_y，f_xy，f_xx及f_yy的值，分别计算所述点P_i的高斯曲率K_Gauss、K-邻域内的局部平均曲率H_i，其中，

S35，重复步骤S31至S34，直至遍历完所述当前点云，并根据当前点云各点的K-邻域内的局部平均曲率计算当前点云的全局平均曲率H_global以及曲率方差σ_H，其中，

公式(10)至(12)中，E＝f_x·f_x，F＝f_x·f_y，G＝f_y·f_y，L＝f_xx·n，M＝f_yx·n，N＝f_yy·n，s为当前点云中点的个数。

基于上述实施例，所述步骤S4进一步包括：

S41，以当前点云中的任一点P_i(x_i,y_i,z_i)作为中心点，搜索所述中心点的K-邻域N_b(P_i)＝{P₁,P₂,…,P_k}，确定该中心点的k-邻域中距离该中心点最近的邻近点P_m(x_m,y_m,z_m)；

S42，以所述最近的邻近点P_m(x_m,y_m,z_m)与中心点P_i(x_i,y_i,z_i)组成的向量为基准向量，其他邻近点P_j(x_j,y_j,z_j)与中心点的向量为计算向量与基准向量间的夹角的余弦并设置标识f⁺用于记录cosθ＞0的邻近点个数，还可设置设置标识f^-，用于记录cosθ＜0的邻近点个数；

S43，设定标识阈值f_σ，当f⁺＞f_σ时，说明向量角中的点个数小于标识阈值，可以认为中心点P_i(x_i,y_i,z_i)的邻域点分布相对集中于中心点的单侧，故标记为边界点进行提取，反之标记为非边界点进行保留。其中，可以通过设定一个距离半径r，在以r为半径的球体内找其k邻域内邻近点的个数，标识阈值则取该个数的一半。

基于上述实施例，所述步骤S5进一步包括：

S51，根据所述全局平均曲率H_global和曲率方差确定曲率阈值的范围，在所述范围内基于预设的点云精简率选取曲率阈值的具体值H_α；

根据所述全局平均曲率H_global和曲率方差σ_H确定曲率阈值的范围为(H_global±σ_H)，在所述范围(H_global±σ_H)内根据预设的点云精简率选取曲率阈值的具体值H_α

S52，对所述当前点云中的任一非边界点，若所述非边界点的K-邻域内的局部平均曲率H_i大于所述曲率阈值H_α，则将所述非边界点标记为保留点；

若H_α<H_i，说明此处为曲面变化相对明显区域，为特征点，反之为非特征点。

S53，计算所述非边界点与其K-邻域内所有邻近点的法向量的点积和，若所述点积和为非正数，则将所述非边界点标记为保留点；

假设非边界点P_i与其邻域点法向量分别为n_p，n_q，若满足点积关系说明局部曲面变化明显，则点P_i为特征点；

S5i4，取步骤S52和步骤S53所获得的保留点的并集，标记为特征点，所述当前点云的非边界点中剩余的点则作为非特征点，可以进行删除。

基于上述实施例，在所述步骤S4之后并在所述步骤S5之前，还包括：

通过设置不同的K值获取多种边界提取效果，并通过综合分析各K值下的相对提取点效率，确定可获得最佳边界提取效果的K值，并将该K值下所获得的当前点云的边界点进行提取。

具体地，可以通过设置不同的K值，获取多种边界提取效果，综合分析相对提取点效率等参数，最终确定边界点云。不同K值的提取效果如图5所示，图5中(a)、(b)、(c)、(d)和(e)分别表示K＝20,K＝24,K＝28,K＝32以及K＝36时边界点提取的效果图，边界点提取的相关参数分析如表1中所示。

通过综合分析，优选K＝32时的边界提取结果作为最终结果，此时相对运算效率较高且具有良好效果。

表1不同K值的提取参数表

最后，将本发明实施例所提供的苹果叶片点云精简方法，与传统的曲率采样法、随机采样法、统一采样法和格栅法进行进行精简效果的对比，如图6所示，为所述各方法的精简效果图，其中，(a)为本发明提出的方法的点云精简结果，(b)为基于Geomagic Studio软件的曲率采样法点云精简结果，(c)为随机采样法点云精简结果，(d)为统一采样法点云精简结果，(e)为格栅采样法点云精简结果，各方法的对比参数分析如表2所示。

表2各方法的对比参数分析

如图7所示，为本发明另一实施例提供的一种苹果叶片点云精简装置的结构示意图，包括：去噪模块71、去离群点模块72、曲率计算模块73、边界提取模块74、特征点提取模块75和非特征点精简模块76，其中，

去噪模块71，用于使用三维激光扫描设备获取苹果叶片的空间三维信息，并对所述空间三维信息进行去噪处理，获得原始点云；

去离群点模块72，用于利用基于包围盒的K-邻域搜索方法计算所述原始点云中的任一点的K-邻域，并通过计算每个点与其K-邻域内各邻近点之间的平均距离确定是否为离群点，将所有离群点从所述原始点云中删除，获得当前点云；

曲率计算模块73，用于基于最小二乘法计算所述当前点云中的任一点的法向量、K-邻域内的局部平均曲率以及所述当前点云的全局平均曲率和曲率方差；

边界提取模块74，用于通过建立k-d树空间拓扑结构并利用所述当前点云中各个点与其K-邻域内各邻近点之间的位置关系，确定当前点云中的边界点；

特征点提取模块75，用于对于当前点云中的任一非边界点，根据所述全局平均曲率、曲率方差以及预设的点云精简率确定曲率阈值，并根据所述非边界点的K-邻域内的局部平均曲率与所述曲率阈值间的大小关系以及所述非边界点与其K-邻域内所有邻近点的法向量的点积和的正负情况，将所述非边界点划分为特征点及非特征点；

非特征点精简模块76，用于根据所述预设的点云精简率，将所有非特征点所组成的点云划分为多个边长等长的子立方体，计算每个所述子立方体的重心点，并将所述边界点、特征点及重心点存储为最终的精简点云结果。

具体地，去噪模块71用于使用三维激光扫描设备(如Artec Spider)获取苹果叶片的空间三维信息即点云，由于扫描设备的限制和随机误差的存在，获取的叶片点云数据在叶片边缘和叶片狭缝处不可避免的会有噪点产生，为了最大程度地还原真实苹果叶片的三维形态，需对叶片三维点云进行去噪处理，可以应用平滑滤波算法如高斯滤波对点云进行去噪处理，在进行滤波处理后叶片边缘仍有噪点，再采用手动交互方法剔除边缘噪点，最后应用孤立点排异法去除其余部分。记原始点云P＝{P_i|i＝1,2,3,…,N}，其中，P_i(x_i,y_i,z_i)为苹果叶片原始点云中的任一点，N为所述原始点云中点的总数。

去离群点模块72采用基于K-邻域搜索的平均距离计算的方法，确定离散点并去除。对经过去噪处理后的所述原始点云中的任一点均执行以下步骤：利用K-邻域搜索方法计算当前点P_i(x_i,y_i,z_i)的K-邻域，记为N_b(P_i)＝{P₁,P₂,…,P_k},计算当前点P_i与其K-邻域内各邻近点之间的距离，并取平均值D(P_i)作为平均距离，预先设定的距离阈值为D_σ，若D(P_i)＞D_σ，则标记P_i为离群点。然后，将所有离群点从所述经过去噪处理后的所述原始点云中删除。

曲率计算模块73基于最小二乘法计算所述当前点云中的任一点的法向量、K-邻域内的局部平均曲率以及所述当前点云的全局平均曲率和曲率方差；所述的最佳曲率计算方法为最小二乘拟合法，具有很好的稳定性且计算简单。

边界提取模块74通过建立散乱点云的kd-tree空间拓扑结构，利用点与其邻域点位置关系信息，确定边界点并提取。主要步骤为：基于步骤S2的快速近邻搜索方法，重新存储K-邻域信息，以距离中心点最近的点与其形成的向量为基准向量，计算其他邻域点与中心点的向量与基准向量之间的夹角θ，设置标识f⁺记录cosθ＞0的邻域点个数，即以向量角作为角的临界；设定标识阈值f_σ，当f⁺＞f_σ，说明向量角中的点个数小于标识阈值，可以认为中心点P_i(x_i,y_i,z_i)的邻域点分布相对集中于中心点的单侧，故标记为边界点提取，反之标记为非边界点保留。本步骤中通过设定的角度阈值以及标识阈值判断是否为边界点，提高了算法的包容性，减小了小部分点的干扰作用导致的结果偏差。

特征点提取模块75对非边界点进行精简处理，为了更好的保留叶片的几何形状和特征，在曲率变化较明显的部分保留较多的点，此部分为特征区域；在曲率变化较小的部分，即平坦区域应删除较多的点。特征点的识别主要利用了两次平行判断，其一为根据所述全局平均曲率H_global和曲率方差σ_H确定曲率阈值的范围为(H_global±σ_H)，在所述范围内根据预设的点云精简率选取曲率阈值的具体值H_α，并将其与当前非边界点的K邻域内的局部平均曲率H_i作比较，若H_α<H_i，说明此处为曲面变化相对明显区域，为特征点，反之为非特征点。其二为计算所述非边界点与其K-邻域内所有邻近点的法向量的点积和若满足点积关系说明局部曲面变化明显，则当前非边界为特征点；取所述两次判断所获得的特征点的并集，则可将所述当前点云的非边界点中剩余的点则作为非特征点，非特征点就可以删除。

非特征点精简模块76对非特征点进行精简处理，首先提取所有非特征点，根据精简率及特征点的数量，确定子立方体的边长，将非边界点云划分为等大的子立方体，计算每个子立方体的重心点，作为每簇立方体内点集的代表点进行保留，子立方体中的其他点则删除。最后将边界点、特征点及每簇点集的重心点存储为最终的精简点云结果。

本发明提出的一种苹果叶片点云精简装置，首先基于K邻域搜索算法对经过去噪处理后的原始点云进行去离群点的处理，然后再进行边界提取，并将非边界点再次划分为特征点与非特征点，能够避免特征信息的流失，减少了空洞现象的发生，可以有效的保留局部曲面变化明显的特征区域，并对非特征点进行基于子立方体重心点保留的精简方法，也进一步地避免了点云的空洞现象的发生，能够为进一步的苹果叶片三维重建研究减少数据优化处理的步骤，可有效提升处理效率。

基于上述实施例，还包括：

边界提取优化模块，通过设置不同的K值获取多种边界提取效果，并通过综合分析各K值下的相对提取点效率，确定可获得最佳边界提取效果的K值，并将该K值下所获得的当前点云的边界点进行提取。

具体地，所述苹果叶片点云精简装置还包括边界提取优化模块，可以通过设置不同的K值，获取多种边界提取效果，综合分析相对提取点效率等参数，最终确定边界点云。

最后，本发明的方法仅为较佳的实施方案，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种苹果叶片点云精简方法，其特征在于，包括：

S1，使用三维激光扫描设备获取苹果叶片的空间三维信息，并对所述空间三维信息进行去噪处理，获得原始点云；

S2，利用基于包围盒的K-邻域搜索方法计算所述原始点云中的任一点的K-邻域，并通过计算每个点与其K-邻域内各邻近点之间的平均距离确定是否为离群点，将所有离群点从所述原始点云中删除，获得当前点云；

S3，基于最小二乘法计算所述当前点云中的任一点的法向量、K-邻域内的局部平均曲率以及所述当前点云的全局平均曲率和曲率方差；

S4，通过建立k-d树空间拓扑结构并利用所述当前点云中各个点与其K-邻域内各邻近点之间的位置关系，确定当前点云中的边界点；

S5，对于当前点云中的任一非边界点，根据所述全局平均曲率、曲率方差以及预设的点云精简率确定曲率阈值，并根据所述非边界点的K-邻域内的局部平均曲率与所述曲率阈值间的大小关系以及所述非边界点与其K-邻域内所有邻近点的法向量的点积和的正负情况，将所述非边界点划分为特征点及非特征点；

S6，根据所述预设的点云精简率，将所有非特征点所组成的点云划分为多个边长等长的子立方体，计算每个所述子立方体的重心点，并将所述边界点、特征点及重心点存储为最终的精简点云结果。

2.根据权利要求1所述的方法，其特征在于，所述步骤S2进一步包括：

S21，建立所述原始点云的最大包围盒；

S22，将所述最大包围盒划分为边长等长的多个子立方体，从所述原始点云中任意选取一个点作为当前点，根据当前点的坐标值，计算当前点所在子立方体的索引并进行存储；

S23，在当前点所在子立方体及相邻的26个子立方体中搜索距离当前点最近的前K个点作为当前点的K-邻域，并在搜索过程中，通过创建两个向量分别存储当前点的K-邻域内各邻近点的索引以及当前点与其K-邻域内各邻近点间的欧氏距离；

S24，根据当前点与其K-邻域内各邻近点间的欧氏距离计算当前点与其K-邻域内各邻近点间的平均距离；

S25，重复步骤S21至S24直至遍历完所述原始点云，根据所述原始点云中的各个点与其K-邻域内各邻近点间的平均距离计算全局平均距离和距离方差；

S26，根据所述全局平均距离和距离方差确定距离阈值，将所获得的各个点与其K-邻域内各邻近点间的平均距离与所述距离阈值逐一进行比较，若所述平均距离大于距离阈值，则将所述平均阈值对应的点标记为离群点；

S27，将所有离群点从所述原始点云中删除，获得当前点云。

3.根据权利要求2所述的方法，其特征在于，所述步骤S3进一步包括：

S31，确定曲面拟和方程，并根据最小二乘拟合原理，使所述曲面拟合方程取最小值，获得下式：

S32，分别对上式中的常系数a,b,c求偏导并取0，获得常系数a，b，c的值；

S33，记

分别记偏导为f_x,f_y,f_xy,f_xx,f_yy，则所述当前点云中的任一点P_i的法向量

S34，根据上述各偏导f_x，f_y，f_xy，f_xx及f_yy的值，分别计算所述点P_i的高斯曲率K_Gauss、K-邻域内的局部平均曲率H_i；

S35，重复步骤S31至S34，直至遍历完所述当前点云，并根据当前点云各点的K-邻域内的局部平均曲率计算当前点云的全局平均曲率H_global以及曲率方差。

4.根据权利要求3所述的方法，其特征在于，所述步骤S4进一步包括：

S41，以当前点云中的任一点P_i(x_i,y_i,z_i)作为中心点，搜索所述中心点的K-邻域N_b(P_i)＝{P₁,P₂,···,P_k}，确定该中心点的k-邻域中距离该中心点最近的邻近点P_m(x_m,y_m,z_m)；

S42，以所述最近的邻近点P_m(x_m,y_m,z_m)与中心点P_i(x_i,y_i,z_i)组成的向量为基准向量，其他邻近点P_j(x_j,y_j,z_j)与中心点的向量为计算向量与基准向量间的夹角的余弦cosθ，并设置标识f⁺用于记录cosθ＞0的邻近点个数；

S43，设定标识阈值f_σ，当f⁺＞f_σ时，将该中心点标记为边界点进行提取，反之标记为非边界点。

5.根据权利要求4所述的方法，其特征在于，所述步骤S5进一步包括：

S51，根据所述全局平均曲率H_global和曲率方差确定曲率阈值的范围，在所述范围内基于预设的点云精简率选取曲率阈值的具体值H_α；

S52，对所述当前点云中的任一非边界点，若所述非边界点的K-邻域内的局部平均曲率H_i大于所述曲率阈值H_α，则将所述非边界点标记为保留点；

S53，计算所述非边界点与其K-邻域内所有邻近点的法向量的点积和，若所述点积和为非正数，则将所述非边界点标记为保留点；

S54，取步骤S52和步骤S53所获得的保留点的并集，标记为特征点，所述当前点云的非边界点中剩余的点则作为非特征点。

6.根据权利要求1-5任一所述的方法，其特征在于，在所述步骤S4之后并在所述步骤S5之前，还包括：

通过设置不同的K值获取多种边界提取效果，并通过综合分析各K值下的相对提取点效率，确定可获得最佳边界提取效果的K值，并将该K值下所获得的当前点云的边界点进行提取。

7.根据权利要求3所述的方法，其特征在于，所述点P_i的高斯曲率K_Gauss计算公式为：

所述点P_i的K-邻域内的局部平均曲率H_i的计算公式为：

上述各式中，E＝f_x·f_x，F＝f_x·f_y，G＝f_y·f_y，L＝f_xx·n，M＝f_yx·n，N＝f_yy·n，s为当前点云中点的个数。

8.根据权利要求2所述的方法，其特征在于，所述步骤S23还包括：

若在所述当前点所在子立方体及相邻的26个子立方体中搜索到的近邻点的个数未达到K个，则将当前点作为离散点从所述原始点云中删除，并开始遍历下一个点。

9.一种苹果叶片点云精简装置，其特征在于，包括：

去噪模块，用于使用三维激光扫描设备获取苹果叶片的空间三维信息，并对所述空间三维信息进行去噪处理，获得原始点云；

去离群点模块，用于利用基于包围盒的K-邻域搜索方法计算所述原始点云中的任一点的K-邻域，并通过计算每个点与其K-邻域内各邻近点之间的平均距离确定是否为离群点，将所有离群点从所述原始点云中删除，获得当前点云；

曲率计算模块，用于基于最小二乘法计算所述当前点云中的任一点的法向量、K-邻域内的局部平均曲率以及所述当前点云的全局平均曲率和曲率方差；

边界提取模块，用于通过建立k-d树空间拓扑结构并利用所述当前点云中各个点与其K-邻域内各邻近点之间的位置关系，确定当前点云中的边界点；

特征点提取模块，用于对于当前点云中的任一非边界点，根据所述全局平均曲率、曲率方差以及预设的点云精简率确定曲率阈值，并根据所述非边界点的K-邻域内的局部平均曲率与所述曲率阈值间的大小关系以及所述非边界点与其K-邻域内所有邻近点的法向量的点积和的正负情况，将所述非边界点划分为特征点及非特征点；

非特征点精简模块，用于根据所述预设的点云精简率，将所有非特征点所组成的点云划分为多个边长等长的子立方体，计算每个所述子立方体的重心点，并将所述边界点、特征点及重心点存储为最终的精简点云结果。

10.根据权利要求9所述的装置，其特征在于，还包括：

边界提取优化模块，通过设置不同的K值获取多种边界提取效果，并通过综合分析各K值下的相对提取点效率，确定可获得最佳边界提取效果的K值，并将该K值下所获得的当前点云的边界点进行提取。