面向东北冷凉区的作物轮作休耕遥感监测方法与装置与流程

本发明涉及农作物遥感监测领域，尤其是涉及一种面向东北冷凉区的轮作休耕遥感监测方法和装置。

背景技术：

农作物种植结构的调整将是中国农业生产可持续发展的国策，轮作休耕是实施″藏粮于地，藏粮于技″的战略布局，对保障我国粮食安全和农业绿色转型具有深远意义。轮作(croprotation)是指在同一田块上不同年度间有顺序地轮换种植不同作物或以复种方式进行的种植方式。休耕(fallow)是指耕地在可种作物的季节只耕不种或不耕不种的方式。轮作休耕按区域分为东北冷凉区、地下水漏斗区、重金属污染区、生态严重退化区等4大区域。目前，东北冷凉区作物种植面积的获取多基于统计数据，其时效性、真实性较差，且无法呈现空间分异特征。利用卫星遥感技术开展耕地轮作监测，具有低成本、便利化、高效率、系统性等明显优势，降低人为因素的干扰，实现对轮作面积的精确监测，实现公平公正。相比于已经在轮作休耕上取得较好成果的国家，中国还处于初级探索阶段，目前主要多从不同轮作复种模式比较、轮作对土壤环境的影响等方面进行研究，对于轮作的遥感核查方面研究较少。另外，大范围的东北地区主要作物的提取多以中低分辨率为主，无法满足精细化的作物提取，难以实现精准的轮休与休耕面积估算。

东北冷凉区轮作与休耕模式并存，轮作作物复杂多样，轮作模式主要是″一主多辅”即玉米与大豆轮作为主、玉米与薯类、籽粒玉米与青贮玉米、苜蓿、黑麦草、谷子、花生、向日葵、等油料作物高梁等轮作为轮作辅。休耕模式是深耕与种植绿肥、种植苜蓿或油菜等肥田养地作物。轮作与休耕试点周期为三年，实施年际间轮作休耕，同一地块可以动态调整，使得东北冷凉区的轮作休耕监测更加复杂。因此，为满足东北冷凉区农作物轮作遥感核查业务需求，亟待需要一种新的遥感监测方法。

技术实现要素：

本发明要解决的技术问题是：为解决现有技术中的不足，从而提供一种面向东北冷凉区的轮作休耕遥感监测方法和装置。

本发明解决其技术问题所采用的技术方案是：

一种面向东北冷凉区的轮作休耕遥感监测方法，包括以下步骤：

s1：获取第n年与第n-1年的监测地区的遥感影像；所述遥感影像包括蓝光、绿光、红光、近红外、3个红边波段、短波红外波段，遥感影像的时间包含5月-9月的影像；

s2：分别识别第n年与第n-1年的遥感影像上各像元的作物类型或者是否为休耕地，并按照相同的作物类型提取仅种植单一作物的地块；

s3：若第n年与第n-1年的遥感影像上同一地块的的识别结果分别为两种不同的作物，则说明进行了轮作；

若第n年与第n-1年的遥感影像上同一地块的识别结果分别为休耕地和作物，则说明进行了休耕。

优选地，本发明的面向东北冷凉区的轮作休耕遥感监测方法，s2步骤中，

遥感影像上的像元符合条件1时识别为玉米，条件1：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08＞t4；

遥感影像上的像元符合条件2时识别为水稻，条件2：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08≤t4，且ndri08≤t5；

遥感影像上的像元符合条件3时识别为甜菜，条件3：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08≤t4，且ndri08＞t5，且gi08≤t6；

遥感影像上的像元符合条件4时识别为大豆，条件4：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08≤t4，且ndvi08＞t5，且gi08＞t6；

遥感影像上的像元识别符合条件5时识别为休耕地，条件5：t1＜ndvi05＜t2，且ndvi08＜t7，且ndbi08＞t8；

其中，ndvi05表示5月份影像的ndvi；

ndvi08表示8月份影像ndvi值；

rep08表示8月份影像的红边位置指数；

ndri08表示8月份影像的归一化水稻指数ndri；

gi08表示8月份影像的绿度指数；

ndbi08表示8月份影像的归一化建筑指数；

阈值t1-t8均由样本数据统计所得。

