CN109946714A - 一种基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法 - Google Patents
一种基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法 Download PDFInfo
- Publication number
- CN109946714A CN109946714A CN201910267303.2A CN201910267303A CN109946714A CN 109946714 A CN109946714 A CN 109946714A CN 201910267303 A CN201910267303 A CN 201910267303A CN 109946714 A CN109946714 A CN 109946714A
- Authority
- CN
- China
- Prior art keywords
- data
- lidar
- forest
- model
- palsar
- 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
- 239000002028 Biomass Substances 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000013461 design Methods 0.000 claims abstract description 7
- 238000011156 evaluation Methods 0.000 claims description 8
- 230000010287 polarization Effects 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 3
- 238000007598 dipping method Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000009897 systematic effect Effects 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 230000003287 optical effect Effects 0.000 claims description 2
- 230000005855 radiation Effects 0.000 claims description 2
- 238000012549 training Methods 0.000 abstract description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 5
- 229910052799 carbon Inorganic materials 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000009877 rendering Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 239000002023 wood Substances 0.000 description 3
- GWEVSGVZZGPLCZ-UHFFFAOYSA-N Titan oxide Chemical compound O=[Ti]=O GWEVSGVZZGPLCZ-UHFFFAOYSA-N 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 239000004408 titanium dioxide Substances 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法，所述的基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法包括分析与设计、LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立步骤，本发明所述的基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法在激光雷达和多角度遥感数据反演森林参数方面，激光雷达的主要优势是能够直接测量树木的高度和垂直结构信息，利用激光雷达提供大量点的精确数据这一优势，然后通过训练多角度遥感数据进行区域推广，就可以实现区域森林垂直结构参数的反演。
Description
技术领域
本发明属于地理科学技术领域，具体涉及一种基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法。
背景技术
森林作为最重要的陆地生态系统，是一个多资源、多功能的综合体，在固定二氧化碳、提供林木资源、调节气候、涵养水源、保护生物多样性、维护生态平衡等方面蕴含着巨大的价值，因此，了解森林生长和消亡的发展规律意义重大，在全球气候变暖的背景下，森林的碳汇功能日益成为研究的热点问题，对森林碳汇量及碳储量进行准确测算，是掌握森林生态系统在全球碳循环中的作用以及碳汇交易的基础，随着科学技术的发展，地理信息系统、遥感技术、全球定位系统和计算机技术在森林资源调查中得到了广泛的应用，为森林资源调查节省了不少人力和物理。尤其是近几年遥感技术的迅速发展，特别是激光雷达（Light Detection And Ranging，LiDAR）技术的出现于应用，更是为我国森林资源调查带来了新的机遇。
发明内容
本发明的目的在于提供一种基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法。
本发明的目的是这样实现的，所述的基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法包括分析与设计、LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立步骤，具体包括：
A、分析与设计：
以待测地森林资源清查数据为基础，选择具有代表性的森林设置样地，对样地内的数目进行每木检尺，分别记录每棵树的胸径、树高、枝下高、冠幅和样地内的相对坐标；
B、LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立：
1）通过机载LiDAR系统测量样地内的点云数据，然后对LiDAR点云数据进行预处理，从中生成DEM、DSM和CHM，提取建立模型所需的相关树高数据；
2）通过对PALSAR数据进行预处理，精确的将地面样地观测数据和SAR影响配准，通过数据生产方提供的定标公式和定标系数将SLC图像的信号转换为用分贝表示的后向散射系数；
3）以样地森林生物量为因变量，以机载LiDAR数据提取的样地林木平均高、极化SAR的敏感信号、极化分解参数为自变量，利用多元数据进行多元后消逐步回归法和强行进入发进行建模。
本发明所述的基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法在激光雷达和多角度遥感数据反演森林参数方面，激光雷达的主要优势是能够直接测量树木的高度和垂直结构信息，利用激光雷达提供大量点的精确数据这一优势，然后通过训练多角度遥感数据进行区域推广，就可以实现区域森林垂直结构参数的反演。
说明书附图
图1为Lidar样地分布示意图；
图2为原始点云数据示意图；
图3为DSM渲染图及图例；
图4为DEM渲染图及图例；
图5为CHM渲染图及图例；
图6为单棵树木位置与数量的估算示意图；
图7为不同最低树高设置下LiDAR平均树高计算示意图；
图8为试验区域实地核查数据示意图；
图9为CHM渲染图及图例；
图10为坡向图及图例；
图11为坡度图及图例。
具体实施方式
下面结合实施例和附图对本发明作进一步的说明，但不以任何方式对本发明加以限制，基于本发明教导所作的任何变换或替换，均属于本发明的保护范围。
本发明所述的基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法包括分析与设计、LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立步骤，具体包括：
A、分析与设计：
以待测地森林资源清查数据为基础，选择具有代表性的森林设置样地，对样地内的数目进行每木检尺，分别记录每棵树的胸径、树高、枝下高、冠幅和样地内的相对坐标；
B、LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立：
1）通过机载LiDAR系统测量样地内的点云数据，然后对LiDAR点云数据进行预处理，从中生成DEM、DSM和CHM，提取建立模型所需的相关树高数据；
2）通过对PALSAR数据进行预处理，精确的将地面样地观测数据和SAR影响配准，通过数据生产方提供的定标公式和定标系数将SLC图像的信号转换为用分贝表示的后向散射系数；
3）以样地森林生物量为因变量，以机载LiDAR数据提取的样地林木平均高、极化SAR的敏感信号、极化分解参数为自变量，利用多元数据进行多元后消逐步回归法和强行进入发进行建模。
B步骤1）中所述的对LiDAR点云数据进行预处理是使用ENVI5.2定义LiDAR点云数据的投影系统和点云类别并进行地理坐标纠正，导入激光点数据后去除高程异常点，将激光点分类，分别从分类后的数据点中提取相应的数据点内插生成数字高程模型（DEM）和数字表面模型（DSM），通过DEM和DSM的计算生成冠层高度模型（CHM）。
B步骤2）中所述的预处理是对PALSAR数据进行辐射定标、斑点噪声滤波、极化分解、正射校正和地形辐射校正。
所述的基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法还包括模型评价和相关数据的精度评价步骤。
所述的模型评价和相关数据的精度评价是通过相关系数检验、F检验、回归系数的t检验、均方根误差、估测精度相关公式、方差膨胀因子进行多重共线性检验的方法来确保构建建模型的显著性和相关数据精度的合理性。
下面以具体实施案例对本发明做进一步说明：
实施例1
样地设置具体见图1，样地坐标分布具体见下表：
样地坐标分布
实施例2
LiDAR点云数据经过系统处理，将最后所得到的激光点进行内插，即可得到较为精确的DEM、DSM：经过DEM提取、DSM提取、IDW插值、CHM生成、平均树高计算等步骤，得到了实验区域的DEM、DSM、CHM、单棵树木位置估计图等，具体见图1~图7。
根据CHM提取平均树高，与核实采样点进行对比分析。通过设置不同的最低树高，平均树高的估计结果如下表所示：
LiDAR点云数据提取平均树高（单位：米/m）
实施例3
选择实验区，区域范围为2平方千米左右，实地共采集高程点48个，单棵树木69株，同时采集了树木的胸径、树高，以及实地照片，如图5所示。该区域的Lidar点云数据经过系统的处理即DEM提取、DSM提取、IDW插值、坡度、坡向生成、CHM生成、平均树高计算等步骤，得到了实验区域的DEM、DSM、CHM、坡度、坡向图，具体见图8~图11。
Claims (5)
1.一种基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法，其特征在于所述的基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法包括分析与设计、LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立步骤，具体包括：
A、分析与设计：
以待测地森林资源清查数据为基础，选择具有代表性的森林设置样地，对样地内的数目进行每木检尺，分别记录每棵树的胸径、树高、枝下高、冠幅和样地内的相对坐标；
B、LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立：
1）通过机载LiDAR系统测量样地内的点云数据，然后对LiDAR点云数据进行预处理，从中生成DEM、DSM和CHM，提取建立模型所需的相关树高数据；
2）通过对PALSAR数据进行预处理，精确的将地面样地观测数据和SAR影响配准，通过数据生产方提供的定标公式和定标系数将SLC图像的信号转换为用分贝表示的后向散射系数；
3）以样地森林生物量为因变量，以机载LiDAR数据提取的样地林木平均高、极化SAR的敏感信号、极化分解参数为自变量，利用多元数据进行多元后消逐步回归法和强行进入发进行建模。
2.根据权利要求1所述的基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法，其特征在于B步骤1）中所述的对LiDAR点云数据进行预处理是使用ENVI5.2定义LiDAR点云数据的投影系统和点云类别并进行地理坐标纠正，导入激光点数据后去除高程异常点，将激光点分类，分别从分类后的数据点中提取相应的数据点内插生成数字高程模型（DEM）和数字表面模型（DSM），通过DEM和DSM的计算生成冠层高度模型（CHM）。
3.根据权利要求1所述的基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法，其特征在于B步骤2）中所述的预处理是对PALSAR数据进行辐射定标、斑点噪声滤波、极化分解、正射校正和地形辐射校正。
4.根据权利要求1所述的基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法，其特征在于所述的基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法还包括模型评价和相关数据的精度评价步骤。
5.根据权利要求4所述的基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法，其特征在于所述的模型评价和相关数据的精度评价是通过相关系数检验、F检验、回归系数的t检验、均方根误差、估测精度相关公式、方差膨胀因子进行多重共线性检验的方法来确保构建建模型的显著性和相关数据精度的合理性。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910267303.2A CN109946714A (zh) | 2019-04-03 | 2019-04-03 | 一种基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910267303.2A CN109946714A (zh) | 2019-04-03 | 2019-04-03 | 一种基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109946714A true CN109946714A (zh) | 2019-06-28 |
Family
ID=67013663
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910267303.2A Pending CN109946714A (zh) | 2019-04-03 | 2019-04-03 | 一种基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109946714A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113009481A (zh) * | 2021-01-15 | 2021-06-22 | 扬州哈工科创机器人研究院有限公司 | 一种基于干涉sar雷达的森林地物成像反演方法 |
| CN114281871A (zh) * | 2022-03-01 | 2022-04-05 | 国家林业和草原局调查规划设计院 | 航空激光雷达估算山地林区chm树高的精度验证方法 |
| CN114924034A (zh) * | 2022-05-06 | 2022-08-19 | 贵州师范大学 | 一种基于生态过程模型的林业碳计量系统 |
| CN115082274A (zh) * | 2022-06-09 | 2022-09-20 | 贵州师范学院 | 基于卫星遥感的地表植被碳汇估算与交易方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104656098A (zh) * | 2015-02-03 | 2015-05-27 | 南京林业大学 | 一种遥感森林生物量反演的方法 |
| WO2016132106A1 (en) * | 2015-02-18 | 2016-08-25 | The University Court Of The University Of Edinburgh | Satellite image processing |
| CN208155273U (zh) * | 2018-05-29 | 2018-11-27 | 海南省林业科学研究所 | 一种基于激光扫描技术的便携式森林蓄积量测定仪 |
| CN108921885A (zh) * | 2018-08-03 | 2018-11-30 | 南京林业大学 | 一种综合三类数据源联合反演森林地上生物量的方法 |
-
2019
- 2019-04-03 CN CN201910267303.2A patent/CN109946714A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104656098A (zh) * | 2015-02-03 | 2015-05-27 | 南京林业大学 | 一种遥感森林生物量反演的方法 |
| WO2016132106A1 (en) * | 2015-02-18 | 2016-08-25 | The University Court Of The University Of Edinburgh | Satellite image processing |
| CN208155273U (zh) * | 2018-05-29 | 2018-11-27 | 海南省林业科学研究所 | 一种基于激光扫描技术的便携式森林蓄积量测定仪 |
| CN108921885A (zh) * | 2018-08-03 | 2018-11-30 | 南京林业大学 | 一种综合三类数据源联合反演森林地上生物量的方法 |
Non-Patent Citations (2)
| Title |
|---|
| 王新云 等: "基于HJ1B和ALOS/PALSAR数据的森林地上生物量遥感估算", 《生态学报》 * |
| 范凤云: "基于机载LIDAR和极化SAR数据的山区森林蓄积量估测方法研究", 《中国优秀博硕士学位论文全文数据库(硕士)农业科技辑》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113009481A (zh) * | 2021-01-15 | 2021-06-22 | 扬州哈工科创机器人研究院有限公司 | 一种基于干涉sar雷达的森林地物成像反演方法 |
| CN114281871A (zh) * | 2022-03-01 | 2022-04-05 | 国家林业和草原局调查规划设计院 | 航空激光雷达估算山地林区chm树高的精度验证方法 |
| CN114924034A (zh) * | 2022-05-06 | 2022-08-19 | 贵州师范大学 | 一种基于生态过程模型的林业碳计量系统 |
| CN114924034B (zh) * | 2022-05-06 | 2024-05-14 | 贵州师范大学 | 一种基于生态过程模型的林业碳计量系统 |
| CN115082274A (zh) * | 2022-06-09 | 2022-09-20 | 贵州师范学院 | 基于卫星遥感的地表植被碳汇估算与交易方法及系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108921885B (zh) | 一种综合三类数据源联合反演森林地上生物量的方法 | |
| Alonzo et al. | Mapping urban forest leaf area index with airborne lidar using penetration metrics and allometry | |
| Nie et al. | Above-ground biomass estimation using airborne discrete-return and full-waveform LiDAR data in a coniferous forest | |
| Tang et al. | Retrieval of vertical LAI profiles over tropical rain forests using waveform lidar at La Selva, Costa Rica | |
| Erdody et al. | Fusion of LiDAR and imagery for estimating forest canopy fuels | |
| Xing et al. | An improved method for estimating forest canopy height using ICESat-GLAS full waveform data over sloping terrain: A case study in Changbai mountains, China | |
| CN109946714A (zh) | 一种基于LiDAR和ALOS PALSAR多元数据的森林生物量模型的建立方法 | |
| Zhang et al. | Mapping forest stand age in China using remotely sensed forest height and observation data | |
| Fang | Canopy clumping index (CI): A review of methods, characteristics, and applications | |
| Chen et al. | A data-model fusion approach for upscaling gross ecosystem productivity to the landscape scale based on remote sensing and flux footprint modelling | |
| Pfeifer et al. | Leaf area index for biomes of the Eastern Arc Mountains: Landsat and SPOT observations along precipitation and altitude gradients | |
| Chu et al. | Fractional vegetation cover | |
| Rautiainen et al. | Coupling forest canopy and understory reflectance in the Arctic latitudes of Finland | |
| He et al. | Using LiDAR data to measure the 3D green biomass of Beijing urban forest in China | |
| Coomes et al. | Airborne laser scanning of natural forests in New Zealand reveals the influences of wind on forest carbon | |
| Wang et al. | Fractional vegetation cover estimation method through dynamic Bayesian network combining radiative transfer model and crop growth model | |
| Tesfamichael et al. | Investigating the impact of discrete-return lidar point density on estimations of mean and dominant plot-level tree height in Eucalyptus grandis plantations | |
| Yuan et al. | Examination of the quantitative relationship between vegetation canopy height and LAI | |
| CN109934109A (zh) | 一种基于遥感的黄土高原水土流失区林草植被信息提取方法 | |
| CN107247809A (zh) | 一种人工林森林年龄空间制图的新方法 | |
| Zhang et al. | Understory biomass measurement in a dense plantation forest based on drone-SfM data by a manual low-flying drone under the canopy | |
| Gonzalez Musso et al. | Applying unmanned aerial vehicles (UAVs) to map shrubland structural attributes in northern Patagonia, Argentina | |
| Hojo et al. | Modeling forest above-ground biomass using freely available satellite and multisource datasets | |
| Bottai et al. | Application of a single-tree identification algorithm to LiDAR data for the simulation of stem volume current annual increment | |
| CN117035174A (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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190628 |