森林地上生物量遥感估测研究进展

合集下载

合成孔径雷达数据估测森林生物量研究综述

合成孔径雷达数据估测森林生物量研究综述
e td. ne K e r :S y wo ds AR ;o e tb o s s etmain f rs ima ; si t o
引 言
近年来 , 森林生物 量的定 量测算一 直是生 态学及 相关领域的一个重 要研究 问题 , 特别是 在全球 气候 变 化、 温室效应加剧 、 界各 国面临温 室气体 减排 、 世 有效 控制温室气体增速 、 履行 《 京都议定书 》 国际协定 的背 景下 , 开展森林生物量的数量 、 时空分 布和动态变化研 究, 具有重要 的理论意义及实践价值 。在 电磁波谱 中, 波 长在 Im m—Im的波段 范围称 为微波 , 微波遥 感具 有穿透树冠 的能力 , 不仅能和树叶发生作用 , 而且 主要 是和森林 生物量 的主体——枝和树干发生作用。特别 是在多云雨和雾 的热带 和亚热带地 区 , 可见光 和红外 遥感受到 了很大限制 , 因而微 波遥感 为森林 生物量全 面 和精确估测 提供 了可行 的工 具 】 。微波 遥感 , 可
测 方面发挥着越 来越重要的作 用。首 先对森林 生物量 S R估 算机 理和方法进 行介绍 , A 然后 分析 主要 的 影 响因素 , 最后 指 出今后 的发展趋 势。
关键词 : 合成孔径 雷达 ; 生物量 ; 估测
中图分类号 : P7 ; 5 T 9 S7 8
文献标识码 : A
a yt n i me, a e o r n r mp ra tt si t oe tb o s . n t i a e ,h c a s a d h s b c me mo e a d moe i ot n o etmae fr s ima s I hs p p r te me h nim n
分为主动 和被 动 2种方 式 , 成孔 径 雷 达 ( ytec 合 Snht i A e ueR dr简称 S R) pr r a a, t A 作为 主动微 波遥感 的方式

运用林分密度和平均高估测思茅松人工林地上生物量

运用林分密度和平均高估测思茅松人工林地上生物量

第49卷第1期东 北 林 业 大 学 学 报Vol.49No.12021年1月JOURNALOFNORTHEASTFORESTRYUNIVERSITYJan.20211)国家自然科学基金项目(31860182、31260156、41571372、42061072),云南省教育厅科学研究基金项目[2018JS330]。

第一作者简介:张国飞,男,1978年12月生,西南林业大学林学院,博士研究生。

E-mail:396316745@qq.com。

通信作者:岳彩荣,西南林业大学林学院,教授。

E-mail:cryue@163.com。

收稿日期:2020年3月30日。

责任编辑:王广建。

运用林分密度和平均高估测思茅松人工林地上生物量1)张国飞 岳彩荣 赵勋 罗洪斌 谷雷(西南林业大学,昆明,650224) 摘 要 以65块云南省普洱地区思茅松人工林圆形样地数据和sentinel-2多光谱影像数据为研究对象,利用林分平均高与林分密度(每公顷株数、林分疏密度、植被覆盖度、叶面积指数)估测思茅松人工林林分地上生物量。

分析思茅松人工林林分地上生物量与林分密度指标的相关性;采用参数模型(不变参数模型和可变参数模型)和非参数模型(包括支持向量机、随机森林和BP神经网络)探索平均高和林分密度等变量估测林分思茅松人工林地上生物量。

结果表明:思茅松人工林林分地上生物量与每公顷株树、林分疏密度、植被覆盖度、叶面积指数呈显著正相关(r>0.5);在构建思茅松人工林地上生物量的所有模型中,每公顷株数-林分平均高构建的可变参数模型(R2=0.9660,RMSE=10.05t·hm-2)效果最优,林分平均高-林分疏密度构建的RF模型(R2=0.9017,RMSE=19.37t·hm-2)次之,林分平均高-植被覆盖度构建的RF模型(R2=0.7484,RMSE=33.36t·hm-2)最差;林分密度-平均高的地上生物量模型与实测地上生物量的相关性较高(R2=0.9660),反演误差值较低(RMSE=10.05t·hm-2);叶面积指数比植被覆盖度对林分地上生物量变动有更好的解释能力,每公顷株数对林分地上生物量变动的解释能力好于林分疏密度。

生态环评中森林植被生物量的估算方法

生态环评中森林植被生物量的估算方法

生态环评中森林植被生物量的估算方法
在生态环评中,估算森林植被生物量的方法可以通过以下几种途径进行:
1. 样地调查法:选择代表性的样地,对其中生长的植被进行调查和测量。

通过统计每
个样地中不同植物组成的生物量,再根据各样地的面积计算总体生物量。

2. 遥感和影像解译法:利用遥感数据和图像解译技术,对森林覆盖面积和植被类型进
行判读和分类,再根据不同的植被类型建立相应的生物量模型,通过计算该地区各植
被类型的面积和对应的生物量模型,估算总体生物量。

3. 生态模型法:根据森林植被的生长特征、环境参数以及人为干扰因素,建立相应的
生态模型。

通过输入相关参数,模型可以预测森林植被的生物量。

4. 无人机遥感法:利用无人机搭载的遥感传感器,对森林植被进行高分辨率的遥感监测。

通过获取的遥感数据,结合地面测量数据进行分析,得出森林植被的生物量估算。

以上是一些常用的森林植被生物量估算方法,在实际应用中可以根据具体的情况选择
适合的方法和技术。

长白山林区森林生物量遥感估测模型

长白山林区森林生物量遥感估测模型
Fa e y L n z Ya gJn n nW n i iMi g e n i mig ( ol efF r t , r es F r t nvrt H ri 5 0 0 C lg oe r Not at oe r U i sy ab 10 4 ) e o sy h sy e i n
9 9 ・m~, .2th 样地 生物 量 真 实值 与 预 测值 相 关 系数 0 803 偏 最 小 二 乘 回归 法 要 优 于 逐 步 回归 法 。 利 用 建 立 的 偏 .6 ,
最 小 二 乘 回归 模 型 计 算 得 到 黑 龙 江 长 白山 2 0 0 7生 物 量 等 级 分 布 图 , 用 2 采 9个 检 验 样 本 对 反 演 结 果 进 行 检 验 , 算 计 得到 2 9个 样 本 的平 均 预 测 精 度 为 8 .3 。 37% 关键词 : T M;森林 生 物 量 ; 步 回归 ; 最 小 二乘 回归 ; otrp 逐 偏 bo t ;长 白山 sa
df rn a d f ryv le(n ldn 1vgtt nid x ,txueifr t n a d e v ome tlf tr.T es p i iee t n so a a f b g u icu igl eeai n e ) etr nomai n n i n na a os h t w s o o r c e e
( 北林业大学林学院 东 哈 尔滨 10 4 ) 5 0 0
摘 要 : 采 用 黑 龙 江长 白 山地 区 T 图像 和 13块 森 林 资 源 连 续 清 查 固定 样 地 数 据 及 野 外 调 查 补 充 样 地 数 据 , M 4 选 择 包 括 各 波 段 灰 度 值 、 同波 段 灰 度 值 之 间 的线 性 和 非 线 性 组 合 ( 括 1 种 植 被 指数 ) 纹 理 信 息 以及 环 境 因子 在 不 包 1 、 内的 7 5个 自变 量 , 别 采 用 逐 步 回归 分 析 法 和偏 最 小 二 乘 回 归 法 建 立 黑 龙 江 长 白 山林 区 森 林 生 物 量 遥 感 估 测 模 分 型 : 步 回归 法 采 用 5个 自变 量 所 建 模 型 平 均 拟 合 精 度 为 7 . % , 方 根 误 差 为 1. 2th 逐 65 均 9 1 ・m~ , 地 生 物 量 真 实 值 样 与 预 测 值 相 关 系 数 为0 5 34 偏 最 小 二 乘 回归 法 采 用 l .4 ; 0个 自变 量 所 建 模 型 平 均 拟 合 精 度 8 .% , 方 根 误 差 58 均

森林生物量遥感估算与应用分析