优选地，本发明的面向东北冷凉区的轮作休耕遥感监测方法，

其中，ρred、ρnir、ρswir分别为红光波段、近红外、短波红外波段像元反射率值；

其中，ρ为拐点反射率，ρ＝(ρ670+ρ780)/2，ρ670、ρ700、ρ740、ρ780分别为在波长670nm，700nm，740nm，780nm处的反射率；

gi＝10000*(ρblue-ρgreen)；

其中，ρred、ρswir、ρblue、ρgreen分别为红光波段、短波红外、蓝光波段、绿光坡段的反射率值。

优选地，本发明的面向东北冷凉区的轮作休耕遥感监测方法，s2步骤中，对地块的提取包括以下步骤：

s201：对遥感影像上的地块进行识别；

s202：对地块内的作物进行识别，若一个地块内含有至少两种作物类型，则标记所述地块为待校正地块；

s203：将包含待校正地块的遥感影像使用高斯滤波器平滑图像；

s204：用一阶微分算子的有限差分计算滤波后图像梯度的幅值和方向；

s205：对梯度幅值进行非极大值抑制，遍历图像，若某个像素的梯度值小于沿着梯度方向上的两个相邻像素的梯度值，则判定此像素点不是边缘点，将其灰度值置为0；

s206：用双阈值算法检测和连接边缘，使用累计直方图计算两个阈值，分别为高阈值和低阈值，分别记thigh和tlow，凡是大于thigh的则判定为边缘，凡是小于tlow的判断为不是边缘；如果检测结果大于等于tlow但又小于等于thigh，则继续判断该像素的邻接像素中有没有超过高阈值的像素点，如果有，则该像素点就是边缘点，否则就不是边缘；

s207：如果检测结果中作为边缘点的像素点不连续而形成断点，则依据边缘局部方向不变性以n×n模板搜索端点；连接搜索到的任意两个端点，评价该端点连线是否为线性，若结果为线性则以连线作为边缘，如果无法搜索端点或者端点不为线性，则由人工进行标记；

s208：将提取出的线性边界线切割耕地地块轮廓边界最终形成耕地地块单元，从而获得仅存在一种作物的地块单元。

优选地，本发明的面向东北冷凉区的轮作休耕遥感监测方法，s3步骤中还按照第n年与第n-1年的遥感影像上同一地块的识别结果不同的变化情况进行统计，分析轮作模式、轮作面积与休耕面积。

本发明还提供一种面向东北冷凉区的轮作休耕遥感监测装置，包括：

影像获取模块：用于获取第n年与第n-1年的监测地区的遥感影像；所述遥感影像包括蓝光、绿光、红光、近红外、3个红边波段、短波红外波段，遥感影像的时间包含5月-9月的影像；

地块类型识别模块：用于分别识别第n年与第n-1年的遥感影像上各像元的作物类型或者是否为休耕地，并按照相同的作物类型提取仅种植单一作物的地块；

结果判断模块：用于判断地块是否为轮作或者休耕，若第n年与第n-1年的遥感影像上同一地块的的识别结果分别为两种不同的作物，则说明进行了轮作；若第n年与第n-1年的遥感影像上同一地块的识别结果分别为休耕地和作物，则说明进行了休耕。

优选地，本发明的面向东北冷凉区的轮作休耕遥感监测装置，地块类型识别模块中，

遥感影像上的像元符合条件1时识别为玉米，条件1：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08＞t4；

遥感影像上的像元符合条件2时识别为水稻，条件2：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08≤t4，且ndri08≤t5；

遥感影像上的像元符合条件3时识别为甜菜，条件3：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08≤t4，且ndri08＞t5，且gi08≤t6；

遥感影像上的像元符合条件4时识别为大豆，条件4：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08≤t4，且ndvi08＞t5，且gi08＞t6;

遥感影像上的像元识别符合条件5时识别为休耕地，条件5：t1＜ndvi05＜t2，且ndvi08＜t7，且ndbi08＞t8；

其中，ndvi05表示5月份影像的ndvi；

ndvi08表示8月份影像ndvi值；

rep08表示8月份影像的红边位置指数；

ndri08表示8月份影像的归一化水稻指数ndri；

gi08表示8月份影像的绿度指数；

ndbi08表示8月份影像的归一化建筑指数；

阈值t1-t8均由样本数据统计所得。

优选地，本发明的面向东北冷凉区的轮作休耕遥感监测装置，

其中，ρred、ρnir、ρswir分别为红光波段、近红外、短波红外波段像元反射率值；

其中，ρ为拐点反射率，ρ＝(ρ670+ρ780)/2，ρ670、ρ700、ρ740、ρ780分别为在波长670nm，700nm，740nm，780nm处的反射率；

gi＝10000*(ρblue-ρgreen)；

其中，ρred、ρswir、ρblue、ρgreen分别为红光波段、短波红外、蓝光波段、绿光坡段的反射率值。

优选地，本发明的面向东北冷凉区的轮作休耕遥感监测装置，地块类型识别模块中，对地块的提取包括以下步骤：

s201：对遥感影像上的地块进行识别；

s202：对地块内的作物进行识别，若一个地块内含有至少两种作物类型，则标记所述地块为待校正地块；

s203：将包含待校正地块的遥感影像使用高斯滤波器平滑图像；

s204：用一阶微分算子的有限差分计算滤波后图像梯度的幅值和方向；

s205：对梯度幅值进行非极大值抑制，遍历图像，若某个像素的梯度值小于沿着梯度方向上的两个相邻像素的梯度值，则判定此像素点不是边缘点，将其灰度值置为0；

s206：用双阈值算法检测和连接边缘，使用累计直方图计算两个阈值，分别为高阈值和低阈值，分别记thigh和tlow，凡是大于thigh的则判定为边缘，凡是小于tlow的判断为不是边缘；如果检测结果大于等于tlow但又小于等于thigh，则继续判断该像素的邻接像素中有没有超过高阈值的像素点，如果有，则该像素点就是边缘点，否则就不是边缘；

s207：如果检测结果中作为边缘点的像素点不连续而形成断点，则依据边缘局部方向不变性以n×n模板搜索端点；连接搜索到的任意两个端点，评价该端点连线是否为线性，若结果为线性则以连线作为边缘，如果无法搜索端点或者端点不为线性，则由人工进行标记；

s208：将提取出的线性边界线切割耕地地块轮廓边界最终形成耕地地块单元，从而获得仅存在一种作物的地块单元。

优选地，本发明的面向东北冷凉区的轮作休耕遥感监测装置，结果判断模块中还按照第n年与第n-1年的遥感影像上同一地块的识别结果不同的变化情况进行统计，分析轮作模式、轮作面积与休耕面积。

本发明的有益效果是：

本申请提供面向东北冷凉区的轮作休耕遥感监测方法和装置，利用归一化植被指数用以区分植被与非植被，归一化建筑指数用以区分建筑和休耕地，红边位置指数用来提取玉米，归一化水稻指数用来提取水稻，绿度指数用以区分大豆和甜菜。具有识别准确度高，能够准确识别出作物，最后根据两年的作物类型进行对比，可以有效地分辨出是否进行了轮作休耕。

附图说明

下面结合附图和实施例对本申请的技术方案进一步说明。

图1是轮作休耕总体技术流程图；

图2是耕地地块提取与轮作休耕模式监测技术流程图；

图3是耕地地块识别结果图；

图4是基于多时相、多指数的作物提取技术

图5是本年度与上一年度作物提取结果；

图6是耕地地块修正示意图；

图7是轮作修耕监测结果图；

具体实施方式

需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。

下面将参考附图并结合实施例来详细说明本申请的技术方案。

实施例1

本实施例提供一种面向东北冷凉区的轮作休耕遥感监测方法，如图1所示，包括以下步骤：

s1：获取第n年与第n-1年的监测地区的遥感影像；

所述遥感影像包括蓝光、绿光、红光、近红外、3个红边波段、短波红外波段，如哨兵2号影像，对遥感影像做波段合成、大气校正、几何校正、裁剪、镶嵌等处理，遥感影像包含5月-9月的影像，覆盖研究区大豆、水稻、玉米等作物的整个生育期；

s2：分别识别第n年与第n-1年的遥感影像上各像元的作物类型或者是否为休耕地，并按照相同的作物类型提取仅种植单一作物的地块，并按照相同的作物类型提取仅种植单一作物的地块；

