专利名称：一种用于检测作物冠层叶片碳氮比的光谱指数构造方法
技术领域：
本发明涉及一种光谱指数构造方法，尤指一种用于检测作物冠层叶片碳氮比的光谱指数构造方法。
背景技术：
碳与氮是作物植株的两个基本元素，它们的代谢状况对作物产量与品质的形成具有重要影响。作物叶片碳氮比(叶片碳含量与氮含量的比值)能够说明碳与氮二者的代谢状况，是一个综合诊断作物植株体内碳氮平衡、营养状况、生长活性和抗病性的重要指标，因此实时、有效和准确地检测作物叶片碳氮比对于田间作物长势诊断和栽培管理具有重要意义。常规的作物冠层叶片碳氮比的检测，通常是采用田间取样、实验室化验的方式实现，检测不仅耗时费力、成本高，而且具有破坏性和滞后性。近年来，随着高光谱遥感技术的快速发展，应用光谱无损探测技术对诸如作物叶片碳、氮含量的生理生化组分进行无损、快速的检测，已成为当前田间作物长势、营养诊断应用的重要内容。当前，应用高光谱技术检测作物叶片氮含量的光谱指数方法很多，但用于检测碳氮比的光谱指数方法的报道却不多见，已报道的少数几个用于检测碳氮比的光谱指数，也多是使用两三个特征波长的原始反射率通过简单地数学运算得到，但由于使用的波长较少，相邻波长原始反射率间存在着高度相关性，使得构造的光谱指数用于碳氮比检测时缺乏应有的稳定性。此外，已有的光谱指数使用的是原始反射率信息，而在实际应用中，田间获取作物冠层光谱时易受天空局部云斑阴影暂时或偶然性遮挡，以及一天中不同时间段太阳高度角变化等因素造成的光照条件差异的影响，使得用原始反射率计算的光谱指数往往携带一定的测量误差，用于作物叶片碳氮比检测时必然会产生不必要的检测误差。

发明内容
(一)要解决的技术问题针对现有技术中的上述缺陷，本发明要解决的技术问题是提供一种用于检测作物冠层叶片碳氮比的光谱指数构造方法，其能够挖掘高光谱数据中丰富的多波段反射率信息，避免已有光谱指数因使用的波长数过少，相邻波长原始反射率间存在高度相关性，导致检测碳氮比时产生的不稳定问题，同时在一定程度上消除由偶然性云斑阴影或不同测量时段太阳高度角变化等因素产生的光照条件差异造成的光谱测量误差问题，以提高所构造的光谱指数对作物冠层叶片碳氮比检测的准确性和可比性。(二)技术方案为解决上述问题，本发明提供了一种用于检测作物冠层叶片碳氮比的光谱指数构造方法，所述方法包括步骤(I)获取作物冠层500nm 910nm波段范围的高光谱原始反射率数据；(2)利用获取的作物冠层原始反射率计算相对反射率；
(3)计算相对反射率曲线中500nm 550nm、680nm 760nm以及760nm 9IOnm
三个子波段范围光谱曲线的斜率；(4)构建用于检测作物冠层叶片碳氮比的光谱指数RCN。进一步,步骤(I)按以下方法获取500nm 910nm波段范围的高光谱原始反射率数据使用可测量作物冠层500nm 910nm波段范围光谱的便携式地物高光谱测量仪，在田间直接测量作物冠层辐亮度或DN值数据，则通过相应的参考板校正参数计算得到作物冠层高光谱原始反射率。进一步,步骤(I)按以下方法获取500nm 910nm波段范围的高光谱原始反射率数据使用包含了 500nm 910nm波段范围的高光谱遥感影像数据，则通过相应的辐射定标与大气校正处理，得到地表作物冠层的高光谱原始反射率数据。进一步，步骤(2)按以下方法计算相对反射率对于获取的由500nm 910nm波段范围内各波长对应的原始反射率形成的光谱曲线，首先进行三点平滑计算得到平滑反射率，然后计算平滑处理后的该光谱曲线反射率的均值，紧接着使用平滑处理后的光谱曲线上各波长对应的反射率除以该均值，所得结果即为500nm 910nm波段范围内各波长对应的相对反射率。进一步，所述平滑反射率的计算公式为Ri=OvAi^iv1)/3，Ri指三点平滑处理后第i个波长所对应的平滑反射率，指500nm 910nm波段范围内第i个波长所对应的原
始反射率；平滑处理后反射率的均值的计算公式为< n指500nm 910nm波段范
围内波长的个数，A指平滑处理后反射率的均值；所述500nm 910nm波段范围内各波长对应的相对反射率的计算公式为R' ^RiA,计算得到的第i个波长所对应的相对反射率。进一步，步骤(3)按以下方法分别计算相对反射率光谱曲线中500nm 550nm、680nm 760nm以及760 910nm三个子波段范围光谱曲线的斜率以各子波段范围内的波长为自变量，以各波长所对应的相对反射率为因变量，通过一次线性拟合的方式计算得到三个子波段范围的光谱曲线的斜率k。进一步，所述一次线性拟合的计算公式为R’ =b+kX入，式中，R’为相对反射率；入为波长；b为线性拟合方程的截距；k为光谱曲线的斜率。进一步，步骤(4)按以下方法计算可用于检测作物冠层叶片碳氮比的光谱指数RCN :记步骤(3)中的500nm 550nm、680nm 760nm以及760 910nm三个波段范围相对反射光谱曲线的斜率分别为kg、kr和kn，则RCN的计算公式为RCN = (kg+kn) / (2 X kr)。进一步，所述500nm 910nm波段范围的高光谱原始反射率数据是指光谱分辨率小于IOnm波段宽度的光谱数据。(三)有益效果本发明构造的用于检测作物冠层叶片碳氮比的光谱指数具有以下优点(I)避免相邻波长反射率间存在的高度相关性问题。本发明构造的光谱指数使用了 500nm 550nm、680nm 760nm以及760 910nm三个波段范围的光谱数据信息，每个光谱范围都包含数十个波长，反射率信息更加丰富，从而可避免现有的只用少数波长反射率所构造的光谱指数，因相邻波长反射率间存在高度相关性而产生检测碳氮比不稳定的问题。(2)消除光照条件差异的影响。由于光谱测量一般要求晴天无云的天气条件，但在实际田间光谱测量过程中，临时或偶然性天空中云斑对光照的遮挡，以及一天不同时间段太阳高度角变化等原因，往往使得实际田间作物冠层的光照条件发生变化，致使从田间获得的光谱反射率数据携带测量误差，使用这些情况下获取的光谱原始反射率在一定程度上缺乏相互间的可比性和同质性，势必造成计算的光谱指数用于碳氮比的检测时会产生不必要的误差。而本发明由于使用的是相对反射率，可在一定程度上消除不同光照差异条件下获取的原始反射率间的异质性，增强检测叶片碳氮比时光谱数据间的可比性。(3)适用性强。本发明所构造的光谱指数主要使用500nm 550nm、680nm 760nm以及760 9 IOnm三个波段范围的高光谱数据，这三个范围涉及绿光、红边和近红外高反射平台光谱波段区域，也是与作物生长紧密相关的特征波段范围，既然碳与氮是作物生长器官的基本组成元素，其代谢变化必然也会在多个反射光谱波段范围有所反应，而不只是表现在两三个特征波长反射率的变化响应上，因此以反射率曲线的斜率作为光谱指数的输入参数，既充分应用了多个光谱波长所蕴含的丰富的光谱响应信息，又以简单明了的表现形式来实现，计算方法简单，使得本发明所构造的光谱指数用于检验作物冠层叶片碳氮比具有更好的扩展性和稳定性，适用性较强。


