Estimation of total nitrogen content in sugarbeet leaves under drip irrigation based on hyperspectral characteristic parameters and vegetation index
LI Zong-Fei1, SU Ji-Xia1, FEI Cong1, LI Yang-Yang1, LIU Ning-Ning1, DAI Yu-Xiang1, ZHANG Kai-Xiang1, WANG Kai-Yong1, FAN Hua,1,*, CHEN Bing,2,*通讯作者:
收稿日期:2019-03-20接受日期:2019-12-26网络出版日期:2020-01-17
基金资助: |
Received:2019-03-20Accepted:2019-12-26Online:2020-01-17
Fund supported: |
作者简介 About authors
E-mail:
摘要
关键词:
Abstract
Keywords:
PDF (653KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
李宗飞, 苏继霞, 费聪, 李阳阳, 刘宁宁, 戴宇祥, 张开祥, 王开勇, 樊华, 陈兵. 基于高光谱数据的滴灌甜菜叶片全氮含量估算[J]. 作物学报, 2020, 46(4): 557-570. doi:10.3724/SP.J.1006.2020.94045
LI Zong-Fei, SU Ji-Xia, FEI Cong, LI Yang-Yang, LIU Ning-Ning, DAI Yu-Xiang, ZHANG Kai-Xiang, WANG Kai-Yong, FAN Hua, CHEN Bing.
氮素是甜菜生长必不可少的营养元素, 是构成植物体内蛋白质、核酸、叶绿素以及多种酶的主要成分, 直接反映植被营养状况, 在植物光合作用与生长发育中发挥着重要作用[1,2], 因此如何根据作物生长发育状况来判断作物适宜的追肥量, 如何运用不同氮素运筹管理以及氮素水平来合理施肥[3], 快速、无损监测甜菜生长状况, 制定适时适宜施肥方案, 对甜菜生产具有重要意义。
传统氮素营养诊断手段以分析测试为主, 因时效性差, 难以在大面积农业生产中有效应用[4]。与传统方法相比, 高光谱遥感可以得到作物生长发育中生理生态与营养特征等微小特性变化, 实现对目标实时、快速监测[5]。近年来, 国内外****研究了在不同施氮处理条件下多种作物的光谱反射特性, 发现了作物不同施氮处理条件下光谱反射率的差异性, 并可用具体波段与波段组合来估测氮素含量。Shibayama等[6]在研究水稻叶片光谱特性时发现, 缺氮水稻叶片与正常水稻叶片的光谱特征显著不同。Smith等[7]研究发现对于同一种植物在不同时期氮素的最佳反演波段存在变化。Stone等[8]基于671 nm和780 nm两个波段反射率组合的氮光谱指数对小麦植株的全氮含量进行了估算。Thomas等[9]研究发现550~675 nm波段与氮素含量具有较好相关性, 可用来估测辣椒的氮素含量。Tian等[10]基于不同氮素水平和多个田间试验, 提出3个蓝光波段构建的植被指数R434/(R496×R401)对叶层氮浓度具有较好预测效果。Wang等[11]对归一化植被指数(NDVI)改进, 提出三波段指数(R924-R703+2×R423)/(R924+ R703-2×R423), 有效提高氮素反演精度。Menesatti等[12]利用高光谱数据建立了锦橙LNC的偏最小二乘法估算模型。Zhang等[13]建立了近红外反射光谱与苹果LNC的多元线性回归方程。
在我国, 利用高光谱技术研究作物氮素营养起步较晚, 但发展较快, 目前已经在小麦[14,15]、水稻[16,17]、棉花[18,19]等作物上得到广泛应用, 赵春江等[20]对小麦不同品种、肥水条件下光谱参数变化规律及与植株生化组分关系研究提出, 可用红边振幅和近红外平台振幅来反演小麦叶片全氮含量。王树文等[21]采用回归分析方法对寒地水稻氮素含量建模分析, 发现用特征光谱参数SI (866)、DI (730、715)、DDI (730、715、685)、DDI (866、685、558)可以较好预测氮素含量。陈兵等[22]研究了基于光谱红边参数的棉花黄萎病叶片氮素含量诊断, 发现以Area672为自变量建立的棉叶LNC的诊断模型能较好地诊断黄萎病棉叶LNC。
尽管前人利用光谱变量估测作物农学参数方面已经开展了大量的研究, 但在甜菜上涉及不多, 相关甜菜氮素营养光谱诊断研究较少, 且氮素估算模型方法较为单一。本文以滴灌甜菜为研究对象, 通过筛选甜菜叶片全氮含量敏感波段, 构建叶片全氮含量敏感参数与植被指数, 建立并对比叶片全氮含量的一元线性、二次项、幂函数以及指数函数估算模型, 以期为滴灌甜菜叶片全氮含量监测和氮肥精准管理提供支持。
1 材料与方法
1.1 试验设计
试验于2017年4月至2017年11月在新疆石河子大学农学院试验站(44°31'N, 88°06'E)进行。试验区域耕层土壤为灌溉灰漠土, pH 7.3, 含有机质13.27 g kg-1、全氮0.76 mg kg-1、速效磷22.56 mg kg-1、速效钾157.97 mg kg-1。共设置N0 (0 kg hm-2)、N75 (75 kg hm-2)、N150 (150 kg hm-2)、N225 (225 kg hm-2) 4个施氮水平。施氮比例按照M0 (0∶0∶0)、M1 (10∶0∶0)、M2 (7∶3∶0)、M3 (5∶3∶2) 3种运筹模式, 即分别在甜菜叶丛快速生长期、块根膨大期与块根糖分积累期按比例施入氮肥。采用随机区组设计, 重复3次, 供试甜菜品种为Beta356, 播种日期为4月15日, 试验田内设置30个小区, 小区面积12 m2, 甜菜播种行距50 cm, 株距20 cm。基肥为重过磷酸钙319.5 kg hm-2, 硫酸钾206 kg hm-2。其他按照大田统一管理。1.2 甜菜冠层光谱数据测定
采用美国Analytical Spectral Device (ASD)公司生产的FileSpecPro FR2500型背挂式野外高光谱辐射仪测定甜菜植株冠层高光谱反射率。该光谱仪波段范围为350~2500 nm, 采样间隔为1 nm, 光谱仪视场角25°, 测量时探头垂直向下, 距离冠层垂直高度约0.5 m, 在测量前均进行白板校正。选择晴朗无云天气, 于当地北京时间12:00—15:00分别在甜菜叶丛快速生长期、块根膨大期、糖分积累期测定光谱反射率, 记录每个光谱测定点5条光谱数据, 每小区测定3次, 取其平均值作为该小区的光谱反射率值。1.3 全氮含量测定
采用凯氏定氮法测定叶片全氮含量。取光谱测定区域叶片, 105℃杀青30 min后, 调至75℃烘干至恒重, 将烘干样品粉碎过60目筛。用分析天平秤取0.5000 g叶片样品放入消煮管, 采用H2O2-浓H2SO4消化处理后, 将冷却的消煮液置于定氮仪中蒸馏, 利用0.01 mol L-1(1/2H2SO4)标准溶液滴定馏出液由蓝绿色至刚变为紫红色, 记录所用标准酸体积, 根据以下公式计算植株全氮含量。叶片全氮(N)量(g kg-1) =
式中, C是0.01 mol L-1(1/2H2SO4)标准溶液浓度, V和V0分别是滴定试液时所用酸标准溶液的体积与滴定空白时所用酸标准溶液的体积(mL), 0.014为氮的摩尔质量kg mol-1, TS为分取倍数, 即定容体积与分取体积之比, 1000为kg与g之间的换算倍数, m为烘干株样的质量(g)。
1.4 数据分析
光谱原数据提取和预处理采用FileSpecPro FR2500型背挂式野外高光谱辐射仪自带软件ViewSpec Program, 用Microsoft Excel 2013与SPSS 12.0软件统计分析农学参数数据和光谱数据, 并利用Origin 2018软件制图, 采用MATLAB软件提取甜菜冠层光谱数据的特征参数, 从中提取常用光谱特征“三边”参数; 通过选择已有光谱植被指数作对比, 结合与本试验中全氮含量相关性较好的敏感性波段, 构建比值(RVI)、差值(DVI)、归一化比值(NDVI)与红边归一化差异(NDI)光谱植被指数(表1)。Table 1
表1
表1高光谱特征参数及植被指数定义
Table 1
光谱特征参数与植被指数 Spectral characteristic parameter and Spectral index | 定义 Definition |
---|---|
红边斜率 Red edge slope (Dr) | 680-760 nm内最大一阶微分光谱值 Maximum first order differential spectrum in 680-760 nm |
红边位置 Red edge position (λr) | 680-760 nm内最大的一阶微分光谱值对应的波长 The wavelength corresponding to the largest first-order differential spectrum in 680-760 nm |
红边面积 Red edge area (SDr) | 680-760 nm内一阶微分光谱值的总和 The sum of the first-order differential spectral values in 680-760 nm |
黄边斜率 Yellow edge slope (Dy) | 560-640 nm内最大一阶微分光谱值 Maximum first order differential spectrum in 560-640 nm |
黄边位置 Yellow edge position (Λy) | 560-640 nm内最大一阶微分光谱值对应的波长 The wavelength corresponding to the largest first-order differential spectrum in 560-640 nm |
黄边面积 Yellow edge area (SDy) | 560-640 nm内一阶微分光谱值的总和 The sum of the first-order differential spectral values in 560-640 nm |
蓝边斜率 Blue edge slope (Db) | 490-530 nm内最大一阶微分光谱值 Maximum first order differential spectrum in 490-530 nm |
蓝边位置 Blue edge position (λb) | 490-530 nm内最大一阶微分光谱值对应的波长 The wavelength corresponding to the largest first-order differential spectrum in 490-530 nm |
蓝边面积 Blue edge area (SDb) | 490-530 nm内一阶微分光谱值的总和 The sum of the first-order differential spectral values in 490-530 nm |
比值植被指数 Ratio vegetation index (RVI) | Rλ1/Rλ2 |
差值植被指数 Difference vegetation index (DVI) | Rλ1- Rλ2 |
归一化比值植被指数 Normalized ratio vegetation index (NDVI) | (RNIR-Rred)/(RNIR+Rred) |
红边归一化差异指数 Red edge normalized difference index (NDI) | (DRλ1-DRλ2)/(DRλ1+DRλ2) |
新窗口打开|下载CSV
通过分析光谱特征参数及植被指数与甜菜冠层叶片全氮含量相关性, 建立全氮含量估测模型。本文主要针对甜菜叶丛快速生长期至糖分积累期3个时期进行测定, 甜菜在叶丛快速生长期后基本达到土壤背景全覆盖, 且甜菜自叶丛快速生长期至糖分积累期冠层结构无明显变化, 因此本文选择总样本量用于甜菜叶片全氮含量估算与检验, 根据测定样本数量大小, 将采集样本随机分为2组, 其中建模样本75个, 模型检验样本26个, 通过对比各模型的建模决定系数(R2)、均方根误差(RMSE)以及相对误差(RE)筛选出最优模型用以研究新疆滴灌条件下甜菜冠层叶片全氮含量的估算。其中R2用以判断模型拟合优度, RE与RMSE用于检验估测模型的可靠程度。
RE =
RMSE =
式中, yi为试验测定全氮含量真实值,
2 结果与分析
2.1 不同氮素水平下甜菜冠层光谱特征分析
通过对不同氮素水平下光谱数据进行平均处理, 由原始光谱曲线图与一阶导数光谱曲线图分析可知(图1), 不同氮素水平下甜菜原始冠层光谱反射率波形趋势相似, 在波段350~670 nm范围内有明显的反射峰与吸收谷, 光谱波形呈现先增高后降低的趋势, 在550 nm附近形成反射小高峰, 均在670 nm附近形成吸收谷, 在波段670~760 nm波段波形变化较为剧烈, 光谱反射率随着波长的增加急剧上升, 形成较大陡坡状波形; 在近红外波段780~1300 nm形成近红外反射平台, 且具有明显的波形变化。施氮处理N75、N150、N225在550 nm处较CK分别降低了7.5%、28.0%、23.0%, 说明各施氮处理相对CK均可提高甜菜植株冠层叶片对红、蓝光的吸收, 提高植株的光合效率, 且施氮处理N150比N75与N225具有更好地提高植株光合效率的作用; 在近红外波段(780~1300 nm)各施氮处理较CK分别增加6.14%~12.6%、14.2%~16.4%、6.82%~18.0%, 各处理之间较容易区分。图1
新窗口打开|下载原图ZIP|生成PPT图1不同氮素水平下高光谱反射率与红边
N0: 施氮0 kg hm-2; N75: 施氮75 kg hm-2; N150: 施氮150 kg hm-2; N225: 施氮225 kg hm-2。
Fig. 1Canopy hyperspectral reflectance and red edge at different nitrogen levels
N0: nitrogen application of 0 kg hm-2; N75: nitrogen application of 75 kg hm-2; N150: nitrogen application of 150 kg hm-2; N225: nitrogen application of 225 kg hm-2.
由图1可知, 红边与原始冠层光谱图曲线近红外(780~1300 nm)波段基本一致, 各施氮处理之间有明显差异, 且具有“双峰”现象, 说明不同施氮处理对甜菜植株有较大影响, 不同施氮处理除N75红边位置与N0相同, 其他施氮处理均出现“红移”现象, 施氮处理N150、N225较N0均向长波方向移动5 nm, 说明随着施氮量的增加, 甜菜植株冠层单位光合效率有所提高, 但在施氮处理N150后红边位置趋于平稳, 这与原始光谱反射率可见光波段“绿峰”处保持一致。
2.2 原始光谱反射率与全氮含量的相关性
结果如图2所示, 冠层光谱反射率与冠层叶片相关性稳定波段主要集中350~1300 nm与1418~ 1809 nm, 由于在1300~2500 nm叶面光谱反射率主要受叶内水含量的控制, 叶内水含量降低, 其光谱反射率升高, 且氮素含量的敏感波段主要集中在短波红外与可见光波段, 所以选择350~1300 nm作为全氮含量估测波段。图2
新窗口打开|下载原图ZIP|生成PPT图2甜菜冠层原始光谱反射率与全氮含量相关性
F0.05: 显著相关(P < 0.05); F0.01: 极显著相关(P < 0.01)。
Fig. 2Correlation between original spectral reflectance and total nitrogen content in sugar beet canopy
F0.05: significant correlation (P < 0.05); F0.01: extremely significant correlation (P < 0.01).
由图2可知, 在350~1300 nm范围内, 在可见光波段与冠层叶片全氮含量呈现极显著正相关的波段主要集中在726~780 nm, 且在波段767 nm处达到最大正相关(r = 0.609, P<0.01), 与全氮含量呈现极显著负相关的波段主要集中在591~615 nm, 且在波段604 nm达到最大负相关(r = -0.267, P<0.01); 在近红外波段与冠层叶片全氮含量呈现极显著正相关的波段主要集中在780~1300 nm, 且在波段1104 nm处达到最大正相关(r = 0.68, P<0.01), 但在此波段范围内没有出现与全氮含量的极显著负相关。
2.3 一阶微分光谱与全氮含量相关性
叶片全氮含量的敏感波段主要集中在短波红外与可见光波段, 因此选择350~1300 nm一阶导数光谱反射率与叶片全氮含量作相关性分析。由图3可知, 在波段389~403、610~620、638~679、709~766、774~805、995~1073、1200~1252 nm范围内呈现极显著正相关, 且在可见光(380~780 nm)范围内在波长762 nm处达到最大正相关(r = 0.703, F<0.01), 在近红外(780~1300 nm)范围内在波长1018 nm处达到最大正相关(r = 0.633, F<0.01); 在波段445~457、486~517、683~693、943~964、1126~1176、1285~1300 nm范围内呈现极显著负相关关系, 且在可见光(380~780 nm)范围内在波长496 nm处达到最大负相关(r = -0.545, F<0.01), 在近红外波段(780~1300 nm)范围内在波长1138 nm处达到最大负相关(r = -0.715, F<0.01); 在红边(680~760 nm)范围内, 最大正相关位于波长747 nm处(r = 0.646, F<0.01), 最大负相关位于波长687 nm处(r = -0.332, F<0.01)。图3
新窗口打开|下载原图ZIP|生成PPT图3甜菜冠层一阶导数光谱与全氮含量相关性
F0.05: 显著相关(P < 0.05); F0.01: 极显著相关(P < 0.01)。
Fig. 3Correlation between the first derivative spectra of sugar beet canopy and the total nitrogen content
F0.05: significant correlation (P < 0.05); F0.01: extremely significant correlation (P < 0.01).
2.4 光谱特征参数、敏感波段植被指数与全氮含量相关性
全氮含量与一阶导数光谱的相关性, 在相应波段内选取相关系数最大的对应波段反射率, 构建原始光谱反射率的比值、差值与归一化比值植被指数, 构建一阶微分光谱反射率的比值、差值、归一化比值与红边比值、差值与归一化差异指数。研究并分析三边参数及敏感波段植被指数与全氮含量的相关性(表2)。Table 2
表2
表2光谱特征参数及敏感波段植被指数与全氮含量相关性
Table 2
光谱特征参数 Spectral characteristic parameter | 相关关系 Correlation | 植被指数 Vegetation index | 相关性Correlation |
---|---|---|---|
Dr | 0.590** | 原始光谱比值植被指数Ratio vegetation index R1104/R767 | 0.358** |
λr | 0.497** | 原始光谱差值植被指数Difference vegetation index R1104-R767 | 0.280** |
SDr | 0.601** | 原始光谱归一化植被指数Normalized ratio vegetation index (R1104-R767)/(R1104+R767) | 0.364** |
Dy | 0.074 | 原始光谱比值植被指数Ratio vegetation index R767/R604 | 0.652** |
Λy | 0.089 | 原始光谱差值植被指数Difference vegetation index R767-R604 | 0.433** |
SDy | -0.086 | 一阶微分比值植被指数Ratio vegetation index Dr1138/Dr762 | -0.169 |
Db | 0.072 | 一阶微分差值植被指数Difference vegetation index Dr1138-Dr762 | -0.744** |
λb | 0.506** | 一阶微分归一化植被指Normalized ratio vegetation index (Dr1138-Dr762)/(Dr1138+Dr762) | 0.500** |
SDb | -0.216 | 一阶微分比值植被指数Ratio vegetation index Dr762/Dr496 | 0.586** |
R604 | -0.267** | 一阶微分差值植被指数Difference vegetation index Dr762-Dr496 | 0.721** |
R1104 | 0.680** | 红边比值植被指数Red edge ratio vegetation index Dr747/Dr687 | 0.431** |
Dr1138 | -0.715** | 红边差值植被指数Red edge difference vegetation index Dr747-Dr687 | 0.589** |
Dr762 | 0.703** | 红边归一化差异指数Red edge normalized difference index (Dr747-Dr687)/( Dr747+Dr687) | 0.545** |
新窗口打开|下载CSV
由表2可知, 所选光谱特征参数中红边参数与叶片全氮含量均达到极显著相关, 黄边参数与蓝边参数中只有蓝边位置与叶片全氮含量达到极显著相关, 在三边参数中, 红边振幅、红边面积、蓝边位置与叶片全氮含量均达到0.5以上的极显著相关, 其中红边面积与之达到最大相关(r = 0.601, F < 0.01); 原始光谱反射率和一阶微分光谱反射率与叶片全氮含量的最大相关性均达到极显著水平, 其中一阶微分最大负相关Dr1138与一阶微分最大正相关Dr762的相关性均达到0.7以上; 所建植被指数与全氮含量的相关性除一阶微分比值植被指数Dr1138/Dr762没有达到极显著水平, 其他植被指数均达到极显著水平, 但整体与叶片全氮含量的相关性较小, 在原始光谱植被指数中, 只有原始比值植被指数R767/ R604与叶片全氮含量的相关性达到了0.6以上; 在一阶微分植被指数与叶片全氮含量的相关性中, 除一阶微分比值植被指数Dr1138/Dr762外均达到0.5以上极显著水平, 其中一阶微分差值植被指数Dr1138-Dr762与一阶微分差值植被指数Dr762- Dr496的相关系数均达到0.7以上; 红边植被指数与叶片全氮含量均达到极显著相关水平, 其中红边差值植被指数Dr747-Dr687与红边归一化差异指数的均达到0.5以上, 且一阶微分光谱植被指数和红边光谱植被指数与叶片全氮含量的相关整体大于原始光谱植被指数的相关性。
2.5 已知高光谱植被指数与全氮含量相关性
在前人研究的基础上, 针对本试验数据对已有植被指数优选, 分析优选高光谱植被指数与全氮含量的相关性, 由表3可知, 所选的15个高光谱植被指数均与全氮含量极显著相关, 其中达到0.5以上相关性的有8个, 植被指数SAVI、OSAVI和MSAVI与全氮含量的相关性均达到0.6以上, 且土壤调节植被指数SAVI与全氮含量达到最大相关性(r = 0.692**, F<0.01)。Table 3
表3
表3已知高光谱植被指数与全氮含量的相关性
Table 3
植被指数 Vegetation index | 计算公式 Calculation formula | 相关关系Correlation | 文献来源 Source of literature |
---|---|---|---|
RSI | R990/R720 | 0.528** | [23] |
NDSI | (R860-R720)/(R860+R720) | 0.582** | [23] |
NDVI | (R790-R670)/(R790+R670) | 0.332** | [24] |
FD-NDVI | (R730-R525)/(R730+R525) | 0.451** | [25] |
RENDVI | (R750-R705/( R750+ R705) | 0.367** | [26] |
mND705 | (R750-R705)/(R750+2R445) | 0.360** | [27] |
GNDVI | (R790-R550)/(R790+R550) | 0.464** | [28] |
SAVI | 1.5(R870-R680)/(R870+R680+0.5) | 0.692** | [29] |
OSAVI | (1+0.16)(R810-R680)/(R810-R680+0.16) | 0.633** | [30] |
MSAVI | 0.5{2R800+1-[(2R800+1)2-8(R800-R670)]0.5} | 0.623** | [31] |
DCNI | (R720-R700/(R700-R670)/(R720-R670+0.03) | 0.460** | [32] |
CIgreen | [(R840-R870)/R550]-1 | -0.493** | [33] |
NINI | [lg(1/R1510)-lg(1/R1680)]/[lg(1/R1510)+lg(1/R1680)] | 0.553** | [34] |
TVI | 0.5[120(R750-R550)-200(R670-R550)] | 0.538** | [35] |
DSI | R800-R680 | 0.510** | [36] |
新窗口打开|下载CSV
2.6 基于高光谱特征参数与敏感植被指数的全氮含量的估算与验证
选取表2、表3相关系数绝对值大于0.65的植被指数用于甜菜叶片全氮含量的估算。其中光谱特征参数与植被指数作为自变量(x), 叶片全氮含量作为因变量(y), 构建叶片全氮含量的一元线性、二次项、幂函数以及指数函数估算模型。通过对比各模型的建模决定系数(R2)、相对误差(RE)以及均方根误差(RMSE)筛选出最优模型用以研究新疆滴灌条件下甜菜叶片全氮含量的估算。其中R2用以判断模型拟合优度, RE与RMSE用于检验估测模型的可靠程度。2.6.1 基于高光谱特征参数的全氮含量的估算与验证 通过对构建光谱特征参数模型与精度检验分析可知(表4), 所选光谱特征参数与全氮含量均具有较好的拟合关系, 同时具有较好的拟合精度(R2 = 0.49~0.78)。通过比较不同函数估算模型发现, 不同光谱特征参数所建模型均表现为二次项函数模型最优, 其中光谱特征参数Dr1138与Dr762的决定系数均达到0.6以上, 说明在光谱短波近红外波段与红光处参数有较好估算全氮含量的能力。通过比较决定系数发现, 光谱特征参数Dr762在所选参数中表现最优, 不同函数构建估算模型决定系数均达到0.7以上, 且以二次项函数所建模型决定系数最高(R2 = 0.772), 其次为线性函数(R2 = 0.770)。
Table 4
表4
表4不同光谱特征参数与全氮含量回归关系模型及验证
Table 4
高光谱特征参数 Characteristic parameter | 拟合模型 Fitting model | 验证模型 Performance model | |||
---|---|---|---|---|---|
拟合方程 Equation | R2 | RE(%) | RMSE | ||
R1104 | y = 63.88x-8.167 | 0.516 | 21.429 | 4.762 | |
y = -30.27x2+96.39-16.64 | 0.519 | 22.299 | 4.829 | ||
y = 58.24x1.296 | 0.513 | 20.898 | 4.766 | ||
y = 7.555e2.277x | 0.493 | 22.904 | 4.898 | ||
Dr1138 | y = 64.73x-1.144 | 0.643 | 19.893 | 4.874 | |
y = -47.08x2+104.7x-9.162 | 0.648 | 20.625 | 4.945 | ||
y = 63.68x1.031 | 0.643 | 19.909 | 4.881 | ||
y = 9.629e2.313x | 0.615 | 23.265 | 5.010 | ||
Dr762 | y = 157.8x-147 | 0.770 | 20.877 | 5.038 | |
y = -192.3x2+578.8x-377.1 | 0.772 | 21.518 | 5.137 | ||
y = 13.96x6.461 | 0.747 | 21.635 | 4.914 | ||
y = 0.04302e5.813x | 0.740 | 22.103 | 4.928 |
新窗口打开|下载CSV
利用检验样本对上述变量建立的估算模型预测精度进行检验, 各模型检验结果如表4所示, 通过对比光谱特征参数Dr762的RE(%)与RMSE发现,虽然二次项函数估算模型决定系数具有更好优势, 但RE(%)与RMSE较高, 且幂函数与指数函数在验证模型中均表现出较低的RE(%)与RMSE, 难以判定出最优估算模型, 因此对估测值与真实值拟合模型再次验证, 通过对光谱特征参数Dr762下不同函数估测值与真实值拟合分析发现(图4), 幂函数估算模型较其他估算模型具有更高的决定系数, 因此光谱特征参数Dr762下幂函数估算模型具有更好的拟合优度。
图4
新窗口打开|下载原图ZIP|生成PPT图4甜菜光谱特征参数Dr762 (一阶微分最大正相关)下全氮含量估测模型验证
Fig. 4Verification of the total nitrogen content estimation model under the spectral characteristic parameters Dr762 (maximum positive correlation of first order differential) of sugar beet
2.6.2 基于高光谱植被指数的叶片全氮含量估算与验证 通过对构建光谱植被指数与精度检验分析可知(表5), 所选光谱植被指数与全氮含量均具有较好的拟合关系, 同时具有较好拟合精度(R2 = 0.415~0.804)。通过比较不同函数估算模型决定系数发现, 不同植被指数所建模型均表现为二次项函数估算模型最优, 其中差值植被指数Dr1138-Dr762与Dr762-Dr496决定系数均达到0.7以上, 且在差值植被指数Dr762-Dr496下二次项估算模型的决定系数达到最高值(R2 = 0.804), 这说明在短波近红外波段与红光波段组合和红光波段与绿光波段组合具有较好估算叶片全氮含量的能力。
Table 5
表5
表5不同植被指数与全氮含量回归关系模型及验证
Table 5
高光谱特征参数 Characteristic parameter | 拟合模型 Fitting model | 验证模型 Performance model | |||
---|---|---|---|---|---|
拟合方程 Equation | R2 | RE(%) | RMSE | ||
SAVI | y = 114.7x-105.7 | 0.541 | 24.492 | 5.092 | |
y = 104.5x2-122.1x+28.18 | 0.544 | 23.153 | 4.926 | ||
y = 12.39x5.208 | 0.539 | 21.878 | 4.830 | ||
y = 0.1489e4.474x | 0.535 | 21.514 | 4.829 | ||
Dr1138-Dr762 | y = 49.36x+0.4913 | 0.732 | 19.708 | 4.709 | |
y = -1.864x2+51.27x+0.04145 | 0.732 | 19.777 | 4.721 | ||
y = 49.73x0.9814 | 0.732 | 19.746 | 4.715 | ||
y = 9.737e1.842x | 0.714 | 20.934 | 4.656 | ||
Dr762-Dr496 | y = 12.85x+5.673 | 0.794 | 23.008 | 5.372 | |
y = -2.828x2+21.68x-0.3212 | 0.804 | 22.996 | 5.608 | ||
y = 18.64x0.755 | 0.799 | 23.025 | 5.447 | ||
y = 11.91e0.4736x | 0.759 | 24.623 | 5.288 |
新窗口打开|下载CSV
利用检验样本对上述变量建立的估算模型预测精度进行检验, 各模型检验结果如表5所示, 通过对比差值植被指数Dr762-Dr496的RE(%)与RMSE发现, 差值植被指数Dr762-Dr496估算模型与光谱特征参数Dr762表现出相似特征, 虽然二次项函数估算模型决定系数最高, 表现出较高拟合优度, 但RMSE在整体函数估算模型中表现最高, 说明偏离真值的程度较大, 虽然差值植被指数Dr1138-Dr762较Dr762-Dr496的决定系数有所降低, 但其表现出较低的RE(%)与较低RMSE, 因此难以对最优全氮估算模型进行判定, 因此对估测值与真实值拟合模型再次验证, 通过对差值植被指数Dr1138-Dr762与Dr762-Dr496下不同函数估测值与真实值拟合分析发现(图5和图6), 在差值植被指数Dr1138-Dr762中指数函数估测值与真实值拟合程度表现最优(R2 = 0.739), 但与差值植被指数Dr762-Dr496的拟合程度相比表现较低。在差值植被指数Dr762-Dr496下线性函数估测值与真实值拟合决定系数和指数函数估测值与真实值拟合决定系数相同(R2 = 0.759), 但线性函数估算模型较指数函数估算模型具有较高决定系数与较低RE(%), 因此差值植被指数Dr762-Dr496下线性函数估算模型具有更好的拟合优度。
图5
新窗口打开|下载原图ZIP|生成PPT图5甜菜植被指数Dr1138-Dr762下全氮含量估测模型验证
Fig. 5Validation of the estimation model of total nitrogen content under vegetation index Dr1138-Dr762
图6
新窗口打开|下载原图ZIP|生成PPT图6甜菜植被指数Dr762-Dr496下全氮含量估测模型验证
Fig. 6Validation of the estimation model of total nitrogen content under vegetation index Dr762-Dr496
3 讨论
3.1 不同氮素水平对甜菜全氮含量及冠层光谱反射率影响
在不同氮素水平下甜菜原始冠层光谱反射率波形相似, 在350~670 nm波段范围内有明显的反射峰与吸收谷, 光谱呈现先增高后降低的趋势, 在波段670~760 nm波形变化较为剧烈, 波形呈现为陡然上升趋势, 在近红外波段780~1300 nm处形成近红外反射平台, 且具有明显波形变化, 这与前人研究基本一致[37,38]。施氮处理N75、N150、N225在550 nm较CK分别降低了7.5%、28.0%、23.0%, 说明各施氮处理相对CK均可提高甜菜植株冠层叶片对红、蓝光的吸收, 提高植株的光合效率, 且施氮处理N150比N75与N225具有更好的提高植株光合效率的作用; 在一阶导数光谱图中不同施氮处理除N75红边位置与N0相同, 其他施氮处理均出现“红移”现象, 施氮处理N150、N225较N0均向长波方向移动5 nm, 说明随着施氮量的增加, 甜菜植株冠层单位光合效率有所提高。3.2 光谱特征参数及植被指数的选择与全氮含量的构建与检验
针对建立光谱特征参数及敏感植被指数进行筛选, 其中“三边”参数中红边面积具有较好估算全氮含量的能力, 在植被指数中R767/R604、Dr1138- Dr762与Dr762-Dr496的相关性均在0.6以上, 其中一阶微分光谱植被指数和红边光谱植被指数与叶片全氮含量的相关性整体大于原始光谱植被指数的相关性, 组合波段相对单一光谱特征参数相关性均有所提高。通过对光谱特征参数筛选建立全氮含量的多种函数估算模型发现, 不同光谱特征参数所建立模型均表现为二次函数拟合效果最好, 利用检验样本对所建立估测模型预测精度进行检验发现, 虽然光谱特征参数Dr762二次函数估测模型决定系数具有最大决定系数(R2 = 0.772),但验证模型中RE(%)与RMSE较高, 通过对Dr762估测值与真实值拟合后发现, 光谱特征参数Dr762幂函数下估算模型具有更好的拟合优度。
通过对敏感植被指数筛选建立全氮含量的多种函数估算模型发现, 不同植被指数所建模型均表现为二次函数估算模型最优, 其中差值植被指数Dr762-Dr496下二次项估算模型的决定系数达到最高值(R2 = 0.804), 这说明在短波近红外波段与红光波段组合和红光波段与绿光波段组合具有较好估算叶片全氮含量的能力。利用检验样本对所建模型进行检验发现, 虽然差值植被指数Dr762-Dr496表现出较高拟合优度, 但验证模型中偏离真值的程度较大, 通过建立Dr762-Dr496估测值与真实值拟合后发现, 差值植被指数Dr762-Dr496下一元线性函数具有更好的拟合优度。
作物冠层叶片氮素含量的高光谱遥感估算较为复杂, 受到作物冠层结构、叶片表面以及内部特征、化学和光学特性、土壤背景值等一系列影响。Knyazikhin等[39]研究认为冠层结构及其光学特性是决定光谱反射率的主要因素, 冠层结构组织具有复杂性, 光子在此结构中会产生碰撞、吸收、散射及反射, 只有一小部分会被截获, 因而叶片反射率与叶片氮素含量间无直接相关性。但叶片、冠层以及生态系统之间的复杂联系可能会导致叶片全氮含量与光谱反射率之间的重复相关性, 虽然可以确定冠层光谱反射率的物理驱动因素是必要的, 模拟冠层中作物生化指标可以起到较好作用, 然而这些物理机制在冠层、叶片内部结构及外部结构等和生态系统之间的一系列生态驱动是联动运作的, 而氮素是作物碳氮循环中一最为重要成分, 氮素有效性会较大程度决定植被冠层丰度[40,41]。碳氮循环作为生态系统中运作的关键, 光合作用为植被提供生长及繁殖的能量及物质, 其冠层氮素含量的定量变化提供了关于植被冠层的直接信息[42]。这也说明氮素含量与冠层结构之间的联动, 并较大程度决定冠层结构, 并与冠层光谱反射率呈现直接或间接的影响。因此对作物的生物学和物理学特性不能孤立地评估, 应更好理解光谱处理方法以及工作原理, 在更多波长范围内进行更广泛的模拟。
本文光谱数据测定主要在甜菜叶丛快速生长期至糖分积累期, 在此阶段范围内甜菜冠层基本达到土壤背景全覆盖, 通过对比表明前人建立植被指数与本文全氮含量均具有较好相关性, 均达到极显著相关, 说明不同纬度地区、不同作物与所建立植被指数均可进行全氮含量的估算, 但由于地区与作物的不同, 所建立植被指数与叶片全氮含量的吻合度存在一定的差异。本研究主要基于石河子滴灌甜菜叶片全氮含量估算得到的初步结果, 但在实际工作中利用高光谱遥感数据对甜菜叶片全氮含量估算时, 易受不同纬度地区以及冠层结构、叶片角度、结构等因素影响, 因此所得结果还有待进一步验证。
4 结论
通过光谱特征参数与叶片全氮含量建立多种函数估算模型, 其中光谱特征参数Dr762幂函数下估算模型具有较好估算甜菜叶片全氮含量的能力, 其决定系数R2 = 0.747, 验证相对误差RE(%)为21.635, 验证均方根误差RMSE为4.914。通过植被指数与叶片全氮含量建立多种函数估测模型, 其中差值植被指数Dr762-Dr496下一元线性函数具有较好估算甜菜叶片全氮含量的能力, 其决定系数R2 = 0.794, 验证相对误差RE(%)为23.008, 验证均方根误差为5.372。参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子
,
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]