一种丘陵山区水田种植模式遥感快速监测方法、系统及计算机可读存储介质技术方案

本发明专利技术公开一种丘陵山区水田种植模式遥感快速监测方法、系统及存储介质，包括：步骤1，对耕作初期获取的多光谱光学遥感影像进行预处理，包括进行辐射定标、大气校正和地形辐射校正；步骤2，根据定量化遥感影像提取初始水体，并对初始水体进行非监督分类，然后统计各类别的光谱均值和方差，最后基于光谱均值和方差建立地表水体指数反演模型，计算地表水体指数影像SWI，监测地表是否含有水体；步骤3，叠加土地利用现状数据和地表水体遥感指数影像SWI，统计水田图斑范围内像元的地表水体指数大小，判别水田图斑的种植模式。本发明专利技术具有人工参与少、自动化程度高，提取结果准确性高，可快速区分水田的种植模式的优点。快速区分水田的种植模式的优点。快速区分水田的种植模式的优点。

【技术实现步骤摘要】
一种丘陵山区水田种植模式遥感快速监测方法、系统及计算机可读存储介质


[0001]本专利技术属于遥感图像处理以及自然资源管理
，具体涉及丘陵山区水田种植模式的监测技术。

技术介绍

[0002]水田是土地利用现状分类体系中耕地类别下的二级类，是指筑有田埂(坎)，可以蓄水，用于种植水稻、莲藕、荸荠、莼菜等水生作物的土地。丘陵山区的水田大多为梯田，主要种植水稻。丘陵山区地形地貌起伏破碎，水田面积小且分散，难以进行规模化、机械化种植水稻，因此在实际耕作中，因种植经济效益低、农村劳动力不足等因素，较多水田没有种植水稻，而改为旱作(种植旱生作物，如玉米)，或撂荒等，影响国家主粮安全。此外，丘陵山区的水田若长期旱作，可能导致田坎破坏、区域地下水位下降等，从而影响水田的蓄水能力，导致耕地的等级下降。因此，水田的种植模式(水作与旱作)监测，是丘陵山区耕地保护工作中的重要内容。
[0003]目前，丘陵山区水田种植模式监测主要有野外调查法和基于高分辨率遥感影像的人工解译法。这两种方法具有各自的不足：
[0004]野外实地调查方法是指调查人员携带定位设备到野外实地查看水田的种植模式，并记录相关信息。这种方法需要投入大量人力和财力，难以应用于大范围的监测，一般只能选择少量样本进行抽样监测。
[0005]基于高分辨率遥感影像的人工解译法是采用一期或多期高分辨率遥感影像，技术人员在专业软件中叠加遥感影像和土地利用现状数据库的水田范围，根据遥感影像的光谱和纹理信息对水田的种植模式进行遥感人工解译。这种方法的效率虽比外业调查高，可用于丘陵山区大范围水田种植模式监测，但仍存在耗费时间长、效率低的问题，同时遥感解译结果的正确性依赖于解译人员的经验，解译人员的经验和主观性差异必然影响结果精度。
[0006]中国专利文献CN112507763A公开了一种基于多源多光谱遥感图像水体提取方法，该方法包括：获取多源遥感图像，并分别对遥感对象进行尺寸预处理；对预处理后的遥感图像进行图像大气校正；对大气校正后的遥感图像进行归一化水体指数NDWI计算；消除归一化水体指数NDWI的异常值，并绘制NDWI图像灰度直方图，根据NDWI图像灰度直方图峰值分布确定水体分割阈值；根据归一化水体指数NDWI与水体分割阈值的预设关系判断遥感图像中的地物是否为水体。但该方法所述的归一化水体指数NDWI仅用了多光谱遥感数据中的两个波段，用于识别河湖库塘等较深的水体较为有效，但难以监测水田中的浅水体。

技术实现思路

