不同播期冬小麦叶面积指数高光谱遥感监测模型

范剑 尤慧 刘凯文 高华东

（1 荆州市气象局，荆州 434020；2 江汉平原生态气象遥感监测技术协同创新中心，荆州 434025；3 湖北省荆州农业气象试验站，荆州 434025）

0 引言

播期是田间管理的一个重要因子，作物不同播期的变化，其光、温、降水等时空位与质量就会有差别，对作物生长发育、产量、品质等会造成一定的影响[1]。聂居超等[2]通过对6个玉米品种在5个不同播期处理下产量、产量构成因素及相关性状和抗病性研究中指出，适时早播能延长玉米生育期，并获得高产。张宁等[3]在对不同播期夏玉米研究中指出，播期每提前1d，灌浆时间延长0.74 d，产量增加9%～13.9%。刘艳阳[4]指出随播期推迟，小麦全生育期所需天数缩短，播期与产量呈极显著二次曲线关系，并对出粉率、籽粒蛋白质含量、湿面筋含量、淀粉含量等有明显影响。

叶面积指数（Leaf area index，LAI ）是反映作物生理生化过程和生产力状况的良好指标，对农作物长势监测和产量估测等应用具有重要意义[5]。传统的作物LAI监测方法（如方格法、激光叶面积法、数码图像处理法等），虽然能相对准确地获得小范围LAI，但成本高、效率低并具有一定的破坏性[6]。高光谱遥感的出现，为实现大范围LAI的快速、高效、无损监测提供了有力的技术手段，并得到了有效的应用。梁亮等[7]选取18种高光谱指数进行比较分析，筛选出敏感高光谱指数模型OSAVI，并结合地面光谱样本数据建立统一的小麦LAI的反演模型。马驿等[6]通过选取12个光谱特征参数和11个植被指数，建立油菜6叶期至角果期LAI的5种线性和非线性定量反演模型，并指出采用统一参数建模难以准确估算全生育期LAI，不同生长时期需进行植被指数分段建模估算LAI。贺佳等[8]通过相关分析、回归分析等统计方法，利用5种高光谱指数OSAVI、 EVI、NVI、mNDVI、mSRI研究不同氮磷水平下冬小麦不同生育时期叶面积指数监测模型。当前国内外相关研究主要集中在针对作物生长发育及品质产量等的影响研究或是不同氮素水平下的作物LAI反演研究，并建立作物统一的LAI监测模型，然而单一的监测模型易受播种时间、小麦品种、气候等因素的影响，导致LAI反演精度不高。

由于各种环境、气候、社会因素（如地理位置、小麦品种、气候条件等）都会造成小麦播种时间的差异，因此通过分析不同播期冬小麦的冠层光谱响应规律以构建不同播期LAI监测模型，具有一定的意义。本文拟通过对不同播期的田间试验，在同步获取分蘖至乳熟期冠层光谱与LAI的基础上，结合不同播期及混合播期的LAI与植被指数之间的定量关系，构建不同播期LAI高光谱遥感监测模型。此外，将4个播期所有生育期的冠层光谱及LAI数据进行组合，模拟出一个混合播期，并进行研究分析。旨在为不同播期水平下冬小麦长势监测和小麦生产的宏观管理调控等提供理论依据和技术支撑。

1 材料与方法

1.1 试验设计

冬小麦分期播种试验于2016年10月—2017年6月，在中国气象局荆州农业气象试验站（30°21′N，112°09′E）进行，播种期为以当地常年大田实际播种期（10月30日）为界，提前10 d播种为播期1，正常播种为播期2，比正常播种晚10 d播种为播期3，晚20 d为播期4，共播种4期，每期分3个重复，每个试验地块面积为20 m2，各地块间留0.5 m保护间隔。小麦种植品种为“郑麦9023”，种植方式采用南北方向条播，保持良好的通风透光，行距统一为20cm，研究区各个试验地块在小麦生长期采用基本一致的施肥和灌溉管理。

1.2 小麦冠层光谱与参数获取

冬小麦冠层光谱反射率的测定采用ASD HandHeld2手持式地物光谱仪，该光谱仪的光谱范围为325～1075 nm，光谱采样间隔约1.5 nm。选择天气晴朗、无云无风的天气进行监测，测定时间控制在北京时间10:00—14:00，分别测定不同播期小麦分蘖期、拔节期、孕穗期、抽穗期、开花期、乳熟期的冠层光谱反射率，每个生长期每隔8 d左右进行一次观测。测量时，传感器距离小麦冠层顶部高度约1m，探头垂直向下。每个试验样本区每次记录5条光谱曲线，重复测定3次，计算该15条曲线的平均值作为该样本小区的冠层反射率值。每半小时用参考板对仪器进行一次校正以消除环境变化所带来的影响。

