一种基于激光点云数据的树干材积获取方法技术

本发明专利技术公开提供了一种基于激光点云数据的树干材积获取方法，包括：步骤1，对树干点云垂直分段；步骤2，逐一计算各个垂直分段的体积；步骤3，计算所有垂直分段体积的累加和。本发明专利技术在获取树干断面轮廓的基础上，通过准确计算树干断面的断面积来准确计算树干每段的树干材积，具有较高的精确性，并为准确提取树干生物量、碳储量与树干量材优化等提供基础数据。

【技术实现步骤摘要】
一种基于激光点云数据的树干材积获取方法
本专利技术涉及一种树干参数提取技术，尤其是涉及一种基于激光点云数据的树干材积获取方法。
技术介绍
地面三维激光扫描技术(TerrestrialLaserScanningTechnology，简称TLS)是一种新型的测绘技术，产生于20世纪90年代。TLS可快速准确的获得目标物体表面的点云数据信息，从而可以快速构建目标物体的三维模型。经过近20多年来的发展，三维激光扫描仪已经可以连续快速地对被观测物体进行非接触式测量，其通过获取物体表面至扫描仪的距离和发射强度获取大量物体表面的三维点云数据。使用TLS提取林业相关参数是近10年来的研究热点。树干材积是测树学中的基本测树因子，是树木经济价值的直接体现，也是计算树木生物量与碳储量的基础。传统上将树干简化为一个绕中心轴旋转的旋转体，以此为基础估算树干材积。然而树干形状是不规则的，旋转体的简化结构过于简化了树干的几何形状，在此基础上计算得到的树干材积必定与树干的真实材积存在着差异，也正因如此，如何准确计算树干材积是林业工作者一直关注的问题。
技术实现思路
有鉴于此，本专利技术的目的是针对现有技术的不足，提供一种基于激光点云数据的树干材积获取方法，用以根据树干点云获取精确的树干材积。为达到上述目的，本专利技术采用以下技术方案：一种基于激光点云数据的树干材积获取方法，包括如下步骤：步骤1，对树干点云垂直分段；步骤2，逐一计算各个垂直分段的体积；步骤3，计算所有垂直分段体积的累加和。优选的，所述步骤1，对树干点云垂直分段，包括：以高度为标准，计算树干点云中的最低点plowest，其中zlowest为最低点plowest的高度值，并以所述最低点plowest为基础计算树干的高度；构建平面π0与π1，其方程分别为z＝zlowest与z＝zlowest+thickness，其中thickness是每一个垂直分段的高度；则位于平面π0与π1间的点云形成树干点云的一个垂直分段Q0；继续向上构建平面π2，其方程为z＝zlowest+thickness*2，则位于平面π1与π2间的树干点云形成另一个垂直分段Q1；则，平面πi的方程可表示为z＝zlowest+thickness*i，位于平面πi与πi+1间的树干点云形成点云垂直分段Qi；其中，平面πi与平面πi+1分别称为垂直分段Qi的下平面与上平面；按照上述垂直分段的构建方法直至树梢部位得到若干个树干点云的垂直分段。优选的，所述步骤2，逐一计算各个垂直分段的体积，包括：计算垂直分段在下平面的投影点集；计算垂直分段的断面轮廓点；构建断面轮廓曲线并计算树干断面的断面积；计算垂直分段的体积。优选的，所述计算垂直分段在下平面的投影点集，包括：计算垂直分段Qi在其下平面的投影点集Qi′；将垂直分段Qi中的树干点云投影至垂直分段Qi的下平面πi得到一个平面点集Qi′；计算平面点集Qi′的凸包多边形的质心点作为平面点集Qi′的几何中心点ci。优选的，所述计算垂直分段的断面轮廓点，包括：设点pm为垂直分段Qi中的一个点，pm′表示点pm在垂直分段Qi下平面的投影，则几何中心点ci与投影点pm′构成的向量与X轴的夹角可表示为：根据角度分区的角度值θ，将平面点集Qi′中的点分为个角度分区，每个角度分区Qij中的点可表示为：计算每一个角度分区Qij中点的重心点，所述重心点为所述角度分区中所有点的几何平均值；并将所述重心点作为所述角度分区的轮廓点，由此得到垂直分段Qi的个轮廓点。优选的，所述构建断面轮廓曲线并计算树干断面的断面积，包括：设垂直分段Qi的断面轮廓点集合用Rprofile表示，首先采用非有理三次B样条曲线构建一条插值于点集Rprofile的闭合曲线，所述闭合曲线即为断面轮廓曲线；断面轮廓曲线所围区域的面积即为树干断面的断面积。优选的，所述构建断面轮廓曲线并计算树干断面的断面积，包括：将Rprofile首部的3个点依次加入到Rprofile的尾部同时也将尾部的3个点依次加入到Rprofile的首部，构成插值点R＝{Rk},k＝0,1,...,n；使用弦长参数法计算插值点R的节点参数值令d为总弦长，节点参数值定义为：根据节点参数值使用平均值法计算节点向量U＝{u0,u0,...,um}，有其中，p＝3表示样条曲线的次数，n＝|R|表示插值点的个数，m＝n+p-1；采用全局曲线插值的方法求解非有理三次B样条曲线的控制点，得到系数矩阵为(n+1)×(n+1)的线性方程组：其中，Ni,p(u)表示为求解方程组得到n+1个控制点{Pi}，根据非有理三次B样条曲线的表达式p＝3；可得到一条连续曲线，设置Rprofile的第一个节点R0为首节点，从首节点处开始构建，当曲线环绕一周到达首节点时停止曲线的构建，即可得到完整的闭合曲线，所述的闭合曲线即为断面轮廓曲线；断面轮廓曲线是一条闭合曲线，根据格林公式将该闭合曲线所围区域的面积转换为第一类曲线积分的计算，从而得到树干断面的断面积Si。优选的，所述计算垂直分段的体积，包括：计算垂直分段Qi的体积Vi，将垂直分段近似为一个柱体，断面形状为断面轮廓曲线所围区域的形状，则根据柱体体积公式，垂直分段Qi的体积Vi可表示为Vi＝Si*thickness。优选的，所述步骤3，计算所有垂直分段体积的累加和，包括：所有垂直分段体积的累加和即为树干材积，树干材积的计算公式可表示为V＝thickness*∑Si，V为树干材积。优选的，所述垂直分段的高度值thickness为0.5厘米。本专利技术的有益效果是：传统上使用区分求积法计算伐倒木的树干材积，其建立在将树干简化为树干横断面是圆，树干中心轴为直线的旋转体的基础上，该方法将树干区分成若干等长或不等长的区分段，使各区段的干形更接近于正几何体，再根据近似求积式计算各分段材积，合计各分段材积得到树干材积。计算立木的树干材积的传统方法是采用经验公式，将立木高度与直径的测量值代入材积经验公式得到立木的树干材积。由于材积经验公式建立在伐倒木材积计算的基础上，因此伐倒木与立木的材积计算都建立在树干横断面为圆，树干中心轴为直线的规则旋转体的几何理论基础之上。然而树干形状是不规则的，树干横断面也不是规则的圆形，因此在此简化基础上计算的材积与树干的真实材积存在着差异。本专利技术建立在精确计算树干断面积的基础上，通过精确获取树干断面的轮廓点，进而构建一条闭合的断面轮廓曲线，闭合曲线所围面积即为断面面积，该段断面积计算方式建立在获取树干断面轮廓的基础上，贴近于树干断面的真实轮廓；而且使用的thickness值为0.5厘米，对于一段长度为0.5厘米的树干，其在各处的断面几乎相同，因此本专利技术计算的垂直分段的体积非常接近于垂直分段的真实体积。相对比传统上伐倒木和立木的不同材积计算方法，本专利技术在获取树干点云的基础上，无需区分伐倒木还是立木，均可使用本专利技术的方法计算树干材积。本专利技术在精确计算树干断面断面积的基础上精确计算树干材积，相比于传统的将树干断面简化为圆来计算树干断面的断面积的方法，本专利技术在获取树干断面轮廓的基础上计算树干断面的断面积，具有较高的精确性，从而为准确提取树干生物量、碳储量与树干量材优化等与树干材积的相关研究与应用提供基础数据。本专利技术的其它特征和优点将在随后的说明本文档来自技高网...

