基于无人机高光谱棉田土壤含水量反演研究.pdf

中 国 煤 炭 地 质COAL GEOLOGY OF CHINA第 35 卷 8 期2023 年 8 月Vol.35 No.8Aug.2023基于无人机高光谱棉田土壤含水量反演研究许将，徐凯磊，翟铄，藏宇晨（中煤航测遥感集团有限公司，陕西西安710199）摘要：以兵团农六师共青团农场示范园棉田为研究区域，采集无人机高光谱数据、土壤含水量及棉花叶片光谱数据，分析土壤含水量与叶片光谱特征相关性，构建地面高光谱与无人机高光谱土壤含水量反演模型。基于两种数据源相同采样点光谱曲线特征对比分析结果，对无人机高光谱土壤含水量反演模型进行优化校正，实现地面研究成果与航天/空遥感数据结合，为无人机高光谱遥感技术在精准农业中的应用提供思路和方法。关键词：土壤含水量；遥感；无人机；高光谱；模型优化中图分类号：S153.6文献标识码：AInversion of Soil Moisture in Cotton Fields Based on UAV Hyperspectral Remote SensingXU Jiang，XU Kailei，ZHAI Shuo，ZANG Yuchen（China Coal Aerial Survey and Remote Sensing Group Co.，Ltd.，Xi an，Shaanxi710199）Abstract:Taking the cotton field in the demonstration garden of the communist youth league farm of the 6th agricultural division of thecorps as the research area,the hyperspectral data of UAV,soil moisture and cotton leaf spectral data were collected,the correlationbetween soil moisture and leaf spectral characteristics was analyzed,and the inversion model of measured hyperspectral and UAVhyperspectral soil moisture was constructed.Based on result of comparative analysis of spectral features of the same sampling pointbetween measured hyperspectral and UAV hyperspectral data,the UAV high-spectral soil moisture content inversion model wasoptimized,that realized ground research results and aerospace/aviation remote sensing data combination and provided ideas andmethods for the application of UAV high spectral remote sensing technology in precision agriculture.Keywords:soil moisture content;remote sensing;UAV;hyperspectral;model optimization0引言土壤水分含量是描述地气能量变换和水循环的重要因素，也是研究地表植被生长发育条件的关键变量1。由于高光谱影像数据难以获取，目前土壤水分高光谱遥感监测方法多数是建立实测土壤光谱反射率及其数学变换形式与土壤水分数学关系进行定量反演监测2-5。随着无人机平台及无人机载高光谱成像设备的研发，解决了高光谱影像数据获取的问题。由于影像数据获取的是地表的信息，多数存在植被覆盖土壤信息的情况，而植被覆盖下的土壤水分遥感监测一直是遥感技术研究的热点和难点。刘万侠、周鹏等均基于雷达数据，利用水-云模型去除植被影响，建立不同极化后向散射系数与土壤含水量之间的关系，开展各自研究区土壤水分遥感监测研究工作6-7；詹志明等分析土壤水分在 ETM的 Nir-Red 光谱特征空间的分布规律，提出了基于Nir-Red光谱特征空间的土壤水分监测模型8，吴春雷等则是在光谱特征分析基础上，提 出 了 PDI 改 进 指 数 植 被 调 整 垂 直 干 旱 指 数（VAPDI）对研究区进行土壤水分监测1；此外还有许多专家学者利用植被指数和温度干旱植被指数开展植被覆盖区土壤水分遥感监测研究。综上所述，目前利用植被光谱特征开展土壤水分高光谱遥感监测研究工作尚少，而将地面研究成果与航天/空遥感数据相结合，实现遥感监测天-地一体化的研究工作就更少。本文利用大疆 M600PRO 无人机搭载 Gaiasky-mini 高光谱成像采集系统，获取兵团农六师共青团农场现代农业园区示范园整块棉田高光谱遥感影像数据，并同步开展地面采样测试工作，分析土壤含水量与棉花叶片光谱特征关系，构建地面实测光谱曲线土壤水分监测模型，开展无人机高光谱数据模型优化研究工作，为基金项目：陕西省重点研发计划（2022ZDLSF06-09），国家自然科学基金项目（62071184）第一作者简介：许将（1991），男，硕士研究生，工程师，主要研究方向为遥感应用。E-mail：收稿日期：2022-07-10责任编辑：张燕生doi:10.3969/j.issn.1674-1803.2023.08.11文章编号：1674-1803（2023）08-0064-068 期许 将，等：基于无人机高光谱棉田土壤含水量反演研究无人机高光谱遥感技术在精准农业中的应用提供思路和方法。1材料与方法1.1研究区概况共青团农场地处天山北麓，准噶尔盆地南缘，昌吉州境内，三屯河流域下游。农场地处欧亚大陆腹地，属温带半荒漠气候，干旱少雨，无霜期短，日照充足，昼夜温差大，光热资源丰富。现代农业园区位于共青团农场东片西部，甘莫公路两侧，南北长7.50km，东西宽6.60km，面积45.4 km2（图1）。图1 研究区地理位置Figure 1 Geographical position of study area1.2高光谱数据采集及预处理1.2.1无人机高光谱数据采集2022年9月利用无人机搭载高光谱成像仪，获取了兵团农六师共青团农场现代农业园区示范园整块棉田无人机高光谱影像数据，飞行期间天气晴朗无云，飞行时间为中午13：0015：00，航向重叠度60%，旁向重叠度40%。1.2.2地面高光谱数据采集土壤和叶片反射率测定采用ASD FieldSpec4光谱仪在室内进行采集，光谱采集时，将样品放置于黑色平板上，采用8视场角探头，探头置于样品表面20cm的垂直上方，先对白板进行校正，获取绝对反射率，并在保证目标与参考白板之间对应相同的光照条件下，每个样品采集5次，取平均值作为该样品的光谱曲线。1.2.3光谱数据预处理为了消除光谱曲线噪音可能引起的误差，利用ENVI软件自带Boxcar Average进行光谱平滑处理。为了减少物理环境因素的影响，提高光谱灵敏度，对平滑后的光谱进行一阶微分、倒数对数一阶微分变换，用于遥感分析建模。1.3模型构建与检验偏最小二乘回归方法（PLSR）是Wold和Albano于 1983 年提出的一种新型的多元统计分析方法。它集多元线性回归分析、主成分分析和典型相关分析的基本功能为一体，可以较好地解决以往用普通多元线性回归分析方法无法解决的问题9。特别是当自变量集合内部存在较高程度的相关性时，其结论比普通多元回归更加可靠。它通过提取概括原数据信息的综合变量（新成分），建立新变量与因变量的回归关系，最后再表达成原变量的回归方程。PLSR算法建模的思想是给出一对数据矩阵X和Y，X为输入自变量矩阵，Y为输出因变量矩阵10，形式如下：X=x11x1pxn1xnp，Y=y11y1qyn1y1nq（1）其中：n是样本数；p和q分别代表自变量和因变量的个数。该算法大体可以分为以下3个步骤：1）将数据矩阵X和Y标准化处理得到E0和F0。E0为矩阵X经过标准化后的矩阵，F0为矩阵Y经过标准化后的矩阵。x*ij=xij-x jSxj，i=1，2，n，j=1，2，py*ij=yij-y jSyj，i=1，2，n，j=1，2，q（2）式（2）为标准化处理公式，其中x j是xj的均值；Sxj是xj的标准差；y j是yj的均值；Syj是yj的标准差；x*ij是xij经过标准化处理后的值，y*ij是yij经过65中国煤炭地质第 35 卷标准化处理后的值，标准化矩阵形式如式（3）所示。E0=x*11x*1px*n1x*np，F0=y*11y*1qy*n1y*nq（3）2）提取E0和F0的主成分t1和u1，式（4）为主成分的基本形式和计算方法，其中w1和c1分别为E0和F0的第一个轴，此处要求t1和u1能够很好地代表数据矩阵X和Y中的变异信息，式（5）为其约束条件，其中Cov（t1，u1）、Var（t1）、Var（u1）、r（t1，u1）分别代表t1和u1的协方差、方差和相关度，s.t.表示约定条件。推导可得w1是矩阵ET0F0FT0E0最大特征值对应的特征向量，c1是矩阵FT0E0ET0F0最大特征值对应的特征向量。t1=w11x*1+w1px*p=E0w1u1=c11y*1+c1py*p=F0c1（4）Cov()t1，u1=Var()t1Var()u1r()t1，u1maxmaxw1，c1E0w1，F0c1s.t.wT1w1=1，cT1c1=1（5）3）构建回归方程。如若达到精度，算法结束；否则就利用X和Y被t1解释后的残余信息进行第二轮成分提取，循环至方程精度最高时为止。假设最终提取了m个成分t1，t2，tm，则PLSR就是通过施行yk对t1，t2，tm的回归，其中k=1，2，q然后再表达成yk关于原自变量x1，x2，xp的回归方程。将 25 个样本随机选取 18 个样本用来构建模型，7 个用来验证模型。模型精度通过决定系数（R2）、均方根误差（RMSE）两个参数进行检验，模型的决定系数介于 01，决定系数越大，模型精度越高，均方根误差越小，模型预测能力越强，精度越高。2结果分析2.1地面棉花叶片光谱特征及其与土壤含水量与相关性分析土壤含水量与叶片反射率相关性曲线在550nm及 750nm 附近呈反射峰，相关系数绝对最大值为0.41；叶片反射率经过不同形式数学变换后，与土壤含水量的相关性有了一定提高，两种变换结果与土壤含水量的相关性曲线形态镜像相似，存在两个负相关和两个正相关波段区间，其中，土壤含水量与反射率一阶微分相关系数绝对最大值为0.69，与反射率倒数对数之一阶微分相关系数绝对最大值为0.68。植被指数是植被光谱特征重要表现形式，利用植被光谱特征进行波段组合，将计算获得的不同植被指数与土壤含水量进行相关性分析，表明土壤含水量与常规植被指数存在明显的正相关性。且与 500670nm 范围内反射率变化趋势存在更高的相关性，在该波段范围内构建的比值植被指数与土壤含水量相关系数达0.80。2.2地面高光谱数据土壤含水量反演模型构建和检验依据土壤含水量与植被光谱特征相关性分析结果，选择反射率变换形式与土壤含水量相关系数绝对最大值对应的波段、常规归一化植被指数（NDVI）、修正的土壤调节植被指数（OSAVI）、归一化 植 被 指 数（NDVI（620670）、比 值 植 被 指 数（SR（620/670）构建土壤水分线性回归模型（表 1），利用验证样本对模型进行检验，如表2、图2所示，新构建植被指数构建的模型精度优于常规植被指数及反射率变换形式所构建的模型；所有模型验证结果显示验证样本接近 1 1 线，个别验证样品偏离1 1线。表1 地面高光谱数据土壤含水量回归模型Table 1 Regression model of soil moisture content from ground hyperspectral data光谱指标NDVIOSAVISR（620670）NDVI（620670）PLSR模型Y=0.30NDVI-0.08Y=0.25OSAVI-0.03Y=0.23*SR（620670）-0.16Y=0.55NDVI（620670）+0.06Y=-0.006-22.32R967+11.81(Log(1/R971)+0.06NDVI+0.05OSAVI+0.04SR（620670）+0.10NDVI（620670）建模集R20.510.550.620.630.68RMSE0.013 40.012 90.011 90.011 70.010 1验证集R20.790.670.680.690.76RMSE0.010 20.011 80.011 70.012 90.008 1注：SR（620670）=B620/B670，NDVI（620670）=（B620-B670）/（B620+B670），其中B620、B670代表波长620nm和670nm处反射率，下同。由于光谱数据各波段以及波段组合存在高度相关性，且建模样本点相对较少，故将土壤含水量作为因变量，选择线性回归模型中的光谱指标作为自变量，利用SIMCA-P软件构建土壤含水量的偏最668 期许 将，等：基于无人机高光谱棉田土壤含水量反演研究小二乘回归模型11，结合表2、图2，偏最小二乘回归模型建模集和验证集的决定系数和均方根误差指标均优于线性回归模型，验证样本与1 1线接近程度最高。图2地面高光谱土壤含水量回归模型验证结果Figure 2 Results of soil moisture regression model validation from ground hyperspectral data2.3基于无人机高光谱数据土壤含水量反演模型构建和检验采集25个调查点对应的无人机高光谱影像上光谱曲线数据，兼顾模型优化校正目的，计算光谱反射率及其数学变换结果以及NDVI、OSAVI、NDVI（620/670）、SR（620670）植被指数，分析上述光谱指标与土壤含水量相关性，结果表明：土壤含水量与高光谱影像数据植被指数呈正相关，与常规植被指数OSAVI 的相关性达 0.68，与新构建比值植被指数SR（620/670）达 0.72，这与地面光谱数据分析结果较一致；而高光谱影像数据光谱反射率及其数学变换结果与土壤含水量相关性分析结果与地面实测光谱数据分析结果存在较大差异，整体相关性均有所降低，且相关性最高对应的波段位置较实测光谱数据发生变化，说明两种数据源在光谱特征细节上存在差异。依据高光谱影像数据与土壤含水量相关性分析结果，选择植被指数构建无人机高光谱数据土壤含水量线性回归、偏最小二乘模型，并利用验证样本对模型进行检验。验证结果显示，常规植被指数 NDVI模型验证样本偏离 1 1线程度最高，其他模型验证样本偏离 1 1线程度接近。结合表1、表2和图2、图3，认为无人机高光谱数据构建土壤含水量模型精度较地面高光谱数据构建模型精度差。2.4无人机高光谱数据土壤含水量模型优化校正为了实现地面研究成果与航天/空遥感数据相结合，提高无人机高光谱模型精度，利用地面光谱数据对无人机高光谱数据模型成果进行优化校正。由于两种数据源构建的植被指数存在高度相关性（图4、图5），其中，常规植被指数NDVI和OSAVI相关系数均为0.90，而新构建的NDVI和SR相关系数均为0.86，建立两种数据源同一植被指数之间的数学关系，将无人机高光谱数据土壤含水量模型（不包括偏最小二乘模型）进行优化校正。将土壤含水量实测值与各模型反演值作对比分析（图6），各模型土壤含水量反演值与实测值分布趋势整体一致，个别点存在明显差异；地面实测高光模型反演数据与无人机高光谱模型反演数据分布趋势一致性高于各自反演值与实测值的趋势一致性，而无人机高光模型优化校正后，提高了与地面高光谱模型反演值的一致性。利用验证集对优化后无人机高光谱土壤含水量模型进行验证，结果显示NDVI指数构建模型的R2由0.55提升到0.79，RMSE由0.015 2提升到0.011 7；OSAVI指数构建模型的R2由0.63提升到 0.67，RMSE 由 0.013 8 提升到 0.012 8；SR（620/670）指数构建模型的R2由 0.62 提升到 0.68，RMSE由0.013 5提升到0.011 8；NDVI（620/670）指数构建模型的R2由0.62提升到0.69，RMSE由0.014 8提升67中国煤炭地质第 35 卷到0.011 8，说明经过模型优化校正后的无人机高光谱土壤含水量反演精度得到提高。表2 无人机高光谱数据土壤含水量回归模型Table 2 Regression model of soil moisture content from UAV hyperspectral data光谱指标NDVIOSAVISR（620670）NDVI（620670）PLSR模型Y=0.39NDVI-0.18Y=0.36OSAVI-0.21Y=0.50SR（620670）-0.49Y=1.89NDVI（620670）-0.07Y=-0.36-0.25NDVI+0.23OSAVI+0.30SR（620670）+0.96NDVI（620670）建模集Y20.420.430.570.530.62RMSE0.014 60.014 40.012 50.013 10.011 8验证集R20.550.630.620.620.63RMSE0.015 20.013 80.013 50.014 80.013 0图3 无人机高光谱土壤含水量回归模型验证结果Figure 3 Results of soil moisture regression model validation from UAV hyperspectral data图4 棉田两种数据源同一调查点NDVI值分布对比趋势图Figure 4 Comparison trend chart of NDVI valuedistribution at the same survey point from two datasources in cotton fields图5 棉田两种数据源同一调查点SR值分布对比趋势图Figure 5 Comparison trend chart of SR valuedistribution at the same survey point from twodata sources in cotton fields688 期许 将，等：基于无人机高光谱棉田土壤含水量反演研究图6 土壤含水量实测值与各模型反演值对比分析Figure 6 Comparison between the measured data and and inversion data of different model of soil moisture3结论本文基于地面高光谱数据和无人机高光谱数据，分析棉花叶片光谱特征与土壤含水量相关性，采用偏最小二乘模型法构建土壤含水量反演模型，并实现无人机高光谱土壤含水量反演模型校正优化研究，得出如下结论：1）综合分析不同土壤含水量棉花叶片光谱特征并计算相应植被指数，发现植被指数与土壤含水量呈明显正相关，且土壤含水量与500670nm范围内新构建植被指数存在更高相关性。2）基于偏最小二乘模型构建无人机高光谱和地面光谱土壤含水量反演模型，结果显示地面高光谱土壤含水量反演模型反演精度优于无人机高光谱土壤含水量反演模型反演精度。3）利用地面高光谱数据对无人机高光谱数据模型进行优化校正。经验证，校正后无人机高光谱土壤含水量估算精度得到有效提高，弥补单一数据源在土壤含水量反演方面的不足。参考文献：1吴春雷，秦其明，李梅，等.基于光谱特征空间的农田植被区土壤湿度遥感监测 J.农业工程学报，2014，30（16）：106-112.2刘伟东，BARET F，张兵，等.高光谱遥感土壤湿度信息提取研究 J.土壤学报，2004，41（5）：700-706.3魏娜.土壤含水量高光谱遥感监测方法研究 M.北京：中国农业科学院，2009.4张锐，李兆富，辛强，等.原状水稻土土壤湿度高光谱特征及反演模型 J.土壤通报，2017，48（1）：32-38.5蔡亮红，丁建丽.基于高光谱多尺度分解的土壤含水量反演J.激光与光电子学进展，2018，55（1）：406-415.6刘万侠，刘旭拢，王娟，等.华南农作物覆盖区土壤水分ENVISAT-ASAR与MODIS数据联合反演算法研究 J.干旱区农业研究，2008，26（3）：39-43.7周鹏，丁建丽，王飞，等.植被覆盖地表土壤水分遥感反演 J.遥感学报，2010，14（5）：959-973.8詹志明，秦其明，阿布都瓦斯提 吾拉木，等.基于NIR-Red光谱特征空间的土壤水分监测新方法 J.地球科学，2006，36（11）：1020-1026.9周倩倩，丁建丽，唐梦迎，等.干旱区典型绿洲土壤有机质的反演及影响因素研究 J.土壤学报，2017，55（2）：313-324.10李萍，赵庚新，高明秀，等.黄河三角洲土壤含水量状况的高光谱估测与遥感反演 J.土壤学报，2015，52（6）：1262-1270.11李相，丁建丽.基于实测高光谱指数与HIS影像指数的土壤含水量监测 J.农业工程学报，2015，31（19）：68-75.69