CN112668404B - Effective identification method for soybean diseases and insect pests - Google Patents
Effective identification method for soybean diseases and insect pests Download PDFInfo
- Publication number
- CN112668404B CN112668404B CN202011455125.5A CN202011455125A CN112668404B CN 112668404 B CN112668404 B CN 112668404B CN 202011455125 A CN202011455125 A CN 202011455125A CN 112668404 B CN112668404 B CN 112668404B
- Authority
- CN
- China
- Prior art keywords
- diseases
- soybean
- value
- insect pests
- plants
- 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.)
- Active
Links
- 201000010099 disease Diseases 0.000 title claims abstract description 301
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 title claims abstract description 301
- 241000607479 Yersinia pestis Species 0.000 title claims abstract description 295
- 241000238631 Hexapoda Species 0.000 title claims abstract description 230
- 244000068988 Glycine max Species 0.000 title claims abstract description 194
- 235000010469 Glycine max Nutrition 0.000 title claims abstract description 150
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000000605 extraction Methods 0.000 claims abstract description 63
- 241000196324 Embryophyta Species 0.000 claims description 237
- 210000000056 organ Anatomy 0.000 claims description 26
- 238000005259 measurement Methods 0.000 claims description 17
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 12
- 239000000284 extract Substances 0.000 claims description 12
- 238000012360 testing method Methods 0.000 claims description 12
- 230000017260 vegetative to reproductive phase transition of meristem Effects 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 3
- 238000012821 model calculation Methods 0.000 claims description 3
- 230000000877 morphologic effect Effects 0.000 claims description 3
- 239000003086 colorant Substances 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 4
- 230000000007 visual effect Effects 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 238000002845 discoloration Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 235000015097 nutrients Nutrition 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 208000035143 Bacterial infection Diseases 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 208000031888 Mycoses Diseases 0.000 description 1
- 229920002472 Starch Polymers 0.000 description 1
- 208000025865 Ulcer Diseases 0.000 description 1
- 241000700605 Viruses Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000012271 agricultural production Methods 0.000 description 1
- 238000010923 batch production Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 208000006278 hypochromic anemia Diseases 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000001338 necrotic effect Effects 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 244000052769 pathogen Species 0.000 description 1
- 230000001717 pathogenic effect Effects 0.000 description 1
- 230000000243 photosynthetic effect Effects 0.000 description 1
- 231100000572 poisoning Toxicity 0.000 description 1
- 230000000607 poisoning effect Effects 0.000 description 1
- 235000019698 starch Nutrition 0.000 description 1
- 239000008107 starch Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 231100000397 ulcer Toxicity 0.000 description 1
- 238000011179 visual inspection Methods 0.000 description 1
- 238000004383 yellowing Methods 0.000 description 1
Images
Classifications
-
- 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
- Image Processing (AREA)
- Catching Or Destruction (AREA)
Abstract
Description
技术领域technical field
本发明涉及一种大豆病虫害的有效鉴定方法。The invention relates to an effective identification method for soybean diseases and insect pests.
背景技术Background technique
大豆作为我国主要的粮食作物和经济作物,在国民经济生产中,尤其是农业生产中的地位非常重要。大豆病虫害不仅是我国,也是全世界范围内最主要的农业灾害之一,截至目前,世界上已有报道中涉及大豆的病虫害种类就高达120多种,包括真菌病害、细菌病害、病毒病害、昆虫病害和寄生植物病害等不同类型,不同的大豆病虫害的方式、种类和影响不尽相同,但不同的病虫害都会严重的影响大豆的生长,使得大豆大幅减产甚至绝产。高程度、大范围的病虫害对我国乃至世界各国的国民经济,特别是以大豆农业的生产种植为主的国家的国民经济造成巨大的损失。大豆的病虫害,具体体现遭受病虫害的农作物平均减产1到2成,严重的甚至绝产;在遭受病虫害后,农作物的口感和品质都会有不同程度的下降,淀粉和蛋白质含量也大幅减小;同时受某些病虫害影响后的大豆会导致食用的人畜中毒、患病甚至死亡;此外,由于病虫害通常具有较大的传播性,往往在同一地区内,由于受到病虫害的影响不能继续对该区域进行大范围种植;另外,受病虫害影响的大豆通常也难于储存和运输,相比于健康的大豆更易腐败,且带有危险性病虫害的品种更容易受到出口的限制。Soybean, as my country's main food crop and cash crop, plays a very important role in national economic production, especially in agricultural production. Soybean diseases and insect pests are not only one of the most important agricultural disasters in my country, but also in the world. Up to now, there have been more than 120 kinds of soybean diseases and insect pests reported in the world, including fungal diseases, bacterial diseases, virus diseases, insect diseases, etc. Different types of diseases and parasitic plant diseases, different soybean diseases and insect pests have different ways, types and impacts, but different diseases and insect pests will seriously affect the growth of soybeans, resulting in a significant reduction in soybean yield or even extinction. High-level and wide-ranging pests and diseases have caused huge losses to the national economy of my country and even all countries in the world, especially the national economy of countries that mainly rely on soybean production and planting. The disease and insect pests of soybeans are reflected in the average yield of crops suffering from pests and diseases by 1 to 20%, and even in severe cases, the production will cease to exist; after suffering from diseases and insect pests, the taste and quality of crops will be reduced to varying degrees, and the starch and protein content will also be greatly reduced; at the same time Soybeans affected by certain diseases and insect pests will cause poisoning, disease and even death of the people and animals they eat; in addition, because the diseases and insect pests usually have a large spread, they are often in the same area, and due to the influence of the disease and insect pests, the area cannot continue to be carried out. In addition, pest-affected soybeans are often difficult to store and transport, more perishable than healthy soybeans, and varieties with dangerous pests and diseases are more susceptible to export restrictions.
病虫害对大豆的影响主要体现在根茎、叶片和果实上。对于大豆的根部而言,病虫害的存在会导致大豆出现发根、黑根、烂根,根瘤数量的减少,从而影响农作物根系对水分、有机质和其他营养物质的吸收,导致大豆植株萎缩、发黄、生长缓慢;对于大豆的茎而言,病虫害的存在会导致茎部出现炭疽、溃疡等局部病斑,从而影响大豆的茎对水分和养分的运输,进而导致大豆植株枯萎、腐烂、倒伏、萎蔫甚至死亡;对于大豆的叶片而言,病虫害的存在会导致作物的叶片出现褪绿、黄化、变紫等变色现象和褪色现象及圆斑、角斑、条纹等病斑现象,从而降低大豆叶片的光合效率,大幅降低农作物的产量和品质;对于大豆的果实而言,病虫害的存在会导致果实出现豆荚表面坏死病斑、大豆籽粒奇特和畸形、大豆籽粒软腐和内部干腐等现象,从而直接影响大豆的品质、降低其口感、减少大豆的产量。The effects of diseases and insect pests on soybean are mainly reflected in the rhizomes, leaves and fruits. For the roots of soybeans, the presence of pests and diseases will cause hairy roots, black roots, and rotten roots in soybeans, and the number of root nodules will decrease, which will affect the absorption of water, organic matter and other nutrients by the root system of crops, causing soybean plants to shrink and yellow. , slow growth; for soybean stems, the presence of pests and diseases will cause local lesions such as anthracnose and ulcers on the stems, which will affect the transportation of water and nutrients by soybean stems, resulting in soybean plants withering, rot, lodging and wilting. For soybean leaves, the presence of diseases and insect pests will cause discoloration and fading phenomena such as chlorosis, yellowing, and purplish discoloration, as well as disease spots such as round spots, angular spots, and stripes, thereby reducing soybean leaves. The photosynthetic efficiency greatly reduces the yield and quality of crops; for soybean fruits, the presence of diseases and insect pests will lead to the appearance of necrotic spots on the surface of the pods, strange and deformed soybean seeds, soft rot and internal dry rot of soybean seeds, etc. It directly affects the quality of soybeans, reduces its taste, and reduces the yield of soybeans.
因此,对于病虫害的识别和病虫害类型的高效、准确地鉴定,对于确定大豆病虫害的病原、病程和由病虫害引起的产量评估、经济损失评估,以及针对病虫害类型所采取的后续防治措施,都具有重要的实际意义和指导价值。而现有的大豆病虫害识别和鉴定方法,往往都是采用柯赫式法则通过对获取染病的植株部分进行目视的田间观察和症状的初步诊断,再结合对染病植株的实验室显微测量和理化参数的分析以确定病虫害类型及产生原因,这种测量往往需要专业的从业人员进行,且判别过程具有一定的滞后性,在测量过程中还可能存在误诊或漏诊等现象,从而贻误最好的治疗时期,对受病虫害影响的大豆造成更大的损失。Therefore, the identification of diseases and insect pests and the efficient and accurate identification of the types of diseases and insect pests are of great importance for determining the pathogen, disease course, yield assessment and economic loss caused by diseases and insect pests, as well as follow-up control measures for the types of diseases and insect pests. practical significance and guiding value. However, the existing methods for identification and identification of soybean diseases and insect pests are often based on the Koch-type rule through visual field observation and preliminary diagnosis of symptoms on the infected plant parts, combined with laboratory microscopic measurements and laboratory microscopic measurements of infected plants. Analysis of physical and chemical parameters to determine the types and causes of pests and diseases, this kind of measurement often requires professional practitioners, and the discrimination process has a certain lag, and there may be misdiagnosis or missed diagnosis in the measurement process, thus delaying the best results. During the treatment period, greater losses are caused to soybeans affected by pests and diseases.
发明内容SUMMARY OF THE INVENTION
基于以上不足之处,本发明提供一种新的大豆病虫害的有效鉴定方法,解决了现有技术基于田间目视观察和实验室分析等现有的对大豆病虫害识别和鉴定方面的滞后性和不准确的问题。Based on the above deficiencies, the present invention provides a new effective identification method for soybean diseases and insect pests, which solves the hysteresis and inconsistency in the existing identification and identification of soybean diseases and insect pests in the prior art based on field visual observation and laboratory analysis. exact question.
本发明所采用的技术如下:一种大豆病虫害的有效鉴定方法,步骤如下:The technology adopted in the present invention is as follows: an effective identification method for soybean diseases and insect pests, the steps are as follows:
步骤1、大豆器官图像的标准化获取Step 1. Standardized acquisition of soybean organ images
对不同生长期内受病虫害影响的大豆进行标准化图像的获取,涉及对受病虫害影响的大豆植株的叶片部分、根系部分和籽粒部分的标准化图像的获取;Acquiring standardized images of soybeans affected by diseases and insect pests in different growth periods involves the acquisition of standardized images of leaf, root and grain parts of soybean plants affected by diseases and insect pests;
步骤2、大豆器官图像的标准化处理Step 2. Standardization of soybean organ images
对不同生长期内大豆器官图像获取的标准化照片进行几何畸变校正,去除拍摄过程对照片带来的几何形变,并针对不同的器官图像选择统一的尺寸进行裁剪操作,包括大豆植株的叶片标准化图像、根系标准化图像、籽粒标准化图像;Geometric distortion correction is carried out on the standardized photos obtained from soybean organ images in different growth periods to remove the geometric distortions brought by the photographing process, and a uniform size is selected for different organ images for cropping operations, including standardized images of soybean plants, leaves, Normalized image of root system, normalized image of grain;
步骤3、大豆病虫害鉴定参数的获取Step 3. Acquisition of soybean disease and insect pest identification parameters
本步骤涉及对不同生长期内,受病虫害影响的大豆器官特征参数的提取,包括对对大豆叶片特征参数的获取、大豆根系特征参数的获取、大豆籽粒特征参数的获取;This step involves the extraction of characteristic parameters of soybean organs affected by diseases and insect pests in different growth periods, including the acquisition of characteristic parameters of soybean leaves, the acquisition of characteristic parameters of soybean roots, and the acquisition of characteristic parameters of soybean grains;
步骤4、大豆病虫害种类数据库的建立Step 4. Establishment of database of soybean pests and diseases
利用已知的不同种类大豆病虫害影响下的大豆植株特征图谱建立数据库,在数据库中对每种病虫害影响下大豆植株在不同生长期下,对步骤3中所涉及的各类特征参数进行提取和保存,根据病虫害种类的差异及其对植株影响的差异对不同的特征参数赋予不同的权重,根据特征参数的提取结果建立病虫害种类的判别模型;A database is established using the known characteristic maps of soybean plants under the influence of different types of soybean diseases and insect pests. In the database, the various characteristic parameters involved in step 3 are extracted and saved for each type of soybean plant under the influence of diseases and insect pests under different growth periods. , assign different weights to different characteristic parameters according to the difference of the types of diseases and insect pests and their effects on plants, and establish the discriminant model of the types of diseases and insect pests according to the extraction results of the characteristic parameters;
步骤5、大豆病虫害的识别和鉴定Step 5. Identification and identification of soybean diseases and insect pests
利用数据库中不同生长期下,受到不同种类大豆病虫害影响的植株特征参数以及数据库中的病虫害种类的判别模型对实际受影响的大豆病虫害种类的识别和鉴定,根据待鉴定的大豆植株的特征参数提取结果与数据库中已有的每种病虫害下植株的特征参数进行对比计算,最后利用待鉴定大豆植株的特征参数的判别模型计算结果与每种病虫害下植株的判别模型计算结果进行对比计算和匹配,从而实现对大豆病虫害类型的识别和鉴定。Using the characteristic parameters of plants affected by different types of soybean diseases and insect pests in the database under different growth periods and the discriminant model of the types of diseases and insect pests in the database to identify and identify the types of soybean diseases and insect pests that are actually affected, extract the characteristic parameters of soybean plants to be identified. The results are compared and calculated with the characteristic parameters of the plants under each disease and insect pest existing in the database. Finally, the calculation results of the discriminant model of the characteristic parameters of the soybean plants to be identified are compared with the calculation results of the discriminant model of the plants under each disease and insect pest to compare, calculate and match. So as to realize the identification and identification of soybean diseases and insect pests.
本发明还具有如下技术特征:The present invention also has the following technical features:
1、如上所述的步骤3:对受病虫害影响的大豆植株在开花期、结荚期、鼓粒期及成熟期的4个大豆关键生长期获取的叶片标准化图像进行进一步处理,针对每个带鉴定病虫害类型的大豆叶片标准化图像提取相应的特征参数,对于某一特定生长期内受到待鉴定病虫害影响大豆植株的叶片,为去除样本间特征参数提取结果的异质性,选择20个植株,每个植株随机选择一个叶片的标准化彩色图像,提取叶片标准化彩色图像的红、绿、蓝三个分量的图像矩阵R(x,y),G(x,y),B(x,y),然后对三个分量按公式F(x,y)=R(x,y)/3+G(x,y)/3+B(x,y)/3将叶片的标准化彩色图像转换为灰度图,对灰度图进行直方图统计,找到大豆叶片绿色部分和受病虫害影响的变色部分两个灰度部分峰值间的频数最低点对应的灰度值t1,将灰度大于阈值t1的所有像元亮度值设为0,灰度小于阈值t1的多有像元亮度值设为1,实现将灰度图转化为二值图;根据叶片图像宽度对应的像元个数值w1和叶片图像对应植株的测量结果w1’厘米,计算每个像元的实际边长z1=w1’/w1厘米,每个像元的实际面积则为z12平方厘米,根据像元的实际边长z1厘米、像元的实际面积z12平方厘米,以及每个斑点不同参数提取结果对应的像元数,计算该特定生长期内待鉴定的受病虫害植株叶片特征参数的均值,包括待鉴定病虫害植株叶片斑点面积比值a1、待鉴定病虫害植株叶片斑点个数值a2,待鉴定病虫害植株叶片斑点面积值a3,待鉴定病虫害植株叶片斑点周长值a4;再次对受病虫害影响植株灰度图进行直方图统计从而进一步提取受病虫害影响叶片的整体形态特征,找到大豆叶片和背景色两个灰度部分峰值间的灰度级频数最低点对应的灰度值t1’,将灰度大于阈值t1’的所有像元亮度值设为0,将灰度小于t1’的所有像元亮度值设为1,实现对大豆整体叶片的二值化处理,根据叶片二值图提取叶片长轴对应的像元数和短轴对应的像元数计算该特定生长期待鉴定病虫害植株叶片长宽比a5和该特定生长期待鉴定病虫害植株叶片椭圆率a6;1. Step 3 as described above: further process the standardized images of leaves obtained during the four key growth stages of soybean plants affected by diseases and insect pests in the flowering stage, pod setting stage, drumming stage and maturity stage. Corresponding characteristic parameters were extracted from the standardized images of soybean leaves that identified the types of diseases and insect pests. For the leaves of soybean plants affected by the diseases and insect pests to be identified in a specific growth period, in order to remove the heterogeneity of the extraction results of characteristic parameters among samples, 20 plants were selected, each Each plant randomly selects a standardized color image of a leaf, extracts the image matrices R(x,y), G(x,y), B(x,y) of the red, green, and blue components of the standardized color image of the leaf, and then Convert the normalized color image of the leaf to a grayscale image according to the formula F(x,y)=R(x,y)/3+G(x,y)/3+B(x,y)/3 for the three components , perform histogram statistics on the grayscale map, find the grayscale value t1 corresponding to the lowest frequency point between the two grayscale peaks of the green part of soybean leaves and the discolored part affected by diseases and insect pests, and set all pixels whose grayscale is greater than the threshold t1. The brightness value is set to 0, and the brightness value of many pixels whose grayscale is less than the threshold t1 is set to 1, so as to convert the grayscale image into a binary image; The measurement result w1' cm, calculate the actual side length of each pixel z1=w1'/w1 cm, the actual area of each pixel is z1 2 square cm, according to the actual side length of the pixel z1 cm, the size of the pixel The actual area z1 is 2 square centimeters, and the number of pixels corresponding to the extraction results of different parameters of each spot, calculate the average value of the leaf characteristic parameters of the plant affected by the disease and insect pest to be identified in this specific growth period, including the ratio of the leaf spot area of the plant to be identified, a1, The number of spots on the leaves of the plants to be identified is a2, the area value of the spots on the leaves of the plants to be identified is a3, and the perimeter value of the spots on the leaves of the plants to be identified is a4; once again, perform histogram statistics on the grayscale images of the plants affected by diseases and insect pests to further extract those affected by diseases and insect pests The overall morphological characteristics of leaves, find the gray value t1' corresponding to the lowest gray level frequency between the two gray peaks of soybean leaf and background color, and set the brightness value of all pixels whose gray value is greater than the threshold t1' to 0 , set the brightness value of all pixels whose gray level is less than t1' to 1, realize the binarization processing of the whole soybean leaf, and extract the number of pixels corresponding to the long axis and the number of pixels corresponding to the short axis of the leaf according to the leaf binary map Calculate the leaf length-to-width ratio a5 of the specific growth expectation identification plant with pests and diseases and the leaf ellipticity a6 of the specific growth expectation identification plant with pests and diseases;
对于某一特定生长期内受病虫害影响大豆植株的根系,为去除样本间指标参数提取结果的异质性,选择20个植株的根系标准化图像,提取叶片彩色图像的红、绿、蓝三个分量的图像矩阵R(x,y),G(x,y),B(x,y),对三个分量按公式F(x,y)=R(x,y)/3+G(x,y)/3+B(x,y)/3将彩色叶片图像转换为灰度图,对根系灰度图像进行直方图统计,找到大豆根系和背景色两个灰度部分峰值间的灰度级频数最低点对应的灰度值t2,将灰度大于阈值t2的所有像元亮度值设为0,将灰度小于t2的所有像元亮度值设为1,实现将根系灰度图转换为二值图,根据根系图像宽度对应的像元个数值w2和根系图像对应植株的测量结果w2’厘米,计算每个像元的实际边长z2=w2’/w2厘米,每个像元的实际面积则为z22平方厘米,根据像元的实际边长z2厘米、像元的实际面积z22平方厘米,以及受病虫害影响的植株根系不同参数提取结果对应的像元数,计算该特定生长期待鉴定病虫害植株根系面积值b1,待鉴定病虫害植株根系根瘤个数值b2,根据灰度图计算该特定生长期内待鉴定的受病虫害植株根系纹理特征参数的均值,根系纹理特征参数的具体计算过程如下,根据公式1计算待鉴定病虫害植株根系对比度纹理特征值b3,根据公式2待鉴定病虫害植株根系一致性纹理特征值b4,根据公式3待鉴定病虫害植株根系熵值纹理特征值b5,根据公式4待鉴定病虫害植株根系能量纹理特征值b6,For the roots of soybean plants affected by diseases and insect pests in a specific growth period, in order to remove the heterogeneity of the extraction results of index parameters between samples, 20 standardized root images of plants were selected, and the red, green and blue components of the leaf color image were extracted. The image matrix R(x,y), G(x,y), B(x,y), for the three components according to the formula F(x,y)=R(x,y)/3+G(x, y)/3+B(x,y)/3 Convert the color leaf image to grayscale, perform histogram statistics on the root grayscale image, and find the grayscale between the two grayscale peaks of soybean root and background color The grayscale value t2 corresponding to the lowest frequency point, set the brightness value of all pixels whose grayscale is greater than the threshold t2 to 0, and set the brightness value of all pixels whose grayscale is less than t2 to 1, so as to convert the root grayscale image into two Value map, according to the pixel value w2 corresponding to the width of the root image and the measurement result of the plant corresponding to the root image w2' cm, calculate the actual side length of each pixel z2=w2'/w2 cm, the actual area of each pixel Then it is z2 2 square centimeters. According to the actual side length of the pixel z2 cm, the actual area of the pixel z2 2 square centimeters, and the number of pixels corresponding to the extraction results of different parameters of the root system of plants affected by pests and diseases, the identification of the specific growth expectation is calculated. The root area value of the plant affected by diseases and insect pests is b1, and the number of root nodules of the plant to be identified is b2. Calculate the average value of the root texture characteristic parameters of the plant affected by diseases and insect pests to be identified in the specific growth period according to the grayscale map. The specific calculation process of the root texture characteristic parameters is as follows: Calculate the root contrast texture feature value b3 of the plant with plant diseases and insect pests to be identified according to formula 1, the texture feature value b4 of the root system consistency of the plants to be identified according to formula 2, the texture feature value b5 of the root system entropy value of the plants to be identified according to formula 3, and the texture feature value b5 of the plant root system to be identified according to formula 4 The eigenvalue b6 of root energy texture of plants with pests and diseases,
其中,p(i,j)表示由根系灰度图像计算的灰度共生矩阵在第i行,第j列位置上的值,表示根系灰度图像中灰度级为i的像元和灰度级为j的像元在固定的方向和像元间隔下同时出现的概率,n为根系灰度图像中i和j两个灰度级的差异,Ng表示根系灰度图像提取的灰度共生矩阵的级数;Among them, p(i,j) represents the value of the grayscale co-occurrence matrix calculated from the root grayscale image at the i-th row and j-th column position, and represents the pixel and gray level of the gray level i in the root grayscale image. The probability that the pixels with level j appear at the same time in a fixed direction and pixel interval, n is the difference between the two gray levels of i and j in the root gray image, and Ng represents the gray co-occurrence matrix extracted from the root gray image series;
对于某一特定生长期内受病虫害影响大豆植株的籽粒,为去除样本间指标参数提取结果的异质性,选择20个植株,提取每个植株籽粒彩色籽粒图像的红、绿、蓝三个分量的图像矩阵R(x,y),G(x,y),B(x,y),对三个分量按公式F(x,y)=R(x,y)/3+G(x,y)/3+B(x,y)/3将彩色籽粒图像转换为灰度图,对灰度图进行直方图统计找到大豆籽粒部分和试验台背景部分两个灰度部分峰值间的频数最低点对应的灰度值t3,将灰度大于阈值t3的所有像元亮度值设为1,灰度小于阈值t3的多有像元亮度值设为0,实现将籽粒灰度图转化为二值图,根据受病虫害影响的籽粒图像宽度对应的像元个数值w3和根系图像对应植株的测量结果w3’厘米,计算每个像元的实际边长z3=w3’/w3厘米,每个像元的实际面积则为z32平方厘米,根据像元的实际边长z3厘米、像元的实际面积z32平方厘米,以及受病虫害影响的植株籽粒不同参数提取结果对应的像元数,计算该特定生长期内待鉴定的受病虫害植株籽粒特征参数的均值,包括待鉴定病虫害植株籽粒个数值c1,待鉴定病虫害植株籽粒面积值c2,待鉴定病虫害植株籽粒周长值c3,待鉴定病虫害植株籽粒长宽比值c4,待鉴定病虫害植株籽粒曲率值c5,待鉴定病虫害植株籽粒椭圆率值c6。For the grains of soybean plants affected by diseases and insect pests in a specific growth period, in order to remove the heterogeneity of the extraction results of index parameters among samples, 20 plants were selected, and the red, green and blue components of the color grain image of each plant's grain were extracted. The image matrix R(x,y), G(x,y), B(x,y), for the three components according to the formula F(x,y)=R(x,y)/3+G(x, y)/3+B(x,y)/3 Convert the color grain image into a grayscale image, and perform histogram statistics on the grayscale image to find the lowest frequency between the two grayscale peaks in the soybean kernel part and the background part of the test bench The gray value t3 corresponding to the point, the brightness value of all pixels whose gray value is greater than the threshold value t3 is set to 1, and the brightness value of many pixels whose gray value is less than the threshold value t3 is set to 0, so as to convert the grayscale image of the grain into two values. Figure, calculate the actual side length of each pixel z3=w3'/w3 cm according to the pixel value w3 corresponding to the width of the grain image affected by the disease and insect pests and the measurement result of the root image corresponding to the plant w3' cm, each pixel The actual area is z3 2 square centimeters, according to the actual side length of the pixel z3 centimeters, the actual area of the pixel z3 2 square centimeters, and the number of pixels corresponding to the extraction results of different parameters of plant grains affected by pests and diseases, calculate the specific The average value of the grain characteristic parameters of the plants affected by diseases and insect pests to be identified during the growth period, including the number of grains of the plants to be identified c1, the value of the grain area of the plants to be identified c2, the value of the grain circumference of the plants to be identified c3, and the length of the grains of the plants to be identified. The width ratio is c4, the grain curvature value of the plants to be identified with pests and diseases is c5, and the grain ellipticity value of the plants to be identified with pests and diseases is c6.
2、如上所述的步骤4:利用开花期、结荚期、鼓粒期及成熟期的不同生长期内受到所有已知的病虫害类型影响的大豆植株图像建立特征参数数据库,对某一特定生长期内已知病虫害影响下大豆叶片特征参数的提取和建库结合步骤3的方法将该特定生长期下受到已知病虫害影响下的彩色叶片图像转换为灰度图,利用步骤3的方法对灰度图进行直方图统计,找到叶片绿色部分和受病虫害影响的变色部分两个灰度部分峰值间的频数最低点对应的灰度值T1,将灰度大于阈值T1的所有像元亮度值设为0,将灰度小于阈值T1的多有像元亮度值设为1,实现将灰度图转化为二值图,根据叶片图像宽度对应的像元个数值W1和叶片图像对应植株的测量结果W1’厘米,计算每个像元的实际边长Z1=W1’/W1厘米,每个像元的实际面积则为Z12平方厘米,根据像元的实际边长Z1厘米、像元的实际面积Z12平方厘米,以及每个斑点不同参数提取结果对应的像元数,计算该特定生长期内已知病虫害类型的植株叶片特征参数的均值,包括已知病虫害植株叶片斑点面积比值A1、已知病虫害植株叶片斑点个数值A2,已知病虫害植株叶片斑点面积值A3,已知病虫害植株叶片斑点周长值A4,为进一步提取受到该已知类型病虫害影响叶片的整体形态特征,再次对受该已知病虫害影响植株灰度图进行直方图统计,找到大豆叶片和背景色两个灰度部分峰值间的灰度级频数最低点对应的灰度值T1’,将灰度大于阈值T1’的所有像元亮度值设为0,将灰度小于T1’的所有像元亮度值设为1,实现对大豆整体叶片的二值化处理,根据叶片二值图提取叶片长轴对应的像元数和短轴对应的像元数计算该特定生长期已知病虫害植株叶片斑点长宽比值A5和已知病虫害植株斑点椭圆率值A6;2. Step 4 as described above: use the images of soybean plants affected by all known types of diseases and insect pests in different growth stages of flowering, pod setting, drumming and maturity to establish a characteristic parameter database, and for a specific growth period. The extraction and database construction of soybean leaf characteristic parameters under the influence of known diseases and insect pests during the period were combined with the method of step 3 to convert the color leaf image under the influence of known diseases and insect pests under the specific growth period into a grayscale image, and the method of step 3 was used to convert the gray scale image. The degree map is used for histogram statistics, and the gray value T1 corresponding to the lowest frequency point between the two gray value peaks of the green part of the leaf and the discolored part affected by diseases and insect pests is found, and the brightness value of all pixels whose gray value is greater than the threshold T1 is set as 0, set the brightness value of many pixels whose grayscale is less than the threshold T1 to 1, to convert the grayscale image into a binary image, according to the pixel value W1 corresponding to the width of the leaf image and the measurement result W1 of the plant corresponding to the leaf image 'cm, calculate the actual side length of each pixel Z1=W1'/W1 cm, the actual area of each pixel is Z1 2 square cm, according to the actual side length of the pixel Z1 cm, the actual area of the pixel Z1 2 square centimeters, and the number of pixels corresponding to the extraction results of different parameters for each spot, calculate the average value of the plant leaf characteristic parameters of the known disease and insect pest types in this specific growth period, including the known disease and insect pest plant leaf spot area ratio A1, known disease and insect pests The number of leaf spots of plants is A2, the area value of leaf spots of plants with known diseases and insect pests is A3, and the perimeter value of leaf spots of plants with known diseases and insect pests is A4. The grayscale map of plants affected by pests and diseases is used for histogram statistics, and the grayscale value T1' corresponding to the lowest point of grayscale frequency between the two grayscale peaks of soybean leaf and background color is found, and all pixels whose grayscale is greater than the threshold T1' The brightness value is set to 0, and the brightness value of all pixels whose grayscale is less than T1' is set to 1, so as to realize the binarization processing of the whole soybean leaf, and extract the number of pixels corresponding to the long axis and the short axis of the leaf according to the leaf binary map. The corresponding pixel number is used to calculate the leaf spot aspect ratio value A5 of the plant with known plant diseases and insect pests and the ellipticity value of the plant spot with known plant diseases and insect pests in the specific growth period A6;
对某一特定生长期内已知病虫害影响下大豆根系特征参数的提取和建库结合步骤3的方法将该特定生长期下受到已知病虫害影响下的彩色根系图像转换为灰度图,利用步骤3的方法对灰度图进行直方图统计,找到大豆根系和背景色两个灰度部分峰值间的灰度级频数最低点对应的灰度值T2,将灰度大于阈值T2的所有像元亮度值设为0,将灰度小于T2的所有像元亮度值设为1,实现将根系灰度图转换为二值图,根据根系图像宽度对应的像元个数值W2和根系图像对应植株的测量结果W2’厘米,计算每个像元的实际边长Z2=W2’/W2厘米,每个像元的实际面积则为Z22平方厘米,根据像元的实际边长Z2厘米、像元的实际面积Z22平方厘米,以及受病虫害影响的植株根系不同参数提取结果对应的像元数,计算所有受到已知病虫害影响植株根系特征参数的均值,包括已知病虫害植株根系面积值B1,已知病虫害植株根系根瘤个数值B2,根据步骤三中的公式1计算已知病虫害植株根系对比度纹理特征值B3,根据公式2计算已知病虫害植株根系一致性纹理特征值B4,根据公式3计算已知病虫害植株根系熵值纹理特征值B5,根据公式4计算已知病虫害植株根系能量纹理特征值B6;The extraction and database establishment of soybean root characteristic parameters under the influence of known diseases and insect pests in a specific growth period are combined with the method of step 3 to convert the color root system images under the influence of known diseases and insect pests in the specific growth period into grayscale images, and the steps The method of 3 performs histogram statistics on the grayscale image, and finds the grayscale value T2 corresponding to the lowest point of the grayscale frequency between the peaks of the two grayscale parts of the soybean root system and the background color. The value is set to 0, the brightness value of all pixels whose grayscale is less than T2 is set to 1, and the root grayscale image is converted into a binary image. According to the pixel value W2 corresponding to the width of the root image and the measurement of the plant corresponding to the root image Result W2' cm, calculate the actual side length of each pixel Z2=W2'/W2 cm, the actual area of each pixel is Z2 2 square cm, according to the actual side length of the pixel Z2 cm, the actual size of the pixel The area Z2 is 2 square centimeters, and the number of pixels corresponding to the extraction results of different parameters of plant roots affected by plant diseases and insect pests. Calculate the mean value of root characteristic parameters of all plants affected by known plant diseases and insect pests, including known plant root system area value B1, known plant diseases and insect pests The number of root nodules in the plant root system is B2. According to the formula 1 in step 3, the root contrast texture characteristic value B3 of the known plant with plant diseases and insect pests is calculated, the root consistency texture characteristic value of the plant with known plant diseases and insect pests is calculated according to the formula 2, and the texture characteristic value of the known plant root system is calculated according to the formula 3. Root entropy value texture feature value B5, calculate root energy texture feature value B6 of known plant diseases and insect pests according to formula 4;
对于某一特定生长期内受到已知病虫害影响大豆植株的籽粒,结合步骤3中的方法将该特定生长期下受到已知病虫害影响下的彩色籽粒图像转换为灰度图,对灰度图进行直方图统计找到大豆籽粒部分和试验台背景部分两个灰度部分峰值间的频数最低点对应的灰度值T3,将灰度大于阈值T3的所有像元亮度值设为1,灰度小于阈值T3的多有像元亮度值设为0,实现将籽粒灰度图转化为二值图,根据受病虫害影响的籽粒图像宽度对应的像元个数值W3和根系图像对应植株的测量结果W3’厘米,计算每个像元的实际边长Z3=W3’/W3厘米,每个像元的实际面积则为Z32平方厘米,根据像元的实际边长Z3厘米、像元的实际面积Z32平方厘米,以及受病虫害影响的植株籽粒不同参数提取结果对应的像元数,计算该特定生长期内已知受病虫害影响的植株籽粒特征参数的均值,包括已知病虫害植株籽粒个数值C1,已知病虫害植株籽粒面积值C2,已知病虫害植株籽粒周长值C3,已知病虫害植株籽粒长宽比值C4,已知病虫害植株籽粒曲率值C5,已知病虫害植株籽粒椭圆率值C6,根据不同类型病虫害对大豆植株器官影响的差异,确定每种已知的大豆病虫害类型叶片特征参数的权重K1,确定每种已知的大豆病虫害类型根系特征参数的权重K2,确定每种已知的大豆病虫害类型籽粒特征参数的权重K3,同时满足权重参数的归一化特征,即K1+K2+K3=1,根据已经提取的待鉴定的大豆病虫害类型的不同器官图像特征参数提取结果和数据库中已知的病虫害种类的不同器官图像特征参数提取结果,建立病虫害判别模型,判别模型如公式5所示:For the grains of soybean plants affected by known diseases and insect pests in a specific growth period, the color grain images under the influence of known diseases and insect pests in the specific growth period are converted into grayscale images in combination with the method in step 3, and the grayscale images are analyzed. The histogram statistics find the gray value T3 corresponding to the lowest frequency point between the two gray value peaks of the soybean grain part and the background part of the test bench, and set the brightness value of all pixels whose gray value is greater than the threshold T3 to 1, and the gray value is less than the threshold value. The brightness value of most pixels in T3 is set to 0, and the grayscale image of the grain is converted into a binary image. , calculate the actual side length of each pixel Z3=W3'/W3 cm, the actual area of each pixel is Z3 2 square cm, according to the actual side length of the pixel Z3 cm, the actual area of the pixel Z3 2 square centimeters, and the number of pixels corresponding to the extraction results of different parameters of plant grains affected by diseases and insect pests, calculate the average value of characteristic parameters of plant grains known to be affected by plant diseases and insect pests in this specific growth period, including the number of known plant grains C1, known The grain area value of plants with known diseases and insect pests is C2, the perimeter value of grains of plants with known diseases and insect pests is C3, the ratio of grain length to width of plants with known diseases and insect pests is C4, the curvature value of grains of plants with known diseases and insect pests is C5, and the ellipticity value of grains of plants with known diseases and insect pests is C6. To determine the difference in the effects on soybean plant organs, determine the weight K 1 of the leaf characteristic parameters of each known soybean pests and diseases, determine the weight K2 of the root characteristic parameters of each known soybean pests and diseases, and determine the weights of each known soybean pests and diseases. The weight K 3 of the characteristic parameter of the type grain, and at the same time satisfy the normalized characteristic of the weight parameter, that is, K 1 +K 2 +K 3 =1, according to the extraction results and Based on the extraction results of different organ image feature parameters of known pest species in the database, a pest and disease discrimination model is established. The discriminant model is shown in Equation 5:
4、如上所述的步骤5:利用步骤4数据库中不同生长期下,受到不同种类已知大豆病虫害影响的植株特征参数以及数据库中已知的病虫害种类判别模型实现对受影响的大豆病虫害种类的识别和鉴定,对某一个待鉴定的受病虫害植株叶片特征参数提取结果a1、a2、a3、a4、a5、a6;待鉴定的受病虫害植株根系特征参数提取结果b1、b2、b3、b4、b5、b6;待鉴定的受病虫害植株籽粒特征参数提取结果c1、c2、c3、c4、c5、c6,将每一个待鉴定的受病虫害影响大豆的特征参数提取结果和该生长期下每一种已知受病虫害影响大豆的特征参数提取结果带入到步骤4中的判别模型中,根据判断模型计算待鉴定病虫害类型与所有已知病虫害类型的模型计算结果P值的大小,P值计算结果越小,证明待鉴定病虫害类型越接近该P值计算结果所对应的那一个已知的病虫害类型,从而实现对大豆病虫害类型的识别和鉴定。4. Step 5 as described above: use the characteristic parameters of plants affected by different types of known soybean diseases and insect pests in the database in step 4 under different growth periods and the known disease and insect pest species discrimination models in the database to realize the identification of affected soybean diseases and insect pests. Identification and identification, extraction results a1, a2, a3, a4, a5, a6 of leaf characteristic parameters of a plant affected by diseases and insect pests to be identified; extraction results b1, b2, b3, b4, b5 of root characteristic parameters of plants affected by diseases and insect pests to be identified , b6; the extraction results c1, c2, c3, c4, c5, and c6 of the grain characteristic parameters of the plants affected by diseases and insect pests to be identified, the extraction results of the characteristic parameters of each soybean affected by diseases and insect pests to be identified and the The extraction results of the characteristic parameters of soybeans known to be affected by diseases and insect pests are brought into the discriminant model in step 4, and the P value of the model calculation results of the types of diseases and insect pests to be identified and all known types of diseases and insect pests is calculated according to the judgment model, and the smaller the calculation result of the P value is , which proves that the pest type to be identified is closer to the known pest and disease type corresponding to the calculation result of the P value, so as to realize the identification and identification of soybean pest and disease types.
本发明所具有的有益效果及优点:本发明涉及的方法与实验室的分析鉴定相比,鉴定成本少、造价低廉,损耗小,对植株的破坏也较小,且具有实时性和准实时性的特点,鉴定人员仅通过图像就可以快速的获取待鉴定病虫害的大豆植株特征参数,并根据提取结果实现对病虫害的鉴定。此外,本发明对鉴别人员的专业知识要求很低,鉴别人员只需要经过简单培训就可以利用本发明涉及的方法,就可以实现利用数据库中的各种已知病虫害类型的大量数据对待鉴定的大豆病虫害类型进行识别和鉴定,而整个鉴定过程不需要鉴别人员自身对病虫害知识的全面掌握。此外,与传统的田间通过目测进行大豆病虫害类型的鉴定方法相比,本发明涉及的方法鉴定速度快,由于待鉴定植株参数的提取过程自动化程度高,本发明还可以通过编程手段同时批量处理大量的待鉴别病虫害类型的大豆植株图像并提取待鉴别病虫害类型的大豆植株参数,实现对大豆植株病虫害类型的快速、高效的批量识别。The beneficial effects and advantages of the present invention: compared with the analysis and identification in the laboratory, the method involved in the present invention has less identification cost, low cost, less loss, less damage to plants, and has real-time and quasi-real-time performance. The identification personnel can quickly obtain the characteristic parameters of soybean plants to be identified with pests and diseases only through images, and realize the identification of pests and diseases according to the extraction results. In addition, the present invention has very low requirements for the professional knowledge of the identification personnel, and the identification personnel can use the method involved in the present invention only through simple training, and can realize the soybeans to be identified by utilizing a large amount of data of various known types of diseases and insect pests in the database. The types of pests and diseases are identified and identified, and the entire identification process does not require the identification personnel to have a comprehensive grasp of the knowledge of pests and diseases. In addition, compared with the traditional method for identifying the types of soybean diseases and insect pests through visual inspection in the field, the method involved in the present invention has a fast identification speed. Due to the high degree of automation in the extraction process of the plant parameters to be identified, the present invention can also batch process a large number of plants simultaneously by programming means. The soybean plant image of the type of diseases and insect pests to be identified is extracted, and the parameters of the soybean plant of the type of diseases and insect pests to be identified are extracted, so as to realize the rapid and efficient batch identification of the types of diseases and insect pests of soybean plants.
附图说明Description of drawings
图1为受病虫害影响植株叶片灰度图;Figure 1 is a grayscale image of leaves of plants affected by diseases and insect pests;
图2为受病虫害影响植株叶片内部二值图;Fig. 2 is a binary map of the interior of leaves of plants affected by diseases and insect pests;
图3受病虫害影响植株叶片整体二值图;Figure 3. The overall binary image of leaves of plants affected by diseases and insect pests;
图4受病虫害影响植株根系灰度图;Fig. 4 Grayscale image of roots of plants affected by pests and diseases;
图5受病虫害影响植株根系二值图;Fig. 5 Binary map of root system of plants affected by diseases and insect pests;
图6受病虫害影响植株籽粒灰度图;Figure 6 Grayscale image of grains of plants affected by diseases and insect pests;
图7受病虫害影响植株籽粒二值图;Fig. 7 Binary image of grains of plants affected by diseases and insect pests;
具体实施方式Detailed ways
下面根据说明书附图举例对本发明做进一步的说明:The present invention will be further described below according to the accompanying drawings of the description:
实施例1Example 1
步骤1:获取大豆在开花期、结荚期、鼓粒期及成熟期等不同生长期受到病虫害影响的典型大豆植株器官标准化图像,Step 1: Obtain the standardized images of typical soybean plant organs affected by diseases and insect pests in different growth stages such as flowering stage, pod setting stage, drumming stage and maturity stage of soybean,
对于某一生长期受病虫害影响的大豆叶片图像,将数码相机安装在固定支架上,将受病虫害影响的叶片放置在试验台上,调整支架位置使相机镜头垂直于试验台,且相机镜头距离大豆叶片20cm的高度,利用光度计确定和调整拍照环境,对整片受病虫害影响的叶片进行拍照获取叶片的标准化图像并拍摄棋盘格网定标板图像;For images of soybean leaves affected by diseases and insect pests in a certain growth period, install a digital camera on a fixed bracket, place the leaves affected by diseases and insect pests on the test bench, adjust the position of the bracket so that the camera lens is perpendicular to the test bench, and the camera lens is far from the soybean leaves At a height of 20cm, use a photometer to determine and adjust the photographing environment, take pictures of the entire leaf affected by diseases and insect pests to obtain a standardized image of the leaf and take a checkerboard grid calibration plate image;
对于某一生长期受病虫害影响的大豆根系图像,将数码相机安装在固定支架上,调整相机的位置,使相机镜头与地表垂直,将大豆植株平铺在地表,保证相机镜头距离地表1m高度,利用光度计确定和调整拍照环境,对大根系进行拍照,用于获取受病虫害影响的大豆根系标准化图像并拍摄棋盘格网定标板图像;For soybean root images affected by diseases and insect pests in a certain growth period, install a digital camera on a fixed bracket, adjust the position of the camera so that the camera lens is perpendicular to the ground, lay the soybean plants on the ground, and ensure that the camera lens is at a height of 1m from the ground. The photometer determines and adjusts the photographing environment, and takes pictures of the large root system, which is used to obtain the standardized image of the soybean root system affected by diseases and insect pests and to take the image of the checkerboard grid calibration plate;
对于某一生长期受病虫害影响的大豆籽粒图像,将受病虫害影响的豆荚放置在试验台上,调整支架位置使相机镜头垂直于试验台,然后将整个植株的大豆籽粒取出并分散放置在试验台上,调整支架位置使相机镜头垂直于试验台,且相机镜头距离大豆籽粒20cm的高度,利用光度计确定和调整拍照环境,拍摄受病虫害影响大豆颗粒的标准化图像并拍摄棋盘格网定标板图像。For an image of soybean grains affected by diseases and insect pests in a certain growth period, place the pods affected by diseases and insect pests on the test bench, adjust the position of the bracket so that the camera lens is perpendicular to the test bench, and then take out the soybean grains of the whole plant and place them on the test bench. , adjust the position of the bracket so that the camera lens is perpendicular to the test bench, and the camera lens is 20 cm from the height of the soybean grains, and the photometer is used to determine and adjust the photographing environment, take the standardized image of soybean particles affected by diseases and insect pests, and take the checkerboard grid calibration plate image.
步骤2:对步骤1中受病虫害影响各类型大豆器官图像提取结果标准化处理,对于不同类型的大豆器官图像,由于标准化图像与相应的定标板图像拍摄条件相同,因此不同类型大豆器官的标准化图像与对应的黑白棋盘定标板图像的几何畸变也是相同的,本实施例对每种类型大豆器官相对应的黑白棋盘格网定标板采用多项式校正方法进行几何畸变校正,记录不同器官类型图像下的多项式公式参数,然后对该类型下大豆器官的标准化图像进行几何校正,并选择合适的尺寸对该类型下的大豆器官图像进行裁剪,从而实现对受病虫害影响的各类型大豆器官图像的标准化预处理操作。Step 2: Standardize the image extraction results of various types of soybean organs affected by diseases and insect pests in step 1. For different types of soybean organ images, since the normalized images and the corresponding calibration plate images have the same shooting conditions, the standardized images of different types of soybean organs are The geometric distortion of the corresponding black-and-white checkerboard calibration plate image is also the same. In this embodiment, the black-and-white checkerboard grid calibration plate corresponding to each type of soybean organ adopts the polynomial correction method to correct the geometric distortion, and records the images of different organ types. Then, the standardized images of soybean organs under this type are geometrically corrected, and the appropriate size is selected to crop the soybean organ images under this type, so as to realize the standardized prediction of various types of soybean organ images affected by diseases and insect pests. Handling operations.
步骤3:对受病虫害影响的大豆植株在开花期、结荚期、鼓粒期及成熟期的4个大豆关键生长期获取的叶片标准化图像进行进一步处理,针对每个带鉴定病虫害类型的大豆叶片标准化图像提取相应的特征参数。Step 3: Further process the standardized images of the leaves obtained in the four key soybean growth stages of the soybean plants affected by diseases and insect pests at the flowering stage, pod setting stage, drumming stage and maturity stage, and for each soybean leaf with the identified disease and insect pest type Normalize the image to extract the corresponding feature parameters.
对于某一特定生长期内受到待鉴定病虫害影响大豆植株的叶片,为去除样本间特征参数提取结果的异质性,选择20个植株,每个植株随机选择一个叶片的标准化彩色图像。提取叶片标准化彩色图像的红、绿、蓝三个分量的图像矩阵R(x,y),G(x,y),B(x,y),然后对三个分量按公式F(x,y)=R(x,y)/3+G(x,y)/3+B(x,y)/3将叶片的标准化彩色图像转换为灰度图,灰度图提取结果如图1所示。对灰度图进行直方图统计,找到大豆叶片绿色部分和受病虫害影响的变色部分两个灰度部分峰值间的频数最低点对应的灰度值t1,将灰度大于阈值t1的所有像元亮度值设为0,灰度小于阈值t1的多有像元亮度值设为1,实现将灰度图转化为二值图,如图2所示。根据叶片图像宽度对应的像元个数值w1和叶片图像对应植株的测量结果w1’厘米,计算每个像元的实际边长z1=w1’/w1厘米,每个像元的实际面积则为z12平方厘米,根据像元的实际边长z1厘米、像元的实际面积z12平方厘米,以及每个斑点不同参数提取结果对应的像元数,计算该特定生长期内待鉴定的受病虫害植株叶片特征参数的均值,包括待鉴定病虫害植株叶片斑点面积比值a1、待鉴定病虫害植株叶片斑点个数值a2,待鉴定病虫害植株叶片斑点面积值a3,待鉴定病虫害植株叶片斑点周长值a4。For the leaves of soybean plants affected by the diseases and insect pests to be identified in a specific growth period, in order to remove the heterogeneity of the extraction results of characteristic parameters between samples, 20 plants were selected, and each plant randomly selected a standardized color image of a leaf. Extract the image matrices R(x,y), G(x,y), B(x,y) of the red, green, and blue components of the leaf standardized color image, and then press the formula F(x,y for the three components )=R(x,y)/3+G(x,y)/3+B(x,y)/3 Convert the normalized color image of the leaf into a grayscale image, and the grayscale image extraction result is shown in Figure 1 . Perform histogram statistics on the grayscale image to find the grayscale value t1 corresponding to the lowest frequency point between the two grayscale peaks of the green part of soybean leaves and the discolored part affected by diseases and insect pests, and set the brightness of all pixels whose grayscale is greater than the threshold t1. The value is set to 0, and the brightness value of many pixels whose grayscale is less than the threshold t1 is set to 1, so as to realize the conversion of the grayscale image into a binary image, as shown in Figure 2. According to the pixel value w1 corresponding to the width of the leaf image and the measurement result of the plant corresponding to the leaf image w1' cm, calculate the actual side length of each pixel z1=w1'/w1 cm, and the actual area of each pixel is z1 2 square centimeters, according to the actual side length of the pixel z1 cm, the actual area of the pixel z1 2 square centimeters, and the number of pixels corresponding to the extraction results of different parameters of each spot, calculate the plant affected by the disease and insect to be identified in this specific growth period The average value of the leaf characteristic parameters, including the ratio a1 of the leaf spot area of the plant to be identified, the number of leaf spots of the plant to be identified a2, the value of the leaf spot area of the plant to be identified a3, and the value of the leaf spot perimeter of the plant to be identified a4.
为进一步提取受病虫害影响叶片的整体形态特征,再次对受病虫害影响植株灰度图进行直方图统计,找到大豆叶片和背景色两个灰度部分峰值间的灰度级频数最低点对应的灰度值t1’,将灰度大于阈值t1’的所有像元亮度值设为0,将灰度小于t1’的所有像元亮度值设为1,实现对大豆整体叶片的二值化处理,如图3为大豆整体叶片的二值化处理结果,根据叶片二值图提取叶片长轴对应的像元数和短轴对应的像元数计算该特定生长期待鉴定病虫害植株叶片长宽比a5和该特定生长期待鉴定病虫害植株叶片椭圆率a6。In order to further extract the overall morphological characteristics of leaves affected by pests and diseases, the histogram statistics of the grayscale images of plants affected by pests and diseases were performed again, and the grayscale corresponding to the lowest point of grayscale frequency between the two grayscale peaks of soybean leaves and background color was found. value t1', set the brightness value of all pixels whose grayscale is greater than the threshold t1' to 0, and set the brightness value of all pixels whose grayscale is less than t1' to 1 to realize the binarization of the whole soybean leaf, as shown in the figure 3 is the result of binarization processing of the whole soybean leaf, extract the number of pixels corresponding to the long axis of the leaf and the number of pixels corresponding to the short axis of the leaf according to the binary map of the leaf to calculate the specific growth expectation and identify the plant leaf aspect ratio a5 of the disease and insect pest and the specific Growth expectation identification of plant leaf ellipticity a6 with pests and diseases.
对于某一特定生长期内受病虫害影响大豆植株的根系,为去除样本间指标参数提取结果的异质性,选择20个植株的根系标准化图像,提取叶片彩色图像的红、绿、蓝三个分量的图像矩阵R(x,y),G(x,y),B(x,y),对三个分量按公式F(x,y)=R(x,y)/3+G(x,y)/3+B(x,y)/3将彩色叶片图像转换为灰度图,如图4所示。For the roots of soybean plants affected by diseases and insect pests in a specific growth period, in order to remove the heterogeneity of the extraction results of index parameters between samples, 20 standardized root images of plants were selected, and the red, green and blue components of the leaf color image were extracted. The image matrix R(x,y), G(x,y), B(x,y), for the three components according to the formula F(x,y)=R(x,y)/3+G(x, y)/3+B(x,y)/3 converts the color leaf image to grayscale, as shown in Figure 4.
对根系灰度图像进行直方图统计,找到大豆根系和背景色两个灰度部分峰值间的灰度级频数最低点对应的灰度值t2,将灰度大于阈值t2的所有像元亮度值设为0,将灰度小于t2的所有像元亮度值设为1,实现将根系灰度图转换为二值图,如图5所示。Perform histogram statistics on the root gray image, find the gray value t2 corresponding to the lowest point of gray level frequency between the two gray peaks of soybean root and background color, and set the brightness value of all pixels whose gray scale is greater than the threshold t2 as If it is 0, set the brightness value of all pixels whose gray level is less than t2 to 1, so as to convert the root gray image into a binary image, as shown in Figure 5.
根据根系图像宽度对应的像元个数值w2和根系图像对应植株的测量结果w2’厘米,计算每个像元的实际边长z2=w2’/w2厘米,每个像元的实际面积则为z22平方厘米,根据像元的实际边长z2厘米、像元的实际面积z22平方厘米,以及受病虫害影响的植株根系不同参数提取结果对应的像元数,计算该特定生长期待鉴定病虫害植株根系面积值b1,待鉴定病虫害植株根系根瘤个数值b2,根据灰度图计算该特定生长期内待鉴定的受病虫害植株根系纹理特征参数的均值,根系纹理特征参数的具体计算过程如下,根据公式1计算待鉴定病虫害植株根系对比度纹理特征值b3,根据公式2待鉴定病虫害植株根系一致性纹理特征值b4,根据公式3待鉴定病虫害植株根系熵值纹理特征值b5,根据公式4待鉴定病虫害植株根系能量纹理特征值b6。According to the pixel value w2 corresponding to the width of the root image and the measurement result of the plant corresponding to the root image w2' cm, calculate the actual side length of each pixel z2=w2'/w2 cm, and the actual area of each pixel is z2 2 square centimeters, according to the actual side length of the pixel z2 cm, the actual area of the pixel z2 2 square centimeters, and the number of pixels corresponding to the extraction results of different parameters of the plant root system affected by the disease and insect pest, calculate the specific growth expectation to identify the root system of the plant with diseases and insect pests The area value b1, the number b2 of the root nodules of the plants to be identified with diseases and insect pests, the average value of the root texture feature parameters of the plants to be identified in the specific growth period is calculated according to the grayscale map, the specific calculation process of the root texture feature parameters is as follows, according to formula 1 Calculate the characteristic value b3 of the contrast texture of the root system of the plant to be identified, according to the formula 2, the texture characteristic value of the root system of the plant to be identified, b4, according to the formula 3, the entropy value of the root system of the plant to be identified, the texture characteristic value b5, and the root of the plant to be identified according to the formula 4 Energy texture feature value b6.
其中,p(i,j)表示由根系灰度图像计算的灰度共生矩阵在第i行,第j列位置上的值,表示根系灰度图像中灰度级为i的像元和灰度级为j的像元在固定的方向和像元间隔下同时出现的概率,n为根系灰度图像中i和j两个灰度级的差异,Ng表示根系灰度图像提取的灰度共生矩阵的级数。Among them, p(i,j) represents the value of the grayscale co-occurrence matrix calculated from the root grayscale image at the i-th row and j-th column position, and represents the pixel and gray level of the gray level i in the root grayscale image. The probability that the pixels with level j appear at the same time in a fixed direction and pixel interval, n is the difference between the two gray levels of i and j in the root gray image, and Ng represents the gray co-occurrence matrix extracted from the root gray image 's level.
对于某一特定生长期内受病虫害影响大豆植株的籽粒,为去除样本间指标参数提取结果的异质性,选择20个植株,提取每个植株籽粒彩色籽粒图像的红、绿、蓝三个分量的图像矩阵R(x,y),G(x,y),B(x,y),对三个分量按公式F(x,y)=R(x,y)/3+G(x,y)/3+B(x,y)/3将彩色籽粒图像转换为灰度图,如图6所示,For the grains of soybean plants affected by diseases and insect pests in a specific growth period, in order to remove the heterogeneity of the extraction results of index parameters among samples, 20 plants were selected, and the red, green and blue components of the color grain image of each plant's grain were extracted. The image matrix R(x,y), G(x,y), B(x,y), for the three components according to the formula F(x,y)=R(x,y)/3+G(x, y)/3+B(x,y)/3 converts the color grain image into a grayscale image, as shown in Figure 6,
对灰度图进行直方图统计找到大豆籽粒部分和试验台背景部分两个灰度部分峰值间的频数最低点对应的灰度值t3,将灰度大于阈值t3的所有像元亮度值设为1,灰度小于阈值t3的多有像元亮度值设为0,实现将籽粒灰度图转化为二值图,如图7所示。Perform histogram statistics on the grayscale image to find the grayscale value t3 corresponding to the lowest frequency point between the peaks of the two grayscale parts of the soybean grain part and the background part of the test bench, and set the brightness value of all pixels whose grayscale is greater than the threshold t3 to 1 , the brightness value of many pixels whose gray level is less than the threshold t3 is set to 0, and the grain gray level image is converted into a binary image, as shown in Figure 7.
根据受病虫害影响的籽粒图像宽度对应的像元个数值w3和根系图像对应植株的测量结果w3’厘米,计算每个像元的实际边长z3=w3’/w3厘米,每个像元的实际面积则为z32平方厘米,根据像元的实际边长z3厘米、像元的实际面积z32平方厘米,以及受病虫害影响的植株籽粒不同参数提取结果对应的像元数,计算该特定生长期内待鉴定的受病虫害植株籽粒特征参数的均值,包括待鉴定病虫害植株籽粒个数值c1,待鉴定病虫害植株籽粒面积值c2,待鉴定病虫害植株籽粒周长值c3,待鉴定病虫害植株籽粒长宽比值c4,待鉴定病虫害植株籽粒曲率值c5,待鉴定病虫害植株籽粒椭圆率值c6。Calculate the actual side length of each pixel z3=w3'/w3 cm according to the pixel number w3 corresponding to the width of the grain image affected by diseases and insect pests and the measurement result w3' cm of the corresponding plant in the root image, and the actual side length of each pixel is calculated. The area is z3 2 square centimeters. According to the actual side length of the pixel z3 cm, the actual area of the pixel z3 2 square centimeters, and the number of pixels corresponding to the extraction results of different parameters of plant grains affected by pests and diseases, the specific growth period is calculated. The average value of the grain characteristic parameters of the plants affected by diseases and insects to be identified, including the number of grains of the plants to be identified c1, the value of the grain area of the plants to be identified c2, the value of the grain circumference of the plants to be identified c3, and the ratio of the grains of the plants to be identified. c4, the grain curvature value c5 of the plants to be identified with pests and diseases, and the grain ellipticity value c6 of the plants to be identified with pests and diseases.
步骤4:利用开花期、结荚期、鼓粒期及成熟期的不同生长期内受到所有已知的病虫害类型影响的大豆植株图像建立特征参数数据库,对某一特定生长期内已知病虫害影响下大豆叶片特征参数的提取和建库结合步骤3中的方法和公式将该特定生长期下受到已知病虫害影响下的彩色叶片图像转换为灰度图,利用步骤3中的方法对灰度图进行直方图统计,找到叶片绿色部分和受病虫害影响的变色部分两个灰度部分峰值间的频数最低点对应的灰度值T1,将灰度大于阈值T1的所有像元亮度值设为0,将灰度小于阈值T1的多有像元亮度值设为1,实现将灰度图转化为二值图,根据叶片图像宽度对应的像元个数值W1和叶片图像对应植株的测量结果W1’厘米,计算每个像元的实际边长Z1=W1’/W1厘米,每个像元的实际面积则为Z12平方厘米,根据像元的实际边长Z1厘米、像元的实际面积Z12平方厘米,以及每个斑点不同参数提取结果对应的像元数,计算该特定生长期内已知病虫害类型的植株叶片特征参数的均值,包括已知病虫害植株叶片斑点面积比值A1、已知病虫害植株叶片斑点个数值A2,已知病虫害植株叶片斑点面积值A3,已知病虫害植株叶片斑点周长值A4,为进一步提取受到该已知类型病虫害影响叶片的整体形态特征,再次对受该已知病虫害影响植株灰度图进行直方图统计,找到大豆叶片和背景色两个灰度部分峰值间的灰度级频数最低点对应的灰度值T1’,将灰度大于阈值T1’的所有像元亮度值设为0,将灰度小于T1’的所有像元亮度值设为1,实现对大豆整体叶片的二值化处理,根据叶片二值图提取叶片长轴对应的像元数和短轴对应的像元数计算该特定生长期已知病虫害植株叶片斑点长宽比值A5和已知病虫害植株斑点椭圆率值A6。Step 4: Use the images of soybean plants affected by all known types of diseases and insect pests in different growth stages of flowering, pod setting, drumming and maturity to establish a feature parameter database, and the impact of known diseases and insect pests in a specific growth period. The extraction and database establishment of soybean leaf characteristic parameters are combined with the methods and formulas in step 3 to convert the color leaf image under the influence of known pests and diseases under the specific growth period into a grayscale image, and the method in step 3 is used to convert the grayscale image. Perform histogram statistics to find the gray value T1 corresponding to the lowest frequency point between the two gray value peaks of the green part of the leaf and the discolored part affected by diseases and insect pests, and set the brightness value of all pixels whose gray value is greater than the threshold T1 to 0, Set the brightness value of many pixels whose grayscale is less than the threshold T1 to 1, and convert the grayscale image into a binary image. , calculate the actual side length of each pixel Z1=W1'/W1 cm, the actual area of each pixel is Z1 2 square cm, according to the actual side length of the pixel Z1 cm, the actual area of the pixel Z1 2 square cm, and the number of pixels corresponding to the extraction results of different parameters of each spot, calculate the average value of the leaf characteristic parameters of plants with known pests and diseases in this specific growth period, including the ratio of leaf spot area of plants with known pests and diseases A1, and the leaves of plants with known pests and diseases. The number of spots is A2, the leaf spot area value of plants with known diseases and insect pests is A3, and the perimeter value of leaf spots of plants with known diseases and insect pests is A4. The histogram statistics of the plant grayscale map are used to find the grayscale value T1' corresponding to the lowest point of the grayscale frequency between the two grayscale peaks of soybean leaf and background color, and the brightness value of all pixels whose grayscale is greater than the threshold T1' is calculated Set to 0, set the brightness value of all pixels whose gray level is less than T1' to 1, realize the binarization processing of the whole soybean leaf, and extract the number of pixels corresponding to the long axis of the leaf and the corresponding short axis according to the leaf binary map. The number of pixels was used to calculate the leaf spot aspect ratio A5 of the plants with known plant diseases and insect pests and the ellipticity value A6 of the plant spots with known plant diseases and insect pests in the specific growth period.
对某一特定生长期内已知病虫害影响下大豆根系特征参数的提取和建库结合步骤3中的方法和公式将该特定生长期下受到已知病虫害影响下的彩色根系图像转换为灰度图,利用步骤3中的方法对灰度图进行直方图统计,找到大豆根系和背景色两个灰度部分峰值间的灰度级频数最低点对应的灰度值T2,将灰度大于阈值T2的所有像元亮度值设为0,将灰度小于T2的所有像元亮度值设为1,实现将根系灰度图转换为二值图,根据根系图像宽度对应的像元个数值W2和根系图像对应植株的测量结果W2’厘米,计算每个像元的实际边长Z2=W2’/W2厘米,每个像元的实际面积则为Z22平方厘米,根据像元的实际边长Z2厘米、像元的实际面积Z22平方厘米,以及受病虫害影响的植株根系不同参数提取结果对应的像元数,计算所有受到已知病虫害影响植株根系特征参数的均值,包括已知病虫害植株根系面积值B1,已知病虫害植株根系根瘤个数值B2,根据步骤三中的公式1计算已知病虫害植株根系对比度纹理特征值B3,根据步骤三中的公式2计算已知病虫害植株根系一致性纹理特征值B4,根据步骤三中的公式3计算已知病虫害植株根系熵值纹理特征值B5,根据步骤三中的公式4计算已知病虫害植株根系能量纹理特征值B6。Extraction and database establishment of soybean root characteristic parameters under the influence of known diseases and insect pests in a specific growth period Combined with the methods and formulas in step 3, the color root system images under the influence of known diseases and insect pests in the specific growth period are converted into grayscale images , use the method in step 3 to perform histogram statistics on the grayscale image, find the grayscale value T2 corresponding to the lowest point of the grayscale frequency between the two grayscale peaks of the soybean root system and the background color, and set the grayscale value greater than the threshold value T2. The brightness value of all pixels is set to 0, and the brightness value of all pixels whose grayscale is less than T2 is set to 1 to convert the root grayscale image into a binary image. According to the pixel value W2 and the root image corresponding to the width of the root image Corresponding to the measurement result of the plant W2' cm, calculate the actual side length of each pixel Z2 = W2'/W2 cm, the actual area of each pixel is Z2 2 square cm, according to the actual side length of the pixel Z2 cm, The actual area of the pixel is Z2 2 square centimeters, and the number of pixels corresponding to the extraction results of different parameters of the root system of the plant affected by diseases and insect pests. Calculate the average value of the characteristic parameters of the root system of all plants affected by known diseases and insect pests, including the root area value of plants affected by known diseases and insect pests B1 , the number B2 of root nodules of plants with known plant diseases and insect pests is known, the characteristic value B3 of root contrast texture of plants with known plant diseases and insect pests is calculated according to the formula 1 in step 3, and the characteristic value B4 of the consistency texture of roots of plants with known plant diseases and insect pests is calculated according to formula 2 in step 3, Calculate the root entropy texture feature value B5 of plants with known plant diseases and insect pests according to formula 3 in step 3, and calculate root energy texture feature value B6 of plants with known plant diseases and insect pests according to formula 4 in step 3.
对于某一特定生长期内受到已知病虫害影响大豆植株的籽粒,结合步骤3中的方法和公式将该特定生长期下受到已知病虫害影响下的彩色籽粒图像转换为灰度图,对灰度图进行直方图统计找到大豆籽粒部分和试验台背景部分两个灰度部分峰值间的频数最低点对应的灰度值T3,将灰度大于阈值T3的所有像元亮度值设为1,灰度小于阈值T3的多有像元亮度值设为0,实现将籽粒灰度图转化为二值图,根据受病虫害影响的籽粒图像宽度对应的像元个数值W3和根系图像对应植株的测量结果W3’厘米,计算每个像元的实际边长Z3=W3’/W3厘米,每个像元的实际面积则为Z32平方厘米,根据像元的实际边长Z3厘米、像元的实际面积Z32平方厘米,以及受病虫害影响的植株籽粒不同参数提取结果对应的像元数,计算该特定生长期内已知受病虫害影响的植株籽粒特征参数的均值,包括已知病虫害植株籽粒个数值C1,已知病虫害植株籽粒面积值C2,已知病虫害植株籽粒周长值C3,已知病虫害植株籽粒长宽比值C4,已知病虫害植株籽粒曲率值C5,已知病虫害植株籽粒椭圆率值C6。根据不同类型病虫害对大豆植株器官影响的差异,确定每种已知的大豆病虫害类型叶片特征参数的权重K1,确定每种已知的大豆病虫害类型根系特征参数的权重K2,确定每种已知的大豆病虫害类型籽粒特征参数的权重K3,同时满足权重参数的归一化特征,即K1+K2+K3=1。根据已经提取的待鉴定的大豆病虫害类型的不同器官图像特征参数提取结果和数据库中已知的病虫害种类的不同器官图像特征参数提取结果,建立病虫害判别模型,判别模型如公式5所示。For the grains of soybean plants affected by known pests and diseases in a specific growth period, the color image of grains affected by known pests and diseases in the specific growth period is converted into a grayscale image by combining the method and formula in step 3. Figure 1. Perform histogram statistics to find the gray value T3 corresponding to the lowest frequency point between the two grayscale peaks of the soybean kernel part and the background part of the test bench. The brightness value of many pixels less than the threshold T3 is set to 0, and the grayscale image of the grain is converted into a binary image. 'cm, calculate the actual side length of each pixel Z3=W3'/W3 cm, the actual area of each pixel is Z3 2 square cm, according to the actual side length of the pixel Z3 cm, the actual area of the pixel Z3 2 square centimeters, and the number of pixels corresponding to the extraction results of different parameters of plant grains affected by diseases and insect pests, calculate the average value of the characteristic parameters of plant grains known to be affected by plant diseases and insect pests in this specific growth period, including the number of known plant grains C1, The grain area value of known plants with diseases and insect pests is C2, the value of grain perimeter of plants with known plant diseases and insect pests is C3, the value of grain aspect ratio of plants with known plant diseases and insect pests is C4, the value of grain curvature of plants with known plant diseases and insect pests is C5, and the ellipticity value of grains of plants with known plant diseases and insect pests is C6. According to the difference in the effects of different types of diseases and insect pests on soybean plant organs, determine the weight K 1 of the leaf characteristic parameters of each known type of soybean diseases and insect pests, determine the weight K 2 of the root characteristic parameters of each known type of soybean diseases and insect pests, and determine the The weight K 3 of the grain characteristic parameters of the known soybean disease and insect pest types, and at the same time satisfies the normalized characteristics of the weight parameters, that is, K 1 +K 2 +K 3 =1. According to the extraction results of different organ image feature parameters of the soybean pests and diseases to be identified and the extraction results of different organ image feature parameters of known pests and diseases in the database, a pest discrimination model is established, and the discrimination model is shown in Equation 5.
步骤5:利用步骤4数据库中不同生长期下,受到不同种类已知大豆病虫害影响的植株特征参数以及数据库中已知的病虫害种类判别模型实现对受影响的大豆病虫害种类的识别和鉴定,对某一个待鉴定的受病虫害植株叶片特征参数提取结果a1、a2、a3、a4、a5、a6;待鉴定的受病虫害植株根系特征参数提取结果b1、b2、b3、b4、b5、b6;待鉴定的受病虫害植株籽粒特征参数提取结果c1、c2、c3、c4、c5、c6。将每一个待鉴定的受病虫害影响大豆的特征参数提取结果和该生长期下每一种已知受病虫害影响大豆的特征参数提取结果带入到步骤4中的判别模型中,根据判断模型计算待鉴定病虫害类型与所有已知病虫害类型的模型计算结果P值的大小,P值计算结果越小,证明待鉴定病虫害类型越接近该P值计算结果所对应的那一个已知的病虫害类型,从而实现对大豆病虫害类型的快速、准确、有效的识别和鉴定工作。Step 5: Use the characteristic parameters of plants affected by different types of known soybean diseases and insect pests in the database in step 4 under different growth periods and the discrimination model of known types of soybean diseases and insect pests in the database to realize the identification and identification of the affected soybean diseases and insect pests. A1, a2, a3, a4, a5, a6 of leaf characteristic parameter extraction results of a plant affected by diseases and insect pests to be identified; extraction results b1, b2, b3, b4, b5, b6 of root system characteristic parameters of a plant affected by diseases and insect pests to be identified; The extraction results of grain characteristic parameters c1, c2, c3, c4, c5, c6 of the plants affected by diseases and insect pests. The extraction results of the characteristic parameters of each soybean affected by diseases and insect pests to be identified and the extraction results of characteristic parameters of each known soybean affected by diseases and insect pests under this growth period are brought into the discriminant model in step 4, and the to-be-determined model is calculated according to the judgment model. The size of the P value of the model calculation result of identifying the type of pest and disease and all known types of pests and diseases. The smaller the calculation result of the P value, the closer the type of pest to be identified is to the known type of pest and disease corresponding to the calculation result of the P value, so as to achieve Fast, accurate and effective identification and identification of soybean pest types.
Claims (3)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011455125.5A CN112668404B (en) | 2020-12-10 | 2020-12-10 | Effective identification method for soybean diseases and insect pests |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011455125.5A CN112668404B (en) | 2020-12-10 | 2020-12-10 | Effective identification method for soybean diseases and insect pests |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112668404A CN112668404A (en) | 2021-04-16 |
CN112668404B true CN112668404B (en) | 2022-08-30 |
Family
ID=75402591
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011455125.5A Active CN112668404B (en) | 2020-12-10 | 2020-12-10 | Effective identification method for soybean diseases and insect pests |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112668404B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113310989B (en) * | 2021-04-26 | 2024-05-14 | 洛阳农林科学院 | Assessment method for severity of soybean disease green |
CN117787453B (en) * | 2023-05-17 | 2024-08-02 | 广东千沣农业科技有限公司 | Corn growth situation change trend prediction method, system and electronic equipment |
CN116310846B (en) * | 2023-05-19 | 2023-08-29 | 成都大学 | Disease and pest identification method and system based on image identification and data mining |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103559511A (en) * | 2013-11-20 | 2014-02-05 | 天津农学院 | Automatic identification method of foliar disease image of greenhouse vegetable |
CN207300299U (en) * | 2017-06-30 | 2018-05-01 | 中国农业科学院农业信息研究所 | A kind of soybean diseases and insect pests monitoring car and monitoring system |
CN108664979A (en) * | 2018-05-10 | 2018-10-16 | 河南农业大学 | The construction method of Maize Leaf pest and disease damage detection model based on image recognition and application |
CN111401245A (en) * | 2020-03-16 | 2020-07-10 | 吉林农业科技学院 | Crop disease and insect pest detection method, system and equipment based on image |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102334422B (en) * | 2010-07-27 | 2013-02-27 | 中国农业科学院蔬菜花卉研究所 | Method and system for real-time diagnosis of vegetable leaf diseases based on machine vision |
US8542939B2 (en) * | 2011-08-04 | 2013-09-24 | Arm Limited | Methods of and apparatus for using tree representations for representing arrays of data elements for encoding and decoding data in data processing systems |
CN103514459A (en) * | 2013-10-11 | 2014-01-15 | 中国科学院合肥物质科学研究院 | Method and system for identifying crop diseases and pests based on Android mobile phone platform |
CN105760871A (en) * | 2014-12-14 | 2016-07-13 | 仲恺农业工程学院 | Novel method for identifying plant leaf spot disease resistance |
TWI656334B (en) * | 2018-02-14 | 2019-04-11 | 國立中山大學 | A system for early detection of orchid pest by hyperspectral imaging techniques |
CN109063816B (en) * | 2018-08-28 | 2024-02-13 | 武汉古奥基因科技有限公司 | Bean pod number statistics device and method |
CN109470639A (en) * | 2018-11-07 | 2019-03-15 | 浙江大学 | Mulberry borer larvae and a rapid identification method for damage to mulberry leaves based on visible and near-infrared hyperspectral imaging |
CN110033015A (en) * | 2019-02-20 | 2019-07-19 | 天津科技大学 | A kind of plant disease detection method based on residual error network |
CN110188635B (en) * | 2019-05-16 | 2021-04-30 | 南开大学 | Plant disease and insect pest identification method based on attention mechanism and multi-level convolution characteristics |
CN112052904A (en) * | 2020-09-09 | 2020-12-08 | 陕西理工大学 | An identification method of pests and diseases based on transfer learning and convolutional neural network |
-
2020
- 2020-12-10 CN CN202011455125.5A patent/CN112668404B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103559511A (en) * | 2013-11-20 | 2014-02-05 | 天津农学院 | Automatic identification method of foliar disease image of greenhouse vegetable |
CN207300299U (en) * | 2017-06-30 | 2018-05-01 | 中国农业科学院农业信息研究所 | A kind of soybean diseases and insect pests monitoring car and monitoring system |
CN108664979A (en) * | 2018-05-10 | 2018-10-16 | 河南农业大学 | The construction method of Maize Leaf pest and disease damage detection model based on image recognition and application |
CN111401245A (en) * | 2020-03-16 | 2020-07-10 | 吉林农业科技学院 | Crop disease and insect pest detection method, system and equipment based on image |
Also Published As
Publication number | Publication date |
---|---|
CN112668404A (en) | 2021-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112668404B (en) | Effective identification method for soybean diseases and insect pests | |
CN109459392B (en) | A method for estimating aboveground biomass of rice based on UAV multispectral imagery | |
CN109916838B (en) | A kind of detection method of rice seed germination ability | |
CN102948282B (en) | Wheatear germination degree detection method | |
CN106525732B (en) | Rapid non-destructive testing method for internal and external quality of apples based on hyperspectral imaging technology | |
CN114049564B (en) | Pine wood nematode disease level prediction model construction method based on hyperspectral remote sensing image | |
US12106547B2 (en) | Method and system for screening spectral indexes of rice resistant to bacterial blight | |
CN102072885A (en) | Machine vision-based paddy neck blast infection degree grading method | |
CN105931223A (en) | Band ratio method based maize embryo segmentation method in high-spectral reflection image | |
CN114092839B (en) | Unmanned aerial vehicle remote sensing-based soybean harvest period maturity judging method | |
CN114494689A (en) | Identification method of tomato drought stress | |
CN107314990A (en) | A kind of spring maize remote sensing recognition method | |
CN115063690B (en) | A vegetation classification method based on NDVI time series characteristics | |
CN112668406B (en) | Effective screening method for soybean varieties | |
CN110221032B (en) | Rice quality detection method based on near infrared spectrum | |
CN109738438B (en) | A rapid method for measuring the senescence degree of wheat leaves | |
Gupta et al. | Applications of RGB color imaging in plants | |
Zhang et al. | Assessing spatio-temporal patterns of sugarcane aphid (Hemiptera: Aphididae) infestations on silage sorghum yield using unmanned aerial systems (UAS) | |
CN105426585B (en) | A kind of potato sprouting method for early warning based on SIN function fitting process | |
Taujuddin et al. | Detection of plant disease on leaves using blobs detection and statistical analysis | |
CN118155096A (en) | Rape maturity monitoring method based on unmanned aerial vehicle multispectral remote sensing | |
CN117253145B (en) | Method and system for determining growth situation of broom corn millet and mung bean in different intercropping modes | |
CN118298330A (en) | Unmanned aerial vehicle rapid diagnosis method and system for corn emasculation progress | |
WO2025138424A1 (en) | Spatio-temporal variation monitoring method for spad of winter wheat | |
CN116579446B (en) | A method for estimating wheat grain yield with high precision using phenotypic characteristics |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |