基于测线图的水稻白叶枯病高光谱图像降维方法
【专利摘要】本发明公开一种基于测线图的水稻白叶枯病高光谱图像降维方法，属于高光谱信息提取领域。该方法包括如下步骤：首先应用连续投影算法在高光谱图像的光谱维选择特征波段图像，然后对特征波段图像建立相应的测线图，计算水稻叶片不同部位测线图的灰度值差异，进一步挑选出能够实现水稻叶片白叶枯病病斑检测所需的最少特征波段，从而降低了病斑检测所需的高光谱图像的维数。本发明方法将连续投影算法和测线图方法相结合，可以实现高维数据的有效降维，得到的特征图像可以准确地识别病斑；结合图像识别，准确计算出病斑面积及病害程度，以少量的特征图像实现水稻白叶枯病病斑的检测，降低了检测模型的复杂程度，有效地缩短了检测所花费的时间。
【专利说明】
基于测线图的水稻白叶枯病高光谱图像降维方法
技术领域
[0001] 本发明属于高光谱信息提取领域，涉及光谱数据分析，具体涉及一种基于测线图 的水稻白叶枯病高光谱图像降维方法。
【背景技术】
[0002] 高光谱图像是一个三维的数据立方体（Image Cube)，二维图像记录样本的空间形 态信息，每一个图像对应一个光谱波段。高光谱图像技术结合了样本光谱信息和图像信息， 实现作物的养分估算和病害检测。当水稻感染白叶枯病后，叶片出现的病斑会逐渐扩大，病 情加重;在农田管理中，需要根据病情的轻重进行农药的喷施。采用高光谱图像进行病情检 测的一个主要问题是获得的数据维数多达512维，由于检测中涉及的波段很多，处理速度 慢，并不是所有的波段的光谱信息及图像信息都可以揭示作物的生长状况，相反，许多数据 可能是一些无关紧要的噪声而将重要的信息掩盖。因此，为了实现病情检测的实时快速处 理，为了提高模型预测精度以及简化模型，需要对采集的数据进行优选，剔除不相关或非线 性变量，得到预测能力强、稳健性好的校正模型。
[0003] 基于光谱维的特征波段选择方法主要有无信息变量的消除法、间隔偏最小二乘 法、遗传算法、连续投影算法等。连续投影算法的工作原理是在光谱矩阵中寻找含有最低限 度的冗余信息的变量组，使变量间的共线性达到最小。
[0004] 基于叶片图像的测线图的方法，通过对叶片图像的不同部分灰度值进行分析，根 据不同波段下的灰度差异度来选取基于图像维的特征波段。
[0005] 由于高光谱图像检测涉及的波段很多，仅仅依靠光谱维的特征波段提取已不易得 到描绘病害特征的最佳组合。
【发明内容】
[0006] 为了克服现有技术的缺点与不足，本发明的目的在于提供一种基于测线图的水稻 白叶枯病高光谱图像降维方法。本发明结合连续投影算法和测线图法可以得到更好的特征 波段组合，降低病斑检测所需的维数，简化检测模型。
[0007] 本发明的目的通过下述技术方案实现：
[0008] -种基于测线图的水稻白叶枯病高光谱图像降维方法，包括如下步骤：
[0009] 首先应用连续投影算法在高光谱图像的光谱维选择特征波段图像，然后对特征波 段图像建立相应的测线图，计算水稻叶片不同部位测线图的灰度值差异，进一步挑选出能 够实现水稻叶片白叶枯病病斑检测所需的最少特征波段，从而降低了病斑检测所需的高光 谱图像的维数。
[0010] 所述的基于测线图的水稻白叶枯病高光谱图像降维方法，具体包括如下步骤：
[0011] 对η个样品的光谱数据的实现预处理，每一个样品的光谱数据为Xl=(Xll，Xl2,...， Xi5i2)T，其中i = l，2, · · ·，n;
[0012] 对样品Xl，对高光谱图像512维光谱信息应用连续投影算法选择m个特征波段，特 征波段集合为T= {ti，t2, · · ·，tm}，m<512;
[0013] 选取样品Xi的特征波段ti(l = l，2, · · ·，m)对应的特征图像料;(/ =私足…、叫；
[0014]在样品Xl的特征图像= 1,2f、、、. ； m)中绘制水稻叶片的测线图；
[0015]根据测线图计算样品Xl的特征图像恥(1 = 1,2，爾)中叶片健康部位(H)、病 斑部位(L)、阴影部位(S)的灰度平均值Gh，Gl，Gs;
[0016] 对样品Xi的特征图像與办=-计算Gh与Gl的差值Chl;当Chl2 100,表 明该特征图像可以区分健康部位和病斑部位，选取ti进入可区分健康部位和病斑部位的特 征波段集合Si;
[0017] 对样品Xi的特征图像抖;(卜? 2、rn)，计算Gh与Gs的差值Chs ;当Chs 2 100，表 明该特征图像可以区分健康部位和阴影部位，选取^进入可区分健康部位和阴影部位的特 征波段集合S2;
[0018] 对样品Xi的特征图像= 1; 2;、、…_m)，计算Gl与Gs的差值Cls;当Cls 2 100，表 明该特征图像可以区分病斑部位和阴影部位，选取^进入可区分病斑部位和阴影部位的特 征波段集合S3;
[0019]检测水稻叶片白叶枯病病斑所需特征波段为S = Si U S2 U S3,SET。病斑检测所需 的特征波段数量进一步降低，因此，可实现病斑检测的高光谱图像的降维。
[0020]对样品的光谱数据Xl进行的预处理由公式1确定，设第k维光谱反射率为Xlk(k=l， 2,...,512):
[0022] 其中，Wj表示采用最小二乘法得到的权值。
[0023] 所述的选择m个特征波段的算法，包括以下步骤：
[0024] 假设初始的迭代向量为Xk(0)，需要提取的变量个数为N，光谱矩阵的列变量数为J 个，则：
[0025] (1)在第1次迭代前(n=l)，任选光谱矩阵的任意1列j，把校正光谱阵第j列赋值给 xj，记为 xk(0);
[0026] (2)将还没有被选入的列向量未知的集合记为s，其中
[0027] ^ {i,!. < i < J\j ^ . , kiji
[0028] (3)分别计算幻对剩下列向量的投影：
[0030] (4)记k(n)=arg[max( | |Pxj| | )，jEs];
[0031] (5)令χ」= Ρχ」，j e s;
[0032] (6)n = n+l，如果n<N，回到第(2)步循环计算;最后提取出的变量：{xk(n) = 0, · · ·， N-l}〇
[0033]所述的叶片健康部位(Η)、病斑部位(L)、阴影部位(S)的灰度平均值Gh，Gl，Gs分别 由公式2、公式3、公式4确定：
[0034]设健康部位长度为1η、病斑部位长度为U、阴影部位长度为Is，则
[0038]所述的差值CHL由公式5确定：
[0039] Chl = Gh-Gl (公式 5)
[0040] 所述的差值Chs由公式6确定：
[0041 ] Chs = Gh-Gs (公式 6)
[0042]所述的差值CLS由公式7确定：
[0043] Cls = Gl-Gs (公式 7)。
[0044] 本发明相对于现有技术，具有如下的优点及效果：
[0045] (1)本发明对高光谱数据进行分析的方法有机结合了光谱维和图像维，对高光谱 数据进行2次降维，相比于仅从光谱维进行特征波段选取的方法，可以选取出较少的特征波 段。
[0046] (2)本发明的方法中采用的测线图方法是对图像的灰度进行运算，数据量小，处理 速度快。
[0047] (3)本发明方法将连续投影算法和测线图方法相结合，可以实现高维数据的有效 降维，得到的特征图像可以准确地识别病斑;结合图像识别，准确计算出病斑面积及病害程 度。
[0048] (4)本发明结合高光谱图像的光谱维及图像维，以少量的特征图像实现水稻白叶 枯病病斑的检测，降低了检测模型的复杂程度，有效地缩短了检测所花费的时间。
【附图说明】
[0049] 图1是本发明的基于连续投影算法和测线图方法对高光谱图像的降维流程图。
[0050] 图2是水稻叶片感染白叶枯病的病斑部位与健康部位的光谱曲线对比图。
[0051 ]图3是水稻叶片的连续投影算法得到的特征波段。
[0052]图4是叶片的特征光谱图像。
[0053]图5是水稻叶片的测线图。
[0054]图6是降维后特征图像的病斑分割结果图。
【具体实施方式】
[0055]下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限 于此。
[0056] 实施例1
[0057] 本发明的基于连续投影算法和测线图方法对高光谱图像的降维流程图，见图1所 示。本发明主要应用于感病水稻叶片的高光谱数据分析。当水稻感染白叶枯病后，叶片出现 的病斑和健康叶片在光谱反射率及图像上都表现出差异。从光谱反射率来看，染病叶片及 健康叶片的光谱曲线如图2所示。根据本发明，首先对叶片的光谱反射率进行处理和提取特 征波段。
[0058] 以多个样品的健康部位及病斑部位的光谱曲线为对象，对光谱数据进行平滑预处 理;每一个样品的光谱数据为Xi=(Xil，Xi2, . . .，Xi512)T，其中i = l，2, . . .，Π ;对高光谱图像 512维光谱信息，采用连续投影算法选取特征波段，实现高光谱图像在光谱维的降维，如图3 所示。选取了m个特征波段，T= {ti，t2,…，tm}，m<512。
[0059] 提取光谱维选取的特征波段对应的特征图像，如图4所示;在叶片上定义一条线贯 穿叶片，计算线对应点的灰度值，得到叶片图像对应的测线图，如图5所示;测线图的横轴为 长度，即该点与起点间的距离，纵轴为该点对应的灰度值。
[0060] 计算叶片^的特征图像中测线图叶片的健康部位(H)、病斑部位(L)、阴影部位(S) 的灰度平均值Gh，Gl，Gs。
[0061 ] 对样品Xi的特征波段ti(l = l，2, · · ·，m)对应的特征图像?二, ,、；计 算Gh与Gl的差值CHL。当100,表明该特征图像可以区分健康部位和病斑部位，选取ti进 入可区分健康部位和病斑部位的特征波段集合Si。
[0062] 对样品Xi的特征波段ti( 1 = 1，2，. . .，m)对应的特征图像2;.…计 算Gh与Gs的差值CHS。当CHS2 100,表明该特征图像可以区分健康部位和阴影部位，选取^进 入可区分健康部位和阴影部位的特征波段集合&。
[0063] 对样品Xi的特征波段ti( 1 = 1，2，. . .，m)对应的特征图像與又1 = ?，计 算&与gs的差值as。当CLS 2 100,表明该特征图像可以区分病斑部位和阴影部位，选取^进 入可区分病斑部位和阴影部位的特征波段集合&。
[0064] 实现水稻叶片白叶枯病病斑检测所需特征波段为S = Si U S2 U S3,SET，应用这些 特征波段实现叶片健康部分、病斑部分和阴影部分的区分，S对应的图像即为特征图像。 [0065]在特征图像中采用0TSU法进行图像分割，如图6所示;根据分割结果，可以准确地 计算病斑面积及病情等级。
[0066]上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的 限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化， 均应为等效的置换方式，都包含在本发明的保护范围之内。
【主权项】
1. 一种基于测线图的水稻白叶枯病高光谱图像降维方法，其特征在于包括如下步骤： 首先应用连续投影算法在高光谱图像的光谱维选择特征波段图像，然后对特征波段图 像建立相应的测线图，计算水稻叶片不同部位测线图的灰度值差异，进一步挑选出能够实 现水稻叶片白叶枯病病斑检测所需的最少特征波段，从而降低了病斑检测所需的高光谱图 像的维数。2. 根据权利要求1所述的基于测线图的水稻白叶枯病高光谱图像降维方法，其特征在 于具体包括如下步骤： 对η个样品的光谱数据的实现预处理，每一个样品的光谱数据为Xl=(Xll，Xl2,...， xi5i2)T，其中i = l，2, · · ·，n; 对样品Xl，对高光谱图像512维光谱信息应用连续投影算法选择m个特征波段，特征波段 集合为T= {ti，t2, · · ·，tm}，m<512; 选取样品xi的特征波段ti( 1 = 1，2，...，m)对应的特征图像'.'二1,2;…，?τ?); 在样品Xl的特征图像與.办二,,,m)中绘制水稻叶片的测线图； 根据测线图计算样品xi的特征图像= H >.、ρπ.)中叶片健康部位、病斑部位、阴 影部位的灰度平均值Gh，Gu Gs; 对样品xi的特征图像=1、2,…< ?川.)，计算GH与GL的差值CHL;当CHL 2 100,选取ti进 入区分健康部位和病斑部位的特征波段集合Si; 对样品Xi的特征图像. , 计算Gh与Gs的差值Chs;当ChsMOO,选取ti进 入区分健康部位和阴影部位的特征波段集合&; 对样品Xi的特征图像代;(〖二1,2,….计算GL与GS的差值C LS;当CLS2 100,选取ti进 入区分病斑部位和阴影部位的特征波段集合&; 检测水稻叶片白叶枯病病斑所需特征波段为S = Si U S2 U S3,SET。3. 根据权利要求2所述的基于测线图的水稻白叶枯病高光谱图像降维方法，其特征在 于： 对样品的光谱数据Xl进行的预处理由公式1确定，设第k维光谱反射率为Xlk(k=l， 2,...,512):其中，Wj表不米用最小二乘法得到的权值。4. 根据权利要求2所述的基于测线图的水稻白叶枯病高光谱图像降维方法，其特征在 于： 所述的选择m个特征波段的算法，包括以下步骤： 假设初始的迭代向量为xk(0)，需要提取的变量个数为N，光谱矩阵的列变量数为J个， 则： (1) 在第1次迭代前(n=l)，任选光谱矩阵的任意1列j，把校正光谱阵第j列赋值给^，记 为χ"0); (2) 将还没有被选入的列向量未知的集合记为s，其中(3) 分别计算&对剩下列向量的投影：(4) 记k(n)=arg[max( | |Pxj | |)，jes]; (5) 令 χ」= Ρχ」，jes; (6) n = n+l，如果n<N，回到第(2)步循环计算;最后提取出的变量：{xk(n) = 0,. . .，N-1}。5. 根据权利要求2所述的基于测线图的水稻白叶枯病高光谱图像降维方法，其特征在 于： 所述的叶片健康部位的灰度平均值Gh由公式2确定： 设健康部位长度为1h，则6. 根据权利要求2所述的基于测线图的水稻白叶枯病高光谱图像降维方法，其特征在 于： 所述的叶片病斑部位的灰度平均值&由公式3确定： 设病斑部位长度为U，贝1J7. 根据权利要求2所述的基于测线图的水稻白叶枯病高光谱图像降维方法，其特征在 于： 所述的叶片阴影部位的灰度平均值Gs由公式4确定： 设阴影部位长度为1 mi18. 根据权利要求2所述的基于测线图的水稻白叶枯病高光谱图像降维方法，其特征在 于： 所述的差值Chl由公式5确定： Chl = Gh-Gl (公式5)〇9. 根据权利要求2所述的基于测线图的水稻白叶枯病高光谱图像降维方法，其特征在 于： 所述的差值Chs由公式6确定： Chs = Gh~Gs (公式6) 〇10. 根据权利要求2所述的基于测线图的水稻白叶枯病高光谱图像降维方法，其特征在 于： 所述的差值as由公式7确定： Cls = Gl-Gs (公式7)。
【文档编号】G06T3/00GK105866040SQ201610177651
【公开日】2016年8月17日
【申请日】2016年3月25日
【发明人】张霞, 兰玉彬, 李继宇, 罗锡文, 周志艳, 魏玉
【申请人】华南农业大学