基于无人机RGB影像的马铃薯植株钾含量估算.pdf

2023年7 月第54卷第7 期农报学业机械doi:10.6041/j.issn.1000-1298.2023.07.019基于无人机RGB影像的马铃薯植株钾含量估算马彦鹏边明博1樊意广1陈志超杨贵军冯海宽1.3（1.北京市农林科学院信息技术研究中心，北京10 0 0 97；2.河南理工大学测绘与国土信息工程学院，焦作4540 0 0；3.南京农业大学国家信息农业工程技术中心，南京2 10 0 95）摘要：马铃薯植株钾含量（Plantpotassiumcontent,PKC）是监测马铃薯营养状况的重要指标，快速准确地获取马铃薯植株钾含量对田间施肥和生产管理具有指导意义。基于无人机遥感平台搭载RCB传感器分别获取马铃薯块茎形成期、块茎增长期和淀粉积累期的RGB影像，并实测马铃薯植株钾含量。首先利用各个生育期的RCB影像提取每个小区冠层平均光谱和纹理特征。然后分别基于冠层光谱和纹理特征构建植被指数和纹理指数（NDTI、RT I 和DTI），并与实测PKC进行相关性分析。最后利用多元线性回归（Multiplelinearregression，M LR）、偏最小二乘（Pa r t i a l l e a s t s q u a r e s r e g r e s s i o n，PLSR）和人工神经网络（Artificial neural networks，A NN）构建马铃薯PKC估算模型。结果表明：各生育期NDTI、RT I 和DTI与马铃薯PKC相关性均高于单一纹理特征，植被指数结合纹理指数均能提高模型的可靠性和稳定性，MLR和PLSR构建的估算模型精度均优于ANN。本研究可为马铃薯PKC监测提供科学参考。关键词：马铃薯；植株钾含量；纹理指数；RGB影像；冠层光谱特征中图分类号：S25文献标识码：A文章编号：10 0 0-12 98（2 0 2 3)0 7-0 196-0 8OSID:Estimation of Potassium Content of Potato Plants Based onUAV RGB ImagesMA YanpengFAN YiguangYANG GuijunFENG Haikuan1BIAN MingboCHEN Zhichao1,3(1.Information Technology Research Center,Beijing Academy of Agriculture and Forestry Sciences,Beijing 100097,China2.School of Geomatics and Land Information Engineering,Henan Polytechnic University,Jiaozuo 454000,China3.National Engineering and Technology Center for Information Agriculture,Nanjing Agricultural University,Nanjing 210095,China)Abstract:Plant potassium content(PKC)of potato plants is an important indicator for monitoring potatonutrition status.Obtaining PKC quickly and accurately has guiding significance for field fertilization andproduction management.RGB images of potato plants during the tuber formation period,tuber growthperiod,and starch accumulation period were obtained by using an unmanned aerial vehicle(UAV)remote sensing platform equipped with an RGB sensor,and PKC was measured.Firstly,the averagespectral and texture features of each plot were extracted from the RGB images of each growth period.Then vegetation indices and texture indices(NDTI,RTI,and DTI)were constructed based on thespectral and texture features of the canopy,and their correlations with the measured PKC were analyzed.Finally,multiple linear regression(MLR),partial least squares regression(PLSR),and artificial neuralnetworks(ANN)were used to construct models for estimating potato PKC.The results showed that thecorrelations between NDTI,RTI,DTI and PKC were higher than those of single texture features duringeach growth period.Combining vegetation and texture indices can improve the reliability and stability ofthe model.MLR and PLSR models were superior to ANN.The research result can provide scientificreferences for monitoring PKC in potato plants.Key words:potato;plant potassium content;texture index;RCB images;canopy spectral features收稿日期：2 0 2 3-0 4-0 3修回日期：2 0 2 3-0 5-0 8基金项目：黑龙江省揭榜挂帅科技攻关项目（2 0 2 1ZXJ05A05）和国家自然科学基金项目（416 0 1346）作者简介：马彦鹏（1997 一），男，实习研究员，主要从事农业定量遥感研究，E-mail：ma y a n p e n g 2 0 2 2 16 3.c o m通信作者：冯海宽（198 2 一），男，高级工程师，主要从事农业定量遥感研究，E-mail：f e n g h a i k u a n 12 3 16 3.c o m第7 期197马彦鹏等：基于无人机RGB影像的马铃薯植株钾含量估算0引言马铃薯是重要的薯类作物，可以粮菜兼用，在我国广泛种植。马铃薯是喜钾作物，钾含量会直接影响马铃薯产量、株高、茎粗、叶面积和生物量等2-3。因此快速准确地获取马铃薯植株钾含量（Pl a n t p o t a s s i u m c o n t e n t,PK C）对监测马铃薯生长发育状况具有重要意义。传统测量PKC的方法以田间破坏性取样为主，测量成本高且费时费力，难以大面积推广使用。无人机遥感技术机动性强，可以用于农作物长势监测、产量估算和病虫害监测等，使田间生产管理更加高效便捷4-6 。基于无人机遥感平台可以搭载RCB、多光谱和高光谱传感器7-8 。LU等9 基于无人机高光谱传感器利用非负矩阵分解提取植被光谱，建立偏最小二乘回归（Partialleast squaresregression，PLSR）模型。结果表明，绿光带和红边带对水稻植株钾积累量估算有显著贡献。THOMSON等【10 研究发现利用无人机高光谱传感器反射率建立的PLSR模型可以有效估算森林植物叶片钾含量。SEVERTSON等 利用无人机多光谱和高光谱反射率有效识别了缺钾油菜。以上研究主要集中在高光谱传感器上，高光谱虽然能获取更多的光谱信息，但因其价格昂贵，处理过程复杂限制了它在农业生产中的应用。RGB相机传感器虽然只能获取3个波段的光谱信息，但因其可以获取高分辨率影像，而且价格便宜，数据处理过程较为简单，成为广受关注的数据获取方式。樊意广等【12 基于无人机RGB相机传感器利用多元线性回归（Multiplelinearregression，M LR）、BP神经网络和Lasso回归结合冠层光谱特征模型有效估算了马铃薯植株氮含量，结果表明MLR构建的模型最优。刘杨等13 基于无人机RCB影像利用MLR、支持向量机和人工神经网络（Artificialneuralnetworks，A NN)结合植被指数估算马铃薯地上生物量，结果表明各生育期利用MLR方法构建的马铃薯地上生物量估算模型效果最佳。陶惠林等14 基于无人机RGB影像有效估算了冬小麦生物量。上述研究结果表明利用无人机RGB相机传感器可以有效估算作物氮素含量和生物量，但是仅利用影像的光谱信息构建的模型容易出现模型精度饱和现象。纹理信息通过像素灰度分布变化来反映目标与环境之间的空间关系，可以提供新的冠层信息15-16 ，在作物生长监测中发挥重要作用17-18 ,有学者提出引人图像的纹理信息结合光谱信息构建模型LU等19 从无人机RCB影像中分别提取光谱信息和纹理信息估算水稻植株钾积累量，结果表明光谱信息和纹理信息可以准确估算水稻植株钾积累量,两者融合能够进一步提高精度。ZHENG等2 0 获取水稻关键生育期冠层多光谱图像，将提取的纹理特征归一化结合植被指数建模估算水稻地上生物量,结果表明，与单独使用光谱信息相比，归一化纹理信息与光谱信息结合显著提高了水稻生物量估算的准确性。MA等2 1 利用棉花收获前冠层RGB影像监测棉花产量，表明从无人机高分辨率RGB影像中提取的植被指数和纹理特征与棉花产量显著相关，且植被指数结合纹理特征估算棉花产量效果最好。马铃薯的生长发育状况与上述作物不同，前期茎叶生长旺盛，但是到生长后期地上干物质不断向地下转移，地上叶片和茎开始死亡。基于无人机RGB影像提取植被指数和灰度共生矩阵（Gray-levelco-occurrence matrix,GLCM）纹理特征可以提高高覆盖度下估算马铃薯地上生物量的准确性2 2】,钾与氮素和生物量机理不同，已有的研究方法是否适用于马铃薯PKC监测还有待进一步研究。因此，本研究以无人机遥感平台搭载RCB传感器获取马铃薯3个关键生育期的RGB影像，并提取各个生育期的冠层光谱特征和纹理信息，结合地面实测的PKC数据，采用3种不同的方法对每个生育期马铃薯PKC进行建模和验证，探究光谱特征结合纹理信息估算马铃薯PKC的适用性。1材料与方法1.1实验设计实验于2 0 19年在北京市昌平区小汤山镇国家精准农业研究示范基地开展。采用的马铃薯品种为Z1（中薯5)和Z2（中薯3）。共设3个实验区，每个实验区3个重复。P区（密度实验区）设置3个密度水平（T0:60000株/hm；T 1:7 2 0 0 0 株/hm；T 2：84000株/hm）；N区（氮肥实验区）设置4个氮肥水平（NO:0 kg/hm；N1:2 44.6 5 k g/h m；N2:489.15kg/hm;N3：7 33.5k g/h m）；K 区（钾肥实验区）设置3个钾肥水平（K0:0kg/hm;K1:970kg/hm；K2:1941kg/hm）。P区和N区钾肥水平均为Kl：970kg/hm。共48 个小区（图1），每个小区面积为5 m 6.5 m。1.2无人机RGB影像获取及预处理采用大疆精灵4Pro无人机遥感平台，分别在块茎形成期、块茎增长期和淀粉积累期选择晴朗无云的天气，在光照强度稳定的12:0 0 14:0 0 获取RGB影像。表1为无人机飞行参数设置。飞行结束后将图像导人三维建模软件AgisoftPhotoScanProfessional(64bit）生成每个生育期数字正射影像。198农2023年机业报械学11600E11700EN重复1重复2重复3重复1重复2重复3T2TOT1T2TOT1K6(P03)(P04)(P09)K5(P10)(P15)(P16)K4T1T2TOT1T2TOP区(P02)(P05)(P08)(P11)(P14)(P17)TOT1T2TOT1T2(PO1)(PO6)(P07)(P12)(P13)(P18)K7K10K3N3NON1N3NON1N(04)N(05)N(12)N(13)N(20)N(21)N,0.0.0tN.0,000tN2N3NON2N3NON(03)N(06)N(11)N(14)N(19)N(22)N区N1N2N3N1N2N3N(02)N(07)N(10)N(15)N(18)N(23)NON1N2NON1N2K8N(01)N(08)N(09)K9N(16)N(17)N(24)K2K2K2K2实验基地K(01)K(02)K(03)K区昌平区地面控制点KOKOKOK11K(06)K(05)K(04)N.,0.0.66N,0,0060.20406080K1北京市各区km11600E11700E图1马铃薯田间位置及实验设计Fig.1Potato field location and experimental design表1无人机飞行参数Tab.1UAVflightparameters参数数值/形式无人机质量/g1391最大飞行速度/(kmh-)50最大起飞海拔高度/m6000飞行时间/min30相机有效像素数2000万波段红(R)、绿(G)、蓝(B)航向重叠度/%85旁向重叠度/%85飞行高度/m20主要步骤包括对齐图像、建立密集点云、生成格网、生成纹理、建立正射、导出数字正射影像1.3地面数据获取在无人机数据采集结束后进行地面数据采集。马铃薯PKC测量方法为田间采样和实验室化学分析相结合。分别于马铃薯块茎形成期、块茎增长期和淀粉积累期在每个实验小区选取3株具有代表性的植株带回实验室进行茎叶分离，用清水冲洗干净后在干燥箱内10 5杀青30 min，然后将温度调至80干燥48 h以上，质量恒定后进行称量，得到各器官的干质量。在实验室内利用电感耦合等离子体发射光谱仪（iCAP6300型）测量各器官钾含量。最后计算PKC，计算公式为CLkMLD+CskMsDPKC=100%(1)MLD+MsD式中Cik叶片钾含量，%CsK地上茎钾含量，%MLD叶片干质量，8MsD地上茎干质量，g1.4光谱指数选取从无人机RGB影像中提取每个小区的冠层光谱特征，然后根据现有的研究进展选取2 4个能够描述作物氮素和钾素营养状况的光谱植被指数对马铃薯块茎形成期、块茎增长期、淀粉积累期的PKC进行估算。植被指数如表2 所示。表2植被指数Tab.2Vegetation indices植被指数表达式文献序号植被指数表达式文献序号RMCRVI(g2-r2)/(g2+r)12GRGBVI(g?-br)/(g?+br)12BNDI(T-g)/(r+g+0.01)12R/(R+G+B)GLA(2g-r+b)/(2g+r+b)12G/(R+G+B)CIVE0.441r-0.881g-0.3856b+18.7874512bB/(R+B+G)(g-b)/(r-g)(g-b)/(r-g)12g+b8+6EXG2g-b-r19g-b8-6EXR1.4r-g19r-bT-bEXCREXG-EXR19r+bR+bIKAW(r-b)/(r+b)19I/gr/gCRVI(g-r)/(g+r)19(r-g-b)/(r+g)(r-g-b)/(r+g)12VARI(g-r)/(g+r-b)19注：R、G、B分别表示红、绿、蓝波段的中心波长。199马彦鹏等：基于无人机RCB影像的马铃薯植株钾含量估算第7 期1.5纹理特征选取遥感影像的纹理特征可以用来描述图像中的空间分布和灰度分布。GLCM是用来描述图像中像素灰度值之间关系的一种统计方法。GLCM纹理特征包括对比度（Contrast，c o n）、相关性（Correlation，cor）、差异性（Dissimilarity，d i s）、熵（Entropy，e n t）、同质性（Homogenetity，h o m）、二阶矩（Secondmoment,sm）、方差（Variance，v a r）和均值（mean）。有研究表明提取纹理时移动窗口大小和方向对结果影响较小，可以忽略2 3-2 4。本研究基于3个生育期的RGB影像选择33的窗口，提取每个波段45方向的8 个纹理特征，共2 4个。为了提高纹理指数与马铃薯PKC之间的相关性，本文在单个纹理特征的基础上按照归一化植被指数（NDVI）、差值植被指数（DVI）和比值植被指数（RVI)方法构建了3个纹理指数（TIS）一NDTI、D T I 和RTI。选出最佳的组合方式用来估算马铃薯PKC，计算公式为T,-T,NDTI:(2)二T,+T,DTI=T,-T,(3)工RTI(4)T式中T、T,3个波段中随机纹理特征值1.6数据分析方法将重复1和重复3的32 组数据作为训练集建立马铃薯PKC的估算模型，将重复2 的16 组数据作为验证集，验证模型精度。采用的数据分析方法为MLR、偏最小二乘（Partial least squares regression,PLSR）和ANN。M LR是一种基本的回归分析方法，它建立了变量和目标之间的线性关系，并且假设误差项是独立同分布的。PLSR是一种经典的线性回归方法，它可以解决多重共线性问题，将原始变量转换为一组新的特征，并且确保每个成分都与目标相关联，通过保留这些成分的不同数量，可以进行多级回归分析，从而实现预测。ANN是一种基于神经元结构的非线性建模方法，通过学习数据中的复杂模式进行建模，并且可以根据数据的结构和分布来自适应地调整模型参数。1.7模型评价指标为评价模型预测结果的准确性，采用决定系数（Co e f f i c i e n t o f d e t e r mi n a t i o n，R）、均方根误差（Rootmeansquare error，RM SE）和标准均方根误差(Normalized root mean square error,NRMSE）作为模型精度评价指标2结果与分析2.1相关性分析采用Pearson相关性分析了不同模型参数与马铃薯PKC之间的相关性。表3为所有模型参数与3个生育期马铃薯PKC的相关系数。块茎形成期除IKAW和（g-b)/（r-g）与马铃薯PKC表现不显著相关外,其余植被指数均与马铃薯PKC表现显著性相关。其中G、b 和r-b表现为0.0 5水平显著性相关,R、B、r、g、g+b、g-b、r +b、r/g、EXG、EXR、EXGR、G RVI、VA RI、（r-g-b)/(r+g)、M G RVI、RGBVI、ND I、G LA 和CIVE共19个植被指数均表现为0.0 1水平显著性相关。相关系数绝对值最大的植被指数是（r-g-b)/（r+g），为0.8 5。块茎增长期G和（g-b)/（r-g）与马铃薯PKC不显著相关，R、B、r、g、b、g+b、g-b、r-b、r+b、r/g、EXG、EXR、EXGR、I K A W、G RV I、V A RI、(r-g-b)/(r+g)、表3马铃薯3个生育期植被指数与PKC的相关系数Tab.3Correlation coefficient between vegetation index and PKC in three growth periods of potatoPKCPKC植被指数植被指数块茎形成期块茎增长期淀粉积累期块茎形成期块茎增长期淀粉积累期R-0.54*0.57*-0.76*EXR0.72*0.82*0.65*G-0.31*-0.19-0.41*EXCR0.67*0.79*0.63*B-0.42*-0.60*-0.64*IKAW0.020.59*0.34*-0.84*-0.77*-0.68*CRVI0.70*0.80*0.64*0.59*0.75*0.59*VARI0.71*0.80*0.64*b0.30*0.65*-0.46*(r-g-b)/(r+g)-0.85*-0.54*-0.65*g+b0.84*0.77*0.68*MGRVI0.70*0.80*0.64*g-b0.48*0.71*0.54*RGBVI0.53*0.72*0.56*r-b-0.33*0.48*0.12NDI-0.70*-0.80*-0.64*I+b-0.59*-0.75*-0.59*CLA0.78*0.85*0.67*/g-0.69*-0.80*-0.64*CIVE0.73*-0.83*-0.65*EXG0.59*0.75*0.59*(g-b)/(r-g)0.17-0.24-0.34*注：*表示0.0 5水平显著性相关，*表示0.0 1水平显著性相关。下同。200农2023年报机业学械MGRVI、RG BVI、ND I、G LA 和CIVE共2 2 个植被指数与马铃薯PKC表现为0.0 1水平显著性相关。相关系数绝对值最大的植被指数是GLA，为0.8 5。淀粉积累期r-b与马铃薯PKC为不显著相关,IKAW和（g-b)/（r-g）表现为0.0 5水平显著性相关,R、G、B、r、g、b、g+b、g-b、r +b、r/g、EXG、EXR、EXGR、G RVI、VA RI、（r -g-b)/(r +g)、M G RVI、RGBVI、ND I、G LA 和CIVE共2 1个植被指数与马铃薯PKC表现为0.0 1水平显著性相关。相关系数绝对值最大的植被指数是R，为0.7 6。如表4所示，单一纹理特征在块茎形成期R-ent和B-hom与马铃薯PKC不显著相关，G-mean和G-ent表现为0.0 5水平显著性相关,R-mean、R-var、R-h o m、R-c o n、R-d i s、R-s e c、R-c o r、G-var、G-h o m、G-c o n、G-d i s、G-s e c、G-c o r、B-mean、B-v a r、B-c o n、B-d i s、B-e n t、B-s e c 和 B-cor共2 0 个纹理特征与马铃薯PKC呈0.0 1水平显著性相关。相关系数绝对值最大的是G-con,为0.79。块茎增长期G-mean表现为不显著相关外，其余纹理指数均表现为显著性相关。其中R-cor和G-cor表现为0.0 5水平显著性相关,R-mean、R-var、R-h o m、R-c o n、R-d i s、R-e n t、R-s e c、G-var、G-h o m、G -c o n、G -d i s、G -e n t、G -s e c、B-mean、B-v a r、B-h o m、B-c o n、B-d i s、B-e n t、B-sec和Bc o r 共2 1个单一纹理特征与马铃薯PKC呈0.0 1水平显著性相关。相关系数绝对值最大的是R-var和R_dis,为0.8 2。淀粉积累期G-cor、B-var、B-h o m、B-c o n、B-d is 和R-cor与马铃薯PKC不显著相关。R-mean、R-v a r、R-h o m、R-con、R-d i s、R-e n t、R-s e c、G -m e a n、G -v a r、G -hom、G-c o n、G-d i s、G-e n t、G -s e c、B-m e a n、B-ent、Bs e c 和B-cor共18 个单一纹理特征表现为0.01水平显著性相关。相关系数绝对值最大的是R-mean,为0.7 6。构建的3个TIS在3个生育期均表现为0.0 1水平显著性相关，相关系数绝对值均大于表现最优的单一纹理特征相关系数绝对值，且3个生育期NDTI、D T I 和RTI相关系数绝对值相差不大，相关系数绝对值最大为0.8 4、0.8 3和0.7 8。表4马铃薯3个生育期纹理特征与PKC的相关系数Tab.4Correlation coefficient of texture characteristics and PKC in three growth stages of potatoPKCPKC纹理特征纹理特征块茎形成期块茎增长期淀粉积累期块茎形成期块茎增长期淀粉积累期R-mean-0.54*0.57*-0.76*G-sec0.45*-0.73*-0.48*R-var0.78*0.82*0.40*G-cor-0.75*-0.29*0.01R-hom-0.41*-0.78*-0.47*B-mean-0.42*-0.60*-0.64*R-con0.78*0.81*0.39*B-var0.66*0.79*-0.17R-dis0.75*0.82*0.43*B-hom-0.11-0.71*0.11R-ent-0.240.74*0.49*B-con0.71*0.79*-0.15R-sec0.41*-0.73*-0.49*B-dis0.68*0.78*0.03R-cor-0.75*-0.29*0.14B-ent0.390.64*-0.35*G-mean-0.31*-0.19-0.41*B-sec0.48*-0.60*0.42*G-var0.79*0.81*0.39*B-cor-0.72*-0.42*-0.44*G-hom-0.38*-0.783*0.47*NDTI-0.82*-0.833*0.77*G-con0.79*0.80*0.38*RTI0.84*0.83*-0.78*G-dis0.75*0.81*0.43*DTI-0.820.82*-0.77*G-ent-0.29*0.74*0.49*注：R-mean表示R的mean纹理特征，其他类推。2.2马铃薯PKC估算与验证分别以3个生育期相关系数绝对值最大的3个植被指数和构建的3个最优的TIS为输人变量,结合MLR、PLSR和ANN构建马铃薯PKC的估算模型。表5、6 分别为3个生育期MLR、PLSR和ANN3种方法以植被指数（VIS）和植被指数结合构建的TIS为输人变量建模和验证的结果。3个生育期利用同种建模方法植被指数结合TIS作为输人变量估算PKC的效果优于仅将植被指数作为输入变量的估算效果。3种建模方法均表现为块茎形成期和块茎增长期的估算效果优于淀粉积累期；同种建模方法从块茎形成期到淀粉积累期模型验证精度逐渐变差。不同的估算模型在3个生育期植被指数结合TIS的建模和验证效果如图2 所示。3个生育期MLR和PLSR的建模和验证结果均优于ANN。不同建模方法中均是块茎形成期和块茎增长期估算效果较好，各样本点分布在1:1线附近。块茎形成期MLR、PLSR和ANN建模R分别为0.8 2、201马彦鹏等：基于无人机RGB影像的马铃薯植株钾含量估算第7 期表53个生育期3种方法马铃薯PKC估算结果对比Tab.5Comparison of potato PKC estimation results of three methods in three growth periodMLRPLSRANN生育期变量R2RMSE/%NRMSE/%R?RMSE/%NRMSE/%R2RMSE/%NRMSE/%VIS0.740.3812.440.670.4212.200.610.5317.35块茎形成期VIS.+TIS0.820.3210.410.800.339.590.750.4013.02VIS0.790.3611.770.760.3912.690.690.4815.72块茎增长期VIS+TIS0.840.3110.070.800.3611.600.710.4815.63VIS0.680.3313.770.680.3313.770.630.3615.34淀粉积累期VIS+TIS0.700.3113.160.700.3213.360.660.3414.15表6 3个生育期3种方法马铃薯PKC验证结果对比Tab.6Comparison of potato PKC verification results of three methods in three growth periodsMLRPLSRANN生育期变量R2RMSE/%NRMSE/%R2RMSE/%NRMSE/%R2RMSE/%NRMSE/%VIS0.860.289.230.840.307.970.680.5518.15块茎形成期VIS+TIS0.860.3110.130.880.297.790.750.4615.26VIS0.640.4715.150.680.4213.720.540.5818.76块茎增长期VIS+TIS0.680.4313.970.720.4113.350.540.5417.75VIS0.560.4720.020.560.4820.020.550.4820.29淀粉积累期VIS+TIS0.600.4519.020.580.4619.330.570.4719.730.80和0.7 5，RMSE分别为0.32%、0.33%和0.40%,NRMSE分别为10.41%、9.59%和13.02%。验证R分别为0.8 6、0.8 8 和0.7 5,RMSE分别为0.31%、0.2 9%和0.46%，NRMSE分别为10.13%、7.7 9%和15.2 6%。块茎增长期MLR建模精度最高R为0.8 4,RMSE为0.31%，NRMSE为10.07%。PLSR建模R、RM SE和NRMSE分别为0.80、0.36%和11.6 0%。ANN建模R、RM SE和NRMSE分别为0.7 1、0.48%和15.6 3%。验证集上PLSR表现最优，R为0.7 2，RMSE为0.41%，NRMSE为13.35%。MLR建模R、RM SE和NRMSE分别为0.6 8、0.43%和13.97%.ANN建模R、RM SE和NRMSE分别为0.54、0.54%和17.75%。淀粉积累期3个模型中MLR估算精度最高，建模R为0.7 0，RMSE为0.31%、NRMSE为13.16%。PLSR建模R、RM SE和NRMSE分别为0.70、0.32%和13.36%。ANN建模R、RM SE和NRMSE分别为0.6 6、0.34%和14.15%。MLR验证R、RM SE和NRMSE分别为0.6 0、0.45%和19.02%。PLSR验证R、RM SE和NRMSE分别为0.58、0.46%和19.33%。ANN验证R、RM SE和NRMSE分别为0.57、0.47%和19.7 3%。3讨论3.1植被指数估算马铃薯PKC本研究基于无人机遥感平台获取了马铃薯块茎形成期、块茎增长期和淀粉积累期的RGB影像，基于冠层光谱特征构建植被指数并与马铃薯3个生育期的PKC进行相关性分析，选取相关系数绝对值最高的3个植被指数建模估算马铃薯PKC，块茎形成期和块茎增长期马铃薯PKC估算效果较好,其原因是马铃薯前期以茎叶生长为主，到块茎形成期马铃薯地上茎叶已经基本发育完整，地上茎叶生长开始缓慢，块茎增长期马铃薯已由单纯的植株营养生长过渡到植株营养生长、生殖生长和物质积累同时进行的时期，此阶段是马铃薯整个生长发育的鼎盛时期，地上茎叶生长茂盛2 5。这两个生育期提取的光谱特征受土壤影响较小，构建的植被指数能够较好地反映马铃薯PKC变化情况。经过对比发现3个生育期通过3种方法建模时淀粉积累期精度稍差。这是因为在淀粉积累期马铃薯地上茎叶停止生长，地下块茎基本停止生长，此阶段以淀粉积累为主，淀粉和蛋白质等营养物质不断增加。地上部营养物质开始向地下转移，地上茎叶开始变黄枯萎，地面开始裸露，导致提取的光谱信息受到土壤的影响，构建的植被指数与马铃薯PKC 相关性降低2 6 3.2植被指数结合TIS估算马铃薯PKC纹理信息可以增加无人机图像的数据维数，经对比发现由单个纹理特征组合得到的TIS与马铃薯PKC之间的相关系数绝对值均大于单个纹理特征。这是因为TIS可以通过平滑冠层结构降低土壤背景、太阳高度角和传感器视角的影响2 7 。本研究基于3个生育期的植被指数结合TIS构建马铃薯PKC的估算模型。发现植被指数结合TIS构建的模型精农2022023年机报学业械5.5验证5.5y=0.80 x+0.71验证5.5y=0.92x+0.35验证y=0.82x+0.64建模R=0.80建模R2=0.75建模5.0R2=0.82-1:1线5.0-1:1线5.0.-:1线%RMSE为0.33%RMSE为0.40%RMSE为0.32%NRMSE为13.0 2%4.5NRMSE为10.41%4.5NRMSE为9.59%4.5上4.04.04.03.53.53.5y=0.76x+0.78y=0.80 x+0.64y=0.80 x+0.483.03.03.0R2=0.86R2=0.88R2-0.752.5RMSE为0.31%2.5RMSE为0.2 9%2.5RMSE为0.46%NRMSE为7.7 9%NRMSE为15.2 6%NRMSE为10.13%2.02.53.03.54.04.55.05.52.02.53.03.54.0 4.55.05.52.02.53.03.54.04.55.05.5实测马铃薯PKC/%实测马铃薯PKC/%实测马铃薯PKC/%(a)块茎形成期,MLR(b)茎形成期,PLSR(c)块茎形成期,ANN验证验证验证5.0y=0.84x+0.47建模5.0y=0.80 x+0.62建模5.0y=0.69x+1.170建模R2=0.84R2=0.714.5.-1:1线4.5R?=0.80-1:1线4.5-1:1线RMSE为0.31%RMSE为0.36%RMSE为0.48%4.0NRMSE为10.0 7%4.0NRMSE为11.6 0%4.0NRMSE为15.6 3%3.53.53.53.03.03.0y=0.59x+1.172.5y=0.53x+1.36y=0.47x+2.012.5R2=0.68R2=0.722.5R2=0.54RMSE为0.43%2.0RMSE为0.41%RMSE为0.54%2.02.0NRMSE为13.97%1.5NRMSE为13.35%NRMSE为17.7 5%1.52.02.53.03.54.04.55.01.52.02.53.03.54.04.55.01.52.02.53.03.54.04.55.0实测马铃薯PKC/%实际马铃薯PKC/%实际马铃薯PKC/%(d)块茎增长期,MLR(e)块茎增长期,PLSR(f)块茎增长期,ANN4.5y=0.70 x+0.939R2=0.70验证建模4.5y=0.70 x+0.967验证建模4.5y=0.66x+1.10R2=0.66验证建模-1:1线R2=0.70-1:1线%RMSE为0.34%-1:1线4.0RMSE为0.31%4.04.0NRMSE为13.16%RMSE为0.32%NRMSE为14.15%NRMSE为13.36%3.53.53.53.03.03.0.y=0.70 x+0.80.y=0.70 x+0.78y=0.70 x+0.792.5R2=0.602.5R2=0.582.5R=0.57RMSE为0.45%RMSE为0.46%RMSE为0.47%2.0NRMSE为19.0 2%2.0NRMSE为19.33%2.0NRMSE为19.7 3%2.02.53.03.5 4.04.52.02.53.03.54.04.52.02.53.03.54.04.5实际马铃薯PKC/%实际马铃薯PKC/%实测马铃薯PKC/%(g)淀粉积累期,MLR(h)淀粉积累期,PLSR(i)淀粉积累期,PLSR图2不同估算模型下3个生育期植被指数结合TIS建模和验证效果Fig.2Modeling and verification