一种不同植株氮含量水平下小麦植株含水率的监测方法

技术领域

本发明属于作物生长监测领域，涉及一种不同植株氮含量水平下小麦植株含水率的监测 方法，具体涉及一种基于高光谱技术确定不同植株氮含量水平下的小麦植株含水率核心波段 共性区域最优光谱指数、建立基于最优光谱指数的小麦PWC监测模型的方法，特别适用于 不同氮营养水平下小麦植株含水率的高光谱监测研究。

背景技术

基于高光谱遥感技术可实现作物水分状况的实时监测和快速诊断，对提高作物灌溉管理 水平和水分利用效率具有重要意义。

植株含水率（PWC）能全面反应作物需水信息。近年来，国内外学者对作物PWC的敏 感波段和特征光谱参数做了大量研究，提出了形式多样的水分光谱指数，但这些光谱指数往 往只考虑了的单一的水分影响。已有研究证明，作物体内生化成分间存在强烈的相关性，其 中氮（N）素通过叶绿素、木质素和纤维素等与作物水分间接相关，因此，不同植株氮含量 （PNC）水平必然影响作物生长发育，从而间接影响作物PWC监测。由于实际大田生产中 往往涉及不同施氮水平，而已有水分光谱指数鲜有考虑不同氮素营养水平的影响，且部分光 谱指数虽对整体PWC监测研究有较高精度，但在某些PNC水平下对PWC的监测精度则较 低，误差较大，因此，迫切需要一种确定适用于不同PNC水平条件下的PWC光谱指数的技 术和方法。

发明内容

本发明的目的是针对上述现有技术的不足，充分考虑不同氮营养水平对小麦水分监测的 影响，通过将所有PWC及对应冠层光谱反射率按照不同PNC大小水平分类，寻找不同 PNC水平下的PWC核心波段共性区域，确定适用于不同PNC水平条件下的PWC最优光谱 指数，从而建立基于最优光谱指数的小麦PWC监测模型，该方法可以快速、无损、精准的 估测不同水氮条件下的小麦PWC。

本发明的目的是通过以下技术方案实现的：

一种不同植株氮含量水平下小麦植株含水率的监测方法，包括以下步骤：

步骤一、取样，采集小麦冠层光谱反射率、PNC和PWC数值；样本点采自不同生育 期、不同水分处理、不同施氮水平和不同年份；

步骤二、按不同PNC水平将全部PWC和冠层光谱反射率分为三个子数据集；

步骤三、获取不同PNC水平子数据集下PWC与NDSI的决定系数（R²）结果集；

步骤四、确定不同PNC水平下的PWC核心波段共性区域；

步骤五、确定基于核心波段共性区域的最优光谱指数；

步骤六、建立基于最优光谱指数的小麦PWC监测模型；

步骤七、检验小麦PWC监测模型的准确性和普适性。

步骤一中，设置不同水分处理和不同施氮水平互作的小麦田间实验，采集拔节至灌浆各 关键生育时期冠层光谱反射率，并在每小区内同步选取可表征小区平均长势的10个单茎进 行破坏性取样，称量鲜重和干重，计算PWC结果，磨样称重后采用凯氏定氮法测定对应 PNC。

采集的冠层高光谱反射率分辨率高、信息量大，光谱波段范围为350-2500nm，波段间 隔为1nm。

PWC的计算方式为：植株含水率(PWC)(%)=(W_F-W_D)/W_F×100，其中W_F为植株鲜重 （g）总和，W_D为植株干重（g）总和。

步骤二中，综合考虑PNC大小范围，及各范围内数据量分布情况，将PNC分为 PNC<1.3%、PNC1.3-1.6%、PNC>1.6%三个范围，按照对应PNC值范围，将全部的PWC和 对应的冠层光谱反射率分为对应不同PNC水平的三个子数据集。

基本分类原则为：保证各水平PNC值的范围存在一定差异的基础上，控制各PNC水平 所属子数据集的数据量无太大差异。总数据量越大越有助于寻找核心波段共性区域。

步骤三中，使用Matlab9.0编程计算，获取不同PNC水平子数据集下PWC与350-2500nm 波段范围内任意两波段组合的NDSI所建对应模型的决定系数（R²）结果集；并将各结果集 中所有波段组合的模型R²矩阵数据绘制成等高线图，以颜色深度变化显示R²大小的变化。

R²结果集具体为：在350-2500nm波段内，以1nm为间隔的两两波段对应组合的NDSI 与PWC所建模型的R²集，呈矩阵排列。

步骤四中，以上述各决定系数（R²）结果集中R²的最大值为标准，定义R²前10%对应的 等高线图区域为核心波段区域；通过Matlab9.0寻找各PNC水平子数据集的核心波段区域的 相同部分，即为核心波段共性区域。

