专利名称：一种受病害胁迫农作物冠层高光谱图像的波段优选方法
技术领域：
本发明涉及高光谱图像数据分析领域，尤其涉及一种受病害胁迫农作物冠层高光谱图像的波段优选方法。
背景技术：
农作物受病虫害胁迫会直接影响其农产品的产量和质量。监测大田作物的冠层是获取作物病情指数信息的最佳方法之一。高光谱成像技术用于农作物冠层图像采集有两个优点(1)像素数据量大，能够分析小尺度空间的病害信息；( 单一像素的波段数据量大， 波段覆盖范围广，能够反映肉眼无法看到的光谱信息。这两个优点使得高光谱成像技术具有将受害像素与健康像素精确分类，从而估计作物病害程度的特性。但是高光谱图像的海量数据往往使得分析耗时较长，对处理器的性能要求也较高。因此寻找一种农作物冠层高光谱图像的波段优选方法来进行数据压缩非常必要。
现有的的高光谱波段优选方法主要有主成分分析法(PCA)、独立成分分析(ICA) 和遗传算法(Ga)波段优选方法。PCA和ICA需要很大的存储空间和计算复杂度；( 波段优选方法通过模拟自然进化过程搜索最优解，算法对新空间的探索能力有限，容易收敛到局部最优解，属于随机类算法，不能得到稳定的解，可靠性差，处理高光谱数据这样的高维数据需要较大的计算时间。发明内容
本发明提供了一种受病害胁迫农作物冠层高光谱图像的波段选择方法，解决了因农作物冠层高光谱图像像素点多、波段数大而需要耗费大量时间进行受病害胁迫和未受病害胁迫的像素点的分类分析，提高了农作物受病害感染的监测效率。
一种受病害胁迫农作物冠层高光谱图像的波段选择方法，包括
(1)采集受病害胁迫农作物冠层的高光谱图像，所述的高光谱图像包含 431-990nm范围内的953个光谱波段；
(2)计算得到所述高光谱图像每个波段的不稳定指数随机选取若干健康像素和病害像素，分为两组，健康像素和病害像素各为一组，不稳定指数ISIi的计算公式如下，Δs + sτ c* τwithin，il,i 2,i
ISI1 = --= ι-γA—. |mu-m2i|
式中ISIi表示第i个波段下的不稳定指数；Awithin, i是第i个波段下病害像素和健康像素的类内偏差，Abetweena是第i个波段下病害像素和健康像素的类间偏差是第 i个波段下病害象素的标准偏差，S2, i是第i个波段下健康像素的标准偏差；！^, i是病害象素在第i个波段下的平均光谱反射值，Hl2ii是健康像素在第i个波段下的平均光谱反射值；
(3)将所有不稳定指数ISIi升序排列，计算相邻不稳定指数的相对距离，RDi = ISIi-ISIw I，RDi表示第i个不稳定指数和第i+Ι个不稳定指数的相对距离；
(4)设定一系列平衡因子，平衡因子在所述相对距离RDi的最小值和最大值之间取
( 抽取其中一个平衡因子，计算每一所述相对距离RDi和抽取的平衡因子之间差值的绝对值，Di表示第i个相对距离与所抽取的平衡因子之间的差值的绝对值；
将Di按顺序累加，用ADi表示第1到第i-Ι个Di值的总和，ADi = D^D2+. . . +D^1， 求所有ADi局部极大值ADj, ADj对应的第j个波段为候选优选波段，其中j = 1，. . .，n，记录该平衡因子对应的候选优选波段组合及其波段数量η ；
(6)提取步骤(5)中的候选优选波段组合，建立最小二乘-支持向量机分类模型， 测算分类模型的分类精度；
(7)将其它所有平衡因子重复步骤(5)-(6)，根据每一个平衡因子对应的候选优选波段组合的波段数量η及分类模型的分类精度选取最优波段组合。
步骤中所述平衡因子在所述相对距离RDi的最小值和最大值之间按平均间隔取值，所述的平均间隔取值的间隔优选为0. 5。
步骤(6)中所述的建立最小二乘-支持向量机分类模型以步骤( 中所选像素的候选优选波段对应的光谱反射值作为输入变量，像素的类标识码作为输出变量。
每一个平衡因子的候选优选波段组合的波段数量η和最小二乘-支持向量机分类模型的分类精度是衡量优选波段质量的两个重要指标，波段数量越低，同时分类模型的分类精度越高，则选出来的优选波段组合质量越高。
本发明的有益效果
(1)实施方便快捷，基本采用简单的算术运算完成，运算速度较快；
(2)可以根据优选波段数量和分类精度的不同要求来选择不同的优选波段组合。
图1是AD局部极大值示意图2是水稻冠层高光谱图像的候选优选波段图。
具体实施方式
以受稻瘟病危害的水稻冠层为例，采集冠层高光谱图像，光谱范围覆盖 431-990nm，共953个光谱波段。
实施例1
(1)对冠层高光谱图像进行标准化矫正、背景分离等，剔除图像中非作物的对象像ο
标准化矫正公式
R= (I-Idark)/(Iwhite-Idark)
其中R是校准后的样本反射值，I是样本的原始反射值，Idark是黑标准象素反射值，Iwhite是白标准象素反射值。
背景分离方法公式[°031] R = R784. 55_R682. 27
其中R是差分运算后的象素值，R78455是784. 55nm波段的反射值，R682.27 ^ 682. 27nm波段的反射值。
按上述公式进行高光谱图像784. 55nm和682. 27nm的波段差分运算，经过差分运算后的单波段图像的植被象素和背景象素的象素值差别较大；取两类象素值的中间值作为阈值，大约为350，此时植被象素大于350，背景象素小于350 ；然后将大于350的象素设为 1，小于350的象素设为0，使图像二值化；然后将二值化后的图像做为掩模来将背景象素屏蔽。
(2)抽取图像中的健康像素和病害像素各50个，分成两组，健康像素和病害像素各为一组，对所选象素的每一个波段计算一个不稳定指数，公式如下
权利要求
1.一种受病害胁迫农作物冠层高光谱图像的波段优选方法，包括(1)采集受病害胁迫农作物冠层的高光谱图像；(2)计算得到所述高光谱图像每个波段的不稳定指数；(3)将所有不稳定指数升序排列，计算相邻不稳定指数的相对距离；(4)设定一系列平衡因子，平衡因子在所述相对距离的最小值和最大值之间取值；(5)抽取一个平衡因子，计算每一所述相对距离和抽取的平衡因子之间的差值的绝对值，Di表示第i个相对距离与所抽取的平衡因子之间的差值的绝对值；将Di按顺序累加，用ADi表示第1到第i-1个Di值的总和，ADi = DjD2+. · · +Di+求所有ADi的局部极大值ADj, ADj对应的第j个波段为候选优选波段，其中j = 1，. . .，n，记录该平衡因子对应的候选优选波段组合及其波段数量η ；(6)提取步骤(5)中的候选优选波段组合，建立最小二乘-支持向量机分类模型，测算分类模型的分类精度；(7)将其它所有平衡因子重复步骤(5)-(6)，根据每一个平衡因子对应的候选优选波段组合的波段数量η及分类模型的分类精度选取最优波段组合。
2.根据权利要求1所述的波段优选方法，其特征在于，步骤(1)中所述的高光谱图像包含431-990nm范围内的953个光谱波段。
3.根据权利要求1所述的波段优选方法，其特征在于，步骤(4)中所述平衡因子在所述相对距离的最小值和最大值之间按平均间隔取值。
4.根据权利要求3所述的波段优选方法，其特征在于，所述的平均间隔取值的间隔为0. 5。
全文摘要
本发明公开了一种受病害胁迫农作物冠层高光谱图像的波段优选方法，包括(1)采集农作物冠层高光谱图像数据；(2)执行不稳定指数计算、平衡因子设置等一系列的波段优选程序；(3)用候选的优选波段组合建立最小二乘-支持向量机分类模型，并估算分类精度；(4)根据优选波段组合内的波段数量和分类模型的分类精度进一步决定最佳的优选波段。本发明方法用于对受病害胁迫的农作物冠层高光谱图像像素进行分类分析，可以大大降低高光谱图像的数据量，操作方便，且可以根据优选波段数量和分类精度的不同要求来选择不同的优选波段组合，大大提高了农作物冠层病害胁迫程度分析的效率和灵活性。
文档编号G06T7/00GK102521830SQ201110391249
公开日2012年6月27日 申请日期2011年11月30日 优先权日2011年11月30日
发明者何勇, 徐兴, 赵芸 申请人:浙江大学, 浙江科技学院