一种无损估算立木材积的低成本树木结构自动重建方法技术

本发明专利技术公开了一种无损估算立木材积的低成本树木结构自动重建方法，具体为基于地面摄影测量点云重建树木三维模型的方法。针对现有立木材积估算的弊端和不足，为了降低树干体积估测成本，在不损伤树木的情况下，本发明专利技术提供了一个基于地面摄影测量点云的立木材积估算的方法，使用该方法重建树木三维几何结构并提取树干体积，为立木材积无损估测提供了一个更加经济的方法。本发明专利技术的研究目标包括：(1)开放AdQSM，开发AdQSM的自动建模方法，有利于更多人使用；(2)开发一种新的低成本树结构建模技术。术。术。

【技术实现步骤摘要】
一种无损估算立木材积的低成本树木结构自动重建方法
一、

[0001]本专利技术涉及一种基于地面摄影测量点云的用于估算立木材积的低成本树木结构自动重建方法。
二、技术背景
[0002]以较高的精度和较低的成本获得立木材积是林业研究与应用的重要问题。从三维点云中重建树木模型对于估算树木体积、森林蓄积量、生物量和碳储量具有非常重要的意义。当前，森林激光雷达点云的获取成本仍然高于森林摄影测量点云。砍伐树木是获取体积最直接和有效的办法，但是具有破坏性并且成本高。通过破坏性采样建立的材积方程是一种估算体积的间接测量方法，但是有时候测量结果的不确定性难以量化甚至未知。本专利技术通过改进的基于激光雷达点云的QSM(Quantitative Structure Model,QSM)重建算法用于从地面摄影测量点云提取树干体积。对已有的QSM方法(AdQSM)进行改进，并提供了一种低成本的树木结构自动重建方法和技术路线。将基于TLS点云的QSM算法用于建模树木的地面近景摄影测量的点云将是一个新的尝试。将来自地面近景摄影测量的点云作为树木定量结构模型的输入数据会降低树木3D重建的成本，因为摄影测量点云更容易获取。基于QSM方法建模来自树木图像的点云的相关报道很少，尤其将地面近景摄影测量点云作为输入数据。激光雷达点云和摄影测量点云都是由一系列的坐标点组成，点云中都包含了X轴,Y轴和Z轴三个方向的数据。本专利技术将QSM算法用于建模树木的地面近景摄影测量的点云，除了可以提高QSM算法的通用性，还为QSM的增加了新的输入数据。本文与已有的方法相比优势如下：
[0003](1)开发了一个基于地面近景摄影测量点云的树木结构自动重建方法，首次基于TLS点云的树木三维重建算法去处理地面摄影测点云。
[0004](2)证明了QSM的数据源可以是很容易获取的地面摄影测量点云，为树干体积无损估测提供了一个更加经济的方法。
[0005](3)在树木结构重建方面，模型更加接近树木真实生长状态，可以以低成本的方式获取立木材积。
三、
技术实现思路

