一种基于光学遥感数据的农田干旱空间分布评价方法与流程

1.本发明涉及农业灾害胁迫监测技术领域，具体为一种基于光学遥感数据的农田干旱空间分布评价方法。


背景技术：

2.传统的干旱监测通常是基于站点的实地观测，最常用的是通过站点测定的土壤水分含量进行插值评价土壤含水情况，或通过气象台站数据构建的降水距平百分率、标准化降水指数、气象干旱综合指数等指数进行插值，进而评价和判断干旱的发生和发展。这些传统方法虽然站点附近的结果真实性较高，但监测精度受控于地面站点的分布密度，较难获得准确及时的农业旱灾发生与发展的时空信息，并且往往要花费大量人力财力。通过定点观测来推测较大区域范围内的农田旱情时空动态变化，会因为站点分布不均使得插值结果推广到更大范围误差较大。而对于相对较小的区域范围，以点推面的农田旱情监测方法由于站点数据分布较少，也难以准确地判断农业旱灾的程度和具体受旱区域。遥感技术的快速发展，尤其是目前在轨的卫星传感器感测的电磁波段涵盖了可见光、近红外、热红外和微波等波段，为区域尺度农业干旱监测提供了新的手段。由于微波遥感受到空间分辨率低的限制，而热红外数据反演的地表温度会受到风速、气温、空气含水量等大气条件的影响，光学遥感数据的监测方法弥补两种数据的不足。
3.近些年来，国内外针对基于遥感技术的农业旱灾进行的监测和评估研究，提出了许多监测方法。农业旱灾遥感监测方法的优势在于它能快速、周期性地对大范围地区进行多波段、多角度的观测，获得连续的有关农田地表特征的面状信息。综合分析这些空间信息，建立相应的遥感反演模型，进一步了解大区域范围的农田土壤水分时空动态变化，结合农田作物长势，判断农业旱灾的发生与发展程度和范围。但是农业是一个复杂的生产系统，在相同干旱胁迫作用下，不同的农作物表现出不同的受旱灾损程度；另一方面，相同的作物，在不同的生长发育阶段对水分的需求也不相同，即不同生长阶段的作物耐旱程度也不相同。因此，虽然对于农田土壤水分以及农业旱灾遥感监测方法的研究已经开展了很多年，但作物生长时期、作物生长的气候、土壤、灌溉条件等因素的差异大都没有被考虑其中，旱灾监测的精度有待提高。


技术实现要素：