森林生物量遥感估算与应用分析
目前 L n st M E M 和 N A / V R 数 据 a da T / T 0 AA H R
展的 3 S技术 为大 尺 度森 林 生物 量估 算 提供 了一 条 快捷 、 济 、 经 方便 和 可靠 的途径[ 7 1 。
本 文对 基 于 遥 感 信 息 的 森林 生 物 量 估算 方 法 与应 用 及其发 展 方 向做 了深 入分 析 . 以期 为我 国森 林 生物量研 究 的发展 提供 参考
1 引 言
森 林 生物 量 约 占全 球 陆地 植 被 生 物 量 9 %f 标 志 . 是 评 估森 林 碳 也
植 被 的遥 感 图像 信 息 是 由其 反 射 光谱 特 征 决
定的 . 植物 的光合 作 用表 现 为 对红 光 和蓝 紫光 的强
优 遥 感光 谱模 型 . 而 确定 了粤西 及 附 近地 区 的森 进 林 生 物 量 ㈣; a e等 结 合 地 面 调 查 和 T A H R Hm M、 V R 数 据 . 欧 洲 森林 生物 量进 行 了成 功估 算 [: e k 对 U Lf y l s 等利 用 雷 达 数 据 对美 国 E s May n at rl d落 叶松 的地 a
维普资讯
第 8 第 4期 卷 20 0 6年 1 2月
地 球 信 息 科 学
GE I ORMAT ON C ENC O— NF I S I E
V 1 ,No4 o. 8 . De .2 0 c, 0 6
森林生物 量遥感估 算与应用分析
徐新良 ,
林生 物量 研究 始 于 2 0世 纪 7 0年 代后 期 . 后 建 立 先 了主要森 林树 种 的生 物量 测定 相 对 生长 方 程 . 算 估
了它们 的生 物 量 和生 产 力 [ 初 步 总 结 了全 国不 同 3 ] . 森林 类 型 的生物 量 与 生产 力及 其 空 间分 布格 局

转载:森林生物量、碳储量的遥感反演与计算

转载:森林生物量、碳储量的遥感反演与计算

转载:森林⽣物量、碳储量的遥感反演与计算 摘要:利⽤地⾯调查获取的森林⽣物量为因变量,以GIS因⼦和遥感因⼦为⾃变量,建⽴多元线性回归模型,再进⾏⽣物量与碳储量的换算。

实现由点到⾯的森林⽣物量及碳储量的获取。

1数据获取因变量的获取:地⾯样地⽣物量数据的获取⽅法是设置森林调查样地,调查样地乔⽊的每⽊胸⾼直径,计算平均胸径,利⽤各树种已建⽴的森林⽣物量相对⽣长式,计算树⼲、树枝、树叶、树根⽣物量,最后得出地⾯样地单位⾯积⽣物量. 利⽤常⽤的碳转换系数乘以森林⽣物量进⾏森林碳储量计算.在南北向⼭体的两侧,从低到⾼,海拔相差约 50 m左右设置调查样地,样地⾯积0.04~0.06 hm2 (20 m ×20 m, 20 m ×30 m) ,调查样地中乔⽊的胸⾼直径,计算平均胸径,利⽤西部常绿阔叶林、温性针叶林、华⼭松林等⽣物量相对⽣长式计算树⼲、树枝、树叶、树根⽣物量,最后得出地⾯样地单位⾯积⽣物量. 样地数量为31个.⾃变量的获取:利⽤GPS获取地⾯样地坐标,并提取与地⾯样地坐标对应的遥感光谱值、波段⽐值、地学因⼦值作为估测模型的⾃变量.遥感数据采⽤2006年1⽉的印度卫星IRS -P6 LISS 3影像. 数据包括B2, B3, B4, B5 4个波段,⼏何分辩率为25 m ×25 m, 1个全⾊波段,⼏何分辩率为5.8m.利⽤1∶5万地形图,选取30个明显地物点,对影像进⾏⼏何校正,校正误差控制在⼀个像元内. 校正后的像元⼤⼩为25 m ×25 m,和地⾯样地⾯积基本吻合,以样地GPS实测坐标为准,在ArcMap下提取4个波段值(B2, B3, B4, B5波段)及3个波段⽐值( (B4 - B3) / (B4 +B3) 、B4 /B3、B4 /B2) ,共7个输⼊变量.利⽤GIS技术提取与⽣物量相关的部分地学因⼦参与⽣物量估测是提⾼估测精度的有效⽅法. 利⽤1∶5万地形图通过建⽴像元为25 m ×25 m的格⽹DEM提取海拔、坡度、坡向值作为3个地学因⼦变量.上述波段值、波段⽐值及地学因⼦共10个因⼦作为模型的候选变量.2建⽴模型以上述波段值、波段⽐值及地学因⼦共10个因⼦作为模型的⾃变量,样地⽣物量为因变量,利⽤Matlab 7.0 的stepwise函数进⾏逐步回归分析,在0.05显著性⽔平下, 通过F 检验筛选变量, 则有B2、B3、B4、B5 4个波段值及海拔、坡度、坡向值7个变量⼊选, 3个波段⽐值被剔除. 对⼊选的7个变量利⽤回归分析regress函数计算,建⽴森林⽣物量最优估测⽅程.Y = - 0.149 4 ×海拔+ 2.620 1 ×坡度- 1.264 5×坡向- 48.227 ×B2 + 34.856 ×B3 +1.569 7 ×B4 – 12.08 ×B5 +2 754.7式中: Y为森林⽣物量( t/ hm2 );B2,B3,B4,B5分别为4个波段值.3⽅差分析及线性回归关系显著性检验根据⽅差分析,模型的F = 3.76, F > F0.05 = 2.44,F > F0.01 =3.54. 模型的7个遥感及地学因⼦和森林⽣物量在0.01⽔平下具有极显著相关关系,相关系数R = 0.730 5. 说明模型可⽤于森林⽣物量估测.⽣物量模型似合效果见图1.(问:应该预留部分数据⽤于检验,如果拿建模的数据来检验说服⼒不强)4 常绿阔叶林⽣物量的估测对IRS - P6数据经⼏何校正、图像增强、假彩⾊合成处理后得到的合成图像,在ArcMap 下采⽤⽬视解释的⽅法区划常绿阔叶林⼩班,并单独提取常绿阔叶林图层. 以该图层为边界,提取和该图层对应的上述7个因⼦.利⽤已建⽴的回归模型计算像元⽣物量,并⽣成⽣物量图层,最后以常绿阔叶林⼩班为边界,取⼩班边界内平均值作为⼩班单位⾯积⽣物量. M. 则有:⼩班⽣物量= . M ×⼩班⾯积.(问:计算结果直接汇总就能算出总⽣物量,为何要⽤均值计算?)5 常绿阔叶林碳储量的估测通常,对植物⽣物量转化为碳储量是按照植物⼲重有机物中碳所占的⽐重进⾏计算. 树种组成、年龄和种群结构不同,转化率也不同,但差异不⼤,⼀般在0.45~0.5之间变化. 由于获取各种植被类型的转化率⽐较困难,所以国际上常⽤的转化率为0.50[ 3, 4 ] . 在此采⽤0.5的碳转化系数计算碳储量.则有:⼩班碳储量=⼩班⽣物量×0.5下⽂将介绍在外业调查数据的基础上,如何计算⽣物量。

6种地形校正方法对雷竹林地上生物量遥感估算的影响

董德 进 周 国模 杜 华 强 徐 小 军 崔 瑞 蕊 沈振 明
临 安 3 10 : 13 0 (.浙 江 省 森 林 生 态 系 统 碳 循 环 与 固碳 减 排 重 点 实 验 室 1 临 安 3 10 ; .浙 江 农 林 大 学 环 境 与资 源 总 站 临 安 3 10 ) . 13 0
o i t p g a h c o r c i n fs x o o r p i c re t mo e s i c u i g Te l t— r g e so o d l , n l d n il e e r s i n, Co i e, C, S sn CS, S CS + C, a d n M i n e t o n ar, n a o e r u d f r s ima s e tma i n o y o t c y r e o b v g o n o e t b o s si to f Ph U sa h s p a c x. Re u t h we s l s o d: E c p o s n n S o r c i n s x e t f r Co i e a d CS c r e t o me h d ,t e r s o r me h d r v d d a if co y c r e to e u t ; S x t p g a h c c r e t n me h d l i r v d t o s h e tf u t o s p o i e s tsa t r o r c i n r s ls i o o r p i o r c i t o s a l mp o e o
中图 分 类 号 :¥ 7 . 7 18 文献标识码 : A 文 章 编 号 : 0 1 7 8 ( 0 1 1 00 — 8 1 0 — 4 8 2 l ) 2— 0 1 0

基于无人机多光谱遥感的春玉米叶面积指数和地上部生物量估算模型比较研究

中国农业科技导报ꎬ2021ꎬ23(9):112-120JournalofAgriculturalScienceandTechnology㊀收稿日期:2020 ̄03 ̄04ꎻ接受日期:2020 ̄05 ̄27㊀基金项目:国家重点研发计划项目(2016YFD0300605 ̄02)ꎻ国家自然科学基金项目(31871519)ꎻ北京市农林科学院2020年度科研创新平台建设项目(PT2020 ̄24)ꎻ国家玉米产业技术体系专项(CARS ̄02 ̄87)ꎮ㊀联系方式:樊鸿叶E ̄mail:1173666393@qq.comꎻ∗通信作者刘玉华E ̄mail:hblyh@126.com基于无人机多光谱遥感的春玉米叶面积指数和地上部生物量估算模型比较研究樊鸿叶1ꎬ2ꎬ3ꎬ㊀李姚姚4ꎬ㊀卢宪菊2ꎬ3ꎬ㊀顾生浩2ꎬ3ꎬ㊀郭新宇2ꎬ3ꎬ㊀刘玉华1∗(1.河北农业大学农学院ꎬ河北保定071000ꎻ2.北京农业信息技术研究中心ꎬ数字植物北京重点实验室ꎬ北京100097ꎻ3.国家农业信息化工程技术研究中心ꎬ北京100097ꎻ4.中国农业科学院作物科学研究所ꎬ北京100081)摘㊀要:明确基于无人机多光谱遥感的玉米叶面积指数(LAI)和地上部生物量的最优估算模型对获取即时㊁无损㊁可靠的长势关键参量具有重要意义ꎮ2018 2019年ꎬ以郑单958(ZD958)和先玉335(XY335)为研究对象ꎬ设置4个施氮处理ꎬ通过无人机搭载多光谱相机获取多光谱影像ꎬ分析两品种LAI和地上部生物量与植被指数相关性ꎬ分别构建了基于植被指数的LAI和地上部生物量预测模型ꎮ结果表明:同一植被指数在两品种中对施氮量的变化响应规律不同ꎻ在吐丝期ꎬ幂函数对ZD958的LAI和地上部生物量估算效果最好ꎬ指数函数对XY335的LAI估算效果好ꎬ幂函数对地上部生物量估算效果好ꎻ在灌浆期ꎬ幂函数对两品种的LAI估算效果最佳ꎬ而指数函数对两品种的地上部生物量估算效果最好ꎮ研究结果为进一步提高春玉米长势监测的精度提供了重要依据ꎮ关键词:春玉米ꎻ叶面积指数ꎻ地上部生物量ꎻ多光谱遥感ꎻ植被指数ꎻ回归模型doi:10.13304/j.nykjdb.2020.0173中图分类号:S127㊀㊀㊀文献标识码:A㊀㊀㊀文章编号:1008 ̄0864(2021)09 ̄0112 ̄09ComparativeAnalysisofLAIandAbove ̄groundBiomassEstimationModelsBasedonUAVMultispectralRemoteSensingFANHongye1ꎬ2ꎬ3ꎬLIYaoyao4ꎬLUXiaju2ꎬ3ꎬGUShenghao2ꎬ3ꎬGUOXinyu2ꎬ3ꎬLIUYuhua1∗(1.CollegeofAgronomyꎬHebeiAgriculturalUniversityꎬHebeiBaoding071000ꎬChinaꎻ2.BeijingKeyLabofDigitalPlantꎬBeijingResearchCenterforInformationTechnologyinAgricultureꎬBeijing100097ꎬChinaꎻ3.NationalEngineeringResearchCenterforInformationTechnologyinAgricultureꎬBeijing100097ꎬChinaꎻ4.InstituteofCropSciencesꎬChineseAcademyofAgriculturalSciencesꎬBeijing100081ꎬChina)Abstract:Determiningtheoptimalestimationmodelofmaizeleafareaindex(LAI)andabove ̄groundbiomassbasedonUAVmultispectralremotesensingplaysasignificantroleinobtainingtimelyꎬnon ̄destructiveꎬandreliablegrowthparameters.ThisstudysetfournitrogentreatmentswithZhengdan958(ZD958)andXianyu335(XY335)asmaterialsfrom2018to2019.MultispectralimageswereacquiredbyadroneequippedwithamultispectralcameraꎬandtherelationshipbetweenLAIandabovegroundbiomassandvegetationindexwasanalyzedꎬandavegetationindex ̄basedLAIandabovegroundbiomasspredictionmodelwasconstructed.Theresultsshowedthattheresponseofthesamevegetationindextothechangeofnitrogenapplicationwasdifferentbetweentwovarieties.AtsilkingstageꎬthepowerequationreceivedthebestrankinginestimatingtheLAIandabove ̄groundbiomassforZD958ꎬandtheexponentialandpowerequationreceivedthebestrankinginestimatingLAIandabove ̄groundbiomassforXY335ꎬrespectively.AtfillingstageꎬthepowerequationreceivedthebestrankingoffittinggoodnessinestimatingLAIꎬwhileexponentialdidinestimatingabove ̄groundbiomassfortwomaizevarieties.Thisstudyprovidedanimportantbasisandtechnicalmethodforfurtherimprovingtheaccuracyofmonitoringgrowthparametersforspringmaize.Keywords:springmaizeꎻleafareaindexꎻabove ̄groundbiomassꎻmultispectralremotesensingꎻvegetationindexꎻregressionmodel㊀㊀无人机遥感具有时空分辨率高㊁成本低㊁灵活性和可重复性强等优势[1]ꎬ能够通过搭载的传感器获取作物冠层反射的电磁波信息ꎬ进而提取与作物长势相关的参量ꎬ为获取田块尺度上即时㊁无损㊁可靠的作物长势信息提供了一种重要手段[2]ꎮ地上部生物量与叶面积指数(leafareaindexꎬLAI)是反映作物生长状况的重要农艺参数[3 ̄4]ꎬ且二者间存在密切的联系[5]ꎮ在作物不同的生长发育阶段ꎬ由于作物冠层结构㊁叶片形态以及生理生态特征的变化ꎬ作物冠层在图像中表现出的光谱特征亦有不同[6]ꎮ植被指数(vegetationindexꎬVI)是两个或多个波段的反射率经过线性或非线性组合运算来增强植被信息ꎬ以削弱环境背景对植被光谱特征的干扰[7]ꎬ如归一化植被指数(normalizeddifferencevegetationindexꎬNDVI)[8]㊁增强植被指数(enhancedvegetationindexꎬEVI)[9 ̄10]和差值植被指数(differencevegetationindexꎬDVI)[11]等ꎮ目前ꎬ光谱植被指数法已经作为一种快速获取地表信息的手段ꎬ被国内外学者用于作物长势监测研究ꎮ杨贵军等[12]通过无人机搭载多传感器ꎬ实现了对作物产量㊁LAI㊁冠层温度等多种信息的监测ꎮ王亚杰[13]基于无人机获取的多光谱数据ꎬ系统比较了玉米叶面积指数的不同监测方法ꎬ结果表明ꎬ在不同水分处理下基于EVI构建的一元线性模型能够较好的预测玉米灌浆期和成熟期的LAIꎮ高林等[14]基于RVI㊁NDVI和DVI等采用经验模型法构建了大豆LAI反演模型ꎬ结果表明ꎬ在鼓粒期基于NDVI构建的大豆LAI效果最佳(R2=0.829ꎬRMSE=0.301)ꎮ孙诗睿等[15]基于无人机多光谱影像ꎬ通过对比不同植被指数建模方法分析冬小麦叶面积指数的反演精度ꎬ发现基于多植被指数的随机森林预测具有更好的拟合效果(R2=0.882ꎬRMSE=1.218)ꎮ邹楠等[16]研究结果表明ꎬ种植密度在6㊁12和18万株 hm-2条件下ꎬ春玉米RVI㊁DVI和DVI分别与LAI的相关性最好ꎬR2分别为0.762㊁0.691和0.648ꎬ可用于监测春玉米LAIꎮ肖武等[17]利用玉米灌浆期地上部生物量与各种植被指数建立一元曲线回归模型ꎬ结果表明ꎬ植被指数NDVI和GNDVI效果较好(R2>0 68)ꎮ庄东英等[18]研究表明ꎬRVI比NDVI对冬小麦生物量的估算具有明显优势ꎮ已有研究大都基于光谱参数对作物LAI和地上部生物量进行模拟估算ꎬ但鲜有研究对作物不同品种LAI和地上部生物量的光谱估算模型进行探讨ꎬ本研究以此为切入点ꎬ以无人机搭载多光谱相机所获多光谱影像数据为基础ꎬ选取7种常用的植被指数与玉米LAI和地上部生生物量进行相关性分析ꎻ综合4种基于光谱植被指数的模拟方法(指数㊁线性㊁对数㊁幂)ꎬ构建2个玉米品种在吐丝期和灌浆期的LAI和地上部生物量估算最优模型ꎬ为进一步提高春玉米长势监测的精度提供依据和方法ꎮ1㊀材料与方法1.1㊀试验材料供试材料为紧凑型郑单958(Zhengdan958ꎬZD958)和半紧凑型先玉335(Xianyu335ꎬXY335)ꎬ由河南豫研种子科技公司和山东登海先锋种业有限公司提供ꎬ施用化肥为尿素(总氮ȡ46.4%)㊁过磷酸钙(有效磷ȡ12%ꎬ水溶性磷ȡ7 0%ꎬ硫ȡ8.0%)和氯化钾(K2Oȡ60%)ꎬ分别由山东润银生物化工股份有限公司㊁江苏省威力磷复肥有限公司和中化化肥有限公司生产ꎮ所用无人机为大疆M600Proꎬ最大载重5kgꎬ最大续航时间30minꎮ搭载于无人机上的多光谱传感器为Rededge ̄MX(MicaSenseꎬUSA)ꎬ该相机有5个波段ꎬ相机焦距为5.5mmꎬ视场角为47.2ʎꎬ图像分辨率为1280pixelsˑ960pixelsꎬ该相机配备了光强传感器及标定板ꎬ其中光强传感器可校正航拍过程中外界光线的变化对光谱影像造成的影响ꎬ标定板具有固定的反射率ꎬ可以利用标定板进行辐射定标ꎮ1.2㊀试验处理设计试验于2018 2019年4 10月在位于吉林省公主岭市的中国农科院公主岭试验基地(43ʎ11ᶄ 44ʎ9ᶄN㊁124ʎ02ᶄ 125ʎ18ᶄE)进行ꎮ该地区属于温带半湿润大陆性季风气候ꎬ长年平均降雨量为595mmꎬ土壤类型主要为黑土ꎬ土壤耕3119期樊鸿叶等:基于无人机多光谱遥感的春玉米叶面积指数和地上部生物量估算模型比较研究层pH6.0ꎬ有机质含量29.68g kg-1㊁全氮1.28g kg-1㊁速效磷28.16mg kg-1㊁速效钾187.74mg kg-1ꎮ设置4个氮肥处理ꎬ分别为无肥N0(0kg hm-2)㊁低肥N1(150kg hm-2全部一次性基肥)㊁中肥N2(225kg hm-2ꎬ基施150kg hm-2ꎬ拔节追施75kg hm-2)和髙肥N3(300kg hm-2ꎬ基施150kg hm-2ꎬ拔节期㊁吐丝期各施75kg hm-2)ꎬ氮㊁磷㊁钾分别为尿素㊁过磷酸钙(P2O542.5kg hm-2)㊁氯化钾(K2O42.5kg hm-2)ꎬ种植密度为67500株 hm-2ꎬ行距0.65mꎬ株距0 22mꎮ除施肥处理不同外ꎬ各处理田间管理方式均保持一致ꎮ1.3㊀玉米地上部生物量及LAI值获取于2018年的7月21日(吐丝期)㊁8月25日(灌浆期)和2019年的7月27日(吐丝期)㊁8月27日(灌浆期)获取田间数据ꎬ将每个处理的试验小区所选取3株具有代表性的玉米植株带回实验室后ꎬ先放置在烘箱中105ħ杀青30minꎬ然后80ħ烘干48h以上ꎬ直至恒重ꎬ最后称取地上部生物量ꎮ测量植株所有叶片的长和最大叶宽ꎬ根据式(1)计算叶面积指数ꎮLAI=0.75ρðni=1liwiN(1)式中ꎬ0.75表示玉米叶面积校正系数ꎻρ表示种植密度ꎬ株 m-2ꎻN表示取样株数ꎬ株ꎻn为总叶片数ꎻli和wi分别表示第i片叶的叶长和最大叶宽ꎬmꎮ1.4㊀无人机多光谱影像数据获取试验采用大疆M600Pro六旋翼小型电动无人机ꎬ搭载多光谱传感器MicasenseRededge ̄MXꎬ上述田间取样当天选择在晴朗㊁无风或微风㊁无云的天气状况下在10:00 14:00采集多光谱影像数据ꎬ无人机飞行高度为40mꎬ航向为70%ꎬ旁向为65%ꎮ1.5㊀多光谱影像数据处理将所获取的四个生育期多光谱影像使用软件Pix4Dmapper进行拼接ꎬ获得覆盖整个试验区的多光谱影像ꎬ拼接前先人为筛选剔除航带外以及无人机调头时的影像ꎬ并输入每个波段对应标定板的反射率ꎬ最终可获得五张试验区单波段的正射图像ꎮ将拼接好的影像导入ENVI5.1中进行几何校正ꎬ以高清数码影像为参考影像ꎬ均匀选取20个点校正所获无人机多光谱影像ꎬ误差控制在0.5个像元以内ꎬ主要是为消除因无人机飞行高度及姿态不稳定对影像造成的影响ꎮ然后将校正后的影像进行波段组合ꎬ将单波段的影像叠加到一起ꎬ最后进行波段运算ꎬ根据需求计算相关植被指数ꎮ将计算所得的植被指数影像在ENVI5.1软件中进行参数提取ꎬ每个小区绘制3个均匀大小的感兴趣区ꎬ绘制时避开小区边界四个方向各一行玉米ꎬ以排除边行长势不均匀植株ꎬ进而提取出每个感兴趣区的平均光谱植被指数ꎬ每个感兴趣区域作为一次重复ꎬ共三次ꎮ1.6㊀光谱植被指数计算本研究选取了与作物长势相关的7个植被指数:NDVI[6]㊁GNDVI[17]㊁RVI[19]㊁GOSAVI[20]㊁VIopt[21]㊁EVI[8]㊁DVI[11]ꎬ计算分别公式如下ꎮNDVI=(RNIR-RRED)/(RNIR+RRED)(2)GNDVI=(RNIR-RGREEN)/(RNIR+RGREEN)(3)RVI=RNIR/RRED(4)GOSAVI=(1+0.16)(RNIR-RGREEN)/(RNIR+RGREEN+0.16)(5)VIopt=1.45(RNIRˑRNIR+1)/(RRED+0.45)(6)EVI=2.5(RNIR-RRED)/(RNIR+6RRED-7.5RRLUE+1)(7)DVI=RNIR-RRED(8)式中ꎬRNIR㊁RRED㊁RGREEN和RBLUE分别表示近红外波段㊁红光波段㊁绿光波段和蓝光波段的反射率ꎮ1.7㊀数据分析采用SPSS22.0软件进行试验数据统计分析ꎬ用Duncan法进行处理间的多重比较(P<0.05)ꎬ用Sigmaplot12.5软件制图ꎮ1.8㊀模型精度评价指标2019年数据用于建模ꎬ2018年数据用于外部独立验证ꎮ采用决定系数(coefficientofdeterminationꎬR2)㊁均方根误差(rootmeansquareerrorꎬRMSE)和归一化均方根误差(normalrootmeansquareerrorꎬNRMSE)评价和检验模型ꎮR2表示回归方程与实测值的拟合程度ꎬ其值越接近1效果越好ꎻRMSE用来衡量模拟值与实测值偏差ꎬ该值越小模拟效果越好ꎻNRMSE<10%模拟效果非常好ꎬ10%~20%效果较好ꎬ20%~30%效果411中国农业科技导报23卷可以接受ꎬ>30%效果差[22]ꎮ本研究评价不同估算模型时综合考虑R2㊁RMSE和NRMSE三个指标[23 ̄24]ꎮR2=1-ðni=1yi-xi()2ðni=1yi-y-()2(9)RMSE=㊀1nðni=1yi-xi()2(10)NRMSE=RMSEy-(11)式中ꎬyi为模拟值ꎬxi为实测值ꎬy-表示实测生长指标的平均值ꎬn为样本数ꎮ2㊀结果与分析2.1㊀不同施氮量对玉米冠层植被指数的影响由图1可见ꎬ两品种在N1㊁N2和N3处理下各植被指数均显著高于N0处理ꎮZD958在吐丝期和灌浆期各植被指数N2处理显著高于N1处理ꎬ说明N2处理下的玉米长势比N1处理好ꎬ在吐丝期N3处理下的GOSAVI㊁DVI和VIopt显著高于N1ꎬ而在灌浆期却显著低于N1(图1A㊁C)ꎮXY335在吐丝期各植被指数除DVI㊁VIopt和RVI外在各施肥处理间无显著性差异(图1B)ꎬ灌浆期的GOSAVI㊁EVI和RVI在N2和N3处理下显著高于N1处理(图1D)ꎮ2.2㊀植被指数与LAI和地上部生物量间的相关性分析LAI和地上部生物量与植被指数间的相关性分析结果(表2)表明ꎬ在吐丝期ꎬ两品种各植被指数与LAI和地上部生物量均呈极显著相关(P<0 01)ꎻ在灌浆期ꎬDVI与ZD958的LAI和地上部生物量相关性不显著(P>0.05)ꎬ而与XY335的LAI和地上部生物量显著相关(P<0.05)ꎻ灌浆期注:A㊁C分别为郑单958吐丝期和灌浆期ꎻB㊁D分别为先玉335吐丝期和灌浆期ꎮ图中不同字母表示氮肥处理间差异在P<0.05水平显著ꎮNote:AandCindicatesilkingandfillingstageofZD958ꎻBandDindicateskillingfillingstageofXY335.DifferentlettersinthefiguredenotessignificantdifferenceatP<0.05levelbetweennitrogenapplicationsusingDuncantest.图1㊀ZD958和XY335在不同施肥处理下的植被指数Fig.1㊀VegetationindexofZD958andXY335underdifferentnitrogenapplications5119期樊鸿叶等:基于无人机多光谱遥感的春玉米叶面积指数和地上部生物量估算模型比较研究表1㊀植被指数与LAI和地上部生物量的相关性Table1㊀CorrelationbetweenvegetationindexandLAIandabove ̄groundbiomass植被指数VegetationindexLAIZD958吐丝期Silkingstage灌浆期FillingstageXY335吐丝期Silkingstage灌浆期Fillingstage地上部生物量Above ̄groundbiomassZD958吐丝期Silkingstage灌浆期FillingstageXY335吐丝期Silkingstage灌浆期FillingstageNDVI0.908∗∗0.912∗∗0.872∗∗0.891∗∗0.844∗∗0.847∗∗0.843∗∗0.843∗∗GNDVI0.965∗∗0.911∗∗0.878∗∗0.932∗∗0.927∗∗0.842∗∗0.862∗∗0.878∗∗RVI0.957∗∗0.903∗∗0.939∗∗0.927∗∗0.890∗∗0.846∗∗0.891∗∗0.858∗∗GOSAVI0.967∗∗0.850∗∗0.902∗∗0.916∗∗0.912∗∗0.800∗∗0.868∗∗0.892∗∗EVI0.964∗∗0.917∗∗0.897∗∗0.903∗∗0.913∗∗0.850∗∗0.833∗∗0.854∗∗DVI0.956∗∗0.5650.927∗∗0.622∗0.860∗∗0.5730.843∗∗0.652∗VIopt0.953∗∗0.618∗0.754∗∗0.734∗∗0.852∗∗0.612∗0.769∗∗0.749∗∗㊀注:∗和∗∗表示在P<0.05水平和P<0.01水平相关性显著ꎮNote:∗and∗∗indicatesignificantcorrelationsatP<0.05andP<0.01levelsꎬrespectively.表2㊀玉米LAI模型精度对比Table2㊀ComparisonoftheaccuracyofdifferentmodelsformaizeLAI(n=12)品种Variety生育期Growthstage植被指数Vegetationindex指数IndexR2RMSE线性LinearR2RMSE对数LogarithmR2RMSE幂PowerR2RMSEZD958吐丝期Silkingstage灌浆期FillingstageNDVI0.9400.4840.9290.3100.9300.2610.9410.267GNDVI0.9430.2590.9310.2580.9310.2580.9430.259RVI0.9260.3010.9160.2850.9220.2760.9320.286GOSAVI0.9450.2530.9360.2500.9360.2490.9460.252EVI0.9410.2670.9290.2610.9300.2600.9420.266DVI0.9170.3250.9130.2900.9210.2770.9260.306VIopt0.9110.3390.9080.2990.9110.2930.9150.334NDVI0.8790.4220.8310.4150.8320.4150.8790.419GNDVI0.8580.4180.8300.4220.8300.4170.8760.418RVI0.8630.4460.8160.4340.8220.4270.8700.418GOSAVI0.7840.5650.7220.5330.7330.5230.7940.552EVI0.8880.620.8420.5090.8420.4830.8880.582XY335吐丝期Silkingstage灌浆期FillingstageNDVI0.8880.2530.8690.2600.8680.2610.8870.254GNDVI0.9220.2200.8940.2350.8900.2390.9190.227RVI0.9470.1790.3020.1900.9190.2050.9400.189GOSAVI0.9460.1750.9310.1890.9280.1940.9440.179EVI0.9500.2650.9380.1800.9360.1820.9450.169DVI0.8810.2600.8910.2390.9000.2300.8910.251VIopt0.8610.2960.8420.2730.8610.2670.8480.293NDVI0.7800.3420.7930.3320.7920.3330.7780.329GNDVI0.8630.2620.8690.2650.8670.2660.8610.261RVI0.8490.7100.8590.2750.8570.2770.8450.275GOSAVI0.8150.2980.8390.2940.8380.2940.8140.296EVI0.8040.3140.8160.3140.8150.3150.8020.313VIopt0.5040.4960.5390.5260.5400.4940.5090.815611中国农业科技导报23卷VIopt与ZD958的LAI和地上部生物量显著相关(P<0.05)ꎬ与XY335的LAI和地上部生物量极显著相关ꎻ其余植被指数与两品种在吐丝期和灌浆期的LAI和地上部生物量均呈极显著相关(P<0 01)ꎬ选择与LAI和地上部生物量分别达极显著水平的植被指数来构建估算模型ꎮ2.3㊀模型比较分析2.3.1㊀玉米LAI模型比较分析㊀从表3中ꎬ选取R2最大ꎬRMSE相对较小的植被指数和拟合方程(线性㊁对数㊁指数和幂)ꎬ如果出现R2相近则选择RMSE最小的拟合方程和植被指数ꎮ在吐丝期ꎬGOSAVI与ZD958的LAI的拟合幂函数效果最好(R2=0.946ꎬRMSE=0.252)ꎻRVI与XY335的LAI拟合的指数函数有较高的R2和最小RMSE(R2=0.947ꎬRMSE=0.179)ꎮ在灌浆期ꎬNDVI与ZD958的LAI拟合的幂函数R2相对较大ꎬRMSE最小(R2=0.879ꎬRMSE=0.419)ꎻGNDVI与XY335的LAI拟合的幂函数效果相对较好(R2=0.861ꎬRMSE=0.261)ꎮ2.3.2㊀玉米地上部生物量模型比较分析㊀吐丝期幂函数对两品种的地上部生物量估算效果较好ꎬ在灌浆期指数函数对两品种的地上部生物量估算效果好(表4)ꎮ在吐丝期ꎬGNDVI与ZD958地上部生物量拟合的幂函数效果最好ꎬR2和RMSE分别为0.867和0.488t hm-2ꎬGOSAVI与XY335地上部生物量拟合的幂函数估算效果最好(R2=0.791ꎬRMSE=0.463t hm-2)ꎮ在灌浆期ꎬEVI与ZD958地上部生物量拟合的指数函数建模效果较好(R2=0.772ꎬRMSE=1.477t hm-2)ꎬGNDVI与XY335地上部生物量拟合的指数效果表3㊀玉米地上部生物量模型精度对比Table3㊀Comparisonofaccuracyofdifferentmodelsformaizeabove ̄groundbiomass(n=12)品种Variety生育期Growthstage植被指数Vegetationindex指数IndexR2RMSE线性LinearR2RMSE对数LogarithmR2RMSE幂PowerR2RMSEZD958吐丝期Silkingstage灌浆期FillingstageNDVI0.8400.5890.8290.5310.8310.5280.8420.553GNDVI0.8650.4910.8600.52710.8600.4770.8670.488RVI0.8060.6200.7930.5850.8130.5560.8240.586GOSAVI0.8410.5490.8320.5260.8380.5170.8460.539EVI0.8440.5460.8350.5230.8370.5190.8460.542DVI0.7590.7080.7410.6530.7610.6280.7770.678VIopt0.7450.7300.7260.6730.7350.6610.7530.717NDVI0.7271.4970.7171.4980.7171.4980.7201.497GNDVI0.7581.5150.7101.5170.7091.5190.7571.517RVI0.7651.5030.7161.5000.7161.5010.7651.502GOSAVI0.6951.7270.6401.6890.6471.6740.7011.709EVI0.7721.4770.7231.4820.7231.4830.7721.479XY3335吐丝期Silkingstage灌浆期FillingstageNDVI0.7560.5080.7210.4990.7220.5000.7560.503GNDVI0.760.4780.7390.4840.7360.4870.7630.480RVI0.7860.4610.7630.4620.7560.4680.7790.463GOSAVI0.7910.4650.7620.460.7630.4610.7910.463EVI0.6940.5510.6670.5470.6670.5460.6940.550DVI0.6710.5780.6510.5600.6630.5500.6830.567VIopt0.5680.6570.5140.6410.5470.6370.5740.653NDVI0.8172.5670.7472.1120.7482.1030.8192.208GNDVI0.8522.0130.7861.9410.7881.9320.8552.192RVI0.7942.2870.7362.1540.7492.1000.8082.192GOSAVI0.8321.9610.7961.8920.7971.8890.8341.942EVI0.7892.2920.7292.1830.7312.1740.7912.273VIopt0.6644.7930.6992.3010.7042.2820.6702.6897119期樊鸿叶等:基于无人机多光谱遥感的春玉米叶面积指数和地上部生物量估算模型比较研究最好(R2=0.852ꎬRMSE=2.0132t hm-2)ꎮ2.4㊀LAI和地上部生物量模型验证为检验模型的可靠性ꎬ利用2018实测数据对所选模型进行验证ꎮ用RMSE和NRMSE两个指标来评价模型精度ꎬ从表5可以看出ꎬ在吐丝期ꎬ与ZD958拟合的LAI和地上部生物量模型估算精度较高ꎬRMSE和NRMSE分别为1.103㊁17.91%和1.485t hm-2㊁16.04%ꎬXY335拟合的模型精度相对较差ꎬRMSE和NRMSE分别为1.369㊁25.08%和1.730t hm-2㊁20.11%ꎮ在灌浆期ꎬZD958的LAI和地上部生物量估算模型精度较高ꎬRMSE和NRMSE分别为0.619ꎬ12.22%和2.482t hm-2ꎬ15 22%ꎻXY335的LAI估算精度较差(RMSE=1 008ꎬNRMSE=25.74%)ꎬ地上部生物量估算模型效果较好(RMSE=2.119t hm-2ꎬNRMSE=12.26%)ꎮXY335的LAI和地上部生物量模型估算效果比ZD958的差ꎬ但总体精度可以接受ꎬ能够用来估算LAI和地上部生物量ꎮ表4㊀LAI和地上部生物量模型验证结果Table4㊀ModelvalidationforLAIandabove ̄groundbiomass(n=12)生长指标Growthtarget品种Variety生育期Growthstage拟合模型FittedmodelRMSENRMSE/%叶面积指数LAIZD958XY335吐丝期SilkingstageLAI=11.447GOSAVI2.4311.10317.91灌浆期FillingstageLAI=11.626NDVI0.8430.62012.22吐丝期SilkingstageLAI=1.262e0.051RVI1.36925.08灌浆期FillingstageLAI=10.273GNDVI4.4371.10825.74地上部生物量Above ̄groundbiomassZD958XY335吐丝期SilkingstageDM=13.346GNDVI2.9321.48516.04灌浆期FillingstageDM=0.249e4.792EVI2.48215.22吐丝期SilkingstageDM=14.958GOSAVI2.4831.73020.11灌浆期FillingstageDM=0.046e7.323GNDVI2.11912.263㊀讨论吐丝期到灌浆期是玉米营养吸收和物质积累的高峰期[25]ꎬ从吐丝期开始玉米植株的养分向籽粒运转[26]ꎬ及时㊁准确的获取吐丝期和灌浆期玉米长势参量并据此进行调优栽培对提高玉米产量具有重要意义ꎮ本研究结果表明ꎬ对于ZD958而言ꎬ在吐丝期和灌浆期各植被指数在N2处理下显著高于N1处理ꎬ说明可能N2处理下玉米叶片叶绿素含量高于N1处理ꎬ冯宗会[27]在冬小麦上的研究表明ꎬ孕穗期施氮量270kg hm-2处理NDVI显著高于施氮量150kg hm-2处理ꎮ经相关性分析发现ꎬ同一植被指数与不同玉米品种间的LAI和地上部生物量的相关性差异较大ꎬ表明不同玉米品种的冠层光谱特征存在差异ꎬ这与张俊华等[28]研究结果一致ꎬ可能是由于品种间株型不同导致冠层光谱差异ꎮ多种植被指数与两玉米品种的LAI和地上部生物量均呈极显著或显著性相关水平ꎬ这与赵巧丽等[29]的研究结果基本吻合ꎮ本研究结果表明ꎬ吐丝期幂函数对ZD958的LAI估算效果最好ꎬ指数函数对XY335的LAI估算效果好ꎬ灌浆期幂函数对两品种的LAI估算效果最佳ꎮ吐丝期幂函数对两品种的地上部生物量估算效果较好ꎬ灌浆期指数函数对两品种的地上部生物量估算效果好ꎮ吐丝期ZD958的GNDVI对地上部生物量的估算效果最佳ꎬ而灌浆期GNDVI对XY335的估算效果最优ꎬShanahan等[30]研究也发现ꎬ抽雄期到灌浆期ꎬGNDVI的估产效应优于NDVI和TSAVIꎮ灌浆期EVI对ZD958的地上部生物量估算精度较高ꎬ韩文霆等[31]发现ꎬ单个生育期EVI对产量的估算效果最佳ꎮXY335在灌浆期GNDVI与LAI和地上部生物量估算效果均较好ꎬ说明在该时期GNDVI对XY335的LAI和地上部生物量的估算可通用ꎮ本研究一方面明确了不同玉米品种长势参量的光谱响应特征和最优预测模型ꎬ另一方面为下一步开展品种分辨率的光谱响应机理解析奠定了基础ꎮ本研究仍存在一些不足ꎬ玉米封垄后覆盖度较高ꎬ易导致植被指数出现 饱和 的现象ꎬ在今后的研811中国农业科技导报23卷究中应在多品种氮肥长期定位试验的基础上ꎬ研究光谱㊁空间和时相特征的融合方法ꎬ面向田间精确栽培的管理需求ꎬ构建适应于不同品种㊁不同生育期的玉米LAI和地上部生物量的遥感精确反演模型ꎮ致谢:感谢中国农业科学院作物科学研究所谢瑞芝研究员为本研究的顺利开展提供了完备的实验条件ꎮ参㊀考㊀文㊀献[1]㊀ZHANGCꎬKOVACSJM.Theapplicationofsmallunmannedaerialsystemsforprecisionagriculture:areview[J].PrecisionAgric.ꎬ13(6):693-712.[2]㊀李冰ꎬ刘镕源ꎬ刘素红ꎬ等.基于低空无人机遥感的冬小麦覆盖度变化监测[J].农业工程学报ꎬ2012ꎬ28(13):160-165.LIBꎬLIURYꎬLIUSHꎬetal..Winterwheatcoverchangemonitoringbasedonlow ̄altitudedroneremotesensing[J].Trans.Chin.Agric.Eng.ꎬ2012ꎬ28(13):160-165. [3]㊀CENHꎬWANLꎬZHUJꎬetal..DynamicmonitoringofbiomassofriceunderdifferentnitrogentreatmentsusingalightweightUAVwithdualimage ̄framesnapshotcameras[J].PlantMethodsꎬ2019ꎬ15(1):15-32[4]㊀吴炳方ꎬ张峰ꎬ刘成林ꎬ等.农作物长势综合遥感监测方法[J].遥感学报ꎬ2004ꎬ8(6):49-514.WUBFꎬZHANGFꎬLIUCLꎬetal..Acomprehensiveremotesensingmethodforcropgrowthmonitoring[J].J.Rem.Sens.ꎬ2004ꎬ8(6):498-514.[5]㊀王登伟ꎬ黄春燕ꎬ马勤建ꎬ等.棉花高光谱植被指数与LAI和地上鲜生物量的相关关系研究[J].新疆农业科学ꎬ2008ꎬ45(5):787-790.WANGDWꎬHUANGCYꎬMAQJꎬetal..StudyonthecorrelationbetweencottonhyperspectralvegetationindexandLAIandabovegroundfreshbiomass[J].J.XinjiangAgric.Sci.ꎬ2008ꎬ45(05):787-790.[6]㊀郑有飞ꎬGUOXꎬOLFERTOꎬ等.高光谱遥感在农作物长势监测中的应用[J].气象与环境科学ꎬ2007ꎬ30(01):10-16.ZHENGYFꎬGUOXꎬOlfertOꎬetal..Applicationofhyperspectralremotesensingincropgrowthmonitoring[J].Meteorol.Environ.Sci.ꎬ2007ꎬ30(01):10-16. [7]㊀GOELNꎬQINW.InfluencesofcanopyarchitectureonrelationshipsbetweenvariousvegetationindicesandLAIandFpar:acomputersimulation[J].RemoteSens.Rev.ꎬ10(4):309-347.[8]㊀林玉英ꎬ胡喜生ꎬ邱荣祖ꎬ等.基于Landsat影像的NDVI对植被与影响因子交互耦合的响应[J].农业机械学报ꎬ2018ꎬ49(10):212-219.LINYYꎬHUXSꎬQIURZꎬetal..ResponseofNDVIbasedonLandsatimagestovegetation ̄interactioncouplingwithimpactfactors[J].Trans.Chin.Soc.Agric.Machineryꎬ2018ꎬ49(10):212-219.[9]㊀米兆荣ꎬ张耀生ꎬ赵新全ꎬ等.NDVI和EVI在高寒草地牧草鲜质量估算和植被动态监测中的比较[J].草业科学ꎬ2010ꎬ27(6):13-19.MIZRꎬZHANGYSꎬZHAOXQꎬetal..ComparisonofNDVIandEVIinfreshqualityestimationandvegetationdynamicsmonitoringinalpinegrassland[J].Pratacult.Sci.ꎬ2010ꎬ27(6):13-19.[10]㊀左丽君ꎬ张增祥ꎬ董婷婷ꎬ等.MODIS/NDVI和MODIS/EVI在耕地信息提取中的应用及对比分析[J].农业工程学报ꎬ2008ꎬ24(3):167-172.ZUOLJꎬZHANGZXꎬDONGTTꎬetal..ApplicationandcomparativeanalysisofMODIS/NDVIandMODIS/EVIinarablelandinformationextraction[J].Trans.Chin.Agric.Engin.ꎬ2008ꎬ24(3):167-172.[11]㊀JORDANFC.Derivationofleaf ̄areaindexfromqualityoflightontheforestfloor[J].Ecologyꎬ1969ꎬ50(4):663-666. [12]㊀杨贵军ꎬ李长春ꎬ于海洋ꎬ等.农用无人机多传感器辅助小麦育种信息获取[J].农业工程学报ꎬ2015ꎬ31(21):184-190.YANGGJꎬLICCꎬYUHYꎬetal..Agriculturaldronemulti ̄sensorremotesensingaidswheatbreedinginformationacquisition[J].Trans.Chin.Agric.Eng.ꎬ2015ꎬ31(21):184-190.[13]㊀王亚杰.基于无人机多光谱遥感的玉米叶面积指数监测方法研究[D].陕西杨凌:西北农林科技大学ꎬ硕士学位论文ꎬ2018.WANGYJ.Astudyonthemonitoringmethodofmaizeleafareaindexbasedonmultispectralremotesensingbydrones[D].ShaanxiYangling:NorthwestUniversityofAgricultureꎬForestryandScienceꎬMasterDissertationꎬ2018.[14]㊀高林ꎬ杨贵军ꎬ王宝山ꎬ等.基于无人机遥感影像的大豆叶面积指数反演研究[J].中国生态农业学报ꎬ2015ꎬ23(7):868-876.GAOLꎬYANGGJꎬWANGBSꎬetal..Aninversionstudyofsoybeanleafareaindexbasedondroneremotesensingimages[J].Chin.J.EcolAgric.ꎬ2015ꎬ23(7):868-876. [15]㊀孙诗睿ꎬ赵艳玲ꎬ王亚娟ꎬ等.基于无人机多光谱遥感的冬小麦叶面积指数反演[J].中国农业大学学报ꎬ2019ꎬ24(11):51-58.SUNSRꎬZHAOYLꎬWANGYJꎬetal..Inversionofwinterwheatleafareaindexbasedondronemultispectralremotesensing[J].J.ChinaAgric.Univ.ꎬ2019ꎬ24(11):51-58. [16]㊀邹楠ꎬ杨文杰ꎬ肖春华ꎬ等.种植密度对玉米冠层高光谱特征的响应研究[J].石河子大学学报(自然科学版)ꎬ2017ꎬ35(6):687-692.ZOUNꎬYANGWJꎬXIAOCHꎬetal..Astudyontheresponseofplantingdensitytohyperspectralcharacteristicsofmaizecanopy[J].J.ShiheziUniv.(Nat.Sci.)ꎬ2017ꎬ35(6):687-692.[17]㊀肖武ꎬ陈佳乐ꎬ笪宏志ꎬ等.基于无人机影像的采煤沉陷区玉米生物量反演分析[J].农业机械学报ꎬ2018ꎬ49(08):169-180.XIAOWꎬCHENJLꎬGAOHZꎬetal..Cornbiomassinversionandanalysisofcoalminingsubsidencezonebasedondroneimagery[J].Trans.Chin.Soc.Agric.Machineryꎬ2018ꎬ49(08):169-180.[18]㊀庄东英ꎬ李卫国ꎬ武立权.冬小麦生物量卫星遥感估测研究[J].干旱区资源与环境ꎬ2013ꎬ27(10):158-162.9119期樊鸿叶等:基于无人机多光谱遥感的春玉米叶面积指数和地上部生物量估算模型比较研究ZHUANGDYꎬLIWGꎬWULQ.Astudyonremotesensingestimationofwinterwheatbiomassbysatellite[J].J.AridLandResour.Environ.ꎬ2013ꎬ27(10):158-162.[19]㊀BUSCHMANNCꎬNAGELE.INvivospectroscopyandinternalopticsofleavesasbasisforremotesensingofvegetation[J].Int.J.RemoteSens.ꎬ14(4):711-722.[20]㊀RONDEAUXGꎬSTEVENMꎬBARETF.Optimizationofsoil ̄adjustedvegetationindices[J].RemoteSens.Environ.ꎬ55(2):95-107.[21]㊀马文玉.基于UAV遥感的冬小麦氮素营养诊断[D].河南焦作:河南理工大学ꎬ硕士学位论文ꎬ2017.MAWY.NutritionaldiagnosisofnitrogeninwinterwheatbasedonUAVremotesensing[D].HenanJiaozuo:HenanUniversityofScienceandTechnologyꎬMasterDissertationꎬ2017.[22]㊀JAMIESONPDꎬPORTERJRꎬWILSONDR.Atestofcomputersimulationmodelarcwheat1onwheatcropsgrowninNewZealand[J].FieldCropsRes.ꎬ1991ꎬ27(4):337-350. [23]㊀BENDIGJꎬBOLTENAꎬBOLTENS.Stimatingbiomassofbarleyusingcropsurfacemodels(CSMs)derivedfromUAV ̄BasedRGBimaging[J].RemoteSens.ꎬ2014ꎬ6(11):10395-10412.[24]㊀YUEJꎬHAIKUANFꎬGUIJUNYꎬetal..Acomparisonofregressiontechniquesforestimationofabove ̄groundwinterwheatbiomassusingnear ̄surfacespectroscopy[J/OL].RemoteSens.ꎬ10(2):66[2021-05-12].https://doi.org/10.3390/rs10010066.[25]㊀刘冰峰ꎬ李军ꎬ贺佳ꎬ等.基于高光谱植被指数的夏玉米地上干物质量估算模型研究[J].农业机械学报ꎬ2016ꎬ47(03):254-262.LIUBFꎬLIJꎬHEJꎬetal..Amodelstudyontheestimationofdrymattermassonsummercorngroundbasedonhyperspectralvegetationindex[J].Trans.Chin.Soc.Agric.Machineryꎬ2016ꎬ47(03):254-262.[26]㊀申丽霞ꎬ王璞ꎬ张软斌.施氮对不同种植密度下夏玉米产量及子粒灌浆的影响[J].植物营养与肥料学报ꎬ2005ꎬ11(03):314-319.SHENLXꎬWANGPꎬZHANGRB.Effectofnitrogenapplicationonsummercornyieldandseedgroutingatdifferentplantingdensities[J].J.PlantNutr.Fert.ꎬ2005ꎬ11(03):314-319.[27]㊀冯宗会.不同水氮管理下NDVI监测及优化施氮研究[D].北京:首都师范大学ꎬ硕士学位论文ꎬ2012.FENGZH.NDVImonitoringandoptimizationofnitrogenapplicationunderdifferentwaterandnitrogenmanagement[D].Beijing:CapitalNormalUniversityꎬMasterDissertationꎬ2012.[28]㊀张俊华ꎬ张佳宝ꎬ贾科利.不同品种夏玉米光谱特征差异及其农学参量相关性研究[J].土壤通报ꎬ2010ꎬ41(2):287-293.ZHANGJHꎬZHANGJBꎬJIAKL.Differencesinthespectralcharacteristicsofsummermaizeandtheircorrelationwithagronomicparameters[J].Chin.J.SoilSci.ꎬ2010ꎬ41(2):287-293.[29]㊀赵巧丽ꎬ郑国清ꎬ段韶芬ꎬ等.基于冠层反射光谱的玉米LAI和地上干物重估测研究[J].华北农学报ꎬ2008ꎬ23(1):219-222.ZHAOQLꎬZHENGGQꎬDUANSFꎬetal..AstudyofcornLAIandabovegrounddrymatterreassessmentbasedoncanopyreflectionspectra[J].NorthChin.J.Agric.Sci.ꎬ2008ꎬ23(1):219-222.[30]㊀SHANAHANJFꎬSCHEPERSJSꎬFRANCISDDꎬetal..Useofremote ̄sensingimagerytoestimatecrongrainyield[J].Agron.J.ꎬ2001ꎬ93(3):583.[31]㊀韩文霆ꎬ彭星硕ꎬ张立元ꎬ等.基于多时相无人机遥感植被指数的夏玉米产量估算[J].农业机械学报ꎬ2020ꎬ51(1):1-15.HANWTꎬPENGXSꎬZHANGLYꎬetal..Summercornyieldestimationbasedonmulti ̄temporaldroneremotesensingvegetationindex[J].Trans.Chin.Soc.Agric.Machineryꎬ2020ꎬ51(1):1-15.(责任编辑:温小杰)021中国农业科技导报23卷。

基于经验模型和机理模型的杉木人工林生物量估测对比研究


The determination coefficient ( R 2 ) of the linear regression equations between the measured and estimated values of each
factor ranged from 0.65 to 0.96. The R 2 of diameter at breast height ( DBH) and tree height was higher than 0.92. The mean
具有较好的表现,但各具特点和局限性。 经验模型参数较易获得,且经验模型预测生物量、林分胸高断面积和林分平均树高的
R2 、MRE 均优于机理模型;但模型对于建模数据内的评价效果较好,对于建模数据外的应用具有局限性,即经验模型更适合模
拟生长期间的某一阶段的林分生物量。 机理模型虽然需要的参数较多,但是考虑了生态学原理,弥补了经验模型的不足,可较
lanceolata plantations based on empirical model and mechanism model.Acta Ecologica Sinica,2024,44(8) :3502⁃ 3514.
基于经验模型和机理模型的杉木人工林生物量估测对
比研究
李佳怡1 ,阮红玉2 ,邱思玉1 ,梁瑞婷1 ,朱兆廷1 ,文 毅1 ,王成雨1 ,王轶夫1,∗ ,孙玉军1
1 北京林业大学森林资源和环境管理国家林业和草原局重点实验室,北京 100083
2 吉水县石阳林场,吉安 331600
摘要:为了揭示间伐干扰下杉木人工林生物量的变化规律,研究利用江西省吉水县石阳林场的 36 块杉木人工林样地的实测数
据和研究区气候数据,通过基于经验的引入地位指数( SI) 的生物量生长方程组和基于机理的 3⁃PG 模型,模拟并预估林分生物

测绘技术中的树木森林测量方法

测绘技术中的树木森林测量方法近年来,随着人们对环境保护意识的提高,森林资源的测量与管理变得愈发重要。

测绘技术在森林资源测量中起着关键作用,尤其是树木森林测量方法的应用。

本文将探讨测绘技术中的树木森林测量方法,展示其在实际应用中的国内外现状和发展趋势。

一、激光雷达测量方法激光雷达技术是树木森林测量领域最为常用的方法之一。

通过向森林区域发射激光束,激光雷达可以测量出激光束与树木之间的距离和高度差。

通过对激光束的反射和回波时间的测量,可以获得树木的高度、形态、密度和位置信息。

激光雷达测量方法具有无侵入性、高精度、高效率等优点,被广泛应用于森林调查和林业资源管理。

国际上已有很多研究通过激光雷达进行树木高度、冠层密度和干体积等参数测量。

二、无人机遥感测量方法随着无人机技术的飞速发展,无人机遥感测量方法在树木森林测量领域也逐渐崭露头角。

通过搭载激光雷达、高分辨率相机等设备的无人机,可以对森林进行高精度、高分辨率的测量和影像获取。

无人机遥感测量方法弥补了传统测量方法的局限性,使得数据获取更加快速、准确。

同时,无人机在树木森林测量中的应用还包括树种分类、林分结构和生物量估算等方面,为森林资源的管理提供了更多的数据支持。

三、全站仪测量方法全站仪是一种通过光学和测距技术进行树木测量的装备。

全站仪测量方法主要通过测量树木的高度、直径和树冠等参数来评估森林资源。

全站仪具有高精度、易操作和快速测量等优点,广泛应用于林业和测绘领域。

在树木森林测量中,全站仪可以用于编制树木清查表、制定采伐计划和评估生态环境等方面,对于森林资源的可持续管理起着重要作用。

四、遥感影像解译方法遥感影像解译方法是树木森林测量中另一重要的测量手段。

通过获取卫星或航空摄影的遥感影像,利用数字图像处理和解译技术,可以获得树木的空间分布、覆盖率和种类等信息。

遥感影像解译方法可以大范围、多时态的获取树木森林数据,为森林资源管理提供了方便快捷的工具。

此外,遥感影像解译方法结合地理信息系统,还可以实现对树木森林资源进行大规模的空间分析和决策支持。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

森林地上生物量遥感估测研究进展娄雪婷;曾源;吴炳方【摘要】Above -ground biomass of forest has great research and application value in the forest ecological system.There are mainly three types of models for estimating above - ground biomass of forest, i. e. , forest measuring method, remote sensing method and integrated method. Remote sensing technique has become an important means for obtaining above - ground biomass of forest at the regional scale. There are mainly four types of remote sensing models, namely empirical, ANN, physical and NPP based models. This paper has analyzed and discussed the present methods for estimating above - ground biomass of forest based on remote sensing as well as their advantages and disadvantages. Finally, this paper points out that the integrated method combining remote sensing technique and forest succession model can be generally used to estimate above - ground biomass of forest at the regional scale in future.%森林生物量是衡量生态系统生产力的重要指标,也是研究森林生态系统物质循环的重要基础,其估测方法可以分为传统地面实测法、遥感监测法和综合模型法.随着生物量估测从样地研究发展到区域应用,空间尺度的增大导致宏观资料和参数的获取存在很多困难.在深入分析目前应用遥感技术估算森林生物量的方法及原理基础上,系统评述了统计模型、物理模型、基于植被净初级生产力(NPP)模型的估算方法以及综合模型法的优缺点,分析了各种方法在不同森林植被及遥感数据源下的适用性及不确定性,探讨了此领域的研究方向.【期刊名称】《国土资源遥感》【年(卷),期】2011(000)001【总页数】8页(P1-8)【关键词】森林地上生物量;遥感;模型估算【作者】娄雪婷;曾源;吴炳方【作者单位】中国科学院遥感应用研究所农业与生态遥感研究室,北京,100101;中国科学院遥感应用研究所农业与生态遥感研究室,北京,100101;中国科学院遥感应用研究所农业与生态遥感研究室,北京,100101【正文语种】中文【中图分类】TP79森林分布面积广,在陆地生物圈层和大气圈层之间的碳循环交换过程中起着重要作用[1]。

森林生物量约占陆地生态系统生物量的 90%[2],是整个森林生态系统运行的能量基础和营养物质来源,对陆地生物圈和其他地表过程有着重要影响。

森林生物量不仅反映森林的生态演替和人为干扰过程,还能充分体现森林生态系统环境的质量,是衡量森林生态系统生产力、生态系统结构功能的最好指标[3]。

森林生物量按来源的部位不同可分为根、干、枝和叶生物量。

其中,干、枝、叶生物量合称为森林地上生物量。

森林地上生物量因其寿命长、体积大而长期、大幅度地影响碳库,是研究森林碳循环、评估森林碳收支的基础,对全球气候变化研究具有重要意义。

本文在总结传统地面实测法的基础上,详细探讨了基于遥感技术估算森林地上生物量的方法,分别讨论了各种方法的适用性及其不确定性,特别提出了综合模型法,为选择估算森林地上生物量方法提供参考。

提取森林地上生物量的研究工作经过简单测定阶段,从 20世纪 60年代起,以统计模型法研究森林地上生物量渐呈主流,即以易测因子作为自变量,生物量为因变量,通过回归分析方法建立起相对生长方程,以此来推算、预测整体的生物量[4-8]。

此后,随着国际生物学计划和人与生物圈计划的开展,遥感监测技术成为研究森林地上生物量的新工具。

虽然该技术在探索影响植物光谱特征的植物生理本质方面还有待完善,但其诸多优势使其成为当前及今后的研究热点。

1.1 传统地面实测法传统地面实测法以测量研究区内单位面积林木重量为主,主要有皆伐实测法、随机抽样法、平均木生物量法、材积源生物量法、相对生长模型法及气体交换法6种。

各种方法的原理、优缺点和适用性如表 1所示。

然而,传统地面实测法通常只能获得点上的数据,且劳动强度大、成本高,不利于研究生物量的空间分布和变化,且在实际应用中存在着林业测量仪器误差、样地空间分布的代表性不足等问题,这些均增加了研究区域尺度森林地上生物量的难度。

尽管如此,将传统地面实测法获得的数据作为生物量遥感提取模型的输入参数来验证数据,仍十分重要。

近年来,传统地面实测法又不断与遥感技术、空间定位技术结合,为遥感监测森林地上生物量奠定了基础。

王志杰[11]尝试利用全站仪和掌上电脑获取监测样地内每棵树木的定位数据,同时采集树木物种信息、形态信息及立地环境数据,建立样地植物监测数据库系统,实时、详细的样地植物信息更易于与高空间分辨率遥感数据接轨。

陈春祥[12]采用 FoxPro及 ArcView GIS改进了大范围林业抽样样地布点方案的修改方法,提高了抽样设计的工作效率。

3S技术发展了传统地面实测法的样地布设和抽样调查方法,为其注入了新活力。

1.2 基于遥感的监测方法遥感作为日趋成熟的对地观测技术,不仅可以为快速、定量获取区域尺度森林地上生物量及其变化提供大范围、多尺度的数据,更可实现对地表参数的实时监测,适用于分析和评价森林生物量的时空格局和变化趋势。

目前主要有统计模型、物理模型和综合模型 3类方法。

1.2.1 统计模型法1.2.1.1 统计模型统计模型法是根据实验区内样地生物量与遥感图像上植被的物理参数、反射光谱特征或不同通道雷达数据间存在的相关关系进行回归拟合,按像元计算生物量的方法。

研究表明,基于植被特殊的光谱特性,多波段光学遥感信息组合成的植被指数能包含90%以上的植被信息[14],以光学遥感的不同波段组合提取生物量已经被广泛应用[14-16]。

高光谱遥感数据能记录下地物完整的光谱曲线,为我们提供了更多的可选波段来表征植被信息。

其中“红边”是植被高光谱曲线最明显的特征,相比植被指数不易受土壤和大气环境的影响[17]。

激光雷达 (LiDAR)在林木高度测量与林分垂直结构信息获取方面具有其他遥感技术无可比拟的优势。

何祺胜[18]等采用高密度L i-DAR小脚印点云数据建立冠层高度模型,估测单株木结构参数树的位置、树高、冠幅,然后采用多元逐步回归分析法建立 LiDAR估测的平均树高、冠幅等与实测森林地上生物量之间的关系,估测方程精度较高,R2=0.738。

文献 Holmgren[19]以 L iDAR数据估测树高、冠幅等参数,计算 10m半径样地的蓄积量,估测精度达到应用要求。

随着遥感数据源的多元化发展,为了集合不同数据源的优势以提高森林生物量的估算精度,越来越多的学者致力于用数学方法和技术工具融合不同传感器、不同波段的数据。

例如高分辨率数据可校正低分辨率数据的误差,提高森林生物量的估算精度。

张宏斌[20]检验了 TM、MOD IS、NOAA数据之间的精度校正关系,并以气候因素 (气温、降雨)和遥感植被指数之间的关系模型监测了 1982~2006年内蒙古草原植被的变化趋势。

应用统计模型的突出优点是可以宏观、连续地监测植被生物量。

在遥感影像光谱或空间分辨率足够高的情况下,此类模型简单易用。

但它的局限性也不可忽视:一方面,由于大气、遥感器定标、太阳几何照明、土壤湿度、亮度和颜色在不同的区域每时每刻变化着,统计关系的适用性常常不稳定[21-23],需要根据应用区域的不同来选择适合的数据源及其波段组合[24,25];另一方面,由于树木的生长受温度、光照及水分等因子限制,不同季节的树干、树枝及树冠生物量之间的统计关系也会发生变化。

总之,基于统计关系估算林地整体地上生物量,增加了结果的不确定性。

1.2.1.2 人工神经网络模型尽管不同的波段组合与森林地上生物量在理论上有明显的相关性,但森林环境的复杂性使各种统计函数无法表达所有的不确定性。

应用人工神经网络估算森林地上生物量,开辟了一条方便快捷且可靠的途径。

人工神经网络 (ANN)系统是以模拟人脑神经系统的结构和功能为基础而建立的一种数据分析处理系统。

在进行知识获取时,由研究者提供样本和相应的解,通过特定的学习算法对样本进行训练,通过网络内部自适应算法不断修改权值分布以达到应用要求[26]。

目前,提取森林生物量神经网络模型很多,如线性网络、回归网络、后向 (BP)神经网络、感知器神经网络 (MLP)及径向基函数神经网络 (RBF)等 [27]。

有实验表明,在生物量估测与实测数据的计算中,人工神经网络法估测的相对误差较小,约小于15%[29],引入海拔、坡度、坡向和立地类型作为神经网络模型的自变量,有助于提高生物量的估算精度[29]。

同时,在多元遥感协同估测生物量方面,韩爱惠[30]结合地理气象因子,采用多元回归和神经网络的方法,研建了以 CBERSWFI数据和 MODIS数据为主要信息源的森林生物量估测模型,估测精度达到了 85%以上。

在森林生物量与遥感因子相关性不显著的情况下,神经网络是建立森林生物量非线性遥感模型的可靠方法。

但其缺点是“黑箱操作”,模型缺乏生物物理意义。

1.2.2 物理模型法物理模型法是指借助二向性反射与生物量之间的关系由遥感信息反演估算生物量的方法。

自然界绝大部分的物体均具有各向异性的反射特性,投射到地物表面上的辐射能量往往由来自太阳的直射辐射与天空散射辐射两部分组成,传感器接收的辐射亮度受辐射环境的影响,称之为二向性反射,可用二向性反射率分布函数 BRDF来描述[31]。

物理模型有明确的物理意义,模型稳定度和普适性好。

目前,应用于提取森林地上生物量的模型主要有辐射传输模型 (Radiative Transfer Model)和几何光学模型(Geometric Optical Model)。

相关文档
最新文档