基于无人机遥感的果园冠层氮素估算及空间分析_李达岁.pdf
《基于无人机遥感的果园冠层氮素估算及空间分析_李达岁.pdf》由会员分享,可在线阅读,更多相关《基于无人机遥感的果园冠层氮素估算及空间分析_李达岁.pdf(9页珍藏版)》请在咨信网上搜索。
1、第42卷 第4期2023年 7月华中农业大学学报Journal of Huazhong Agricultural UniversityVol.42 No.4July 2023,158166基于无人机遥感的果园冠层氮素估算及空间分析李达岁1,阮思奇1,胡青青1,张金智1,2,张亚昊3,佃袁勇1,胡春根1,2,刘永忠1,2,雷宏伟4,周靖靖11.华中农业大学园艺林学学院,武汉 430070;2.果蔬园艺作物种质创新与利用全国重点实验室,武汉 430070;3.武汉市洪山区园林局林业防护中心,武汉 430070;4.江西绿萌科技控股有限公司,赣州 341600摘要 为快速准确地获取植株冠层氮素含量及
2、空间分布特征,对大尺度的果园进行精准动态的管理,以宽行窄株小冠模式、宽行窄株篱壁模式和传统栽培模式3种栽培模式的120棵柑橘树为研究对象,通过测定冠层氮素含量并提取无人机遥感影像多光谱数据中的纹理指数和植被指数,运用随机森林算法(RF)建立基于植被指数、纹理指数以及融合植被指数和纹理指数的柑橘冠层氮素反演模型,并比较融合植被指数和纹理指数的支持向量机(SVM)、BP神经网络算法(BP)和RF的模型反演精度。结果显示:在随机森林算法中,融合植被指数和纹理指数比单独的植被指数或纹理指数更能准确预测柑橘冠层氮素含量;植被指数训练集R2为0.710,测试集R2为0.430;纹理指数训练集R2为0.76
3、1,测试集R2为0.349;融合植被指数和纹理指数训练集R2为0.775,测试集R2为0.533。融合植被指数和纹理指数在SVM算法训练集R2为0.511,测试集R2为0.371;BP神经网络训练集R2为0.651,测试集R2为0.204。用融合植被指数和纹理指数的RF模型对3种栽培模式的柑橘园进行氮素反演,得到宽行窄株小冠模式的柑橘冠层平均氮素含量最高,其次为宽行窄株篱壁模式,传统栽培模式最低,氮素含量均值分别为31.33、30.20和27.82 mg/g。结合无人机遥感与融合植被指数和纹理指数的随机森林算法能够有效预测柑橘冠层氮素含量,可为大尺度柑橘果园定量施肥提供参考。关键词 无人机遥感
4、;柑橘;精准施肥;氮素;多光谱;植被指数;纹理指数中图分类号 S127;S666 文献标识码 A 文章编号 1000-2421(2023)04-0158-09柑橘是世界第一大类水果,也是我国南方栽培面积最广、经济地位最重要的果树1。氮素是柑橘生长发育的必需元素,其盈亏直接影响柑橘的产量及果实品质2,但柑橘对氮肥的利用率相对较低3,因此,实时、快捷、准确地检测果树氮素含量并进行科学施肥,是果树优质丰产的重要保障4。传统的果树氮素含量测量方法费时、费力,不能满足大面积、快速、实时监测的需要5,且植物氮素的测定需要在专门的实验室中对叶片样品进行化学分析,费用昂贵且会产生对环境有负面影响的有毒化学物质
5、6。基于无人机(unmanned aerial vehicle,UAV)获取多光谱遥感数据作为一种无损、干净、快速的方法,可用于监测植物的氮含量7。基于无人机遥感反演农作物长势参数(形态参数、生理生化参数、胁迫参数、产量参数)的相关研究近年来发展迅速,已成为国内外农业遥感领域的热点8。基于遥感光谱数据估算植物氮含量的基本方法有 2种:基于作物地面实测参数和其遥感光谱反射率相关关系的统计模型和基于模拟光在叶片或冠层辐射传输过程的物理模型9。基于统计模型估算植物氮含量的研究包含光谱指数、原始光谱、包络线去除光谱、导数光谱等不同参数。然而,仅使用植被指数进行氮素含量估算,精度较低,植被指数结合机器学
6、习算法可适当提高果园冠层氮素反演精度10-11。高分辨率遥感影像除了具备光谱信息外,还具有丰富的纹理信息,能够有效反映植被内部的结构信息,有效缓和光谱饱和现象12。而利用纹理信息进行森林参数、作物(如冬小麦和水稻等)生物量和氮素营养诊断已有报道13-14,但在果树冠层营养元素收稿日期:2022 12 15基金项目:国家重点研发计划项目(2019YFD1000104);国家自然科学基金项目(31901963);国家柑橘产业技术体系(CARS-26)李达岁,E-mail:lds_通信作者:周靖靖,E-mail:李达岁,阮思奇,胡青青,等.基于无人机遥感的果园冠层氮素估算及空间分析 J.华中农业大学
7、学报,2023,42(4):158166.DOI:10.13300/ki.hnlkxb.2023.04.018反演及诊断方面鲜有尝试。本研究利用大疆精灵4多光谱版无人机采集江西赣州宽行窄株小冠模式、宽行窄株篱壁模式、传统栽培模式的柑橘园多光谱遥感数据,提取植被指数和纹理指数,利用随机森林算法(random forest,RF)、支持向量机(support vector machine,SVM)、BP 神经网络算法(back propagation,BP)3种机器语言算法对柑橘冠层氮素含量进行反演,遴选出适合柑橘冠层氮素反演的高精度模型,以期为大尺度柑橘果园高效精准定量施肥提供技术支撑。1材料与
8、方法1.1材料研究区位于江西省赣州市信丰县绿萌柑橘基地(2429 2709 N、11354 11638 E),属亚热带丘陵山区湿润季风气候。本研究选择宽行窄株小冠模式、宽行窄株篱壁模式和传统栽培模式的柑橘树作为研究对象15,采用分层抽样和随机抽样相结合的方式选取 120 株柑橘树,每种栽培模式 40 株,采用A8 Plus小型化RTKGPS差分定位仪对样本植株进行定位,采取每株柑橘树的冠层叶片 35 片,放入-20 冰箱保存。1.2遥感数据采集和辐射定标于 2020年 11月 16日上午 10:00晴朗无云的天气下,利用大疆精灵4多光谱版无人机对柑橘果园进行遥感影像采集,设置飞行高度为100
9、m,镜头焦距为5.74 mm。遥感影响空间分辨率为5 cm,包含5个波段:蓝波段(blue,B)、绿波段(green,G)、红波段(red,R)、红边波段(rededge band,RE)和近红外波段(near infrared band,NIR)。辐射定标:DN(digital number)值与地面表观反射率之间的回归方程为:y=0.01743x-0.23758(x 为地球表观反射率,y为DN值),得到表观反射率影像。1.3柑橘冠层氮素含量测定采用凯氏定氮法16测定柑橘叶片氮素含量,将每株柑橘树的35片叶的氮素含量均值作为其冠层叶片氮素含量。1.4植被指数计算提取17个植被指数,构建基于多
10、光谱的柑橘冠层氮素含量反演模型(表1)。表1植被指数计算公式Table 1Calculation formula of vegetation index注:R为光谱反射率:RNIR为近红外反射率、RR为红光反射率、RG为绿光反射率、RB为蓝光反射率。Note:R is spectral reflectance:RNIR is near-infrared reflectance,RR is the red light reflectivity,RG is the green light reflectivity,RB is the blue light reflectivity.第 4 期李达岁
11、 等:基于无人机遥感的果园冠层氮素估算及空间分析反演及诊断方面鲜有尝试。本研究利用大疆精灵4多光谱版无人机采集江西赣州宽行窄株小冠模式、宽行窄株篱壁模式、传统栽培模式的柑橘园多光谱遥感数据,提取植被指数和纹理指数,利用随机森林算法(random forest,RF)、支持向量机(support vector machine,SVM)、BP 神经网络算法(back propagation,BP)3种机器语言算法对柑橘冠层氮素含量进行反演,遴选出适合柑橘冠层氮素反演的高精度模型,以期为大尺度柑橘果园高效精准定量施肥提供技术支撑。1材料与方法1.1材料研究区位于江西省赣州市信丰县绿萌柑橘基地(242
12、9 2709 N、11354 11638 E),属亚热带丘陵山区湿润季风气候。本研究选择宽行窄株小冠模式、宽行窄株篱壁模式和传统栽培模式的柑橘树作为研究对象15,采用分层抽样和随机抽样相结合的方式选取 120 株柑橘树,每种栽培模式 40 株,采用A8 Plus小型化RTKGPS差分定位仪对样本植株进行定位,采取每株柑橘树的冠层叶片 35 片,放入-20 冰箱保存。1.2遥感数据采集和辐射定标于 2020年 11月 16日上午 10:00晴朗无云的天气下,利用大疆精灵4多光谱版无人机对柑橘果园进行遥感影像采集,设置飞行高度为100 m,镜头焦距为5.74 mm。遥感影响空间分辨率为5 cm,包
13、含5个波段:蓝波段(blue,B)、绿波段(green,G)、红波段(red,R)、红边波段(rededge band,RE)和近红外波段(near infrared band,NIR)。辐射定标:DN(digital number)值与地面表观反射率之间的回归方程为:y=0.01743x-0.23758(x 为地球表观反射率,y为DN值),得到表观反射率影像。1.3柑橘冠层氮素含量测定采用凯氏定氮法16测定柑橘叶片氮素含量,将每株柑橘树的35片叶的氮素含量均值作为其冠层叶片氮素含量。1.4植被指数计算提取17个植被指数,构建基于多光谱的柑橘冠层氮素含量反演模型(表1)。表1植被指数计算公式T
14、able 1Calculation formula of vegetation index名称 Name归一化差值植被指数 Normalized difference vegetation index差值植被指数 Difference vegetation index增强植被指数 Enhanced vegetation index土壤修正植被指数 Soil adjusted vegetation index非线性植被指数 None liner index绿色归一化植被指数 Green normalized difference vegetation index绿色比植被指数 Green rat
15、io vegetation index归一化叶绿素指数 Normalized pigment chlorophyll index作物氮反射指数 Nitrogen reflectance index土壤调节植被指数 Optimized soil adjusted vegetation index植物衰老反射率指数 Plant senescing reflectance index比植被指数 Ratio vegetation index冠层结构不敏感植被指数 Structure insensitive pigment index三角植被指数 Triangle vegetation index可见光
16、大气阻抗植被指数 Visible atmospherically resistant vegetation index宽范围动态植被指数 Wide dynamic vegetation index转化叶绿素吸收反射指数Transformation chlorophyll absorption reflection index公式 FormulationINDV=(RNIR-RR)/(RNIR+RR)IDV=(RNIR-RR)IEV=2.5()RNIR-RR/(RNIR-7.5RB+6RR+1)ISAV=1.5(RNIR-RR)(RNIR+RR+0.5)INL=(RNIR2-RR)/(RNIR2
17、+RR)IGNDV=(RNIR-RG)/(RNIR+RG)IGRV=RNIR/RGINPC=(RR-RB)/(RR+RB)INR=RR/(RR+RG+RB)IOSAV=1.16(RNIR-RR)/(0.16+RNIR+RR)IPSR=(RR-RB)/RNIRIRV=RNIR/RRISIP=(RNIR-RB)/(RNIR-RR)ITV=0.5 120(RR-RG)-200(RR-RG)IVAR=(RG-RR)/(RR+RG-RB)IWDRV=(0.12RNIR-RR)/(0.12RNIR+RR)ITCAR=3(RNIR-RR)-0.2(RNIR-RG)(RNIR/RR)参考文献Reference
18、1717171718171719171720172021171821注:R为光谱反射率:RNIR为近红外反射率、RR为红光反射率、RG为绿光反射率、RB为蓝光反射率。Note:R is spectral reflectance:RNIR is near-infrared reflectance,RR is the red light reflectivity,RG is the green light reflectivity,RB is the blue light reflectivity.159第 42 卷 华 中 农 业 大 学 学 报1.5纹理指数提取采用灰度共生矩阵方法分别提取红光
19、波段、绿光波段、蓝光波段、近红外波段和红边波段的8种纹理指数22,分别为:同质性(IHOM)、平均值(IMEAN)、对比度(ICON)、熵(IEnt)、非相似度(IDIS)、变化量(IVAR)、角二阶矩阵(ISM)和相关性(ICOR)(表2)。1.6数据处理根据决定系数R2、均方误差(mean square error,MSE)、平均绝对值误差(mean absolute error,MAE)来评价反演精度。1.7模型建立利用Python软件中的RF、SVM、BP算法程序包选取70%的样本作为训练集,30%样本作为测试集验证模型精度,分别建立基于植被指数、纹理指数、融合植被指数和纹理指数的柑橘
20、冠层氮素反演模型。2结果与分析2.1不同栽培模式下柑橘冠层氮素含量分布特征由表3可见,3种栽培模式的柑橘冠层的氮素含量由高到低为:宽行窄株小冠宽行窄株篱壁传统栽培模式,平均值分别为31.57、30.10、27.64 mg/g。2.2基于植被指数的随机森林柑橘冠层氮素反演利用基于特征递归消除的机器学习算法,对提取的17个植被指数进行重要性排序,递归重要性最低的指数,筛选至最佳的3个指数并以此确定模型。模型的训练集 R2为 0.710,MSE 和 MAE 分别为0.311和0.449;测试集R2为0.430,MSE和MAE分别为0.467和0.555。重要性高到低排序分别为归一化叶绿素指数(INP
21、C)、冠层结构不敏感植被指数(ISIP)、作物氮反射指数(INR)。2.3基于纹理指数的随机森林柑橘冠层氮素反演对提取的40个纹理指数进行重要性排序,递归消减重要性最低的指数,筛选至最佳的7个指数并以此确定模型(图1)。模型的训练集R2为0.761,MSE和 MAE 分 别 为 0.256 和 0.389;测 试 集 的 R2为0.349,MSE 和 MAE 分别为 0.533和 0.599。重要性最高的为绿光波段的平均值(IMEAN_G),其次分别是近红外波段的相关性(ICOR_NIR)、红边波段的相关性(ICOR_RE)。2.4融合植被指数和纹理指数的随机森林柑橘冠层氮素反演对提取的57个
22、融合指数进行重要性排序,递归消减重要性最低的指数,筛选至最佳的10个指数数量并以此确定模型(图2)。模型训练集R2为0.775,MSE 和 MAE分别为 0.241和 0.387;测试集的 R2为0.533,MSE 和 MAE 分别为 0.382 和 0.497。其中重要性最高的绿光波段的平均值(IMEAN_G),其次是归一 化 叶 绿 素 指 数(INPC)、红 边 波 段 的 相 关性(ICOR_RE)。2.5融合植被指数和纹理指数的SVM和BP神经网络柑橘冠层氮素反演使用 SVM 和 BP 算法以叶片氮素和筛选的 10个重要性排序靠前的参数建立模型。结果显示,表3不同模式柑橘冠层的氮素含
23、量分布Table 3Distribution of nitrogen content in canopy of different patternsmg/g 模式Pattern宽行窄株篱壁 Wide row and narrow plant fence cultivation宽行窄株小冠 Wide row and narrow plant cultivation传统栽培模式 Traditional cultivation 平均值Mean30.1031.5727.64最大值Max35.6835.9936.02最小值Min23.7625.6618.51中值Median30.4032.2527.93
24、标准差Standard deviation2.182.013.07表2纹理指数计算公式Table 2Calculation formula of texture features index名称 Name同质性 Homogeneity(IHOM)平均值 Mean(IMEAN)对比度 Contrast(ICON)熵 Entropy(IEnt)非相似度 Dissimilarity(IDIS)变化量 Variance(IVAR)二阶矩阵 Second moment(ISM)相关性Correlation(ICOR)公式 FormulationIHOM=i,j=0N-1Pi,j/1+(i-j)2IMEA
25、N=i,pN-1ipICON=i,j=0N-1iPi,j(i-j)2IEnt=i,j=0N-1iPi,j(-lnPi,j)IDIS=i,j=0N-1iPi,j|i-j|IVAR=i,j=0N-1(Pi,j-)2/(N-1)ISM=i,j=0N-1iP2i,jICOR=i,j=0N-1ijPi,j-121.1=i=0N-1ii=0N-1Pi,j2.2=j=0N-1ij=0N-1Pi,j注:i,j为像素灰度、N 为灰度级数、Pi,j是在给定空间距离和方向时,灰度以i(行)为起始点,出现灰度级j(列)的概率。Note:i,j is pixel gray,N is gray level,Pi,j is
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 无人机 遥感 果园 氮素 估算 空间 分析 李达岁
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。