图1为本发明所述用于检测作物冠层叶片碳氮比的光谱指数构造方法流程图；图2为应用光谱指数RCN分别估测冬小麦与大麦叶片碳氮比的效果；图3为不考虑作物类别时RCN检测作物叶片碳氮比的效果。
具体实施例方式下面结合附图和实施例，对本发明的具体实施方式
作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。如图1所示，本发明所述用于检测作物冠层叶片碳氮比的光谱指数构造方法，包括步骤(I)获取作物冠层500nm 9IOnm波段范围的高光谱原始反射率数据如果使用可测量作物冠层500nm 910nm波段范围光谱的便携式地物高光谱测量仪，在田间直接测量作物冠层辐亮度或DN值数据，则通过相应的参考板校正参数计算得到作物冠层高光谱原始反射率；如果使用包含了 500nm 910nm波段范围的高光谱遥感影像数据，则通过相应的辐射定标与大气校正处理，得到地表作物冠层的高光谱原始反射率数据。例如，构造用于检测冬小麦和大麦作物冠层叶片碳氮比的光谱指数RCN，本步骤中应用美国ASD (Analytical Spectral Device)公司生产的便携式地物野外地物高光谱仪，于2010年7月在内蒙古自治区海拉尔农垦大麦田间测量大麦拔节至抽穗期的冠层辐亮度或DN值数据；于2012年4月至5月间在北京昌平小汤山测量了冬小麦拔节到灌浆期间的冠层辐亮度或DN值数据，并通过相应的参考板校正参数各自计算获得两个作物冠层500nm 9IOnm波段范围高光谱原始反射率。(2)利用获取的作物冠层原始反射率计算相对反射率；
本步骤中，对于已获取的冬小麦和大麦500nm 910nm波段范围内各波长对应的原始反射率形成的光谱曲线，首先进行三点平滑计算得到平滑反射率(如公式(I)所示)，然后计算平滑处理后的该光谱曲线反射率的均值(如公式(2)所示)，紧接着使用平滑处理后的光谱曲线上各波长对应的反射率除以该均值，所得结果即为500nm 910nm波段范围内各波长对应的相对反射率(如公式(3)所不)。
权利要求
1.一种用于检测作物冠层叶片碳氮比的光谱指数构造方法，其特征在于，所述光谱指数构造方法包括以下步骤 (O获取作物冠层500nm 910nm波段范围的高光谱原始反射率数据； (2)利用获取的作物冠层原始反射率计算相对反射率； (3)计算相对反射率曲线中500nm 550nm、680nm 760nm以及760nm 910nm三个子波段范围光谱曲线的斜率； (4)构建用于检测作物冠层叶片碳氮比的光谱指数RCN。
2.如权利要求1所述的光谱指数构造方法，其特征在于，步骤(I)按以下方法获取500nm 910nm波段范围的高光谱原始反射率数据使用可测量作物冠层500nm 910nm波段范围光谱的便携式地物高光谱测量仪，在田间直接测量作物冠层辐亮度或DN值数据，则通过相应的参考板校正参数计算得到作物冠层高光谱原始反射率。
3.如权利要求1所述的光谱指数构造方法，其特征在于，步骤(I)按以下方法获取500nm 910nm波段范围的高光谱原始反射率数据使用包含了 500nm 910nm波段范围的高光谱遥感影像数据，则通过相应的辐射定标与大气校正处理，得到地表作物冠层的高光谱原始反射率数据。
4.如权利要求1所述的光谱指数构造方法，其特征在于，步骤(2)按以下方法计算相对反射率对于获取的由500nm 910nm波段范围内各波长对应的原始反射率形成的光谱曲线，首先进行三点平滑计算得到平滑反射率，然后计算平滑处理后的该光谱曲线反射率的均值，紧接着使用平滑处理后的光谱曲线上各波长对应的反射率除以该均值，所得结果即为500nm 910nm波段范围内各波长对应的相对反射率。
5.如权利要求4所述的光谱指数构造方法，其特征在于，所述平滑反射率的计算公式为Ri= (IVJrdiV1) /3，Ri指三点平滑处理后第i个波长所对应的平滑反射率，巧指500nm 910nm波段范围内第i个波长所对应的原始反射率；平滑处理后反射率的均值的计算公式为A = 4 R1，n指500nm 9IOnm波段范围内波长的个数，A指平滑处理后反射率的均值；所述500nm 910nm波段范围内各波长对应的相对反射率的计算公式为R' ^RiA, Ri i指计算得到的第i个波长所对应的相对反射率。
6.如权利要求1所述的光谱指数构造方法，其特征在于，步骤(3)按以下方法分别计算相对反射率光谱曲线中500nm 550nm、680nm 760nm以及760 910nm三个子波段范围光谱曲线的斜率以各子波段范围内的波长为自变量，以各波长所对应的相对反射率为因变量，通过一次线性拟合的方式计算得到三个子波段范围的光谱曲线的斜率k。
7.如权利要求6所述的光谱指数构造方法，其特征在于，所述一次线性拟合的计算公式为R’ =b+kX λ，式中，R’为相对反射率；λ为波长；b为线性拟合方程的截距；k为光谱曲线的斜率。
8.如权利要求1所述的光谱指数构造方法，其特征在于，步骤(4)按以下方法计算可用于检测作物冠层叶片碳氮比的光谱指数RCN :记步骤(3)中的500nm 550nm、680nm 760nm以及760 910nm三个波段范围相对反射光谱曲线的斜率分别为kg、kr和kn，则RCN的计算公式为 RCN = (kg+kn)/(2Xkr)。
9.如权利要求1所述的光谱指数构造方法，其特征在于，所述500nm 910nm波段范围的高光谱原 始反射率数据是指光谱分辨率小于IOnm波段宽度的光谱数据。
全文摘要
本发明公开了一种用于检测作物冠层叶片碳氮比的光谱指数构造方法，所述方法包括步骤(1)获取作物冠层500nm～910nm波段范围的高光谱原始反射率数据；(2)利用获取的作物冠层原始反射率计算相对反射率；(3)计算相对反射率光谱曲线中500nm～550nm、680nm～760nm以及760nm～910nm三个子波段范围光谱曲线的斜率；(4)构建用于检测作物冠层叶片碳氮比的光谱指数RCN。采用本发明构造的光谱指数进行作物冠层叶片碳氮比的检测，可以消除光照条件差异对碳氮比检测的影响，同时还可避免已有光谱指数所存在的相邻波长反射率间存在高度相关的问题。
文档编号G01N21/25GK103048276SQ201210545190
公开日2013年4月17日 申请日期2012年12月14日 优先权日2012年12月14日
发明者徐新刚, 赵春江, 宋晓宇, 杨贵军, 杨小冬, 顾晓鹤, 杨浩, 龙慧灵, 董燕生, 张竞成, 冯海宽 申请人:北京农业信息技术研究中心