一种基于高分卫星遥感数据的耕地变化监测方法及系统与流程

本发明涉及农业信息技术领域，具体来说，涉及一种基于高分卫星遥感数据的耕地变化监测方法及系统。

背景技术：

近年来，由于我国经济发展、人口增长以及城市扩张等因素，大量耕地被占用，我国的耕地面积在持续减少，为了进一步保护耕地资源，减少耕地的非法侵占，国家制定了一系列关于耕地保护有关的政策和法规，以放缓耕地减少的速率，然而依靠传统的人工调查方法，无法对耕地变化情况进行实时准确的监测，耕地侵占情况时有发生，耕地保护效果并不乐观。遥感技术由于具有监测范围广、数据采集快速，节约成本等优点，现在已经被广泛应用于农作物信息快速提取中，为耕地面积变化监测提供了十分强大的技术支撑。因此本发明一种基于高分卫星遥感数据的耕地变化监测方法及系统，利用不同时相的高分遥感影像实现大范围、大尺度、短周期的耕地面积变化监测，为土地违规行为提供执法依据，也为土地合理利用、生态环境保护提供决策支持。

技术实现要素：

本发明提出了一种基于高分卫星遥感数据的耕地变化监测方法及系统。

第一方面，提供了一种基于高分卫星遥感数据的耕地变化监测方法，包括以下步骤：

步骤a，遥感数据获取以及遥感影像预处理，根据研究区农作物生长的物候特征选择适宜时相的高分遥感数据，对不同时相的高分卫星遥感影像进行影像预处理，分别包括对影像的辐射定标、大气校正和正射校正，得到预处理后的不同时相的全色和多光谱遥感影像；

步骤b，地面实测数据获取，采用实地调查的形式收集和存储野外调查样地数据，包括研究区主要土地类型的名称以及地理坐标信息，根据获取的土地类型实测数据结合人工目视解译选取样本数据和验证数据；

步骤c，提取不同时相的耕地信息，运用面向对象分类方法和机器学习算法，结合上一步选取的地面实测样本数据对不同时相遥感影像进行土地类型分类识别，重点提取研究区耕地信息,并根据地面实测验证数据进行精度检验；

步骤d，分析统计耕地信息的变化数据，根据提取的不同时相耕地信息，采用逻辑运算的方法分析不同时相耕地面积变化情况，将其分为耕地面积增加、耕地面积减少以及耕地面积未变化三种类型；

步骤e,生成耕地信息分布专题图，根据提取的不同时相耕地信息，得出耕地不断减少；根据统计的耕地信息变化数据，生成对应的耕地变化监测专题图，供用户参考；

进一步的，所述步骤c中的机器学习算法包括最大似然法、支持向量机、神经网络、随机森林、决策树以及k最近邻算法。

进一步的所述步骤d中，逻辑运算的算法过程为：将遥感影像识别的土地类型分为耕地与非耕地，其中耕地均赋值为1，其余非耕地赋值为0；将前后两个时相的遥感影像分类结果进行减法运算，得到-1、0和1三种结果，分别对应于-1非耕地变为耕地即耕地增加、0即耕地未变化、1耕地变为非耕地即耕地减少。

第二方面，提供了一种基于高分卫星遥感数据的耕地变化监测系统，包括以下模块：

第一模块，遥感数据获取及影像预处理模块，该模块由高分遥感卫星信号收集器和数据存储处理系统组成，用于收集和存储目标区域不同时相的高分遥感影像，并对不同时相的遥感影像进行预处理，包括辐射定标、大气校正、正射校正和影像融合，得到预处理后的遥感影像；

第二模块，地面实测数据获取模块，地面实测数据获取模块由地理信息数据采集器组成，用于获取野外调查样地数据，研究区主要土地类型的名称以及地理坐标信息，作为下一个耕地信息提取模块的样本数据和验证数据；

第三模块，耕地信息提取模块，用于接收影像预处理模块预处理后的遥感影像以及地面实测数据获取模块的样本数据和验证数据，运用面向对象分类方法和机器学习算法对不同时相遥感影像进行土地类型分类识别和精度检验，提取研究区不同时相耕地信息；

第四模块，耕地信息变化分析统计模块，用于接收耕地信息提取模块的不同时相耕地信息，采用逻辑运算分析不同时相耕地面积变化情况，统计耕地面积增加、耕地面积减少以及耕地面积不变三种类型的耕地信息变化数据；

第五模块,耕地信息分布显示模块，用于接收耕地信息变化分析统计模块得到的耕地信息变化结果是不断减少；

进一步，所述第三模块中的机器学习算法包括最大似然法、支持向量机、神经网络、随机森林、决策树以及k最近邻算法。

进一步，所述第四模块中，逻辑运算的算法过程为：将遥感影像识别的土地类型分为耕地与非耕地，其中耕地均赋值为1，其余非耕地赋值为0；将前后两个时相的遥感影像分类结果进行减法运算，得到-1、0和1三种结果，分别对应于-1非耕地变为耕地即耕地增加、0即耕地未变化、1耕地变为非耕地即耕地减少。

本发明一种基于高分卫星遥感数据的耕地变化监测方法及系统，利用不同时相的高分遥感影像，采用面向对象和机器学习算法对耕地信息进行大范围、大尺度的面积变化自动监测，有效提高了耕地信息变化的识别精度，为土地违规行为提供执法依据，也为土地合理利用、生态环境保护提供决策支持。

