本发明涉及遥感技术应用领域，具体涉及基于高光谱的水稻穗期氮营养监测与诊断方法及应用。
背景技术：
随着城镇化建设不断推进，从事农业劳动人口显著降低，稻作经营主体和稻作方式均较传统水稻生产发生极大变化，以土地流转形式为主的规模化经营主体以及机械化轻简栽培模式呈逐年增长趋势。
氮素营养作为作物生长发育和产量、品质的重要影响因子，是植物需求量最大的矿质营养元素。氮素营养的亏缺对水稻的生长尤为重要，基于高光谱的氮素营养诊断手段就是一种表观诊断和化学诊断相结合的现代监测技术，其通过光谱数据快速无损的了解作物氮素营养状况。随着遥感技术的应用推广，作物氮素的遥感监测在作物遥感领域研究已经得到一定的发展，利用光谱影像、光谱反射率及其衍生出的植被指数快速获取农作物的农艺指标的研究实例已有很多。薛利红等人发现群体含氮量与光谱指数显著相关，杨燕穷等人应用3ccd多光谱成像技术估算生物量等农艺指标。综上所述，基于光谱技术的研究大多为水稻氮素诊断某一环节提供理论依据，未针对水稻生产实践提供一系列诊断方案，尤其是在特定时期对水稻给出准确的氮素施肥方案。因此，提出水稻穗期氮素亏缺诊断方法及应用。
技术实现要素：
本发明所要解决的技术问题在于：如何解决水稻生长过程中对产量影响较大的特定时期中存在的氮肥亏缺问题，提供了水稻穗期氮素亏缺诊断方法，构建了基于归一化植被指数(ndsi(1136，1587))的水稻氮浓度模型和基于土壤调节优化植被指数(osavi)的水稻地上部干物质量模型，并利用由水稻地上部干物质量建立的临界氮浓度模型对水稻穗期氮素营养状况进行估算，并结合各生态点穗期氮肥利用率，给出施肥方案。
本发明是通过以下技术方案解决上述技术问题的，本发明包括以下步骤：
s1：测量水稻冠层光谱、地上部干物质量和氮浓度
利用地物光谱仪器测量生态点水稻冠层光谱；进行区域采样，计算单位面积下水稻地上部干物质量；将水稻地上部干物质烘干至恒重，测定并计算水稻植株氮浓度；
s2：构建地上部干物质量模型与水稻临界氮浓度稀释模型
利用冠层土壤调节优化植被指数(osavi)对水稻地上部干物质量进行相关性分析，通过估计标准误差(se)和拟合决定系数(r2)优化方程，建立水稻地上部干物质量模型；通过分析对比不同施氮水平对不同生长阶段的水稻地上部干物质量与植株氮浓度的动态变化来构建水稻的临界氮浓度稀释模型；
s3：构建水稻氮浓度模型
利用减量精细采样法，分析植株氮浓度与归一化植被指数(ndsi(1136，1587))之间的相关性，筛选出敏感波段，建立氮浓度模型，利用水稻地上部干物质量模型与氮浓度模型对水稻氮素营养状况进行估算；
s4：计算并决定氮肥追施策略
通过光谱遥感数据，结合临界氮浓度稀释模型和实际氮浓度对水稻氮肥需求规律的直观表现，计算出氮亏缺量(nand)；并根据水稻穗期氮肥利用效率，决定氮肥追施策略。
更进一步地，在所述步骤s1中，利用地物光谱仪器测量水稻冠层光谱时，地物光谱仪器的探头垂直于冠层，距离冠层1m，每种样本选择10个均匀分布且不同行的采样点(避开边行)，每点连续记录10次，对测量结果取平均值以消除误差，对于水稻冠层光谱数据需要进行原始光谱反射值的提取。
更进一步地，在所述步骤s1中，进行水稻地上部干物质量测量时，根据各区域平均茎蘖数，从各小区内选取代表性水稻植株3～5穴，连根拔起，带回实验室，分离水稻各器官，然后高温杀青后烘干至恒重，称量地上部干物质量，折算为单位面积下水稻地上部干物质量(w)。
更进一步地，在所述步骤s1中，测量水稻植株氮浓度时，将水稻地上部干物质烘干至恒重，粉碎后过筛，通过瑞士foss全自动凯氏凯式定氮仪对样品氮浓度进行测定，根据水稻地上部干物质量计算水稻氮浓度(nt)。
更进一步地，在所述步骤s21中，水稻地上部干物质重模型如下：
w＝0.04e6.56osavi
其中，w表示水稻地上部干物质量，osavi为水稻冠层土壤调节优化植被指数。
更进一步地，在所述步骤s22中，水稻临界氮浓度稀释模型如下：
nc＝3.62w-0.29
其中，3.46为水稻地上部生物量为1tha-1时的临界氮浓度值，w为水稻地上部生物量。
更进一步地，在所述临界氮浓度稀释模型构建完毕后，采用均方根误差(rmse)和标准化均方根误差(n-rmse)对所述临界氮浓度稀释模型的可行性进行检验。
更进一步地，在所述步骤s3中，氮浓度模型如下：
nt＝6.50ndsi(1136，1587) 4.70
其中，nt为水稻实际氮浓度，ndsi(1136，1587)为1136nm和1587nm的光谱反射值的归一化植被指数。
更进一步地，在所述步骤s4中，计算出水稻穗期氮亏缺量过程包括以下步骤：
s41：根据临界氮浓度稀释模型和氮浓度模型，计算得到临界氮浓度下水稻氮积累量和实际氮浓度下水稻氮积累量；
s42：计算临界氮浓度下水稻氮积累量与实际氮浓度下水稻氮积累量的差值，得到氮亏缺量。
更进一步地，在所述步骤s4中，决定氮肥追施策略时利用区域内不同生态点的植株氮浓度和氮肥施用量，确定各生态点的水稻穗期氮素利用效率，得出水稻氮素追施量。
本发明还提供了水稻穗期氮营养监测与诊断方法的应用，将所述水稻穗期氮素亏缺诊断方法应用于水稻穗肥亏缺诊断工作中。
本发明相比现有技术具有以下优点：该水稻穗期氮营养监测与诊断方法，充分利用了遥感技术可以在田间规模化种植中快速精准无损的优势，在降低监测成本的同时，提高了监测结果的客观性，避免了传统人工调查时的误差和以点代面的不确定性；能够利用多年数据增加模型的普适性，可以通过当年的遥感数据实时监测水稻穗期氮素营养情况，同时利用氮营养指数模型获取氮亏缺状态，使得氮素施肥方案更加具体和可操作；还利用幼穗分化初期和颖花分化期的水稻遥感信息和农艺信息对诊断方法进行验证，使监测结果更客观，值得被推广用于水稻穗期氮素亏缺诊断工作中。
附图说明
图1是氮素亏缺诊断方法的流程示意图；
图2是各植被指数与水稻地上部生物量指数拟合结果比较图；
图3是基于水稻地上部生物量构建的氮浓度稀释曲线图；
图4是两波段组合的归一化植被指数ndsi(1136，1587)与水稻氮浓度的决定系数等势图；
图5是水稻地上部干物质量预测值与实测值比较；
图6是水稻氮浓度的预测值与实测值比较；
图7是基于无人机高光谱遥感采用本发明建立的水稻穗期氮素亏缺计算方法，估算的位于长江中下游研究区水稻穗分化阶段氮素营养状况分布图。
具体实施方式
下面对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
如图1所示，本实施例提供一种技术方案：水稻穗期氮营养监测与诊断方法，包括以下步骤：
s1：测量水稻穗分化阶段冠层光谱、地上部干物质量和氮浓度
利用地物光谱仪器测量生态点水稻冠层光谱；进行区域采样，计算单位面积下水稻地上部干物质量；将水稻地上部干物质烘干至恒重，测定并计算水稻氮浓度；
s2：构建地上部干物质量模型与水稻临界氮浓度稀释模型
利用冠层土壤调节优化植被指数(osavi)对水稻地上部干物质量进行相关性分析，通过估计标准误差(se)和拟合决定系数(r2)优化方程，建立水稻地上部干物质量模型；通过分析对比不同施氮水平对不同生长阶段的水稻地上部干物质量与植株氮浓度的动态变化来构建水稻临界氮浓度稀释模型；
s3：构建水稻氮浓度模型
利用减量精细采样法，分析植株氮浓度与植被指数(ndsi(1136，1587))之间的相关性，筛选出敏感波段，建立氮浓度模型，结合水稻地上干物质量模型与氮浓度模型对水稻氮素营养状况进行估算；
s4：计算并决定氮肥追施策略
通过光谱遥感数据，结合临界氮浓度稀释模型和实际氮浓度对水稻氮肥需求规律的直观表现，计算出氮亏缺量(nand)；并根据水稻穗分化阶段肥料利用效率，决定氮肥追施策略。
在所述步骤s1中，利用地物光谱仪器测量水稻冠层光谱时，地物光谱仪器的探头垂直于冠层，距离冠层1m，每种样本选择10个均匀分布且不同行的采样点(避开边行)，每点连续记录10次，对测量结果取平均值以消除误差，对于水稻冠层光谱数据需要进行原始光谱反射值提取。
在所述步骤s1中，进行水稻地上部干物质量测量时，根据各区域平均茎蘖数，从各小区内选取代表性水稻植株3～5穴，连根拔起，带回实验室，分离水稻各器官，然后高温杀青后烘干至恒重，称量地上部干物质量，折算为单位面积下水稻地上部干物质量。
在所述步骤s1中，测量水稻植株氮浓度时，将水稻地上部干物质烘干至恒重，粉碎后过筛，通过瑞士foss全自动凯氏凯式定氮仪对样品氮浓度进行测定，根据水稻地上部干物质量计算植株氮浓度(nt)。
在所述步骤s21中，利用现有的植被指数来估算水稻地上部干物质量，采用线性、指数、对数、幂函数来分析所选与生物量相关的植被指数与水稻地上部干物质量之间的关系。利用rondeaux等提出的水稻冠层土壤调节优化植被指数(osavi)与水稻地上部干物质量进行相关性分析，利用估计标准误差(se)和拟合决定系数(r2)优化方程。基于土壤调节优化植被指数(osavi)构建本地区水稻地上部干物质量模型如下：
w＝0.04e6.56osavi
其中，w表示水稻地上部干物质量，osavi为水稻冠层土壤调节优化植被指数。
表一文献中出现的用于水稻地上部干物质量估算的指数
图2为各植被指数与水稻地上部生物量指数拟合结果比较图，冠层土壤调节优化植被指数(osavi)呈现最高的决定系数值r2＝0.80。相比于其他指数，冠层土壤调节优化植被指数(osavi)表现出对水稻穗分化阶段地上部干物质量较好的敏感性。因此，基于冠层土壤调节优化植被指数(osavi)构建水稻穗分化阶段地上部干物质量模型是可行的。
在所述步骤s22中，临界氮浓度(nc)是一个广泛使用的能够描述作物最大地上部干物质量所需要的最小氮浓度，其计算公式为：
nc＝aw-b
其中，nc为物地上部干物质量对应的临界氮浓度(％)，w为作物地上部干物质量(tha-1)，a为作物地上部干物质量为1tha-1时的临界氮浓度，b为控制模型曲线斜率的数学参数。greenwood等提出了关于c3和c4植物临界氮浓度与地上部干物质量的通用定量模型方程(c3：nc＝5.7dm-0.5，c4：nc＝4.1dm-0.5)，c3供试作物有马铃薯、小麦、卷心菜、豆类等，c4供试作物有玉米、高粱等。然而以上模型是在作物生长不受氮素制约条件下构建，临界氮浓度稀释模型的构建应考虑作物生长受氮素制约情况。由于不同生态区域和品种类型，甚至同为籼稻亚种所构建的模型存在较大差异，说明该模型不具备普适性，具有明显的本地验参性。因此，需要构建区域化的临界氮浓度模型。
因此，本发明构建本地区水稻临界氮浓度稀释模型如下：
nc＝3.62w-0.29
其中，3.46为作物地上部生物量为1tha-1时水稻的临界氮浓度值，w为水稻地上部干物质量。在所述临界氮浓度稀释模型构建完毕后，采用均方根误差(rmse)和标准化均方根误差(n-rmse)对所述临界氮浓度稀释模型的可行性进行检验。
图3为基于水稻地上部生物量构建的氮浓度稀释曲线图，受氮限制处理的数据点低于或者接近临界氮浓度稀释曲线(位于曲线一侧)，不受氮限制处理的数据点高于或者接近临界氮浓度稀释曲线(位于曲线另一侧)。利用水稻每个取样日观察到植株氮浓度最大值与最小值(nmax、nmin)及对应水稻植株地上部生物量，得到两个氮浓度稀释边界模型(nmax＝3.74w-0.29、nmin＝3.48w-0.28)，同样也符合作物临界氮浓度稀释模型。因此，本发明构建的本地区水稻临界氮浓度稀释模型是可行的。
在所述步骤s3中，由于现有的植被指数对于水稻穗期植株氮浓度预测效果较差，需进一步挖掘冠层光谱敏感波段进一步提高植株氮浓度的预测效果。本发明通过减量精细采样法，系统分析植株氮浓度与全波段的双波段组合归一化植被指数(ndsi)之间的相关性，从而筛选出敏感波段。以10nm波段间隔从350nm-2500nm波段分析相关性，发现敏感双波段组合区域为800-1300nmvs1400-1800nm，再以1nm波段间隔从敏感双波段组合区域对归一化植被指数(ndsi)与水稻氮素浓度进行相关性分析，得到ndsi(1136，1587)，利用估计标准误差(se)和拟合决定系数(r2)优化方程。基于归一化植被指数(ndsi(1136，1587))构建本地区水稻氮浓度模型如下：
nt＝6.50ndsi(1136，1587) 4.70
其中，nt为水稻实际氮浓度，ndsi(1136，1587)为1136nm和1587nm的光谱反射值的归一化植被指数。
图4为两波段组合的归一化植被指数ndsi(1136，1587)与水稻氮浓度的决定系数等势图，归一化植被指数ndsi(1136，1587)在所有双波段重为最优组合。相比于其他组合，归一化植被指数ndsi(1136，1587)表现出对水稻穗分化阶段氮浓度最好的敏感性，因此，基于归一化植被指数ndsi(1136，1587)构建水稻穗分化阶段氮浓度模型是可行的。
在所述步骤s4中，利用本发明中构建的氮浓度、干物质重以及临界氮浓度稀释模型，计算出水稻穗分化阶段追氮亏缺量。通过收集不同生态点的水稻氮浓度和氮肥施用量，得到水稻穗肥氮素利用效率，按照下列计算方法得到水稻穗分化阶段追施氮肥用量。
w＝0.04e6.56osavi(1)
nc＝3.62w-0.29(2)
nt＝6.50ndsi(1136，1587) 4.70(3)
nna＝nt×w×10(4)
ncna＝nc×w×10(5)
nand＝ncna-nna(6)
nni＝nt/nc(7)
nr＝nand/nre(8)
式(1)-(8)中：
nr：水稻氮素追施用量(kgnha-1)；
nre：水稻穗肥氮素利用效率(％)，本实施例中取值60％(kgnha-1)；
nand：水稻氮亏缺量(kgnha-1)；
ncna：临界氮浓度下水稻氮积累量(kgnha-1)；
nna：实际氮浓度下水稻氮积累量(kgnha-1)；
nc：生产当前生物量所需的最低氮浓度(％)，采用公式(6)对本地区中籼稻临界氮浓度稀释曲线计算；
nt：当前生物量下实际氮浓度(％)，采用公式(7)进行光谱参数反演；
w：当前地上部生物量(tha-1)，采用公式(8)进行光谱参数反演；
nni：氮营养指数；
ndsi(1136，1587)：1136nm和1587nm的冠层原始反射光谱的归一化光谱指数；
osavi：水稻冠层土壤调节优化植被指数。
具体实施案例：
本发明提供了水稻穗分化阶段氮素亏缺诊断方法，还将所述水稻穗期氮素亏缺诊断方法应用于水稻幼穗分化初期和颖花分化期氮素亏缺诊断具体工作中。
本发明具体应用于中国安徽省皖中地区的研究区。皖中地区地处于安徽省淮河以南、大别山—巢湖以北(31°28′n，117°13′e)，属亚热带季风气候。该地区主要河流为长江和淮河，在其与巢湖的交互作用下形成了肥沃的农田土壤资源，该地区主要种植的粮食作物为水稻，包括单季稻和双季稻。因此，研究区是诊断水稻穗分化阶段氮素亏缺的典型地区，本发明以幼穗分化初期和颖花分化期的水稻为实验对象，其具体实施方案如下：
1.1、水稻冠层光谱测量
水稻冠层反射光谱采用美国asdfieldspec4地物光谱仪测量，该地物光谱仪波段范围为350-2500nm，带宽为1nm，裸光纤探头视场角为25°。实验中使用的参考板为99％的朗伯体白板。水稻冠层光谱应于晴朗无云微风天气，10:00-14:00采集数据，采集数据前地物光谱仪预热30min以上，遇多变天气则增加白板矫正次数，以确保试验精度；测量时，探头垂直于冠层，距离冠层1m，每种样本选择10个均匀分布且不同行的采样点(避开边行)，每点连续记录10次，取平均值消除误差。
1.2水稻地上部干物质重模型和临界氮浓度稀释模型的性能评估
本发明构建的水稻地上部干物质量模型为：
w＝0.04e6.56osavi
其中，w表示水稻地上部干物质量，osavi为水稻冠层土壤调节优化植被指数。利用幼穗分化初期和颖花分化期的水稻数据检验模型的可行性，采用均方根差(rmse)、预测相关性(p-r2)来评价。
图5为水稻地上部干物质量预测值与实测值比较。预测相关性(p-r2)为0.57，均方根差(rmse)为2.54。冠层土壤调节优化植被指数(osavi)能够较好地反映水稻地上部干物质量(w)变化情况，因此，基于冠层土壤调节优化植被指数(osavi)估算水稻穗分化阶段地上部干物质量是可行的。
本发明构建的本地区水稻临界氮浓度稀释模型为：
nc＝3.62w-0.29
其中，3.46为水稻地上部干物质量为1tha-1时的临界氮浓度值，w为水稻地上部干物质量。利用幼穗分化初期和颖花分化期的水稻数据检验模型的可行性，本发明构建的临界氮浓度稀释曲线可以很好地将氮限制和非氮限制组区幼穗分化初期和颖花分化期的水稻数据分开。
1.3水稻氮浓度模型性能评估
本发明构建的水稻氮浓度模型为：
nt＝6.50ndsi(1136，1587) 4.70
其中，nt为水稻实际氮浓度，ndsi(1136，1587)为1136nm和1587nm的光谱反射值的归一化植被指数。利用幼穗分化初期和颖花分化期的水稻数据检验模型的可行性，采用均方根差(rmse)、预测相关性(p-r2)来评价。
图6为水稻氮浓度的预测值与实测值比较，预测相关性(p-r2)为0.51，均方根差(rmse)为0.52。归一化植被指数ndsi(1136，1587)能够较好地反映水稻地上部干物质量(w)变化情况，因此，基于归一化植被指数ndsi(1136，1587)估算水稻穗分化阶段氮浓度模型是可行的。
1.4计算并决定氮肥追施策略
通过光谱遥感数据，结合临界氮浓度稀释模型和实际氮浓度对水稻氮肥需求规律的直观表现，计算出氮亏缺量(nand)；并根据水稻幼穗分化初期和颖花分化期阶段肥料利用效率，决定氮肥追施策略。
图7是基于无人机高光谱遥感采用本发明建立的水稻穗期氮素亏缺计算方法，估算的位于长江中下游研究区水稻穗分化阶段氮素营养状况分布图。将水稻穗分化阶段氮素亏缺值分为5个级别：0-25、25-50、50-75、75-100、大于100(图中数值单位为kgha-1)。由于施氮水平的差异，图中水稻穗分化阶段氮亏缺值的空间分布特征和异质性差异较为明显。前期施氮较低的实验区块比其他实验区块具有较高的氮素亏缺。
综上所述，上述实施例的水稻穗分化期氮亏缺诊断方法，充分利用了遥感技术可以在田间规模化种植中快速精准无损的优势，在降低监测成本的同时，提高了监测结果的客观性，避免了传统人工调查时的误差和以点代面的不确定性；能够利用多年数据增加模型的普适性，可以通过当年的遥感数据实时监测水稻穗分化期的氮素营养情况，同时利用上述计算方法得到水稻氮亏缺值，使得氮素施肥方案更加具体和可操作；还利用幼穗分化初期和颖花分化期的水稻遥感信息和农艺信息对诊断方法进行验证，结果表明，水稻穗分化期氮素亏缺诊断方法准确性较高，预测能力较强，操作便捷。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
1.基于高光谱的水稻穗期氮营养监测与诊断方法，其特征在于，包括以下步骤：
s1：测量水稻冠层光谱、地上部干物质量和氮浓度
利用地物光谱仪器测量生态点水稻冠层光谱；然后进行区域采样，计算单位面积下水稻地上部干物质量；将水稻地上部干物质烘干至恒重，测定并计算水稻氮浓度；
s2：构建水稻地上部干物质量模型与临界氮浓度稀释模型
通过分析对比不同施氮水平对不同生长阶段的水稻地上部干物质量与氮浓度的动态变化构建水稻的临界氮浓度稀释模型；对冠层土壤调节优化植被指数osavi与水稻地上部干物质量进行相关性分析，通过估计标准误差和拟合决定系数优化方程，建立水稻地上部干物质量模型；
s3：构建水稻氮浓度模型
利用减量精细采样法，分析植株氮浓度与植被指数之间的相关性，筛选出敏感波段，建立水稻氮浓度模型，利用水稻地上部干物质量模型与氮浓度模型对水稻穗期氮素营养状况进行估算；
s4：计算并决定氮肥追施策略
通过光谱遥感数据，结合临界氮浓度稀释模型和实际氮浓度对水稻氮肥需求规律的直观表现，计算出氮亏缺量；并根据水稻穗期氮肥利用效率，决定氮肥追施策略。
2.根据权利要求1所述的基于高光谱的水稻穗期氮营养监测与诊断方法，其特征在于：在所述步骤s1中，利用地物光谱仪器测量水稻冠层光谱时，地物光谱仪器的探头垂直于冠层，距离冠层指定高度，每种样本选择多个均匀分布且不同行的采样点，每点连续记录多次，对测量结果取平均值。
3.根据权利要求1所述的基于高光谱的水稻穗期氮营养监测与诊断方法，其特征在于：在所述步骤s1中，进行水稻地上部干物质量测量时，根据各区域平均茎蘖数，从各小区内选取代表性水稻植株3～5穴，连根拔起，带回实验室，分离水稻各器官，然后高温杀青后烘干至恒重，称量地上部干物质量，折算单位面积下水稻地上部干物质量。
4.根据权利要求1所述的基于高光谱的水稻穗期氮营养监测与诊断方法，其特征在于：在所述步骤s1中，测量水稻植株氮浓度时，将水稻地上部干物质烘干至恒重，粉碎后过筛，通过瑞士foss全自动凯氏凯式定氮仪对样品氮浓度进行测定，根据水稻地上部干物质量计算植株氮浓度nt。
5.根据权利要求1所述的基于高光谱的水稻穗期氮营养监测与诊断方法，其特征在于：在所述步骤s2中，包括以下步骤：
s21水稻地上部干物质量模型如下：
w＝0.04e6.56osavi
其中，w表示水稻地上部干物质量，osavi为水稻冠层土壤调节优化植被指数；
s22临界氮浓度稀释模型如下：
nc＝3.62w-0.29
其中，3.46为作物地上部生物量为1tha-1时的水稻临界氮浓度值，w为水稻地上部生物量。
6.根据权利要求5所述的基于高光谱的水稻穗期氮营养监测与诊断方法，其特征在于：在所述地上部干物质量模型与临界氮浓度稀释模型构建完毕后，采用均方根误差rmse和标准化均方根误差n-rmse对所述临界氮浓度稀释模型的可行性进行检验。
7.根据权利要求1所述的基于高光谱的水稻穗期氮营养监测与诊断方法，其特征在于：在所述步骤s3中，氮浓度模型如下：
nt＝6.50ndsi(1136，1587) 4.70
其中，nt为水稻实际氮浓度，ndsi(1136，1587)为1136nm和1587nm的光谱反射值的归一化植被指数。
8.根据权利要求1所述的基于高光谱的水稻穗期氮营养监测与诊断方法，其特征在于：在所述步骤s4中，计算出水稻穗期氮亏缺量过程包括以下步骤：
s41：根据临界氮浓度稀释模型和氮浓度模型，计算得到临界氮浓度下水稻氮积累量和实际氮浓度下水稻氮积累量；
s42：计算临界氮浓度下水稻氮积累量与实际氮浓度下水稻氮积累量的差值，得到氮亏缺量。
9.根据权利要求8所述的基于高光谱的水稻穗期氮营养监测与诊断方法，其特征在于：在所述步骤s4中，决定氮肥追施策略时利用区域内不同生态点的植株氮浓度和氮肥施用量，确定各生态点的水稻穗肥氮素利用效率，得出水稻氮素追施量。
10.一种如权利要求1～9任一项所述的基于高光谱的水稻穗期氮营养监测与诊断方法的应用，其特征在于：将所述水稻穗期氮素营养诊断方法应用于水稻穗肥亏缺诊断工作中。
技术总结