发明名称 基于小波包变换和灰色预测模型的高光谱图像分类方法
摘要 一种基于小波包变换和灰色预测模型的新型高光谱图像分类方法,属于高光谱图像处理领域,其步骤如下:第一,获得待处理的高光谱数据;第二,应用小波包变换分解每个像元的高光谱响应曲线;第三,应用灰色预测模型处理分解结果;第四,使用特征构造结果对高光谱数据进行监督分类;第五,输出高光谱图像地物分类结果。该方法是一种自动的高光谱图像分类方法,能够有效的去除波段相关性,降低数据冗余度,避免维数灾难对分类精度的负面影响,适用范围广。
申请公布号 CN102663420A 申请公布日期 2012.09.12
申请号 CN201210078647.7 申请日期 2012.03.22
申请人 北京航空航天大学 发明人 尹继豪;徐胤;高超;顾则通;孙建颖;李辉
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 代理人
主权项 1.基于小波包变换和灰色预测模型的高光谱图像分类方法,其特征在于该方法包括如下步骤:(1)获得待处理的高光谱数据,对任一高光谱像元X,记为X=(x<sub>1</sub>,…,x<sub>i</sub>,…,x<sub>K</sub>),其中K是高光谱数据的波段总数,x<sub>i</sub>,i=1,…,K代表第i个波段的光谱响应数值;(2)应用小波包变换分解每个像元的高光谱响应曲线:给定小波母函数ψ及最大分解深度j,应用ψ对X=(x<sub>1</sub>,…,x<sub>i</sub>,…,x<sub>K</sub>)进行j层小波包变换,一共得到2<sup>j</sup>个分量,其中包括1个近似分量A和2<sup>j</sup>-1个细节分量<img file="FDA0000146048150000011.GIF" wi="255" he="56" />它们对应的能量系数分别是a和<img file="FDA0000146048150000012.GIF" wi="234" he="58" />并且满足关系式(1):<maths num="0001"><![CDATA[<math><mrow><msub><mi>d</mi><mn>1</mn></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>d</mi><mrow><msup><mn>2</mn><mi>j</mi></msup><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><mn>1</mn><mo>-</mo><mi>a</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,j是最大分解深度;(3)应用灰色预测模型处理分解结果:计算细节分量的能量系数序列<img file="FDA0000146048150000014.GIF" wi="210" he="58" />的1阶累积和序列<img file="FDA0000146048150000015.GIF" wi="233" he="66" /><maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>d</mi><mo>&OverBar;</mo></mover><mi>k</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>d</mi><mi>i</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msup><mn>2</mn><mi>j</mi></msup><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,j是最大分解深度;令<maths num="0003"><![CDATA[<math><mrow><mi>Y</mi><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>d</mi><mn>2</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>d</mi><mrow><msup><mn>2</mn><mi>j</mi></msup><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mi>B</mi><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mo>-</mo><mrow><mo>(</mo><msub><mover><mi>d</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>+</mo><msub><mover><mi>d</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>-</mo><mrow><mo>(</mo><msub><mover><mi>d</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mo>+</mo><msub><mover><mi>d</mi><mo>&OverBar;</mo></mover><mn>3</mn></msub><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>-</mo><mrow><mo>(</mo><msub><mover><mi>d</mi><mo>&OverBar;</mo></mover><mrow><msup><mn>2</mn><mi>j</mi></msup><mo>-</mo><mn>2</mn></mrow></msub><mo>+</mo><msub><mover><mi>d</mi><mo>&OverBar;</mo></mover><mrow><msup><mn>2</mn><mi>j</mi></msup><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>则:<maths num="0005"><![CDATA[<math><mrow><mfenced open='(' close=')'><mtable><mtr><mtd><mi>b</mi></mtd></mtr><mtr><mtd><mover><mi>b</mi><mo>^</mo></mover></mtd></mtr></mtable></mfenced><msup><mrow><mo>(</mo><msup><mi>B</mi><mi>T</mi></msup><mi>B</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>B</mi><mi>T</mi></msup><mi>Y</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,b称为发展系数,<img file="FDA00001460481500000110.GIF" wi="20" he="44" />是灰作用量;保留近似分量的能量系数a和发展系数b,则任一像元的特征构造结果由能量系数a和发展系数b组成;(4)使用特征构造结果对高光谱数据进行监督分类;(5)输出高光谱图像地物分类结果。
地址 100191 北京市海淀区学院路37号