附图说明

图1为本发明实施例的方法技术路线图；

图2为本发明监测到的研究区耕地违规侵占影像图；

图3为本发明系统模块功能示意图。

具体实施方式

下面结合附图和具体实施实例对本发明做进一步说明。

本发明的实施实例是基于高分一号遥感数据，对2015年和2017年研究区的耕地变化情况进行监测，技术路线参照图1，具体步骤如下：

步骤a，研究区遥感数据获取以及遥感影像预处理，根据研究区农作物生长的物候特征选择2015年6月和2017年6月时相的高分一号遥感数据，对不同时相的高分一号遥感影像进行影像预处理，分别包括对影像的辐射定标、大气校正和正射校正，得到预处理后的不同时相的全色和多光谱遥感影像；

步骤b，研究区地面实测数据获取，采用实地调查的形式收集和存储野外调查样地数据，包括研究区主要土地类型的名称(耕地、林地、建筑用地、道路、水域)以及地理坐标信息，根据获取的土地类型实测数据结合人工目视解译选取样本数据和验证数据；

步骤c，提取研究区不同时相的耕地信息，运用面向对象法和最大似然分类法，结合上一步选取的地面实测样本数据对2015年6月和2017年6月时相的遥感影像进行土地类型分类识别，将研究区土地类型划分为耕地、林地、建筑用地、道路、水域和其它六大类，并根据地面实测验证数据进行精度检验，总体分类精度均在95％以上；接着重点提取研究区耕地信息,并将林地、建筑用地、道路、水域和其它五类统一划为非耕地，不同时相的研究区土地类型均被划分为耕地和非耕地两种类型，以便进行下一步的耕地信息变化分析；

步骤d，分析统计研究区耕地信息的变化数据，根据提取的2015年和2017年耕地信息，采用逻辑运算的方法分析2015年和2017年耕地面积变化情况，将其分为耕地面积增加、耕地面积减少以及耕地面积未变化三种类型。逻辑运算的算法过程为：将遥感影像识别的土地类型分为耕地与非耕地，其中耕地均赋值为1，其余非耕地赋值为0；将前后两个时相的遥感影像分类结果进行减法运算，得到-1、0和1三种结果，分别对应于-1非耕地变为耕地即耕地增加、0即耕地未变化、1耕地变为非耕地即耕地减少；

步骤e,得到耕地面积减少，可以为耕地违规侵占行为提供执法依据，如图2；

具体实施时，本发明所提供方法是可以通过内置于计算机、服务器内的软件、硬件或软硬件结合的系统以模块化的方式实现，本发明还提供了一种基于高分卫星遥感数据的耕地变化监测系统，系统模块功能示意图如图3，具体实施包括以下模块：

第一模块，遥感数据获取及影像预处理模块，该模块由高分遥感卫星信号收集器和数据存储处理系统组成，在本发明实施例中用于收集和存储研究区2015年6月和2017年6月时相的高分一号遥感影像，并对不同时相的遥感影像进行预处理，包括辐射定标、大气校正、正射校正和影像融合，得到预处理后的遥感影像；

第二模块，地面实测数据获取模块，地面实测数据获取模块由地理信息数据采集器组成，在本发明实施例中用于获取研究区野外调查样地数据，研究区主要土地类型的名称(耕地、林地、建筑用地、道路、水域)以及地理坐标信息，作为下一个耕地信息提取模块的样本数据和验证数据；

第三模块，耕地信息提取模块，在本发明实施例中用于接收影像预处理模块预处理后的高分一号遥感影像以及地面实测数据获取模块的研究区样本数据和验证数据，运用面向对象方法和最大似然分类法对2015年6月和2017年6月时相的遥感影像进行土地类型分类识别，提取耕地信息将研究区划分为耕地和非耕地两种类型，对分类结果进行精度检验，总体分类精度均在95％以上；

第四模块，耕地信息变化分析统计模块，在本发明实施例中用于接收耕地信息提取模块的2015年6月和2017年6月时相的耕地信息，采用逻辑运算分析不同时相耕地面积变化情况，统计耕地面积增加、耕地面积减少以及耕地面积不变三种类型的耕地信息变化数据。逻辑运算的算法过程为：将遥感影像识别的土地类型分为耕地与非耕地，其中耕地均赋值为1，其余非耕地赋值为0；将前后两个时相的遥感影像分类结果进行减法运算，得到-1、0和1三种结果，分别对应于-1非耕地变为耕地即耕地增加、0即耕地未变化、1耕地变为非耕地即耕地减少；

第五模块,耕地信息分布显示模块，在本发明实施例中用于接收耕地信息变化分析统计模块得到的研究区耕地信息变化结果并生成专题图呈现给用户不断减少，可以为耕地违规侵占行为提供执法依据，如图2。

本发明一种基于高分卫星遥感数据的耕地变化监测方法及系统，利用不同时相的高分遥感影像，采用面向对象和机器学习算法对耕地信息进行大范围、大尺度的面积变化自动监测，有效提高了耕地信息变化的识别精度，为土地违规行为提供执法依据，也为土地合理利用、生态环境保护提供决策支持。

以上所述，仅为本发明的具体实施方式，但本发明并不仅限于说明书和实施方式中所例运用，任何本领域的技术人员在不脱离本发明揭露的技术范围内，相对于上述实施例进行各种推演、修改和替换仍属于本发明所保护的范围，因此本发明的保护范围以权利要求书所界定的为准。