光学与SAR遥感协同反演植被覆盖区土壤含水量.pdf

收稿日期：基金项目：山东省自然科学基金资助项目（）作者简介：杨晶晶（），女，陕西富平人，硕士研究生，研究方向为遥感科学与技术 通信作者：李莎（），女，山东荣成人，高级工程师，主要从事水工环地质方面的研究工作：【水环境与水生态】光学与 遥感协同反演植被覆盖区土壤含水量杨晶晶，邓清海，李 莎，张丽萍，陈 桥，孙桂宗，孙振洲（山东科技大学 地球科学与工程学院，山东 青岛；山东省地质矿产勘查开发局 第七地质大队，山东 临沂）摘 要：在进行土壤含水量反演时，单纯使用传统遥感反演模型很难有效消除干扰因素。以山东省东营市为研究区，基于光学遥感与合成孔径雷达（）数据，采用植被光谱指数修正水云模型中的植被含水量，并将修正后的水云模型与高级积分方程模型（）耦合，以消除植被含水量和土壤粗糙度对土壤含水量反演结果的影响，从而达到提高遥感模型反演土壤含水量精度的目的。结果表明：基于比值植被指数（）的二次函数修正水云模型后，与 模型耦合反演土壤含水量的精度最高，决定系数大于，均方根误差（）小于等于；土壤含水量在空间上呈现西北部大，向南逐渐减小的连续空间分布特征，该耦合模型具有普适性。关键词：植被光谱指数；植被含水量；模型；多源遥感协同反演中图分类号：；文献标志码：引用格式：杨晶晶，邓清海，李莎，等光学与 遥感协同反演植被覆盖区土壤含水量人民黄河，（）：，（，；，）：，（），（）（），（），（），：；（）；旱情监测是农作物产品估值研究等的基础。近几十年来，多源遥感对地监测技术迅速发展，突破了传统土壤含水量监测的局限性，使得大面积连续土壤含水量监测成为可能。微波遥感技术可以穿透地表植被覆盖对地进行监测，并利用微波信号与土壤物理性质之间的高度相关性，灵敏地探测土壤含水量的动态变化。基于此，有学者利用可见光、短红外、近红外图像结合微波探测技术，研究得出了一系列反演土壤含水量模型，如常用的水云模型（）、模型以及高级积分方程模型（）等，其中 模型被广泛应用于植被覆盖区土壤含水量反演、第 卷第 期 人 民 黄 河 ，年 月 ，模型可以去除地表粗糙度对土壤含水量反演效果的影响。传统模型反演土壤含水量受多种干扰因素影响，使用一种模型很难保证反演精度，且模型反演参数在不同气候条件下受干扰程度不同，逐一剔除干扰因素提高反演精度的难度较大。因此，可根据光学遥感影像采集植被层信息，选择合适的植被光谱指数修正水云模型中的植被含水量，并耦合高级积分方程模型，去除基于雷达后向散射系数计算的地表粗糙度的影响，从而准确获取土壤含水量情况，最后结合光学遥感连续对地观测优势，验证修正后的水云模型耦合 模型反演土壤水分的普适性。黄河下游地区是我国粮食生产核心区，近年来城市扩张给当地生态系统造成了巨大影响。土壤含水量直接影响地区生态环境和农作物生长，实时监测土壤含水量分布，对保护黄河下游地区生态系统以及农业发展具有重要意义。因此，本文以黄河下游山东省东营市为研究对象，基于遥感技术探究植被覆盖区土壤含水量，以期为其生态环境保护和农业发展规划编制提供参考。数据预处理山东省东营市位于黄河入海口，属雨热同季的暖温带大陆性季风气候区。年降水量 ，降水集中在 月。该地区地形起伏不大，地势沿黄河自西南向东北倾斜。黄河自西南向东北贯穿全市，在垦利区东北注入渤海。光学遥感数据 属于光学遥感卫星，搭载陆地成像仪（）和热红外传感器（），其中 对波段重新调整，排除了 处水汽吸收特征，具有丰富的波段信息。山东省东营市 月降水量少，因此选择该地区 年 月 日、年 月 日、年 月 日三景 影像解译东营市植被覆盖情况。原始数据使用 软件进行图像融合、影像裁剪、辐射定标和大气校正等预处理，以增强波段特征，其空间分辨率为 。影像数据本文采用的 影像数据为来自欧洲太空局（）的 地距影像（）数据，其空间分辨率为 。对照 数据的日期进行下载，使用哨兵数据应用 对原始数据进行预处理，并计算雷达后向散射系数。预处理后的 影像的雷达后向散射特征更为明显，能提供基础计算参数，其中雷达入射角为。实测数据研究区实测地表含水量信息从 土壤湿度数据集中获取。该数据集通过同化地表土壤湿度、植被叶面积指数和植被光学深度信息，使用国际土壤水分观测网和现有各区土壤湿度数据集对生产目标范围内的数据。此次下载使用的实测土壤水分数据为 版本，数据空间分辨率为 。研究方法 植被含水量反演模型通过对比分析，选择合适的植被光谱指数构建植被含水量反演模型。基于光学遥感的多光谱信息，常用植被持水度（）、叶片含水量（）和相对含水量（）模拟植被含水量，其中 与植被含水量相关性最强，常用于构建相关函数、反演植被含水量，因此选用 修正水云模型中的植被含水量。现阶段，根据光学遥感数据的波段特征，计算集中在红外、近红外波段的植被光谱指数，是基于光学遥感数据表达植被含水量的较为成熟的方法。本文选择的植被光谱指数有增强植被指数（）、归一化植被指数（）、归一化水体指数（）、比值植被指数（）、红边叶绿素指数（）和归一化差值水分指数（），各指数反演公式见表（其中、为 卫星的可见光波段反射率，对应波长分别为 、；为近红外波段反射率，其波长为 ；、为短波红外波段反射率，对应波长分别为、）。表 植被光谱指数及其反演公式植被光谱指数符号计算公式增强植被指数（）（）归一化植被指数（）（）归一化水体指数（）（）比值植被指数 红边叶绿素指数（）归一化差值水分指数（）（）初步构建植被光谱指数与植被含水量之间的关系时，发现基于植被光谱指数的二次多项式反演的植被含水量与实测植被含水量相关性最强，因此植被含水人 民 黄 河 年第 期量反演模型最高次数为二次。模型高级积分方程模型可以准确捕获地表辐射的雷达后向散射系数信息，相较于其他模型能够更真实地模拟较广范围的地表粗糙度。由于获取建立 模型的完整参数比较困难，因此通过模拟数据建立土壤含水量和雷达后向散射系数之间的半经验关系，进而使用查找法确定 模型中地表粗糙度系数。耦合模型采用植被光谱指数表达植被含水量修正水云模型，从而剔除植被含水量对土壤含水量反演结果的影响，再耦合 模型进一步剔除地表粗糙度的影响，最终得出反演该地区土壤含水量的耦合模型。本文使用的 数据为一级 数据，包括、两种极化方式，雷达入射角为。使用 模型反演地表粗糙度，将模拟波段和地表粗糙度等相关参数输入模型，根据研究区实际情况选择参数并修正模型，得到关于地表粗糙度的雷达后向散射系数：（）（）式中：和 分别为 数据的、极化参数。水云模型将雷达后向散射系数定义为植被后向散射系数和受植被遮挡的后向散射系数两部分。本文采用东营市 月植被覆盖度普遍较低时的遥感数据，根据 等于 年提出的经验参数，即令水云模型中的植被参数 、，修正水云模型。使用 表达植被含水量，并与 模型耦合反演得到土壤含水量：（）（）（）结果与讨论 植被含水量影响分析为了确保修正的水云模型中植被含水量满足随机验证精度，将 年的样本数据分为 组，每组 个样本点，进行植被含水量初步计算，并与实测数据进行对比，验证不同植被光谱指数表达植被含水量的精度和稳定性，进而确定修正的水云模型中植被含水量的计算公式。本文采用计算数据拟合的均方根误差（）和决定系数（）评价耦合模型的反演精度、反演结果的代表性和稳定性，并据此评价模型参数的适用性，计算结果见表（表示植被光谱指数）。表 基于植被光谱指数的植被含水量计算公式及精度植被光谱指数植被含水量计算公式 （）（）（）（）（）（）由表 可知，除了 外，使用二次多项式反演植被含水量的决定系数最大，其中由 反演植被含水量的决定系数最大，但均方根误差也较大，说明基于 构建的二次多项式反演植被含水量精度高，但稳定性差，普适性不强；由 构建二次多项式反演植被含水量的决定系数较大，同时均方根误差最小，说明该二次多项式能够在保证反演精度的前提下，确保模型的稳定性和普适性。因此基于 的二次多项式表达植被含水量修正水云模型，并耦合 模型反演土壤含水量，即将 代入式（）求得土壤含水量。验证模型精度根据耦合模型反演 年 月山东省东营市土壤含水量，并与实测土壤含水量数据进行对比，验证耦合模型精度。为突出模型修正方法的适用性，对样本点进行均匀分区，每一区的样本点都能包含东营市所有的植被覆盖特征，并适当调整耦合模型中植被含水量表达式的系数，共进行 组试验，耦合模型中植被含水量计算公式以及土壤含水量反演精度见表（土壤含水量公式略）。该耦合模型反演的土壤含水量较稳定，决定系数均大于，并且均方根误差仅为，说明模型具有较强普适性。耦合模型的应用根据耦合模型，选取 年典型植被覆盖区内一组样本点反演土壤含水量（以土壤体积含水量表示），见图。其中，耦合模型植被含水量修正公式为 人 民 黄 河 年第 期。西北部沿海种植区土壤含水量最大，为 ；受黄河水补给与灌溉影响，黄河沿岸水稻种植区土壤含水量较大，为 ；南部内陆地区地形起伏较明显，受地表径流、距海距离远等因素影响，土壤含水量最小，仅为 。西部内陆地区部分土壤含水量反演值与实测值相比偏高。原因是，下载的土壤含水量实测数据与遥感数据的分辨率不一致，且研究区周遭城市土壤含水量实测值对该区域实测数据的裁剪和收集产生了一定影响，但反演结果整体较好，说明模型适用性强。表 耦合模型反演 年东营市 个样本点土壤含水量结果统计序号耦合模型中植被含水量计算公式土壤含水量反演精度 图 东营市 年随机样本反演精度验证与空间分布 为了避免验证结果的偶然性，对植被覆盖区选取的所有样本点进行逐年统计，通过对比最终反演结果的空间分布，选取最具代表性的 个样本点展示该地区土壤含水量空间分布情况，进行土壤含水量空间变化分析。利用该耦合模型反演 年和 年土壤含水量，对耦合模型中植被含水量参数进行调整以确保反演结果的精度和可靠性，反演结果见图、图。其中，年、年反演土壤含水量的耦合模型中，植被含水量计算公式分别为 ，。图 东营市 年随机样本反演精度验证与空间分布图 东营市 年随机样本反演精度验证与空间分布 用该耦合模型反演 年和 年土壤含水人 民 黄 河 年第 期量，决定系数大于，精度较高，但较 年的反演精度有所降低，且均方根误差增大。结合整体土壤含水量反演结果的空间分布情况得知，年东营市北部地区土壤干旱情况加重，其中土壤含水量普遍较高（）的地区，土壤含水量以 的速率降低，在整体空间范围内变化最为突出。此外，经过土壤含水量反演发现，年东营市南部土壤干旱情况持续北移，表明该地区土壤含水量对黄河水补给的依赖加重。对土壤含水量反演结果进行分析发现，年和 年的反演误差与 年相比有所增大。原因是，采用的 影像分辨率为 ，数据的分辨率为 ，数据差异较大且分辨率不高。综上分析，以后应基于更高分辨率的数据进行研究，以增强半经验模型的普适性。此外，对比不同年份土壤含水量反演结果，发现使用该耦合模型时决定系数不低于，说明反演值与实测值之间相关性较强。对 年土壤含水量空间分布情况进行对比发现，土壤含水量在空间上呈现连续性变化，与实际情况相符，说明该耦合模型可以反映该地区土壤含水量的实际情况。基于耦合模型反演的土壤含水量具有实际应用价值，能够在一定程度上为农业播种、浇灌等提供依据。结论以山东省东营市为研究区，基于光学遥感数据与 协同反演植被覆盖区土壤含水量。首先利用光学遥感数据的波段特征修正水云模型中的植被含水量，然后根据 的雷达后向散射系数和雷达入射角，引入 模型消除地表粗糙度对土壤含水量反演精度的影响，最终耦合修正后的水云模型和 模型，达到提高传统遥感反演土壤含水量模型精度的目的，并通过连续年份反演讨论了耦合模型的精度和适用性。耦合模型反演土壤含水量与实测值之间的关系表明：使用植被光谱指数表达植被含水量时，基于 的二次函数表达植被含水量的效果最好；修正的水云模型与 模型构建的耦合模型反演土壤含水量精度较高，且效果稳定，其决定系数在 以上，小于等于，该耦合模型具有较强普适性。参考文献：湛江，李志萍，于小朋黄河下游冲积平原典型区包气带土壤水分特征曲线的模型优选人民黄河，（）：李向龙，赵洪丽，赵红莉，等极限学习机模型的土壤含水量反演研究测绘科学，（）：李伯祥，陈晓勇基于 多源遥感数据的河北省景县农田土壤水分协同反演生态与农村环境学报，（）：张新乐，秦乐乐，郑兴明，等基于主动微波遥感的典型黑土区土壤水分反演 东北农业大学学报，（）：谢媛，严颂华，马秦生，等融合合成孔径雷达和光学影像的植被覆盖区土壤湿度反演科学技术与工程，（）：张鹏岩，耿文亮，杨丹，等黄河下游地区土地利用和生态系统服务价值的时空演变农业工程学报，（）：徐辉，丁祖栋，武玲玲黄河下游沿黄城市生态系统健康评价人民黄河，（）：，黄志蕾，朱纹君，齐善忠近 年黄河下游城市后备可利用土地资源时空演变分析：以济南市为例水土保持研究，（）：屈凤莲，姚西文，毛泽秦淤地坝对沟坡土壤水分及作物生长影响的研究人民黄河，（）：，（）：，（）：郭华东，王为民，朱博勤多模态传感器辐射模态数据在 飞船模拟验证研究 全国国土资源与环境遥感技术应用交流会论文文集 武汉：中国地质学会，：郭文，马梦梦，孙培彦基于 数据和 神经网络的裸露地表土壤含水量反演研究中国农村水利水电，（）：王睿馨，宋小宁，马建威，等基于 全极化数据的张掖地区土壤水分的反演中国科学院大学学报，（）：沙莎，胡蝶，王丽娟，等基于 数据的黄土高原植被含水量的估算模型研究遥感技术与应用，（）：，（）：【责任编辑 吕艳梅】人 民 黄 河 年第 期