基于神经网络的灌溉空间分布监测方法研究

摘要  利用神经网络对光学数据和雷达数据样本进行训练，以山东省胶州市和桓台县为研究区域，建立以哨兵1号和哨兵2号数据为数据源的灌溉监测模型，将监测结果与实地考察数据进行对比验证，并评价所建立的监测模型。结果表明，该模型在研究区域的灌溉空间分布监测能得到质量较高的结果，与使用单一雷达数据相比在高植被覆盖区域也能够有较高的监测精度，为今后的灌溉空间分布监测研究提供了新方向。

关键词  哨兵1号;哨兵2号;神经网络;灌溉空间分布

中图分类号  S 127   文献标识码  A    文章编号  0517-6611（2022）15-0183-05

doi： 10.3969/j.issn.0517-6611.2022.15.047

开放科学（资源服务）标识码（OSID）：

Research on Monitoring Method of Irrigation Spatial Distribution Based on Neural Network

ZHAO Xu-dong1，2，ZHUANG Hui-bo3，CHEN Gan-qin3 et al

（1.Shandong University of Science and Technology，Qingdao，Shandong 266500;2.Qingdao Star-rising Technologies Co.，Ltd.，Qingdao，Shandong 266500;3.Shandong Hydrological Center，Jinan，Shandong 250000）

Abstract  Using neural network to train optical data and radar data samples，taking Jiaozhou City and Huantai County in Shandong Province as the research area，an irrigation monitoring model with Sentinel No.1 and Sentinel No.2 data as data sources was established，and the monitoring results with on-site inspections data were compared，and the established monitoring model was evaluated.The results showed that the model could obtain higher quality results in the monitoring of irrigation spatial distribution in the study area.Compared with the use of single radar data，it can also have higher monitoring accuracy in areas with high vegetation coverage.The study provided a new direction for future research on irrigation spatial distribution monitoring.

Key words  Sentinel-1;Sentinel-2;Neural network;Irrigation spatial distribution

山东省农业用水量大，农业用水占经济社会用水比例高，对农业灌溉空间分布进行监测有助于提高农业用水量核算精度，对用水总量控制和最严格水资源管理具有重要意义。早期的灌溉空间分布统计主要通过人工获取，该方法人力物力消耗大，时效性差，空间分布信息缺失，导致灌溉信息难以及时获取。

近年来，遥感技术的快速发展为地面信息的快速获取提供了新的技术手段，且越来越多地应用于灌溉面积空间分布提取。遥感具有周期短、覆盖范围大和信息获取速度快、信息采集实时性强等特性，可有效弥补地面观测系统成本高、空间代表性差等缺点，为灌溉空间分布监测提供了基于面上的实时有效的监测手段和方法。目前基于遥感手段进行灌溉空间分布提取的应用主要有2种，一种是根据灌区的土地利用和农作物种植结构信息，结合灌溉工程控制范围来获取可能的灌溉分布范围;刘逸竹等[1]利用时间序列NDVI进行灌溉面积统计，研制了我国2010年、空间分辨率250 m的灌溉耕地空间分布地图。

另一种是基于遥感土壤含水量变化信息获取实际灌溉面积的方法，通过分析和对比土壤含水量的变化特征，获得轮次实际灌溉面积[2-3]。目前利用遥感进行土壤含水量的监测主要分为3类：热红外法[4]、可见光/近红外法（反射率法、植被指数法）和微波方法[5]。针对不同的地表类型可以选择不同的模型计算土壤水分，近年来，在裸露地表通常使用一些物理经验模型（oh模型[6]、Dubois 模 型[7]）进行土壤含水量估测，但需要现场测量进行校准;在低植被覆盖度地区，通常选择将半经验模型（Shi模型[8] 、水云模型[9]等）应用于水分反演中，植被区土壤含水量仍被过低估计;对于植被高覆盖度地区，植被冠层双向散射衰减了土壤后向散射的影响，需要采用密歇根微波植被散射模型[10]（michigan microwave canopy scattering，MIMICS）及其改进模型（Bi-MIMICS模型 [11]）计算裸土后向散射系数，由于参数太过烦琐不易于使用。研究表明当微波遥感数据协同光学植被指数[12]的方式  反演土壤含水量时，用光学数据计算植被生物物理参数，能更好地补充植被的影响。

合成孔径雷达具有全天时、全天候的观测能力，且对植被和土壤表层具有一定的穿透能力，已被广泛应用于地表土壤水分的监测和反演中[13]。但雷达后向散射系数除受到土壤介电常数的影响外，还受到地表粗糙度、植被覆盖等的影响，本质上土壤水分反演是一个病态问题。为了解决病态反演的问题，提出利用神经网络深入挖掘雷达观测数据结合光谱数据与地表是否灌溉之间的关系，从而得到地面灌溉结果。人工神经网络具有连续的学习能力，在无法建立完整理论模型的应用中具有显著优势，尤其是处理非线性关系时。研究表明，神经网络在经过大量数据集的训练后，由于其有效的插值能力，可以提供良好的估计性能。

