高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于高光谱激光雷达的绿萝叶绿素3维重建

汪慧民 邵慧 孙龙 李伟 王程 陈杰 朱家兵

引用本文:
Citation:

基于高光谱激光雷达的绿萝叶绿素3维重建

    作者简介: 汪慧民(1997-),男,硕士研究生,现主要从事激光雷达信号处理的研究.
    通讯作者: 邵慧, shaohui@ahjzu.edu.cn
  • 基金项目:

    安徽省高校省级自然科学研究项目 KJ2021JD16

    安徽省高校协同创新项目 GXXT-2021-028

    安徽省高校省级自然科学研究项目 KJ2021A0622

    安徽省高校省级自然科学研究项目 KJ2020A0471

    安徽省科技厅面上项目 2008085MF182

    安徽省住房城乡建设科学技术计划资助项目 2020-YF22

    安徽省住房城乡建设科学技术计划资助项目 2022-YF077

    安徽省科技重大专项基金资助项目 202003a05020031

  • 中图分类号: TP721.2

3-D reconstruction of chlorophyll content of epipremnum aureum based on hyperspectral LiDAR

    Corresponding author: SHAO Hui, shaohui@ahjzu.edu.cn ;
  • CLC number: TP721.2

  • 摘要: 为了准确重建正反面叶片叶绿素3维分布,利用高光谱激光雷达,采集了不同生长状态的绿萝叶片与植株的空间-光谱域点云数据,设计了一种基于分类预测的重建方法。通过偏最小二乘回归构建叶片正面与反面光谱的叶绿素含量预测模型,采用光谱自适应阈值选择方法实现植株点云中叶片正反面的分类,并根据类别标签选择模型计算叶绿素含量,重建植株的叶绿素3维分布。结果表明,该方法得到的植株叶绿素3维分布更接近真实值,决定系数达到0.69,均方根误差为4.97。这一结果可为植物表型研究提供新的数据基础和理论方法。
  • 图 1  HSL系统结构图

    Figure 1.  Schematic of HSL system

    图 2  绿萝实验样本

    Figure 2.  Experimental samples of epipremnum aureum

    图 3  样本扫描点

    Figure 3.  Scanning points of sample

    图 4  叶绿素含量3维重建框图

    Figure 4.  Schematic diagram of 3-D reconstruction of chlorophyll content

    图 5  植被指数与叶绿素含量相关矩阵图

    Figure 5.  Correlation matrix diagram of vegetation indices and chlorophyll content

    图 6  叶片正反面平均光谱与分布

    Figure 6.  Average spectrum and range of adaxial and abaxial leaves

    图 7  不同波长的差异比值F指数

    Figure 7.  F-value of different wavelength

    图 8  不同植株分类结果

    Figure 8.  Classification results of different plants

    图 9  两种重建方法对比

    a~c—不分类重建  d~f—分类重建

    Figure 9.  Comparison of two reconstruction methods

    a~c—reconstruction without classification  d~f—reconstruction with classification

    图 10  两种方法预测值与真实值关系

    a—不分类预测  b—分类预测

    Figure 10.  Relationship between the predicted and true values for two methods

    a—prediction without classification  b—prediction with classification

    表 1  各种最优植被指数的波长组合与相关系数

    Table 1.  Wavelength combinations and correlation coefficients of various optimal vegetation indices

    adaxial abaxial
    rmax position of the wavelength (λi, λj)/nm rmax position of the wavelength (λi, λj)/nm
    iRVI 0.83 (720, 840) 0.87 (725, 825)
    iDVI 0.83 (720, 905) 0.89 (720, 820)
    iNDVI 0.84 (705, 885) 0.87 (725, 825)
    下载: 导出CSV

    表 2  PLSR对叶片正反面建立的叶绿素含量预测模型

    Table 2.  Predictive model of chlorophyll content developed by PLSR on adaxial and abaxial leaves

    dataset model modeling dataset validation dataset
    R2e eRMSE R2 eRMSE
    adaxial y=12.19x1+1.48x2+30.91x3-8.09 0.75 4.31 0.80 4.49
    abaxial y=43.37x1+2.27x2+100.12x3-33.88 0.82 3.96 0.80 3.67
    下载: 导出CSV
  • [1]

    GITELSON A A, VIÑA A, VERMA S B, et al. Relationship between gross primary production and chlorophyll content in crops: Implications for the synoptic monitoring of vegetation productivity[J]. Journal of Geophysical Research: Atmospheres, 2006, 111: D08S11.
    [2] 杨敏文. 快速测定植物叶片叶绿素含量方法的探讨[J]. 光谱实验室, 2002, 19(4): 478-481. doi: 10.3969/j.issn.1004-8138.2002.04.015

    YANG M W. Study on rapid determination of chlorophyll content of leaves[J]. Chinese Journal of Spectroscopy Laboratory, 2002, 19(4): 478-481(in Chinese). doi: 10.3969/j.issn.1004-8138.2002.04.015
    [3]

    SIMS D A, GAMON J A. Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages[J]. Remote Sensing of Environment, 2002, 81(2-3): 337-354. doi: 10.1016/S0034-4257(02)00010-X
    [4]

    di VITTORIO A V. Enhancing a leaf radiative transfer model to estimate concentrations and in vivo specific absorption coefficients of total carotenoids and chlorophylls a and b from single-needle reflectance and transmittance[J]. Remote Sensing of Environment, 2009, 113(9): 1948-1966. doi: 10.1016/j.rse.2009.05.002
    [5]

    MATTHEWS M W, BERNARD S, ROBERTSON L. An algorithm for detecting trophic status (chlorophyll-a), cyanobacterial-dominance, surface scums and floating vegetation in inland and coastal waters[J]. Remote Sensing of Environment, 2012, 124: 637-652. doi: 10.1016/j.rse.2012.05.032
    [6] 刘爽, 于海业, 张郡赫, 等. 基于最优光谱指数的大豆叶片叶绿素含量反演模型研究[J]. 光谱学与光谱分析, 2021, 41(6): 1912-1919.

    LIU Sh, YU H Y, ZHANG J H, et al. Study on inversion model of chlorophyll content in soybean leaf based on optimal spectral indices[J]. Spectroscopy and Spectral Analysis, 2021, 41(6): 1912-1919(in Chinese). 
    [7] 王芳东, 严志雁, 赵小敏, 等. 油茶叶片叶绿素含量高光谱估测的偏最小二乘模型参数选择[J]. 江西农业大学学报, 2022, 44(1): 86-96.

    WANG F D, YAN Zh Y, ZHAO X M, et al. Parameter selection of partial least squares model for hyperspectral estimation of chlorophyll content in Camellia oleifera leaves[J]. Acta Agriculturae Universitatis Jiangxiensis, 2022, 44(1): 86-96(in Chinese). 
    [8] 张建, 李勇, 谢静, 等. 针对水稻植株三维叶绿素(SPAD) 信息获取的最优波段选择问题研究[J]. 光谱学与光谱分析, 2017, 37(12): 3749-3757.

    ZHANG J, LI Y, XIE J, et al. Research on optimal near-infrared band selection of chlorophyll (SPAD) 3-D distribution about rice plant[J]. Spectroscopy and Spectral Analysis, 2017, 37(12): 3749-3757(in Chinese). 
    [9] 王鑫, 王梓橦, 尤文强, 等. 利用叶片正反面反射光谱估算叶绿素含量[J]. 光谱学与光谱分析, 2018, 38(8): 2524-2528.

    WANG X, WANG Z T, YOU W Q, et al. Estimation of chlorophyll content by reflectance spectra of the positive and negative blades[J]. Spectroscopy and Spectral Analysis, 2018, 38(8): 2524-2528(in Chinese). 
    [10]

    EITEL J U H, VIERLING L A, LONG D S. Simultaneous measurements of plant structure and chlorophyll content in broadleaf saplings with a terrestrial laser scanner[J]. Remote Sensing of Environment, 2010, 114(10): 2229-2237. doi: 10.1016/j.rse.2010.04.025
    [11]

    ZHANG Zh J, WANG Ch, ZHANG H, et al. Analysis of permafrost region coherence variation in the Qinghai-Tibet Plateau with a high-resolution TerraSAR-X image[J]. Remote Sensing, 2018, 10(2): 298. doi: 10.3390/rs10020298
    [12]

    CALDERS K, ADAMS J, ARMSTON J, et al. Terrestrial laser scanning in forest ecology: Expanding the horizon[J]. Remote Sensing of Environment, 2020, 251: 112102. doi: 10.1016/j.rse.2020.112102
    [13]

    HOSOI F, OMASA K. Estimation of vertical plant area density profiles in a rice canopy at different growth stages by high-resolution portable scanning LiDAR with a lightweight mirror[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2012, 74: 11-19. doi: 10.1016/j.isprsjprs.2012.08.001
    [14]

    EITEL J U H, VIERLING L A, LONG D S. Simultaneous measurements of plant structure and chlorophyll content in broadleaf saplings with a terrestrial laser scanner[J]. Remote Sensing of Environment, 2010, 114(10): 2229-2237. doi: 10.1016/j.rse.2010.04.025
    [15]

    BI K, XIAO S, GAO S, et al. Estimating vertical chlorophyll concentrations in maize in different health states using hyperspectral LiDAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2020, 58(11): 8125-8133. doi: 10.1109/TGRS.2020.2987436
    [16]

    DU L, JIN Z, CHEN B, et al. Application of hyperspectral LiDAR on 3-D chlorophyll-nitrogen mapping of rohdea japonica in laboratory[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2021, 14: 9667-9679. doi: 10.1109/JSTARS.2021.3111295
    [17]

    SHAO H, CAO Z, LI W, et al. Feasibility study of wood-leaf separation based on hyperspectral LiDAR technology in indoor circumstances[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2021, 15: 729-738.
    [18]

    UDDLING J, GELANG-ALFREDSSON J, PⅡKKI K, et al. Evaluating the relationship between leaf chlorophyll concentration and SPAD-502 chlorophyll meter readings[J]. Photosynthesis Research, 2007, 91(1): 37-46. doi: 10.1007/s11120-006-9077-5
    [19] 金秀, 朱先志, 李绍稳, 等. 基于梯度提升树的土壤速效磷高光谱回归预测方法[J]. 激光与光电子学进展, 2019, 56(13): 131102.

    JIN X, ZHU X Zh, LI Sh W, et al. Predicting soil available phosphorus by hyperspectral regression method based on gradient boosting decision tree[J]. Laser & Optoelectronics Progress, 2019, 56 (13): 131102(in Chinese). 
    [20]

    GAULTON R, DANSON F M, RAMIREZ F A, et al. The potential of dual-wavelength laser scanning for estimating vegetation moisture content[J]. Remote Sensing of Environment, 2013, 132: 32-39.
    [21] 尼加提·卡斯木, 师庆东, 王敬哲, 等. 基于高光谱特征和偏最小二乘法的春小麦叶绿素含量估算[J]. 农业工程学报, 2017, 33(22): 208-216.

    NEJAT K, SHI Q D, WANG J Zh, et al. Estimation of spring wheat chlorophyll content based on hyperspectral features and PLSR model[J]. Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(22): 208-216(in Chinese). 
    [22] 龚威, 史硕, 陈必武, 等. 对地观测高光谱激光雷达发展及展望[J]. 遥感学报, 2021, 25(1): 501-513.

    GONG W, SHI Sh, CHEN B W, et al. Development and prospect of hyperspectral LiDAR for earth observation[J]. National Remote Sensing Bulletin, 2021, 25(1): 501-513(in Chinese). 
    [23] 张翠红, 张小娟, 朱大洲, 等. 基于激光散射图像小麦叶片叶绿素检测研究[J]. 激光技术, 2012, 36(4): 459-462.

    ZHANG C H, ZHANG X J, ZHU D Zh, et al. Detection of chlorophyll content of wheat leaves based on laser scattering images[J]. Laser Technology, 2012, 36(4): 459-462(in Chinese). 
  • [1] 刘璐邵慧孙龙陈杰徐恒胡玉霞肖晓 . 利用高光谱激光雷达检测木材的霉变与含水量. 激光技术, 2023, 47(5): 620-626. doi: 10.7510/jgjs.issn.1001-3806.2023.05.007
    [2] 薛芳谭勇邬志强王玉诏 . 基于激光遥感的酒驾遥感探测技术研究. 激光技术, 2019, 43(1): 93-97. doi: 10.7510/jgjs.issn.1001-3806.2019.01.019
    [3] 郝晓剑孙永凯 . 激光诱导击穿光谱用于煤中多元素同步检测. 激光技术, 2020, 44(1): 119-124. doi: 10.7510/jgjs.issn.1001-3806.2020.01.021
    [4] 宗蓉郭伟仲凌志 . 毫米波雷达和激光雷达探测云边界的对比分析. 激光技术, 2016, 40(2): 171-177. doi: 10.7510/jgjs.issn.1001-3806.2016.02.005
    [5] 何艳坤白玉杰 . 基于残差偏置和查找表的高光谱图像无损压缩. 激光技术, 2014, 38(5): 643-646. doi: 10.7510/jgjs.issn.1001-3806.2014.05.014
    [6] 刘翠连陶于祥罗小波李青妍 . 混合卷积神经网络的高光谱图像分类方法. 激光技术, 2022, 46(3): 355-361. doi: 10.7510/jgjs.issn.1001-3806.2022.03.009
    [7] 齐永锋马中玉 . 邻域谱概率协同表示的高光谱图像分类方法. 激光技术, 2019, 43(4): 448-452. doi: 10.7510/jgjs.issn.1001-3806.2019.04.003
    [8] 范琪朱晓林周鼎富周杰罗雄陈春利童文皓杨泽后 . 激光测风雷达分析典型高原机场风场特征. 激光技术, 2020, 44(5): 525-531. doi: 10.7510/jgjs.issn.1001-3806.2020.05.001
    [9] 孙越娇雷武虎胡以华赵楠翔任晓东 . 基于视觉显著模型的遥感图像舰船快速检测. 激光技术, 2018, 42(3): 379-384. doi: 10.7510/jgjs.issn.1001-3806.2018.03.017
    [10] 吴建军李磊方平凯孟小前谭均铭 . 电力巡线直升机激光扫描数据的高效组织与显示. 激光技术, 2019, 43(3): 318-323. doi: 10.7510/jgjs.issn.1001-3806.2019.03.006
    [11] 方俊龙郭宝峰沈宏海杨名宇 . 一种改进的基于自动形态学的端元提取算法. 激光技术, 2017, 41(1): 106-112. doi: 10.7510/jgjs.issn.1001-3806.2017.01.022
    [12] 王晓维梁应剑李翔龙彦志梅运桥 . 基于光学多普勒频移的低空速测量方法研究. 激光技术, 2016, 40(5): 629-632. doi: 10.7510/jgjs.issn.1001-3806.2016.05.003
    [13] 王琪杨桄张俭峰向英杰田张男 . 结合K-L散度和互信息的无监督波段选择算法. 激光技术, 2018, 42(3): 417-421. doi: 10.7510/jgjs.issn.1001-3806.2018.03.024
    [14] 张昌赛刘正军杨树文左志权 . 基于LiDAR数据的布料模拟滤波算法的适用性分析. 激光技术, 2018, 42(3): 410-416. doi: 10.7510/jgjs.issn.1001-3806.2018.03.023
    [15] 吴香伟郭宝峰陈春种沈宏海 . 基于加权组合核RX算法异物检测及其参量选择. 激光技术, 2015, 39(6): 745-750. doi: 10.7510/jgjs.issn.1001-3806.2015.06.003
    [16] 李峰吴燕雄卫爱霞刘文龙米晓楠 . 机载激光雷达3维建筑物模型重建的研究进展. 激光技术, 2015, 39(1): 23-27. doi: 10.7510/jgjs.issn.1001-3806.2015.01.005
    [17] 房垚鑫郭宝峰马超 . 基于改进点扩散函数的遥感图像超分辨率重建. 激光技术, 2019, 43(5): 713-718. doi: 10.7510/jgjs.issn.1001-3806.2019.05.024
    [18] 何宁谭智诚吕杏利蒋红艳 . 基于语音识别控制的激光超声水下遥感研究. 激光技术, 2017, 41(3): 391-395. doi: 10.7510/jgjs.issn.1001-3806.2017.03.017
    [19] 鲍鸿曾海涛白玉磊胡忠向志聪周延周申作春 . 基于概率密度最小二乘拟合的叶片后缘轮廓. 激光技术, 2016, 40(4): 555-559. doi: 10.7510/jgjs.issn.1001-3806.2016.04.021
    [20] 欧江霞刘伟诚 . 基于IGGⅢ方案的加权总体最小二乘点云球面拟合. 激光技术, 2017, 41(5): 749-753. doi: 10.7510/jgjs.issn.1001-3806.2017.05.026
  • 加载中
图(10) / 表(2)
计量
  • 文章访问数:  1502
  • HTML全文浏览量:  1028
  • PDF下载量:  13
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-07-25
  • 录用日期:  2022-09-08
  • 刊出日期:  2023-09-25

基于高光谱激光雷达的绿萝叶绿素3维重建

    通讯作者: 邵慧, shaohui@ahjzu.edu.cn
    作者简介: 汪慧民(1997-),男,硕士研究生,现主要从事激光雷达信号处理的研究
  • 1. 安徽建筑大学 电子与信息工程学院,合肥 230601
  • 2. 安徽建筑大学 安徽省古建筑智能感知与高维建模国际联合研究中心,合肥 230601
  • 3. 北京航空航天大学 无人机系统研究院,北京 100191
  • 4. 淮南师范学院 电子工程学院,淮南 232038
基金项目:  安徽省高校省级自然科学研究项目 KJ2021JD16安徽省高校协同创新项目 GXXT-2021-028安徽省高校省级自然科学研究项目 KJ2021A0622安徽省高校省级自然科学研究项目 KJ2020A0471安徽省科技厅面上项目 2008085MF182安徽省住房城乡建设科学技术计划资助项目 2020-YF22安徽省住房城乡建设科学技术计划资助项目 2022-YF077安徽省科技重大专项基金资助项目 202003a05020031

摘要: 为了准确重建正反面叶片叶绿素3维分布,利用高光谱激光雷达,采集了不同生长状态的绿萝叶片与植株的空间-光谱域点云数据,设计了一种基于分类预测的重建方法。通过偏最小二乘回归构建叶片正面与反面光谱的叶绿素含量预测模型,采用光谱自适应阈值选择方法实现植株点云中叶片正反面的分类,并根据类别标签选择模型计算叶绿素含量,重建植株的叶绿素3维分布。结果表明,该方法得到的植株叶绿素3维分布更接近真实值,决定系数达到0.69,均方根误差为4.97。这一结果可为植物表型研究提供新的数据基础和理论方法。

English Abstract

    • 叶绿素是植物进行光合作用合成养分的重要色素,准确测定植物叶绿素含量可以估计植物总初级生产力和监测植物生长状态[1]

      传统的叶绿素含量的测定方法主要有分光光度法、光声光谱法,其中分光光度法应用最为广泛,但是此类方法需现场破坏性采集样品,耗时耗力且受实验条件限制[2]。大量研究表明, 叶片的光谱特征与叶绿素含量之间密切相关[3-5],由于速度快、效率高、成本低、无损检测等优点,被动多光谱与高光谱方法被广泛用于分析、估计和反演叶片叶绿素含量。LIU等人使用地物光谱仪构建了大豆叶片叶绿素含量的预测模型,为大面积监测大豆的生长状态提供了参考[6]。WANG等人利用光谱仪分析了偏最小二乘回归(partial least squares regression, PLSR)中不同输入参数对油茶叶片叶绿素含量估算的影响[7]。ZHANG等人利用多光谱图像采集系统获取并重建了水稻叶绿素含量3维分布多视角图像[8]。WANG等人利用被动光谱仪在实验室环境中对植物叶片正反面光谱采集,通过PLSR构建了对叶片正反面光谱不敏感的叶绿素含量预测模型[9],由于光谱仪视场角问题,每次采集到的是整个视场范围内的平均光谱,完成整棵植株的单点精确扫描所需时间成本太高。被动遥感在外场实验中容易受到光源和背景条件的限制以及冠层结构的影响[10-11],因此主动探测遥感越来越受到学者们的关注[12]。激光雷达(light detection and ranging, LiDAR)作为一种主动遥感技术,不依赖于外界光源,可以高效获得植被3维结构信息[13]。EITEL等人使用地面激光扫描仪采集阔叶树苗532 nm的后向散射强度并反演了其叶绿素含量[14],但单波长激光雷达能够提供的光谱信息有限,在反演精度上受到限制。

      随着超连续谱激光器的发展,人们设计和发展了能够获得目标物体丰富光谱信息和精确空间信息的高光谱激光雷达(hyperspectral LiDAR, HSL),使得植被生化参数反演与3维重建成为可能。BI等人借助HSL光谱信息构建的PLSR模型,结合空间信息比较了不同施氮条件下玉米植株垂直方向上叶绿素浓度[15]。DU等人通过构建宽带植被指数,将HSL数据应用于万年青叶绿素-氮含量的3维重建[16]。以上研究考虑的都是叶片正面光谱与叶绿素含量的关系,在实际的遥感测量中,叶片的反面光谱不可避免地被接收到。仅依赖叶片正面光谱所建立的叶绿素含量预测模型不适用于叶片反面情况,造成预测精度的降低,影响到植物叶绿素3维分布模型准确建立。

      本文作者利用自研的百通道HSL,以绿萝为实验样本,利用最优植被指数结合偏最小二乘法构建叶片正面与反面光谱预测模型。分析叶片正反面的光谱差异,并提出自适应阈值选择方法实现绿萝植株正反面叶片的分类。在分类条件下,将预测模型应用于多基元植株点云,实现叶绿素含量3维分布重建。

    • 利用百通道地基HSL系统[17]采集样本空间-光谱域点云数据。系统由发射单元、扫描控制单元和信号采集单元3个部分组成,如图 1所示。通过超连续谱激光器和声光可调谐滤波器(acousto-optic tunable filters, AOTF)组成的发射单元可在不同时刻发射不同波长的激光脉冲,其光谱分辨率为5 nm,有效工作波长为550 nm~1050 nm。激光发射信号经光纤耦合到准直器,再经折转发射镜将不同波长的激光信号发射出去,然后由接收望远镜头将回波信号汇聚到探测器上,通过雪崩光电二极管(avalanche photodiode, APD)和高速采集卡记录不同时刻接收的不同波长的激光脉冲。同时利用2维转台记录目标的3维空间坐标,获得目标空间-光谱的点云信息。

      图  1  HSL系统结构图

      Figure 1.  Schematic of HSL system

    • 使用柯尼卡手持叶绿素仪(SPAD-502Plus)测量叶片叶绿素含量,该仪器可以测量叶片的相对叶绿素含量, 即土壤与作物分析开发(soil and plant analyzer development, SPAD)值,SPAD测量范围为0~200, 分辨率为0.1,准确度在±1以内,重复性在±3以内。由于SPAD值与叶绿素含量具有高度相关性,且能反映叶绿素含量的相对高低[18],本文中以其表征叶绿素含量。

    • 实验样本为绿萝叶片与植株,绿萝属于天南星科麒麟叶属植物,叶片硕大、叶茎呈绿色,如图 2所示。图 2a图 2b分别为绿萝叶片的正面与反面,将其固定在叶面夹上,共采集15枚绿萝叶片的数据;图 2c~图 2e为不同长势与不同生长阶段的绿萝植株,均为土栽法种植, 图 2c发育良好,叶片较大呈深绿色; 图 2d由于缺肥,底部叶片叶绿素合成受阻,开始出现泛黄现象; 图 2e处于生长初期,叶片较小。将不同样本置于红色泡沫板上,背景为黑布,其中叶片反面由红色方框框出。

      图  2  绿萝实验样本

      Figure 2.  Experimental samples of epipremnum aureum

    • 使用HSL采集每枚叶片的正面与反面各6个单点,用手持叶绿素仪测量相同位置的叶绿素含量,取3次测量的均值作为该位置的叶绿素含量。实验中共获取叶片正面和反面数据各90组,划分为正面数据集与反面数据集,每种数据集以2∶1比例划分为建模集和验证集。

      绿萝植株放置于距离HSL系统4.8 m处,HSL系统发射的激光光束从样本的左上角起始点开始从左至右扫描,当扫描至最右侧换行从右至左扫描,重复扫描过程直至完成整个样本扫描,其中水平扫描间距与垂直扫描间距均为0.05°。如图 3a所示,共获得9964个HSL扫描点,以固定距离值进行分割,消除背景回波信号的干扰,提取仅包含植株样本扫描点,再采用高斯滤波器去除随机噪声,共获得3400个点,如图 3b所示。

      图  3  样本扫描点

      Figure 3.  Scanning points of sample

      在相同条件下,采集反射率为99%标准漫反射参照板回波强度信号,利用下式,计算目标点的光谱反射率:

      $ \rho_{\mathrm{t}}\left(\lambda_i\right)=\frac{V_{\mathrm{t}}\left(\lambda_i\right)}{V_{\mathrm{b}}\left(\lambda_i\right)} \rho_{\mathrm{b}}\left(\lambda_i\right) $

      (1)

      式中, ρt(λi)为目标点在波长λi处的反射率,ρb(λi)为参照白板在波长λi处的反射率, Vt(λi)和Vb(λi)分别为目标点和参照板在波长λi处的回波信号电压峰值。

    • 图 4为绿萝叶绿素3维重建方法框图,包括预处理、分类、建模、叶绿素含量3维重建4个部分。

      图  4  叶绿素含量3维重建框图

      Figure 4.  Schematic diagram of 3-D reconstruction of chlorophyll content

      预处理部分包括3维空间坐标解算、反射率计算以及点云滤波,已在1.4节中讨论。建模部分通过单点测量的叶片反射率,提取最优植被指数,结合实测SPAD值,利用PLSR建立叶绿素预测模型。分类部分包括两个部分,首先计算各个波长差异比值指数选择特征波长,然后利用十折交叉验证确立最佳阈值,完成对植株点云的分类。依据植株样本每个扫描点的类别选择模型计算SPAD值,结合3维空间坐标完成叶绿素含量的3维重建。

    • 选择PLSR作为构建叶绿素含量预测模型的方法。PLSR是一种新型多元统计数据分析方法,能够有效克服自变量之间的相关性与冗余性,广泛应用于高光谱遥感定量反演中[19]。将叶片层次上构建的叶绿素反演模型应用到3维冠层结构上,消除角度与距离因素的影响是首要条件。相较于原始光谱,通过两个波段反射率数学组合的植被指数能够一定程度上消除这些影响[20]。因此,选择3种形式植被指数:比值植被指数(ratio vegetation index, RVI)iRVI、差值植被指数(difference vegetation index, DVI)iDVI和归一化差分植被指数(normalized difference vegetation index, NDVI)iNDVI作为PLSR的输入变量,3种植被指数的计算公式如下:

      $ i_{\mathrm{RVI}}=\frac{\rho\left(\lambda_j\right)}{\rho\left(\lambda_i\right)} $

      (2)

      $ i_{\mathrm{DVI}}=\rho\left(\lambda_j\right)-\rho\left(\lambda_i\right) $

      (3)

      $ i_{\mathrm{NDVI}}=\frac{\rho\left(\lambda_j\right)-\rho\left(\lambda_i\right)}{\rho\left(\lambda_j\right)+\rho\left(\lambda_i\right)} $

      (4)

      式中, ρ(λi)与ρ(λj)分别表示λiλj波长的反射率。

      为了挑选出效果最佳的波段组合,计算所有可能的波段组合与叶绿素含量间的相关系数,并以相关矩阵形式呈现,选择相关系数最高的波段组合作为最优波段组合。

    • 采用单阈值法实现植株扫描点中叶片正反面的分类。为了从高光谱丰富的特征中寻找最合适的特征波长与最佳阈值,以正反面叶片数据集作为训练样本,定义差异比值指数F(λi),计算公式如下:

      $ F\left(\lambda_i\right)=\frac{\left|\overline{\rho_{\text {ada }}}\left(\lambda_i\right)-\overline{\rho_{\text {aba }}}\left(\lambda_i\right)\right|}{\hat{\rho}_{\text {ada }}\left(\lambda_i\right) \cdot \hat{\rho}_{\text {aba }}\left(\lambda_i\right)} $

      (5)

      式中,ρada(λi)与$\hat{\rho}_{\text {ada }}\left(\lambda_i\right)$)分别表示正面数据集所有样本在λi波长的反射率平均值与标准差,ρaba(λi)与$\hat{\rho}_{\text {ada }}\left(\lambda_i\right)$)分别表示反面数据集所有样本在λi波长的反射率平均值与标准差。F(λi)越大, 表明在波长λi处叶片正反面反射率类间差异很大,且类内离散程度很低,越利于选取合适阈值进行分类,因此将最大F值对应波长作为特征波长λ

      特征波长确立后,采用十折交叉验证方法选取最佳阈值。将数据集等划分为10份,轮流将9组作为训练集,剩余1组作为测试集。在一次交叉验证中,依次将训练集每个样本在波长λ处的反射率作为阈值,进行分类并计算准确率,选择准确率最高一次对应阈值,完成对验证集的分类并计算准确率。最终,选择平均准确率最高一次交叉验证对应的阈值作为最佳阈值T

      λT确立后,遍历植株扫描点,若该点在λ波长处的反射率高于阈值T,则判定为叶片反面,反之为叶片正面。

    • 植株点云分类完成后,依据扫描点的类别标签选择对应模型完成SPAD值的计算。以HSL发射中心为原点建立空间直角坐标系,利用扫描点的3维坐标完成空间上的重建。将每个扫描点的SPAD值映射成不同颜色,将SPAD值融合到空间坐标上,实现叶绿素含量的3维重建。

    • 为了评价模型的好坏与植株整体预测的精度,采取决定系数R2与均方根误差(root mean square error, RMSE)eRMSE作为评价标准。

      $ R^2=1-\frac{\sum\limits_{i=1}^n\left(y_i-y_i{ }^{\prime}\right)^2}{\sum\limits_{i=1}^n\left(y_i-\bar{y}\right)^2} $

      (6)

      $ e_{\mathrm{RMSE}}=\sqrt{\frac{\sum\limits_{i=1}^n y_i-y_i{ }^{\prime}}{n}} $

      (7)

      式中, yi为第i个样本的叶绿素含量实测值,yi′为第i个样本的叶绿素含量预测值,y为叶绿素含量平均值,n为该数据集样本总数。

    • 图 5为不同植被指数与叶绿素含量的相关矩阵图。可以看出,相同的波段组合的不同植被指数与叶绿素含量间的相关系数r不同,选择出相关系数最大的波段组合才能确定最佳的植被指数。

      图  5  植被指数与叶绿素含量相关矩阵图

      Figure 5.  Correlation matrix diagram of vegetation indices and chlorophyll content

      提取与叶绿素相关系数最大的波段组合,如表 1所示。3种植被指数与叶绿素之间高度相关,相关系数均大于0.8。各种植被指数用到的波段组合虽然不同,但均利用了红边(670 nm~760 nm)内的波长,红边是绿色植物的特性,主要表现为反射率快速上升,叶绿素的含量高低直接影响红边反射率[21]

      表 1  各种最优植被指数的波长组合与相关系数

      Table 1.  Wavelength combinations and correlation coefficients of various optimal vegetation indices

      adaxial abaxial
      rmax position of the wavelength (λi, λj)/nm rmax position of the wavelength (λi, λj)/nm
      iRVI 0.83 (720, 840) 0.87 (725, 825)
      iDVI 0.83 (720, 905) 0.89 (720, 820)
      iNDVI 0.84 (705, 885) 0.87 (725, 825)
    • 将上述最优植被指数作为PLSR的输入变量,叶绿素含量作为因变量,建立回归模型,结果如表 2所示。表中,x1x2x3分别表示iRVIiDVIiNDVI,两种模型在验证集中均有较好的精度。

      表 2  PLSR对叶片正反面建立的叶绿素含量预测模型

      Table 2.  Predictive model of chlorophyll content developed by PLSR on adaxial and abaxial leaves

      dataset model modeling dataset validation dataset
      R2e eRMSE R2 eRMSE
      adaxial y=12.19x1+1.48x2+30.91x3-8.09 0.75 4.31 0.80 4.49
      abaxial y=43.37x1+2.27x2+100.12x3-33.88 0.82 3.96 0.80 3.67
    • 图 6为正反面数据集中叶片的平均光谱曲线与分布。可以看出,在750 nm~900 nm叶片的正反面平均反射率差异较大,且光谱波动较大,有部分重合,难以选取合适的阈值进行分类。而在700 nm~750 nm,均表现出叶片反面反射率大于叶片正面反射率。

      图  6  叶片正反面平均光谱与分布

      Figure 6.  Average spectrum and range of adaxial and abaxial leaves

      选择600 nm~1000 nm范围计算各个波长的F值,其大小如图 7所示。可以看出,在700 nm~750 nm的F值高于其它波长,这与观测结果是一致的,且Fmax=F(705),故选择705 nm作为分类的特征波长。其次,阈值设定为14%时,正反面数据集的分类准确率最高。最后,将植株扫描点中所有点在特征波长705 nm处反射率高于14%的判定为叶片反面,低于14%的判定为叶片正面,瓷盆未参与分类,结果如图 8所示。

      图  7  不同波长的差异比值F指数

      Figure 7.  F-value of different wavelength

      图  8  不同植株分类结果

      Figure 8.  Classification results of different plants

      将不同样本的分类结果与实物图对比,除叶片反面的边缘点误判成正面之外,总体分类准确较高。边缘点的误判主要原因是在激光扫描过程中,光斑没有完全打在叶片上,一部分被后方的黑布所吸收,导致回波信号强度值的降低,从而反射率相较于正常点较低,造成分类的错误[22]

    • 利用已构建的SPAD预测模型完成叶绿素含量的3维重建。考虑两种情况下的模型预测,对绿萝样本叶片不进行分类处理,即利用正面光谱模型反演样本的SPAD值。分类处理是在分类基础上,根据类别标签选择对应模型反演SPAD值。不同样本的HSL点云中包含叶片点、瓷盆点、边缘点和茎秆点。由于在构建模型的训练样本中,未考虑瓷盆等不具有叶片特性的点,也未考虑边缘点的影响,依据模型计算这些点存在极端的SPAD值,影响样本3维分布图的整体效果,将超出叶片正常SPAD值范围的异常值置为0。绿萝茎秆呈绿色,理论上存在一定含量叶绿素,但由于茎秆厚度较大,无法通过手持叶绿素仪获取稳定的读数[23],不能作为参考标准,因此不考虑茎秆点的预测效果,主要对叶片点进行分析。

      发育良好绿萝植株培育过程中受光较为均匀且处于健康状态,因此,各叶片间的SPAD值差异较小,集中在40~60之间。在不分类条件下进行预测的重建结果如图 9a~图 9c所示,框中各叶片的SPAD预测值区间集中在7~35之间,和实测值差距较大。而在分类条件下的重建结果如图 9d~图 9f所示,对应位置的叶片SPAD预测值集中在35~64之间,符合植株的真实情况。状态不良绿萝重建结果如图 9b图 9e所示,两种方法能对底部泛黄叶片有较好预测效果,但框中正常叶片被错误预测为低SPAD含量,而分类预测的重建结果能较为准确预测泛黄叶片与正常叶片的SPAD值。生长初期绿萝处于叶绿素合成阶段,各叶片SPAD在20~45之间。不分类重建如图 9c所示,框中叶片SPAD预测值接近于0,分类重建如图 9f所示,较为符合植株的真实叶绿素分布。

      图  9  两种重建方法对比

      Figure 9.  Comparison of two reconstruction methods

      为了定量分析两种方法植株整体的预测精度,在各个样本点云的叶片中心处共选择47个单点,利用手持叶绿素仪实测相同点作为叶绿素含量的真实值。两种方法预测值与真实值的关系如图 10a图 10b所示,图中散点越接近于1∶1的直线表示预测效果越好。可见,使用不分类的预测方法,叶片反面点的叶绿素含量预测值相较于真实值误差较大,表现为预测值远低于真实值。由于正面光谱所构建的叶绿素含量预测模型不适合反面叶片,因此造成样本整体预测精度过低,R2<0,eRMSE=14.45。使用分类的预测方法,样本整体预测精度提升到一个较好的程度,R2=0.69,eRMSE=4.97。因此,在HSL获得的样本点云中同时拥有正反面光谱时,基于先分类后预测的方法,可以重建更为精确的叶绿素3维分布。

      图  10  两种方法预测值与真实值关系

      Figure 10.  Relationship between the predicted and true values for two methods

    • 以绿萝为研究对象,测量了绿萝叶片正反面单点的高光谱数据与叶绿素含量,提取3种植被指数最优波长组合。将最优波长组合的植被指数作为偏最小二乘回归的输入变量,构建叶片正反面的叶绿素含量预测模型。提出了基于差异比值指数结合十折交叉验证的阈值选择策略,利用阈值法完成绿萝点云的正反面叶片分类,依据分类标签选择模型完成叶绿素含量的反演,将叶绿素含量融合到3维空间坐标上,实现植株叶绿素3维分布重建。研究结果表明,在HSL获取到包含植株正反面叶片点云时,基于分类预测的重建方法可以有效获取植株真实状态的叶绿素3维分布。

参考文献 (23)

目录

    /

    返回文章
    返回