4.为解决现有的农业旱灾遥感监测方法没有考虑作物生长时期、作物生长的气候、土壤、灌溉条件等因素，旱灾监测的精度有待提高的问题，本发明公开一种基于光学遥感数据的农田干旱空间分布评价方法，该方法综合考虑了影响农业旱情发展的多种因素，并且能够针对不同作物以及作物所处的不同生长时期采取针对性的监测参数和监测标准，使农业干旱监测方法具有更好的普适性。
5.为实现上述目的，本发明采取如下的技术方案：
6.一种基于光学遥感数据的农田干旱空间分布评价方法，包括以下步骤：
7.步骤1：剔除旱情监测区域卫星影像中的非作物区域，根据作物种植的分布情况、作物类型及作物所处生长时期，将同一生长时期下的同一作物类型所在种植区域作为一个监测单元，将旱情监测区域卫星影像分成若干个监测单元；
8.步骤2：构建植被-土壤水分特征空间模型，计算监测区域内每个监测单元的模型参数，获得监测单元内每个像元的土壤含水量m；
9.步骤3：根据农业干旱性质考虑土壤含水量与相应农田种植区域内田间持水系数和永久萎蔫点系数之间的差异，并综合多种影响农业旱情的因素，构建农业旱情空间分布监测模型dsdi；
10.步骤4：区分作物以及作物的生长时期划分阈值，并划分所述阈值对应的旱情等级，将所述农业旱情空间分布监测模型dsdi的计算结果与阈值进行比较，得到不同类型作物、不同生长时期下的旱情监测区域农田的干旱程度及旱情空间分布评价结果，并绘制旱情监测区域的农业旱情空间分布图。
11.与现有技术相比，本发明所提供的基于光学遥感数据的农田干旱空间分布评价方法具有如下特点和有益效果：
12.(1)该方法适用于含有短波红外波段的卫星数据，卫星数据源更加丰富，不同分辨率的卫星遥感数据源使得模型适用于各种监测尺度；
13.(2)该方法在监测土壤含水量的基础上，还考虑了农作物种类、农作物所处生长时期、农作物所处的生长环境等影响作物发生旱情的因素，最终获得的农业旱情空间分布评价结果精度更高。
附图说明
14.图1为本发明提出的基于光学遥感数据的农田干旱空间分布评价方法流程图；
15.图2为植被-土壤水分特征空间模型原理图；
16.图3为本发明实例中得到的吉林省旱情监测区域的农业旱情空间分布图。
具体实施方式
17.下面将结合附图及较佳实施例对本发明的技术方案进行详细描述。
18.在其中一个实施例中，本发明提供一种基于光学遥感数据的农田干旱空间分布评价方法，如图1所示，该方法包括如下步骤：
19.步骤1：剔除旱情监测区域卫星影像中的非作物区域，根据作物种植的分布情况、作物类型及作物所处生长时期，将同一生长时期下的同一作物类型所在种植区域作为一个监测单元，将旱情监测区域卫星影像分成若干个监测单元。
20.本发明为主要针对农田旱情空间分布的监测与评价，评价区域确定为农田，因此在获取旱情监测区域卫星影像后，需要将卫星影像中的建筑、水体、裸地等非作物区域进行剔除，提取农作物种植区域，排除非农田地物对监测的干扰；另外，由于卫星影像受到云、雨等因素的限制，影像中可能存在的云及云阴影也需要作为异常区域剔除。剔除非作物区域和云及云阴影异常区域后的卫星影像作为后续处理的输入图像。
21.不同农作物由于种植方式、自身特性的不同对土壤水分的需求以及对旱情的敏感性也不同，因此本发明从卫星影像提取作物种植的分布情况(如水田作物和旱地作物)、作
物类型(如水稻、玉米、大豆、花生、小麦等)及作物所处生长时期(如拔节期、开花期等)。
22.根据上述提取结果，将同一生长时期下的同一种作物所在的种植区域确定为一个监测单元。如在提取到的农田种植区域的基础上，该种植区域内种植了玉米、大豆两种作物，且玉米处于拔节期、大豆处于开花期，则该种植区域可认为存在“拔节期-玉米”、“开花期-大豆”两个监测单元。
23.步骤2：构建植被-土壤水分特征空间模型，计算监测区域内每个监测单元的模型参数，获得监测单元内每个像元的土壤含水量m。
24.根据农田种植周期内影像的表现特点，农作物从播种、出苗到抽穗成熟表现为从裸土到植被全覆盖，植被的影响逐渐展现。从植被+土壤水分两个角度综合评价，利用遥感数据构建植被-土壤水分特征空间模型，监测土壤含水量。具体步骤为：
25.步骤21：利用已经剔除建筑、水体、裸地、云等非作物区域的卫星遥感数据计算每个监测单元内各像元的植被指数evi和基于短波红外波段的变换反射率str；根据模型参数构建特点，主要选用含有红光波段、蓝光波段、近红外波段以及短波红光波段的卫星影像。各像元的植被指数evi和基于短波红外波段的变换反射率str的计算公式如下：
[0026][0027][0028]
其中，r
red
、r
blue
、r
nir
、r
swir
分别为卫星影像的红光波段、蓝光波段、近红外波段和短波红外波段的地表反射率；
[0029]
步骤22：利用步骤21计算得到的evi和str构建植被-土壤水分特征空间，根据监测单元内作物所处的生长时期，确定该监测单元的特征空间的监测范围。例如，作物所处的生长时期为播种期、出苗期等低植被覆盖期，evi在0.1～0.5范围内与对应的str构建特征空间，然后再求取特征空间对应的干边拟合方程strd和湿边拟合方程strw；作物所处的生长时期为抽穗期、拔节期或成熟期等高植被覆盖期，evi在0.4～0.8范围内与对应的str构建特征空间，然后再求取特征空间对应的干边拟合方程strd和湿边拟合方程strw；
[0030]
步骤23：对于特征空间的干边拟合方程strd和湿边拟合方程strw的计算，基于evi和str的像素值在特征空间的分布，采用分位数回归的方法进行非线性指数函数数据拟合，计算方法为：
[0031][0032][0033]
其中，sw、iw为湿边拟合方程参数，sd、id为干边拟合方程参数；
[0034]
步骤24：根据干边拟合方程strd和湿边拟合方程strw的计算公式，计算监测单元内每个像元的土壤含水量m，计算公式为：
[0035][0036]
即：
[0037][0038]
步骤3：在土壤含水量监测结果的基础上，根据农业干旱性质考虑土壤含水量与相应农田种植区域内田间持水系数和永久萎蔫点系数之间的差异，并综合多种影响农业旱情的因素，构建农业旱情空间分布监测模型dsdi。
[0039]
农业干旱通常是指由于土壤水分持续不足而造成的作物体内水分亏缺，影响作物正常生长发育的现象。当蒸散损失超过植物有效土壤水分时，可用田间持水系数与永久萎蔫点系数之间的差异表征农业干旱。根据步骤2中土壤含水量反演结果，获得农业旱情空间分布监测模型dsdi的过程为：
[0040]
步骤31：利用rosetta土壤转换函数计算监测单元内每个像元的土壤的田间持水系数θ
fc
和永久萎蔫点系数θ
pwp
，土壤中砂粒含量、粉粒含量和黏粒含量来源于世界土壤数据库(hwsd，harmonizedworld soil database)；
[0041]
步骤32：根据监测单元内各像元的土壤含水量m与田间持水系数θ
fc
和永久萎蔫点系数θ
pwp
间的差异，构建农田土壤水分盈亏模型，该模型的公式为：
[0042][0043]
步骤33：构建农业旱情空间分布监测模型dsdi。
[0044]
具有良好的灌溉条件的农田，通常作物需水能够得到充分满足，不易发生农田干旱；无灌溉条件的农田，依靠降水补充作物需水，降水减少的情况下易发生干旱。坡度c、海拔d等地形因素，因为坡度较大海拔过高而易发生水分的流失，导致农田更易失水，相同降水条件下更易发生干旱。因此将灌溉条件、海拔、坡度因素综合加入农田土壤水分盈亏模型中，构建农业旱情空间分布监测模型dsdi，其公式为：
[0045][0046]
其中，b为swdi归一化后的值，c为农田种植区域的归一化后坡度值，d为农田种植区域的归一化后海拔值。
[0047]
步骤4：区分作物以及作物的生长时期划分阈值，并划分阈值对应的旱情等级，将农业旱情空间分布监测模型dsdi的计算结果与阈值进行比较，得到不同类型作物、不同生长时期下的旱情监测区域农田的干旱程度及旱情空间分布评价结果，并绘制旱情监测区域的农业旱情空间分布图。
[0048]
在相同干旱胁迫作用下，不同的农作物表现出不同的受旱灾损程度，另外，相同的作物，在不同的生长发育阶段对水分的需求也不相同，即不同生长阶段的作物耐旱程度也不相同。对于大多数作物而言，干旱胁迫对作物生产影响最大的两个阶段是播种出苗期和开花授粉期。前者是影响作物的适时播种与出苗率，没有适时种植和所需要的出苗率，就极难获得正常的收成。而在开花授粉灌浆早期，因干旱而导致授粉成功率和灌浆不足，直接影响最后的收成。因此，本发明通过阈值的划定以及旱情等级的划分，来评价不同作物、不同生长时期下的旱情监测区域农田的干旱程度及旱情空间分布情况，具体步骤为：
[0049]
步骤41：为评价监测单元干旱程度，区分作物、区分作物的生长时期来进行阈值的
划分，阈值分级标准参考《农业干旱等级》(gbit 32136-2015)，并将阈值对应的旱情等级划分为正常或湿润、轻度干旱、中度干旱、重度干旱以及特大干旱五个等级，利用作物特性进行对应的不同等级干旱程度阈值的划分；划分阈值以及对应的旱情等级后，将步骤3中农业旱情空间分布监测模型dsdi的计算结果与阈值进行比较，得到监测单元的干旱程度；
[0050]
步骤42：一个监测区域被分成若干个监测单元，将监测区域内多个监测单元合并，绘制该监测区域的农业旱情空间分布图。
[0051]
下面以吉林省2022年5月的农作物种植区域为例，对本发明的技术方案及技术效果作进一步说明。
[0052]
1、将旱情监测区域卫星影像分成若干个监测单元。
[0053]
利用吉林省2022年5月的terra/modis卫星遥感数据对吉林省进行耕种区域的提取以及农作物种植结构的提取，剔除建筑、水体、裸地、云及云阴影等非农田区域。为方便说明展示，主要选取了吉林省种植的大宗作物：水稻、玉米和大豆进行实例说明。5月的玉米和大豆农作物主要处于播种出苗期，水稻处于移栽分蘖期。将同一生长时期下的同一种作物所在种植区域确定为一个监测单元。因此，监测单元共分为三个：播种出苗期-玉米、播种出苗期-大豆以及移栽分蘖期-水稻。
[0054]
2、构建植被-土壤水分特征空间模型，计算监测区域内每个监测单元的模型参数，模型参数包括监测单元特征空间的干边、湿边拟合方程strd和strw，获得监测单元内每个像元的土壤含水量m。
[0055]
利用terra/modis卫星遥感数据计算每个监测单元内各像元的植被指数evi和基于短波红外波段的变换反射率str；根据三个监测单元内作物的生长时期，确定evi在0.1～0.5范围内与对应的str构建特征空间求取strd和strw；采用分位数回归的方法进行非线性指数函数数据拟合获取strd和strw；计算监测单元内每个像元的土壤含水量m。
[0056]
3、在土壤含水量监测结果的基础上，根据农业干旱性质考虑土壤含水量与该地块内田间持水系数和永久萎蔫点系数之间的差异，并综合多种影响农业旱情的因素，构建农业旱情空间分布监测模型dsdi。
[0057]
根据世界土壤数据库(hwsd，harmonizedworld soil database)中土壤中砂粒含量、粉粒含量和黏粒含量，利用rosetta土壤转换函数计算监测单元内每个像元土壤的田间持水系数θ
fc
和永久萎蔫点系数θ
pwp
；获取吉林省的dem高程数据并计算每个像元的坡度值；吉林省内水稻拥有良好的灌溉措施，玉米和大豆灌溉措施缺乏。
[0058]“移栽分蘖期-水稻”监测单元具有良好的灌溉设施，为水稻正常生长提供所需水分，dsdi大于0，不受农田干旱影响；
[0059]“播种出苗期-玉米”和“播种出苗期-大豆”两个监测单元通过以下公式获得干旱评价结果：dsdi＝0.65*b+0.15*c+0.1*d。
[0060]
4、评价不同作物、不同生长时期下农田的干旱程度及旱情空间分布情况，绘制吉林省旱情监测区域的农业旱情空间分布图。
[0061]“移栽分蘖期-水稻”监测单元由于具有良好的灌溉条件，无旱情发生，处于正常即湿润等级；“播种出苗期-玉米”和“播种出苗期-大豆”监测单元，两种作物处于的相同生长时期，且同样是需求不耐旱的作物，阈值分级标准适用相同的标准，分级标准如下：
[0062]
dsdi等级
(0，+∞)正常或湿润(-2，0)轻度干旱(-5，-2)中度干旱(-10，-5)重度干旱(-∞，-10)特大干旱
[0063]
通过上述过程，将监测区域内三个监测单元合并，得到吉林省水稻、大豆玉米种植区域的农业旱情空间分布结果，最终绘制出吉林省旱情监测区域的农业旱情空间分布图，如图3所示。
[0064]
本发明所提出的一种基于光学遥感数据的农田干旱空间分布评价方法所得到的农田旱情空间分布结果综合考虑了影响农业旱情发展的多种因素，并且方法能够针对不同作物以及作物所处的不同生长时期采取针对性的监测参数和监测标准，农业干旱空间分布评价精度更高，效果更好。
[0065]
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0066]
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

