先选取的至少一种模型曲线（比如预先选取的至少一种模型 曲线包括线性曲线、复合曲线、生长曲线、指数曲线、戈珀兹曲线和逻辑斯蒂曲线）对玉米 大斑病不同发病时期的病情指数进行拟合，得到不同的模型曲线相应的以发病时期相应的 调查时间为自变量，玉米大斑病的病情指数为因变量的预测方程，并根据各个预测方程确 定出各个预测方程用于显著性检验的统计量（决定系数R2、残差平方和SSE、F值和检验概 率P);
[0058] 利用所述统计量对相应的预测方程进行显著性检验，从所述各个预测方程中筛选 出显著性最高的预测方程相应的预测模型作为最佳预测模型。
[0059] 本发明实施例中，利用SPSS筛选最佳预测模型一具体实施例的过程如下：首先选 择线性模型、复合模型、生长模型、指数模型、戈珀兹模型和逻辑斯蒂模型；然后，利用不同 时期的玉米叶部大斑病病情分级结果估计各个模型的参数，并输出各个模型回归方程显著 性检验的决定系数R2、残差平方和SSE、F值和检验概率P等统计量；最后，以各个统计量为 依据选择其中的最优模型（以R2越大、残差平方和SSE越小、F值越大、检验概率P越大为 原则筛选最优模型）。假设利用SPSS筛选出的最佳预测模型为逻辑斯蒂模型（Logistic)， 模型的回归方程为
其中：^是t时刻大斑病的病情指数；a，b是拟合系 数；c是病情指数增长率（拟定值）。将Logistic公式变形为t = ln((l/xt-a)/b)/c，公 式经过求导后得到 X't= -bcExp (ct) (a+bExp (ct)) 2，当 a = bExp (ct)，即
=-c/4a取最大值。因此，计算结果表明在调查时间t为
时病情发展比较快。
[0060] 如图3所示，本实施例公开一种计算玉米叶部大斑病病斑面积比例的装置，包括：
[0061] 图像采集单元1、图像处理单元2和面积比例计算单元3 ;其中，
[0062] 所述图像采集单元1，用于以白板为背景，采集玉米叶部的二维彩色图像；
[0063] 所述图像处理单元2,用于对所述二维彩色图像依次进行超绿特征提取、反色处 理、大津法二值化，得到玉米叶部图像，对所述二维彩色图像进行中值滤波、超绿特征提取、 大津法二值化、开运算，得到玉米叶部大斑病病斑图像，其中所述玉米叶部图像的目标为白 色，背景为黑色，所述玉米叶部图像的目标为白色，背景为黑色；
[0064] 所述面积比例计算单元3,用于对所述玉米叶部图像进行遍历，获取每一个像素点 的R通道值，计算R通道值为255的像素点数量，确定该数量为玉米叶片像素点数量，对所 述玉米叶部图像相对应的玉米叶部大斑病病斑图像进行遍历，获取每一个像素点的R通道 值，计算R通道为255的像素点数量，确定该数量为玉米叶部大斑病病斑处像素点数量，计 算玉米叶部大斑病病斑占玉米植株叶片面积的比例，计算公式为：面积比例=玉米叶部大 斑病病斑处像素点数量/玉米植株的叶片像素点数量。
[0065] 本发明实施例的计算玉米叶部大斑病病斑面积比例的装置，应用计算机图像处理 技术对玉米叶部大斑病病斑进行自动分割及面积比例计算，具有快速、准确、低成本的特 点；另外，可以长期存储和重复利用样本图像数据，有利于对宝贵样本的统一标准化管理。
[0066] SPSS软件集数据录入、资料编辑、数据管理、统计分析、报表制作、图形绘制为一 体，包括常规的集中量数和差异量数、相关分析、回归分析、方差分析、卡方检验、t检验和非 参数检验，也包括多元统计技术。
[0067] 虽然结合附图描述了本发明的实施方式，但是本领域技术人员可以在不脱离本发 明的精神和范围的情况下做出各种修改和变型，这样的修改和变型均落入由所附权利要求 所限定的范围之内。
【主权项】
1. 一种计算玉米叶部大斑病病斑面积比例的方法，其特征在于，包括： 51、 以白板为背景，采集玉米植株叶部的二维彩色图像； 52、 基于图像处理技术，根据所述二维彩色图像分别提取玉米叶部图像和玉米叶部大 斑病病斑图像，并根据所述玉米叶部图像计算玉米叶片像素点数量，根据所述玉米叶部大 斑病病斑图像计算玉米叶部大斑病病斑处像素点数量； 53、 根据所述玉米叶片像素点数量和玉米叶部大斑病病斑处像素点数量计算玉米叶部 大斑病病斑占玉米植株叶片面积的比例。2. 根据权利要求1所述的计算玉米叶部大斑病病斑面积比例的方法，其特征在于，所 述S2,包括： 521、 对所述二维彩色图像进行超绿特征提取； 522、 对所述步骤S21得到的图像进行反色处理，并进行大津法二值化，得到玉米叶部 图像，其中所述玉米叶部图像的目标为白色，背景为黑色。3. 根据权利要求1所述的计算玉米叶部大斑病病斑面积比例的方法，其特征在于，所 述S2,包括： 523、 对所述二维彩色图像进行中值滤波； 524、 对所述步骤S23得到的图像进行超绿特征提取； 525、 对所述步骤S24得到的图像进行大津法二值化； 526、 对所述步骤S25得到的图像进行开运算，得到玉米叶部大斑病病斑图像，其中所 述玉米叶部大斑病病斑图像的目标为白色，背景为黑色。4. 根据权利要求1所述的计算玉米叶部大斑病病斑面积比例的方法，其特征在于，所 述S2,包括： 527、 对所述玉米叶部图像进行遍历，获取每一个像素点的R通道值，计算R通道值为 255的像素点数量，确定该数量为玉米叶片像素点数量； 528、 对所述玉米叶部大斑病病斑图像进行遍历，获取每一个像素点的R通道值，计算R 通道为255的像素点数量，确定该数量为玉米叶部大斑病病斑处像素点数量； 529、 计算玉米叶部大斑病病斑占玉米植株叶片面积的比例，计算公式为：面积比例= 玉米叶部大斑病病斑处像素点数量/玉米植株的叶片像素点数量。5. 根据权利要求1所述的计算玉米叶部大斑病病斑面积比例的方法，其特征在于，还 包括： 根据所述玉米叶部大斑病病斑占叶片面积的比例对玉米叶部大斑病进行病情分级。6. 根据权利要求5所述的计算玉米叶部大斑病病斑面积比例的方法，其特征在于，还 包括： 获取玉米叶部大斑病病斑图像相应的玉米叶片的位置；其中， 所述根据所述玉米叶部大斑病病斑占叶片面积的比例对玉米叶部大斑病进行病情分 级，包括： 按照预设的病情分级标准，根据所述玉米叶部大斑病病斑占叶片面积的比例和玉米 叶部大斑病病斑图像相应的玉米叶片的位置对玉米叶部大斑病进行病情分级，所述病情分 级标准由大斑病病斑面积比例范围、病斑所在的玉米叶片的位置范围及对应的病情级数组 成。7. 根据权利要求5所述的计算玉米叶部大斑病病斑面积比例的方法，其特征在于，还 包括： 按照预设的抗性评价标准，根据所述玉米叶部大斑病病情分级结果对玉米叶部大斑病 进行抗性评价。8. 根据权利要求5所述的计算玉米叶部大斑病病斑面积比例的方法，其特征在于，还 包括： 根据所述玉米叶部大斑病病情分级结果建立玉米叶部大斑病预测模型，并利用所述玉 米叶部大斑病预测模型对玉米叶部大斑病疫情进行预测； 其中，所述S1，包括： 以白板为背景，采集玉米大斑病不同发病时期被调查玉米植株叶部的二维彩色图像。9. 根据权利要求8所述的计算玉米叶部大斑病病斑面积比例的方法，其特征在于，所 述根据所述玉米叶部大斑病病情分级结果建立玉米叶部大斑病预测模型，包括： 根据玉米大斑病同一发病时期玉米叶部大斑病病情分级结果计算该发病时期的病情 指数，计算公式为：其中，RV为玉米大斑病发病植株叶部大斑病病 情级数，N为病情级数是RV的发病植株数，NT代表调查玉米植株总数，RVmax为大斑病病情 级数的最大值； 以发病时期相应的调查时间为自变量，玉米大斑病的病情指数为因变量，利用SPSS中 曲线估计程序，采用预先选取的至少一种模型曲线对玉米大斑病不同发病时期的病情指数 进行拟合，得到不同的模型曲线相应的以发病时期相应的调查时间为自变量，玉米大斑病 的病情指数为因变量的预测方程，并根据各个预测方程确定出各个预测方程用于显著性检 验的统计量； 利用所述统计量对相应的预测方程进行显著性检验，从所述各个预测方程中筛选出显 著性最高的预测方程相应的预测模型作为最佳预测模型。10. -种计算玉米叶部大斑病病斑面积比例的装置，其特征在于，包括： 图像采集单元、图像处理单元和面积比例计算单元；其中， 所述图像采集单元，用于以白板为背景，采集玉米叶部的二维彩色图像； 所述图像处理单元，用于对所述二维彩色图像依次进行超绿特征提取、反色处理、大津 法二值化，得到玉米叶部图像，对所述二维彩色图像进行中值滤波、超绿特征提取、大津法 二值化、开运算，得到玉米叶部大斑病病斑图像，其中所述玉米叶部图像的目标为白色，背 景为黑色，所述玉米叶部图像的目标为白色，背景为黑色； 所述面积比例计算单元，用于对所述玉米叶部图像进行遍历，获取每一个像素点的R 通道值，计算R通道值为255的像素点数量，确定该数量为玉米叶片像素点数量，对所述玉 米叶部图像相对应的玉米叶部大斑病病斑图像进行遍历，获取每一个像素点的R通道值， 计算R通道为255的像素点数量，确定该数量为玉米叶部大斑病病斑处像素点数量，计算玉 米叶部大斑病病斑占玉米植株叶片面积的比例，计算公式为：面积比例=玉米叶部大斑病 病斑处像素点数量/玉米植株的叶片像素点数量。
【专利摘要】本发明公开一种计算玉米叶部大斑病病斑面积比例的方法及装置，能够快速、准确、低成本地对玉米叶部大斑病病斑面积比例进行计算。所述方法包括：S1、以白板为背景，采集玉米植株叶部的二维彩色图像；S2、基于图像处理技术，根据所述二维彩色图像分别提取玉米叶部图像和玉米叶部大斑病病斑图像，并根据所述玉米叶部图像计算玉米叶片像素点数量，根据所述玉米叶部大斑病病斑图像计算玉米叶部大斑病病斑处像素点数量；S3、根据所述玉米叶片像素点数量和玉米叶部大斑病病斑处像素点数量计算玉米叶部大斑病病斑占玉米植株叶片面积的比例。
【IPC分类】G06T1/00, G01B11/28
【公开号】CN105180850
【申请号】CN201510416395
【发明人】马钦, 苏琛, 刘峻明, 朱德海, 范梦扬, 王越
【申请人】中国农业大学
【公开日】2015年12月23日
【申请日】2015年7月15日