您当前的位置:首页 > 论文范文 > 建筑论文 > 水利论文 > 正文

去除土壤背景的冬麦根土壤含水率遥感反演方法

  摘要:为剔除无人机多光谱图像中的土壤背景、提高作物根域土壤含水率反演精度,以不同水分处理的拔节期冬小麦为研究对象,利用无人机多光谱相机分别在09:00、11:00、13:00、15:00和17:00等5个时刻获取高分辨率多光谱图像,采用改进的植被指数阈值法快速确定植被像元与土壤像元的分类阈值,通过阈值划分剔除土壤背景,并根据阈值变化研究土壤背景对冬小麦冠层反射率的影响,建立了剔除土壤背景前后基于植被指数的土壤含水率反演模型。结果表明,应用改进的植被指数阈值法可有效剔除多光谱图像中的土壤背景,其中基于植被指数RDVI的剔除精度最高,总体精度在91.32%以上;土壤背景对冬小麦冠层近红外波段的反射率影响较大,红边波段次之,而对可见光波段的反射率影响较小;剔除土壤背景前后的植被指数与土壤含水率均呈线性关系,剔除土壤背景对反演土壤含水率的精度有显著提高,其中NGRDI反演深度10~20cm的冬小麦根域土壤含水率效果最好,建模集R2和RMSE分别为0.739和2.0%,验证集R2和RMSE分别为0.787和2.1%。

  关键词:冬小麦;土壤含水率;土壤背景;多光谱图像;植被指数

  0引言

  土壤水分是作物生长发育的重要条件,也是指导灌溉的理论依据。土壤水分包括表层土壤水分(保持在10cm以上的土壤中)和根域土壤水分(保持在10cm以下的土壤中)[1]。表层土壤水分是解决地面与大气之间水分和能量交换的关键组成部分,而作物生长发育的水分主要来源于根域土壤[2]。土壤含水率测定方法有干燥法、土壤湿度传感器测量法、土壤水量平衡法和遥感技术测量法等。传统的干燥法一直被用作土壤含水率测定的标准技术方法,但是该方法相对耗时、耗力,且具有破坏性;土壤湿度传感器测量法和土壤水量平衡法则存在设备成本高、准确性差等缺点[3];而无人机遥感技术具有操作灵活、时空分辨率较高等优点,为快速无损测定土壤含水率提供了可能[4-5]。

  早期利用热红外遥感与现场测量相结合的方法测量或定量估算土壤含水率,并且将光学、热红外遥感技术和微波传感器相结合,同样为大规模监测土壤含水率提供了可能[6-9]。尽管目前已取得了一些研究成果[10-11],但主要是反演表层土壤含水率,利用高分辨率的无人机多光谱遥感数据对作物覆盖下根域土壤含水率的研究很少。

  研究表明,根域土壤含水率对作物生理状态影响很大,理论上可以通过作物冠层反射率或植被指数间接测定根域土壤含水率[12-14]。由于高分辨率无人机多光谱图像受到土壤背景的干扰,导致根域土壤含水率的监测精度普遍较低。因此,有学者对降低土壤背景的干扰进行了研究[15-17]。尽管在剔除土壤背景研究上有了很大的改进,但土壤线和各种指数仍然对土壤背景和覆盖度变化有一定的敏感性[18-19]。采用植被指数可以在最大程度上剔除土壤背景,而植被指数阈值法是基于植被指数有效识别作物与非作物像元之间的阈值,通过阈值划分来剔除土壤背景,而应用植被指数阈值法剔除土壤背景的关键是确定土壤与作物之间的植被指数阈值。本文以4种水分处理下的拔节期冬小麦为研究对象,利用无人机搭载多光谱传感器获取高分辨率多光谱图像,采用改进的植被指数阈值法快速确定多光谱图像中土壤与冬小麦的植被指数阈值,以实现土壤背景的剔除,研究不同土壤含水率和作物阴影对土壤背景剔除精度的影响,以及土壤背景对提取冬小麦冠层反射率的影响,利用剔除土壤背景和未剔除土壤背景的植被指数与土壤含水率建立关系模型,探究剔除土壤背景对冬小麦根域土壤含水率反演精度的影响。

  1材料与方法

  1.1试验地概况

  试验地位于陕西省杨凌示范区西北农林科技大学的中国旱区节水农业研究院(34°17'42.17″N,108°4'20″E,海拔525m)。该地区年平均降雨量635mm,且主要集中在7—9月,属于暖温带季风半湿润气候区,试验地深度0~60cm的土壤类型为中壤土,平均田间持水率为26%,凋萎系数为8.6%(均为质量含水率),土壤干容重1.44g/cm3,土壤pH值为8.1,有机质质量比为13.3g/kg,总氮质量比为0.82g/kg,试验地概况如图1所示。

  1.2试验设计

  试验小区布设采用完全随机设计,如图1b所示,设置4个水分处理,灌水上限分别为T1(50%田间持水量)、T2(65%田间持水量)、T3(80%田间持水量)和T4(对照组,95%田间持水量),每个水分处理3个重复,如小区编号T1-1代表T1水分处理的第1个重复,共计12个小区。试验地播种前施氮肥和磷肥各240kg/hm2作为基肥,冬小麦品种选用的是陕西中部地区推广的小偃22,播种时间为2018年10月15日,每个小区设置为4m×4m,种植行距为0.25m,播种密度为100万株/hm2。本研究于2019年3月24、25、28、31日进行试验,24日为灌水处理前的对照试验,25、28、31日为灌水处理后的试验。此时冬小麦正处于拔节期,天气晴朗,而且试验期间对每个小区都配有移动式遮雨棚,防止降雨干扰。

  1.3无人机多光谱遥感数据采集

  试验期间,使用大疆创新公司(DJI)生产的经纬M600型六旋翼无人机搭载多光谱传感器获取冬小麦冠层多光谱图像。如图2所示,多光谱传感器为Tetracam公司生产的Micro-MCA型多光谱相机(简称μ-MCA)。μ-MCA多光谱相机有490、550、680、720、800、900nm共6个波长的光谱采集通道,能够搭载于M600型六旋翼无人机。试验期间每日采集5次多光谱影像(采集时间为09:00、11:00、13:00、15:00和17:00),无人机飞行高度为15m,多光谱相机镜头垂直向下,对应的地面分辨率为8cm,并利用放置在试验田间的白色参考板(100%反射参考目标)进行影像的辐射校正。

  0.4地面数据采集与处理

  试验期间,在每日13:00光谱数据采集完成后,及时在每个水分处理小区的中心区域打钻取土,根据拔节期冬小麦根系活动层所在深度,每个小区取土深度分别为0~10cm、10~20cm、20~30cm和30~40cm,共4个深度的土壤样本,并采用干燥法测定4个深度下的冬小麦土壤含水率(体积含水率,%),每个深度土壤样本共48个。随机抽取3/4数据用于建模,1/4数据用于验证,对4个深度土壤含水率的样本进行统计分析,样本统计特征如表1所示。在拔节期冬小麦试验结束后(2019年3月31日)对每个小区中间区域的冬小麦取样,并用干燥法测定冬小麦的植株含水率(相对含水率,%)。

  1.5数据分析

  1.5.1土壤背景剔除

  基于无人机多光谱图像的波段范围,选择植被指数DWSI[20]、NGRDI[21]、NDVI[22]、EXG[23]、OSAVI[24]、RDVI[25]和MSR[26]来剔除土壤背景并用于反演土壤含水率。使用PixelWrench2对μ-MCA多光谱相机采集的原始图像进行预处理,并用软件ENVI5.3对预处理后的各个小区的多光谱图像进行波段运算,得到多种植被指数的灰度图像。为实现准确、自动的植被检测,以及正确的分割参数,需要一种简便、高效的方法来寻找设置植被与裸地之间断点的植被指数阈值,而植被指数阈值法[27]可以有效确定作物与土壤像元之间的阈值。本文利用改进的植被指数阈值法将OTSU算法和高斯滤波结合,假设图像包含两类像素(土壤和植被),然后基于最小化组合传播(类内方差)和高斯滤波线性减噪过程,可以快速确定多光谱图像中作物与土壤像元之间的阈值。在Python3.6.3中运行该算法得到小区的各种植被指数阈值,在ENVI软件中通过植被指数阈值对各个小区多光谱图像的冬小麦和土壤进行掩膜处理,实现土壤背景的剔除。

  1.5.2土壤含水率反演及精度评价

  利用混淆矩阵精度验证法评价植被指数阈值法的分类精度,评价指标为Kappa系数和总体精度[28]。用分类精度最高的植被指数剔除土壤,并通过该植被指数确定的阈值来提取冬小麦冠层反射率,分析土壤背景中土壤含水率变化和阴影干扰对反射率的影响,进一步判别植被指数用于剔除土壤背景的准确性。

  用未剔除土壤和剔除土壤的反射率计算不同植被指数,通过这些植被指数反演不同深度的土壤含水率,以7个植被指数为自变量和土壤含水率为因变量建立线性模型,比较剔除土壤背景前后模型的决定系数(R2)和均方根误差(RMSE),分析土壤背景对反演土壤含水率的影响。R2越接近1,RMSE越接近0,则该反演模型的精度和准确性越好[11,29]。

  2结果与分析

  2.1土壤背景对植被指数阈值法分类精度的影响

  基于1.5节描述的数据处理过程,应用7种植被指数处理不同时刻(09:00、11:00、13:00、15:00和17:00)获取的冬小麦多光谱图像,分析作物阴影对植被指数阈值法剔除背景的影响。通过7个植被指数剔除背景后,进行精度评价得到的Kappa系数及图像分类精度如图3所示。

  由图3可看出,利用植被指数阈值法进行分类的拔节期冬小麦多光谱图像效果较好,总体精度和Kappa系数整体呈正态分布。其中RDVI对冬小麦多光谱图像的分类精度最高,在所有图像的分类结果中总体精度达到了91.35%,大部分都集中于97.25%,而且Kappa系数均在0.83以上,均值为0.91。而基于植被指数NDVI和OSAVI的分类效果也相对较好,Kappa系数整体都大于0.7,总体精度也都在88.59%以上。但是基于EXG、NRGDI和DWSI这3种植被指数分类效果相对较差,总体精度都出现80%甚至低于80%的情况,而且对应的Kappa系数出现了低于0.6的情况。由于太阳高度角的变化,导致09:00、11:00、15:00和17:00采集多光谱数据时的太阳没有完全正射冬小麦冠层,造成土壤背景中不只有裸土,还存在较多冬小麦投影到土壤中的阴影,这使得土壤存在阳光直射的裸土和有阴影的土壤。所以,为避免作物阴影的干扰,对24、25、28、31日13:00(此时阳光正射试验地冬小麦冠层)的冬小麦多光谱图像进行分类,分析土壤含水率变化对植被指数与阈值法剔除背景的影响,得到的精度评价结果如图4所示。

  从图4可以发现,RDVI的分类效果依然最好,总体精度达到91.32%以上,Kappa系数都在0.85以上。NDVI和OSAVI的分类精度相似,仅比RDVI略低一点。DWSI、NGRDI和EXG的分类效果依然较差,但是基于植被指数NGRDI在土壤含水率变化下的图像分类效果明显高于作物阴影干扰下的图像分类效果,总体精度普遍都在96.23%左右,但是也会出现分类精度降低的情况。

  造成这种现象的原因可能是计算DWSI、EXG和NGRDI这3个植被指数的中心波长(490、550、680nm)都在可见光波段,冬小麦和土壤背景之间的反射率差别较小,而在近红外波段(中心波长为720、800、900nm)二者之间的反射率差别较大[20,30]。NDVI和OSAVI基于冬小麦和土壤在红色、近红外波段反射率的变化差异,通过归一化扩大土壤与冬小麦差异,同时最小化了光照强度变化的影响,增强了对冬小麦的识别能力,但它们对土壤亮度敏感。而RDVI有效消除了土壤反射率的影响,更好地适应于不同水分处理导致覆盖度差异悬殊的试验区域,使得对冬小麦的识别能力更强[31]。因此含有近红外波段的植被指数在一定程度上更能响应冬小麦生长环境的变化,使对应的分类精度更高。

  2.2土壤背景对提取冬小麦冠层反射率的影响

  土壤背景中土壤含水率变化和裸土上的冬小麦阴影会导致土壤亮度发生变化,同时也对提取的冬小麦冠层反射率有很大影响。

  2.2.1不同土壤含水率对提取冬小麦冠层反射率的影响

  尽管植被指数RDVI的分类精度最好,但是同一类型的土壤在不同含水率下的光谱反射率变化差异较大,这导致了用植被指数阈值法确定的阈值受土壤含水率变化的影响较大,这也会影响提取冠层反射率。为避免阴影和土壤水分变化的交互影响,本研究以13:00获取的冬小麦多光谱图像为研究对象,图5给出了试验期间植被指数阈值法确定13:00多光谱图像的RDVI阈值与不同深度土壤含水率的相关性矩阵图。

  图5相关性矩阵图表明,在水分处理前后基于RDVI确定的阈值和土壤含水率均呈正相关,并与0~10cm的土壤含水率相关性均最高,而与10cm以下的土壤含水率相关性普遍较低。其中水分处理前的24日阈值与深度0~10cm的土壤含水率相关性最高(r=0.914),而水分处理后的25、28、31日阈值与深度0~10cm土壤含水率的相关性呈减小趋势,深度30~40cm的土壤含水率与RDVI阈值的相关性表现出同样的变化趋势,10~20cm和20~30cm深度的土壤含水率与阈值之间的相关性变化不大。这表明0~10cm处土壤含水率对植被指数阈值法确定作物与土壤像元之间的阈值影响最大,而水分处理导致土壤含水率与阈值之间的相关性降低,这对提取冬小麦冠层反射率影响较大。

  图6所示(图中b1~b6分别表示中心波长为490、550、680、720、800、900nm的光谱反射率,下同)为水分处理前后基于RDVI阈值提取的冬小麦冠层反射率。T1处理的冬小麦由于土壤含水率没有变化,因此冠层波段反射率变化很微弱。

作者:张智韬,周永财,杨帅,谭丞轩,劳聪聪,许崇豪 单位:西北农林科技大学旱区农业水土工程教育部重点实验室 西北农林科技大学水利与建筑工程学院