CN116432873B - 树木修剪废弃物碳排放预测方法、系统、电子设备和介质 - Google Patents
树木修剪废弃物碳排放预测方法、系统、电子设备和介质 Download PDFInfo
- Publication number
- CN116432873B CN116432873B CN202310700457.2A CN202310700457A CN116432873B CN 116432873 B CN116432873 B CN 116432873B CN 202310700457 A CN202310700457 A CN 202310700457A CN 116432873 B CN116432873 B CN 116432873B
- Authority
- CN
- China
- Prior art keywords
- tree
- branches
- trimming
- volume
- branch
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000013138 pruning Methods 0.000 title claims abstract description 119
- 238000000034 method Methods 0.000 title claims abstract description 65
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 61
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 61
- 239000002699 waste material Substances 0.000 title claims abstract description 44
- 238000009966 trimming Methods 0.000 claims abstract description 192
- 241000282806 Rhinoceros Species 0.000 claims abstract description 69
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 230000011218 segmentation Effects 0.000 claims abstract description 16
- 238000004891 communication Methods 0.000 claims description 16
- 238000007781 pre-processing Methods 0.000 claims description 7
- 238000012952 Resampling Methods 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 239000002028 Biomass Substances 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000000750 progressive effect Effects 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009264 composting Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/005—Tree description, e.g. octree, quadtree
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/188—Vegetation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/84—Greenhouse gas [GHG] management systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Tourism & Hospitality (AREA)
- Software Systems (AREA)
- General Business, Economics & Management (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Water Supply & Treatment (AREA)
- Entrepreneurship & Innovation (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请属于绿地生物量计算领域，本申请实施例提供了一种树木修剪废弃物碳排放预测方法、系统、电子设备和介质，树木修剪废弃物碳排放预测方法包括获取由地面激光雷达扫描仪采集到的目标区域的树木群落点云数据；对点云数据进行点云分割得到目标树木的点云文本文件，并利用定量结构模型将其转换成3D模型文件格式的目标树木的三维结构模型；根据在Rhinoceros 3D建模软件中预先设置的修剪条件以及修剪参数模拟出目标树木的修剪后的形态并计算出目标树木的修剪体积；以单棵树木为单位遍历目标区域，得到树木群落的预测修剪体积并计算修剪废弃物产生的碳排放，本申请的方法能较为准确的预测出树木修剪废弃物的碳排放，且整个预测过程简单，容易实现。
Description
技术领域
本申请属于绿地生物量计算领域，尤其涉及一种树木修剪废弃物碳排放预测方法、系统、电子设备和介质。
背景技术
城市树木是城市绿地最重要的组成部分，然而城市树木的隐含碳排放有可能抵消城市树木的碳汇效益，其中树木修剪下来的废弃枝条是城市树木重要的碳排放来源。然而目前研究树木修剪产生的碳排放时，通常使用固定的修剪比例预测树木修剪体积从而计算树木修剪废弃物产生的碳排放，这往往忽视了树木年龄、形态和修剪方式的差异，造成树木修剪体积的预测不准确进而导致树木修剪废弃物碳排放预测不准确。此外，目前也有研究直接采用对待修剪的枝条实际测量的方式来预测修剪体积进而预测出树木修剪废弃物碳排放，但该方法成本较高，效率较低，且具有实施难度。
因此，如何准确高效的预测出树木修剪废弃物碳排放，这是急需解决的问题。
发明内容
有鉴于此，本申请实施例所解决的技术问题之一在于提供一种树木修剪废弃物碳排放预测方法、系统、电子设备和介质，用以尽量克服或避免上述问题。
基于上述目的，本申请提供一种树木修剪废弃物碳排放预测方法，包括：获取由地面激光雷达扫描仪采集到的目标区域的树木群落点云数据；对所述点云数据进行点云分割，得到目标树木的点云文本文件，并利用定量结构模型将所述点云文本文件转换成3D模型文件格式的所述目标树木的三维结构模型；利用Rhinoceros 3D建模软件读取所述目标树木的三维结构模型；根据在所述Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系模拟出所述目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积；将所述主干、二级枝条、三级枝条的体积之和作为所述目标树木的整树体积，将所述主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积之和作为所述目标树木的修剪后的体积；将所述目标树木的整树体积与修剪后的体积作差，得到所述目标树木的修剪体积；以单棵树木为单位，遍历所述目标区域，对所述目标区域内的每棵树木的修剪体积进行计算并求和，得到所述目标区域内树木群落的预测修剪体积；根据所述目标区域内树木群落的预测修剪体积计算所述目标区域内树木群落的修剪废弃物产生的碳排放。
可选地，在所述对所述点云数据进行点云分割之前，还包括：对所述点云数据进行预处理，所述预处理至少包括点云重采样、点云去噪、地面点分类、点云归一化。
可选地，在所述利用Rhinoceros 3D建模软件读取所述目标树木的三维结构模型之后，还包括：在所述Rhinoceros 3D建模软件中根据所述目标树木的三维结构模型构建所述目标树木的Mesh三角网格面模型；根据所述目标树木的树木结构将所述Mesh三角网格面模型划分为若干个封闭的三角网格体，所述三角网格体用于区分出所述目标树木的一种树木结构的整体，所述树木结构至少包括主干、二级枝条、三级枝条；所述根据在所述Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系模拟出所述目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积，包括：根据所述三角网格体区分出来的所述树木结构的整体以及在所述Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系，模拟出所述目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积。
可选地，所述根据在所述Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系模拟出所述目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积，包括：将体积小于预设的有效枝条体积阈值的所述三角网格体进行删除，得到所述目标树木的有效枝条；在所述有效枝条中划分出所述目标树木的主干并计算出所述主干的体积；对于划分出主干后的所述有效枝条，判断每个所述有效枝条是否与所述主干相交；根据枝条间的相交关系，将判断结果为“是”的所述有效枝条判定为二级枝条并计算出所述二级枝条的体积；根据在所述Rhinoceros 3D建模软件中预先设置的二级枝条对应的修剪条件以及修剪参数对所述二级枝条进行修剪，得到二级枝条修剪后剩余部分的三角网格体并计算所述二级枝条修剪后剩余部分的体积；根据在所述Rhinoceros 3D建模软件中预先设置的主干对应的修剪条件以及修剪参数对所述主干进行修剪，得到主干修剪后剩余部分的三角网格体并计算所述主干修剪后剩余部分的体积；对于划分出主干以及二级枝条后的所述有效枝条，判断每个所述有效枝条是否与所述二级枝条相交；根据枝条间的相交关系，将判断结果为“是”的所述有效枝条判定为三级枝条并计算出所述三级枝条的体积；根据在所述Rhinoceros 3D建模软件中预先设置的三级枝条对应的修剪条件以及修剪参数对所述三级枝条进行修剪，得到三级枝条修剪后剩余部分的三角网格体并计算所述三级枝条修剪后剩余部分的体积；根据所述主干修剪后剩余部分、所述二级枝条修剪后剩余部分、所述三级枝条修剪后剩余部分的三角网格体模拟出所述目标树木的修剪后的形态。
可选地，在所述根据在所述Rhinoceros 3D建模软件中预先设置的二级枝条对应的修剪条件以及修剪参数对所述二级枝条进行修剪之前，还包括：计算有效二级枝条体积阈值，所述有效二级枝条体积阈值为预设的有效二级枝条系数与体积最大的所述二级枝条的体积的乘积；判断每个所述二级枝条的体积是否小于所述有效二级枝条体积阈值；将判断结果为“是”的所述二级枝条进行删除。
可选地，所述根据在所述Rhinoceros 3D建模软件中预先设置的三级枝条对应的修剪条件以及修剪参数对所述三级枝条进行修剪之前，还包括：计算每一个三级枝条对应的有效三级枝条体积阈值，所述三级枝条对应的有效三级枝条体积阈值为预设的有效三级枝条系数与和所述三级枝条相连接的二级枝条的体积的乘积；判断每个所述三级枝条的体积是否小于所述三级枝条对应的有效三级枝条体积阈值；将判断结果为“是”的所述三级枝条进行删除。
可选地，所述碳排放的计算公式如下：
C=V×ρ×i×EF
其中，C表示所述目标区域内树木群落的修剪废弃物产生的碳排放量，V表示所述目标区域内树木群落的预测修剪体积，ρ表示被修剪的枝条的平均密度，i表示被修剪的枝条的含碳量系数，EF表示碳排放因子。
本申请实施例还提供了一种树木修剪废弃物碳排放预测系统，包括：数据获取模块，用于获取由地面激光雷达扫描仪采集到的目标区域的树木群落点云数据；三维结构模型生成模块，用于对所述点云数据进行点云分割，得到目标树木的点云文本文件，并利用定量结构模型将所述点云文本文件转换成3D模型文件格式的所述目标树木的三维结构模型；单棵树木修剪体积计算模块，利用Rhinoceros 3D建模软件读取所述目标树木的三维结构模型；根据在所述Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系模拟出所述目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积；将所述主干、二级枝条、三级枝条的体积之和作为所述目标树木的整树体积，将所述主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积之和作为所述目标树木的修剪后的体积；将所述目标树木的整树体积与修剪后的体积作差，得到所述目标树木的修剪体积；目标区域树木修剪体积之和计算模块，用于以单棵树木为单位，遍历所述目标区域，对所述目标区域内的每棵树木的修剪体积进行计算并求和，得到所述目标区域内树木群落的预测修剪体积；树木修剪废弃物碳排放预测模块，用于根据所述目标区域内树木群落的预测修剪体积计算所述目标区域内树木群落的修剪废弃物产生的碳排放。
本申请实施例还提供了一种电子设备，包括：处理器、存储器、通信接口和通信总线，所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信；所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行如上所述方法对应的操作。
本申请实施例还提供了一种计算机存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述方法。
由以上技术方案可见，首先获取由地面激光雷达扫描仪采集到的目标区域的树木群落点云数据；对点云数据进行点云分割，得到目标树木的点云文本文件，并利用定量结构模型将点云文本文件转换成3D模型文件格式的目标树木的三维结构模型；再将三维结构模型导入Rhinoceros 3D建模软件中；根据在Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系模拟出目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积；将主干、二级枝条、三级枝条的体积之和作为目标树木的整树体积，将主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积之和作为目标树木的修剪后的体积；将目标树木的整树体积与修剪后的体积作差，得到目标树木的修剪体积；以单棵树木为单位，遍历目标区域，对目标区域内的每棵树木的修剪体积进行计算并求和，得到目标区域内树木群落的预测修剪体积，最后根据目标区域内树木群落的预测修剪体积计算目标区域内树木群落的修剪废弃物产生的碳排放。本申请通过地面激光雷达扫描仪采集到的点云数据计算修剪体积，大大提高了计算效率，采用Rhinoceros 3D建模软件对树木模型进行处理，能够直观的展示给用户树木的修剪后的形态，用户可以根据该树木形态对树木进行修剪，同时Rhinoceros 3D建模软件能够根据修剪条件以及修剪参数进行遍历直接得出每个树木的各个枝条的体积，用户可以根据这些数据预测出修剪废弃物产生的碳排放，整个预测过程简单高效，容易实现，且预测结果较为准确。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请实施例中记载的一些实施例，对于本领域普通技术人员来讲，还可以根据这些附图获得其他的附图。
图1为本申请实施例一的一种树木修剪废弃物碳排放预测方法的流程图；
图2为本申请实施例二的一种树木修剪方法的流程图；
图3为本发明实施例三的一种树木修剪废弃物碳排放预测系统的示意图；
图4为本申请实施例四的一种电子设备的结构示意图。
具体实施方式
为了使本领域的人员更好地理解本申请实施例中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本申请实施例一部分实施例，而不是全部的实施例。基于本申请实施例中的实施例，本领域普通技术人员所获得的所有其他实施例，都应当属于本申请实施例保护的范围。
需要说明的是，实施本申请实施例的任一技术方案必不一定需要同时达到以上的所有优点。下面结合本申请实施例附图进一步说明本申请实施例具体实现。
实施例一
图1为本申请实施例的一种树木修剪废弃物碳排放预测方法的流程图，如图1所示，该树木修剪废弃物碳排放预测方法包括如下步骤：
步骤S101、获取由地面激光雷达扫描仪采集到的目标区域的树木群落点云数据。
具体地，目标区域为用户需要对树木进行修剪体积预测的区域，树木群落点云数据是指在三维坐标系中树木的空间特征数据，树木群落点云数据具有几何位置信息以及颜色信息等，为了使得采集到的树木的较小枝条的较为精确的点云数据，可以选择相对精度为1CM以内的地面激光雷达扫描仪采集目标区域的树木群落点云数据。
需要说明的是，在进行树木群落点云数据采集时，为了避免树木叶片影响后续的树木体积计算以及分析等，可以选择目标区域的树木在无叶期时进行点云数据的采集。
步骤S102、对点云数据进行点云分割，得到目标树木的点云文本文件，并利用定量结构模型将点云文本文件转换成3D模型文件格式的目标树木的三维结构模型。
具体地，这里点云分割是将树木群落的点云数据根据树木的空间、几何等特征点进行划分，使得同一划分内的点云拥有相似的特征，将目标树木从树木群落中划分出来，从而便于对其单独进行处理。树木定量结构模型(Quantitative Structure Model,QSM)是基于地基激光雷达点云数据对目标树木的结构精细建模方法，是以层次结构顺序描述完整的地上树木木质部位的单棵树木的几何模型。点云文本文件为XYZ格式，以目标树木生长的地面水平面为XY面，垂直于地面的树干方向为法线方向（即地面向上方向为Z轴正方向）。该定量结构模型能够根据目标树木的点云文本文件实现对目标树木的快速建模，将点云文本文件转换成3D模型文件格式的目标树木的三维结构模型，这里的3D模型文件格式可以为OBJ文件格式，即标准3D模型文件格式，定量结构模型还可以从目标树木的三维结构中准确的计算出目标树木建模体积，包括主干体、枝条等的体积和长度等。
可选地，在对点云数据进行点云分割之前，还可以对点云数据进行预处理，预处理包括点云重采样、点云去噪、地面点分类、点云归一化等。
具体地，点云重采样可以采用5%的采样率将点云进行重采样处理，随机保留指定数量的点，这里保留的点数为点云总采样数与采样率之积，点云重采样能够减少点云数量，精简数据，有助于节省储存空间和计算时间。点云数据常见噪声包括高位粗差和低位粗差，高位粗差通常是因为机载雷达系统在采集数据的过程中受到低飞的飞行物，比如鸟类或者飞机的影响，误将这些物体反射回来的信号当作被测树木的反射信号记录下来，低位粗差则是由于测量过程中的多路径误差或者激光测距仪的误差导致产生的极低点，地面激光雷达一般多产生低位粗差。这些噪点会对树木体积的计算产生影响，因此需要移除噪点来提高数据质量。可以对每一个点搜索指定邻域点个数的相邻点，计算点到相邻点的距离平均值以及这些平均值距离的中值和标准差，如果这个点的平均值距离大于最大距离(其中，最大距离=中值+标准差倍数×标准差)，则认为是噪点，应将该点去除。地面点分类是对点云数据进行过滤的过程，通常采用改进的渐进加密三角网滤波算法分类地面点，该算法的具体步骤可以参考现有的渐进加密三角网滤波算法实施步骤，本申请对此不作限制。点云归一化是指将点云数据进行标准化处理，以消除地形起伏对点云的影响，生成归一化的点云数据。
本实施方式中，在对点云数据进行点云分割之前，对点云数据进行预处理，预处理包括点云重采样、点云去噪、地面点分类、点云归一化等，能够节省储存空间和计算时间，提高整个预测过程的效率，同时提高数据质量，消除地形起伏对点云的影响，使得预测结果更为准确。
步骤S103、利用Rhinoceros 3D建模软件读取目标树木的三维结构模型。
具体地，Rhinoceros 3D建模软件是一款功能强大的高级建模软件，可以创建、编辑、分析、提供、渲染、动画与转换 NURBS 线条、曲面、实体与多边形网格，且不受精度、复杂阶数或是尺寸的限制，这里将目标树木的三维模型导入Rhinoceros 3D建模软件中，Rhinoceros 3D建模软件读取目标树木的三维结构模型，并将该三维模型结构数据拾取进入Grasshopper中，Grasshopper是一款基于Rhinoceros 3D建模软件运行的可视化编程语言，并利用Grasshopper平台中GhPython Script来进行后续的树木修剪体积的核算、修剪后树木的三维可视化模拟以及修剪废弃物碳排放的预测等工作。
步骤S104、根据在Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系模拟出目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积。
具体地，修剪条件用于确定目标树木进行修剪的枝条、枝条上需要进行修剪的枝条位置以按照何种方式对枝条进行修剪等，修剪参数用于规定修剪枝条的修剪量。例如，对二级枝条进行修剪时，可以先确定该所在的坐标系中所述二级枝条Z方向的最高点与最低点的XYZ坐标以及该二级枝条的体心XYZ坐标，计算最高点与最低点的距离并将计算结果作为所述二级枝条修剪标准距离，假设该修剪参数为1/3，则可以将该体心进行移动，直至该体心与最低点的距离为所述二级枝条修剪标准距离的1/3，创建原点为上述移动后的体心、法线方向为最低点至最高点的矢量方向的平面，以该平面作为该二级枝条的修剪面并对该二级枝条进行修剪。枝条间的相交关系可以判断出枝干间的分级拓扑关系，例如，与主干相交的枝条可以判断其为二级枝条。
步骤S105、将主干、二级枝条、三级枝条的体积之和作为目标树木的整树体积，将主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积之和作为目标树木的修剪后的体积。
具体地，除主干、二级枝条、三级枝条以外的枝条体积过小，对目标树木的整树体积的计算结果影响不大，故可以将这些过小枝条进行删除，不计算在整树枝条的范围内，可以减少计算过程，提高树木修剪体积的预测效率，进而提高修剪废弃物碳排放的预测效率。
步骤S106、将目标树木的整树体积与修剪后的体积作差，得到目标树木的修剪体积。
步骤S107、以单棵树木为单位，遍历目标区域，对目标区域内的每棵树木的修剪体积进行计算并求和，得到目标区域内树木群落的预测修剪体积。
具体地，以单棵树木为单位，以此对目标区域的每棵树木进行点云分割，得到单棵树木的点云文本文件，并利用定量结构模型将点云文本文件转换成3D模型文件格式的单棵树木的三维结构模型，再利用Rhinoceros 3D建模软件对单棵树木进行修剪体积进行计算，遍历目标区域的所有树木，即可得到目标区域内树木群落的预测修剪体积。
步骤S108、根据目标区域内树木群落的预测修剪体积计算目标区域内树木群落的修剪废弃物产生的碳排放。
可选地，碳排放的计算公式可以是：
C=V×ρ×i×EF
其中，C表示目标区域内树木群落的预测修剪体积修剪废弃物产生的碳排放量，单位为Kg；V表示目标区域内树木群落的预测修剪体积，单位为;ρ表示被修剪的枝条的平均密度，单位kg//>;i表示被修剪的枝条的含碳量系数，EF表示被修剪枝条的处理方式的碳排放因子,单位为kg-C。
具体地，树木修剪部分的生物量转化为碳排放量是通过生物量的含碳量系数和修剪后废弃物处理的碳排放因子确定的，这里的被修剪的枝条的含碳量系数通常可以设置为0.5。由于树木修剪部分有多种处理方式，不同处理方式的碳排放因子通常是不同的，主要的树木被修剪枝条的处理方式是焚烧、填埋、生物质发电和堆肥，可以根据不同的处理方式在碳排放因子相应的表格中进行查询，从而较为准确的计算出修剪废弃物所产生的碳排放量。
由以上技术方案可见，首先获取由地面激光雷达扫描仪采集到的目标区域的树木群落点云数据；对点云数据进行点云分割，得到目标树木的点云文本文件，并利用定量结构模型将点云文本文件转换成3D模型文件格式的目标树木的三维结构模型；再将三维结构模型导入Rhinoceros 3D建模软件中；根据在Rhinoceros 3D建模软件中预先设置的修剪条件以及修剪参数模拟出目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积；将主干、二级枝条、三级枝条的体积之和作为目标树木的整树体积，将主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积之和作为目标树木的修剪后的体积；将目标树木的整树体积与修剪后的体积作差，得到目标树木的修剪体积；以单棵树木为单位，遍历目标区域，对目标区域内的每棵树木的修剪体积进行计算并求和，得到目标区域内树木群落的预测修剪体积，最后根据目标区域内树木群落的预测修剪体积和计算目标区域内树木群落的修剪废弃物产生的碳排放。本申请通过地面激光雷达扫描仪采集到的点云数据计算修剪体积，大大提高了计算效率，采用Rhinoceros 3D建模软件对树木模型进行处理，能够直观的展示给用户树木的修剪后的形态，用户可以根据该树木形态对树木进行修剪，同时Rhinoceros 3D建模软件能够根据修剪条件以及修剪参数进行遍历直接得出每个树木的各个枝条的体积，用户可以根据这些数据预测出修剪废弃物碳排放，整个预测过程简单高效，容易实现，且预测结果较为准确。
实施例二
图2为本申请实施例的一种树木修剪方法的流程图，如图2所示，该树木修剪方法包括如下步骤：
步骤S201、将体积小于预设的有效枝条体积阈值的三角网格体进行删除，得到目标树木的有效枝条。
可选地，在利用Rhinoceros 3D建模软件读取目标树木的三维结构模型之后，还可以在Rhinoceros 3D建模软件中根据目标树木的三维结构模型构建目标树木的Mesh三角网格面模型；根据目标树木的树木结构将Mesh三角网格面模型划分为若干个封闭的三角网格体，三角网格体用于区分出目标树木的一种树木结构的整体，树木结构至少包括主干、二级枝条、三级枝条，根据在Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系模拟出目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积，包括：根据三角网格体区分出来的树木结构的整体以及在Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系，模拟出目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积。
具体地，Mesh三角网格面模型是计算机图形学中用于为各种不规则物体建立模型的一种数据结构，现实世界中的物体表面直观上看都是由曲面构成的，而在计算机世界中，由于只能用离散的结构去模拟现实中连续的事物，所以现实世界中的曲面实际上在计算机里可以是由无数个小的多边形面片组成，Mesh三角网格面模型由若干个网格面封闭而成，Rhinoceros 3D建模软件根据目标树木的树木结构自动将Mesh三角网格面模型划分为若干个封闭的三角网格体，三角网格体用于区分出目标树木的一种树木结构的整体，例如代表主干、二级枝条、三级枝条或其他较小枝条等。
本实施方式中，根据目标树木的三维结构模型构建目标树木的Mesh三角网格面模型，再根据目标树木的树木结构自动将Mesh三角网格面模型划分为若干个封闭的三角网格体，三角网格体用于区分出目标树木的一种树木结构的整体，树木结构至少包括主干、二级枝条、三级枝条，构建目标树木的Mesh三角网格面模型能够利用较少的数据规模表示出精细的三维模型细节，且不同的树木结构利用不同的三角网格体进行表示，方便用户对树木结构直接进行删除。
具体地，目标树木上存在很多体积过小的无效枝条，且无效枝条数量较多计算量较大，但这些无效枝条对树木体积的影响又较小，因此可以直接对其进行删除，因此可以设置有效枝条体积阈值，将无效枝条进行删除，例如，有效枝条体积阈值可以设置为，判断目标树木中的三角网格体的体积是否小于有效枝条体积阈值，对体积小于/>的三角网格体进行删除，以得到目标树木的有效枝条。
步骤S202、在有效枝条中划分出目标树木的主干并计算出主干的体积。
具体地，上述有效枝条中对目标树木的主干进行筛选，得到目标树木的主干，并利用Rhinoceros 3D建模软件自动计算出主干的体积。
步骤S203、对于划分出主干后的有效枝条，判断每个有效枝条是否与主干相交。
具体地，对于划分出主干后的有效枝条，可以判断枝干间的分级拓扑关系，判断除主干外的每个有效枝条是否与主干相交。
步骤S204、根据枝条间的相交关系，将判断结果为“是”的有效枝条判定为二级枝条并计算出二级枝条的体积。
具体地，将与主干相交的有效枝条判定为二级枝条，判断结果为“不相交”的有效枝条判定为其他有效枝条，并利用Rhinoceros 3D建模软件自动计算出二级枝条的体积。
步骤S205、根据在Rhinoceros 3D建模软件中预先设置的二级枝条对应的修剪条件以及修剪参数对二级枝条进行修剪，得到二级枝条修剪后剩余部分的三角网格体并计算二级枝条修剪后剩余部分的体积。
可选地，在根据在Rhinoceros 3D建模软件中预先设置的二级枝条对应的修剪条件以及修剪参数对二级枝条进行修剪之前，还可以计算有效二级枝条体积阈值，有效二级枝条体积阈值为预设的有效二级枝条系数与体积最大的二级枝条的体积的乘积，并判断每个二级枝条的体积是否小于有效二级枝条体积阈值，将判断结果为“是”的二级枝条进行删除。
需要说明的是，在进行枝条修剪时，存在一些体积较小的二级枝条，因枝条体积过小，无需对其进行修剪，可以直接对这些枝条整体进行删除，也方便了对枝条进行修剪以及同时比较少了二级枝条修剪的计算数量，保证修剪结果的准确性，
具体地，在计算出二级枝条的体积之后，基于二级枝条集合，可以对其进行排序并从中筛选出体积最大的二级枝条，然后计算有效二级枝条体积阈值，有效二级枝条系数可以根据实际树木情况进行设置，例如，这里的有效二级枝条系数可以设置为0.1，此处对具体设置数值不做限制，再根据该有效二级枝条系数与体积最大的二级枝条的体积的乘积计算有效二级枝条体积阈值，对于体积小于有效二级枝条体积阈值的二级枝条进行删除，再对删除掉无效二级枝条后的二级枝条进行枝条修剪。
本实施方式中，根据最大体积的二级枝条以及有效二级枝条系数对二级枝条进行筛选，删除体积较小的无效二级枝条，减少了二级枝条的修剪数量并提高了修剪结果的准确性。
具体地，对二级枝条进行修剪时，应确定需要修剪的二级枝条在该目标树木所在的坐标系中位置，确定该二级枝条在Z方向的最高点与最低点的XYZ坐标以及该二级枝条的体心XYZ坐标，该体心即为该二级枝条的重心，计算最高点与最低点的坐标距离并将计算结果作为所述二级枝条修剪标准距离，并以该二级枝条最低点至最高点的矢量方向的反方向作为体心移动的移动方向，根据预设的修剪参数确定该二级枝条的修剪面，例如，设置修剪参数为1/3，将该二级枝条的体心沿着移动方向进行移动，并实时判断该体心与二级枝条最低点坐标距离是否大于标准距离的1/3，若大于，则继续移动体心，直至该体心与二级枝条最低点坐标距离是等于标准距离的1/3，停止移动体心，并创建原点为上述移动后的体心，法线方向为最低点至最高点的矢量方向的平面，以该平面作为该二级枝条的修剪面并对该二级枝条进行修剪，并将修剪后的三角网格体进行填洞转化，转化为可计算体积的封闭的三角网格体，然后根据Rhinoceros 3D建模软件计算出二级枝条修剪后剩余部分的体积。
步骤S206、根据在Rhinoceros 3D建模软件中预先设置的主干对应的修剪条件以及修剪参数对主干进行修剪，得到主干修剪后剩余部分的三角网格体并计算主干修剪后剩余部分的体积。
具体的，对目标树木的主干进行修剪时，可以将二级枝条的最低点中的最高点作为主干修剪面的坐标原点，以目标树木的Z方向作为法线方向来确定主干修剪面，并对其进行修剪，并将修剪后的三角网格体进行填洞转化，转化为可计算体积的封闭的三角网格体，然后根据Rhinoceros 3D建模软件计算出主干修剪后剩余部分的体积。
步骤S207、对于划分出主干以及二级枝条后的有效枝条，判断每个有效枝条是否与二级枝条相交。
步骤S208、根据枝条间的相交关系，将判断结果为“是”的有效枝条判定为三级枝条并计算出三级枝条的体积。
具体地，将判断结果为“是”的有效枝条判定为三级枝条，将判断结果为“否”的枝条进行删除，并根据Rhinoceros 3D建模软件计算出三级枝条的体积。
步骤S209、根据在Rhinoceros 3D建模软件中预先设置的三级枝条对应的修剪条件以及修剪参数对三级枝条进行修剪，得到三级枝条修剪后剩余部分的三角网格体并计算三级枝条修剪后剩余部分的体积。
可选地，根据在Rhinoceros 3D建模软件中预先设置的三级枝条对应的修剪条件以及修剪参数对三级枝条进行修剪之前，还可以计算每一个三级枝条对应的有效三级枝条体积阈值，三级枝条对应的有效三级枝条体积阈值为预设的有效三级枝条系数与和三级枝条相连接的二级枝条的体积的乘积。判断每个三级枝条的体积是否小于三级枝条对应的有效三级枝条体积阈值。将判断结果为“是”的三级枝条进行删除。
具体地，与二级枝条类似，也存在一些体积较小的三级枝条，因枝条体积过小，可以直接对这些枝条整体进行删除，在对三级枝条中的过小枝条进行删除时，可以利用与该三级枝条连接的二级枝条的体积与预设的有效三级枝条系数的乘积作为有效三级枝条体积阈值，有效三级枝条系数可以根据实际树木情况进行设置，例如可以将其设置为0.2，此处对具体设置数值不做限制，再判断每个三级枝条的体积是否小于三级枝条对应的有效三级枝条体积阈值，对于判断结果为“是”的三级枝条进行删除，再对删除掉无效三级枝条后的三级枝条进行枝条修剪。
本实施方式中，根据与和三级枝条相连接的二级枝条的体积以及有效三级枝条系数对三级枝条进行筛选，删除体积较小的无效三级枝条，减少了三级枝条的修剪数量并提高了修剪结果的准确性。
具体地，对三级枝条时，可参考上述二级枝条的修剪过程来确定三级枝条的修剪面，此处不再赘述。对三级枝条进行修剪后，还需要将修剪后的三级枝条的三角网格体进行填洞转化，转化为可计算体积的封闭的三角网格体，然后根据Rhinoceros 3D建模软件计算出三级枝条修剪后剩余部分的体积。
步骤S210、根据主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的三角网格体模拟出目标树木的修剪后的形态。
具体地，Rhinoceros 3D建模软件可以根据主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的三角网格体模拟出目标树木的修剪后的形态，同时也可以将修剪后的目标树木的形态与修剪前的目标树木的形态进行对比，直观的呈现给用户目标树木修剪后的三维可视化结果。
由以上技术方案可见，将体积小于预设的有效枝条体积阈值的三角网格体进行删除，得到目标树木的有效枝条，在有效枝条中划分出目标树木的主干并计算出主干的体积，对于划分出主干后的有效枝条，判断每个有效枝条是否与主干相交，将判断结果为“是”的有效枝条判定为二级枝条并计算出二级枝条的体积，根据在Rhinoceros 3D建模软件中预先设置的二级枝条对应的修剪条件以及修剪参数对二级枝条进行修剪，得到二级枝条修剪后剩余部分的三角网格体并计算二级枝条修剪后剩余部分的体积，根据在Rhinoceros 3D建模软件中预先设置的主干对应的修剪条件以及修剪参数对主干进行修剪，得到主干修剪后剩余部分的三角网格体并计算主干修剪后剩余部分的体积，对于划分出主干以及二级枝条后的有效枝条，判断每个有效枝条是否与二级枝条相交，将判断结果为“是”的有效枝条判定为三级枝条并计算出三级枝条的体积，根据在Rhinoceros 3D建模软件中预先设置的三级枝条对应的修剪条件以及修剪参数对三级枝条进行修剪，得到三级枝条修剪后剩余部分的三角网格体并计算三级枝条修剪后剩余部分的体积，根据主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的三角网格体模拟出目标树木的修剪后的形态，通过设定的修剪条件以及修剪参数实现了对目标树木枝条的修剪，且修剪过程简单易于实现，计算结果较为准确，直观的呈现给用户树木修剪后的形态，提高了用户体验。
实施例三
图3为本发明实施例的一种树木修剪废弃物碳排放预测系统的示意图，如图3所示，该树木修剪废弃物碳排放预测系统包括：
数据获取模块301，用于获取由地面激光雷达扫描仪采集到的目标区域的树木群落点云数据；
三维结构模型生成模块302，用于对点云数据进行点云分割，得到目标树木的点云文本文件，并利用定量结构模型将点云文本文件转换成3D模型文件格式的目标树木的三维结构模型；
单棵树木修剪体积计算模块303，利用Rhinoceros 3D建模软件读取目标树木的三维结构模型；根据在Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系模拟出目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积；将主干、二级枝条、三级枝条的体积之和作为目标树木的整树体积，将主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积之和作为目标树木的修剪后的体积；将目标树木的整树体积与修剪后的体积作差，得到目标树木的修剪体积；
目标区域树木修剪体积之和计算模块304，用于以单棵树木为单位，遍历目标区域，对目标区域内的每棵树木的修剪体积进行计算并求和，得到目标区域内树木群落的预测修剪体积；
树木修剪废弃物碳排放预测模块305，用于根据目标区域内树木群落的预测修剪体积计算目标区域内树木群落的修剪废弃物产生的碳排放。
由以上技术方案可见，首先获取由地面激光雷达扫描仪采集到的目标区域的树木群落点云数据；对点云数据进行点云分割，得到目标树木的点云文本文件，并利用定量结构模型将点云文本文件转换成3D模型文件格式的目标树木的三维结构模型；再将三维结构模型导入Rhinoceros 3D建模软件中；根据在Rhinoceros 3D建模软件中预先设置的修剪条件以及修剪参数模拟出目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积；将主干、二级枝条、三级枝条的体积之和作为目标树木的整树体积，将主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积之和作为目标树木的修剪后的体积；将目标树木的整树体积与修剪后的体积作差，得到目标树木的修剪体积；以单棵树木为单位，遍历目标区域，对目标区域内的每棵树木的修剪体积进行计算并求和，得到目标区域内树木群落的预测修剪体积，最后根据目标区域内树木群落的预测修剪体积计算目标区域内树木群落的修剪废弃物产生的碳排放。本申请通过地面激光雷达扫描仪采集到的点云数据计算修剪体积，大大提高了计算效率，采用Rhinoceros 3D建模软件对树木模型进行处理，能够直观的展示给用户树木的修剪后的形态，用户可以根据该树木形态对树木进行修剪，同时， Rhinoceros 3D建模软件能够根据修剪条件以及修剪参数进行遍历直接得出每个树木的各个枝条的体积，用户可以根据这些数据预测出修剪废弃物碳排放，整个预测过程简单高效，容易实现，且预测结果较为准确。
实施例四
图4为本申请实施例的一种电子设备的结构示意图，本申请具体实施例并不对电子设备的具体实现做限定。
该电子设备可以包括：处理器(processor)402、通信接口(CommunicationsInterface)404、存储器(memory)406、以及通信总线408。
其中：
处理器402、通信接口404、以及存储器406通过通信总线408完成相互间的通信。
通信接口404，用于与其它电子设备或服务器进行通信。
处理器402，用于执行程序410，具体可以执行上述方法实施例中的相关步骤。
具体地，程序410可以包括程序代码，该程序代码包括计算机操作指令。
处理器402可能是CPU，或者是特定集成电路ASIC（Application SpecificIntegrated Circuit），或者是被配置成实施本申请实施例的一个或多个集成电路。智能设备包括的一个或多个处理器，可以是同一类型的处理器，如一个或多个CPU；也可以是不同类型的处理器，如一个或多个CPU以及一个或多个ASIC。
存储器406，用于存放程序410。存储器406可能包含高速RAM存储器，也可能还包括非易失性存储器（non-volatile memory），例如至少一个磁盘存储器。
程序410具体可以用于使得处理器402执行前述多个方法实施例中任一实施例所描述的方法对应的操作。
程序410中各步骤的具体实现可以参见上述方法实施例中的相应步骤和单元中对应的描述，并具有相应的有益效果，在此不赘述。所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的设备和模块的具体工作过程，可以参考前述方法实施例中的对应过程描述，在此不再赘述。
本申请实施例最后还提供了一种计算机存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述多个方法实施例中的任一方法对应的操作。
需要指出，根据实施的需要，可将本申请实施例中描述的各个部件/步骤拆分为更多部件/步骤，也可将两个或多个部件/步骤或者部件/步骤的部分操作组合成新的部件/步骤，以实现本申请实施例的目的。
上述根据本申请实施例的方法可在硬件、固件中实现，或者被实现为可存储在记录介质（诸如CD、ROM、RAM、软盘、硬盘或磁光盘）中的软件或计算机代码，或者被实现通过网络下载的原始存储在远程记录介质或非暂时机器可读介质中并将被存储在本地记录介质中的计算机代码，从而在此描述的方法可被存储在使用通用计算机、专用处理器或者可编程或专用硬件（诸如ASIC或FPGA）的记录介质上的这样的软件处理。可以理解，计算机、处理器、微处理器控制器或可编程硬件包括可存储或接收软件或计算机代码的存储组件（例如，RAM、ROM、闪存等），当软件或计算机代码被计算机、处理器或硬件访问且执行时，实现在此描述的方法。此外，当通用计算机访问用于实现在此示出的方法的代码时，代码的执行将通用计算机转换为用于执行在此示出的方法的专用计算机。
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及方法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请实施例的范围。
以上实施方式仅用于说明本申请实施例，而并非对本申请实施例的限制，有关技术领域的普通技术人员，在不脱离本申请实施例的精神和范围的情况下，还可以做出各种变化和变型，因此所有等同的技术方案也属于本申请实施例的范畴，本申请实施例的专利保护范围应由权利要求限定。
Claims (10)
1.一种树木修剪废弃物碳排放预测方法，其特征在于，包括：
获取由地面激光雷达扫描仪采集到的目标区域的树木群落点云数据；
对所述点云数据进行点云分割，得到目标树木的点云文本文件，并利用定量结构模型将所述点云文本文件转换成3D模型文件格式的所述目标树木的三维结构模型；
利用Rhinoceros 3D建模软件读取所述目标树木的三维结构模型；
根据在所述Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系模拟出所述目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积；
将所述主干、二级枝条、三级枝条的体积之和作为所述目标树木的整树体积，将所述主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积之和作为所述目标树木的修剪后的体积；
将所述目标树木的整树体积与修剪后的体积作差，得到所述目标树木的修剪体积；
以单棵树木为单位，遍历所述目标区域，对所述目标区域内的每棵树木的修剪体积进行计算并求和，得到所述目标区域内树木群落的预测修剪体积；
根据所述目标区域内树木群落的预测修剪体积计算所述目标区域内树木群落的修剪废弃物产生的碳排放。
2.根据权利要求1所述的方法，其特征在于，在所述对所述点云数据进行点云分割之前，还包括：
对所述点云数据进行预处理，所述预处理至少包括点云重采样、点云去噪、地面点分类、点云归一化。
3.根据权利要求1所述的方法，其特征在于，在所述利用Rhinoceros 3D建模软件读取所述目标树木的三维结构模型之后，还包括：
在所述Rhinoceros 3D建模软件中根据所述目标树木的三维结构模型构建所述目标树木的Mesh三角网格面模型；
根据所述目标树木的树木结构将所述Mesh三角网格面模型划分为若干个封闭的三角网格体，所述三角网格体用于区分出所述目标树木的一种树木结构的整体，所述树木结构至少包括主干、二级枝条、三级枝条；
所述根据在所述Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系模拟出所述目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积，包括：
根据所述三角网格体区分出来的所述树木结构的整体以及在所述Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系，模拟出所述目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积。
4.根据权利要求3所述的方法，其特征在于，所述根据在所述Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系模拟出所述目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积，包括：
将体积小于预设的有效枝条体积阈值的所述三角网格体进行删除，得到所述目标树木的有效枝条；
在所述有效枝条中划分出所述目标树木的主干并计算出所述主干的体积；
对于划分出主干后的所述有效枝条，判断每个所述有效枝条是否与所述主干相交；
根据枝条间的相交关系，将判断结果为“是”的所述有效枝条判定为二级枝条并计算出所述二级枝条的体积；
根据在所述Rhinoceros 3D建模软件中预先设置的二级枝条对应的修剪条件以及修剪参数对所述二级枝条进行修剪，得到二级枝条修剪后剩余部分的三角网格体并计算所述二级枝条修剪后剩余部分的体积；
根据在所述Rhinoceros 3D建模软件中预先设置的主干对应的修剪条件以及修剪参数对所述主干进行修剪，得到主干修剪后剩余部分的三角网格体并计算所述主干修剪后剩余部分的体积；
对于划分出主干以及二级枝条后的所述有效枝条，判断每个所述有效枝条是否与所述二级枝条相交；
根据枝条间的相交关系，将判断结果为“是”的所述有效枝条判定为三级枝条并计算出所述三级枝条的体积；
根据在所述Rhinoceros 3D建模软件中预先设置的三级枝条对应的修剪条件以及修剪参数对所述三级枝条进行修剪，得到三级枝条修剪后剩余部分的三角网格体并计算所述三级枝条修剪后剩余部分的体积；
根据所述主干修剪后剩余部分、所述二级枝条修剪后剩余部分、所述三级枝条修剪后剩余部分的三角网格体模拟出所述目标树木的修剪后的形态。
5.根据权利要求4所述的方法，其特征在于，在所述根据在所述Rhinoceros 3D建模软件中预先设置的二级枝条对应的修剪条件以及修剪参数对所述二级枝条进行修剪之前，还包括：
计算有效二级枝条体积阈值，所述有效二级枝条体积阈值为预设的有效二级枝条系数与体积最大的所述二级枝条的体积的乘积；
判断每个所述二级枝条的体积是否小于所述有效二级枝条体积阈值；
将判断结果为“是”的所述二级枝条进行删除。
6.根据权利要求4所述的方法，其特征在于，所述根据在所述Rhinoceros 3D建模软件中预先设置的三级枝条对应的修剪条件以及修剪参数对所述三级枝条进行修剪之前，还包括：
计算每一个三级枝条对应的有效三级枝条体积阈值，所述三级枝条对应的有效三级枝条体积阈值为预设的有效三级枝条系数与和所述三级枝条相连接的二级枝条的体积的乘积；
判断每个所述三级枝条的体积是否小于所述三级枝条对应的有效三级枝条体积阈值；
将判断结果为“是”的所述三级枝条进行删除。
7.根据权利要求1所述的方法，其特征在于，所述碳排放的计算公式如下：
其中，C表示所述目标区域内树木群落的修剪废弃物产生的碳排放量，V表示所述目标区域内树木群落的预测修剪体积，ρ表示被修剪的枝条的平均密度，i表示被修剪的枝条的含碳量系数，EF表示被修剪枝条的处理方式的碳排放因子。
8.一种树木修剪废弃物碳排放预测系统，其特征在于，包括：
数据获取模块，用于获取由地面激光雷达扫描仪采集到的目标区域的树木群落点云数据；
三维结构模型生成模块，用于对所述点云数据进行点云分割，得到目标树木的点云文本文件，并利用定量结构模型将所述点云文本文件转换成3D模型文件格式的所述目标树木的三维结构模型；
单棵树木修剪体积计算模块，利用Rhinoceros 3D建模软件读取所述目标树木的三维结构模型；根据在所述Rhinoceros 3D建模软件中预先设置的修剪条件、修剪参数以及枝条间的相交关系模拟出所述目标树木的修剪后的形态并分别计算出主干、二级枝条、三级枝条、主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积；将所述主干、二级枝条、三级枝条的体积之和作为所述目标树木的整树体积，将所述主干修剪后剩余部分、二级枝条修剪后剩余部分、三级枝条修剪后剩余部分的体积之和作为所述目标树木的修剪后的体积；将所述目标树木的整树体积与修剪后的体积作差，得到所述目标树木的修剪体积；
目标区域树木修剪体积之和计算模块，用于以单棵树木为单位，遍历所述目标区域，对所述目标区域内的每棵树木的修剪体积进行计算并求和，得到所述目标区域内树木群落的预测修剪体积；
树木修剪废弃物碳排放预测模块，用于根据所述目标区域内树木群落的预测修剪体积计算所述目标区域内树木群落的修剪废弃物产生的碳排放。
9.一种电子设备，其特征在于，包括：处理器、存储器、通信接口和通信总线，所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信；
所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行如权利要求1-7中任一项所述的方法对应的操作。
10.一种计算机存储介质，其特征在于，所述计算机存储介质上存储有计算机程序，该程序被处理器执行时实现如权利要求1-7中任一所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310700457.2A CN116432873B (zh) | 2023-06-14 | 2023-06-14 | 树木修剪废弃物碳排放预测方法、系统、电子设备和介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310700457.2A CN116432873B (zh) | 2023-06-14 | 2023-06-14 | 树木修剪废弃物碳排放预测方法、系统、电子设备和介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116432873A CN116432873A (zh) | 2023-07-14 |
| CN116432873B true CN116432873B (zh) | 2023-09-19 |
Family
ID=87080137
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310700457.2A Active CN116432873B (zh) | 2023-06-14 | 2023-06-14 | 树木修剪废弃物碳排放预测方法、系统、电子设备和介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116432873B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105608739A (zh) * | 2016-03-16 | 2016-05-25 | 福州大学 | 数据和规则综合驱动的三维树木精细建模方法 |
| KR20200022262A (ko) * | 2018-08-22 | 2020-03-03 | 고려대학교 산학협력단 | 국가 산림탄소 인벤토리를 위한 산림탄소를 모델링하는 방법 |
| CN114136208A (zh) * | 2021-11-22 | 2022-03-04 | 北京林业大学 | 一种无损估算立木材积的低成本树木结构自动重建方法 |
| CN115659815A (zh) * | 2022-10-31 | 2023-01-31 | 江苏中发建筑设计有限公司 | 老旧建筑拆除的碳排放数据处理方法、系统、设备及介质 |
-
2023
- 2023-06-14 CN CN202310700457.2A patent/CN116432873B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105608739A (zh) * | 2016-03-16 | 2016-05-25 | 福州大学 | 数据和规则综合驱动的三维树木精细建模方法 |
| KR20200022262A (ko) * | 2018-08-22 | 2020-03-03 | 고려대학교 산학협력단 | 국가 산림탄소 인벤토리를 위한 산림탄소를 모델링하는 방법 |
| CN114136208A (zh) * | 2021-11-22 | 2022-03-04 | 北京林业大学 | 一种无损估算立木材积的低成本树木结构自动重建方法 |
| CN115659815A (zh) * | 2022-10-31 | 2023-01-31 | 江苏中发建筑设计有限公司 | 老旧建筑拆除的碳排放数据处理方法、系统、设备及介质 |
Non-Patent Citations (1)
| Title |
|---|
| 果树虚拟修剪研究进展;吕萌萌 等;系统仿真学报(第03期);14-26 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116432873A (zh) | 2023-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106815847B (zh) | 基于激光雷达点云的树木分割方法及单棵树提取方法 | |
| CN106529469B (zh) | 基于自适应坡度的无人机载LiDAR点云滤波方法 | |
| CN110363861B (zh) | 基于激光雷达点云的田地作物三维重构方法 | |
| CN111985496B (zh) | 一种基于植被三维凸包的树障隐患快速检测方法和终端 | |
| CN108038908B (zh) | 基于人工智能的空间对象识别及建模方法和系统 | |
| CN110223314A (zh) | 一种基于树冠三维点云分布的单木分割方法 | |
| Sithole et al. | Report: ISPRS comparison of filters | |
| CN107016725A (zh) | 一种顾及LiDAR点云数据分布差异的植被三维实景建模方法 | |
| CN104298998A (zh) | 一种3d点云的数据处理方法 | |
| CN102708587A (zh) | 一种快速获取三维建筑信息的方法及系统 | |
| CN110765962A (zh) | 一种基于三维点云轮廓分维值的植物识别分类方法 | |
| Zhao et al. | A comparison of LiDAR filtering algorithms in vegetated mountain areas | |
| CN110322497A (zh) | 一种基于三维可视化的交互式点云对象提取方法 | |
| CN115600307B (zh) | 一种从城市场景的Mesh模型中生成单体化建筑的方法 | |
| Hu et al. | A fast and simple method of building detection from LiDAR data based on scan line analysis | |
| Özdemir et al. | Automatic extraction of trees by using multiple return properties of the lidar point cloud | |
| CN115830234A (zh) | 一种输电线路建模用点云处理方法及处理系统 | |
| CN116432873B (zh) | 树木修剪废弃物碳排放预测方法、系统、电子设备和介质 | |
| Sarıtaş et al. | Enhancing Ground Point Extraction in Airborne LiDAR Point Cloud Data Using the CSF Filter Algorithm | |
| CN110033459B (zh) | 顾及地物完整性的大规模点云快速分块方法 | |
| Korzeniowska et al. | Generating DEM from LiDAR data–comparison of available software tools | |
| CN115410036A (zh) | 一种高压架空输电线路关键要素激光点云自动分类方法 | |
| CN116030350A (zh) | 一种基于点云技术的城市绿地植被碳汇量计算方法及系统 | |
| CN112381029B (zh) | 一种基于欧氏距离的机载LiDAR数据建筑物提取方法 | |
| CN115655098A (zh) | 高密度激光点云技术用于电网工程挖填土石方测算方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |