CN114241155A - 一种基于车载激光点云数据的城市树木三维可视化方法 - Google Patents
一种基于车载激光点云数据的城市树木三维可视化方法 Download PDFInfo
- Publication number
- CN114241155A CN114241155A CN202111472677.1A CN202111472677A CN114241155A CN 114241155 A CN114241155 A CN 114241155A CN 202111472677 A CN202111472677 A CN 202111472677A CN 114241155 A CN114241155 A CN 114241155A
- Authority
- CN
- China
- Prior art keywords
- tree
- points
- point
- point cloud
- crown
- 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.)
- Pending
Links
- 238000007794 visualization technique Methods 0.000 title claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 15
- 230000011218 segmentation Effects 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000000605 extraction Methods 0.000 claims abstract description 4
- 238000012800 visualization Methods 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 9
- 238000010606 normalization Methods 0.000 claims description 4
- 239000002023 wood Substances 0.000 claims description 4
- 230000000877 morphologic effect Effects 0.000 abstract description 5
- 238000000926 separation method Methods 0.000 abstract description 4
- 230000007547 defect Effects 0.000 abstract description 3
- 238000011835 investigation Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
Images
Classifications
-
- 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
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明属于车载激光雷达点云数据处理技术领域，具体涉及一种基于车载激光点云数据的城市树木三维可视化方法。基于原始点云数据，经树木分割、树冠树干分离后，对树干和树冠单独建模，最后构建城市场景下的树木三维模型。相比于传统的树木三维建模方法，能够根据树干和树冠的不同形态特征进行单独建模，避免了单树建模形态特征表达不精确的缺点，能够简单、快速、准确地构建出城市树木三维模型，可以应用于林业资源调查、3D城市建模，树木参数提取等诸多领域。
Description
技术领域
本发明属于车载激光雷达点云数据处理技术领域，具体涉及一种基于车载激光点云数据的城市树木三维可视化方法。
背景技术
激光雷达是一种主动遥感设备，突破了传统的单点测量方法，通过主动发射激光脉冲测定传感器与目标物体之间的距离。具有高精度、高分辨率、高效率的独特优势，特别能穿透地表覆盖的森林植被快速获取地形信息，为森林结构参数的提取以及树木三维可视化提供了全新的技术手段。
树木等植被模型是城市场景建模中不可或缺的重要组成部分。同时，三维树木模型也是虚拟场景中极具代表性的自然景观模型，已被广泛应用于虚拟旅游、虚拟城市、虚拟生态景观等方面。三维点云数据可以完美地呈现物体表面结构信息、空间信息以及物体的模型细节，能够有效避免传统的基于图像的三维重建无法直接提供空间信息的问题。构建三维树木模型可以定量提取树高、胸径等植被参数，同时也是推动林业数字化进程的基础和关键。
目前常用的点云三维重建方法主要有贪婪投影三角化算法、移动立方体以及泊松曲面重建算法。贪婪投影三角化算法首先通过法线将点云投影到某一平面，在二维平面上建立点云间的拓扑关系，最后得到三维的三角网格曲面模型，但该方法不能同时进行曲面平滑和孔洞修复；移动立方体算法是在三维离散数据场中通过线性插值来逼近等值面，但重建后会包含大量的三角形面片数据且信息丢失严重；泊松重建通过平滑滤波指示函数，构建柏松方程，建立表面滤波后的指示函数梯度场和表面向内法向量场的关系，求解等值面，但不能很好捕捉重建表面的局部细节。
发明内容
针对上述存在问题或不足，本发明提供了一种基于车载激光点云数据的树木三维可视化方法，基于原始点云数据，经树木分割、树冠树干分离后，对树干和树冠单独建模，最后构建城市场景下的树木三维场景。
本发明具体技术方案如下：
步骤1、去噪与滤波：
由于激光雷达扫描设备的精度以及周围环境的影响，点云数据中经常包含噪声点，因此第一步工作是去除点云中的噪声点。然后针对去噪后的点云数据进行数据滤波，得到地面点与非地面点数据。
步骤2、单木分割与树干树冠提取：
将滤波后的数据进行归一化处理。首先使用滤波得到的地面点生成数字高程模型(DEM，Digital Elevation Model)，之后通过非地面点高程减去DEM高度得到冠层高度模型(CHM，Canopy Height Model)。使用基于CHM单木分割算法提取出单木。之后利用圆柱拟合算法从单木中分离树干和树冠。
步骤3、表面重构：
使用alpha-shape算法对树干和树冠分别构建三维模型，步骤如下：
3-1、设置距离阈值α
3-2、构建待建模点集Q＝{qi}，(i＝1,2,…n)，从中任选一个点q1，与之距离小于2α的点构成新的点集Q1，从点集Q1中任取一组点q2和q3，求出过点q1、q2、q3且半径为α的球以及球心o1、o2；
3-3、遍历点集Q1，依次求出其他点到球心o1、o2的距离集合L1、L2。如果L1和L2中有一个集合的距离均大于α，则q1、q2、q3为边缘轮廓点，连接三点构成边界三角面片；反之，则不是边缘轮廓点，停止遍历，执行下一步；
3-4、选择点集Q1中的下一组点按步骤(2)、(3)进行判断，直到点集Q1中的所有点判断结束。输出三角面片集合，集合中暴露的三角面片形成局部凸包；
3-5、选择点集Q中下一个点按步骤(2)、(3)、(4)进行判断，直到点集Q中的所有点判断结束，多个凸包重构表面S(Q)。
步骤4、三维可视化：
将单株树木对应的树干树冠模型组合，完成城市树木三维可视化。
步骤3所涉及的原理：
alpha-shape的原理：
假设有一点集P，其alpha-shape是一个由点集P和半径参数α决定的唯一多边形。可设想成存在一个半径为α的圆，从点集P外某一点开始滚动。若α足够大，那么这个圆就不会滚到点集P的内部，则这个圆滚动的痕迹就是点集P的边界线。若α很大，则检测出的边界线是点集P的凸包。同理，三维空间的alpha-shape算法则是通过三点作半径为α的球来判断边界点，并在得到的边界点处建立三角面片，重构出曲面。
本发明首先对点云进行去噪和滤波处理，从原始车载激光雷达点云中分离出地面点和非地面点；再利用非地面点得到的数字高程模型对非地面点进行高程归一化处理，得到冠层高度模型；再使用基于CHM的单木分割算法提取单株木，使用圆柱拟合算法分离单株木的树冠和树干；然后使用alpha-shape算法将树冠树干分别建模，最后得到城市场景下的树木三维可视化模型，其流程如图1所示。
本发明是通过车载激光雷达获取城市树木三维激光点云数据，并根据树干和树冠点云不同的三维分布特征，经过去噪与滤波、点云归一化、单木分割、树干树冠分离，完成树木三维可视化，建立基于原始激光点云从城市场景中提取树木并进行三维可视化的方法。
综上所述，本发明相比于传统的树木三维建模方法，能够根据树干和树冠的不同形态特征进行单独建模，避免了单树建模形态特征表达不精确的缺点，能够简单、快速、准确地构建出城市树木三维模型。
附图说明
图1为本发明的流程示意图；
图2为单木分割结果示意图；
图3为树干、树冠分离结果图；a.plot1树干集结果图；b.plot2树冠集结果图；
图4为alpha-shape原理图；
图5为树木三维可视化结果图；a.plot1单株木结果图；b.plot2所有树木结果图。
具体实施方式
以下通过一个样方的实施并结合附图对本发明作进一步详细说明。本实例开发环境Microsoft Visual Studio 2019、外部库PCL1.12.0、编程语言C++。
一种基于车载激光点云数据的城市树木三维可视化方法，包括以下步骤：
步骤1，以某校车载激光点云数据为例开展城市树木三维可视化表达，数据使用SSW-2型车载激光扫描仪获取，点云中包括车载系统获取的所有地物信息，如地面、行道树、建筑物、绿化带、基础公共设施等，其中行道树包括独立行道树、相连行道树，点云密度为平均201点/米2。使用cloudcompare软件对获取的点云数据进行去噪处理，使用CSF滤波算法将去噪后的点云分为地面点与非地面点。
步骤2，使用非地面点数据插值生成DEM模型，分辨率设置为0.1m。非地面点云高程值减去所对应的DEM高程值得到CHM模型，使用基于CHM模型的单木分割算法得到249棵单木。使用基于圆柱拟合的树干提取算法，分离出树干和树冠。圆环高度h设置为0.1m，圆环拟合误差阈值为0.25，半径变化百分比阈值为0.3，移动步长设置为0.1m。
步骤3，从树冠集中，每次取出单一树木的树冠进行表面重建，α值设为0.35，rgb值设为0,90,0；从树干集中，每次取出单一树木的树干，α值设为0.1，rgb值设为140,90,50，拼接得到单株树木的三维模型。对所有单株树木的树冠和树干建模完成后，得到城市场景下的树木三维可视化模型。
步骤4，精度检验。经过以上几步，最终得到树木三维模型如图所示。从模型精度方面进行检验，如图所示，和三维点云对比可以看出，该方法准确地从原始点云中提取出了单木点云并分离出树干和树冠，树木三维模型较好地还原了树木的形态特征。但是由于点云密度较低，对树冠内部更精细程度的建模还存在一定的优化空间。
该实例根据本发明提出的方法，对样地地面激光雷达数据进行处理，依照方法步骤，依次对点云进行了去噪、滤波、单木分割以及树干树冠分离，都取得了较好地精度，说明该方法在城市场景下的树木三维可视化方面具有较好的效果。
Claims (2)
1.一种基于车载激光点云数据的城市树木三维可视化方法，包括以下步骤：
步骤1、去噪与滤波：
将场景点云中的噪声剔除并分离地面点与非地面点；
步骤2、单木分割与树干树冠提取：
将滤波后的数据进行归一化处理，得到冠层高度模型(CHM，Canopy Height Model)。使用基于CHM单木分割算法提取出单木。之后利用基于圆柱拟合的方法从单木中分离树干和树冠。
步骤3、表面重构：
使用alpha-shape算法对树干和树冠分别构建三维模型，步骤如下：
3-1、设置距离阈值α
3-2、构建待建模点集Q＝{qi}，(i＝1,2,…n)，从中任选一个点q1，与之距离小于2α的点构成新的点集Q1，从点集Q1中任取一组点q2和q3，求出过点q1、q2、q3且半径为α的球以及球心o1、o2；
3-3、遍历点集Q1，依次求出其他点到球心o1、o2的距离集合L1、L2。如果L1和L2中有一个集合的距离均大于α，则q1、q2、q3为边缘轮廓点，连接三点构成边界三角面片；反之，则不是边缘轮廓点，停止遍历，执行下一步；
3-4、选择点集Q1中的下一组点按步骤(2)、(3)进行判断，直到点集Q1中的所有点判断结束。输出三角面片集合，集合中暴露的三角面片形成局部凸包；
3-5、选择点集Q中下一个点按步骤(2)、(3)、(4)进行判断，直到点集Q中的所有点判断结束，多个凸包重构表面S(Q)。
步骤4、三维可视化：
将单株木对应的树干树冠模型组合，完成城市树木三维可视化。
2.如权利要求1所述基于车载激光点云数据的城市树木三维可视化方法，其特征在于：所述步骤1中采用CSF算法分离场景点云中的地面点与非地面点。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111472677.1A CN114241155A (zh) | 2021-12-06 | 2021-12-06 | 一种基于车载激光点云数据的城市树木三维可视化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111472677.1A CN114241155A (zh) | 2021-12-06 | 2021-12-06 | 一种基于车载激光点云数据的城市树木三维可视化方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114241155A true CN114241155A (zh) | 2022-03-25 |
Family
ID=80753106
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111472677.1A Pending CN114241155A (zh) | 2021-12-06 | 2021-12-06 | 一种基于车载激光点云数据的城市树木三维可视化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114241155A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115937454A (zh) * | 2023-01-09 | 2023-04-07 | 航天宏图信息技术股份有限公司 | 一种在大规模城市场景中自动摆放树木模型的方法和装置 |
| CN116129391A (zh) * | 2023-04-18 | 2023-05-16 | 山东省国土测绘院 | 一种从车载激光点云提取行道树的方法及系统 |
-
2021
- 2021-12-06 CN CN202111472677.1A patent/CN114241155A/zh active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115937454A (zh) * | 2023-01-09 | 2023-04-07 | 航天宏图信息技术股份有限公司 | 一种在大规模城市场景中自动摆放树木模型的方法和装置 |
| CN116129391A (zh) * | 2023-04-18 | 2023-05-16 | 山东省国土测绘院 | 一种从车载激光点云提取行道树的方法及系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR101165534B1 (ko) | 수관 식물점 그룹에 대해 시뮬레이트된 나무 줄기 및 나무 가지를 제공하는 지리공간 모델링 시스템 | |
| CN113034689A (zh) | 基于激光点云的地形三维模型及地形图构建方法和系统、存储介质 | |
| CN106815847A (zh) | 基于激光雷达点云的树木分割方法及单棵树提取方法 | |
| Gobeawan et al. | Modeling trees for virtual Singapore: From data acquisition to CityGML models | |
| CN106529469A (zh) | 基于自适应坡度的无人机载LiDAR点云滤波方法 | |
| CN106503060A (zh) | 一种输电线路三维点云数据处理及交跨物获取方法 | |
| WO2015149302A1 (zh) | 基于点云与数据驱动的树木模型重建方法 | |
| CN114241155A (zh) | 一种基于车载激光点云数据的城市树木三维可视化方法 | |
| CN103258345A (zh) | 一种基于地面激光雷达三维扫描的树木枝干参数提取方法 | |
| CN107832849B (zh) | 一种基于知识库的电力线廊道三维信息提取方法及装置 | |
| CN111754618B (zh) | 一种面向对象的实景三维模型多层次解译方法及系统 | |
| Hui et al. | A mean shift segmentation morphological filter for airborne LiDAR DTM extraction under forest canopy | |
| CN110009745B (zh) | 根据平面基元和模型驱动对点云提取平面的方法 | |
| Chen et al. | A mathematical morphology-based multi-level filter of LiDAR data for generating DTMs | |
| CN113362359A (zh) | 融合高度和光谱信息的倾斜摄影数据建筑物自动提取方法 | |
| CN114723907A (zh) | 面向无人机倾斜摄影三维模型数据的水面区域重建方法 | |
| CN112906719A (zh) | 基于消费级深度相机的立木因子测量方法 | |
| CN115018982A (zh) | 一种基于地基激光雷达点云的树木数字孪生方法 | |
| CN110580468B (zh) | 一种基于影像匹配点云的单木结构参数提取的方法 | |
| JP2004093632A (ja) | 地形形状抽出方法、地形形状抽出システム及びプログラム | |
| Yue et al. | A method for extracting street trees from mobile LiDAR point clouds | |
| CN115861532A (zh) | 基于深度学习的植被地物模型重建方法与系统 | |
| CN113724381B (zh) | 一种基于高分遥感影像的动态三维场景快速重建方法 | |
| CN112509142B (zh) | 基于面向表型精确鉴定的豆株快速三维重建方法 | |
| Cárdenas et al. | Reconstruction of tree branching structures from UAV-LiDAR data |
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 |