遥感影像上的像元符合条件1时识别为玉米，条件1：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08＞t4；

遥感影像上的像元符合条件2时识别为水稻，条件2：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08≤t4，且ndri08≤t5；

遥感影像上的像元符合条件3时识别为甜菜，条件3：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08≤t4，且ndri08＞t5，且gi08≤t6；

遥感影像上的像元符合条件4时识别为大豆，条件4：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08≤t4，且ndvi08＞t5，且gi08＞t6；

遥感影像上的像元识别符合条件5时识别为休耕地，条件5：t1＜ndvi05＜t2，且ndvi08＜t7，且ndbi08＞t8；

其中，ndvi05表示5月份影像的ndvi，阈值t1与t2由样本数据统计所得，如t1＝0，t2＝0.5；

ndvi08表示8月份影像ndvi值，阈值t3与t7由样本数据统计所得，例如t3＝0.6，t7＝0.2；

rep08表示8月份影像的红边位置指数，阈值t4由样点数据统计而得，例如t4＝728；

ndri08表示8月份影像的归一化水稻指数ndri，阈值t5由样点数据统计而得，例如t5＝0.5；

gi08表示8月份影像的绿度指数，阈值t6由样点数据统计而得，t6＝52；

ndbi08表示8月份影像的归一化建筑指数，阈值t8由样点数据统计而得，例如t8＝-0.1，

样本数据统计阈值是指在影像上创建随机点，目视判别其类型，分别提取各波段的反射率，计算作物分类的阈值。

其具体识别过程如下：

s21：利用第n-1年或者第n年的5月份影像的ndvi剔除水体和林地，判断ndvi05的取值范围是否符合t1＜ndvi05＜t2，若条件为否则像元标识为水体或林地，若条件为真则进行s22判断；

s22：然后利用8月份影像ndvi剔除建筑或其他人造地物。判断像元ndvi08是否大于t3，若为否则像元标识为建筑或其他人造地物，若为是则进行s23判断；

可进一步判断ndvi08是否小于t7，其中ndvi08表示8月份影像ndvi值，阈值t7由样本数据统计所得，例如t7＝0.2，若ndvi08＞t7则像元标识为浓密植被，若ndvi08＜t7，则设定条件ndbi08＞t8，其中，ndbi08表示8月份影像的归一化建筑指数，阈值t8由样点数据统计而得，例如t8＝-0.1，若条件为否则像元则像元标识为建筑，若条件为真则像元标识为休耕地。

s23：8月份玉米的红边位置指数比大豆、水稻等作物的大，rep利于提取玉米。判断像元的rep08是否小于t4，若条件为否则像元标识为玉米，若条件为真则进行s4判断；

s24：短波红外光谱反射率随土壤含水量的降低而升高，其对植被水分变化敏感。基于短波红外与红波段构建归一化水稻指数ndri，定义为短波红外波段的反射值与红光波段的反射值之差与两者之和的比值。在步骤s3后判断像元的ndri08是否大于t5，若条件为否则像元标识为水稻，若条件为真则进行s25判断；

s25：利用颜色不同区分大豆甜菜，8月份大豆和甜菜表现出两种不同的颜色，大豆墨绿色，甜菜为亮绿色。基于蓝光与绿光波段构建绿度指数gi，定义为蓝光与绿光波段的反射值之差。在步骤s24后判断像元的gi08＞t6，若条件为否则像元标识为甜菜，若条件为真则像元标识为大豆；

本步骤中，各种指数的定义如下：

归一化植被指数(ndvi)用以区分植被与非植被；归一化建筑指数(ndbi)用以区分建筑和休耕地，红边位置指数(rep)用来提取玉米；归一化水稻指数ndri用来提取水稻；绿度指数gi用以区分大豆和甜菜。

ndvi定义为近红外波段与红波段两者反射率之差与之和比值，ndbi定义为短波红外波段的反射率与近红外波段的反射率之差与两者之和的比值，计算公式定义如下：

其中，ρred、ρnir、ρswir分别为红光波段、近红外、短波红外波段像元反射率值。

红边位置指数(rep)，计算公式为

ρ＝(ρ670+ρ780)/2