【技术保护点】
一种基于激光点云数据的树干材积获取方法，其特征在于：包括如下步骤：步骤1，对树干点云垂直分段；步骤2，逐一计算各个垂直分段的体积；步骤3，计算所有垂直分段体积的累加和。

【技术特征摘要】
1.一种基于激光点云数据的树干材积获取方法，其特征在于：包括如下步骤：步骤1，对树干点云垂直分段；步骤2，逐一计算各个垂直分段的体积；步骤3，计算所有垂直分段体积的累加和。2.根据权利要求1所述的一种基于激光点云数据的树干材积获取方法，其特征在于：所述步骤1，对树干点云垂直分段，包括：以高度为标准，计算树干点云中的最低点plowest，其中zlowest为最低点plowest的高度值，并以所述最低点plowest为基础计算树干的高度；构建平面π0与π1，其方程分别为z＝zlowest与z＝zlowest+thickness，其中thickness是每一个垂直分段的高度值；则位于平面π0与π1间的点云形成树干点云的一个垂直分段Q0；继续向上构建平面π2，其方程为z＝zlowest+thickness*2，则位于平面π1与π2间的树干点云形成另一个垂直分段Q1；则，平面πi的方程可表示为z＝zlowest+thickness*i，位于平面πi与πi+1间的树干点云形成点云垂直分段Qi；其中，平面πi与平面πi+1分别称为垂直分段Qi的下平面与上平面；按照上述垂直分段的构建方法直至树梢部位得到若干个树干点云的垂直分段。3.根据权利要求2所述的一种基于激光点云数据的树干材积获取方法，其特征在于：所述步骤2，逐一计算各个垂直分段的体积，包括：计算垂直分段在下平面的投影点集；计算垂直分段的断面轮廓点；构建断面轮廓曲线并计算树干断面的断面积；计算垂直分段的体积。4.根据权利要求3所述的一种基于激光点云数据的树干材积获取方法，其特征在于：所述计算垂直分段在下平面的投影点集，包括：计算垂直分段Qi在其下平面的投影点集Qi′；将垂直分段Qi中的树干点云投影至垂直分段Qi的下平面πi得到一个平面点集Qi′；计算平面点集Qi′的凸包多边形的质心点作为平面点集Qi′的几何中心点ci。5.根据权利要求4所述的一种基于激光点云数据的树干材积获取方法，其特征在于：所述计算垂直分段的断面轮廓点，包括：设点pm为垂直分段Qi中的一个点，pm′表示点pm在垂直分段Qi下平面的投影，则几何中心点ci与投影点pm′构成的向量与X轴的夹角可表示为：根据角度分区的角度值θ，将平面点集Qi′中的点分为个角度分区，每个角度分区Qij中的点可表示为：计算每一个角度分区Qij中点的重心点，所述重心点为所述角度分区中所有点的几何平均值；并将所述重心点作为所述角度分区的轮廓点，由此得到垂直分段Qi的个轮廓点。6.根据权利要求5所述的一种基于激光点云数据的树干材积获取方法，其特征在于：构建断面轮廓曲线并计算树干断面的断面积，包括：设垂直分段Qi的断面轮廓点集合用Rprofile表示，首先采用非...

【专利技术属性】
技术研发人员：尤磊，宋新宇，冯岩，郭颂，李小可，
申请(专利权)人：信阳师范学院，
类型：发明
国别省市：河南,41