CN116416524A - Early-stage asymptomatic detection method for bacterial leaf blight of rice - Google Patents
Early-stage asymptomatic detection method for bacterial leaf blight of rice Download PDFInfo
- Publication number
- CN116416524A CN116416524A CN202310330513.8A CN202310330513A CN116416524A CN 116416524 A CN116416524 A CN 116416524A CN 202310330513 A CN202310330513 A CN 202310330513A CN 116416524 A CN116416524 A CN 116416524A
- Authority
- CN
- China
- Prior art keywords
- hyperspectral
- characteristic
- 3dcnn
- rice
- asymptomatic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/188—Vegetation
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/0464—Convolutional networks [CNN, ConvNet]
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/764—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/7715—Feature extraction, e.g. by transforming the feature space, e.g. multi-dimensional scaling [MDS]; Mappings, e.g. subspace methods
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/774—Generating sets of training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/80—Fusion, i.e. combining data from various sources at the sensor level, preprocessing level, feature extraction level or classification level
- G06V10/806—Fusion, i.e. combining data from various sources at the sensor level, preprocessing level, feature extraction level or classification level of extracted features
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/82—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using neural networks
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/194—Terrestrial scenes using hyperspectral data, i.e. more or other wavelengths than RGB
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Multimedia (AREA)
- Software Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Computing Systems (AREA)
- Databases & Information Systems (AREA)
- Medical Informatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Molecular Biology (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Remote Sensing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
Description
技术领域technical field
本发明涉及水稻白叶枯病检测技术领域,具体涉及一种水稻白叶枯病早期无症状检测方法。The invention relates to the technical field of rice bacterial blight detection, in particular to an early asymptomatic detection method for rice bacterial blight.
背景技术Background technique
随着农业生产需求和检测技术的发展,对于农作物病害早期无症状的检测日益成为研究的热点。当下基于计算机视觉的农作物病害检测和分类主要依赖于病害爆发后展现的颜色、纹理等特征构建的模型来实现。研究人员已经提出了多种基于RGB图像检测感病水稻叶片的外部变化的模型,并且取得了很好的效果。With the development of agricultural production requirements and detection technology, early asymptomatic detection of crop diseases has increasingly become a research hotspot. At present, the detection and classification of crop diseases based on computer vision mainly rely on the model constructed by the color, texture and other features displayed after the outbreak of the disease. Researchers have proposed a variety of models for detecting external changes of susceptible rice leaves based on RGB images, and achieved good results.
其中,基于深度学习技术的农作物病害检测模型越来越引起研究人员的注意。深度学习可以理解为“特征学习”,它是通过多层处理,逐渐将初始的“低层”特征表示转化为“高层”特征表示后,用“简单模型”完成复杂的分类与回归等任务。自2015年起,随着计算机数据处理能力的大幅提升,深度卷积神经网络取得了飞速发展,在农作物病害检测中的应用研究越来越广泛。深度学习工具中各类经典卷积神经网络(例如ResNet、GoogLeNet、VGG等)和自定义卷积神经网络已被广泛应用于农作物病害检测。Among them, the crop disease detection model based on deep learning technology has attracted more and more attention of researchers. Deep learning can be understood as "feature learning". It gradually converts the initial "low-level" feature representation into a "high-level" feature representation through multi-layer processing, and uses "simple models" to complete complex tasks such as classification and regression. Since 2015, with the substantial improvement of computer data processing capabilities, deep convolutional neural networks have achieved rapid development, and their application research in crop disease detection has become more and more extensive. Various classic convolutional neural networks (such as ResNet, GoogLeNet, VGG, etc.) and custom convolutional neural networks in deep learning tools have been widely used in crop disease detection.
目前尽管深度神经网络在从RGB图像检测水稻病害方面取得了巨大成功,但值得注意的是,这些网络无法基于RGB图像检测农作物早期无症状病害。如果不提前给出标签,由于其相似的视觉纹理,很难通过目视检测识别叶片是健康状态还是无症状病害状态。而且就目前的高光谱成像技术的发展状况来看,高光谱图像中海量的光谱信息数据未被充分处理和挖掘,信息处理远远不能满足实时需要,以及当前基于深度学习技术的高光谱病害检测模型忽略了目标尺度的问题(在给定尺寸的图像中,有的病斑比较小,有的病斑比较大),而使用固定尺寸的卷积核会导致较低的识别准确率。因此当前技术中存在无法基于RGB图像检测农作物早期无症状病害,海量的高光谱图像数据的未被充分处理和挖掘难以满足检测实时性需求,使用固定尺寸的卷积核会导致较低的识别准确率的问题。Although deep neural networks have achieved great success in detecting rice diseases from RGB images, it is worth noting that these networks cannot detect early asymptomatic diseases of crops based on RGB images. Without prior labeling, it is difficult to identify leaves in a healthy state or asymptomatic disease state by visual inspection due to their similar visual textures. Moreover, as far as the current development status of hyperspectral imaging technology is concerned, the massive spectral information data in hyperspectral images has not been fully processed and mined, and information processing is far from meeting real-time needs, and the current hyperspectral disease detection based on deep learning technology The model ignores the problem of the target scale (in an image of a given size, some lesions are relatively small, and some lesions are relatively large), and the use of a fixed-size convolution kernel will result in lower recognition accuracy. Therefore, in the current technology, it is impossible to detect early asymptomatic diseases of crops based on RGB images. The massive hyperspectral image data has not been fully processed and mined to meet the real-time detection requirements. The use of fixed-size convolution kernels will lead to lower recognition accuracy. rate problem.
发明内容Contents of the invention
本发明的目的在于提供一种水稻白叶枯病早期无症状检测方法,以解决现有技术中无法基于RGB图像检测农作物早期无症状病害,高光谱图像中海量的光谱信息数据未被充分处理和挖掘难以满足检测实时性需求,使用固定尺寸的卷积核会导致较低的识别准确率的技术问题。The purpose of the present invention is to provide a method for early asymptomatic detection of rice bacterial blight, to solve the problem that in the prior art, it is impossible to detect early asymptomatic diseases of crops based on RGB images, and the massive spectral information data in hyperspectral images is not fully processed. Mining is difficult to meet the real-time requirements of detection, and the use of fixed-size convolution kernels will lead to technical problems of low recognition accuracy.
为解决上述技术问题,本发明具体提供下述技术方案:In order to solve the above technical problems, the present invention specifically provides the following technical solutions:
一种水稻白叶枯病早期无症状检测方法,包括以下步骤:A method for early asymptomatic detection of rice bacterial blight, comprising the following steps:
获取健康水稻叶片和感染白叶枯病害的水稻叶片在预设波段内各个光谱波长处的高光谱图像组成高光谱图像数据集,并利用高光谱图像数据集构建健康水稻叶片和感染白叶枯病害的水稻叶片的高光谱曲线确定出特征波长区间,其中,位于所述特征波长区间内的高光谱图像表征为可区分健康水稻叶片和感染白叶枯病害的水稻叶片的图像特征;Obtain hyperspectral images of healthy rice leaves and rice leaves infected with bacterial blight at various spectral wavelengths within the preset band to form a hyperspectral image dataset, and use the hyperspectral image dataset to construct healthy rice leaves and rice leaves infected with bacterial blight The hyperspectral curve of the rice leaf determines the characteristic wavelength interval, wherein the hyperspectral image located in the characteristic wavelength interval is characterized as an image feature that can distinguish healthy rice leaves from rice leaves infected with bacterial blight;
利用随机森林算法基于位于特征波长区间内的所述高光谱图像对特征波长区间内各个光谱波长进行重要性评分,并在特征波长区间中选取出高重要性评分的光谱波长作为特征敏感波长;Using a random forest algorithm to score the importance of each spectral wavelength in the characteristic wavelength interval based on the hyperspectral image located in the characteristic wavelength interval, and select a spectral wavelength with a high importance score in the characteristic wavelength interval as the characteristic sensitive wavelength;
将特征敏感波长处的高光谱图像作为用于区分高光谱图像所属叶片类别的敏感图像特征,以实现高光谱图像数据集的数据降维;The hyperspectral image at the characteristic sensitive wavelength is used as a sensitive image feature for distinguishing the leaf category to which the hyperspectral image belongs, so as to achieve data dimensionality reduction of the hyperspectral image dataset;
基于敏感图像特征对3DCNN模型训练得到用于水稻白叶枯病早期无症状检测的3DCNN无症状检测模型;The 3DCNN asymptomatic detection model for early asymptomatic detection of rice bacterial blight was obtained by training the 3DCNN model based on sensitive image features;
在3DCNN无症状检测模型中引入多尺度光谱空洞卷积模块进行精度优化得到MS-SDC-3DCNN模型,以实现水稻白叶枯病早期无症状的高精度检测。The MS-SDC-3DCNN model was obtained by introducing a multi-scale spectral atrous convolution module into the 3DCNN asymptomatic detection model for precision optimization, so as to achieve high-precision detection of early asymptomatic rice bacterial blight.
作为本发明的一种优选方案,所述利用高光谱图像数据集构建健康水稻叶片和感染白叶枯病害的水稻叶片的高光谱曲线确定出特征波长区间,包括:As a preferred solution of the present invention, the hyperspectral curves of healthy rice leaves and rice leaves infected with bacterial blight are determined by using the hyperspectral image data set to determine the characteristic wavelength range, including:
将高光谱图像数据集中属于健康水稻叶片的所有高光谱图像绘制在光谱坐标系中得到健康水稻叶片的高光谱曲线,其中,所述光谱二维坐标系中纵坐标为反射强度,横坐标为光谱波长;Draw all hyperspectral images belonging to healthy rice leaves in the hyperspectral image data set in the spectral coordinate system to obtain the hyperspectral curve of healthy rice leaves, wherein the ordinate in the spectral two-dimensional coordinate system is the reflection intensity, and the abscissa is the spectrum wavelength;
将高光谱图像数据集中属于感染白叶枯病害的水稻叶片中无症状水稻叶片的所有高光谱图像绘制在光谱坐标系中得到无症状水稻叶片的高光谱曲线;Draw all the hyperspectral images of the asymptomatic rice leaves in the hyperspectral image data set belonging to the rice leaves infected with bacterial blight in the spectral coordinate system to obtain the hyperspectral curve of the asymptomatic rice leaves;
将高光谱图像数据集中属于感染白叶枯病害的水稻叶片中轻度症状水稻叶片的所有高光谱图像绘制在光谱坐标系中得到轻度症状水稻叶片的高光谱曲线;Draw all hyperspectral images of rice leaves with mild symptoms in the hyperspectral image data set belonging to rice leaves infected with bacterial blight in the spectral coordinate system to obtain the hyperspectral curve of rice leaves with mild symptoms;
将高光谱图像数据集中属于感染白叶枯病害的水稻叶片中中度症状水稻叶片的所有高光谱图像绘制在光谱坐标系中得到中度症状水稻叶片的高光谱曲线;Draw all hyperspectral images of rice leaves with moderate symptoms in the hyperspectral image data set belonging to rice leaves infected with bacterial blight in the spectral coordinate system to obtain a hyperspectral curve of rice leaves with moderate symptoms;
将高光谱图像数据集中属于感染白叶枯病害的水稻叶片中重度症状水稻叶片的所有高光谱图像绘制在光谱坐标系中得到重度症状水稻叶片的高光谱曲线;Draw all hyperspectral images of rice leaves with severe symptoms in the hyperspectral image data set belonging to rice leaves infected with bacterial blight in the spectral coordinate system to obtain hyperspectral curves of rice leaves with severe symptoms;
在健康水稻叶片的高光谱曲线、无症状水稻叶片的高光谱曲线、轻度症状水稻叶片的高光谱曲线、中度症状水稻叶片的高光谱曲线和重度症状水稻叶片的高光谱曲线中进行逐点测算相似度,并将相似度低于预设阈值的所有曲线点进行逐点连接得到特征区分曲线;Point-by-point in hyperspectral curves of healthy rice leaves, hyperspectral curves of asymptomatic rice leaves, hyperspectral curves of mildly symptomatic rice leaves, hyperspectral curves of moderately symptomatic rice leaves, and hyperspectral curves of severely symptomatic rice leaves Calculate the similarity, and connect all the curve points whose similarity is lower than the preset threshold point by point to obtain the characteristic distinction curve;
将特征区分曲线位于的光谱波长区间作为所述特征波长区间。The spectral wavelength interval where the characteristic distinguishing curve is located is taken as the characteristic wavelength interval.
作为本发明的一种优选方案,所述利用随机森林算法基于位于特征波长区间内的所述高光谱图像对特征波长区间内各个光谱波长进行重要性评分,包括:As a preferred solution of the present invention, the random forest algorithm is used to score the importance of each spectral wavelength in the characteristic wavelength interval based on the hyperspectral image located in the characteristic wavelength interval, including:
利用随机森林算法依次对位于特征波长区间内的每个高光谱图像进行随机替换,并计算替换前后的区分高光谱图像所属叶片类别的分类误差作为高光谱图像的重要性评分;Use the random forest algorithm to randomly replace each hyperspectral image in the characteristic wavelength interval in turn, and calculate the classification error before and after the replacement to distinguish the leaf category of the hyperspectral image as the importance score of the hyperspectral image;
将高光谱图像的重要性评分作为高光谱图像对应的光谱波长的重要性评分,并依据光谱波长的重要性评分对光谱波长进行由高到低排序。The importance score of the hyperspectral image is used as the importance score of the spectral wavelength corresponding to the hyperspectral image, and the spectral wavelengths are sorted from high to low according to the importance score of the spectral wavelength.
作为本发明的一种优选方案,所述在特征波长区间中选取出高重要性评分的光谱波长作为特征敏感波长,包括:As a preferred solution of the present invention, selecting a spectral wavelength with a high importance score in the characteristic wavelength interval as a characteristic sensitive wavelength includes:
设定选取阈值,对特征波长区间中排序后的光谱波长的重要性评分由高到低进行求和直到满足所选定阈值为止,将此累加过程中对应的光谱波长作为所述特征敏感波长;Setting the selection threshold, summing the importance scores of the sorted spectral wavelengths in the characteristic wavelength interval from high to low until the selected threshold is met, and using the corresponding spectral wavelength in the accumulation process as the characteristic sensitive wavelength;
其中,选取阈值为0.9。Among them, the selected threshold is 0.9.
作为本发明的一种优选方案,所述3DCNN无症状检测模型的获取,包括:As a preferred solution of the present invention, the acquisition of the 3DCNN asymptomatic detection model includes:
将高光谱图像数据集中位于特征敏感波长处的健康水稻叶片的高光谱图像、高光谱图像数据集中位于特征敏感波长处的感染白叶枯病害的水稻叶片的高光谱图像作为3DCNN模型的输入,将高光谱图像所属叶片类别作为3DCNN模型的输出;The hyperspectral images of healthy rice leaves located at characteristic sensitive wavelengths in the hyperspectral image dataset and the hyperspectral images of rice leaves infected with bacterial blight at characteristic sensitive wavelengths in the hyperspectral image dataset are used as input to the 3DCNN model. The leaf category to which the hyperspectral image belongs is used as the output of the 3DCNN model;
利用3DCNN模型基于所述3DCNN模型的输入和3DCNN模型的输出进行训练得到所述3DCNN无症状检测模型。Using the 3DCNN model to train based on the input of the 3DCNN model and the output of the 3DCNN model to obtain the 3DCNN asymptomatic detection model.
作为本发明的一种优选方案,所述在3DCNN无症状检测模型中引入多尺度光谱空洞卷积模块得到MS-SDC-3DCNN模型,包括:As a preferred solution of the present invention, the introduction of the multi-scale spectral hole convolution module in the 3DCNN asymptomatic detection model to obtain the MS-SDC-3DCNN model includes:
通过选取不同光谱空洞率的三维光谱空洞卷积模块组合为多尺度光谱空洞卷积模块,分别对所述高光谱图像进行光谱特征提取,其中,不同光谱空洞率的三维光谱空洞卷积模块所提取到的光谱特征信息在光谱维度上具有不同的特征尺度;By selecting three-dimensional spectral hole convolution modules with different spectral hole ratios and combining them into multi-scale spectral hole convolution modules, the spectral features of the hyperspectral images are extracted respectively, wherein the three-dimensional spectral hole convolution modules with different spectral hole ratios extract The obtained spectral feature information has different characteristic scales in the spectral dimension;
所述多尺度光谱空洞卷积模块将获得的不同特征尺度的光谱特征信息进行融合,以提高高光谱图像中光谱特征信息的获取;The multi-scale spectral atrous convolution module fuses the obtained spectral feature information of different feature scales to improve the acquisition of spectral feature information in hyperspectral images;
将多尺度光谱空洞卷积模块输出的光谱特征信息作为3DCNN无症状检测模型的输入,以提高3DCNN无症状监测模型的检测性能;The spectral feature information output by the multi-scale spectral hole convolution module is used as the input of the 3DCNN asymptomatic detection model to improve the detection performance of the 3DCNN asymptomatic detection model;
所述三维光谱空洞卷积模块的感受野与光谱空洞率的关系式为:The relationship between the receptive field of the three-dimensional spectral hole convolution module and the spectral hole rate is:
Rf=2×(rSDR-1)×(k-1)+k;R f =2×(r SDR −1)×(k−1)+k;
其中Rf表示三维光谱空洞卷积模块中单个卷积核的感受野;rSDR表示光谱空洞率,k表示卷积核的大小,这里k默认设置为3;Where R f represents the receptive field of a single convolution kernel in the three-dimensional spectral hole convolution module; r SDR represents the spectral hole rate, and k represents the size of the convolution kernel, where k is set to 3 by default;
所述多尺度光谱空洞卷积模块中的特征提取操作由3D卷积来实现,其中,计算3D卷积网络的第i层中第j个特征图上位置(x,y,z)处的值的计算公式为:The feature extraction operation in the multi-scale spectral atrous convolution module is implemented by 3D convolution, wherein the value at the position (x, y, z) on the jth feature map in the i-th layer of the 3D convolutional network is calculated The calculation formula is:
其中,为第i层第j个特征图上位置(x,y,z)处的值,Hi、Wi和Di是第i层卷积核的高度、宽度和深度,f为第i层的激活函数,/>是连接到上一层第m个特征图的卷积核位置(h,w,d)处的权重参数,/>是上一层第m个特征图上位置(x+h,y+w,z+d)处的值,bij是第i层第j个特征图的偏差;in, is the value at position (x, y, z) on the jth feature map of the i-th layer, H i , W i and D i are the height, width and depth of the convolution kernel of the i-th layer, and f is the value of the i-th layer activation function, /> is the weight parameter at the convolution kernel position (h, w, d) connected to the mth feature map of the previous layer, /> is the value at position (x+h, y+w, z+d) on the mth feature map of the previous layer, b ij is the deviation of the jth feature map of the i-th layer;
将多尺度光谱空洞卷积模块的输出作为3DCNN无症状检测模型的输入得到所述MS-SDC-3DCNN模型。The MS-SDC-3DCNN model is obtained by using the output of the multi-scale spectral atrous convolution module as the input of the 3DCNN asymptomatic detection model.
作为本发明的一种优选方案,所述利用MS-SDC-3DCNN模型基于敏感图像特征进行训练,包括:As a preferred solution of the present invention, said utilizing MS-SDC-3DCNN model to train based on sensitive image features, including:
将高光谱图像数据集中位于特征敏感波长处的健康水稻叶片的高光谱图像、高光谱图像数据集中位于特征敏感波长处的感染白叶枯病害的水稻叶片的高光谱图像混合作为训练图像数据集;Mixing hyperspectral images of healthy rice leaves at characteristic sensitive wavelengths in the hyperspectral image data set and hyperspectral images of rice leaves infected with bacterial blight at characteristic sensitive wavelengths in the hyperspectral image data set are used as training image data sets;
将训练图像数据集以8:1:1的比例分为训练集、验证集和测试集,并利用所述训练集、验证集和测试集对所述MS-SDC-3DCNN模型进行训练得到所述检测模型;The training image data set is divided into a training set, a verification set and a test set in a ratio of 8:1:1, and the MS-SDC-3DCNN model is trained using the training set, verification set and test set to obtain the described detection model;
其中,MS-SDC-3DCNN模型使用交叉熵作为损失函数,使用随机梯度下降优化器进行训练,MS-SDC-3DCNN模型的具体参数如下:学习率设置为1×10-3,权重衰减系数设置为1×10-6,动量设为0.95,ε设为1×10-5,epoch设置为50,dropout设置为0.45,MS-SDC-3DCNN输入的高光谱图像尺寸为(9,9,10)。Among them, the MS-SDC-3DCNN model uses cross entropy as the loss function, and uses the stochastic gradient descent optimizer for training. The specific parameters of the MS-SDC-3DCNN model are as follows: the learning rate is set to 1×10 -3 , and the weight decay coefficient is set to 1×10 -6 , momentum is set to 0.95, ε is set to 1×10 -5 , epoch is set to 50, dropout is set to 0.45, and the hyperspectral image size input by MS-SDC-3DCNN is (9,9,10).
作为本发明的一种优选方案,所述MS-SDC-3DCNN模型中设置有残差块,以实现残差块来避免梯度消失问题。As a preferred solution of the present invention, the MS-SDC-3DCNN model is provided with a residual block to implement the residual block to avoid the problem of gradient disappearance.
作为本发明的一种优选方案,所述高光谱图像由31680个像素组成,尺寸为(132,240,10),其中,高光谱图像像素分为三类:健康像素、无症状病害像素和有症状病害像素。As a preferred solution of the present invention, the hyperspectral image consists of 31680 pixels with a size of (132,240,10), wherein the hyperspectral image pixels are divided into three categories: healthy pixels, asymptomatic disease pixels and symptomatic disease pixels pixels.
作为本发明的一种优选方案,所述特征波长区间为450-950nm,所述特征波长区间内包含232个光谱波长。As a preferred solution of the present invention, the characteristic wavelength range is 450-950 nm, and the characteristic wavelength range includes 232 spectral wavelengths.
本发明与现有技术相比较具有如下有益效果:Compared with the prior art, the present invention has the following beneficial effects:
本发明利用随机森林算法对高光谱图像进行降维,确定对病害敏感的波长,降低水稻叶片的无症状检测模型的训练时间,构建多尺度光谱空洞卷积模块,实现在不增加计算量的前提下,扩大神经网络的感受野,以提高特征提取能力,采用多尺度光谱空洞卷积模块对无症状检测模型进行优化,无症状检测模型利用经提取和融合后多个波长分辨率的特征,更有效地使用重要的波长信息,以提高无症状检测模型的检测性能。The invention uses the random forest algorithm to reduce the dimension of the hyperspectral image, determines the wavelength sensitive to the disease, reduces the training time of the asymptomatic detection model of rice leaves, and constructs a multi-scale spectral hole convolution module, which can be achieved without increasing the amount of calculation. Next, expand the receptive field of the neural network to improve the feature extraction ability, and optimize the asymptomatic detection model by using the multi-scale spectral hole convolution module. Efficient use of important wavelength information to improve detection performance of asymptomatic detection models.
附图说明Description of drawings
为了更清楚地说明本发明的实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是示例性的,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图引伸获得其它的实施附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following will briefly introduce the accompanying drawings that are required in the description of the embodiments or the prior art. Apparently, the drawings in the following description are only exemplary, and those skilled in the art can also obtain other implementation drawings according to the provided drawings without creative work.
图1为本发明实施例提供水稻白叶枯病早期无症状检测方法的结构示意图的结构示意图;Fig. 1 provides the structural representation of the structural representation of rice bacterial blight early stage asymptomatic detection method for the embodiment of the present invention;
图2为健康水稻叶片的RGB高光谱图像;Figure 2 is an RGB hyperspectral image of healthy rice leaves;
图3为感染白叶枯病害的无症状水稻叶片的RGB高光谱图像;Figure 3 is an RGB hyperspectral image of asymptomatic rice leaves infected with bacterial blight;
图4为高光谱曲线图;Fig. 4 is hyperspectral graph;
图4(A)为232个波长的重要性评分;Figure 4(A) is the importance score of 232 wavelengths;
图4(B)为232个波长中前50个波长的重要性分数的排名;Figure 4(B) is the ranking of the importance scores of the top 50 wavelengths among the 232 wavelengths;
图5为不同光谱空洞率的三维光谱空洞卷积模块示意图;Fig. 5 is a schematic diagram of a three-dimensional spectral hole convolution module with different spectral hole ratios;
图5(A)为光谱空洞率为1的三维光谱空洞卷积模块示意图;Figure 5 (A) is a schematic diagram of a three-dimensional spectral void convolution module with a spectral void rate of 1;
图5(B)为光谱空洞率为2的三维光谱空洞卷积模块示意图;Figure 5(B) is a schematic diagram of a three-dimensional spectral void convolution module with a spectral void rate of 2;
图5(C)为光谱空洞率为3的三维光谱空洞卷积模块示意图;Figure 5(C) is a schematic diagram of a three-dimensional spectral void convolution module with a spectral void rate of 3;
图6为多尺度光谱空洞卷积模块示意图;Fig. 6 is a schematic diagram of a multi-scale spectral hole convolution module;
图7为MS-SDC-3DCNN模型的神经网络框架图。Figure 7 is a neural network framework diagram of the MS-SDC-3DCNN model.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some, not all, embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts belong to the protection scope of the present invention.
深度神经网络在从RGB图像检测水稻病害方面取得了巨大成功,但值得注意的是,这些网络可能无法为基于RGB图像的早期无症状BLB病害检测产生正确的结果。水稻白叶枯病害胁迫下健康和无症状叶片的高光谱图像的RGB图示如图2和图3所示。如果不提前给出标签,由于其相似的视觉纹理,很难通过目视检测识别叶片是健康的(图2)还是无症状的(图3)。Deep neural networks have achieved great success in detecting rice diseases from RGB images, but it is worth noting that these networks may not produce correct results for early asymptomatic BLB disease detection based on RGB images. RGB illustrations of hyperspectral images of healthy and asymptomatic leaves under rice bacterial blight stress are shown in Figures 2 and 3. Without prior labeling, it is difficult to identify by visual inspection whether leaves are healthy (Fig. 2) or asymptomatic (Fig. 3) due to their similar visual textures.
然而,病害会使得水稻叶片内部结构发生变化,水稻叶片的这种内部变化导致叶片对于不同波段的反射率也发生变化,使得利用高光谱成像技术来检测水稻病害的无症状感染得以实现。因此本发明基于水稻叶片的高光谱图像结合深度学习模型进行水稻白叶枯病害早期无症状检测方法(一种水稻白叶枯病早期无症状检测方法),利用高光谱图像实现无症状检测,有利于预防和控制白叶枯病害。However, the disease will change the internal structure of rice leaves, and this internal change of rice leaves will lead to changes in the reflectance of the leaves for different bands, making it possible to use hyperspectral imaging technology to detect asymptomatic infections of rice diseases. Therefore, the present invention is based on the hyperspectral image of rice leaves combined with a deep learning model for early asymptomatic detection of rice bacterial blight (a method for early asymptomatic detection of rice bacterial blight), using hyperspectral images to achieve asymptomatic detection, there is Conducive to the prevention and control of bacterial blight.
如图1所示,本发明提供了一种水稻白叶枯病早期无症状检测方法,包括以下步骤:As shown in Figure 1, the present invention provides a kind of asymptomatic detection method of rice bacterial blight early stage, comprises the following steps:
获取健康水稻叶片和感染白叶枯病害的水稻叶片在预设波段内各个光谱波长处的高光谱图像组成高光谱图像数据集,并利用高光谱图像数据集构建健康水稻叶片和感染白叶枯病害的水稻叶片的高光谱曲线确定出特征波长区间,其中,位于所述特征波长区间内的高光谱图像表征为可区分健康水稻叶片和感染白叶枯病害的水稻叶片的图像特征;Obtain hyperspectral images of healthy rice leaves and rice leaves infected with bacterial blight at various spectral wavelengths within the preset band to form a hyperspectral image dataset, and use the hyperspectral image dataset to construct healthy rice leaves and rice leaves infected with bacterial blight The hyperspectral curve of the rice leaf determines the characteristic wavelength interval, wherein the hyperspectral image located in the characteristic wavelength interval is characterized as an image feature that can distinguish healthy rice leaves from rice leaves infected with bacterial blight;
利用随机森林算法基于位于特征波长区间内的所述高光谱图像对特征波长区间内各个光谱波长进行重要性评分,并在特征波长区间中选取出高重要性评分的光谱波长作为特征敏感波长;Using a random forest algorithm to score the importance of each spectral wavelength in the characteristic wavelength interval based on the hyperspectral image located in the characteristic wavelength interval, and select a spectral wavelength with a high importance score in the characteristic wavelength interval as the characteristic sensitive wavelength;
将特征敏感波长处的高光谱图像作为用于区分高光谱图像所属叶片类别的敏感图像特征,以实现高光谱图像数据集的数据降维;The hyperspectral image at the characteristic sensitive wavelength is used as a sensitive image feature for distinguishing the leaf category to which the hyperspectral image belongs, so as to achieve data dimensionality reduction of the hyperspectral image dataset;
基于敏感图像特征对3DCNN模型训练得到用于水稻白叶枯病早期无症状检测的3DCNN无症状检测模型;The 3DCNN asymptomatic detection model for early asymptomatic detection of rice bacterial blight was obtained by training the 3DCNN model based on sensitive image features;
在3DCNN无症状检测模型中引入多尺度光谱空洞卷积模块进行精度优化得到MS-SDC-3DCNN模型,以实现水稻白叶枯病早期无症状的高精度检测。The MS-SDC-3DCNN model was obtained by introducing a multi-scale spectral atrous convolution module into the 3DCNN asymptomatic detection model for precision optimization, so as to achieve high-precision detection of early asymptomatic rice bacterial blight.
病害会使得水稻叶片内部结构发生变化,水稻叶片的这种内部变化导致叶片对于不同波段的反射率也发生变化,因此在同一光谱波长处健康水稻叶片和感染白叶枯病害的水稻叶片的高光谱图像对应的反射率不同,反射率的不同就可以区分出高光谱图像是属于健康水稻叶片还是感染白叶枯病害的水稻叶片,因此通过对比各种水稻叶片的高光谱曲线能够筛选出区分出高光谱图像所属水稻叶片类别的光谱波段,即表明该光谱波段内的高光谱图像包含有区分水稻叶片类别的图像特征,因此,通过高光谱曲线的绘制,能够获取可区分健康水稻叶片和感染白叶枯病害的水稻叶片的图像特征,具体如下:The disease will change the internal structure of rice leaves, and this internal change of rice leaves will lead to changes in the reflectance of the leaves for different wavelength bands. Therefore, the hyperspectral data of healthy rice leaves and rice leaves infected with bacterial blight at the same spectral wavelength The reflectance corresponding to the image is different, and the difference in reflectance can distinguish whether the hyperspectral image belongs to healthy rice leaves or rice leaves infected with bacterial blight. The spectral band of the rice leaf category that the spectral image belongs to indicates that the hyperspectral image in this spectral band contains image features that distinguish the rice leaf category. Therefore, through the drawing of the hyperspectral curve, it is possible to obtain a The image features of rice leaves with blight disease are as follows:
所述利用高光谱图像数据集构建健康水稻叶片和感染白叶枯病害的水稻叶片的高光谱曲线确定出特征波长区间,包括:The hyperspectral curves of healthy rice leaves and rice leaves infected with bacterial blight are determined by using the hyperspectral image data set to determine the characteristic wavelength range, including:
将高光谱图像数据集中属于健康水稻叶片的所有高光谱图像绘制在光谱坐标系中得到健康水稻叶片的高光谱曲线,其中,所述光谱二维坐标系中纵坐标为反射强度,横坐标为光谱波长;Draw all hyperspectral images belonging to healthy rice leaves in the hyperspectral image data set in the spectral coordinate system to obtain the hyperspectral curve of healthy rice leaves, wherein the ordinate in the spectral two-dimensional coordinate system is the reflection intensity, and the abscissa is the spectrum wavelength;
将高光谱图像数据集中属于感染白叶枯病害的水稻叶片中无症状水稻叶片的所有高光谱图像绘制在光谱坐标系中得到无症状水稻叶片的高光谱曲线;Draw all the hyperspectral images of the asymptomatic rice leaves in the hyperspectral image data set belonging to the rice leaves infected with bacterial blight in the spectral coordinate system to obtain the hyperspectral curve of the asymptomatic rice leaves;
将高光谱图像数据集中属于感染白叶枯病害的水稻叶片中轻度症状水稻叶片的所有高光谱图像绘制在光谱坐标系中得到轻度症状水稻叶片的高光谱曲线;Draw all hyperspectral images of rice leaves with mild symptoms in the hyperspectral image data set belonging to rice leaves infected with bacterial blight in the spectral coordinate system to obtain the hyperspectral curve of rice leaves with mild symptoms;
将高光谱图像数据集中属于感染白叶枯病害的水稻叶片中中度症状水稻叶片的所有高光谱图像绘制在光谱坐标系中得到中度症状水稻叶片的高光谱曲线;Draw all hyperspectral images of rice leaves with moderate symptoms in the hyperspectral image data set belonging to rice leaves infected with bacterial blight in the spectral coordinate system to obtain a hyperspectral curve of rice leaves with moderate symptoms;
将高光谱图像数据集中属于感染白叶枯病害的水稻叶片中重度症状水稻叶片的所有高光谱图像绘制在光谱坐标系中得到重度症状水稻叶片的高光谱曲线;Draw all hyperspectral images of rice leaves with severe symptoms in the hyperspectral image data set belonging to rice leaves infected with bacterial blight in the spectral coordinate system to obtain hyperspectral curves of rice leaves with severe symptoms;
在健康水稻叶片的高光谱曲线、无症状水稻叶片的高光谱曲线、轻度症状水稻叶片的高光谱曲线、中度症状水稻叶片的高光谱曲线和重度症状水稻叶片的高光谱曲线中进行逐点测算相似度,并将相似度低于预设阈值的所有曲线点进行逐点连接得到特征区分曲线;Point-by-point in hyperspectral curves of healthy rice leaves, hyperspectral curves of asymptomatic rice leaves, hyperspectral curves of mildly symptomatic rice leaves, hyperspectral curves of moderately symptomatic rice leaves, and hyperspectral curves of severely symptomatic rice leaves Calculate the similarity, and connect all the curve points whose similarity is lower than the preset threshold point by point to obtain the characteristic distinction curve;
将特征区分曲线位于的光谱波长区间作为所述特征波长区间。The spectral wavelength interval where the characteristic distinguishing curve is located is taken as the characteristic wavelength interval.
如图4,图4(A),图4(B)所示为白叶枯病害胁迫下不同病害状态的水稻叶片378到1033nm波长之间的高光谱曲线图,因此区分水稻叶片所述类别的光谱特征信息涵盖在为378-1033nm光谱波长之间。实验所使用的高光谱成像系统的光谱分辨率为2.14nm,拍摄获取的高光谱波长区间为378.28nm-1033.05nm,为了减少光谱两端因仪器、环境而产生的噪声的影响,对378-1033nm光谱波长区间进行进一步截取为451.276nm-949.664nm,因此本发明将特征波长区间设置为450-950nm,既包含区分水稻叶片所述类别的光谱特征信息,又减少光谱两端因仪器、环境而产生的噪声的影响。As shown in Fig. 4, Fig. 4 (A), Fig. 4 (B) shows the hyperspectral curve graph between 378 to 1033nm wavelengths of rice leaves in different disease states under bacterial blight stress, so distinguish the described category of rice leaves Spectral characteristic information is covered between 378-1033nm spectral wavelength. The spectral resolution of the hyperspectral imaging system used in the experiment is 2.14nm, and the hyperspectral wavelength range obtained by shooting is 378.28nm-1033.05nm. The spectral wavelength interval is further intercepted to 451.276nm-949.664nm, so the present invention sets the characteristic wavelength interval to 450-950nm, which not only includes the spectral characteristic information for distinguishing the categories of rice leaves, but also reduces the spectrum at both ends due to instruments and environments. the influence of noise.
在特征波长区间378-1033nm中获取所有高光谱图像,截取450至950nm光谱波长处的高光谱图像以去除特征波长区间两端的噪声,同时更好的区分水稻叶片所述类别的光谱特征信息。所获得的高光谱图像数据其特点是相邻波长之间具有高维冗余。过多的冗余光谱信息给检测方法和计算复杂性带来了巨大挑战。因此,有必要通过降维方法压缩数据量,在不丢弃有效特征光谱信息的基础上降低后续处理的成本,本发明利用随机森林算法进行区分水稻叶片所属类别的光谱特征信息降维,降低原始高光谱数据的维数,从而降低计算复杂度,具体步骤如下:All hyperspectral images were obtained in the characteristic wavelength range of 378-1033nm, and hyperspectral images at spectral wavelengths from 450 to 950nm were intercepted to remove noise at both ends of the characteristic wavelength range and to better distinguish the spectral characteristic information of the rice leaf category. The acquired hyperspectral image data is characterized by high dimensional redundancy between adjacent wavelengths. Excessive redundant spectral information brings great challenges to detection methods and computational complexity. Therefore, it is necessary to compress the amount of data through a dimensionality reduction method and reduce the cost of subsequent processing without discarding effective characteristic spectral information. The dimensionality of the spectral data, thereby reducing the computational complexity, the specific steps are as follows:
所述利用随机森林算法基于位于特征波长区间内的所述高光谱图像对特征波长区间内各个光谱波长进行重要性评分,包括:The random forest algorithm is used to score the importance of each spectral wavelength in the characteristic wavelength interval based on the hyperspectral image located in the characteristic wavelength interval, including:
利用随机森林算法依次对位于特征波长区间内的每个高光谱图像进行随机替换,并计算替换前后的区分高光谱图像所属叶片类别的分类误差作为高光谱图像的重要性评分;Use the random forest algorithm to randomly replace each hyperspectral image in the characteristic wavelength interval in turn, and calculate the classification error before and after the replacement to distinguish the leaf category of the hyperspectral image as the importance score of the hyperspectral image;
将高光谱图像的重要性评分作为高光谱图像对应的光谱波长的重要性评分,并依据光谱波长的重要性评分对光谱波长进行由高到低排序。The importance score of the hyperspectral image is used as the importance score of the spectral wavelength corresponding to the hyperspectral image, and the spectral wavelengths are sorted from high to low according to the importance score of the spectral wavelength.
设定选取阈值,对特征波长区间中排序后的光谱波长的重要性评分由高到低进行求和直到满足所选定阈值为止,将此累加过程中对应的光谱波长作为所述特征敏感波长;Setting the selection threshold, summing the importance scores of the sorted spectral wavelengths in the characteristic wavelength interval from high to low until the selected threshold is met, and using the corresponding spectral wavelength in the accumulation process as the characteristic sensitive wavelength;
其中,选取阈值为0.9。Among them, the selected threshold is 0.9.
利用RF对450至950nm中原始232个波长的重要性得分进行排序(图4(A))。可以看出,重要性得分最高的波长主要分布在500–700nm。图4(B)清楚地表明,重要性得分最高的前10个波长为547.2nm、534.5nm、551.4nm、566.2nm、697.4nm、530.3nm、693.0nm、543.0nm、538.7nm和568.4nm,将重要性得分最高的前10个波长为547.2nm、534.5nm、551.4nm、566.2nm、697.4nm、530.3nm、693.0nm、543.0nm、538.7nm和568.4nm作为敏感特征波长,即区分水稻叶片类别的最少数量的特征,并能保证检测模型的检测准确度。将前10个敏感特征波长对应的高光谱图像作为3DCNN深度学习模型的输入,训练3DCNN检测模型,对水稻白叶枯早期无症状病害进行检测。The importance scores of the original 232 wavelengths from 450 to 950 nm were ranked using RF (Fig. 4(A)). It can be seen that the wavelengths with the highest importance scores are mainly distributed in 500–700nm. Figure 4(B) clearly shows that the top 10 wavelengths with the highest importance scores are 547.2nm, 534.5nm, 551.4nm, 566.2nm, 697.4nm, 530.3nm, 693.0nm, 543.0nm, 538.7nm and 568.4nm, and the The top 10 wavelengths with the highest importance scores are 547.2nm, 534.5nm, 551.4nm, 566.2nm, 697.4nm, 530.3nm, 693.0nm, 543.0nm, 538.7nm and 568.4nm as the sensitive characteristic wavelengths, that is, the wavelengths for distinguishing rice leaf categories. The minimum number of features can guarantee the detection accuracy of the detection model. The hyperspectral images corresponding to the first 10 sensitive characteristic wavelengths were used as the input of the 3DCNN deep learning model, and the 3DCNN detection model was trained to detect the early asymptomatic disease of rice bacterial blight.
所述3DCNN无症状检测模型的获取,包括:The acquisition of the 3DCNN asymptomatic detection model includes:
将高光谱图像数据集中位于特征敏感波长处的健康水稻叶片的高光谱图像、高光谱图像数据集中位于特征敏感波长处的感染白叶枯病害的水稻叶片的高光谱图像作为3DCNN模型的输入,将高光谱图像所属叶片类别作为3DCNN模型的输出;The hyperspectral images of healthy rice leaves located at characteristic sensitive wavelengths in the hyperspectral image dataset and the hyperspectral images of rice leaves infected with bacterial blight at characteristic sensitive wavelengths in the hyperspectral image dataset are used as input to the 3DCNN model. The leaf category to which the hyperspectral image belongs is used as the output of the 3DCNN model;
利用3DCNN模型基于所述3DCNN模型的输入和3DCNN模型的输出进行训练得到所述3DCNN无症状检测模型。Using the 3DCNN model to train based on the input of the 3DCNN model and the output of the 3DCNN model to obtain the 3DCNN asymptomatic detection model.
在水稻叶片病害高光谱图像分析过程中,由于病害入侵的程度和水稻植株对病害的抵抗作用不同,选定的感兴趣区域中,水稻叶片高光谱图像中包含大量的光谱信息,这使得采用固定尺寸的卷积核构建的检测模型难以获得充分的光谱信息,从而导致识别准确率降低,因此本发明在3DCNN无症状检测模型中引入了三维光谱空洞卷积模块,以扩展光谱维度的感受野,利用三维光谱空洞卷积模块分别提取不同尺度的光谱特征,最后再进行融合计算,以获得更为丰富的光谱信息,提升模型检测性能,具体如下:In the process of rice leaf disease hyperspectral image analysis, due to the degree of disease invasion and the resistance of rice plants to disease, in the selected region of interest, the rice leaf hyperspectral image contains a large amount of spectral information, which makes the use of fixed It is difficult to obtain sufficient spectral information for a detection model constructed with a convolution kernel of a large size, which leads to a reduction in recognition accuracy. Therefore, the present invention introduces a three-dimensional spectral hole convolution module in the 3DCNN asymptomatic detection model to expand the receptive field of the spectral dimension. Use the three-dimensional spectral atrous convolution module to extract spectral features of different scales, and finally perform fusion calculations to obtain richer spectral information and improve model detection performance, as follows:
所述在3DCNN无症状检测模型中引入多尺度光谱空洞卷积模块得到MS-SDC-3DCNN模型,包括:The MS-SDC-3DCNN model is obtained by introducing a multi-scale spectral hole convolution module in the 3DCNN asymptomatic detection model, including:
通过选取不同光谱空洞率的三维光谱空洞卷积模块组合为多尺度光谱空洞卷积模块,分别对所述高光谱图像进行光谱特征提取,其中,不同光谱空洞率的三维光谱空洞卷积模块所提取到的光谱特征信息在光谱维度上具有不同的特征尺度;By selecting three-dimensional spectral hole convolution modules with different spectral hole ratios and combining them into multi-scale spectral hole convolution modules, the spectral features of the hyperspectral images are extracted respectively, wherein the three-dimensional spectral hole convolution modules with different spectral hole ratios extract The obtained spectral feature information has different characteristic scales in the spectral dimension;
所述多尺度光谱空洞卷积模块将获得的不同特征尺度的光谱特征信息进行融合,以提高高光谱图像中光谱特征信息的获取;The multi-scale spectral atrous convolution module fuses the obtained spectral feature information of different feature scales to improve the acquisition of spectral feature information in hyperspectral images;
将多尺度光谱空洞卷积模块输出的光谱特征信息作为3DCNN无症状检测模型的输入,以提高3DCNN无症状监测模型的检测性能;The spectral feature information output by the multi-scale spectral hole convolution module is used as the input of the 3DCNN asymptomatic detection model to improve the detection performance of the 3DCNN asymptomatic detection model;
所述三维光谱空洞卷积模块的感受野与光谱空洞率的关系式为:The relationship between the receptive field of the three-dimensional spectral hole convolution module and the spectral hole rate is:
Rf=2×(rSDR-1)×(k-1)+k;R f =2×(r SDR −1)×(k−1)+k;
其中Rf表示三维光谱空洞卷积模块中单个卷积核的感受野;rSDR表示光谱空洞率,k表示卷积核的大小,这里k默认设置为3;Where R f represents the receptive field of a single convolution kernel in the three-dimensional spectral hole convolution module; r SDR represents the spectral hole rate, and k represents the size of the convolution kernel, where k is set to 3 by default;
图5(A)中黑色cube为尺寸为(3,3,3)的卷积核,光谱空洞率SDR=1说明不添加空洞,感受野尺寸也为(3,3,3)。图5(B)中光谱空洞率SDR=2,黑色cube为卷积核添加空洞之后的位置,白色cube为添加空洞之后卷积核感受野的尺寸(7,7,7)。图5(C)中光谱空洞率SDR=3,黑色cube为卷积核添加空洞之后的位置,白色cube为添加空洞之后卷积核感受野的尺寸(11,11,11)。黑色立方体代表卷积核,白色立方体覆盖感受野。即当SDR设置为1、2和3时,三维光谱空洞卷积模块中卷积核的感受野分别为3×3×3、7×7×7和11×11×11。不同光谱空洞率的三维光谱空洞卷积模块组合为多尺度光谱空洞卷积模块,实现分别对所述高光谱图像进行光谱特征提取。In Figure 5(A), the black cube is a convolution kernel with a size of (3,3,3), and the spectral hole ratio SDR=1 indicates that no hole is added, and the size of the receptive field is also (3,3,3). In Figure 5(B), the spectral hole ratio SDR=2, the black cube is the position of the convolution kernel after the hole is added, and the white cube is the size of the convolution kernel receptive field after the hole is added (7,7,7). In Figure 5(C), the spectral hole ratio SDR=3, the black cube is the position of the convolution kernel after the hole is added, and the white cube is the size of the convolution kernel receptive field after the hole is added (11, 11, 11). Black cubes represent convolution kernels, and white cubes cover receptive fields. That is, when the SDR is set to 1, 2, and 3, the receptive fields of the convolution kernels in the three-dimensional spectral atrous convolution module are 3×3×3, 7×7×7, and 11×11×11, respectively. The three-dimensional spectral hole convolution modules with different spectral hole ratios are combined into a multi-scale spectral hole convolution module, which realizes the spectral feature extraction of the hyperspectral image respectively.
具体的,所述多尺度光谱空洞卷积模块中的特征提取操作由3D卷积来实现,其中,计算3D卷积网络的第i层中第j个特征图上位置(x,y,z)处的值的计算公式为:Specifically, the feature extraction operation in the multi-scale spectral atrous convolution module is implemented by 3D convolution, wherein the position (x, y, z) on the jth feature map in the i-th layer of the 3D convolutional network is calculated The formula for calculating the value at is:
其中,为第i层第j个特征图上位置(x,y,z)处的值,Hi、Wi和Di是第i层卷积核的高度、宽度和深度,f为第i层的激活函数,/>是连接到上一层第m个特征图的卷积核位置(h,w,d)处的权重参数,/>是上一层第m个特征图上位置(x+h,y+w,z+d)处的值,bij是第i层第j个特征图的偏差。in, is the value at position (x, y, z) on the jth feature map of the i-th layer, H i , W i and D i are the height, width and depth of the convolution kernel of the i-th layer, and f is the value of the i-th layer activation function, /> is the weight parameter at the convolution kernel position (h, w, d) connected to the mth feature map of the previous layer, /> is the value at position (x+h, y+w, z+d) on the mth feature map of the previous layer, b ij is the deviation of the jth feature map of the i-th layer.
将多尺度光谱空洞卷积模块的输出作为3DCNN无症状检测模型的输入得到所述MS-SDC-3DCNN模型。The MS-SDC-3DCNN model is obtained by using the output of the multi-scale spectral atrous convolution module as the input of the 3DCNN asymptomatic detection model.
所述利用MS-SDC-3DCNN模型基于敏感图像特征进行训练,包括:Described utilizing MS-SDC-3DCNN model to carry out training based on sensitive image features, including:
将高光谱图像数据集中位于特征敏感波长处的健康水稻叶片的高光谱图像、高光谱图像数据集中位于特征敏感波长处的感染白叶枯病害的水稻叶片的高光谱图像混合作为训练图像数据集;Mixing hyperspectral images of healthy rice leaves at characteristic sensitive wavelengths in the hyperspectral image data set and hyperspectral images of rice leaves infected with bacterial blight at characteristic sensitive wavelengths in the hyperspectral image data set are used as training image data sets;
将训练图像数据集以8:1:1的比例分为训练集、验证集和测试集,并利用所述训练集、验证集和测试集对所述MS-SDC-3DCNN模型进行训练得到所述检测模型;The training image data set is divided into a training set, a verification set and a test set in a ratio of 8:1:1, and the MS-SDC-3DCNN model is trained using the training set, verification set and test set to obtain the described detection model;
其中,MS-SDC-3DCNN模型使用交叉熵作为损失函数,使用随机梯度下降优化器进行训练,MS-SDC-3DCNN模型的具体参数如下:学习率设置为1×10-3,权重衰减系数设置为1×10-6,动量设为0.95,ε设为1×10-5,epoch设置为50,dropout设置为0.45,MS-SDC-3DCNN输入的高光谱图像尺寸为(9,9,10)。Among them, the MS-SDC-3DCNN model uses cross entropy as the loss function, and uses the stochastic gradient descent optimizer for training. The specific parameters of the MS-SDC-3DCNN model are as follows: the learning rate is set to 1×10 -3 , and the weight decay coefficient is set to 1×10 -6 , momentum is set to 0.95, ε is set to 1×10 -5 , epoch is set to 50, dropout is set to 0.45, and the hyperspectral image size input by MS-SDC-3DCNN is (9,9,10).
所述MS-SDC-3DCNN模型中设置有残差块,以实现残差块来避免梯度消失问题。The MS-SDC-3DCNN model is provided with a residual block to implement the residual block to avoid the problem of gradient disappearance.
所述高光谱图像由31680个像素组成,尺寸为(132,240,10),其中,高光谱图像像素分为三类:健康像素、无症状病害像素和有症状病害像素。将所述高光谱图像裁剪成尺寸为(9,9,10)的图像块作为MS-SDC-3DCNN输入的高光谱图像尺寸,并将每个图像块中心像素所对应的分类标签分配给该图像块。The hyperspectral image consists of 31680 pixels, and the size is (132, 240, 10). Among them, the hyperspectral image pixels are divided into three categories: healthy pixels, asymptomatic disease pixels and symptomatic disease pixels. Crop the hyperspectral image into image blocks of size (9,9,10) as the hyperspectral image size input by MS-SDC-3DCNN, and assign the classification label corresponding to the central pixel of each image block to the image piece.
所述特征波长区间为450-950nm,所述特征波长区间内包含232个光谱波长。实验所使用的高光谱成像系统的光谱分辨率为2.14nm,选取的特征波长区间为450-950nm,把波长区间除以光谱分辨率可计算得到波长区间内共包含232个波长。The characteristic wavelength range is 450-950nm, and the characteristic wavelength range includes 232 spectral wavelengths. The spectral resolution of the hyperspectral imaging system used in the experiment is 2.14nm, and the selected characteristic wavelength range is 450-950nm. Dividing the wavelength range by the spectral resolution can be calculated to include 232 wavelengths in the wavelength range.
本发明利用随机森林算法对高光谱图像进行降维,确定对病害敏感的波长,降低水稻叶片的无症状检测模型的训练时间,构建多尺度光谱空洞卷积模块,实现在不增加计算量的前提下,扩大神经网络的感受野,以提高特征提取能力,采用多尺度光谱空洞卷积模块对无症状检测模型进行优化,无症状检测模型利用经提取和融合后多个波长分辨率的特征,更有效地使用重要的波长信息,以提高无症状检测模型的检测性能。The invention uses the random forest algorithm to reduce the dimension of the hyperspectral image, determines the wavelength sensitive to the disease, reduces the training time of the asymptomatic detection model of rice leaves, and constructs a multi-scale spectral hole convolution module, which can be achieved without increasing the amount of calculation. Next, expand the receptive field of the neural network to improve the feature extraction ability, and optimize the asymptomatic detection model by using the multi-scale spectral hole convolution module. Efficient use of important wavelength information to improve detection performance of asymptomatic detection models.
以上实施例仅为本申请的示例性实施例,不用于限制本申请,本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内,对本申请做出各种修改或等同替换,这种修改或等同替换也应视为落在本申请的保护范围内。The above embodiments are only exemplary embodiments of the present application, and are not intended to limit the present application, and the protection scope of the present application is defined by the claims. Those skilled in the art may make various modifications or equivalent replacements to the present application within the spirit and protection scope of the present application, and such modifications or equivalent replacements shall also be deemed to fall within the protection scope of the present application.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310330513.8A CN116416524A (en) | 2023-03-30 | 2023-03-30 | Early-stage asymptomatic detection method for bacterial leaf blight of rice |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310330513.8A CN116416524A (en) | 2023-03-30 | 2023-03-30 | Early-stage asymptomatic detection method for bacterial leaf blight of rice |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116416524A true CN116416524A (en) | 2023-07-11 |
Family
ID=87057603
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310330513.8A Pending CN116416524A (en) | 2023-03-30 | 2023-03-30 | Early-stage asymptomatic detection method for bacterial leaf blight of rice |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116416524A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118608877A (en) * | 2024-08-08 | 2024-09-06 | 四川省林业和草原调查规划院(四川省林业和草原生态环境监测中心) | A tree screening management system based on image technology |
| CN119360237A (en) * | 2024-09-25 | 2025-01-24 | 宁波大学 | Inversion method of rice bacterial blight severity based on UAV hyperspectral images |
-
2023
- 2023-03-30 CN CN202310330513.8A patent/CN116416524A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118608877A (en) * | 2024-08-08 | 2024-09-06 | 四川省林业和草原调查规划院(四川省林业和草原生态环境监测中心) | A tree screening management system based on image technology |
| CN119360237A (en) * | 2024-09-25 | 2025-01-24 | 宁波大学 | Inversion method of rice bacterial blight severity based on UAV hyperspectral images |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ghaderizadeh et al. | Hyperspectral image classification using a hybrid 3D-2D convolutional neural networks | |
| Zhang et al. | Classification modeling method for near‐infrared spectroscopy of tobacco based on multimodal convolution neural networks | |
| CN113095409B (en) | Hyperspectral Image Classification Method Based on Attention Mechanism and Weight Sharing | |
| US20230334829A1 (en) | Hyperspectral image classification method based on context-rich networks | |
| CN112446388A (en) | Multi-category vegetable seedling identification method and system based on lightweight two-stage detection model | |
| CN111291826B (en) | A pixel-by-pixel classification method for multi-source remote sensing images based on correlation fusion network | |
| CN116416524A (en) | Early-stage asymptomatic detection method for bacterial leaf blight of rice | |
| CN108052911A (en) | Multi-modal remote sensing image high-level characteristic integrated classification method based on deep learning | |
| CN114429564B (en) | A collaborative classification method for hyperspectral and LiDAR data based on dual branches | |
| CN114694178A (en) | Method and system for monitoring safety helmet in power operation based on fast-RCNN algorithm | |
| CN111652273B (en) | Deep learning-based RGB-D image classification method | |
| CN107977683A (en) | Joint SAR target identification methods based on convolution feature extraction and machine learning | |
| CN114937182B (en) | Image emotion distribution prediction method based on emotion wheel and convolutional neural network | |
| CN112434662B (en) | Tea leaf scab automatic identification algorithm based on multi-scale convolutional neural network | |
| CN111222545B (en) | Image classification method based on linear programming incremental learning | |
| CN112364974A (en) | Improved YOLOv3 algorithm based on activation function | |
| CN114694042A (en) | A camouflaged person target detection method based on improved Scaled-YOLOv4 | |
| CN115527193A (en) | Chinese medicinal material type identification method | |
| CN115908907A (en) | A hyperspectral remote sensing image classification method and system | |
| CN118247558A (en) | Plant disease quantitative inversion method, device and equipment for multi-mode deep learning | |
| CN118135392A (en) | Remote sensing image detection method based on dual-temporal interactive enhanced CNN-Transformer | |
| Zhao et al. | BIHAF-Net: Bilateral interactive hierarchical adaptive fusion network for collaborative classification of hyperspectral and LiDAR data | |
| CN112766161A (en) | Hyperspectral target detection method based on integrated constraint multi-example learning | |
| CN115206455B (en) | Deep neural network-based rare earth element component content prediction method and system | |
| CN116563600A (en) | A multi-label recognition method for Chinese food based on multi-scale fusion and attention mechanism |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |