芮婷婷，徐云飞，程 琦，杨 斌，冯志军，周 涛，张世文

(1.安徽理工大学空间信息与测绘工程学院，安徽淮南 232001;2.安徽理工大学矿山采动灾害空天地协同监测与预警安徽普通高校重点实验室，安徽淮南 232001;3.安徽理工大学矿区环境与灾害协同监测煤炭行业工程研究中心，安徽淮南 232001；4.安徽理工大学地球与环境学院,安徽淮南 232001)

冬小麦在我国是人们日常生活中主要的粮食作物，快速、准确地获取冬小麦水分含量对农田灌溉管理、旱情状况监测、作物长势监测等具有重要的科学意义。含水量影响着作物的生理特征和健康状况，在缺水时作物叶片颜色、纹理特征和形态结构等会发生变化，如叶片出现蜷缩、枯黄和色斑等，进而影响作物的生长和产量。作物水分含量可以由根、枝、叶、冠层的生理特性来反映，其中作物叶片的新陈代谢最旺盛，生理特性最敏感，可以通过监测作物叶片水分含量来评估作物整体水分状况。

近年来，国内外学者基于遥感技术对作物水分含量的快速监测进行了广泛探讨和深入研究。遥感技术主要分为卫星遥感、地物光谱仪和无人机遥感。大多数研究结合卫星遥感影像监测作物水分含量，但是卫星运行时会受到云层的影响，使得利用卫星遥感监测精度受限。利用地物光谱仪可以获取作物的冠层光谱信息，但所得信息是点状的，难以用于区域性水分含量的动态监测。而无人机遥感平台具有移动性强、运营成本低、获取影像分辨率高、作业周期短等优点，可快速地获取农田高精度影像，更好地满足农田尺度上精准监测的要求。胡珍珠等筛选了与叶片含水量相关的光谱水分指数，构建了核桃果实不同生育时期的叶片含水量反演模型。吾木提·艾山江等通过灰色关联度分析，选择与叶片含水量关联度较高的5种植被指数，构建了春小麦叶片含水量的偏最小二乘模型和bp神经网络反演模型。彭要奇等采用改进的神经网络模型对玉米叶片含水量进行反演。陈秀青等在叶片和冠层两个尺度上，采用两波段植被指数，构建冬小麦叶片含水量的偏最小二乘回归和竞争自适应重加权采样-偏最小二乘回归模型。wang等分析了不同水分胁迫条件下春小麦冠层反射率与冠层含水量的关系。chen等利用无人机数据，构建了植被补给水指数，筛选出棉花不同部位的敏感波段。以上研究大多基于单一的敏感波段或高光谱植被指数监测作物含水量，而将无人机多光谱数据与敏感波段、植被指数结合反演作物叶片含水量的研究尚未见报道。

本研究利用无人机获得多光谱影像，同步采集了冬小麦叶片含水量数据，分析了叶片含水量与光谱反射率和植被指数之间的相关性，进而筛选出与叶片含水量相关的敏感波段组及植被指数组，再利用得到的敏感波段组和植被指数组分别构建了冬小麦叶片含水量的pls、elm和pso-elm预测模型，通过选定模型评价指标比较3种模型的预测精度，最后将最优反演模型应用于研究区，获得研究区冬小麦叶片含水量的空间分布图，以期为作物叶片水分含量的快速监测和田间管理提供科学依据和理论支持。

研究区位于安徽省淮北市杜集区和烈山区的交界区域(116.85°e，33.97°n)，地处中纬度地区，属于暖温带半湿润季风气候，全年平均温度为17 ℃，平均降水量为849.6 mm，日照时数为 4 430.2 h，雨水多集中于6-8月，土壤类型为壤土。研究区采取冬小麦和夏玉米轮作的种植方式，灌溉方式为雨养。供试小麦品种为淮麦22，小麦实行机播，行距15 cm，播种密度约为每公顷100万株。研究区位置和采样点分布情况如图1所示。

图1 研究区位置和采样点分布图

1.2.1 无人机多光谱数据

本研究使用的无人机为大疆精灵4多光谱版植保无人机，其光谱相机配备6个1/2.9英寸cmos影像传感器，其中1个彩色传感器用于常规可见光(rgb)成像，5个单色传感器用于包含蓝光波段(中心波长450 nm，带宽32 nm)、绿光波段(中心波长560 nm，带宽32 nm)、红光波段(中心波长650 nm，带宽32 nm)、红边波段(中心波长730 nm，带宽32 nm)和近红外波段(中心波长840 nm，带宽52 nm)的多光谱成像。

无人机多光谱影像的采集时间为2021年1月7日，冬小麦处于越冬期，飞行时间为11：00-14：00，飞行时天空晴朗无云，风力较小。起飞前手动控制无人机飞行至校准白板的正上方约1 m处，采用相机单摄模式来拍摄标准白板。本次飞行设置航线为s型，飞行高度为100 m，航向和旁向重叠度分别设置为70%和60%，相机镜头与地面呈90°，拍照模式为等时间间隔。

1.2.2 地面实测数据

根据研究区范围及冬小麦种植面积，利用均匀布点原则布设采样点，并利用gps记录每个采样点的位置信息。冬小麦叶片样品均采用五点法混合采集，然后用密封袋对冬小麦叶片样品进行密封保存并放置于阴凉处，以此来保证叶片水分不受到损失。实验室处理时，冬小麦叶片含水量的测定采用烘干法，将冬小麦植株的叶片和茎分离后称取每一个样品的叶片鲜重，然后将所有样品放入烘箱中，在105 ℃下杀青30 min，再在 80 ℃下烘干至恒重，称取叶片干重。冬小麦叶片含水量(lwc)按公式“lwc=(鲜重-干重)/鲜重×100%”进行计算。

利用pix4dmapper软件对无人机获取的多光谱影像进行拼接处理。拼接后多光谱影像的主要预处理包括几何校正和辐射校正。结合arcmap 10.6软件，以无人机高清数码正射影像为参考影像，分别在5个单波段影像上均匀选取30个参考点进行几何校正，几何校正误差在0.5个像元之内。在多光谱影像中计算标准白板在各波段的dn值，利用公式“=dn/dn×”进行多光谱影像的辐射校正，其中为目标地物的反射率，dn为目标地物的dn均值，dn为标准白板的dn均值，为标准白板的反射率值。

2) 当不确定加箱或减箱事件发生时见图2b)，进入下一阶段t=t+1，更新集装箱序列和船舶贝内箱位序列，减箱事件对应的集装箱c4和箱位p6从序列中被删除，加箱事件对应的集装箱c7被添加到集装箱序列。

将预处理后的无人机多光谱影像及采样点的gps点位信息导入到envi 5.3软件中。以地面实测取样点为中心，在图像上裁剪出200×200(像素)的光谱影像，将感兴趣区域(roi)范围内冬小麦叶片样本的平均反射光谱作为该取样点的光谱反射率，获得不同波段的光谱反射率数据。

作物因其内在生化参数存在差异而表现出不同的光谱反射率。在可见光波段内,健康绿色作物的主要吸收峰在红光波段和蓝光波段的附近形成，主要反射峰在绿光波段形成；而在近红外波段，作物的光谱特征为高反射率和低吸收率；在红边波段，作物的光谱反射率增长较快。将这些特征波段范围内的光谱反射率通过线性或非线性的组合构成植被指数，可用来诊断作物生长状态以及反演各种作物参数。通过以上对植被光谱特征的分析，本研究获取的多光谱影像波段特征，从现有的多光谱植被指数中选取了20个多光谱植被指数用于模型构建(表1)。

表1 植被指数及其相关计算公式table 1 vegetation index and its related calculation formula

偏最小二乘模型是一种具有广泛适用性的线性多元统计数据分析方法，具有主成分分析、典型相关分析和多元线性回归分析的特点，可以有效地解决变量之间存在多重共线性或者样本量较少的问题。因此,其在处理因变量较多或者因变量内部信息高度线性相关，而样本量较少的数据时具有明显的优势，具有简化数据结构、降低数据维度、减少噪声干扰等特点。

极限学习机是一种针对单隐含层前馈神经网络的学习算法。该算法随机产生输入层与隐含层的连接权值和隐含层神经元的阈值，且在训练过程中无需调整，只需要设置隐含层神经元个数，便可以得到全局最优解。与传统的训练方法相比，极限学习机具有学习速度快和训练误差小等优点。elm模型的主要步骤为：①确定样本集；②确定隐含层神经元个数，随机生成输入层与隐含层的连接权值和隐含层神经元阈值；③确定隐含层神经元的激活函数，计算隐含层输出矩阵和输出层权值矩阵，根据权值矩阵和激活函数，最终计算极限学习机模型的预测值。

粒子群算法是模拟群体智能所建立的一种优化算法，是基于鸟类的觅食行为而提出的全局寻优的一种方法。粒子群算法具有收敛速度快、可调参数少等优点，广泛应用于神经网络的优化。为了提高模型的预测精度，可用粒子群算法优化极限学习机的初始输入权值和阈值，其优化流程如图2所示。

图2 pso-elm算法流程图

采用决定系数()和均方根误差(rmse)两个指标来评价模型的精度。表示了预测值与实测值的拟合程度，其值越接近1，表明该模型的拟合效果越高；rmse反映了预测值与实测值的偏离程度，其值越小，表明该模型的预测精度 越高。

经皮尔逊相关性分析，选取的20个植被指数与冬小麦叶片含水量之间的相关程度不同，其中gndvi、evi、nrbdi与叶片水分含量的相关性达到显著水平 (