其中，ρ为拐点反射率，ρ670、ρ700、ρ740、ρ780分别为在波长670nm，700nm，740nm，780nm处的反射率。700与740是700nm～740nm区间进行内插产生的常数。

除了以上常用的植被指数外，本发明定义了两种新的植被指数分别为归一化水稻指数ndri和绿度指数gi。ndri定义为短波红外波段的反射值与红光波段的反射值之差与两者之和的比值，计算如下：

绿度指数gi，计算如下：

gi＝10000*(ρblue-ρgreen)

其中，ρred、ρswir、ρblue、ρgreen分别为红光波段、短波红外、蓝光波段、绿光坡段的反射率值。

轮作与休耕在地块基础上监测，地块是基于深度学习获取的地块边缘矢量结果。耕地地块提取结果未考虑地面种植情况，往往获取的地块上存在多种作物种植，为满足轮作休耕精准监测需要与进一步提高地块提取结果的准确性，需要针对异常耕地地块修正。异常耕地地块一般表现为地块内含有至少两种作物类型，作物由于颜色不同具有明显的界限，这种界限宽度较小，通常可视作直线。异常耕地地块的修正就是检测出这种直线并基于此直线切割耕地地块。切割后的耕地地块内仅包含一种作物。

本实施例中东北冷凉区主要种植大豆、水稻、玉米，在8月份玉米的红边位置达到最大，此时设定红边位置的阈值，便于玉米的识别；

现有技术中的水稻提取大多是基于ndvi或是ndwi，利用水稻的移栽期和抽穗期之间ndvi变化差异实现的。这类方法需要多景影像参与，增加了计算量。另外，这种方法往往将大豆错分为水稻。短波红外在植被水分、矿物类型识别上具有较为明显的反射谱特征。8月份影像近红外、短波红外、红波段组成的假彩色影像具有明显特征，表现为水稻为玫瑰红色，大豆和玉米呈现黄色。ndri是短波红外和红波段的反射率之和与反射率之差的比值，充分利用了短波红外对土壤含水量敏感的特征，该指数更加利于提取水稻。大豆和甜菜在遥感影像上表现形式不同，基于蓝光与绿光波段的反射值之差构建的绿度指数，用以区分大豆和甜菜。

本实施例组合ndvi、rep、ndri构建的决策树模型，利用较少的影像数据，充分利用作物的ndvi特征、红边特征以及短波红外水稻的特征，减少了数据量，提取方便、精度较准。

s3：若第n年与第n-1年的遥感影像上同一地块的的识别结果分别为两种不同的作物，则说明进行了轮作；

若第n年与第n-1年的遥感影像上同一地块的识别结果分别为休耕地和作物，则说明进行了休耕。

进一步地，还可以统计玉米变为大豆、大豆变为玉米、玉米变为甜菜、玉米变为水稻、水稻变为休耕地、休耕地变为水稻等情况，根据以上变化情况分析轮作和休耕模式，标记出轮作地块与轮作类型以及休耕地块，统计分析轮作模式、轮作面积与休耕面积。

步骤s2中，可以采用以下地块提取方法：

步骤一：首先将遥感影像的真彩色三通道与全色波段进行融合，得到分辨率为0.8米，数据格式为tiff，r-g-b三通道存储的影像。在原影像利用arcgis软件手工画出耕地地块的轮廓边界用以制作标签数据，基于轮廓边界的矢量数据创建与原数据具有相同投影信息单波段存储的标注影像，影像内容为0、1，分别表示非耕地(背景类)、耕地地块。将原影像与标注影像同时进行随机裁剪，裁剪大小为256×256像素，并对裁剪的图像进行数据增强处理，经随机裁剪和数据增强后共获得10000张256×256像素的.png格式的图片，其中25％即2500张作为验证样本，7500张为训练样本。

步骤二：利用tensorflow和keras深度学习框架搭建的u-net网络模型，进行模型训练。在训练之前设置学习率参数为0.001，训练批次大小为32，验证批次大小为16，训练次数为2000次，激活函数为sigmod函数。

步骤三：使用训练出的设施农业提取模型预测设施农业空间分布与类型，将待分类影像裁剪为固定大小(256×256)像素的图像块分别进行预测，采用滑动窗口像素为32的方式获取有一定重叠区域的图像块，然后对每一张预测的图像块保留中间某一区域的分类结果，舍弃边缘不准确的结果，再依次拼接。