步骤五中，在步骤四已确定的核心波段共性区域内计算两两波段组合的所有NDSI与 PWC所建模型所对应的建模精度（R²）和检验误差（RRMSE），选取R²-RRMSE最大值对应 的光谱指数为最优光谱指数，确定最优光谱指数为NDSI(1302,1190)。

步骤六中，基于最优光谱指数NDSI(1302,1190)的小麦PWC监测模型为：Y=- 1289.4X+78.29，其中，模型R²为0.876，SE为2.702。

步骤七中，利用独立年份小麦试验数据验证监测模型的准确性和普适性，采用复相关系 数R²、相对根均方差RRMSE对监测模型进行综合评价（检验监测模型R²为0.7979，RRMSE 为0.0662），并验证监测模型对PNC的敏感性（PNC检验模型R²为0.3242，RRMSE为 45.7011）。

本发明相对传统技术的优点：

本发明基于不同PNC分类水平下PWC核心波段共性区域筛选出的NDSI最优波段组 合，充分考虑了低、中、高不同氮素营养条件对PWC监测的影响，所选光谱指数适用于不 同氮素营养下的小麦PWC高光谱监测。而传统光谱指数构建筛选技术则往往只是针对PWC 的研究，很少考虑不同氮营养条件的影响，特别是对不同氮营养水平条件下的PWC监测广 适性鲜有提到。本发明可以快速、无损、精准的估测不同水氮条件下的小麦PWC，将对精 确农业中不同氮素营养条件下的小麦植株水分含量的高光谱监测提供重要技术支撑。

附图说明

图1为不同PNC水平下的小麦PWC与任意两波段组合的NDSI决定系数前10%等高线 图（A:PNC<1.3%;B:PNC1.3-1.6%;C:PNC>1.6%）及其共性区域图（D）。

图2为基于不同PNC水平核心波段共性区域NDSI(1302,1190)的小麦PWC模型构建 （A）和模型检验（B）。

图3为基于总数据集的小麦PWC与任意两波段组合的NDSI决定系数等高线图。

图4为基于总数据集NDSI(1727,1539)的小麦PWC模型构建（A）和模型检验（B）。

图5为本发明不同植株氮含量水平下小麦植株含水率的监测方法流程图。

具体实施方式

本发明通过实施不同年份、不同施氮水平和不同水分处理互作的两个小麦试验，结合附 图，考察了基于高光谱技术的不同PNC水平下的PWC核心波段共性区域最优光谱指数的应 用优势。

如图5所示，不同植株氮含量水平下小麦植株含水率的监测方法，包括以下步骤：

S101：取样，采集小麦冠层光谱反射率、PNC和PWC数值；样本点采自不同生育期、 不同水分处理、不同施氮水平和不同年份；

S102：按不同PNC水平将PWC和冠层光谱反射率数据分为三个子数据集；

S103：获取不同PNC水平子数据集下PWC与NDSI的决定系数（R²）结果集；

S104：确定不同PNC水平下的PWC核心波段共性区域；

S105：确定核心波段共性区域的最优光谱指数；

S106：建立基于最优光谱指数的小麦PWC监测模型；

S107：检验小麦PWC监测模型的准确性和普适性。

本发明利用不同年份（连续两季2010.11-2011.06、2011.11-2012.06）、不同水分处理、 不同施氮水平和不同生育期的小麦田间试验，采集不同氮含量水平下的小麦植株含水率及对 应冠层光谱反射率，按照不同PNC大小水平分类，寻找不同PNC水平下的PWC核心波段 共性区域，确定适用于不同PNC水平条件下的PWC最优光谱指数，从而基于最优光谱指数 构建面向小麦拔节到灌浆期的监测模型。

小麦品种为扬麦18，种植方式采用条播，行距为25cm，基本苗为每亩12万株，小区 面积为10m²。

试验1（2010.11-2011.06）中设置4个水分处理水平，为9.5-10.5%、15.5-16.5%、21.5- 22.5%和29.5-30.5%（以体积含水率表示）；设置2个施氮水平，为150kg/hm²、 300kg/hm²。试验2（2011.11-2012.06）中设置3个水分处理为13.5-14.5%、21.5-22.5%、 29.5-30.5%；设置3个施氮水平，为90kg/hm²、180kg/hm²、270kg/hm²。试验1用于建立监 测模型，试验2用于检验监测模型。

具体包括以下步骤：

S101：取样，采集小麦冠层光谱反射率、PNC和PWC数值。实施不同水分处理和不同 施氮水平互作的小麦田间实验，设置4个水分处理水平，为9.5-10.5%、15.5-16.5%、21.5- 22.5%和29.5-30.5%（以体积含水率表示），2个施氮水平，为150kg/hm²、300kg/hm²。

采集2010.11-2011.06扬麦18从拔节至灌浆各关键生育时期的冠层光谱反射率，冠层光 谱反射率采用野外高光谱辐射仪，波段范围为350-2500nm，波段间隔为1nm；冠层光谱反 射率测定：选择晴朗无云的天气，10:00-14:00，距冠层上方1m测试。

于各小区同步选取可表征本小区平均长势的10株单茎进行破坏性取样，称量鲜重和干 重，计算PWC，磨样称重后采用凯氏定氮法测试计算对应PNC。

PWC计算方式：植株含水率(PWC)(%)=(W_F-W_D)/W_F×100，其中W_F为植株鲜重（g） 总和，W_D为植株干重（g）总和。

S102：按不同PNC水平将PWC和冠层光谱反射率数据分为三个子数据集。将采集的 PWC及对应冠层光谱反射率的总数据集按照PNC值大小范围及数据量分布情况，分为 PNC<1.3%、PNC1.3-1.6%、PNC>1.6%三个范围，按照对应PNC值范围，将全部的PWC和 对应的冠层光谱反射率分为对应不同PNC水平的三个子数据集。

基本分类原则：保证各水平PNC值的范围存在一定差异的基础上，控制各PNC水平所 属子数据集的数据量无过大差异。其中，试验1中PNC<1.3%、PNC1.3-1.6%、PNC>1.6%对 应的三个子数据集数据量分别为：32、27、21，试验2中子数据集数据量分别为：23、38、 46。

S103：获取不同PNC水平子数据集下PWC与NDSI的决定系数（R²）结果集。使用 Matlab9.0编程计算，获取不同PNC水平子数据集下PWC与350-2500nm波段范围内任意两波 段组合的NDSI所建对应模型的决定系数（R²）结果集；并将各结果集中所有波段组合的模 型R²矩阵数据绘制成等高线图输出，以颜色深度变化显示R²大小的变化。

R²结果集具体为：在350-2500nm波段内，以1nm为间隔的两两波段对应组合的NDSI 与PWC所建模型的R²集，呈矩阵排列。

S104：确定不同PNC水平下的PWC核心波段共性区域。以上述各决定系数（R²）结果 集中R²的最大值为标准，定义R²前10%对应的等高线图区域为核心波段区域（如图1中A、 B、C）；通过Matlab9.0寻找各PNC水平子数据集的核心波段区域的相同部分，即为核心波段 共性区域，以等高线图形式输出（如图1中D）。

S105：确定基于核心波段共性区域的最优光谱指数。在上述确定的核心波段共性区域内 综合计算两两波段组合的所有NDSI与PWC所建模型所对应的建模精度（R²）和检验误差 （相对根均方差RRMSE），综合考虑R²和RRMSE大小，选取R²-RRMSE最大值对应的核心波 段共性区域波段组合NDSI(1302,1190)为最优光谱指数，其中1190nm位于水分特征吸收波 段1200nm附近，1302nm为PWC的最优敏感波段1300nm附近。

S106：如图2A，建立基于最优光谱指数NDSI(1302,1190)的小麦PWC监测模型：Y=- 1289.4X+78.29，其中，检测模型R²为0.876，SE为2.702。

S107：如图2B，利用独立年份（试验2，2011.11-2012.06）小麦试验数据验证监测模型 的准确性和普适性，采用复相关系数R²、相对根均方差RRMSE对监测模型进行综合评价 （检验监测模型R²为0.7979，RRMSE为0.0662），并验证监测模型对PNC的敏感性（PNC检 验模型R²为0.3242，RRMSE为45.7011）。

其中，RRMSE的计算公式如下：

n为样本数，P_i为模型预测值，O_i为试验观察值，为观察平均值。

若按照传统的光谱指数确定方式，即不分不同PNC水平的总数据集下计算PWC最优光谱 指数，则结果如图3所示：PWC的最优光谱指数波段组合为NDSI（1727,1539），其中核心 波段均位于短波红外范围内。该光谱指数对总数据集有较高的建模和检验精度，但对不同 PNC水平的监测和预测性能不稳定，如：PNC>1.6时模型的监测精度很低（R²=0.5064），且 存在较明显的饱和现象（见图4A、B虚线圆圈内区域）；PNC1.3-1.6水平下的预测能力也偏 低（R²=0.6702）（见图4B）。

本发明中基于不同PNC水平核心波段共性区域最优光谱指数NDSI(1302,1190)所建立 的小麦PWC监测模型，对不同PNC水平下的PWC均显示出较高的监测预测能力，同时对 PNC不敏感，且相对传统光谱指数构建的监测模型而言极大程度地克服了高PNC水平下 PWC模型的易饱和现象。进一步与前人提出的光谱指数比较发现，NDSI(1302,1190)的总监 测预测能力远优于前人光谱指数。说明本发明基于高光谱技术确定的不同PNC水平下的小 麦植PWC核心波段共性区域最优光谱指数NDSI(1302,1190)具有较好的应用前景。该方法 克服了传统方法对试验因素考虑不全面的缺陷，同样适用于其他多因素互作的农作物生产监 测研究。

表1基于不同光谱指数的小麦植株含水率建模（n=271）及检验（n=308）效果