发明名称 基于条件互信息的遥感高光谱图像波段选择方法
摘要 本发明提出了一种基于条件互信息的遥感高光谱图像波段选择方法,包括下述步骤:A、打开高光谱遥感图像,人机交互标注待分类的样本;B、分组:根据步骤A得到的待分类样本,利用给定类别条件下各相邻波段之间的条件互信息进行波段分组;C、搜索:利用支持向量机和遗传算法相结合的搜索算法对步骤B得到的分组波段进行搜索计算,找到最优的波段组合;并在此基础上使用自适应分支定界法进行剪枝,通过条件互信息波段分组和自适应分支定界法剪枝的结合使用,避免了因为噪音扰动而出现冗余和噪音分组,减少了波段分组的数量,提高了波段组合的分类准确率。
申请公布号 CN101853392A 申请公布日期 2010.10.06
申请号 CN201010152915.6 申请日期 2010.04.21
申请人 河海大学 发明人 李士进;吴昊;余宇峰;朱跃龙;冯钧;万定生;郑伏广;王继民
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 许方
主权项 一种基于条件互信息的遥感高光谱图像波段选择方法,包括下述步骤:A、打开高光谱遥感图像,人机交互标注待分类的样本;B、分组:根据步骤A得到的待分类样本,将图像所有波段分为若干组;C、搜索:利用支持向量机和遗传算法相结合的搜索算法对步骤B得到的分组波段进行搜索计算,找到最优的波段组合;其特征在于:所述步骤B是使用条件互信息进行波段分组,具体包括下述步骤:B1、利用下述方法计算给定类别条件下各相邻波段之间的条件互信息:相邻两个波段Bi和Bi+1之间的条件互信息计算方法如下: <mrow> <mi>I</mi> <mrow> <mo>(</mo> <mi>C</mi> <mo>;</mo> <msub> <mi>B</mi> <mi>i</mi> </msub> <mo>|</mo> <msub> <mi>B</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>I</mi> <mrow> <mo>(</mo> <mi>C</mi> <mo>;</mo> <msub> <mi>B</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mrow> <mi>I</mi> <mrow> <mo>(</mo> <mi>C</mi> <mo>;</mo> <msub> <mi>B</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mi>H</mi> <mrow> <mo>(</mo> <msub> <mi>B</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mi>I</mi> <mrow> <mo>(</mo> <msub> <mi>B</mi> <mi>i</mi> </msub> <mo>;</mo> <msub> <mi>B</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow>其中,i是波段序号I(C;Bi)=H(C)+H(Bi)-H(C;Bi)是波段Bi和类别之间的互信息;I(C;Bi+1)=H(C)+H(Bi+1)-H(C;Bi+1)是波段Bi+1和类别之间的互信息;I(Bi;Bi+1)=H(Bi)+H(Bi+1)-H(Bi;Bi+1)是波段Bi和波段Bi+1之间的互信息; <mrow> <mi>H</mi> <mrow> <mo>(</mo> <mi>C</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>c</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> </munderover> <mi>P</mi> <mrow> <mo>(</mo> <mi>C</mi> <mo>=</mo> <mi>c</mi> <mo>)</mo> </mrow> <msub> <mi>log</mi> <mn>2</mn> </msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>C</mi> <mo>=</mo> <mi>c</mi> <mo>)</mo> </mrow> </mrow>是样本类别的熵,Nc是用户标注的待分类样本的类别总数; <mrow> <mi>H</mi> <mrow> <mo>(</mo> <msub> <mi>B</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mn>255</mn> </munderover> <mi>P</mi> <mrow> <mo>(</mo> <mi>X</mi> <mo>=</mo> <mi>j</mi> <mo>)</mo> </mrow> <msub> <mi>log</mi> <mn>2</mn> </msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>X</mi> <mo>=</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow>是划分样本后波段Bi的熵,其中X表示波段Bi中图象灰度的随机变量,j表示图像中灰度的值; <mrow> <mi>H</mi> <mrow> <mo>(</mo> <msub> <mi>B</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mn>255</mn> </munderover> <mi>P</mi> <mrow> <mo>(</mo> <mi>Y</mi> <mo>=</mo> <mi>j</mi> <mo>)</mo> </mrow> <msub> <mi>log</mi> <mn>2</mn> </msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>Y</mi> <mo>=</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow>是划分样本后波段Bi+1的熵,其中Y表示波段Bi+1中图象灰度的随机变量,j表示图像中灰度的值;B2、以波段序号为横坐标,条件互信息为纵坐标绘制条件互信息-波段曲线;B3、以上述条件互信息-波段曲线中的各极大值点为边界,将相邻两极大值点之间的连续波段作为一组。
地址 210098 江苏省南京市鼓楼区西康路1号
您可能感兴趣的专利