一种冬小麦识别方法
【专利摘要】本发明公开了一种冬小麦识别方法。利用冬小麦生育期内的遥感植被指数时序数据和物候历数据，通过分类类别设定、训练样本选取、样本时序曲线提取、冬小麦识别特征选择及参量化、最佳特征阈值确定和冬小麦识别模型构建，提取出区域范围内冬小麦的空间分布信息。本方法的特点是可基于未经去噪声处理的遥感植被指数时序数据进行冬小麦空间分布信息提取，具有较强的抗时序数据噪声能力；方法的稳定性和普适性较高，可应用于具有一定物候差异的大范围区域的冬小麦分布信息遥感提取。
【专利说明】
-种冬小麦识别方法
技术领域
[0001] 本发明设及农业遥感技术领域，尤其设及一种冬小麦空间分布信息的遥感提取方 法。
【背景技术】
[0002] 遥感影像数据因具有覆盖面积大、现势性强等优势，已被广泛应用于农作物空间 分布信息的提取工作。目前，基于遥感影像的作物识别方法可分为单期影像源方法和时间 序列影像源方法。
[0003] 单期影像源方法是基于最佳时期的单期或少数几期影像，利用待识别作物与背景 地物的光谱差异来实现作物空间分布信息的提取。此类方法的效率高、操作性强，但也存在 明显不足。一方面，此类方法多采用中高空间分辨率遥感数据，然而，此类数据的卫星重访 周期一般较长，且受云雨天气的影响，往往难W提供最佳时期的影像数据;另一方面，在某 几个时相内，此类方法很难区分光谱相似的地物类型。因此，此类方法的应用通常要受到遥 感数据源的限制。
[0004] 时间序列影像源方法是基于植被指数时间序列数据，利用作物的物候变化规律特 征W实现作物空间分布信息的提取。不同作物类型的物候特征具有独特性，对不同地物类 型具有可区分性，能够进一步提高作物的识别精度。植被指数时序数据由于能够充分反映 不同作物的物候节律信息，从而使基于植被指数时序数据的作物空间信息提取方法得到了 不断发展，并逐渐成为目前最为主流的技术方法之一。
[0005] 传统的基于时间序列的作物识别方法多利用作物在几个关键生长节点上的影像 特征及其与其它节点上影像特征之间的逻辑关系来构建作物识别模型。然而，受时序数据 噪声的影响，某些待识别作物像元的时序曲线在关键节点上的影像特征易发生突变，甚至 与其他节点上影像特征之间的逻辑关系也发生了改变，致使传统方法的作物识别精度较 低。因此，传统的基于时序数据的作物识别方法往往需要依靠高质量的植被指数时间序列 数据才能获得较高的作物分布信息提取精度，但受大气条件(云雨天气等）的影响，我们有 时无法获取质量普遍较高的时序数据，通常需要在应用之前对其进行去噪声处理。
[0006] 去噪声处理一方面确实能够去掉一些比较明显的噪声，使时序数据变得更为平 滑，但另一方面，去噪声处理也使时序数据所反映的农作物物候特征信息被修改或作为噪 声信息而被去除（如冬小麦的双峰特征不明显），并同时会混进新的噪声，从而影响对后续 作物识别特征的选择。此外，去噪声处理也增加了作物识别工作的复杂度，降低了工作效 率。因此，在基于时序数据进行作物空间信息提取时，应尽量避免对时序数据进行去噪声处 理。
【发明内容】
[0007] 本发明正是针对传统的基于时间序列的作物识别方法的一些不足，提出了一种抗 遥感时序数据噪声的冬小麦识别方法，旨在利用未经滤波处理的植被指数时序数据，快速 提取区域范围内冬小麦的空间分布信息。
[000引为达此目的，本发明采用W下技术方案：
[0009] -种冬小麦识别方法，包括W下步骤：
[0010] A、分类类别设定
[0011] 根据各地物类型与冬小麦的混分状况，设定区域范围内的分类类别。
[0012] B、各地物类型训练样本选取
[0013] 针对设定的分类类别，基于植被指数时序数据或符合应用需求的其它数据，利用 目视解译的方法，遵循基本的采样原则，分别采集各地物类型的训练样本;采样的基本原则 如下：
[0014] (1)各地物类别的样本数量一般要与区域范围内各地物类型的分布面积成比例。
[0015] (2)各地物类型的样本要在整个区域范围内均匀分布，使各个样本在局部范围内 具有区域代表性。
[0016] (3)各地物类型的样本应位于纯像元区域，且其周围像元也应尽量为纯像元。
[0017] C、样本时序曲线提取
[0018] 基于植被指数时序数据，利用采集的训练样本，提取各地物类型对应的样本时序 曲线。
[0019] D、冬小麦识别特征选择及参量化
[0020] 结合区域的物候历数据，根据提取的样本时序曲线，确定冬小麦的识别特征，并对 其进行参量化。确定的识别特征和参量化结果为：
[0021] (1)识别特征。冬小麦的识别特征为冬小麦分葉期内的生长状态、冬小麦出苗期- S叶期和开花期-乳熟期内各自的总体生长趋势。
[0022] (2)参量化结果。分别将上述冬小麦的识别特征参量化为冬小麦分葉期内植被指 数的平均值、冬小麦出苗期叶期和开花期-乳熟期内各自植被指数的线性拟合斜率。
[0023] E、最佳特征阔值确定
[0024] 首先基于训练样本，对参量化后各地物类型的识别特征分别进行频数统计，然后 根据统计结果，将冬小麦与其它各地物类型区分程度最高处所对应的特征取值范围作为相 应识别特征的最佳阔值选择区间，最后在确定的区间范围内利用等步长逐步捜索的方式， 选定冬小麦识别精度最高时各识别特征所对应的取值作为其各自的最佳阔值。
[0025] F、冬小麦识别模型构建
[0026] 根据参量化后的冬小麦识别特征及其对应的最佳特征阔值，构建冬小麦的识别模 型。其具体的模型构建步骤为：
[0027] (1)利用冬小麦分葉期内植被指数的平均值排除与冬小麦物候特征差异较大的非 植被和其它作物类型。其数学模型为：
[002引
[0029] (2)利用冬小麦出苗期-S叶期和开花期-乳熟期两个时间段内植被指数的线性拟 合斜率进一巧排除其它植被类型，最终实现冬小麦的提取。其数学模型为：
[0030]
[0031]
[0032] 上述数学模型中，tia = l ,2,...6)分别代表冬小麦出苗期普遍开始时间、S叶期 普遍结束时间、分葉期普遍开始时间、分葉期普遍结束时间、开花期普遍开始时间、乳熟期 普遍结束时间；巧,表示某一像元在ti至ti + i时间段内包含的所有植被指数数据； WC0"仍J表示巧,的平均值；s/we(F/,,_,,+i)表示巧,_,,,1线性拟合后获得的斜率;Ta、Tb 和Tc分别代表冬小麦各识别特征对应的最佳阔值。将上述3个数学公式得到的结果求交集 之后的空间分布区即为识别的冬小麦分布范围。
[0033] G、冬小麦空间分布信息提取
[0034] 基于植被指数时序数据，利用构建的冬小麦识别模型，提取出区域范围内冬小麦 的空间分布信息。
[0035] 本发明具有W下特点：
[0036] (1)原理简单，易实施，工作效率高。
[0037] (2)可基于未经去噪声处理的遥感植被指数时序数据进行冬小麦空间分布信息的 提取，具有较强的抗时序数据噪声能力，且提取结果的稳定性较高。
[0038] (3)可应用于具有一定物候差异的大范围区域的冬小麦空间分布信息遥感提取。
【附图说明】
[0039] 图1是河南省冬小麦的空间分布信息提取结果。
【具体实施方式】
[0040] 下面结合附图对本发明的技术实施方案作进一步描述。
[0041] A、分类类别设定
[0042] 根据各地物类型与冬小麦的混分状况，设定区域范围内的分类类别。
[0043] 本案例中，设定河南省的分类类别为冬小麦、林地、草地、非植被和其它作物5种。
[0044] B、各地物类型训练样本选取
[0045] 针对设定的采样类别，基于植被指数时序数据或符合应用需求的其它数据，利用 目视解译的方法，遵循基本的采样原则，分别采集各地物类别的训练样本;采样的基本原则 如下：
[0046] (1)各地物类别的样本数量一般要与区域范围内各地物类型的分布面积成比例。
[0047] (2)各地物类型的样本要在整个区域范围内均匀分布，使各个样本在局部范围内 具有区域代表性。
[004引（3)各地物类型的样本应位于纯像元区域，且其周围像元也应尽量为纯像元。
[0049] 在本案例中，基于8天合成的MOD I S (Modera t e Re S O 1U t i on Imag i ng SpectrometerWDVKNormalized Difference Vege1:ation Index)时序数据，同时参考了 高分辨率遥感影像数据和物候历数据，目视解译选取了河南省各地物类别的训练样本。
[0050] C、样本时序曲线提取
[0051] 基于植被指数时序数据，利用采集的训练样本，提取各地物类型对应的样本时序 曲线。
[0052] 本案例中，基于NDVI时序数据，利用各地物类型训练样本的坐标数据，提取出其对 应像元的NDVI时序曲线作为各地物类型的样本时序曲线。
[0053] D、冬小麦识别特征选择及参量化
[0054] 结合冬小麦的物候历数据，根据提取的冬小麦样本时序曲线，确定冬小麦的识别 特征，并对其进行参量化。确定的识别特征和参量化结果为：
[0055] (1)识别特征。冬小麦的识别特征为冬小麦分葉期内的生长状态、冬小麦出苗期- S叶期和开花期-乳熟期内各自的总体生长趋势。
[0056] (2)参量化结果。分别将上述冬小麦的识别特征参量化为冬小麦分葉期内植被指 数的平均值、冬小麦出苗期叶期和开花期-乳熟期内各自植被指数的线性拟合斜率。
[0057] E、最佳特征阔值确定
[0058] 首先基于训练样本，对参量化后各地物类型的识别特征分别进行频数统计，然后 根据统计结果，将冬小麦与其它各地物类型区分程度最高处所对应的特征取值范围作为相 应识别特征的最佳阔值选择区间，最后在确定的区间范围内利用等步长逐步捜索的方式， 选定冬小麦识别精度最高时各识别特征所对应的取值作为其各自的最佳阔值。
[0059] 本案例中，初步确定分葉期NDVI平均值、出苗期-S叶期NDVI线性拟合斜率和开花 期-乳熟期NDVI线性拟合斜率S个冬小麦识别特征的最佳阔值选择区间为[0.32，0.3引、[- 0.02，0.03]和[-0.04，-0.Ol]。Wo.Ol为步长，通过等步长逐步捜索的方式确定上述识别特 征的最佳阔值依次为0.36、0.00和-0.03。
[0060] F、冬小麦识别模型构建
[0061] 根据参量化后的冬小麦识别特征及其对应的最佳特征阔值，构建冬小麦的识别模 型。其具体的模型构建步骤为：
[0062] (1)利用冬小麦分葉期内植被指数的平均值排除与冬小麦物候特征差异较大的非 植被和其它化物类巧。巧敬学模型为：
[0063]
[0064] (2)利用冬小麦出苗期-S叶期和开花期-乳熟期两个时间段内植被指数的线性拟 合斜率进一步排除其它植被类型，最终实现冬小麦的提取。其数学模型为：
[00 化]
[0066]
[0067] 上述数学模型中，tia = l，2，...6)分别代表冬小麦出苗期普遍开始时间、S叶期 普遍结束时间、分葉期普遍开始时间、分葉期普遍结束时间、开花期普遍开始时间、乳熟期 普遍结束时间；表示某一像元在ti至ti + i时间段内包含的所有植被指数数据； we。柳7,,J表示巧的平均值；J表示1^/,线性拟合后获得的斜率;Ta、Tb 和Tc分别代表冬小麦各识别特征对应的最佳阔值。将上述3个数学公式得到的结果求交集 之后的空间分布区即为识别的冬小麦分布范围。
[0068] 本案例中，参考NDVI时序数据的合成周期、样本时序曲线W及河南省物候历，选定 ti(i = l，2，. . .6)分别为：10月24日、12月3日、12月11日、12月27日、次年的4月23日和次年 的6月10日。
[0069] G、冬小麦空间分布信息提取
[0070] 基于植被指数时序数据，利用构建的冬小麦识别模型，提取出区域范围内冬小麦 的空间分布信息。
[0071] W上所述，仅为本发明较佳的【具体实施方式】，但本发明的保护范围并不局限于此， 任何熟悉该技术的人在本发明所掲露的技术范围内，可轻易想到的变化或替换，都应涵盖 在本发明的保护范围之内。因此，本发明的保护范围应该W权利要求的保护范围为准。
【主权项】
1. 一种冬小麦识别方法，其特征在于包括如下步骤： A、 分类类别设定，根据各地物类型与冬小麦的混分状况，设定区域范围内的分类类别； B、 各地物类型训练样本选取，针对设定的分类类别，基于植被指数时序数据或符合应 用需求的其它数据，利用目视解译的方法，遵循基本的采样原则，分别采集各地物类型的训 练样本； C、 样本时序曲线提取，基于植被指数时序数据，利用采集的训练样本，提取各地物类型 对应的样本时序曲线； D、 冬小麦识别特征选择及参量化，结合区域的物候历数据，根据提取的样本时序曲线， 确定冬小麦的识别特征，并对其进行参量化； E、 最佳特征阈值确定，首先基于训练样本，对参量化后各地物类型的识别特征分别进 行频数统计，然后根据统计结果，将冬小麦与其它各地物类型区分程度最高处所对应的特 征取值范围作为相应识别特征的最佳阈值选择区间，最后在确定的区间范围内利用等步长 逐步搜索的方式，选定冬小麦识别精度最高时各识别特征所对应的取值作为其各自的最佳 阈值。 F、 冬小麦识别模型构建，根据参量化后的冬小麦识别特征及其对应的最佳特征阈值， 构建冬小麦识别模型； G、 冬小麦空间分布信息提取，基于植被指数时序数据，利用构建的冬小麦识别模型，提 取出区域范围内冬小麦的空间分布信息。2. 如权利要求1步骤D中所说的冬小麦识别特征选择及参量化，其特征在于将冬小麦分 蘖期内的生长状态、冬小麦出苗期-三叶期和开花期-乳熟期内各自的总体生长趋势作为冬 小麦的识别特征，并分别将其参量化为冬小麦分蘖期内植被指数的平均值、冬小麦出苗期-三叶期和开花期-乳熟期内各自植被指数的线性拟合斜率。3. 如权利要求1步骤F中所说的冬小麦识别模型构建，其特征在于将以下3个数学公式 得到的结果求交集之后的空间分布区即为识别的冬小麦分布范围。3个数学公式分别为：式中，ti(i = 1，2，... 6)分别代表冬小麦出苗期普遍开始时间、三叶期普遍结束时间、分 蘖期普遍开始时间、分蘖期普遍结束时间、开花期普遍开始时间、乳熟期普遍结束时间； 7人,-，,+1表示某一像元在ti至ti+i时间段内包含的所有植被指数(Vegetation Index，VI)数 据;)表示F/,,~的平均值;AM%-,,+,)表示对％ a,线性拟合后所获得的斜率； Ta、Tb和Tc分别代表冬小麦各识别特征对应的最佳阈值。
【文档编号】G06K9/00GK106022224SQ201610307697
【公开日】2016年10月12日
【申请日】2016年5月12日
【发明人】朱文泉, 姜涛, 唐珂, 詹培
【申请人】北京师范大学