本发明专利技术公开了一种土地整治后农田作物生长缺陷区遥感快速诊断方法，其包括：1、获取土地整治后农作物生长季的高分辨率多光谱影像；2、预处理相应的遥感影像，并结合工程施工图裁剪得到整治项目区内所有田块的遥感影像；3、计算各个田块内每个像元的归一化植被指数NDVI；4、进行空间聚类分析；5、分离低值聚类、低值异常两种类型的空间范围，并将其作为土地整治后农作物生长缺陷区。该方法可以快速便捷地实现地整治后农田作物生长缺陷区的检测，为完善土地整治的规划方案和精细后期养护管理提供依据。
【技术实现步骤摘要】
一种土地整治后农田作物生长缺陷区遥感快速诊断方法
本专利技术涉及基于遥感的快速检测领域，具体为一种土地整治后农田作物生长缺陷区遥感快速诊断方法。
技术介绍
为提高农田生产稳定性和提高粮食生产能力，保障国家粮食安全，我国正在大规模开展土地整治工程。根据全国土地整治规划(2016-2020年)，我国将建设4-6亿亩高标准农田，新补充耕地2000万亩。土地整治是一项工程密集、对生态环境扰动大的特殊活动。经过土地平整、土壤改良、灌排沟渠建设后，一些地块耕地地力得到提升，作物长势较好，但部分地块可能因为整治工程扰动后土体构型破坏、填埋压不实沉降、肥力短期难恢复、排水条件不良等原因，农作物长势较差，呈现为农作物生长的缺陷区。这些缺陷区需要后期追加和补充相应的工程措施，从而保证土地整治项目区农田产能得到全面提升。因此，识别土地整治后农田作物生长缺陷区十分必要。一般而言，土地整治要求监理或验收单位对土地整治工程质量、耕地地力进行检测。现行的检测方法主要是采用现场查验的方法，即在土地整治后现场布点采集作物长势、产量数据，然后通过比较分析来识别土地整治的有效性和耕地产能。近年也有一些单位尝试开发流动性检测车，例如专利一种适用于土地整治工程的流动检测车，专利公开号：CN203666474U；现场调查系统，例如一种土地整治工程实地调查系统，专利公开号：CN104731924A；来实现土地整治各项工程的检测。从现有的技术来看，还存在以下缺点：(1)效率低，通过现场查验、取样调查、统计分析需要大量的时间，在技术上很难保障农作物生长缺陷区的及时发现和再整治；(2)成本高，田间实地查验时需要耗费大量的人力、财力，为土地整治工程和后期管护增加了额外的成本；(3)缺乏空间可视化，在土地整治后，需要识别农田作物生长缺陷区的空间范围和位置，制备相关图件，供二次土地整治参考，但传统方法难以提供有空间显示功能的图件。
技术实现思路
本专利技术旨在提供一种基于遥感和空间统计的土地整治后农田作物生长缺陷区遥感快速诊断方法，提供准确、可靠、可利用性较好的检测成果。本专利技术解决上述技术问题采用的技术方案是：一种土地整治后农田作物生长缺陷区遥感快速诊断方法，包括以下步骤：步骤1)在土地整治后的农作物生长季内，获取覆盖土地整治区的基于卫星、航空摄影或无人机的高分辨率多光谱遥感影像；步骤2)预处理所述遥感影像，包括对遥感影像进行大气校正和几何配准，将土地整治项目区的工程施工图与预处理后的遥感影像叠加在一起，然后利用工程施工图中的所有田块边界对预处理后的遥感影像进行裁剪，得到土地整治项目区内所有田块的遥感影像；步骤3)计算所有田块内每个像元的归一化植被指数NDVI，NDVI是近红外与红色通道反射比值的一种变化形式，是农作物生长状况的最佳指示因子，计算公式为：NDVI＝(NIR-R)/(NIR+R)式中，NIR代表遥感影像近红外波段的反射率，R代表红光波段的反射率，NDVI取值范围在-1到1之间，NDVI值越高，表明农作物生长状况越好；在ENVI即TheEnvironmentforVisualizingImages软件中采用bandmath功能实现对各个田块每个像元的NDVI值；步骤4)进行空间聚类分析，计算各个像元NDVI值空间关联的LocalMoran’sI指数、各个像元NDVI值与所有田块NDVI平均值的差值Di、各个像元NDVI值标准化统计量Zi，计算公式分别为：式中，Xi、Xj分别为第i、j个像元的NDVI值，为所有田块NDVI的均值，n为像元总数，Wij为像元i和j之间的空间权重，VAR(Ii)是Ii的方差；当Zi取值在小于-1.96或大于+1.96时，该像元的空间聚类在置信度为95％的水平下没有显著性，属于“不显著”类型；当Zi取值在-1.96和+1.96区间时，该像元的空间聚类在置信度为95％的水平下有显著性，对于空间聚类有显著性的像元，Ii和Di均大于0时，该像元属于高值聚类HH，此时该像元包含同样高的NDVI的邻近像元；Ii大于0且Di小于0时，该像元属于低值聚类LL，此时该像元包含同样低的NDVI的邻近像元；Ii和Di均小于0时，该像元属于低值被高值围绕的异常值LH，此时该像元的邻近像元的NDVI值都较高；Ii小于0且Di大于0，该像元属于高值被低值围绕的异常值HL，此时该像元的邻近像元的NDVI值都较低；在上述5种类型中，低值聚类LL、低值被高值围绕的异常值LH区域长势显著低于其他地区，是作物生长的缺陷区；步骤5、分离低值聚类LL、低值被高值围绕的异常值LH两种类型的空间范围，并将其作为土地整治后农田作物生长缺陷区，将其单独显示在地图上，供二次土地整治工程施工时查阅、定位。本专利技术与现有技术相比的优点在于：(1)本专利技术可实现土地整治后农田作物生长缺陷区的快速检测。本专利技术基于遥感数据和空间分析方法，全过程基于计算机处理，可以快速地、便捷地完成检测工作。(2)本专利技术可实现空间显示功能。土地整治的规划、设计、施工需要明确缺陷区的位置、范围和面积。本专利技术最终提供一个地图来显示缺陷区，具有较高的实用价值。附图说明图1是土地整治后农田作物生长缺陷区遥感快速检测的技术流程图。图2是本专利技术一个实施例中覆盖土地整治区的原始高分遥感影像。图3是本专利技术一个实施例中土地整治区所有田块的高分遥感影像。图4是本专利技术一个实施例中土地整治区所有田块的NDVI空间分布。图5是本专利技术一个实施例中土地整治区所有田块的NDVI空间聚类模式。图6是本专利技术一个实施例中土地整治后农田作物生长缺陷区的检测结果。具体实施方式实施例1以江苏省金坛市某地基本农田土地整治为例，对本专利技术实施作进一步详细说明。该区土地整治工程完成于2017年5月，随后种植水稻。以下实例用于说明本专利技术，但不用来限制本专利技术的范围。如图1，本专利技术提供了一种土地整治后农田作物生长缺陷区遥感快速诊断方法，具体实施步骤如下：步骤1)确定检测区域，在土地整治后第一季水稻生长季(抽穗期)获取1期基于Geo-eye-01卫星的高分遥感影像，多光谱空间分辨率为1.65米，有蓝、绿、红、近红外四个波段，遥感数据获取时间为2017年8月23日，原始遥感影像如图2所示。步骤2)预处理所述遥感影像，包括：对Geo-eye-01遥感影像进行大气校正、几何配准、研究区裁剪的预处理操作。然后将土地整治项目区的工程施工图与预处理后的遥感影像叠加在一起，利用工程施工图中所有田块的边界对预处理后的遥感影像进行裁剪，得到土地整治项目区内所有田块的遥感影像，以上为成熟的遥感影像处理技术，此处不再详细说明，结果如图3所示。步骤3)利用ENVI(TheEnvironmentforVisualizingImages)软件中采用bandmath功能计算得到所有田块的归一化植被指数NDVI值，如图4所示，NDVI取值在0.89到0.32之间。步骤4)对检测区域的NDVI进行空间聚类，聚类使用ArcGIS的聚类分析模块实现。总计得到不显著(区内各个像元的Zi取值全部小于-1.96或大于+1.96时)、低值聚集(区内各个像元的Zi取值全部介于-1.96到+1.96之间，且Ii大于0，且Di小于0)、高值聚集(区内各个像元的Zi取值全部介于-1.96到+1.96之间，且Ii本文档来自技高网...
【技术保护点】
1.一种土地整治后农田作物生长缺陷区遥感快速诊断方法，其特征是，包括以下步骤：步骤1)在土地整治后的农作物生长季内，获取覆盖土地整治区的基于卫星、航空摄影或无人机的高分辨率多光谱遥感影像；步骤2)预处理所述遥感影像，包括对遥感影像进行大气校正和几何配准，将土地整治项目区的工程施工图与预处理后的遥感影像叠加在一起，然后利用工程施工图中的所有田块边界对预处理后的遥感影像进行裁剪，得到土地整治项目区内所有田块的遥感影像；步骤3)计算所有田块内每个像元的归一化植被指数NDVI，NDVI是近红外与红色通道反射比值的一种变化形式，是农作物生长状况的最佳指示因子，计算公式为：NDVI＝(NIR‑R)/(NIR+R)式中，NIR代表遥感影像近红外波段的反射率，R代表红光波段的反射率，NDVI取值范围在‑1到1之间，NDVI值越高，表明农作物生长状况越好；在ENVI即The Environment for Visualizing Images软件中采用band math功能实现对各个田块每个像元的NDVI值；步骤4)进行空间聚类分析，计算各个像元NDVI值空间关联的Local Moran’s I指数、各个像元NDVI值与所有田块NDVI平均值的差值Di、各个像元NDVI值标准化统计量Zi，计算公式分别为：...
【技术特征摘要】
1.一种土地整治后农田作物生长缺陷区遥感快速诊断方法，其特征是，包括以下步骤：步骤1)在土地整治后的农作物生长季内，获取覆盖土地整治区的基于卫星、航空摄影或无人机的高分辨率多光谱遥感影像；步骤2)预处理所述遥感影像，包括对遥感影像进行大气校正和几何配准，将土地整治项目区的工程施工图与预处理后的遥感影像叠加在一起，然后利用工程施工图中的所有田块边界对预处理后的遥感影像进行裁剪，得到土地整治项目区内所有田块的遥感影像；步骤3)计算所有田块内每个像元的归一化植被指数NDVI，NDVI是近红外与红色通道反射比值的一种变化形式，是农作物生长状况的最佳指示因子，计算公式为：NDVI＝(NIR-R)/(NIR+R)式中，NIR代表遥感影像近红外波段的反射率，R代表红光波段的反射率，NDVI取值范围在-1到1之间，NDVI值越高，表明农作物生长状况越好；在ENVI即TheEnvironmentforVisualizingImages软件中采用bandmath功能实现对各个田块每个像元的NDVI值；步骤4)进行空间聚类分析，计算各个像元NDVI值空间关联的LocalMoran’sI指数、各个像元NDVI值与所有田块NDVI平均值的差值Di、各个像元NDVI值标准化统计量Zi，计算公式分别为：式中，Xi、Xj分别为第i、j个像元的NDVI值，为所有田块NDVI的均值，n为像元...
【专利技术属性】
技术研发人员：杨永均，陈浮，李肖肖，张绍良，朱凤武，许桃元，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：江苏,32
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。