[0007]针对现有方法效率低、成本高、主观性大等问题，本专利技术的目的是提供一种基于多光谱光学遥感数据的丘陵山区水田种植模式定量遥感监测方法、系统及计算机可读存储介质，减少人工参与程度、提高自动化程度，提高提取结果的准确性，快速区分水田的种植模
式，为丘陵山区耕地保护工作提供技术支撑。
[0008]本专利技术方法是采用耕作初期获取的多光谱光学遥感数据和土地利用现状数据，建立的一种基于多光谱光学遥感数据的丘陵山区水田种植模式定量遥感监测方法，利用光学遥感影像获得多光谱光学遥感数据，其遥感影像是通过航天(卫星)或航空(飞机)成像方式，拍摄的地表的影像(照片)，其实质是获取地面目标的反射、散射或辐射的电磁波信息。光学遥感影像是指获取可见光到红外线这部分电磁波信息进行成像，得到的影像。
[0009]本专利技术的技术方案如下：
[0010]一种丘陵山区水田种植模式遥感快速监测方法，所述方法包括如下步骤：
[0011]步骤1，多光谱光学遥感影像定量化预处理
[0012]对耕作初期获取的多光谱光学遥感影像进行预处理，包括进行辐射定标、大气校正和地形辐射校正，获得地表反射率遥感影像SR，使得遥感影像能定量反映地表真实信息。
[0013]步骤2，地表水体指数反演
[0014]根据定量化遥感影像提取初始水体，并对初始水体进行非监督分类，然后统计各类别的光谱均值和方差，最后基于光谱均值和方差建立地表水体指数反演模型，然后计算地表水体指数影像SWI，根据地表水体指数监测地表是否含有水体。
[0015]步骤3，水田图斑种植模式的自动判别
[0016]叠加土地利用现状数据和地表水体遥感指数影像，统计水田图斑范围内像元的地表水体指数大小，利用水体指数影像SWI对土地利用现状数据中的水田图斑进行计算，判别水田图斑的种植模式。
[0017]进一步，所述步骤1中，辐射定标是利用多光谱光学遥感影像原始数据的像元亮度值，结合遥感影像元数据中的增益和偏差系数转换为辐亮度数据L，然后利用公式计算星上反射率
[0018][0019]式中d为日地距离，L为辐亮度数据，单位是w/m2/sr/μm，E0为太阳常数，β为成像时太阳天项角。
[0020]进一步，所述步骤1中，大气校正是采用辐射传输模型法，公式如下
[0021][0022]式中SR为地表反射率，TOA为星上反射率，E0为太阳常数，L
path
为大气的路径辐射，β、v分别为太阳天顶角和卫星观测天顶角，T(β)、T(v)为分别为大气下行辐射总透过率和上行辐射总透过率，S是大气底层向下的半球反照率。
[0023]进一步，所述步骤1中，地形辐射校正是采用C辐射校正模型结合DEM数据来消除地形的影响，校正公式为：
[0024][0025]c＝a/b
ꢀꢀꢀꢀ
(6)
[0026]其中，SR
T
是倾斜地表(即校正前)的反射率，a和b为拟合方程的系数，i为太阳入射角，β为太阳天顶角。
[0027]进一步，所述步骤2的地表水体指数反演具体包括如下步骤：
[0028]步骤2.1，初始水体识别：采用归一法或人工交互解译方法提取初始水体。
[0029]步骤2.2，初始水体非监督分类：将提取的所有初始水体像元利用k
‑
均值算法、核k
‑
均值算法或迭代自组织数据分析技术(ISODATA)等进行非监督分类。
[0030]步骤2.3，初始水体参数统计：结合初始水体非监督分类结果和地表反射率影像，逐波段、逐类别统计初始水体的光谱反射率均值Mean
i，j
和均方差STD
i，j
，其中i为不同的波段，j为不同的类别。
[0031]步骤2.4，计算地表水体指数影像SWI
[0032][0033][0034]式中i为不同的波段，j为不同的类别，SR
i
为波段i的地表反射率值，m为波段数，n为类别数。
[0035]更进一步，所述步骤2.1，采用归一法提取初始水体具体如下：
[0036]基于地表反射率影像，逐像元计算归一化水体指数MNDWI，红、绿、近红外三波段最小值Min_RGN，红、绿、蓝、近红外和短波红外最大值Max_Val，公式如下
[0037][0038]Min_RGN＝min{SR
R
，SR
G
，SR
NIR本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种丘陵山区水田种植模式遥感快速监测方法，其特征在于，所述方法包括如下步骤：步骤1，多光谱光学遥感影像定量化预处理对耕作初期获取的多光谱光学遥感影像进行预处理，包括进行辐射定标、大气校正和地形辐射校正，获得地表反射率遥感影像SR，使得遥感影像能定量反映地表真实信息；步骤2，地表水体指数反演根据定量化遥感影像提取初始水体，并对初始水体进行非监督分类，然后统计各类别的光谱均值和方差，最后基于光谱均值和方差建立地表水体指数反演模型，并应用于地表反射率影像SR计算地表水体指数影像SWI，根据地表水体指数监测地表是否含有水体；步骤3，水田图斑种植模式的自动判别叠加土地利用现状数据和地表水体遥感指数影像，统计水田图斑范围内像元的地表水体指数大小，利用水体指数影像SWI对土地利用现状数据中的水田图斑进行计算，判别水田图斑的种植模式。2.根据权利要求1所述的丘陵山区水田种植模式遥感快速监测方法，其特征在于，所述步骤1中，辐射定标是利用多光谱光学遥感影像原始数据的像元亮度值，结合遥感影像元数据中的增益和偏差系数转换为辐亮度数据L，然后利用公式计算星上反射率式中d为日地距离，L为辐亮度数据，单位是w/m2/sr/μm，E0为太阳常数，β为成像时太阳天顶角。3.根据权利要求1所述的丘陵山区水田种植模式遥感快速监测方法，其特征在于，所述步骤1中，大气校正是采用辐射传输模型法，公式如下式中SR为地表反射率，TOA为星上反射率，E0为太阳常数，L
path
为大气的路径辐射，β、v分别为太阳天顶角和卫星观测天顶角，T(β)、T(ν)为分别为大气下行辐射总透过率和上行辐射总透过率，S是大气底层向下的半球反照率。4.根据权利要求1所述的丘陵山区水田种植模式遥感快速监测方法，其特征在于，所述步骤1中，地形辐射校正是采用C辐射校正模型结合DEM数据来消除地形的影响，校正公式为：c＝a/b
ꢀꢀꢀꢀ
(6)其中，SR
T
是倾斜地表即校正前的反射率，a和b为拟合方程的系数，i为太阳入射角，β为太阳天顶角。5.根据权利要求1所述的丘陵山区水田种植模式遥感快速监测方法，其特征在于，所述步骤2的地表水体指数反演具体包括如下步骤：步骤2.1，初始水体识别：采用归一法或人工交互解译方法提取初始水体；步骤2.2，初始水体非监督分类：将提取的所有初始水体像元利用k
‑
均值算法、核k
‑
均
值算法或迭代自组织数据分析技术(ISODATA)进行非监督分类；步骤2.3，初始水体参数统计：结合初始水体非监督分类结果和地表反射率影像，逐波段、逐类别统计初始水体的光谱反射率均值Mean
i，j
和均方差STD
i，j
，其中i为不同的波段，j为不同的类别；步骤2.4，计算地表水体指数影像SWI步骤2.4，计算地表水体指数影像SWI式中i为不同的波段，j为不同的类别，SR
i
为波段i的...

【专利技术属性】
技术研发人员：胡勇，邓琳，熊先才，李晓俊，周川，黄健，周志跃，
申请(专利权)人：重庆市规划和自然资源调查监测院，
类型：发明
国别省市：