[0006]针对现有立木材积估算的弊端和不足，为了降低树干体积估测成本，在不损伤树木的情况下，本专利技术提供了一个基于地面摄影测量点云的立木材积估算的方法，使用该方法重建树木三维几何结构并提取树干体积，为立木材积无损估测提供了一个更加经济的方法。
[0007]主要
技术实现思路
：
[0008]一种基于地面摄影测量点云重建树木三维模型的方法。本专利技术的技术方案主要包含以下步骤：
[0009]1.采集立木地面摄影测量点云。使用普通的数码相机采集立木，获取所有树木的坐标位置；然后生成基于图像的三维点云，采用运动重构法构造密集点云(SfM)，使用软件
自动实现SfM过程并进行人工编辑和修正；最后进行单木分割，从样地点云中提取单棵树木，进行地面点滤波、重采样、生成DEM、归一化等操作步骤实现从样地级别的点云中提取单木。分割后可对单棵点云树木进行逐一检查，对分类不正确的树木进行人工重新分割，手工去除噪点。
[0010]2.基于AdQSM从单棵树木的地面摄影测量点云中重建树木结构并提取体积。树木定量结构模型是基于地基激光雷达点云的单棵树木枝干几何结构的3D模型。利用AdQSM重建的地上单木三维几何结构去测量树干或者树枝的体积，可以量化材积方程带来的不确定性。
[0011](1)S树木主干自动识别
[0012]本专利技术采用改进的AdQSM自下而上重建树枝几何体的圆柱体，需要精确计算树干底部的初始圆柱体半径。通过自动选择树干的点云，使用不同颜色标记点云，在树干底部安装第一个圆柱体。为了减少用户的干预，本专利技术实现了自动重建方法。本专利技术可以准确识别和分割树干点云。基于树木骨架节点的分布，结合枝下高的分布，识别出树干点云变化相对稳定的位置。
[0013]算法的设计和实现过程如下：
[0014]①
从主干骨架的每个分支节点向上取20厘米，将整个主干的点云水平切成5厘米厚的切片。每个分支节点将有四个切片。骨架分支节点的搜索范围为树高的10
‑
55％。得到该范围内所有分支顶点的切片，并计算所有切片的边界框面积。
[0015]②
由于主干骨架的第一个分叉节点更有可能是主干的顶端，所以首先判断第一个分叉节点。以第一个分叉节点处的第一个切片的边界框区域作为参考点。如果其他三个切片的边界框面积变化显著(1.3倍以上)，则使用第一个分叉节点作为主干的顶部断点P。
[0016]③
如果不符合
②
的情况，遵循“先比较自己，后比较他人”的原则。以第i个节点(i＝2,3,4，
…
)的第一个切片的边界框面积为参考，如果其他三个切片的边界框面积变化显著(1.3倍以上)，则以第i个节点为主干的顶端端点P。
[0017]④
如果第i个节点不满足
③
条件，将第i个节点的1片边界框面积与第i
‑
1个节点的第1片边界框面积进行比较。如果有较大的变化(1.3倍或更多)，则认为第i个节点是主干的顶部端点P。
[0018]⑤
如果所有分叉节点都不满足
④
条件，则取第一个节点为分叉点，即主干分叉的顶端端点P。
[0019]⑥
为了保证参与拟合初始圆柱体的树木主干点云更可靠，将P点位置向下移动底部由分开。
[0020](2)初始圆柱半径聚类
[0021]为了拟合位于主干底部的第一个圆柱体，(1)设计并实现了一种算法，自动选择主干点云的相对稳定段。识别并标记点云变化相对稳定的树干，找到树干顶端端点P。对初始圆柱体进行拟合，以精确计算其半径。在本专利技术中，识别出的树干点云被不同长度(10cm，20cm，30cm，40cm，
…
，100cm)等分。根据标记点云每段拟合的圆柱体半径分布规律，对初始圆柱体半径进行聚类。由计算得出，聚类中参与初始聚类圆柱体半径的圆柱体数目
为N。L为树干点云的总长度。
[0022]树枝结构拟合和优化过程下：。
[0023]①
采用(1)的方法，从下到上搜索树干的顶部端点P，将树干的点云用固定长度均匀分割，得到无间隔的N段点云。
[0024]②
在分段点云(用于拟合初始圆柱体)中，将z值最大的点设置为顶面中心，将z值最小的点设置为底面中心。以被分割的云的最小包络盒的长度或宽度的一半(两者的最大值)为半径，确定初始圆柱体，包括圆柱体顶面中心、圆柱体底面中心和半径。
[0025]③
确定圆柱体段的中心轴线向下延伸后与地面的交点，取交点的平均值作为树的底部位置。在得到圆柱的粗半径后，利用非线性最小二乘法计算圆柱各分段的半径。
[0026]④
将N段等长点云拟合成N个圆柱体。采用K
‑
means方法对所有圆柱体半径进行聚类，尽可能保证拟合树干初始圆柱体半径的准确性。我们将圆柱体半径的大小分为大、中、小三个层次。聚本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于地面摄影测量点云重建树木三维模型的方法，其特征在于：降低了树干体积估测成本，在不损伤树木的情况下本发明提供了一个基于地面摄影测量点云的立木材积估算的方法，使用该方法重建树木三维几何结构并提取树干体积，为立木材积无损估测提供了一个更加经济的方法。2.根据权利要求1所述的一种基于地面摄影测量点云重建树木三维模型的方法，其特征在于，所述该种基于地面摄影测量点云重建树木三维模型的方法的技术方案：(1)采集立木地面摄影测量点云，使用普通的数码相机采集立木，获取所有树木的坐标位置；然后生成基于图像的三维点云，采用运动重构法构造密集点云(SfM)，使用软件自动实现SfM过程并进行人工编辑和修正；最后进行单木分割，从样地点云中提取单棵树木，进行地面点滤波、重采样、生成DEM、归一化等操作步骤实现从样地级别的点云中提取单木，分割后可对单棵点云树木进行逐一检查，对分类不正确的树木进行人...

【专利技术属性】
技术研发人员：范光鹏，苏晓慧，许诺，
申请(专利权)人：北京林业大学，
类型：发明
国别省市：