将待分类影像根据影像大小和计算机内存进行分块裁剪，将分块的影像切割成大小为256×256像素的小图，切割后的数据被输入到训练好的模型中进行分类，然后将分类的小图拼接。对于分类结果进行众数滤波、腐蚀、膨胀、开闭运算等分类后处理操作。众数滤波剔除个数较少的像元，保留像元数量较多的目标。腐蚀操作腐蚀可以使目标区域范围″变小″，其实质造成图像中耕地地块的边界收缩，可以用来消除小且无意义的目标物。

步骤四：将耕地地块结果转为矢量。深度学习结果获得是栅格数据，利用arcgls软件栅格转矢量操作，分别经过二值化、细化等转为矢量，然后经平滑面操作，获取最初的耕地地块边界。耕地地块提取结果如图3所示。研究区内耕地地块共计6196块，面积共计42.52km²，最小地块面积为0.45亩，最大耕地地块面积为455亩。

步骤s2中，可以采用以下地块提取方法：

选择关键期影像的灰度图像或纹理图像。灰度图像是单波段影像，纹理图像是利用灰度共生矩阵计算的均值、方差、熵值等。

s201：对遥感影像上的地块进行识别；

s202：对地块内的作物进行识别，若一个地块内含有至少两种作物类型，则标记所述地块为待校正地块；

s203：将包含待校正地块的遥感影像使用高斯滤波器平滑图像；

s204：用一阶微分算子的有限差分计算滤波后图像梯度的幅值和方向；

s205：对梯度幅值进行非极大值抑制，遍历图像，若某个像素的梯度值小于沿着梯度方向上的两个相邻像素的梯度值，则判定此像素点不是边缘点，将其灰度值置为0；

s206：用双阈值算法检测和连接边缘，使用累计直方图计算两个阈值，分别为高阈值和低阈值，分别记t_high和t_low，凡是大于t_high的则判定为边缘，凡是小于t_low的判断为不是边缘；如果检测结果大于等于t_low但又小于等于t_high，则继续判断该像素的邻接像素中有没有超过高阈值的像素点，如果有，则该像素点就是边缘点，否则就不是边缘；

尽管canny能够提取出部分边缘信息，通过分析canny边缘检测的结果，发现在提取后的边缘中间存在单像素间断现象。断线边缘无法到达耕地地块的边界，即无法分割耕地地块。

s207：如果检测结果中作为边缘点的像素点不连续而形成断点，则依据边缘局部方向不变性以n×n(3×3像素)模板搜索端点；连接搜索到的任意两个端点，评价该端点连线是否为线性，若结果为线性则以连线作为边缘，如果无法搜索端点或者端点不为线性，则由人工进行标记；一种常用的线性评价方法是长宽比指数即线性兴趣区的最小外接矩形的长宽比值；

s208：将提取出的线性边界线切割耕地地块轮廓边界最终形成耕地地块单元，从而获得仅存在一种作物的地块单元。

哨兵2号是截止到申请日时唯一一个在红边范围含有三个波段的数据，其分辨率为10米，利于作物类型的监测，提供了一种新的遥感监测数据。高分辨遥感影富含丰富的纹理特征，利用高分辨率遥感影像结合深度学习技术能够确定地块边界，为轮作核查提供准确地块矢量。本发明提出了一种基于高分辨率遥感影像与哨兵影像的东北冷凉区作物轮作休耕监测技术。首先，基于高分辨遥感影像利用深度学习技术获取地块耕地边界，其次，基于哨兵影像结合作物物候特征、红边参数位置、归一化植被指数以及归一化水稻指数的提取大豆、水稻、玉米、小麦、甜菜等作物，将作物地块结果修正耕地地块，分析年际间作物变化情况，得出轮作休耕作物与面积。

实施例2

本实施例中还提供一种面向东北冷凉区的轮作休耕遥感监测装置，包括：

影像获取模块：用于获取第n年与第n-1年的监测地区的遥感影像；所述遥感影像包括蓝光、绿光、红光、近红外、3个红边波段、短波红外波段，遥感影像的时间包含5月-9月的影像；

