高级搜索

留言板

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

高光谱成像的水稻冠层穗颈瘟早期识别

,  ,  ,  , 

上一篇

下一篇

袁建清, 仇逊超, 贾银江, 等. 高光谱成像的水稻冠层穗颈瘟早期识别[J]. 西南大学学报（自然科学版）, 2023, 45(6): 57-65. doi: 10.13718/j.cnki.xdzk.2023.06.006
引用本文: 袁建清, 仇逊超, 贾银江, 等. 高光谱成像的水稻冠层穗颈瘟早期识别[J]. 西南大学学报（自然科学版）, 2023, 45(6): 57-65. doi: 10.13718/j.cnki.xdzk.2023.06.006
shu
YUAN Jianqing, QIU Xunchao, JIA Yinjiang, et al. Early Identification of Rice Panicle Blast in Canopy Based on Hyperspectral Imaging[J]. Journal of Southwest University Natural Science Edition, 2023, 45(6): 57-65. doi: 10.13718/j.cnki.xdzk.2023.06.006
Citation: YUAN Jianqing, QIU Xunchao, JIA Yinjiang, et al. Early Identification of Rice Panicle Blast in Canopy Based on Hyperspectral Imaging[J]. Journal of Southwest University Natural Science Edition, 2023, 45(6): 57-65. doi: 10.13718/j.cnki.xdzk.2023.06.006
shu

高光谱成像的水稻冠层穗颈瘟早期识别

  • 1.

    哈尔滨金融学院 计算机系，哈尔滨 150030

  • 2.

    东北农业大学 电气与信息学院，哈尔滨 150030

  • 基金项目: 科技创新2030-“新一代人工智能”重大项目(2021ZD0110904)；黑龙江省省属本科高校基本科研业务费项目(2020-KYYWF-E009)；黑龙江省高等教育教学改革重点委托项目(SJGZ20200067)
详细信息
    作者简介:

    袁建清，博士，副教授，主要从事农业遥感与光谱分析研究 .

    通讯作者: 苏中滨，教授，博士研究生导师

  • 中图分类号: S435.111.1; S127

Early Identification of Rice Panicle Blast in Canopy Based on Hyperspectral Imaging

  • 1.

    Department of Computer Science, Harbin Finance University, Harbin 150030, China

  • 2.

    College of Electrical and Information, Northeast Agricultural University, Harbin 150030, China

  • 摘要: 为实现田间水稻冠层穗颈瘟的早期识别，利用室外高光谱成像系统采集早期自然发病大田的水稻冠层穗颈瘟图像，提取、分析反射率光谱特征.对预处理后的高光谱数据，采用主成分分析(Principal Component Analysis，PCA)、植被指数(Vegetation Index，VI)和竞争性自适应重加权法(Competitive Adaptive Reweighted Sampling，CARS)3种方法提取特征变量，结合支持向量机(Support Vector Machine，SVM)和线性判别分析(Linear Discriminant Analysis，LDA)分类算法构建识别模型.结果显示：以CARS特征波长和植被指数构建的模型，从分类结果看都取得了不错的效果，但是特征波长数量较多，可能存在过拟合的风险；单独使用PCA获得的主成分构建水稻冠层识别模型，没有明显效果.为此，研究尝试对选取的植被指数和提取的CARS特征使用PCA进一步降维，得到4个VI-PCs特征和5个CARS-PCs特征用于建模，取得了很好的效果.基于VI-PCs特征的SVM模型和LDA模型的总体分类精度分别为94%和95%；基于CARS-PCs特征的SVM模型和LDA模型总体分类精度分别为95%和97%，实现用较少变量获得较好的区分效果.从模型构建算法来看，LDA算法模型均优于SVM算法模型，说明LDA方法更适合于水稻冠层穗颈瘟识别模型的构建.研究可为航空、航天大面积的作物病虫害遥感监测提供理论依据.
    Abstract: In order to realize the early recognition of rice panicle blast in canopy in the field, the outdoor hyperspectral imaging system was used to collect the images of early stage natural occurring rice canopy panicle blast in the fields, then the spectral characteristics of reflectance were extracted and analyzed. For preprocessed hyperspectral data, PCA, vegetation index and CARS were used to extract feature variables, and build recognition models in combining with LDA and SVM classification algorithms. The results showed that the model constructed by CARS characteristic wavelength and vegetation index achieved good classification results. However, the number of characteristic wavelengths was large, there may be a risk of over-fitting. The principal components obtained by PCA alone were used to construct the rice canopy recognition model, which did not reflect the obvious effect. So, the study tried to use PCA to further reduce the dimension of the selected vegetation index and the extracted CARS features, and obtained 4 VI-PCs features and 5 CARS-PCs features for modeling, which achieved good results. The overall classification accuracy of SVM model and LDA model based on VI-PCs features was 94% and 95%, respectively. The overall classification accuracy of SVM model and LDA model based on CARS-PCs features was 95% and 97%, respectively, which achieved a better discrimination effect with fewer variables. From the view of model construction algorithm, the LDA algorithm model was superior to the SVM algorithm model, showing that LDA method is more suitable for the construction of rice canopy panicle blast identification model. The study can provide a theoretical basis for the aeronautical and astronautical remote sensing monitoring crop diseases and pests over a large area.
  • 加载中

  • 图 1  水稻冠层穗颈瘟样本RIO选取示意图

    下载: 全尺寸图片 幻灯片

    图 2  水稻冠层全样本光谱

    下载: 全尺寸图片 幻灯片

    图 3  水稻冠层穗颈瘟的平均光谱曲线

    下载: 全尺寸图片 幻灯片

    图 4  CARS特征波长选择过程

    下载: 全尺寸图片 幻灯片

    图 5  LDA模型的建模分类结果图

    下载: 全尺寸图片 幻灯片

    表 1  水稻冠层穗颈瘟前8个主成分的累积贡献率

    主成分 累积贡献率% 主成分 累积贡献率%
    PC1 70.12 PC5 92.73
    PC2 80.53 PC6 95.65
    PC3 83.48 PC7 96.97
    PC4 89.22 PC8 97.78
    下载: 导出CSV

    表 2  冠层穗颈瘟植被指数

    名称 定义 描述
    NRI (R570-R670)/(R570+R670) 氮反射指数
    PRI (R531-R570)/(R531+R570) 光化学植被指数
    GNDVI (R786-R566)/(R786+R566) 绿光归一化植被指数
    DVI R786-R682 差值植被指数
    SGI $\operatorname{mean}\left(\sum\limits_{500}^{600} R_i\right)$ 绿波段总和指数
    SAVI $1.5\left(\frac{R_{786}-R_{682}}{R_{786}+R_{682}+0.5}\right)$ 土壤调整植被指数
    VARI $\frac{\left(R_{566}-R_{682}\right)}{\left(R_{566}+R_{682}-R_{450}\right)}$ 可见光大气阻抗指数
    Db R′530 蓝边(490~530 nm)最大一阶导数值
    SDb SR′490~530 蓝边一阶导数总和
    Dy R′550 黄边(550~582 nm)最大一阶导数值
    SDy SR′550~582 黄边一阶导数总和
    Rr R′670 红边(640~680 nm)最大一阶导数值
    SRr SR′640~680 红边一阶导数总和
    下载: 导出CSV

    表 3  不同特征变量建模结果

    建模算法 特征变量 参数个数 建模精度/% 预测精度/% 总体精度/% 建模算法 特征变量 参数个数 建模精度/% 预测精度/% 总体精度/%
    SVM 全谱 173 90.91 88.23 90 LDA 全谱 173 93.39 88.23 92
    PCs 8 92.42 88.23 91 PCs 8 93.39 88.23 92
    植被指数 13 95.45 91.18 94 植被指数 13 98.48 91.18 96
    VI-PCs参数 4 96.97 88.23 94 VI-PCs参数 4 96.97 91.18 95
    CARS特征 14 96.97 94.12 96 CARS特征 14 98.48 94.12 97
    CARS-PCs参数 5 96.97 91.18 95 CARS-PCs参数 5 98.48 94.12 97
    下载: 导出CSV
  • [1] 杨燕, 何勇. 基于高光谱图像的稻瘟病抗氧化酶值早期预测[J]. 农业工程学报, 2013, 29(20): 135-141. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201320019.htm
    [2] 吴习宇, 曾凯芳, 郭启高, 等. 基于高光谱成像的枇杷果实品质检测[J]. 西南大学学报(自然科学版), 2022, 44(9): 52-60. doi: http://xbgjxt.swu.edu.cn/article/doi/10.13718/j.cnki.xdzk.2022.09.006
    [3] XIE C Q, CHU B Q, HE Y. Prediction of Banana Color and Firmness Using a Novel Wavelengths Selection Method of Hyperspectral Imaging[J]. Food Chemistry, 2018, 245(4): 132-140.
    [4] 孙俊, 靳海涛, 芦兵, 等. 基于高光谱图像及深度特征的大米蛋白质含量预测模型[J]. 农业工程学报, 2019, 35(15): 295-303. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201915036.htm
    [5] 郭琳, 王化宏, 程道德, 等. 施氮量对覆膜旱作水稻产量形成和冠层光谱特性的影响[J]. 西南大学学报(自然科学版), 2021, 43(1): 67-74. doi: http://xbgjxt.swu.edu.cn/article/doi/10.13718/j.cnki.xdzk.2021.01.008
    [6] WANG J Q, WU W M, WANG T W, et al. Estimation of Leaf Chlorophyll Content and Density in Populus Euphratica Based on Hyperspectral Characteristic Variables[J]. Spectroscopy Letters, 2018, 51(9): 485-495. doi: 10.1080/00387010.2018.1523194
    [7] 秦占飞, 常庆瑞, 谢宝妮, 等. 基于无人机高光谱影像的引黄灌区水稻叶片全氮含量估测[J]. 农业工程学报, 2016, 32(23): 77-85. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201623011.htm
    [8] WANG Y, JIANG F, GUPTA B B, et al. Variable Selection and Optimization in Rapid Detection of Soybean Straw Biomass Based on CARS[J]. IEEE Access, 2017, 6(10): 5290-5299.
    [9] LU J Z, ZHOU M C, GAO Y W, et al. Using Hyperspectral Imaging to Discriminate Yellow Leaf Curl Disease in Tomato Leaves[J]. Precision Agriculture, 2018, 19(3): 379-394.
    [10] JIN X, JIE L, WANG S A, et al. Classifying Wheat Hyperspectral Pixels of Healthy Heads and Fusarium Head Blight Disease Using a Deep Neural Network in the Wild Field[J]. Remote Sensing, 2018, 10(3): 395.
    [11] ASHOURLOO D, AGHIGHI H, MATKAN A A, et al. An Investigation into Machine Learning Regression Techniques for the Leaf Rust Disease Detection Using Hyperspectral Measurement[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(9): 4344-4351.
    [12] 秦立峰, 张熹, 张晓茜. 基于高光谱病害特征提取的温室黄瓜霜霉病早期检测[J]. 农业机械学报, 2020, 51(11): 212-220. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-NYJX202011023.htm
    [13] 黄双萍, 孙超, 齐龙, 等. 基于深度卷积神经网络的水稻穗瘟病检测方法[J]. 农业工程学报, 2017, 33(20): 169-176. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201720021.htm
    [14] 袁建清, 苏中滨, 贾银江, 等. 基于高光谱成像的寒地水稻叶瘟病与缺氮识别[J]. 农业工程学报, 2016, 32(13): 155-160. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201613022.htm
    [15] ZHANG G S, XU T Y, TIAN Y W, et al. Assessment of Rice Leaf Blast Severity Using Hyperspectral Imaging during Late Vegetative Growth[J]. Australasian Plant Pathology, 2020, 49(5): 571-578.
    [16] 朱梦远, 杨红兵, 李志伟. 高光谱图像和叶绿素含量的水稻纹枯病早期检测识别[J]. 光谱学与光谱分析, 2019, 39(6): 1898-1904. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201906045.htm
    [17] 康丽, 袁建清, 高睿, 等. 高光谱成像的水稻稻瘟病早期分级检测[J]. 光谱学与光谱分析, 2021, 41(3): 898-902. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN202103042.htm
    [18] KOBAYASHI T, SASAHARA M, KANDA E, et al. Assessment of Rice Panicle Blast Disease Using Airborne Hyperspectral Imagery[J]. The Open Agriculture Journal, 2016, 10(1): 28-34.
    [19] 孔繁昌, 刘焕军, 于滋洋, 等. 高寒地区粳稻穗颈瘟的无人机高光谱遥感识别[J]. 农业工程学报, 2020, 36(22): 68-75. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU202022009.htm
    [20] 张竞成. 多源遥感数据小麦病害信息提取方法研究[D]. 杭州: 浙江大学, 2012.
    [21] GITELSON A, MERZLYAK M N. Spectral Reflectance Changes Associated with Autumn Senescence of Aesculus Hippocastanum L. and Acer Platanoides L. Leaves. Spectral Features and Relation to Chlorophyll Estimation[J]. Journal of Plant Physiology, 1994, 143(3): 286-292.
    [22] 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.
    [23] GITELSON A A, ZUR Y, CHIVKUNOVA O B, et al. Assessing Carotenoid Content in Plant Leaves with Reflectance Spectroscopy[J]. Photochemistry and Photobiology, 2002, 75(3): 272-281.
    [24] HUETE A R, LIU H Q, BATCHILY K, et al. A Comparison of Vegetation Indices over a Global Set of TM Images for EOS-MODIS[J]. Remote Sensing of Environment, 1997, 59(3): 440-451.
    [25] LOBELL D B, ASNER G P. Hyperion Studies of Crop Stress in Mexico[C]. Pasadena: Proceedings ofthe 12th Annual JPL Airborne Earth Science Workshop, 2004: 1-6.
  • 加载中
图( 5) 表( 3)
计量
  • 文章访问数:  3316
  • HTML全文浏览数:  3316
  • PDF下载数:  666
  • 施引文献:  0
出版历程
  • 收稿日期:  2023-03-09
  • 刊出日期:  2023-06-20

高光谱成像的水稻冠层穗颈瘟早期识别

    通讯作者: 苏中滨，教授，博士研究生导师
    作者简介: 袁建清，博士，副教授，主要从事农业遥感与光谱分析研究
  • 1. 哈尔滨金融学院 计算机系，哈尔滨 150030
  • 2. 东北农业大学 电气与信息学院，哈尔滨 150030
  • 收稿日期:  2023-03-09
基金项目:  科技创新2030-“新一代人工智能”重大项目(2021ZD0110904)；黑龙江省省属本科高校基本科研业务费项目(2020-KYYWF-E009)；黑龙江省高等教育教学改革重点委托项目(SJGZ20200067)

摘要: 为实现田间水稻冠层穗颈瘟的早期识别，利用室外高光谱成像系统采集早期自然发病大田的水稻冠层穗颈瘟图像，提取、分析反射率光谱特征.对预处理后的高光谱数据，采用主成分分析(Principal Component Analysis，PCA)、植被指数(Vegetation Index，VI)和竞争性自适应重加权法(Competitive Adaptive Reweighted Sampling，CARS)3种方法提取特征变量，结合支持向量机(Support Vector Machine，SVM)和线性判别分析(Linear Discriminant Analysis，LDA)分类算法构建识别模型.结果显示：以CARS特征波长和植被指数构建的模型，从分类结果看都取得了不错的效果，但是特征波长数量较多，可能存在过拟合的风险；单独使用PCA获得的主成分构建水稻冠层识别模型，没有明显效果.为此，研究尝试对选取的植被指数和提取的CARS特征使用PCA进一步降维，得到4个VI-PCs特征和5个CARS-PCs特征用于建模，取得了很好的效果.基于VI-PCs特征的SVM模型和LDA模型的总体分类精度分别为94%和95%；基于CARS-PCs特征的SVM模型和LDA模型总体分类精度分别为95%和97%，实现用较少变量获得较好的区分效果.从模型构建算法来看，LDA算法模型均优于SVM算法模型，说明LDA方法更适合于水稻冠层穗颈瘟识别模型的构建.研究可为航空、航天大面积的作物病虫害遥感监测提供理论依据.

English Abstract

    全文HTML

  • 开放科学（资源服务）标志码（OSID）：

  • 稻瘟病传染速度快，危害严重，是水稻重要病害之一[1]. 水稻稻瘟病最为普遍的是水稻拔节期的叶瘟病和灌浆期的穗颈瘟. 穗颈瘟病斑位置主要集中在穗颈、穗轴或穗粒上，直接影响水稻的质量和产量. 在水稻安全生产中，加强穗颈瘟病害的防治是很重要的环节. 目前实际生产中的主要防治措施是政府组织集中利用无人机进行统一喷药或由农户凭经验提前喷药预防，农药的过度施用，使得农产品药残超标，从而破坏生态，污染环境. 因此，水稻田间发病情况的准确监测和早期识别，对水稻稻瘟病的提前防治、精准施药具有实际指导意义.

    高光谱成像技术具有无损、高效的特点，被广泛应用于作物品质检测[2-4]、作物养分及产量监测[5-8]和病虫害识别[9-12]等. 在水稻病虫害方面，黄双萍等[13]使用一种改进的深度卷积神经网络算法，提取水稻穗株的穗颈瘟病斑，建立了水稻穗颈瘟的检测模型；袁建清等[14]采用不同光谱预处理方法，结合不同特征提取方法和分类算法构建模型，实现水稻缺氮叶片和不同感病程度水稻叶瘟病叶片的判别；Zhang等[15]构建了基于叶片光谱比值的支持向量机评价模型，实现了水稻叶片尺度的叶瘟病分级检测；朱梦远等[16]采用多种预处理方法和不同分类算法建立模型，实现水稻纹枯病的早期检测；康丽等[17]提取特征变量作为输入，建立了支持向量机多级分类模型，实现水稻叶瘟病早期检测. Kobayashi等[18]利用航空高光谱遥感图像，根据特征波段对水稻稻瘟病田间感病区域进行评估；孔繁昌等[19]利用无人机高光谱平台获取水稻穗颈瘟冠层数据，使用随机森林方法构建了水稻穗颈瘟分类模型. 综上所述，研究者们对于实验室条件下的水稻病害机理做了深入的研究，对田间稻瘟病的外业监测和识别还处于探索阶段，需要进行更深入的研究.

    近地冠层水稻稻瘟病的监测同航空、航天遥感环境接近，可以衔接机理研究和实际应用研究. 本研究在之前初探性研究的基础之上，进行新数据采集，扩充数据量，对水稻冠层穗颈瘟的早期识别进行深入研究. 研究利用高光谱成像平台获取自然感病水稻大田的近地冠层图像，提取并分析其光谱反射率特征，采用不同特征提取方法，分别结合判别分析和支持向量机分类算法，构建水稻冠层穗颈瘟早期识别模型，以期为稻瘟病的无人机航空、卫星航天大面积遥感监测和识别提供一定的理论和技术支持.

1.   材料与方法

    1.1.   高光谱图像采集

  • 近地冠层水稻穗颈瘟数据的采集在黑龙江省哈尔滨市方正县水稻研究院和稻作博物馆进行. 方正县位于黑龙江省中南部，松花江中游南岸、长白山支脉张广才岭北段西北麓，地理位置为128°13′41″-129°33′20″E，45°32′46″-46°09′00″N，面积约为3 000 km2. 方正县以水稻生产为主，水稻为一年一熟，一般在每年的5月中旬进行移栽插秧，10月初进行收割. 本研究的信息采集于2019年8月16日和8月17日的11：00-14：00进行，天气晴朗少风. 水稻穗颈瘟为自然感病，处于感病初期. 共采集5块自然感病地块图像，水稻品种均为“黑绥1166”.

    研究使用Head Wall高光谱成像系统的室外平台进行数据采集，室外采集系统由高光谱相机、图像采集仪、推扫平台、笔记本电脑、三脚架和白板构成，另有发电机、稳压器、线轴和电源等设备. 每次采集数据时，首先在推扫平台上把高光谱成像仪固定好，进行高光谱成像仪的推扫设置，根据实景情况设置倾斜角度和左右摆扫角，角速度设置为0.5 deg/s，接着进行相机镜头的对焦调试，再进行暗电流校正和亮电流校正，用镜头盖遮住镜头进行暗电流校正，获得暗场高光谱图像，用白板进行亮电流校正，获得标准光谱图像，对待校正光谱图像经公式(1)计算后，最终得到校正后的高光谱图像，即为真实的水稻冠层穗颈瘟高光谱图像.

    式中：R为校正后高光谱图像的光谱反射率；Rr为待校正高光谱图像的光谱反射率；Rd为暗电流光谱反射率；Rw为亮电流光谱反射率.

    为了消除拍摄时可能产生的几何变形和失真现象，消除光照等产生的图像误差和辐射噪声，使用RPC正射校正法和暗像素法对水稻穗颈瘟冠层影像进行几何校正和辐射校正. 对校正之后的影像，使用ENVI软件的RIO工具进行感兴趣区的选取，如图 1所示，每个感兴趣区大小为40×25 Pix，这1 000个像素的平均反射率光谱即为一个冠层光谱样本. 根据5块水稻地块穗颈瘟感病情况，结合图像采集时的标记，共提取得到中度感病样本36个，轻度感病和无病样本各32个(图像采集为穗颈瘟感病早期，经植保专家鉴定，没有重度感病情况)，共计100个光谱样本用于研究.

    • 1.2.   高光谱数据预处理与裁剪

  • 本研究采用5幅图像，由于图像采集时间和地点不同可能会引起偏差，因此对提取的水稻冠层穗颈瘟反射率光谱通过公式(2)进行标准化处理，得到标准化光谱[20].

    式中：Ref′i为标准化光谱的光谱反射率；Refi为原光谱的光谱反射率；n为波段数.

    然后对标准化后的光谱进行5点、3项式Savitzky-Golay卷积平滑预处理，以消除样本间散射所引起的误差，降低仪器和环境等产生的噪声影响.

    Head Wall高光谱成像系统采集的光谱范围为400~1 000 nm，采样间隔为3 nm. 910~1 000 nm波段范围光谱波动剧烈，说明受环境、机器噪音等的影响较大，失真严重，因此对这部分数据进行裁剪，保留400~910 nm波段的光谱用于后续分析.

    所有样本平滑后的反射率光谱如图 2所示.

    • 1.3.   特征变量提取

  • 高光谱数据中的大量冗余、多重相关性的信息，会增加模型复杂度，降低模型精度和运算速度. 本研究采用主成分分析、竞争性自适应重加权算法和植被指数方法提取水稻冠层穗颈瘟的光谱特征参数.

    主成分分析法(Principal Component Analysis，PCA)是通过线性变换的方式将多变量转换为信息不重叠、相互正交的新变量，即为主成分. PCA能够尽可能多地保留有效信息，并解决了大量原始数据信息的重叠性问题，在高维数据降维、特征提取等方面得到广泛应用.

    竞争性自适应重加权法(Competitive Adaptive Reweighted Sampling，CARS)是基于蒙特卡洛采样法和偏最小二乘PLS模型回归系数的特征变量选择方法. 每次通过自适应加权采样技术保留PLS模型中回归系数绝对值权重较大的变量，去掉权重值较小的变量，建立基于新子集的PLS模型，反复多次运算后，选择交互验证均方根误差最小子集中的变量作为特征波长.

    植被指数法(Vegetation Index，VI)是由两个或多个光谱数据，经线性和非线性变换得到的具有一定植被指示意义的各种数值. 植被指数被广泛地应用于外业监测中，研究参考文献，选取能够反映植被胁迫和病害识别的植被指数、适合于冠层遥感观测的植被指数以及三边一阶导数参数，用于水稻穗颈瘟的判别.

    • 1.4.   分类建模算法

  • 研究主要采用支持向量机(Support Vector Machine，SVM)和线性判别分析(Linear Discriminant Analysis，LDA)两种分类算法进行水稻冠层穗颈瘟早期识别模型的构建.

    支持向量机是一种具有监督机制的机器学习方法，它通过结构风险最小化原理、最优化理论和核方法来提高泛化能力，适用于小样本、非线性以及高维特征空间等复杂问题的处理，在很多领域得到广泛应用. 线性判别分析的主要思想是将一个高维空间中的数据根据特征投影到一个低维的空间向量上，使得投影后各类别的类内方差小而类间均值差别大，即相同类别的样本投影点聚在一起，异类样本的投影点相距较远. LDA是一种经典的线性学习方法，广泛地被应用于模式分类等领域.

    本研究按照2∶1的比例随机划分建模集和预测集，随机取出中度感病样本24个、轻度感病和无病各21个，共计66个样本用于模型的建立，剩余34个样本用于模型的验证.

    光谱数据提取与处理、数据预处理、光谱特征提取和建模用到的软件有ENVI 5.1，Excel 2010，SPSS 21，The Unscrambler 10.4和Matlab 7.1.

2.   结果与分析

    2.1.   近地冠层水稻穗颈瘟光谱特征分析

  • 分别计算水稻冠层中度感病、轻度感病和无病样本的光谱平均值，得到3类样本的平均高光谱反射率光谱曲线. 如图 3所示，感病水稻冠层反射率光谱整体变得平缓，感病越严重越平缓. 比较3类样本存在着3个明显差异区间：首先是绿波段反射蜂处，这里是以566 nm为中心的波段区域，感病水稻低于健康水稻，感病越严重反射率越低. 其次是红光吸收谷处，感病水稻在这个区间反射率高于健康水稻，且中度感病高于轻度感病，说明感病越严重在红光区域吸收作用越弱. 最后是近红外高反射平台，感病水稻反射率低于健康水稻，感病越严重反射率越低. 分析出现这些特征的原因，主要是感病水稻因稻瘟病病菌的侵入破坏了水稻叶片、穗颈、穗轴等的细胞结构，引起其间隙发生变化，使得叶绿素等生化含量、持水量降低，光合作用强度降低，导致感病水稻在绿光区间和近红外区间反射能力降低，在红光区域吸收能力减弱. 通过上述分析得到水稻冠层穗颈瘟反射率光谱的3个敏感波段区域，分别是以566 nm为中心的绿波段区域、以786 nm为中心的近红外波段区域及以682 nm为中心的红光吸收谷处的波段区域.

    • 2.2.   PCA降维与植被指数选取

  • 用PCA方法对水稻冠层光谱进行降维处理，提取主成分，前8个主成分(PCs)累积贡献率为97.78%，各个主成分贡献率如表 1所示.

    研究参阅文献[21-25]，选取能够反映植被受病害胁迫、表征生物量变化的植被指数，包括氮反射指数NRI、光化学植被指数PRI等5个植被指数，以及适合冠层观测的土壤调整植被指数SAVI和可见光大气阻抗指数VARI，加上植被病害分析中常用的蓝边、黄边和红边3边参数共计13个参数. 通过单因素方差分析检验这13个参数对水稻冠层穗颈瘟均具有显著区分能力(显著水平为sig=0.05). 植被指数的定义和描述见表 2. 根据2.1中水稻冠层穗颈瘟的光谱反射率特征分析得到的敏感波段，计算宽波段植被指数时各参数取值为：蓝光为450 nm，绿光为566 nm，红光为682 nm，近红外光线为786 nm.

    13个植被指数数量比较多，参数间有可能存在数据共线性，因此本研究尝试使用PCA法对13个植被参数做进一步降维处理，PCA处理显示前4个主成分的累积贡献率达到99.13%，将这4个植被指数主成分(记为VI-PCs)用于模型构建分析.

    • 2.3.   CARS特征选取

  • 用CARS对水稻冠层全部样本进行降维，蒙特卡洛采样次数设为50，特征波长的筛选过程如图 4所示. 图 4a显示，随着采样次数的增加，特征波长的个数逐步减小，图 4b显示，PLS模型的RMSECV随着采样次数的增加先是缓慢下降后又增大，图 4c表示变量回归系数的变化趋势，“*”是RMSECV的最小位置，即采样次数为24次时，RMSECV值最小，得到最优的特征变量组合. 共筛选出14个CARS特征波长，分别是412.77，534.27，590.57，620.21，688.37，709.12，729.86，753.57，765.42，800.98，812.84，827.65，851.36，860.25 nm.

    CARS筛选出14个特征波长，与训练样本数量相比较，特征数量偏多，可能存在过拟合的风险，本研究用PCA法对14个CARS特征波长进行降维处理，PCA处理得到前5个主成分的累积贡献率达到99.29%，将这5个CARS主成分(记为CARS-PCs)也用于后续模型构建分析.

    • 2.4.   模型评估与分析

  • 分别将全谱、PCA提取的8个主成分因子、13个植被指数、14个CARS特征波长、4个VI-PCs特征和5个CARS-PCs特征作为输入，采用LDA线性判别分析和SVM支持向量机算法构建水稻冠层穗颈瘟早期识别模型. LDA算法选用Fisher线性判别；SVM算法采用RBF核函数，利用GRID网格搜索法确定分类机参数惩罚系数C和gamma的最优组合. 建模结果如表 3所示.

    从表 3可以看出，特征变量选取方法中PCA主成分分析的降维方法在水稻冠层穗颈瘟建模中没有体现出明显的效果，与全谱建模相比基本没有多大改善，分析原因可能PCA方法不能够很好地剔除室外光线和环境等噪声的影响. 作为外业遥感特别是航空航天遥感常用的植被指数法在近地冠层水稻穗颈瘟的识别中取得了较好的效果，基于植被指数的SVM模型和LDA模型总体分类精度分别达到了94%和96%，说明将不同波段进行组合、变换确实可以从不同角度反映植被受病害胁迫而引起的生化指标的变化，从而反演植被受病虫害胁迫的情况. 用CARS方法进行特征波长的筛选取得效果也很好，以CARS特征波长作为输入的SVM模型和LDA模型总体分类精度分别达到了96%和97%.

    然而，与训练样本数量相比较，选取的植被指数特征和筛选的CARS特征数量偏多，可能会存在过拟合的风险，因此本研究将选取的植被指数和CARS方法筛选出的特征波长再结合PCA方法进一步降维，得到4个VI-PCs特征和5个CARS-PCs特征用于建模：基于VI-PCs特征的SVM模型和LDA模型的总体分类精度分别达到94%和95%；基于CARS-PCs特征的SVM模型和LDA模型总体分类精度分别达到95%和97%. 实现了用较少变量获得较好区分的效果，同时避免了过拟合的风险.

    从模型构建算法上看，LDA线性判别的方法更适合于水稻冠层穗颈瘟分类模型的构建，无论以哪种特征变量作为输入，模型的分类精度都优于SVM模型. 图 5是基于VI-PCs参数和CARS-PCs参数的LDA模型建模分类结果图. 从图中可以清晰看到，基于VI-PCs参数和CARS-PCs参数的LDA模型在建模时分别出现两个中度感病被误判为轻度感病和一个中度感病被误判为轻度感病的情况.

3.   结论

  • 本研究以大田自然发病为基础，用室外高光谱成像系统采集早期发病大田的近地冠层水稻穗颈瘟图像，提取、分析反射率光谱. 对光谱进行标准化处理和Savitzky-Golay卷积平滑预处理，分析水稻冠层的光谱反射率特征，采用不同数据降维、特征提取方法和LDA，SVM分类算法，构建水稻冠层穗颈瘟早期识别模型. 主要结论有：

    1) 研究水稻灌浆期中度感病、轻度感病和无病3类样本的冠层反射率光谱，发现感病水稻的光谱反射率较无病水稻整体平缓，感病越严重水稻在红光区域吸收能力越低，在绿光区域和近红外区域反射能力越低. 分析得出水稻冠层穗颈瘟反射率光谱的3个敏感波段区域，即以566 nm为中心的绿波段区域、以786 nm为中心的近红外波段区域、以682 nm为中心的红光吸收谷区域.

    2) 本研究分别采用PCA、植被指数和CARS 3种方法提取光谱特征，构建识别模型. 单独使用PCA降维方法，获得主成分因子构建水稻冠层识别模型，没有取得明显效果. 选用植被指数和CARS方法筛选特征波长进行模型构建，模型都获得较高的区分精度. 为了避免所用特征波长数量较多、训练样本数量不多而可能存在过拟合的风险，本研究尝试对筛取的植被指数和CARS特征波长，用PCA方法进一步降维，获得4个VI-PCs参数和5个CARS-PCs参数用于建模，模型都获得较好的区分效果：基于VI-PCs参数的SVM模型和LDA模型总体分类精度分别为94%和95%；基于CARS-PCs参数的SVM模型和LDA模型总体分类精度分别为95%和97%.

    3) 从模型构建算法看，LDA方法更适合于水稻冠层穗颈瘟识别模型的构建，基于各种特征模型的分类精度均优于SVM算法模型.

    本研究实现了田间水稻冠层穗颈瘟病害的早期识别，可以为无人机航空遥感和卫星航天遥感的大面积病虫害监测与识别提供一定的理论和技术依据.

参考文献 (25)

目录

版权所有：西南大学学报（自然科学版）

地址：重庆市北碚区西南大学 电话：023-68254576, 68252540, 68366165,68366917 电子邮箱：lkxb@swu.edu.cn

本系统由北京仁和汇智信息技术有限公司开发

/

返回文章
返回