植物叶片受害症状的自动识别方法技术

一种植物叶片受害症状的自动识别方法，它包括以下步骤：1)获取目标叶片图像；2)目标叶片图像处理；3)进行图像微分和区域标准色值分析；4)根据区域图像标准色值进行二次函数回归分析和回归参数计算；5)用方差分析法进行阈值参数的计算；6)用优化参数进行目标叶片受害症状判别分析。本发明专利技术通过应用图象测量和计算机识别技术实现了植物叶片“∧”和“∨”型受害症状的自动判别，通过受害症状来确定其受害机理，对准确定量地判别“∧”和“∨”型受害症状及其受害根源具有重要的指导作用；有利于解决目前植物叶片受害症状在科研和生产中的不确定性和争议；对于自动化判别植物叶片气象等灾害和病害具有重要理论和实际意义。

【技术实现步骤摘要】

本专利技术涉及植物气象灾害和病害等症状识别
，具体地说是一种植物叶片 Λ 和 V 型受害症状的自动识别方法。
技术介绍
纵观植物和树木等学科领域的研究历史，较为常见的植物/树木叶片受害症状有 叶斑、溃疡、花叶等等，而且相关的计算机识别技术已有报道。然而，至今尚未见到定量地 描述植物叶片 Λ 和 V 型受害症状的相关技术。 在自然生态系统中，植物/树木常遭受生物和非生物的危害而成灾。由可侵染性 的生物有机体造成的危害称为病害，而由非侵染性的环境因素，如干旱、冻害、营养匮乏等 等造成的危害为生理病害。无论侵染性病害还是生理病害皆表现出一定的受害症状。本 申请的专利技术人在研究植物/树木遭受气象灾害、机械伤害、病虫害等侵袭后表现出来的不 同症状中，分离出两种症状相近且容易混淆的 Λ 和 V 型植物叶片受害症状，并用数字 图像法进行了相应的描述和观测（WANG F. and Omasa Κ.，2012. Image measurements of leaf scorches on landscape trees subjected to extreme meteorological events. Ecological Informatics，12(2012) 16 - 22)。Λ和V型植物叶片受害症状在诸多的叶片受害症状 中较为常见且形态也很相似，其中， V 型叶片局部枯萎症状常表现为生物性病原的侵染 与叶片组织或细胞的抗性反应。由于侵染病原菌遭受叶片活组织或细胞的阻抗使得受害部 位从叶尖叶缘到叶内叶基渐窄。Λ叶片局部枯萎症状常表现为叶片内部系统性生理脱水 而始发于远离中央叶脉的叶尖叶缘部位，这为区分和定量观测此类症状奠定了基础。大量 观测发现，受害症状的内角（Internal angle of injured area简称IAIA)是区分二者的适 宜指标，IAIA是指枯萎部位的叶缘对应的圆心角，且大多生理性病害的IAIA大于等于180 度，而生物性侵染病害的IAIA往往小于180度，通过这些受害症状的类型的判别可自动识 别一些生理性干旱、冻害、营养匮乏等症状，以区别于生物性侵染病害，如炭疽病、腐烂病和 叶斑病等。如图1所示，IAIA是以枯萎部位在中央主脉上到达的最远端为顶点（图la和图 lb中的TP)，以该点到枯委部位在两侧叶缘上到达的最远端为边（图la和图lb中的BD)， 而测得的圆心角。如果枯萎部位远离中央主脉，则以枯萎部位距叶缘的最远端为IAIA顶点 (图lc中的TP)，以该点到枯萎部位在叶缘上到达的最远端为边（图lc中的BD)，由此测 得的圆内角为IAIA。目前IAIA的测定可以使用角度测量工具，如UTHSCSA Image T〇〇13. 00 进行测定，也可以使用目视法估测，以180度为判别标准，大于等于180度的为Λ型受害 症状（如图la所示），而小于180度的为V型受害症状（如图lb所示）但是，还没有出 现通过计算机自动识别IAIA的技术。
技术实现思路
针对上述现有技术的不足，本专利技术提供了一种植物叶片受害症状的自动识别方 法，其能够应用数字图像自动判别植物叶片的Λ和V型受害症状，以便及时准确地判 别受害的根源。 本专利技术解决其技术问题采取的技术方案是：， 其特征是，包括以下步骤： 1)获取目标叶片图像 以白色复印纸为背景，在光线良好的顺光条件下离体或活体拍摄而得到的高分辨 率、叶尖向上的公正RGB图像； 2)目标叶片图像处理 首先运行识别程序并赋值微分块数、微分方向、边缘处理方式等参数，并对参数有 效性判别后打开图像，然后分析计算目标叶片图像每一像素的R、G、B和L值最后检查目标 叶片图像的有效色值范围； 3)进行图像微分和区域图像标准色值分析 首先对目标叶片图像从左向右按等距离进行纵向微分得Λ X，然后根据公式（1) 计算每个Λ X段的R，G，Β以及L色彩平均值和G/L值， SDGi = 100 X (G/LrG/Lmin) / (G/Lma-G/Lmin) (1) 式中，SDGi为第i分段的标准色值，G/Li为第i单元的色值比，G/Lmin为单元色 值之极小值，G/Lmax为单元色值的极大值； 4)根据区域图像标准色值进行二次函数回归分析和回归参数计算 以标准色值为y轴，叶片微分单元为X轴构建二次曲线函数： y = aX2+bX+c (2) 式中，a为二次项系数、b为一次项系数、c为常数项，且a尹0， 求解式（2)的导函数式： 2aX+b = 0 (3) 并判别式（3)在等于零时的极值： X = -b/2a (4) 本文档来自技高网...

【技术保护点】
植物叶片受害症状的自动识别方法，其特征是，包括以下步骤：1)获取目标叶片图像以白色复印纸为背景，在光线良好的顺光条件下离体或活体拍摄而得到的高分辨率、叶尖向上的公正RGB图像；2)目标叶片图像处理首先运行识别程序并赋值微分块数、微分方向、边缘处理方式等参数，并对参数有效性判别后打开图像，然后分析计算目标叶片图像每一像素的R、G、B和L值最后检查目标叶片图像的有效色值范围；3)进行图像微分和区域图像标准色值分析首先对目标叶片图像从左向右按等距离进行纵向微分得△X,然后根据公式(1)计算每个△X段的R,G,B以及L色彩平均值和G/L值，SDGi＝100×(G/Li‑G/Lmin)/(G/Lmax‑G/Lmin)    (1)式中，SDGi为第i分段的标准色值，G/Li为第i单元的色值比，G/Lmin为单元色值之极小值，G/Lmax为单元色值的极大值；4)根据区域图像标准色值进行二次函数回归分析和回归参数计算以标准色值为y轴，叶片微分单元为X轴构建二次曲线函数：y＝aX2+bX+c    (2)式中，a为二次项系数、b为一次项系数、c为常数项，且a≠0，求解式(2)的导函数式：2aX+b＝0    (3)并判别式(3)在等于零时的极值：X＝‑b/2a    (4)式中，5)用方差分析法进行阈值参数的计算对目标叶片图像从上向下按等距离进行横向微分并按步骤3)进行图像微分和区域图像标准色值计算，对起始点进行位移若干像素后重复计算3到5次，以此为基础，进行方差阈值分析；目标叶片图像横向微分后的阈值Ts由式(5)进行判别：Ts(i)＝Minimum{(Gx(1),Gx(2),……Gx(n)}    (5)式中,n是微分的块数，通过式(6)可求得Gx(1),Gx(2),……Gx(n)：Gx(1)＝Fx(1)+Fy(n),Gx(2)＝Fx(2)+Fy(n‑1),……Gx(n)＝Fx(n)+Fy(1)    (6)而Fx(m)和Fy(m)可由式(7)和式(8)计算得到：Fx(m)=(mΣi=1k(x‾i-x‾)2/k-1)/(Σi=1kΣj=1n(xij-x‾i)2/k(m-1))---(7)]]>Fy(m)=(mΣi=1k(x‾i-x‾)2/k-1)/(Σi=1kΣj=n1(xij-x‾i)2/k(m-1))---(8)]]>式中，m＝1,2……n，Fx(m)是以阈值为分界点叶片上部和下部标准色值间方差分析结果的F值，Fx(m)是叶片上部和下部标准色值方差分析结果的F值，Xij是在i重复和j分块的标准色值，是i重复的标准色值平均值，是标准色值的总体平均值，k为重复次数；6)用优化参数进行目标叶片受害症状判别分析构建三组分段判别函数，如式(9)所示：fl(a,R,-b/2ax‾,Ts,Fv,Mv)=lnql+Cl0+Σi=13Clixi---(9)]]>式中，a为二次函数的二次项系数、R为二次函数的相关系数、为极值比、是标准色值的总体平均值、Ts为方差分析法获得的阈值、Fv为阈值两侧标准色值的方差分析F值、Mv为阈值两侧平均标准色值比；ql为样叶频率项、Cl0为判别函数的常数项、Cl1为判别函数的系数项、xi判别函数的自变量向量；根据式(9)进行判别：如果f1(a,R,-b/2ax‾,Ts,Fv,Mv)=max1<=l<=mfl(a,R,-b/2ax‾,Ts,Fv,Mv)}]]>则归为“∧”型受害症状组，如果f2(a,R,-b/2ax‾,Ts,Fv,Mv)=max1<=l<=mfl(a,R,-b/2ax‾,Ts,Fv,Mv)}]]>则归为“∨”型受害症状组，如果f3(a,R,-b/2ax‾,Ts,Fv,Mv)=max1<=l<=mfl(a,R,-b/2ax‾,Ts,Fv,Mv)}]]>则归为其他症状组；其中，m＝3。...

【技术特征摘要】
1.植物叶片受害症状的自动识别方法，其特征是，包括以下步骤： 1) 获取目标叶片图像 以白色复印纸为背景，在光线良好的顺光条件下离体或活体拍摄而得到的高分辨率、 叶尖向上的公正RGB图像； 2) 目标叶片图像处理 首先运行识别程序并赋值微分块数、微分方向、边缘处理方式等参数，并对参数有效性 判别后打开图像，然后分析计算目标叶片图像每一像素的R、G、B和L值最后检查目标叶片 图像的有效色值范围； 3) 进行图像微分和区域图像标准色值分析 首先对目标叶片图像从左向右按等距离进行纵向微分得Λ X，然后根据公式（1)计算 每个Λ X段的R，G，Β以及L色彩平均值和G/L值， SDGi = 100 X (G/U-G/L^) / (G/Lma-G/Lmin) (1) 式中，SDGi为第i分段的标准色值，G/Li为第i单元的色值比，G/Lmin为单元色值之 极小值，G/Lmax为单元色值的极大值； 4) 根据区域图像标准色值进行二次函数回归分析和回归参数计算 以标准色值为y轴，叶片微分单元为X轴构建二次曲线函数： y = aX2+bX+c (2) 式中，a为二次项系数、b为一次项系数、c为常数项，且a尹0, 求解式（2)的导函数式： 2aX+b = 0 (3) 并判别式（3)在等于零时的极值： X = -b/2a (4)5) 用方差分析法进行阈值参数的计算 对目标叶片图像从上向下按等距离进行横向微分并按步骤3)进行图像微分和区域图 像标准色值计算，对起始点进行位移若干像素后重复计算3到5次，以此为基础，进行方差 阈值分析； 目标叶片图像横向微分后的阈值Ts由式（5)进行判别： Ts(i) = Minimum {(Gx (1), Gx (2),......Gx (n)} (5) 式中，n是微分的块数， 通过式（6)可求得 Gx(l), Gx(2),......Gx(n): Gx(l) = Fx (1)+Fy (η), Gx (2) = Fx (2)+Fy ...

【专利技术属性】
技术研发人员：王斐，宋磊，大政鎌次，王静，
申请(专利权)人：山东省林业科学研究院，王斐，
类型：发明
国别省市：山东;37