[发明专利]一种作物分布制图方法

说明书

本发明公开了一种作物分布制图，依据遥感数据提取的栅格尺度作物种植模式信息以及相应行政单元该种植模式上各类作物的面积比例，计算每个栅格上各类作物的面积，从而制作作物空间分布图。本发明将遥感识别的作物物候特征与作物实际物候结合，以确定作物空间分布状况，能够通过遥感宏观、快速获取、多尺度等特征，实现大范围、多尺度、且具有较高空间异质性的作物空间分布制图。

技术领域

本发明涉及数据分析领域，具体地说，是涉及一种作物分布制图。

背景技术

作物时空分布反映了不同时期人类在空间上利用土地从事农业生产的状况，它是开展农田生态系统格局与功能、陆地生态系统循环、全球变化等研究的重要基础数据。了解作物时空分布对保障国家粮食安全和资源环境可持续发展具有重要意义。

目前，作物时空分布信息获取途径主要包括行政统计和遥感监测。其中，统计方法仅能反映一定行政单元内的作物的数量变化，难以体现统计单元内部作物的空间差异性；并且，统计数据获取消耗大量人力、物力和财力，也受到人为因素干扰。单纯的遥感技术存在大气干扰、尺度转换和混合像元等因素的影响，难以获取大区域尺度全农作物系列的空间分布信息。因此，统计数据和遥感数据进行深度融合，开展多作物、甚至全作物类型的空间分布制图显得尤为必要。

发明内容

本发明的目的在于提供一种作物分布制图方法 ，

为实现上述目的，本发明采用的技术方案如下：

本发明包括以下步骤

A获取遥感数据信息，构建像元尺度的植被指数时间序列曲线，通过所述像元尺度时间序列曲线波峰数量的提取复种信息，依据波峰发生时间，对所述像元尺度反映的种植季节信息进行提取和归类细化，得到各个栅格作物种植模式信息；

B获取行政单元作物播种面积数据，将所述行政单元不同季相作物进行归并，求得各种植模式不同季相作物播种面积信息；

C依据遥感提取的栅格尺度作物种植模式及其对应的耕地面积分布以及所述不同种植模式上各季相作物的面积信息，计算每个栅格上各类作物的面积，从而制作作物空间分布图。

进一步地，所述特征参数包括时序曲线波峰波峰数和时序曲线波峰时间点，根据所述特征参数推定得到所述栅格种植模式信息：

具体地，所述预处理包括采用时间序列谐波分析法对数据进行去噪处理使时间序列曲线平滑。

进一步地，所述像元尺度提取的具体流程如下：

设像元i各时相像素值为a₁，a₂…a₂₃，

①比较两相邻像素之间的EVI的大小关系，并用Δa_i序列数据来表示：

②计算相邻Δa_i之差：

ΔΔa_i＝Δa_i+1-Δa_i (3-8)

其中，若ΔΔa_i0，则i-1为局部波峰所在时间点，a_i-1为局部波峰EVI值；若ΔΔa_i0，则i-1为局部波谷所在时间点，a_i-1为局部波谷EVI值。