技术特征：
1.一种基于光学遥感数据的农田干旱空间分布评价方法，其特征在于，包括以下步骤：步骤1：剔除旱情监测区域卫星影像中的非作物区域，根据作物种植的分布情况、作物类型及作物所处生长时期，将同一生长时期下的同一作物类型所在种植区域作为一个监测单元，将旱情监测区域卫星影像分成若干个监测单元；步骤2：构建植被-土壤水分特征空间模型，计算监测区域内每个监测单元的模型参数，获得监测单元内每个像元的土壤含水量m；步骤3：根据农业干旱性质考虑土壤含水量与相应农田种植区域内田间持水系数和永久萎蔫点系数之间的差异，并综合多种影响农业旱情的因素，构建农业旱情空间分布监测模型dsdi；步骤4：区分作物以及作物的生长时期划分阈值，并划分所述阈值对应的旱情等级，将所述农业旱情空间分布监测模型dsdi的计算结果与阈值进行比较，得到不同类型作物、不同生长时期下的旱情监测区域农田的干旱程度及旱情空间分布评价结果，并绘制旱情监测区域的农业旱情空间分布图。2.根据权利要求1所述的一种基于光学遥感数据的农田干旱空间分布评价方法，其特征在于，步骤2包括以下步骤：步骤21：利用含有红光波段、蓝光波段、近红外波段以及短波红光波段的卫星影像，计算每个所述监测单元内各像元的植被指数evi和基于短波红外波段的变换反射率str，计算公式如下：公式如下：其中，r
red
、r
blue
、r
nir
、r
swir
分别为卫星影像的红光波段、蓝光波段、近红外波段和短波红外波段的地表反射率；步骤22：利用步骤21计算得到的各像元的植被指数evi和基于短波红外波段的变换反射率str构建植被-土壤水分特征空间，根据监测单元内作物所处的生长时期确定特征空间的监测范围；步骤23：采用分位数回归的方法进行非线性指数函数数据拟合，得到特征空间的干边拟合方程str
d
和湿边拟合方程str
w
：：其中，s
w
、i
w
为湿边拟合方程参数，s
d
、i
d
为干边拟合方程参数；步骤24：根据干边拟合方程str
d
和湿边拟合方程str
w
计算监测单元内每个像元的土壤含水量m，计算公式如下：3.根据权利要求1或2所述的一种基于光学遥感数据的农田干旱空间分布评价方法，其特征在于，步骤3包括以下步骤：
步骤31：利用rosetta土壤转换函数计算监测单元内每个像元的田间持水系数θ
fc
和永久萎蔫点系数θ
pwp
；步骤32：根据各像元的土壤含水量m、田间持水系数θ
fc
和永久萎蔫点系数θ
pwp
构建农田土壤水分盈亏模型，模型公式如下：步骤33：将灌溉条件、海拔、坡度因素综合加入所述农田土壤水分盈亏模型中，构建如下的农业旱情空间分布监测模型dsdi：其中，b为swdi归一化后的值，c为农田种植区域的归一化后坡度值，d为农田种植区域的归一化后海拔值。

技术总结
本发明公开一种基于光学遥感数据的农田干旱空间分布评价方法，包括：剔除旱情监测区域卫星影像中的非作物区域，将影像分成若干个监测单元；构建植被-土壤水分特征空间模型，计算监测单元内各像元的土壤含水量；根据农业干旱性质考虑土壤含水量与田间持水系数和永久萎蔫点系数之间的差异，并综合多种影响农业旱情的因素，构建农业旱情空间分布监测模型DSDI；根据DSDI评价不同作物、不同生长时期下的农田干旱程度及旱情空间分布情况，绘制旱情监测区域的农业旱情空间分布图。本发明综合考虑了影响农业旱情发展的多种因素，并能够针对不同作物、不同生长时期采取针对性的监测参数和监测标准，干旱空间分布评价精度更高，效果更好。更好。更好。


技术研发人员：武传昱 朱瑞飞 巩加龙 许照鑫
受保护的技术使用者：长光卫星技术股份有限公司
技术研发日：2023.03.15
技术公布日：2023/7/31