1.3 叶面积指数测定

在冬小麦冠层光谱反射率监测的相应位置处，连续取10株长势较均匀的小麦，采用面积系数法计算冬小麦叶面积指数。通过测定单株上各叶片长宽乘积之和与校正系数之积，根据单株叶面积与1 m2密度之积/10000，求得该区小麦叶面积指数。

1.4 植被指数选取

研究采用如表1所示具有代表性的8种植被指数进行不同播期冬小麦LAI的监测，并以不同植被指数为自变量，小麦实测LAI为因变量，进行冬小麦LAI高光谱遥感监测模型的建立。其中播期1、播期2、播期3及播期4参与计算的采样数据量（n）分别为66，64，58及64组。

表1 本文采用的植被指数Table 1 The vegetation indices

1.5 数据处理与分析

利用ViewSpec软件对采集的冠层光谱数据进行预处理，首先对每个采样点采集的光谱数据进行筛选，剔除变化差异较大的误差数据后计算平均，作为该样本点的光谱反射率数据。其次，为了减弱噪声对实际光谱数据的影响，将原始光谱反射率在325～400 nm范围内的数据进行剔除，以便于进行数学分析及监测模型的建立。再基于SPSS和Excel软件分别对4个播期和基于该4个播期的同生育期数据分别进行平均处理后模拟出的混合播期数据进行相关统计计算分析。

2 结果与分析

2.1 不同播期冬小麦光谱变化分析

随着播期的不同，冬小麦冠层反射率光谱存在不同的响应。图1为不同播期小麦在同一时段观测的冠层光谱均值曲线，此时播期1处于孕穗末期，播期2为孕穗始期，播期3和播期4均为拔节期。可以看出在550 nm附近出现小的反射峰，不同播期差异较小；在760～1075 nm，不同播期光谱反射率出现明显差异，反射率随着播期的延后而减小；在970 nm附近光谱反射率出现吸收谷，并随着播期的延后，吸收谷有加深的趋势。

图1 不同播期下冠层光谱反射率的变化Fig. 1 Changes in spectral reflectance of canopy at different sowing dates

2.2 冬小麦冠层反射率与叶面积指数的相关性

将4个播期及利用所有播期模拟的混合播期的冬小麦冠层反射率分别与其对应LAI进行相关分析（图2）。从图中可以看出，4个播期处理及混合播期的LAI随冠层光谱反射率的变化规律相似。从播期1（图2a）中可以看出，在400～727 nm冬小麦冠层光谱反射率与LAI呈显著负相关（P＜0.05），在743～1075 nm呈极显著正相关（P＜0.01）；播期2（图2b）在400～721nm为显著负相关（P＜0.05），在749～1075 nm为极显著正相关（P＜0.01）；播期3（图2c）在400～728 nm为极显著负相关（P＜0.01），在745～1075 nm呈极显著正相关（P＜0.01）；播期4（图2d）在400～717nm为极显著负相关（P＜0.01），在732～1075 nm呈极显著正相关（P＜0.01）；混合播期（图2e）在400～725 nm冬小麦冠层光谱反射率与LAI呈极显著负相关（P＜0.01），在735～1075 nm呈极显著正相关（P＜0.01）。

通过分析，在467 nm蓝光波段处有最大负相关，相关系数约为—0.81；红光波段最大相关在651 nm，相关系数为—0.8，绿光波段最大相关在491 nm，相关系数为—0.801；近红外波段824 nm处有最大正相关，相关系数为0.762。因而植被指数计算时红光波段波长选用651 nm，绿光波段波长为491 nm，蓝光波段波长为467 nm，近红外波段波长为824 nm。

2.3 不同播期及混合播期植被指数相关性分析

图2 不同播期下冬小麦冠层反射率与LAI的相关性（a～e分别为不同播期及混合播期冬小麦冠层光谱反射率与LAI的相关性）Fig.2 Correlation between canopy reflectance and LAI of winter wheat in different sowing dates(a-e shows the correlation between spectral reflectance and LAI of winter wheat canopy at different sowing dates and all sowing dates together)

通过表1计算8种植被指数，分析不同播期及混合播期冬小麦LAI与对应植被指数的相关性（表2）。从表2中可以看出，混合播期不同植被指数与LAI均达到极显著性相关水平（P＜0.01），相关系数（r）范围为0.798～0.891，按照统计学规定，植被指数与LAI 相关系数均达到中度到高度相关水平；播期1的植被指数与LAI为显著相关（P＜0.05），r为0.506～0.799，达中度相关水平；播期2的8种植被指数与LAI为极显著相关水平（P＜0.01），r为0.615～0.898，达中度相关水平；播期3的植被指数与LAI达到极显著相关水平（P＜0.01），r为0.828～0.952，植被指数与LAI为高度相关水平；播期4中植被指数与LAI为极显著相关（P＜0.01），r为0.863～0.894，植被指数与LAI为高度相关；因此所选的8种植被指数均可用来监测不同播期的冬小麦LAI，以作为构建冬小麦LAI监测模型参数。

2.4 不同播期叶面积指数监测模型的拟合与验证

将LAI与8种植被指数分别进行拟合，建立不同播期及混合播期冬小麦LAI分期监测模型（表3），并以试验田实测LAI与反演LAI数据（n=60）进行拟合模型的验证。可以发现，在混合播期中植被指数NDGI与LAI之间线性关系极显著，R2为0.798，F统计量为243.192，实测LAI与反演LAI的均方根误差（RMSE）为0.91；播期1基于EVI2建立的监测模型具有较高的拟合精度，R2和F统计量分别为0.803和76.187，通过对EVI2监测模型的验证，RMSE为0.81，监测效果较好；播期2、播期3均能基于NDGI建立较高精度的拟合模型，R2分别为0.823和0.907，F统计量分别为86.278和127.057，RMSE分别为0.78和0.63，通过验证，NDGI监测模型对播期2及播期3具有较好的验证效果；播期4基于mNDVI建立的监测模型具有较高的拟合精度，R2和F统计量分别为0.819和75.698，RMSE为0.82。

表2 不同播期植被指数与叶面积指数的相关系数Table 2 Correlation coefficients between vegetation indices and leaf area index at different sowing dates

从表3中可知，混合播期及播期2和播期3冬小麦叶面积指数反演模型参数中均以NDGI为最佳拟合光谱参数，播期1和播期4分别以植被指数EVI2和mNDVI拟合效果最好。通过综合比较分析来看，混合播期的拟合模型虽然具有一定的普适性，反演效果较好，但不同播期的分期监测模型较混合播期统一拟合模型效果更佳，验证精度也更高，并且播期3的拟合效果最佳。因此对于不同播期的冬小麦LAI的监测时，应建立不同播期的分期监测模型，其中播期1的监测模型为LAI＝10.401×EVI2－0.456；播期2监测模型为LAI＝63.276×NDGI＋0.5884；播期3为LAI＝66.174×NDGI＋0.5712；播期4为LAI＝14.181×mNDVI－4.8695。

3 结论

1）随着播期的不同，冬小麦冠层反射率光谱存在不同的响应。冬小麦LAI与冠层光谱反射率在近红外波段、红光波段、绿光波段和蓝光波段的最大相关波长分别为824，651，491和467 nm。通过分析不同播期及混合播期8种植被指数与LAI的相关性，发现不同类型的植被指数与LAI呈极显著相关。在混合播期及播期2、播期3中植被指数NDGI与LAI之间相关性最高，播期1和播期4分别以EVI2和mNDVI为最佳。

2）通过监测模型拟合及验证可知，混合播期的拟合模型虽然具有一定的普适性，但是播期1至播期4的分期监测模型较混合播期统一监测模型的拟合效果更好，精度更高，并且播期3的拟合效果及精度最高。因此，反演不同播期的LAI时通过建立统一监测模型的精度较差，应根据不同播期建立LAI分期监测模型更合理，并且LAI估测的最佳播期为播期3。

高光谱数据波段数目多，数据量庞大，且不同植被指数具有各自的特点和优势，但在如何进行植被指数的筛选及监测模型的建立上，仍需不断探索和研究。由于本研究是基于田间小区试验条件下得到的，气候及相关影响因素同质性大，异质性低，所建立的不同播期分期监测模型虽然精确性较高，但适用范围可能较窄，存在一定的局限性。因此，未来的研究方向应集中在不同播期的不同耕作类型、品种、生态条件下的LAI大面积、高精度光谱监测。

表3 不同播期冬小麦LAI与植被指数的定量关系Table 3 Quantitative relationship between LAI and the vegetation indices of winter wheat with different sowing dates