（1）
式中：
θ=arccos[（R-G）+R-B]/2（R-G）2+（R-B）（G-B）。
在本试验中，利用Matlab中的Imread函数分别提取彩色叶片图像中R、G、B3个通道的图像，再利用式
（1）将其转换为H、I、S模式。
然后，采用Otsu阈值法分别将H、I、S图像转换为二值化图像，以提取叶片病斑区域。
不妨将得到的病斑二值化图像仍记为H、I、S。
为了得到能区分作物病害的颜色特征值，分别计算病害叶片图像R、G、B分量和分割后的病斑图像的H、I、S分量的均值、方差、偏度、峰值、能量、熵，共6×6=36个统计特征参数[9]。
计算公式如下：
p（b）=h（b）/s
b=∑b=1bp（b）
bk=1/δ3∑（b-b）3p（b）
δ2=∑b=1（b-b）2p（b）
bF=1/δ4∑b=1（b-b）4p（b）-3
bN=∑b=1[p（b）]2
bE=∑b=1p（b）lg[p（b）]。
（2）
式中，p（b）为灰度级，h（b）为直方图，s为图像的面积，b为灰度均值，b为b的均值、bk为偏度，bF为峰值，bN为能量，bE为熵。
若直接把提取的36个病害图像特征参数用于构建病害识别模型，则识别效果的稳定性较差，因为各个特征对识别结果的贡献不同，有些特征甚至会影响识别结果。
因此，须要对得到的特征进行选择。
笔者使用基于区分矩阵的属性约简方法对得到的36个特征进行选择。
2基于区分矩阵的属性约简方法
基于区分矩阵的属性约简方法基本思想：
将属性在区分矩阵中出现的次数作为目标属性的重要性判断函数。
计算各属性的重要性判断函数f，由此可以确定第1个重要属性，此时重要属性的判断函数f值最大。
在第1个重要属性简单求出后，进一步快速简单地求出后续的重要属性。
设决策表S的核值集合为Core，其属性约简过程如下：
输入S的属性约简矩阵RM；输出S的属性约简Redu；
（1）初始化，令Redu=Core，H=；
（2）令Q={RM（aj）|aj∈Redu或AFI（aj）=0}，则RM=RM-Q，B=A-Redu-H；
（3）对aj∈B，计算f（aj）=Max{ak}，（k=1，2，…，m），其中f（x）为属性重要性判断函数；
（4）Redu←Redu∪{aj}；
（5）对每个aj∈B，令H2（aj）∩RM（ak），H3（aj）=|H2（aj）|；
（6）对所有aj∈B，令RM（aj）←RM（aj）-H2（aj），f（aj）←f（aj）-H3（aj）；若f（ak）=0，则H1=∪{ak}；
（7）重复步骤
（2）～（6），直到判断函数AFI为零向量；
（8）输出Redu。
Redu为信息表的一个属性约简。
3基于病害叶片的作物病害识别方法
由以上分析可得作物病害识别步骤如下：
（1）将采集到的每幅病害叶片图像转换为HIS模式；
（2）采用Otsu阈值法分别将叶片图像的H、I、S分量进行病斑分割，再转换为二值化图像，不妨还记为H、I、S；
（3）利用式
（2）分别提取H、I、S的36个分类特征，得到一个特征向量，然后生成包括条件属性集合和结论属性集合的满足粗糙集数据处理要求的二维关系规则表；
（4）将关系表中的每类特征属性进行归一化和离散化处理；
（5）利用基于区分矩阵的属性约简方法对关系表进行属性约简，依次消去可省略的属性（列）和合并重复的对象（行），再对每一个对象进行简化，消去冗余的属性值；
（6）根据一定的评选准则选取有效识别规则的属性简化表，最终获得优选的最简单判定规则；
（7）利用最近邻分类器对植物病害进行识别。
该分类器的基本思路是计算待测试样本与训练样本之间的欧式距离的最小值，由此确定待识别测试样本的类别。
4结果与分析
为了验证本研究提出的作物病害识别方法，采集温室黄瓜常见的细菌性角斑病、黄瓜炭疽病和黄瓜褐斑病的病害叶片各100幅，其中50幅图像作为训练集，其余50幅作为测试集。
本研究所使用的病害样本图像均是在陕西杨凌农业示范区科学园黄瓜温室采集得到的，黄瓜品种为中农26号。
以白色为背景色，在自然光照的非强光条件下对自然发病的黄瓜叶片进行图像采集得到黄瓜病害叶片图像。
以Matlab7.X软件中自带的图像处理工具箱为图像处理和分析平台，计算以上特征参数，数据统计分析采用常用统计分析SAS软件。
采用Otsu阈值法分别将病害叶片进行病斑分割，分割定位待识别的黄瓜叶片病害部分病斑，利用式
（2）提取病害叶片病斑的36个分类特征。
图1为3种黄瓜病害叶片图像及对应的分割病斑图像的H分量。
首先将得到的36个特征进行离散化，即把每个特征划分成有限个区域，本试验中取5个区域，使得每个区域中对象的决策值相同。
然后利用基于区分矩阵的属性约简方法对这36个特征进行属性约简，得到6个重要的属性，分别为病斑分量H的均值、方差、能量，病斑分量I的均值，病斑分量S的方差和能量。
再将这6个属性组成一个特征向量来表示这个样本。
最后利用最近邻分类器对病害种类进行分类，结果见表1。
为了说明本研究所提出方法的有效性，表1中给出了基于全部36个特征和基于其他2种方法[10-11]的识别结果。
由表1看出，本研究提出的方法的识别率最高，对黄瓜褐斑病的识别率高达94.26%，表明该方法用于黄瓜叶部病害图像识别是有效可行的。
表1不同黄瓜病害识别方法的识别结果
方法识别率（%）36个特征识别法78.59文献[10]82.88文献[11]90.35本研究所提出的方法94.26
叶片图像采集和识别试验结果显示：
（1）病害识别中并非特征越多越好，应对识别特征进行优化组合或特征选择；
（2）即使同一病害样本在不同成像环境下拍摄，图像特征向量也会存在一定差异；（3）不同黄瓜品种、不同发病时期的病症表现差异较大。
这些差异应与识别效果相结合，可以通过制定图像的采集规范予以解决。
5结论
本研究利用最大类间方差法分割黄瓜病斑图像，再提取病斑图像的36个特征，然后利用粗糙集对36个特征进行属性约简，得到6个重要特征，最后利用最近邻分类器对病害进行分类。
该方法提取病斑具有操作简单、分割效果好的特点。
本试验结果表明该方法是有效可行的。
由于病害叶片拍摄状态和光照等环境因素对颜色特征存在影响，发病程度及病斑的典型性也可能对识别准确率有影响，这些因素对拍摄照片中颜色特征的定量关系有待进一步研究。
对于某些颜色、形状特征非常相似的其他黄瓜病症，还应进一步提取病斑的颜色、形状和纹理等特征，综合考虑决定黄瓜病害症状的物理特征，以实现对病害的有效识别。
参考文献：
[1]李旺，唐少先.基于图像处理的农作物病害识别研究现状[J].湖南农机，2012，39
（1）：
176-178.
[2]QianYH，LiangJY，PedryczW，etal.Anefficientacceleratorforattributereductionfromincompletedatainroughsetframework[J].PatternRecognition，2011，44（8）：
1658-1670.
[3]胡可云，陆玉昌，石纯一.粗糙集理论及其应用进展[J].清华大学学报：
自然科学版，2001，41
（1）：
64-68.
[4]樊持杰，陈丽敏，夏春艳.基于RS的属性约简算法[J].微计算机信息，2010，26
（2）：
222-223.
[5]李智玲，胡?
?
.一种改进的区分矩阵属性约简算法[J].计算机系统应用，2008，17（10）：
51-55.
[6]王祥科，郑志强.Otsu多阈值快速分割算法及其在彩色图像中的应用[J].计算机应用，2006，26（增刊1）：
14-15.
[7]虎晓红，李炳军，席磊.基于多示例图的小麦叶部病害分割方法[J].农业工程学报，2012，28（13）：
154-159.
[8]TusharHJ，RavindraDB，PrashantGP.Cropdiseasedetectionusingimagesegmentation[J].WorldJournalofScienceandTechnology，2012，2（4）：
190-194.
[9]RevathiP，RevathiR.Knowledgediscoveryindiagnoseofcropdiseasesusingmachinelearningtechniques[J].InternationalJournalofEngineeringScienceandTechnology，2011，3（9）：
7187-7190.
[10]岑?
傣危?
李宝聚，石延霞，等.基于彩色图像颜色统计特征的黄瓜炭疽病和褐斑病的识别研究[J].园艺学报，2007，34（6）：
1425-1430.
[11]贾建楠，吉海彦.基于病斑形状和神经网络的黄瓜病害识别[J].农业工程学报，2013，29（增刊1）：
115-121.于翰林，卢泽民，朱咏莉.基于嵌入式Web服务器的pH值检测系统的设计[J].江苏农业科学，2015，43
（1）：
390-393.
升级会员 