• 《中国科学引文数据库(CSCD)》来源期刊
  • 中国科技期刊引证报告(核心版)期刊
  • 《中文核心期刊要目总览》核心期刊
  • RCCSE中国核心学术期刊

融合无人机多光谱和纹理特征的马铃薯LAI估算

李健, 江洪, 罗文彬, 麻霞, 张雍

李健, 江洪, 罗文彬, 等. 融合无人机多光谱和纹理特征的马铃薯LAI估算[J]. 华南农业大学学报, 2023, 44(1): 93-101. DOI: 10.7671/j.issn.1001-411X.202201002
引用本文: 李健, 江洪, 罗文彬, 等. 融合无人机多光谱和纹理特征的马铃薯LAI估算[J]. 华南农业大学学报, 2023, 44(1): 93-101. DOI: 10.7671/j.issn.1001-411X.202201002
LI Jian, JIANG Hong, LUO Wenbin, et al. Potato LAI estimation by fusing UAV multi-spectral and texture features[J]. Journal of South China Agricultural University, 2023, 44(1): 93-101. DOI: 10.7671/j.issn.1001-411X.202201002
Citation: LI Jian, JIANG Hong, LUO Wenbin, et al. Potato LAI estimation by fusing UAV multi-spectral and texture features[J]. Journal of South China Agricultural University, 2023, 44(1): 93-101. DOI: 10.7671/j.issn.1001-411X.202201002

融合无人机多光谱和纹理特征的马铃薯LAI估算

基金项目: 福建省科技计划引导性项目(2021Y0005);国家重点研发计划(2017YFB0504203)
详细信息
    作者简介:

    李健,硕士研究生,主要从事无人机遥感技术与应用研究,E-mail: 1432492120@qq.com

    通讯作者:

    江 洪, 副研究员,博士,主要从事遥感信息处理与应用研究,E-mail: jh9l0@fzu.edu.cn

  • 中图分类号: S532;S127

Potato LAI estimation by fusing UAV multi-spectral and texture features

  • 摘要:
    目的 

    研究融合无人机遥感影像多光谱信息和纹理特征估算马铃薯Solanum tuberosum叶面积指数(Leaf area index,LAI)方法,提高马铃薯LAI反演精度。

    方法 

    利用大疆P4M无人机采集2021年2-4月南方冬种马铃薯幼苗期、现蕾期、块茎膨大期多光谱影像,用LAI-2000冠层分析仪实测LAI数据。提取影像光谱、纹理等信息,分析植被指数、纹理特征与LAI的相关性,基于R2adj的全子集分析优选特征变量。采用主成分分析,融合光谱和纹理特征,用PCA-MLR(Principal component analysis-multiple linear regression)模型估算马铃薯LAI。

    结果 

    从幼苗期到块茎膨大期,PCA-MLR估算模型优于T-MLR(Texture multiple linear regression)和VI-MLR(Vegetation index multiple linear regression)模型,R2分别为0.73、0.59和0.66。

    结论 

    本研究提出一种估算马铃薯LAI的PCA-MLR方法,为马铃薯的长势监测和田间管理提供数据支持。

    Abstract:
    Objective 

    Develop a method to improve the potato (Solanum tuberosum) leaf area index (LAI) estimation accuracy using the UAV multiple spectral wavebands and texture information.

    Method 

    The DJI P4M drone was used to collect multispectral images of the southern winter potato at seedling period, budding period and tuber swelling period from February to April 2021. LAI data were measured by LAI-2000 canopy analyzer. The spectral and texture characteristics of images were extracted. The correlations between vegetation index, texture characteristics and LAI were analyzed. The selected characteristic variables were analyzed based on subset of adjusted R2adj. The principal component analysis was used to fuse spectrum and texture features, and the principal component analysis-multiple linear regression (PCA-MLR) model was used to estimate potato LAI.

    Result 

    From the seedling period to the tuber swelling period, the PCA-MLR estimation model was better than texture multiple linear regression (T-MLR) and vegetation index multiple linear regression (VI-MLR) model, with R2 of 0.73, 0.59 and 0.66 respectively.

    Conclusion 

    This study proposed a method of PCA-MLR to estimate the potato LAI and improve the levels of the potato growth monitoring and field management.

  • 2015年以来,我国启动了马铃薯主粮化战略,马铃薯成为我国又一主粮作物[1]。叶面积指数(Leaf area index,LAI)作为作物生长的重要参数,与地上生物量积累和地下块茎产量有密切的关系,是监测作物长势、指导田间管理的重要指标[2]。当前,地面实测与遥感估算是测量农作物LAI的主要方式。地面实地测量分为直接测量法与间接测量法[3]。直接测量法是利用直接测定的植被叶面积计算LAI,间接测量法是用光学仪器测量有关参数得到LAI[4]。卫星遥感估算法主要采用光学遥感影像估测农作物LAI,估算面积大,但受天气、时空分辨率的影响。相比卫星遥感,无人机遥感具有空间分辨率高、能准确反映地块作物长势、监测灵活机动等优势,适用于田块尺度作物的信息获取与LAI估算。国内外学者利用光学遥感影像计算植被指数(Vegetation index,VI)来估算LAI已有大量的研究[5],研究表明VI能较好地估算LAI[6]。其中,归一化植被指数(Normalized difference vegetation index,NDVI)被广泛用于建立LAI估算模型[7],但茂密作物的NDVI会出现饱和现象[8],在估算LAI时存在一定的局限性。因此,综合利用遥感影像光谱、纹理特征估算作物生长参数的研究逐渐增多。VI与纹理特征融合可以提高生物量[9]、叶绿素含量[10]以及氮含量[11]估算精度,将两者结合估算作物氮含量、生物量的研究较多,如玉米、小麦、水稻等[12-14],对马铃薯LAI估算的研究较少。前人研究大多从植被光谱指数或者单一的纹理特征入手,忽略了无人机多光谱影像光谱、纹理信息的融合,本研究基于高分辨率的大疆P4M无人机多光谱影像,比较马铃薯不同生育期光谱、纹理单一指标以及光谱和纹理融合方法,进行马铃薯LAI估算。

    研究区位于福建省宁德市霞浦县沙江镇(E119°59′,N26°47′)粮食产能区“五星”示范片区,面积约37.33 hm2。该地属于亚热带季风气候,年平均气温为17°~23°,春季雨水丰沛,夏有台风,冬季少见霜雪。研究区是福建省农科院南方马铃薯研究基地之一,田间管理完善,进行2个施肥试验,分别为肥料联合筛选试验和氮肥分期施用试验。

    2种试验各6个处理,3次随机重复试验。肥料联合筛选试验处理有:不施肥、常用有机肥+常规化肥、沃尔田生物有机肥、沃尔田生物有机肥+缓释高钾肥、生命源黄腐酸生物有机肥、生命源黄腐酸生物有机肥+缓释高钾肥;氮肥分期施用试验中,氮、磷、钾用量分别为225、150和300 kg/hm2,分为基肥(播种前施用)和追肥(播种后施用),试验处理及各种肥料质量分数为:不施肥、100%基肥+无追肥、75%基肥+25%追肥、50%基肥+50%追肥、25%基肥+75%追肥、无基肥+100%追肥。

    于马铃薯幼苗期(2021年2月20日)、现蕾期(2021年3月15日)、块茎膨大期(2021年4月1日) 3个关键生育期分别采集马铃薯LAI数据,测量样本均来自2种试验处理地块,每次随机重复试验中至少采集1个样本,测量样地在田间有相应的标记牌,每个生育期45个样本,3个生育期LAI值的箱线图如图1a所示。LAI的测量使用LAI-2000冠层分析仪,测量时间在17:00—18:30,排除太阳直射对仪器测量的影响。每次测量测1次A值(天空光)、4次B值(冠层下方)。同时,选择样区内多个相同植株,标记多个样点,在正午日照充足且无云层遮盖的时候,用ASD 光谱仪测量不同生育期的光谱(图1b)。

    图  1  实测LAI值的箱线图(a)和实测光谱曲线(b)
    Figure  1.  Box plot of the measured LAI value (a) and measured spectral curve (b)

    利用大疆P4M多光谱无人机在无风且阳光充足的中午时段采集影像,飞行航高100 m,速度50 km/h,航向和侧向重叠率在80%以上,像片分辨率为1600像素×1300像素,实际地面分辨率约5 cm,相机具体参数见表1。3次采集的数据运用大疆智图软件拼接,利用反射率为0.5的标准定标版进行辐射校正,得到反射率数据,用ENVI软件裁剪拼接的影像。

    表  1  大疆P4M多光谱相机波段信息
    Table  1.  Band information of DJI P4M multispectral camera
    波段 Band 中心波长/nm Center wavelength 波长宽度/nm Wavelength width
    蓝波段 Blue band 450 32
    绿波段 Green band 560 32
    红波段 Red band 650 32
    红边波段 Red edge band 730 32
    近红外波段 Near-infrared band 840 52
    下载: 导出CSV 
    | 显示表格

    通过预处理的多光谱影像计算植被指数和纹理特征,提取马铃薯植被。将实测的LAI与植被指数、纹理特征进行相关性分析,基于R2adj全子集回归分析筛选最优特征变量,用主成分分析(Principal component analysis,PCA)融合光谱、纹理得到主要成分进行马铃薯LAI估算。流程如图2所示。

    图  2  LAI估算流程图
    Figure  2.  LAI estimation flowchart

    1)马铃薯提取。用支持向量机(Support vector machine,SVM)分类方法提取马铃薯植被像元,去除土壤、阴影、地膜等背景。2)植被指数计算。根据大疆P4M多光谱波段、实测的不同生育期光谱曲线和已有的研究,选取与叶面积指数相关的植被指数(表2)。3)纹理特征计算。计算8类基于灰度共生矩阵(Gray-level co-occurrence matrix,GLCM)的纹理特征值有:均值(Mean)、协同性(Homogeneity, Hom)、相异性(Dissimilarity, Dis)、信息熵(Entropy, En)、二阶矩(Second moment, Sm)、相关性(Correlation, Cor)、对比度(Contrast, Contr)方差(Variance, Var),共40个纹理特征值。4)特征优选。将LAI与植被指数、纹理特征进行相关性分析,挑选出相关性较高的特征,并引入调整的R2(R2adj)对植被指数、纹理特征进行全子集回归分析,对相关性较高的特征进行二次优选。5)光谱、纹理特征融合。利用PCA融合二次优选特征,用少量的主要成分表示原始光谱、纹理信息,融合成新的指标。6)模型比较及验证。将PCA方法融合的主要成分进行多元回归建模,并与植被指数、纹理特征回归模型进行比较。在所测试样本中,分别由30个建模样本和15个验证样本组成。利用R2R2adj和均方根误差(RMSE)作为模型的评价指标。7)选较优模型进行马铃薯LAI估算。

    表  2  多光谱植被指数
    Table  2.  Multispectral vegetation index
    植被指数 Vegetation index 计算公式1) Calculating formula 文献来源 Literature source
    NDVI NDVI= (NIR−Red)/ (NIR−Red) [15]
    GNDVI GNDVI=(NIR−Green)/ (NIR−Green) [16]
    NDRE NDRE=(NIR−Rededge)/ (NIR+Rededge) [17]
    LCI LCI=(NIR−Rededge)/ (NIR+Red) [18]
    OSAVI OSAVI=(NIR−Red)/ (NIR+Red+0.16) [19]
    DVI DVI=NIR−Red [20]
    DVI_GRE DVI_GRE=NIR−Green
    DVI_EDG DVI_EDG=NIR−Rededge
    RVI RVI=NIR/Red [21]
    RVI_GRE RVI_GRE=NIR/Green
    RVI_EDG RVI_EDG=NIR/Rededge
    RDVI RDVI=(NDVI)1/2 [22]
    RDVI_GRE RDVI=(GNDVI)1/2
    RDVI_EDG RDVI_EDG=(NDRE)1/2
    G G=Green
    R R=Red
    EDG EDG=Rededge
    NIR NIR=NIR
     1) Green、Red、Rededge和NIR分别表示绿、红、红边及近红外波段反射率  1) Green, Red, Rededge and NIR represent the reflectance of green, red, red edge and near-infrared bands respectively
    下载: 导出CSV 
    | 显示表格

    马铃薯叶片、地膜、土壤、阴影分类结果如图3所示,每个生育期总体精度和用户精度达到95%以上。将马铃薯纯植被像元区域转换为矢量,裁剪出各地块马铃薯植被区域。

    图  3  马铃薯提取前后的马铃薯植被区域对比
    Figure  3.  Comparison of potato planting area before and after potato extraction

    幼苗期、现蕾期、块茎膨大期马铃薯多光谱植被指数与LAI的相关性见表3。幼苗期除Green、Rededge、NIR外,其余植被指数与LAI相关系数绝对值均为0.50以上,P<0.01;现蕾期,除Green、Rededge外,植被指数与LAI相关系数绝对值大多为0.60~0.80,P<0.01;块茎膨大期,NDVI、GNDVI、RVI、RDVI、OSAVI、Red和RDVI_GRE相关系数绝对值约为0.50~0.65,P<0.01。

    表  3  光谱特征与LAI的相关系数绝对值
    Table  3.  Absolute value of correlation coefficient between spectral feature and LAI
    植被指数 Vegetation index 幼苗期 Seedling period 现蕾期 Budding period 块茎膨大期 Tuber swelling period
    Green 0.49 0.46 0.30
    Red 0.73 0.58 0.64
    Rededge 0.30 0.48 0.15
    NIR 0.43 0.61 0.24
    NDVI 0.77 0.76 0.73
    GNDVI 0.66 0.70 0.65
    NDRE 0.53 0.63 0.30
    OSAVI 0.69 0.70 0.66
    LCI 0.58 0.66 0.33
    RVI 0.76 0.79 0.73
    RVI_GRE 0.54 0.64 0.54
    RVI_EDG 0.67 0.74 0.24
    DVI 0.54 0.64 0.35
    DVI_GRE 0.50 0.65 0.32
    GVI_REDED 0.51 0.70 0.30
    RDVI 0.77 0.76 0.73
    RDVI_GRE 0.66 0.69 0.54
    RDVI_EDG 0.52 0.62 0.30
    下载: 导出CSV 
    | 显示表格

    马铃薯各生育期纹理特征与LAI的相关性如表4所示。幼苗期,b1_Contr、b1_Hom、b2_Cor、b2_Mean、b3_Dis、b3_Ent、b3_Mean、b4_Cor和b5_Mean的相关系数绝对值为0.40~0.60,P<0.01。现蕾期,b1_Contr、b1_Hom、b2_Cor、b3_Mean、b3_Dis、b4_Cor和b5_Mean相关系数绝对值为0.40~0.60,P<0.01。块茎膨大期,b1_Mean、b3_Cor、b3_Ent、b3_Hom、b3_Mean、b4_Hom、b5_Hom和b5_Mean相关系数绝对值为0.40~0.60,P<0.01。

    表  4  纹理特征与LAI的相关系数绝对值(|r|)
    Table  4.  Absolute value of correlation coefficient (|r|) between texture feature and LAI
    生育期 Growth period 纹理特征1)Texture feature |r|
    幼苗期 Seedling period b1_Contr 0.57
    b1_Hom 0.63
    b2_Cor 0.58
    b2_Mean 0.42
    b3_Dis 0.55
    b3_En 0.44
    b3_Mean 0.30
    b3_Var 0.38
    b4_Cor 0.56
    b5_Mean 0.56
    现蕾期 Budding period b1_Contr 0.53
    b1_En 0.50
    b1_Mean 0.36
    b2_Cor 0.46
    b2_Con 0.57
    b2_Dis 0.52
    b2_En 0.65
    b3_Mean 0.49
    b4_Cor 0.62
    b5_Mean 0.54
    块茎膨大期 Tuber swelling period b1_Mean 0.63
    b2_Var 0.24
    b3_Com 0.53
    b3_En 0.47
    b3_Hom 0.43
    b3_Sm 0.62
    b3_Mean 0.69
    b4_Hom 0.63
    b5_Mean 0.45
    b5_Hom 0.53
     1)b1、b2、b3、b4和b5分别代表蓝、绿、红、红边、近红外波段  1) b1, b2, b3, b4 and b5 represent blue, green, red, red edge and near-infrared bands respectively
    下载: 导出CSV 
    | 显示表格

    依据VI、纹理特征与LAI相关性分析结果,基于R2adj全子集回归构建LAI多元线性估算模型,计算模型R2adjR2、RMSE。由图4光谱特征可知,全子集优化植被指数效果较好的指标,在幼苗期有RDVI、RVI_GRE、Red、RDVI_EDG、GVI_EDG、LCI和GNDVI,在现蕾期有RDVI、NDRE、LCI、GNDVI和Green,在块茎膨大期有NDVI、OSAVI、RVI_EDG、LCI、NIR 和 GNDVI。优选的纹理特征重要度效果较好指标,在幼苗期有b2_ Hom、b2_Dis、b1_Mean和b3_Contr,在现蕾期有b3_Mean、b4_Cor和b5_Mean,在块茎膨大期有b3_Mean、b5_Mean、b5_Hom、b3_Hom、b3_En和b1_Mean(图5)。

    图  4  基于R2adj全子集回归优选的光谱特征重要度
    Figure  4.  Importance of spectral features based on R2adj full subset regression optimization
    图  5  基于R2adj全子集回归优选的纹理特征重要度
    b1、b2、b3、b4和b5分别代表蓝、绿、红、红边、近红外波段
    Figure  5.  Importance of texture features based on R2adj full subset regression optimization
    b1, b2, b3, b4 and b5 represent blue, green, red, red edge and near-infrared bands respectively

    利用PCA方法对植被指数和纹理特征进行降维,将累计贡献率大于99%的成分用来反映原始信息。确定3个生育期最佳主成分为6、3和4个,基于PCA主成分多元线性回归(Principal component analysis-multiple linear regression,PCA-MLR)建模,并与纹理特征多元回归(Texture multiple linear regression,T-MLR)、植被指数多元回归(Vegetation index multiple linear regression,VI-MLR)模型进行对比。在马铃薯的3个不同生育期内,PCA-MLR模型精度更高,幼苗期、现蕾期与块茎膨大期的R2分别约为0.73、0.59和0.66,RMSE分别为0.426、0.540和0.432,模型稳定性较好,其R2R2adj都较高,RMSE明显降低(表5)。

    表  5  马铃薯各生育期LAI估算建模比较1)
    Table  5.  Comparison of LAI estimation modeling for potatoes at different growth stages
    生育期 Growth Period 自变量 Independent variable 建模方式 Modeling method R2 R2adj RMSE 优化的回归模型 Optimized regression model
    幼苗期 Seedling period RDVI、 RVI_GRE、 Red、 RDVI_EDG、 GVI_EDG、 LCI GNDVI VI-MLR 0.647 0.587 0.490
    b2_Hom、b2_Dis、 b1_Mean、b3_Contr T-MLR 0.637 0.581 0.511
    PC1、PC2、PC3、PC4、PC5、PC6 PCA-MLR 0.739 0.674 0.426 LAI=PC6×4.16+PC4×0.27+ PC2×0.89−PC1×0.19− PC3×3.06−PC5×3.17+0.9
    现蕾期 Budding period RDVI、NDRE、LCI、GNDVI、Green VI-MLR 0.483 0.470 0.571
    b3_Mean、b4_Cor、b5_Mean T-MLR 0.465 0.417 0.609
    PC1、PC2、PC3 PCA-MLR 0.592 0.558 0.542 LAI=PC3×1.17−PC2×0.84− PC1×0.11+0.34
    块茎膨大期 Tuber swelling period NDVI、OSAVI、RVI_EDG、LCI、NIR、GNDVI VI-MLR 0.608 0.539 0.540
    b3_Mean、b5_Mean、b5_Hom、 b3_Hom、b3_En、b1_Mean T-MLR 0.594 0.561 0.536
    PC1、PC2、PC3、PC4 PCA-MLR 0.659 0.592 0.432 LAI=PC1×0.53+PC2×0.26− PC3×0.72−PC4×0.15−0.08
     1)PC1~PC6表示对应的主成分  1)PC1−PC6 represent the corresponding principal components respectively
    下载: 导出CSV 
    | 显示表格

    PCA-MLR模型较优,利用其进行LAI估算作图(图6),实测值与预测值如图7所示。从幼苗期到块茎膨大期,LAI逐渐增大,图中颜色逐渐加深。幼苗期叶子较少,块茎膨大期马铃薯茎叶增长都达到了峰值,是叶面积最大的时期[23],叶子相对茂密,LAI较高,不同处理的空间分布更为明显。

    图  6  施肥处理地块LAI空间分布
    Figure  6.  Spatial distribution of LAI in fertilized plots
    图  7  马铃薯各生育期LAI预测值与实测值
    Figure  7.  Estimated and measured values of LAI of potato at different growth stages

    2种施肥试验统计如表6所示,不同处理地块的LAI在3个时期皆大于不施肥处理。肥料联合筛选试验地块的3个生育期LAI均值分别为2.89、4.88和5.21,不同组合肥料施用地块LAI普遍高于单种肥料施用地块LAI。常用有机肥+常规化肥处理LAI在块茎膨大期最大,为5.67,比不施肥处理同期LAI增加约20%。单种肥料施用时,沃尔田生物有机肥处理块茎膨大期LAI最大,为5.00,相比不施肥处理同期LAI增加约5%;在氮肥分期施用试验中,3个生育期LAI均值分别为1.40、4.13和4.56。75%基肥+25%追肥处理地块块茎膨大期LAI最大,为4.92,与同期不施肥处理相比,LAI增加约14%。

    表  6  施肥试验各处理地块LAI统计结果
    Table  6.  Statistical results of LAI of each treatment plot in fertilization experiment
    施肥试验类型 Type of ferilization experiment 处理 Treatment 幼苗期 Seedling period 现蕾期 Budding period 块茎膨大期 Tuber swelling period
    肥料联合筛选试验 Fertilizer joint screening experiment 常用有机肥+常规化肥 Commonly used organic fertilizer + conventional chemical fertilizer 2.90 4.95 5.67
    生命源黄腐酸生物有机肥+缓释高钾肥 Life source fulvic acid bio-organic fertilizer + slow-release high potassium fertilizer 2.91 5.08 5.63
    沃尔田生物有机肥+缓释高钾肥 Waltian bio-organic fertilizer + slow-release high potassium fertilizer 2.99 5.05 5.46
    沃尔田生物有机肥 Waltian bio-organic fertilizer 2.82 4.78 5.00
    不施肥 No fertilization 2.81 4.61 4.75
    生命源黄腐酸生物有机肥 Life source fulvic acid bio-organic fertilizer 2.88 4.79 4.72
    均值 Mean 2.89 4.88 5.21
    氮肥分期施用试验 Nitrogen fertilizer application experiment by stages 25%基肥+75%追肥 25% basic fertilization+75% additional fertilization 1.37 4.09 4.68
    50%基肥+50%追肥 50% basic fertilization+50% additional fertilization 1.36 4.12 4.40
    100%基肥+无追肥 100% basic fertilization+no additional fertilization 1.36 4.11 4.44
    不施肥 No fertilization 1.33 4.03 4.25
    无基肥+100%追肥 No basic fertilization+100% additional fertilization 1.55 4.24 4.69
    75%基肥+25%追肥 75% basic fertilization+25% additional fertilization 1.45 4.16 4.92
    均值 Mean 1.4 4.13 4.56
    下载: 导出CSV 
    | 显示表格

    从大疆P4M无人机影像的光谱、纹理信息及二者之间的融合入手,探究其对马铃薯LAI的估算能力。模型中马铃薯植被的植被指数和纹理特征的选取只考虑了不同生育期马铃薯冠层反射差异(即在红、蓝波段吸收,在绿波段反射以及在红边和近红外波段高反射的特性)。植被指数选用大多来自波段的比值、差值组合,而纹理特征选用仅能表现单波段植被冠层结构,没有进行组合,下一步研究中,还可以根据不同处理构建新的指标,更精细地反应施肥、水分等特征。在构造指标时,本文只采用了PCA主成分融合的方法,后面的研究还可考虑层次分析法、加权法等,比较不同方法的精度,加入其他纹理特征提取方法,增加样本量提升模型稳定性,以及探究将3个生育期结合起来估算LAI。也可以进一步探讨将马铃薯LAI估算结果应用于产量估算、作物表型特征研究、肥料筛选及施药指导等。

    基于大疆P4M无人机多光谱数据,通过融合其光谱和纹理信息估算马铃薯LAI:1)在现蕾期、幼苗期和块茎膨大期,基于植被指数的LAI估算模型比纹理特征模型表现得更为优越;2)基于PCA-MLR的马铃薯LAI估算模型优于基于多光谱植被指数和纹理特征的多元回归模型,其决定系数有所提高,RMSE也明显降低;3) 3个不同生育期,马铃薯LAI值不断变化,从幼苗期LAI均值小于2,到现蕾期LAI均值为4~5,块茎膨大期LAI达到峰值,接近6,这与原始测量数据相吻合。综上,通过光谱、纹理信息融合的PCA-MLR模型,能够较好地估算马铃薯LAI,可以为马铃薯长势监测和参数估算提供一种可行的方法。

  • 图  1   实测LAI值的箱线图(a)和实测光谱曲线(b)

    Figure  1.   Box plot of the measured LAI value (a) and measured spectral curve (b)

    图  2   LAI估算流程图

    Figure  2.   LAI estimation flowchart

    图  3   马铃薯提取前后的马铃薯植被区域对比

    Figure  3.   Comparison of potato planting area before and after potato extraction

    图  4   基于R2adj全子集回归优选的光谱特征重要度

    Figure  4.   Importance of spectral features based on R2adj full subset regression optimization

    图  5   基于R2adj全子集回归优选的纹理特征重要度

    b1、b2、b3、b4和b5分别代表蓝、绿、红、红边、近红外波段

    Figure  5.   Importance of texture features based on R2adj full subset regression optimization

    b1, b2, b3, b4 and b5 represent blue, green, red, red edge and near-infrared bands respectively

    图  6   施肥处理地块LAI空间分布

    Figure  6.   Spatial distribution of LAI in fertilized plots

    图  7   马铃薯各生育期LAI预测值与实测值

    Figure  7.   Estimated and measured values of LAI of potato at different growth stages

    表  1   大疆P4M多光谱相机波段信息

    Table  1   Band information of DJI P4M multispectral camera

    波段 Band 中心波长/nm Center wavelength 波长宽度/nm Wavelength width
    蓝波段 Blue band 450 32
    绿波段 Green band 560 32
    红波段 Red band 650 32
    红边波段 Red edge band 730 32
    近红外波段 Near-infrared band 840 52
    下载: 导出CSV

    表  2   多光谱植被指数

    Table  2   Multispectral vegetation index

    植被指数 Vegetation index 计算公式1) Calculating formula 文献来源 Literature source
    NDVI NDVI= (NIR−Red)/ (NIR−Red) [15]
    GNDVI GNDVI=(NIR−Green)/ (NIR−Green) [16]
    NDRE NDRE=(NIR−Rededge)/ (NIR+Rededge) [17]
    LCI LCI=(NIR−Rededge)/ (NIR+Red) [18]
    OSAVI OSAVI=(NIR−Red)/ (NIR+Red+0.16) [19]
    DVI DVI=NIR−Red [20]
    DVI_GRE DVI_GRE=NIR−Green
    DVI_EDG DVI_EDG=NIR−Rededge
    RVI RVI=NIR/Red [21]
    RVI_GRE RVI_GRE=NIR/Green
    RVI_EDG RVI_EDG=NIR/Rededge
    RDVI RDVI=(NDVI)1/2 [22]
    RDVI_GRE RDVI=(GNDVI)1/2
    RDVI_EDG RDVI_EDG=(NDRE)1/2
    G G=Green
    R R=Red
    EDG EDG=Rededge
    NIR NIR=NIR
     1) Green、Red、Rededge和NIR分别表示绿、红、红边及近红外波段反射率  1) Green, Red, Rededge and NIR represent the reflectance of green, red, red edge and near-infrared bands respectively
    下载: 导出CSV

    表  3   光谱特征与LAI的相关系数绝对值

    Table  3   Absolute value of correlation coefficient between spectral feature and LAI

    植被指数 Vegetation index 幼苗期 Seedling period 现蕾期 Budding period 块茎膨大期 Tuber swelling period
    Green 0.49 0.46 0.30
    Red 0.73 0.58 0.64
    Rededge 0.30 0.48 0.15
    NIR 0.43 0.61 0.24
    NDVI 0.77 0.76 0.73
    GNDVI 0.66 0.70 0.65
    NDRE 0.53 0.63 0.30
    OSAVI 0.69 0.70 0.66
    LCI 0.58 0.66 0.33
    RVI 0.76 0.79 0.73
    RVI_GRE 0.54 0.64 0.54
    RVI_EDG 0.67 0.74 0.24
    DVI 0.54 0.64 0.35
    DVI_GRE 0.50 0.65 0.32
    GVI_REDED 0.51 0.70 0.30
    RDVI 0.77 0.76 0.73
    RDVI_GRE 0.66 0.69 0.54
    RDVI_EDG 0.52 0.62 0.30
    下载: 导出CSV

    表  4   纹理特征与LAI的相关系数绝对值(|r|)

    Table  4   Absolute value of correlation coefficient (|r|) between texture feature and LAI

    生育期 Growth period 纹理特征1)Texture feature |r|
    幼苗期 Seedling period b1_Contr 0.57
    b1_Hom 0.63
    b2_Cor 0.58
    b2_Mean 0.42
    b3_Dis 0.55
    b3_En 0.44
    b3_Mean 0.30
    b3_Var 0.38
    b4_Cor 0.56
    b5_Mean 0.56
    现蕾期 Budding period b1_Contr 0.53
    b1_En 0.50
    b1_Mean 0.36
    b2_Cor 0.46
    b2_Con 0.57
    b2_Dis 0.52
    b2_En 0.65
    b3_Mean 0.49
    b4_Cor 0.62
    b5_Mean 0.54
    块茎膨大期 Tuber swelling period b1_Mean 0.63
    b2_Var 0.24
    b3_Com 0.53
    b3_En 0.47
    b3_Hom 0.43
    b3_Sm 0.62
    b3_Mean 0.69
    b4_Hom 0.63
    b5_Mean 0.45
    b5_Hom 0.53
     1)b1、b2、b3、b4和b5分别代表蓝、绿、红、红边、近红外波段  1) b1, b2, b3, b4 and b5 represent blue, green, red, red edge and near-infrared bands respectively
    下载: 导出CSV

    表  5   马铃薯各生育期LAI估算建模比较1)

    Table  5   Comparison of LAI estimation modeling for potatoes at different growth stages

    生育期 Growth Period 自变量 Independent variable 建模方式 Modeling method R2 R2adj RMSE 优化的回归模型 Optimized regression model
    幼苗期 Seedling period RDVI、 RVI_GRE、 Red、 RDVI_EDG、 GVI_EDG、 LCI GNDVI VI-MLR 0.647 0.587 0.490
    b2_Hom、b2_Dis、 b1_Mean、b3_Contr T-MLR 0.637 0.581 0.511
    PC1、PC2、PC3、PC4、PC5、PC6 PCA-MLR 0.739 0.674 0.426 LAI=PC6×4.16+PC4×0.27+ PC2×0.89−PC1×0.19− PC3×3.06−PC5×3.17+0.9
    现蕾期 Budding period RDVI、NDRE、LCI、GNDVI、Green VI-MLR 0.483 0.470 0.571
    b3_Mean、b4_Cor、b5_Mean T-MLR 0.465 0.417 0.609
    PC1、PC2、PC3 PCA-MLR 0.592 0.558 0.542 LAI=PC3×1.17−PC2×0.84− PC1×0.11+0.34
    块茎膨大期 Tuber swelling period NDVI、OSAVI、RVI_EDG、LCI、NIR、GNDVI VI-MLR 0.608 0.539 0.540
    b3_Mean、b5_Mean、b5_Hom、 b3_Hom、b3_En、b1_Mean T-MLR 0.594 0.561 0.536
    PC1、PC2、PC3、PC4 PCA-MLR 0.659 0.592 0.432 LAI=PC1×0.53+PC2×0.26− PC3×0.72−PC4×0.15−0.08
     1)PC1~PC6表示对应的主成分  1)PC1−PC6 represent the corresponding principal components respectively
    下载: 导出CSV

    表  6   施肥试验各处理地块LAI统计结果

    Table  6   Statistical results of LAI of each treatment plot in fertilization experiment

    施肥试验类型 Type of ferilization experiment 处理 Treatment 幼苗期 Seedling period 现蕾期 Budding period 块茎膨大期 Tuber swelling period
    肥料联合筛选试验 Fertilizer joint screening experiment 常用有机肥+常规化肥 Commonly used organic fertilizer + conventional chemical fertilizer 2.90 4.95 5.67
    生命源黄腐酸生物有机肥+缓释高钾肥 Life source fulvic acid bio-organic fertilizer + slow-release high potassium fertilizer 2.91 5.08 5.63
    沃尔田生物有机肥+缓释高钾肥 Waltian bio-organic fertilizer + slow-release high potassium fertilizer 2.99 5.05 5.46
    沃尔田生物有机肥 Waltian bio-organic fertilizer 2.82 4.78 5.00
    不施肥 No fertilization 2.81 4.61 4.75
    生命源黄腐酸生物有机肥 Life source fulvic acid bio-organic fertilizer 2.88 4.79 4.72
    均值 Mean 2.89 4.88 5.21
    氮肥分期施用试验 Nitrogen fertilizer application experiment by stages 25%基肥+75%追肥 25% basic fertilization+75% additional fertilization 1.37 4.09 4.68
    50%基肥+50%追肥 50% basic fertilization+50% additional fertilization 1.36 4.12 4.40
    100%基肥+无追肥 100% basic fertilization+no additional fertilization 1.36 4.11 4.44
    不施肥 No fertilization 1.33 4.03 4.25
    无基肥+100%追肥 No basic fertilization+100% additional fertilization 1.55 4.24 4.69
    75%基肥+25%追肥 75% basic fertilization+25% additional fertilization 1.45 4.16 4.92
    均值 Mean 1.4 4.13 4.56
    下载: 导出CSV
  • [1] 卢肖平. 马铃薯主粮化战略的意义、瓶颈与政策建议[J]. 华中农业大学学报(社会科学版), 2015(3): 1-7. doi: 10.13300/j.cnki.hnwkxb.2015.03.001
    [2]

    LUO S, HE Y, LI Q, et al. Nondestructive estimation of potato yield using relative variables derived from multi-period LAI and hyperspectral data based on weighted growth stage[J]. Plant Methods, 2020, 16(1): 150. doi: 10.1186/s13007-020-00693-3.

    [3] 常好雪, 蔡晓斌, 陈晓玲, 等. 基于实测光谱的植被指数对水稻叶面积指数的响应特征分析[J]. 光谱学与光谱分析, 2018, 38(1): 205-211.
    [4] 孙越, 顾祝军, 李栋梁. 无人机与卫星影像的叶面积指数遥感反演研究[J]. 测绘科学, 2021, 46(2): 106-112. doi: 10.16251/j.cnki.1009-2307.2021.02.016
    [5] 李月, 何宏昌, 王晓飞, 等. 农作物冠层光谱分析及反演技术综述[J]. 测绘通报, 2019(9): 13-17. doi: 10.13474/j.cnki.11-2246.2019.0277
    [6]

    HUNT E R, HIVELY W D, FUJIKAWA S J, et al. Acquisition of NIR-green-blue digital photographs from unmanned aircraft for crop monitoring[J]. Remote Sensing, 2010, 2(1): 290-305. doi: 10.3390/rs2010290

    [7]

    KAMAL M, SIDIK F, PRANANDA A R A, et al. Mapping leaf area index of restored mangroves using WorldView-2 imagery in Perancak Estuary, Bali, Indonesia[J]. Remote Sensing Applications: Society and Environment, 2021, 23: 100567. doi: 10.1016/j.rsase.2021.100567.

    [8] 宋开山, 张柏, 王宗明, 等. 基于人工神经网络的大豆叶面积高光谱反演研究[J]. 中国农业科学, 2006, 39(6): 1138-1145. doi: 10.3321/j.issn:0578-1752.2006.06.007
    [9] 刘畅, 杨贵军, 李振海, 等. 融合无人机光谱信息与纹理信息的冬小麦生物量估测[J]. 中国农业科学, 2018, 51(16): 3060-3073. doi: 10.3864/j.issn.0578-1752.2018.16.003
    [10] 陈鹏, 冯海宽, 李长春, 等. 无人机影像光谱和纹理融合信息估算马铃薯叶片叶绿素含量[J]. 农业工程学报, 2019, 35(11): 63-74. doi: 10.11975/j.issn.1002-6819.2019.11.008
    [11] 杨福芹, 冯海宽, 肖天豪, 等. 融合无人机影像光谱与纹理特征的冬小麦氮营养指数估算[J]. 农业现代化研究, 2020, 41(4): 718-726. doi: 10.13872/j.1000-0275.2020.0061
    [12]

    RYU C, SUGURI M, UMEDA M. Multivariate analysis of nitrogen content for rice at the heading stage using reflectance of airborne hyperspectral remote sensing[J]. Field Crops Research, 2011, 122(3): 214-224. doi: 10.1016/j.fcr.2011.03.013

    [13] 刘昌华, 王哲, 陈志超等. 基于无人机遥感影像的冬小麦氮素监测[J]. 农业机械学报, 2018, 49(6): 207-214. doi: 10.6041/j.issn.1000-1298.2018.06.024
    [14]

    YUE J, YANG G, TIAN Q, et al. Estimate of winter-wheat above-ground biomass based on UAV ultrahigh-ground-resolution image textures and vegetation indices[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2019, 150: 226-244. doi: 10.1016/j.isprsjprs.2019.02.022

    [15]

    TANAKA S, KAWAMURA K, MAKI M, et al. Spectral index for quantifying leaf area index of winter wheat by field hyperspectral measurements: A case study in Gifu prefecture, central Japan[J]. Remote Sensing, 2015, 7(5): 5329-5346. doi: 10.3390/rs70505329

    [16] 王来刚, 徐建华, 贺佳, 等. 基于无人机遥感的玉米叶面积指数与产量估算[J]. 玉米科学, 2020, 28(6): 88-93. doi: 10.13597/j.cnki.maize.science.20200613
    [17]

    ZHANG J, WANG C, YANG C, et al. Assessing the effect of real spatial resolution of in situ UAV multispectral images on seedling rapeseed growth monitoring[J]. Remote Sensing, 2020, 12(7): 1207. doi: 10.3390/rs12071207.

    [18] 李宗南. 冬小麦长势遥感监测指标研究[D]. 北京: 中国农业科学院, 2010.
    [19]

    LIANG L, DI L, ZHANG L, et al. Estimation of crop LAI using hyperspectral vegetation indices and a hybrid inversion method[J]. Remote Sensing of Environment, 2015, 165: 123-134. doi: 10.1016/j.rse.2015.04.032

    [20] 夏天, 吴文斌, 周清波, 等. 冬小麦叶面积指数高光谱遥感反演方法对比[J]. 农业工程学报, 2013, 29(3): 139-147.
    [21] 姚雄, 余坤勇, 刘健. 基于无人机多光谱遥感的马尾松林叶面积指数估测[J]. 农业机械学报, 2021, 52(7): 213-221. doi: 10.6041/j.issn.1000-1298.2021.07.022
    [22]

    ROUJEAN J L, BREON F M. Estimating PAR absorbed by vegetation from bidirectional reflectance measurements[J]. Remote Sensing of Environment, 1995, 51(3): 375-384. doi: 10.1016/0034-4257(94)00114-3

    [23] 周岑岑. 马铃薯生育期及形态建成的模拟研究[D]. 武汉: 华中农业大学, 2015.
  • 期刊类型引用(4)

    1. 齐浩,孙海芳,吕亮杰,李偲,闵家楠,侯亮. 基于无人机多光谱信息与纹理特征融合的小麦叶面积指数估测. 农业机械学报. 2025(03): 334-344 . 百度学术
    2. 汪煜,谈俊忠,张梅,朱振宇. 基于无人机遥感的互花米草长势参数提取与相关性分析. 海洋环境科学. 2024(01): 100-108+118 . 百度学术
    3. 刘昕哲,武璐,陈李金,马宇帆,李涛,吴婷婷. 基于半经验半机理建模的冬小麦LAI反演及长势评估. 农业工程学报. 2024(01): 162-170 . 百度学术
    4. 吴路华,杨东妮,陈丹. 贵州武陵山区植被LAI时空演变及其未来持续性特征. 科学技术创新. 2023(22): 31-36 . 百度学术

    其他类型引用(8)

图(7)  /  表(6)
计量
  • 文章访问数:  653
  • HTML全文浏览量:  37
  • PDF下载量:  355
  • 被引次数: 12
出版历程
  • 收稿日期:  2022-01-03
  • 网络出版日期:  2023-05-17
  • 刊出日期:  2023-01-09

目录

/

返回文章
返回