地块类型识别模块：用于分别识别第n年与第n-1年的遥感影像上各像元的作物类型或者是否为休耕地，并按照相同的作物类型提取仅种植单一作物的地块；

结果判断模块：用于判断地块是否为轮作或者休耕，若第n年与第n-1年的遥感影像上同一地块的的识别结果分别为两种不同的作物，则说明进行了轮作；若第n年与第n-1年的遥感影像上同一地块的识别结果分别为休耕地和作物，则说明进行了休耕。

进一步地，地块类型识别模块中，

遥感影像上的像元符合条件1时识别为玉米，条件1：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08＞t4；

遥感影像上的像元符合条件2时识别为水稻，条件2：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08≤t4，且ndri08≤t5；

遥感影像上的像元符合条件3时识别为甜菜，条件3：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08≤t4，且ndri08＞t5，且gi08≤t6；

遥感影像上的像元符合条件4时识别为大豆，条件4：t1＜ndvi05＜t2，且ndvi08＞t3，且rep08≤t4，且ndvi08＞t5，且gi08＞t6；

遥感影像上的像元识别符合条件5时识别为休耕地，条件5：t1＜ndvi05＜t2，且ndvi08＜t7，且ndbi08＞t8；

其中，ndvi05表示5月份影像的ndvi；

ndvi08表示8月份影像ndvi值；

rep08表示8月份影像的红边位置指数；

ndri08表示8月份影像的归一化水稻指数ndri；

gi08表示8月份影像的绿度指数；

ndbi08表示8月份影像的归一化建筑指数；

阈值t1-t8均由样本数据统计所得。

优选地，本发明的面向东北冷凉区的轮作休耕遥感监测装置，

其中，ρred、ρnir、ρswir分别为红光波段、近红外、短波红外波段像元反射率值；

其中，ρ为拐点反射率，ρ＝(ρ670+ρ780)/2，ρ670、ρ700、ρ740、ρ780分别为在波长670nm，700nm，740nm，780nm处的反射率；

gi＝10000*(ρblue-ρgreen)；

其中，ρred、ρswir、ρblue、ρgreen分别为红光波段、短波红外、蓝光波段、绿光坡段的反射率值。

进一步地，地块类型识别模块中，对地块的提取包括以下步骤：

s201：对遥感影像上的地块进行识别；

s202：对地块内的作物进行识别，若一个地块内含有至少两种作物类型，则标记所述地块为待校正地块；

s203：将包含待校正地块的遥感影像使用高斯滤波器平滑图像；

s204：用一阶微分算子的有限差分计算滤波后图像梯度的幅值和方向；

s205：对梯度幅值进行非极大值抑制，遍历图像，若某个像素的梯度值小于沿着梯度方向上的两个相邻像素的梯度值，则判定此像素点不是边缘点，将其灰度值置为0；

s206：用双阈值算法检测和连接边缘，使用累计直方图计算两个阈值，分别为高阈值和低阈值，分别记thigh和tlow，凡是大于thigh的则判定为边缘，凡是小于tlow的判断为不是边缘；如果检测结果大于等于tlow但又小于等于thigh，则继续判断该像素的邻接像素中有没有超过高阈值的像素点，如果有，则该像素点就是边缘点，否则就不是边缘；

s207：如果检测结果中作为边缘点的像素点不连续而形成断点，则依据边缘局部方向不变性以n×n模板搜索端点；连接搜索到的任意两个端点，评价该端点连线是否为线性，若结果为线性则以连线作为边缘，如果无法搜索端点或者端点不为线性，则由人工进行标记；一种常用的线性评价方法是长宽比指数即线性兴趣区的最小外接矩形的长宽比值；

s208：将提取出的线性边界线切割耕地地块轮廓边界最终形成耕地地块单元，从而获得仅存在一种作物的地块单元。

进一步地，结果判断模块中还按照第n年与第n-1年的遥感影像上同一地块的识别结果不同的变化情况进行统计，分析轮作模式、轮作面积与休耕面积。

以上述依据本申请的理想实施例为启示，通过上述的说明内容，相关工作人员完全可以在不偏离本项申请技术思想的范围内，进行多样的变更以及修改。本项申请的技术性范围并不局限于说明书上的内容，必须要根据权利要求范围来确定其技术性范围。

本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。