1 研究区概况与遥感数据来源

1.1 研究区概况

为建立神经网络提取耕种地区灌溉空间分布情况，以山东省胶州市（119°37′～120°12′E、36°00′～36°30′N）和淄博市桓台县（117°50′～118°11′E、36°51′～37°06′N） 2处作为试验区域。2个研究区域以平原和丘陵为主，由图1、2可知，2个区域的地表覆盖类型主要包括农业用地和人造地表2种，大范围的农业用地用于灌溉分布的提取。

1.2 遥感数据

1.2.1    哨兵1号。

为确定灌溉的空间分布情况，通过估算试验区的土壤含水量从而确定灌溉情况，试验中使用Sentinel-1数据。哨兵1号（Sentinel-1）卫星是欧洲航天局哥白尼计划（GMES）中的地球观测卫星，由2颗极轨卫星（Sentinel-1A和Sentinel-1B）组成，载有C波段合成孔径雷达。根据工作模式，Sentinel-1 可在全天时、全天候条件下获取5～40 m分辨率的图像，重访周期由1颗星的12 d减到2颗星的6 d。 携带的合成孔径雷达（SAR）有4种成像模式（表1），该研究所用的数据是干涉宽模式（Interferometric wide swath，IW）下的GRD数据。

1.2.2    哨兵2号。

哨兵二号是高分辨率多光谱成像卫星，高度为786 km用于陆地监测，可提供植被、土壤和水覆盖、内陆地区及海岸区域等图像，还可以用于紧急救援服务。与哨兵一号相同，哨兵二号也拥有2颗卫星（Sentinel-2A和Sentinel-2B），高度为786 km，可覆盖13个光谱波段，从可见光和近红外到短波红外，幅宽达290 km。地面分辨率分别为10、20和60 m，一颗卫星的重访周期为10 d，2颗互补，重访周期为5 d。

1.3 数据预处理及训练样本制作

1.3.1    数据预处理。

数据预处理主要包括对光学数据和雷达数据的处理，目的是获得模型输入需要的反射率和后向散射系数。为了获得地面的后向散射系数，需要对雷达数据（哨兵1号）进行地理编码和辐射定标，利用SNAP软件批量完成。辐射定标的公式：

σoij=10log  DN2ij A2σ

其中，σo为后向散射系数;i和j分别表示第i行和第j列;DN为雷达图像的灰度值;A为定标参数。

光学数据选择哨兵2号的L2A数据，该数据已完成大气校正、几何校正，只需进行波段合成、裁剪和镶嵌的操作，也通过SNAP软件批量操作完成（表2）。

1.3.2    训练样本。

神经网络模型的精度严重依赖训练样本，依靠大量可靠性高的训练样本为网络模拟结果做支撑，因此在选择样本时，为了提高网络模型的适应性，利用地面灌溉信息选择遥感数据，得到可靠的灌溉样本信息。

利用灌区的灌溉周期信息，内蒙古黄河灌区的解放闸灌域从2021年4月10日开始一水的灌溉，5月23日前后结束二水灌溉。利用灌区灌溉信息判断灌溉区域与非灌溉区域，对覆盖灌区的哨兵1号数据和哨兵2号数据进行预处理获得反射率和后向散射系数，提取训练数据的输入数据。

由经验知识以及试验证明，后向散射系数与土壤水分含量相关，土壤水分含量的多少往往能表示是否进行灌溉，反射率和NDVI值则能反映地面植被覆盖情况。最终选择后向散射系数、可见光波段反射率、近红外波段反射率以及NDVI值作为模型输入特征向量，NDVI值用哨兵2号的4波段和8波段计算。利用大型灌区的灌溉时间和灌溉分布情况，获取相应时间和区域的数据，对遥感影像进行预处理，包括辐射校正、几何校正和波段计算等处理过程，选择灌溉样本制作灌溉类型样本数据集。输入向量通过非灌溉样本选择云、水体、建筑用地、森林以及非灌溉的农业用地，根据输入特征向量制作相应的非灌溉样本数据集。

NDVI=（ρband8-ρband4）/（ρband8+ρband4）

1.4 灌溉分布监测模型设计

1.4.1    神经网络。

人工神经网络是由众多神经元使用可调节的连接权值连接而成。一个神经网络通常由一个输入层，一个或多个隐藏层和一个输出层连接构成[14]。其中隐藏层的层数往往决定神经网络的数据挖掘能力，输入层神经元个数为影响结果的特征向量个数，输出层的神经元个数由输出结果决定（图3）。

神经网络具有高处理速度、高容错能力，能够较好地解决非线性问题的优化解，对数据的依赖性较强。神经网络在遥感中应用广泛，主要应用范围有图像识别、目标检测、遥感图像分类、参数定量反演等领域。神经网络能够以较短的计算时间和较高的精度模拟各变量之间的非线性关系，因此该研究利用神经网络深入挖掘灌溉与非灌溉像元之间的特征差异从而实现对灌溉区域的提取。