CN102538673B - 一种机载小光斑LiDAR波形数据提取林分垂直结构参数的方法 - Google Patents
一种机载小光斑LiDAR波形数据提取林分垂直结构参数的方法 Download PDFInfo
- Publication number
- CN102538673B CN102538673B CN201110424549.XA CN201110424549A CN102538673B CN 102538673 B CN102538673 B CN 102538673B CN 201110424549 A CN201110424549 A CN 201110424549A CN 102538673 B CN102538673 B CN 102538673B
- Authority
- CN
- China
- Prior art keywords
- forest
- waveform data
- vertical structure
- waveform
- lidar
- 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
Images
Landscapes
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明的目的在于提供一种机载小光斑LiDAR波形数据提取林分垂直结构参数的方法,包括:基于机载小光斑LiDAR波形数据生成激光点云数据集;对该激光点云数据集进行地面点与植被点的滤波处理,利用滤波处理后的地面点生成地面高程模型,将所述机载小光斑LiDAR波形数据中的高程减去相应位置的DEM值得到高程归一化后的波形数据;对波形数据进行林分尺度的波形合成处理,生成林分的合成波形数据;将合成波形数据分解为一系列子波形的组合;根据所述一系列子波形的组合,获取林分垂直结构参数。本发明提供的方法,比使用点云数据对森林垂直结构的刻画更加准确,并可以确定林分垂直结构参数。
Description
技术领域
本发明涉及激光雷达数据处理和信息提取技术,尤其涉及一种机载小光斑LiDAR波形数据提取林分垂直结构参数的方法。
背景技术
激光雷达(Light Detection And Ranging,LiDAR)是一项通过由传感器所发出的激光脉冲来测定传感器与目标物之间距离的主动遥感技术。它可以应用于林业研究,尤其可以用于遥感森林垂直结构并估测林木高度。
随着林分年龄的增大和林木的生长,树冠组分的垂直分布相对于幼林时的状态也发生变化。成过熟天然林的特点是林隙多、林龄结构和高度变化大,相对同龄林而言在垂直结构上表现得更复杂,而同龄林树顶附近的冠层物质占了绝大部分。激光雷达的波形对这些林分连续性垂直结构的变化可以很好地刻画。现有的激光雷达反演林分垂直结构可以通过分析小光斑LiDAR的点云数据来实现。但由于机载小光斑LiDAR数据仅能获取少数离散的点云,对林分下层的垂直结构很难得到很好的刻画,对林分垂直结构的刻画精度较低。
发明内容
本发明的目的是提供一种机载小光斑LiDAR波形数据提取林分垂直结构参数的方法,从而解决上述技术问题。
本发明的目的通过以下技术方案来实现。
一种机载小光斑LiDAR波形数据提取林分垂直结构参数的方法,包括:
步骤1、基于机载小光斑LiDAR波形数据生成激光点云数据集;
步骤2、对所述的激光点云数据集进行地面点与植被点的滤波处理,利用滤波处理后的地面点生成地面高程模型,将所述机载小光斑LiDAR波形数据中的高程减去相应位置的DEM值得到高程归一化后的波形数据;
步骤3、对所述高程归一化后的波形数据进行林分尺度的波形合成处理,得到林分尺度上的合成波形数据;
步骤4、对所述的林分尺度上的合成波形数据进行处理,分解为子波形的集合;
步骤5、根据所述的子波形的集合,获取林分垂直结构信息。
本发明提供的方法,比使用点云数据对林分垂直结构的刻画更加准确,可以确定林分的垂直结构参数。
附图说明
图1本发明实施例提供的方法流程图。
具体实施方式
本发明提供了一种机载小光斑LiDAR波形数据提取林分垂直结构参数的方法,其实现方式如图1所示,具体包括如下操作:
步骤1、基于机载小光斑LiDAR波形数据生成激光点云数据集;
对飞行过程中获取的机载小光斑LiDAR波形数据通过峰值检测、高斯分解或反卷积方法可以提取每个波形数据对应的反映地物特征点的文件,这些反映地物特征点的文件构成激光点云数据集。
步骤2、对所述的激光点云数据集进行地面点与植被点的滤波处理,利用滤波处理后的地面点生成地面高程模型,将所述机载小光斑LiDAR波形数据中的高程减去相应位置的地面高程模型值得到高程归一化后的波形数据;
本发明中,地面高程模型是指DEM(Digital Elevation Model,数字高程模型)。
步骤3、对所述高程归一化后的波形数据进行林分尺度的波形合成处理,得到林分尺度上的合成波形数据;
步骤3的具体实现方式可以是:基于林分尺度对高程归一化后的波形数据进行抽取、背景噪声去除、垂直剖面构建和波形合成,最终得到林分尺度上的合成波形数据。
步骤4、对所述的林分尺度上的合成波形数据进行处理,分解为子波形的集合。
步骤4的实现方式可以但不仅限于以下两种:
(一)对所述的林分尺度上的合成波形数据通过高斯分解处理分解为子波形的集合,具体实现方式是,通过高斯分解,将所述的林分尺度上的合成波形数据分解成为一系列子波形,每个子波形代表不同垂直分层目标(其中,垂直分层目标是指在一定高度上植被组分的集合)的回波信号,可以用高斯分量和噪声之和的形式表示,其公式表示如下:
其中,f(x)为一个林分内合成的回波信号;b为噪声,n为所述的林分内垂直分层目标的数量;fj(x)表示第j个垂直分层位置目标的回波信号;x为可能的植物组分。
通过对林分尺度上的合成波形数据进行高斯分解,得到一系列的子波形,以便后续步骤中,可以通过子波形的个数来确定回波信号的数量,从而获取垂直分层目标的特征信息。所述的目标特征信息是指林层位置和林层高度,林层位置和林层高度通过高斯分解得到的子波形的振幅和展宽体现。
(二)对所述的林分尺度上的合成波形数据通过反卷积处理分解为子波形的集合。
激光雷达的发射脉冲与目标物(即垂直分层目标)的作用过程相当于是一个卷积的过程,也就是说,激光雷达接收到的回波信号是发射脉冲和目标物响应的卷积。在已知发射脉冲及回波信号的情况下,要了解目标物的特征就相当于求解反卷积的过程。因此,本发明实施例中,可以对所述的林分尺度上的合成波形数据进行反卷积处理,以便后续步骤中获取垂直分层目标的特征信息,用公式表示如下:
f(x)=s(x)*h(x)+n(x) 公式2
其中,f(x)为一个林分的合成回波信号;s(x)表示激光雷达的发射脉冲;n(x)为噪声,h(x)则为所需要的垂直分层目标的特征信息。
步骤5、根据所述的子波形的集合,获取林分垂直结构信息。
上述林分尺度上的合成波形数据通过高斯分解或反卷积等方法的处理,被分解为一系列子波形的集合。其中,子波形的个数代表林层数量,一个子波形代表了林分内一个林层的回波信号,通过子波形坐标的解算可以获取林层位置信息,波形宽度信息可以用于提取林层高度信息,从而获取林分的垂直结构信息(即垂直分层目标的特征信息)。
本发明实施例提供的技术方案,直接将机载小光斑LiDAR波形数据在林分尺度上进行波形合成,然后使用高斯分解和反卷积从合成的波形数据中提取林分垂直结构的参数,获取各林层信息(林层数量、林层位置、林层高度等参数)的方法。通过本发明,可以在传统的激光雷达测绘应用生成的DEM的基础上,自动获得林分的垂直结构信息,为森林结构监测、森林质量评价、森林经营决策、生物多样性评价等方面提供了可靠的遥感估测手段。随着机载小光斑LiDAR波形数据日益广泛地获取,本发明具有很强的实用价值。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发发明的保护范围应该以权利要求书的保护范围为准。
Claims (4)
1.一种机载小光斑LiDAR波形数据提取林分垂直结构参数的方法,其特征在于,包括:
步骤1、基于机载小光斑LiDAR波形数据生成激光点云数据集;
步骤2、对所述的激光点云数据集进行地面点与植被点的滤波处理,利用滤波处理后的地面点生成地面高程模型,将所述机载小光斑LiDAR波形数据中的高程减去相应位置的DEM值得到高程归一化后的波形数据;
步骤3、基于林分尺度对高程归一化后的波形数据进行抽取、背景噪声去除、垂直剖面结构和波形合成,得到林分尺度上的合成波形数据;
步骤4、将所述合成波形数据通过处理分解为一系列子波形的集合;
步骤5、根据所述的子波形的集合,获取林分垂直结构信息。
2.根据权利要求1所述的方法,其特征在于,所述步骤4具体包括:
对所述的林分尺度上的合成波形数据进行高斯分解,得到一系列子波形的集合,每个子波形代表不同垂直分层目标的回波信号。
3.根据权利要求1所述的方法,其特征在于,所述步骤4具体包括:
对所述的林分尺度上的合成波形数据进行反卷积处理,得到一系列子波形的集合。
4.根据权利要求2或3所述的方法,其特征在于,所述林分垂直结构信息包括:林层数量、林层位置、林层高度;所述步骤5具体包括:
获取子波形的个数,所述子波形的个数代表林层数量;
通过对子波形坐标进行解算获取林层位置;
获取子波形的波形宽度信息,从所述波形宽度信息中获取林层高度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110424549.XA CN102538673B (zh) | 2011-12-19 | 2011-12-19 | 一种机载小光斑LiDAR波形数据提取林分垂直结构参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110424549.XA CN102538673B (zh) | 2011-12-19 | 2011-12-19 | 一种机载小光斑LiDAR波形数据提取林分垂直结构参数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102538673A CN102538673A (zh) | 2012-07-04 |
CN102538673B true CN102538673B (zh) | 2014-04-09 |
Family
ID=46346182
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110424549.XA Active CN102538673B (zh) | 2011-12-19 | 2011-12-19 | 一种机载小光斑LiDAR波形数据提取林分垂直结构参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102538673B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104867180B (zh) * | 2015-05-28 | 2017-09-15 | 南京林业大学 | 一种集成uav与lidar的森林林分特征反演方法 |
CN105676205A (zh) * | 2016-01-27 | 2016-06-15 | 武汉大学 | 一种机载LiDAR波形数据高斯分解方法 |
CN106199627B (zh) * | 2016-09-14 | 2018-07-10 | 中国农业科学院农业资源与农业区划研究所 | 一种无人机lidar反演草地植被参数的精度改进方法 |
CN112241661A (zh) * | 2019-07-17 | 2021-01-19 | 临沂大学 | 一种结合机载LiDAR点云数据和航空影像的城市地物精细化分类方法 |
CN111077532A (zh) * | 2019-11-22 | 2020-04-28 | 同济大学 | 一种基于反卷积和高斯分解的地物空间信息获取方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FI112402B (fi) * | 1999-10-28 | 2003-11-28 | Diware Oy | Menetelmä puustotunnusten määrittämiseksi sekä tietokoneohjelma menetelmän suorittamiseksi |
CN101196562B (zh) * | 2007-12-14 | 2012-01-04 | 武汉大学 | 一种基于改进的em算法的激光雷达波形数据分解的方法 |
US8577611B2 (en) * | 2010-03-30 | 2013-11-05 | Weyerhaeuser Nr Company | System and method for analyzing trees in LiDAR data using views |
CN102103202B (zh) * | 2010-12-01 | 2013-01-09 | 武汉大学 | 一种融合影像的机载激光雷达数据半监督分类方法 |
-
2011
- 2011-12-19 CN CN201110424549.XA patent/CN102538673B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN102538673A (zh) | 2012-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104656098B (zh) | 一种遥感森林生物量反演的方法 | |
CN109031344B (zh) | 一种全波形激光雷达和高光谱数据联合反演森林结构参数的方法 | |
CN102538673B (zh) | 一种机载小光斑LiDAR波形数据提取林分垂直结构参数的方法 | |
CN108921885B (zh) | 一种综合三类数据源联合反演森林地上生物量的方法 | |
Ørka et al. | Effects of different sensors and leaf-on and leaf-off canopy conditions on echo distributions and individual tree properties derived from airborne laser scanning | |
Sumnall et al. | Comparison of small-footprint discrete return and full waveform airborne lidar data for estimating multiple forest variables | |
CN113204998B (zh) | 一种基于单木尺度的机载点云森林生态估测方法及*** | |
Razak et al. | Airborne laser scanning of forested landslides characterization: Terrain model quality and visualization | |
CN106990401B (zh) | 基于全波形机载激光雷达数据二类高程误差修正方法 | |
KR20200122897A (ko) | 라이다 측량자료를 이용한 숲틈 모니터링 시스템 및 방법 | |
CN104808191B (zh) | 一种基于全波形LiDAR单木冠层容积分解的树种分类方法 | |
CN107679476B (zh) | 一种海冰类型遥感分类方法 | |
CN104020475A (zh) | 一种基于机载LiDAR数据的电力线提取及建模方法 | |
CN103398957A (zh) | 基于高光谱与激光雷达提取叶面积垂直分布的方法 | |
CN104820830A (zh) | 一种基于全波形LiDAR冠层剖面模型的树种识别方法 | |
Rosette et al. | Lidar remote sensing for biomass assessment | |
CN104502919A (zh) | 利用机载激光雷达点云提取城市植被三维覆盖的方法 | |
CN108872964A (zh) | 基于无人机LiDAR数据的银杏人工林郁闭度提取方法 | |
CN108959705A (zh) | 一种预测大面积亚热带森林生物量的混合效应模型 | |
Chen et al. | Site quality assessment of a Pinus radiata plantation in Victoria, Australia, using LiDAR technology | |
CN105866792A (zh) | 一种新的星载激光雷达树高提取方法 | |
Wack et al. | Forest inventory for eucalyptus plantations based on airborne laserscanner data | |
Feng et al. | PyFLEXTRKR: A flexible feature tracking Python software for convective cloud analysis | |
CN101833087B (zh) | 基于小波分析的星载激光雷达波形数据分解的方法 | |
Reitberger et al. | Benefit of airborne full waveform lidar for 3D segmentation and classification of single trees |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |