CN105866792B - 一种星载激光雷达树高提取方法 - Google Patents

一种星载激光雷达树高提取方法 Download PDF

Info

Publication number
CN105866792B
CN105866792B CN201610378509.9A CN201610378509A CN105866792B CN 105866792 B CN105866792 B CN 105866792B CN 201610378509 A CN201610378509 A CN 201610378509A CN 105866792 B CN105866792 B CN 105866792B
Authority
CN
China
Prior art keywords
forest canopy
forest
height
canopy height
tree
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
Application number
CN201610378509.9A
Other languages
English (en)
Other versions
CN105866792A (zh
Inventor
王成
聂胜
习晓环
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Remote Sensing and Digital Earth of CAS
Original Assignee
Institute of Remote Sensing and Digital Earth of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Remote Sensing and Digital Earth of CAS filed Critical Institute of Remote Sensing and Digital Earth of CAS
Priority to CN201610378509.9A priority Critical patent/CN105866792B/zh
Publication of CN105866792A publication Critical patent/CN105866792A/zh
Application granted granted Critical
Publication of CN105866792B publication Critical patent/CN105866792B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88—Lidar systems specially adapted for specific applications

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

星载激光雷达GLAS的波形数据可以分解为两个混合广义高斯波形，分别代表激光光斑森林冠层部分和地面部分。利用这些地面部分波形的波峰位置精确确定地面的位置，然后根据GLAS波形开始位置可以确定森林冠层顶部的位置，最终得到森林冠层高度。但是森林冠层高度受到光斑形状、大小以及坡度、坡向的影响，提出一种新的树高几何物理改正模型。该模型充分考虑上述因素对树高的影响，并考虑发射脉冲宽度对树高的影响，最终精确提取森林冠层高度。该发明专利利用GLAS数据精确提取森林冠层高度，弥补了以往遥感方法在森林植被高度研究上的不足，提高了森林冠层高度估测精度，为全球碳循环和全球气候变化研究开辟了新的途径。

Description

一种星载激光雷达树高提取方法
所属技术领域
本发明专利是一种对地观测领域的技术，涉及一种利用星载激光雷达数据提取森林冠层高度的方法。该方法具有一定的普适性，是一项具有理论研究和实际应用价值的方法。
背景技术
“碳”不仅是环境问题，而且蕴含着更多的经济和政治问题。全球碳循环研究是气候变化的核心研究内容之一，碳储量的时空变化依然是世界各国需要解决的共性问题。搞清中国在全球碳循环中的碳源/碳汇问题，关系到我国未来能源政策的制定以及如何履行联合国气候变化框架公约。
森林高度是反映碳储量的重要参数，传统获取森林冠层高度的方法是实地测量，采用测高仪器对单木进行测量。该类仪器依据空间几何原理进行测量树高，操作简单、结果精度高，适用于单木森林冠层高度的精准测量。但是，大范围森林高度测量则需要投入巨大的人力和物力，并且在一些条件苛刻的区域，人员难以到达且不易获取大面积测量数据。随着遥感技术的发展，光学、微波技术成为估测森林植被高度的重要手段(Brown andSchroeder,1999)，但在测量精度上均存在一定的缺陷，激光雷达技术成为其中重要的补充手段。
激光雷达(Light Detection And Ranging，简称LiDAR)是一种集激光、全球定位系统(Global Positioning System，简称GPS)和惯性导航系统(Intertial NavigationSystem，简称INS)技术于一身的用于快速获取地面及地面目标三维信息的主动式观测系统。LiDAR技术与成像光谱、合成孔径雷达一起被列为对地观测系统(EOS)计划最核心的信息获取与处理技术，该技术能精确探测地物垂直方向上的信息。其中星载激光雷达(GLAS)数据已经被广泛应用于森林植被结构参数的反演(Sun et al.,2000&2008；Lefsky etal., 2010)。但是由于GLAS为大光斑(每个光斑为70m的椭圆)，而森林植被区域地形起伏较大，影像森林冠层高度估算精度，一些情况下导致地形和植被波形难以分离。
发明内容
本发明专利就是利用星载激光雷达数据(GLAS)提供的三维空间信息，提取森林在垂直方向上的高度差即森林冠层高度。为精确提取森林冠层高度，建立基于星载激光雷达的树高坡度改正模型。该模型充分考虑光斑形状、大小以及坡度、坡向对树高的耦合影响，建立适用性广的几何物理模型。该模型能够精确地估测森林冠层高度，为全球碳循环研究提供基础。
附图说明
下面结合附图和实例对本发明专利进一步说明。
图1是工作流程图
图2是GLAS建筑物高度计算示意图
具体实施方式
本发明专利总体思路为：利用GLAS数据提取粗略的森林冠层高度，建立树高的几何物理模型，该模型充分考虑光斑大小、朝向以及坡度、坡向因素对树高的影响，并分析其影响规律，最终得到精确的森林冠层高度(图1)。
森林冠层高度H通过计算激光雷达第一个冠层回波信号与最后一个回波(地面回波)峰值的差得到(图 2)。星载GLAS波形开始为系统噪声，冠层回波从超过噪声阈值开始以一定的时间间隔进行记录，第一个回波信号为信号的开始，即激光雷达信号接触冠层顶部并反射的开始同时也表示冠层到传感器的距离，最后一个回波的峰值点表示从地面到传感器的距离，二者之差即为森林冠层高度。
通常情况下，GLAS回波的起始位置(最高点)到最后一个波峰位置(地面)之间的距离可以代表光斑内森林冠层高度。但是，GLAS数据为大光斑数据(直径70m)，大光斑激光雷达回波信号是激光脉冲与植被冠层、地面共同作用的结果，其波形受地形起伏/坡度影响显著，主要表现为波形的展宽和重叠。随着坡度的增加，波形的展宽会增加，当坡度增加一定程度时，地面数据会与植被的波形数据混合，大大增加了树高提取的难度。在这种情况下，为了精确提取树高，必须要能够精确地确定地面位置和树高改正模型。
为了精确确定地面位置，本发明专利首先提出了地面位置确定方法，通过波形模型得知，GLAS波形分为两个部分：植被部分和地面部分。这两个部分都能用广义高斯函数表示，地面部分的广义高斯函数波峰所在位置可以确定为地面位置。其具体表示为：
其中μground表示为地面位置。
另外，本发明专利还建立一种新的树高的几何物理改正模型。该模型充分考虑光斑大小、朝向以及坡度坡向对树高的影响，并在此基础上分析上述因素的影响规律。此外，还需要考虑光斑发射脉冲宽度对树高的影响，最终精确提取森林冠层高度。树高改正模型表示为式(2)，精确提取表达为式(3)。
H＝Wfstart-FWHM/2-GroundExtent/2-μground (3)
其中a表示为光斑长半轴、b为光斑的短半轴，θ是光斑朝向和坡向的夹角，η是地形坡度， GroundExtent是树高改正高度，Wfstart表示森林冠层高度顶部，FWHM为发射脉冲宽度，H为改正后的森林冠层高度。

Claims (3)

1.基于星载激光雷达数据精确提取森林树高方法，其特征主要包括以下几个步骤：
(1)将GLAS波形数据分解为两个混合广义高斯函数，确定地面位置；
(2)利用GLA01数据确定森林冠层顶部位置，结合地面位置得到森林冠层高度；
(3)建立树高几何物理改正模型，该模型充分考虑光斑大小、朝向以及坡度、坡向等对树高的影响；
(4)结合树高几何改正模型，得到精确的森林冠层高度；
树高改正模型为：
精确提取的森林冠层高度为：
H＝Wfstart-FWHM/2-GroundExtent/2-μground；
其中，a表示为光斑长半轴、b为光斑的短半轴，θ是光斑朝向和坡向的夹角，η是地形坡度，GroundExtent是树高改正高度，Wfstart表示森林冠层高度顶部，FWHM为发射脉冲宽度，H为改正后的森林冠层高度；
所述森林冠层高度通过计算激光雷达第一个冠层回波信号与最后一个回波峰值的差得到。
2.根据权利要求1所述的方法，其特征在于：利用星载激光雷达GLA01、GLA05、GLA14数据并利用广义高斯分解计算出粗略的森林冠层高度。
3.根据权利要求1所述的方法，其特征在于：建立树高的几何物理模型，分析光斑大小、朝向以及坡度坡向因素对树高的影响规律，并考虑发射脉冲宽度对树高的影响，最终得到精确的森林冠层高度。
CN201610378509.9A 2016-05-31 2016-05-31 一种星载激光雷达树高提取方法 Active CN105866792B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610378509.9A CN105866792B (zh) 2016-05-31 2016-05-31 一种星载激光雷达树高提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610378509.9A CN105866792B (zh) 2016-05-31 2016-05-31 一种星载激光雷达树高提取方法

Publications (2)

Publication Number Publication Date
CN105866792A CN105866792A (zh) 2016-08-17
CN105866792B true CN105866792B (zh) 2019-04-26

Family

ID=56641830

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610378509.9A Active CN105866792B (zh) 2016-05-31 2016-05-31 一种星载激光雷达树高提取方法

Country Status (1)

Country Link
CN (1) CN105866792B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107479065B (zh) * 2017-07-14 2020-09-11 中南林业科技大学 一种基于激光雷达的林窗立体结构量测方法
CN110070488B (zh) * 2019-04-25 2023-01-03 北京工业大学 一种基于卷积神经网络的多角度遥感影像森林高度提取方法
CN111781605B (zh) * 2020-07-24 2022-07-22 江苏大学 一种基于超声回波信号的树体冠层厚度探测方法
CN112560661A (zh) * 2020-12-10 2021-03-26 首都师范大学 叶面积指数的计算方法、装置、电子设备及可读存储介质
CN117520733B (zh) * 2024-01-05 2024-03-19 云南师范大学 森林冠层高度与地理环境协变量关系确定方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103760565A (zh) * 2014-02-10 2014-04-30 中国科学院南京地理与湖泊研究所 一种区域尺度森林冠层高度遥感反演方法
CN104062644A (zh) * 2013-11-22 2014-09-24 董立新 一种从激光雷达高斯回波数据中提取树高的方法
CN105005047A (zh) * 2015-07-17 2015-10-28 武汉大学 后向散射优化的森林复杂地形校正及树高反演方法、系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104062644A (zh) * 2013-11-22 2014-09-24 董立新 一种从激光雷达高斯回波数据中提取树高的方法
CN103760565A (zh) * 2014-02-10 2014-04-30 中国科学院南京地理与湖泊研究所 一种区域尺度森林冠层高度遥感反演方法
CN105005047A (zh) * 2015-07-17 2015-10-28 武汉大学 后向散射优化的森林复杂地形校正及树高反演方法、系统

Also Published As

Publication number Publication date
CN105866792A (zh) 2016-08-17

Similar Documents

Publication Publication Date Title
CN105866792B (zh) 一种星载激光雷达树高提取方法
Gong et al. ICEsat GLAS data for urban environment monitoring
Popescu et al. Seeing the trees in the forest
Mallet et al. Full-waveform topographic lidar: State-of-the-art
JP6381100B2 (ja) ３次元路面下診断システムおよび３次元路面下診断方法
CN109031344B (zh) 一种全波形激光雷达和高光谱数据联合反演森林结构参数的方法
Pirotti Analysis of full-waveform LiDAR data for forestry applications: a review of investigations and methods
Wang et al. Separation of ground and low vegetation signatures in LiDAR measurements of salt-marsh environments
Wallace Assessing the stability of canopy maps produced from UAV-LiDAR data
KR101080985B1 (ko) 항공 라이다 성과를 이용한 수고 및 개체 수 확인방법
JP2003344048A (ja) 森林情報処理システム
CN103398957A (zh) 基于高光谱与激光雷达提取叶面积垂直分布的方法
CN108959705B (zh) 一种预测亚热带森林生物量的方法
Cui et al. Retrieving crop fractional cover and LAI based on airborne Lidar data
Coluzzi et al. On the LiDAR contribution for landscape archaeology and palaeoenvironmental studies: the case study of Bosco dell'Incoronata (Southern Italy)
CN104062644A (zh) 一种从激光雷达高斯回波数据中提取树高的方法
Apostol et al. Height extraction and stand volume estimation based on fusion airborne LiDAR data and terrestrial measurements for a Norway spruce [Picea abies (L.) Karst.] test site in Romania
CN106443665A (zh) 基于辅助距离的雷达与esm航迹关联方法
Fu et al. Prediction of subtropical forest parameters using airborne laser scanner
Nie et al. A continuous wavelet transform based method for ground elevation estimation over mountainous vegetated areas using satellite laser altimetry
Hyyppä et al. Forest inventory using laser scanning
CN104457756B (zh) 一种基于双机测距的海面物体定位方法
CN109902686B (zh) 一种森林单木参数提取方法
Grasso et al. The use of SLAM and UAV technology in geological field for monitoring: the case study of the Bossea Cave
CN102313537A (zh) 基于单片InSAR正射影像的调绘方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant