[0038] 上述数学模型中,Mean()表示遥感像元在时间上的求均值操作,NDVI、RED和NIR分别代表归一化差值植被指数、红光反射率和近红外反射率影像;下标文字代表春玉米的物候期;T1到T6分别代表春玉米各识别特征对应的阈值。将上述6个数学模型得到的结果求交集之后的空间分布区即为识别的春玉米分布范围。
[0039] G、春玉米空间分布信息提取
[0040] 基于区域的遥感红光反射率、近红外反射率以及归一化差值植被指数时间序列数据,利用构建的春玉米识别模型,提取区域内春玉米的空间分布信息。
[0041] 本发明具有以下特点:
[0042] (1)原理简单,易实施,工作效率高。
[0043] (2)可基于遥感红光和近红外反射率、NDVI时间序列数据进行春玉米空间分布信息的提取,且提取结果的稳定性较高。
[0044] (3)可应用于具有一定物候差异的大范围区域的春玉米空间分布信息遥感提取。
附图说明
[0045] 图1是辽宁省春玉米的空间分布信息提取结果。
具体实施方式
[0046] 下面结合附图对本发明的技术实施方案作进一步描述。
[0047] A、分类体系构建
[0048] 根据区域范围内的地物类型结构和春玉米与其他地物类型的遥感光谱混分情况,构建相应的分类体系。
[0049] 本案例设定辽宁省的地物类型为春玉米、水稻、林地(含灌木)、草地(含其他类型的作物,如蔬菜、大豆等)、建设用地、水域六类。
[0050] B、各地物类型训练样本选取
[0051] 根据所构建的分类体系和当地主要农作物物候历,基于高空间分辨率遥感影像和实地调查数据,遵循训练样本选取的基本原则,选取各地物类型的训练样本。训练样本选取的基本原则为:
[0052] (1)各地物类型的样本数量要足够多,且一般要与区域范围内各地物类型的面积成比例。
[0053] (2)各地物类型的训练样本要在整个区域范围内均匀分布,使各训练样本在局部范围内均具有代表性。
[0054] (3)各地物类型的训练样本应位于纯净像元区域,且其周围像元也应尽量为纯净像元。
[0055] 本案例同时参考了高分辨率遥感影像数据、物候历数据和实地调查数据等,目视解译选取了辽宁省各地物类型的训练样本。
[0056] C、各地物类型遥感属性指标时间序列曲线均值提取
[0057] 基于各地物类型训练样本数据,提取各地物类型训练样本在遥感红光反射率、近红外反射率以及归一化差值植被指数数据上的时间序列曲线均值。
[0058] 本案例主要使用8天合成的MODIS(Moderate Resolution Imaging Spectrometer)地表反射率数据产品(MOD09Q1),分别提取各地物类型训练样本对应像元的红光反射率、近红外反射率和NDVI时间序列曲线,并取其均值作为各地物类型的样本时序曲线。
[0059] D、识别特征选择
[0060] 结合区域的物候历数据,根据提取的各地物类型遥感属性指标时间序列曲线均值,选择各地物类型遥感属性指标时间序列曲线均值在春玉米特定物候期的差异作为识别特征,具体包括以下两个步骤:
[0061] (1)物候特征选择。分析各地物类型在遥感红光反射率、近红外反射率和NDVI时间序列数据的变化规律,确定特征差异明显的物候阶段为春玉米拔节期前期、抽雄期至乳熟期和成熟期前期。本案例确定辽宁省春玉米的拔节期前期为年日序第129–153天、抽雄期至乳熟期为年日序第193–257天、成熟期前期为年日序第257–265天。
[0062] (2)识别特征选择。分别将上述春玉米物候特征所对应的遥感属性指标参量化,确定的识别特征为:春玉米抽雄期至乳熟期、拔节期前期、成熟期前期的NDVI均值,拔节期前期的红光反射率均值,拔节期前期和成熟期前期的近红外反射率均值。本案例确定辽宁省的识别特征为:春玉米抽雄期至乳熟期年日序第193–257天,拔节期前期年日序第129–153天,成熟期前期年日序第257–265天的NDVI均值,拔节期前期年日序第129–153天的红光反射率均值,拔节期前期年日序第129–153天和成熟期前期年日序第257–265天的近红外反射率均值。
[0063] E、特征阈值确定
[0064] 根据选择的识别特征,统计各识别特征所对应的各地物类型训练样本的频数分布,根据频数分布图,将各地物类型区分程度最大处所对应的特征取值作为相应识别特征的阈值。
[0065] 本案例通过对训练样本的频数分布进行统计,确定辽宁省春玉米抽雄期至乳熟期、拔节期前期、成熟期前期的NDVI均值,拔节期前期的红光反射率均值,拔节期前期和成熟期前期的近红外反射率均值阈值依次为0.40、0.33、0.57、0.10、0.18和0.30。
[0066] F、春玉米识别模型构建
[0067] 根据识别特征及其对应的特征阈值,构建春玉米识别模型。其具体的模型构建步骤为:
[0068] (1)根据植被地物类型在春玉米抽雄期至乳熟期的NDVI平均值明显高于非植被,剔除非植被地物类型。本案例春玉米抽雄期至乳熟期为年日序第193–257天,该物候期对应的NDVI特征阈值为0.40,故其数学模型为:
[0069] Mean(NDVI193-257)>0.40
[0070] (2)根据春玉米拔节期前期农作物的NDVI平均值低于非农作物,剔除非农作物植被类型。本案例春玉米拔节期前期为年日序第129–153天,该物候期对应的NDVI特征阈值为0.33,故其数学模型为:
[0071] Mean(NDVI129-153)<0.33
[0072] (3)根据春玉米成熟期前期的NDVI平均值低于水稻,剔除部分水稻并提取玉米。本案例春玉米成熟期前期为年日序第257–265天,该物候期对应的NDVI特征阈值为0.57,故其数学模型为:
[0073] Mean(NDVI257-265)<0.57
[0074] (4)根据春玉米拔节期前期的红光反射率平均值高于水稻,进一步剔除部分水稻并提取玉米。本案例春玉米拔节期前期为年日序第129–153天,该物候期对应的红光反射率特征阈值为0.10,故其数学模型为:
[0075] Mean(RED129-153)>0.10
[0076] (5)根据春玉米拔节期前期的近红外反射率平均值高于水稻,进一步剔除部分水稻并提取玉米。本案例春玉米拔节期前期为年日序第129–153天,该物候期对应的近红外反射率特征阈值为0.18,故其数学模型为:
[0077] Mean(NIR129-153)>0.18
[0078] (6)根据春玉米成熟期前期的近红外反射率平均值低于水稻,进一步剔除部分水稻并提取玉米。本案例春玉米成熟期前期为年日序第257–265天,该物候期对应的近红外反射率特征阈值为0.30,故其数学模型为:
[0079] Mean(NIR257-265)<0.30
[0080] 上述数学模型中,Mean()表示遥感像元在时间上的求均值操作,NDVI、RED和NIR分别代表归一化差值植被指数、红光反射率和近红外反射率影像;下标数字代表年日序,即一年内的第几天。将上述6个数学模型得到的结果求交集之后的空间分布区即为识别的春玉米分布范围。
[0081] G、春玉米空间分布信息提取
[0082] 基于区域的遥感红光反射率、近红外反射率以及归一化差值植被指数时间序列数据,利用构建的春玉米识别模型,提取区域内春玉米的空间分布信息。
[0083] 本案例为辽宁省MODIS遥感数据的红光反射率、近红外反射率以及归一化差值植被指数时间序列数据。
[0084] 以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。