发明名称 基于粒子群优化的高光谱遥感影像自适应波段选择方法
摘要 本发明公开了一种基于粒子群优化的高光谱遥感影像自适应波段选择方法,主要包括两个嵌套在一起的PSO算法,其中外部PSO负责自适应地优化波段数目,而内部PSO则负责优化指定数目的具体波段;本发明方法同时还改进了原始PSO算法中的目标函数。本发明能够解决高光谱遥感传统波段选择算法效果不好,以及不能事先确定波段数目的问题。本发明方法具有智能选择波段、适应性广等特点。
申请公布号 CN103150577B 申请公布日期 2017.02.08
申请号 CN201310101517.5 申请日期 2013.03.26
申请人 河海大学 发明人 苏红军
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 夏雪
主权项 一种基于粒子群优化的高光谱遥感影像自适应波段选择方法,其特征在于,包括以下步骤:步骤1,选择需进行降维的高光谱影像S;步骤2,随机初始化M个粒子;步骤3,对于外部PSO,其波段数目依据<img file="FDA0001126402750000011.GIF" wi="262" he="70" />进行迭代优化,x为待优化的波段数目,k为外部PSO的粒子序号,id为当前波段的数目;步骤4,令k=k+1,随机初始化波段数目<img file="FDA0001126402750000012.GIF" wi="158" he="63" />随机初始化内部嵌套PSO的波段<img file="FDA0001126402750000013.GIF" wi="282" he="71" />l为内部嵌套PSO的粒子序号,Y为内部嵌套PSO所对应的波段,<img file="FDA0001126402750000014.GIF" wi="155" he="79" />为当前波段数为id时内部嵌套PSO的粒子序号为l的粒子所处的位置;步骤5,对于内部嵌套PSO的每次迭代,执行以下步骤:(a)对于每一个内部嵌套PSO的粒子,评估其目标函数;(b)确定所有步骤(a)所述粒子的全局最优解<img file="FDA0001126402750000015.GIF" wi="123" he="70" />(c)对于每一个步骤(a)所述粒子,确定其历史局部最优解<img file="FDA0001126402750000016.GIF" wi="123" he="63" />(d)利用粒子速度更新公式和粒子位置更新公式更新所有内部嵌套PSO的粒子;其中,粒子速度更新公式和粒子位置更新公式分别由以下两个公式实现:<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>v</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>i</mi><mi>n</mi><mi>n</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>=</mo><mi>&omega;</mi><mo>&times;</mo><msubsup><mi>v</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>i</mi><mi>n</mi><mi>n</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>r</mi><mn>1</mn></msub><mo>&times;</mo><mrow><mo>(</mo><msubsup><mi>p</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>i</mi><mi>n</mi><mi>n</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>-</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>d</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>r</mi><mn>2</mn></msub><mo>&times;</mo><mrow><mo>(</mo><msubsup><mi>p</mi><mrow><mi>g</mi><mi>d</mi></mrow><mrow><mi>i</mi><mi>n</mi><mi>n</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>-</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>i</mi><mi>n</mi><mi>n</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001126402750000017.GIF" wi="1191" he="79" /></maths><maths num="0002"><math><![CDATA[<mrow><msubsup><mi>Y</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>i</mi><mi>n</mi><mi>n</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>=</mo><msubsup><mi>Y</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>i</mi><mi>n</mi><mi>n</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>+</mo><msubsup><mi>v</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>i</mi><mi>n</mi><mi>n</mi><mi>e</mi><mi>r</mi></mrow></msubsup></mrow>]]></math><img file="FDA0001126402750000018.GIF" wi="533" he="79" /></maths>式中,<img file="FDA0001126402750000019.GIF" wi="94" he="63" />为当前波段数为id时内部嵌套PSO的粒子速度;c<sub>1</sub>和c<sub>2</sub>分别为认知学习因子和社会学习因子;r<sub>1</sub>和r<sub>2</sub>均为独立随机变量,其取值范围为[0,1];ω为惯性权重,决定粒子速度的更新比例;<img file="FDA00011264027500000110.GIF" wi="106" he="63" />为内部嵌套PSO的所有粒子上一步迭代的局部最优解,<img file="FDA00011264027500000111.GIF" wi="108" he="71" />为内部嵌套PSO的所有粒子上一步迭代的全局最优解,<img file="FDA00011264027500000112.GIF" wi="126" he="70" />为当前波段数为id时内部嵌套PSO的粒子所处的位置,<img file="FDA00011264027500000113.GIF" wi="93" he="63" />为当前波段数为id时内部PSO的粒子所处的位置,x<sub>id</sub>为当前波段数为id时PSO的粒子所处的位置;(e)重复执行步骤a)–d),直到内部嵌套PSO收敛,将<img file="FDA0001126402750000021.GIF" wi="99" he="70" />的值作为波段数目<img file="FDA0001126402750000022.GIF" wi="130" he="63" />的最优值;步骤6,如果k&lt;M,转向步骤4;如果k=M,则检查外部PSO是否收敛,如果收敛则结束,此时外部PSO所有粒子的全局最优解<img file="FDA0001126402750000023.GIF" wi="107" he="71" />为最优化的波段数目,其对应的波段为最优化的波段;如果不收敛则执行以下步骤:a)对于每一个外部PSO的粒子,检索基于内部嵌套PSO得到的波段所对应的目标函数;b)确定步骤6的步骤a)中所述粒子的全局最优解<img file="FDA0001126402750000024.GIF" wi="131" he="77" />c)对于每一个外部PSO的粒子,确定其历史局部最优解<img file="FDA0001126402750000025.GIF" wi="131" he="62" />d)利用粒子速度更新公式和粒子位置更新公式更新所有外部PSO的粒子;其中,粒子速度更新公式和粒子位置更新公式分别由以下两个公式实现:<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>v</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>o</mi><mi>u</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>=</mo><mi>&omega;</mi><mo>&times;</mo><msubsup><mi>v</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>o</mi><mi>u</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>r</mi><mn>1</mn></msub><mo>&times;</mo><mrow><mo>(</mo><msubsup><mi>p</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>o</mi><mi>u</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>-</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>d</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>r</mi><mn>2</mn></msub><mo>&times;</mo><mrow><mo>(</mo><msubsup><mi>p</mi><mrow><mi>g</mi><mi>d</mi></mrow><mrow><mi>o</mi><mi>u</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>-</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>o</mi><mi>u</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001126402750000026.GIF" wi="1198" he="78" /></maths><maths num="0004"><math><![CDATA[<mrow><msubsup><mi>x</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>o</mi><mi>u</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>=</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>o</mi><mi>u</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msubsup><mo>+</mo><msubsup><mi>v</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>o</mi><mi>u</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msubsup></mrow>]]></math><img file="FDA0001126402750000027.GIF" wi="486" he="77" /></maths>式中,<img file="FDA0001126402750000028.GIF" wi="100" he="63" />为当前波段数为id时外部PSO的粒子速度;<img file="FDA0001126402750000029.GIF" wi="101" he="63" />为外部PSO的所有粒子上一步迭代的局部最优解,<img file="FDA00011264027500000210.GIF" wi="100" he="70" />为外部PSO的所有粒子上一步迭代的全局最优解,<img file="FDA00011264027500000211.GIF" wi="94" he="62" />为当前波段数为id时外部PSO的粒子所处的位置,x<sub>id</sub>为当前波段数为id时PSO的粒子所处的位置;e)令k=1,转向步骤4。
地址 210098 江苏省南京市鼓楼区西康路1号