Processing math: 100%

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于蚁群算法的光谱分解方法剔除植被干扰信息

夏浩东 薛云 邓会娟 刘凤姣

夏浩东, 薛云, 邓会娟, 等, 2012. 基于蚁群算法的光谱分解方法剔除植被干扰信息. 地质力学学报, 18 (1): 72-78.
引用本文: 夏浩东, 薛云, 邓会娟, 等, 2012. 基于蚁群算法的光谱分解方法剔除植被干扰信息. 地质力学学报, 18 (1): 72-78.
XIA Hao-dong, XUE Yun, DENG Hui-juan, et al., 2012. ELIMINATING THE DISTURBANCE OF VEGETATION INFORMATION BY SPECTRAL MIXTURE ANALYSIS BASED ON ANT COLONY ALGORITHM. Journal of Geomechanics, 18 (1): 72-78.
Citation: XIA Hao-dong, XUE Yun, DENG Hui-juan, et al., 2012. ELIMINATING THE DISTURBANCE OF VEGETATION INFORMATION BY SPECTRAL MIXTURE ANALYSIS BASED ON ANT COLONY ALGORITHM. Journal of Geomechanics, 18 (1): 72-78.

基于蚁群算法的光谱分解方法剔除植被干扰信息

详细信息
    作者简介:

    夏浩东(1968-), 男, 高级工程师, 主要从事地质矿产勘查与实物地质资料管理工作。E-mail:xiahaodong68@163.com

  • 中图分类号: P627

ELIMINATING THE DISTURBANCE OF VEGETATION INFORMATION BY SPECTRAL MIXTURE ANALYSIS BASED ON ANT COLONY ALGORITHM

  • 摘要: 针对遥感图像预处理工作中,光谱分解方法处理速度慢而蚁群算法识别目标速度快的特点,结合蚁群算法和线性光谱混合模型,建立基于蚁群搜索的光谱分解模型,以剔除植被干扰信息。选取青海黄南州吉地地区为研究区,首先确定蚂蚁移动规则,然后建立基于蚁群算法的光谱分解模型,最后根据模型重构不含有植被信息的新的多波段图像,通过残差图分析以及原图与剔除植被后影像对比分析,初步验证了基于蚁群算法的光谱分解方法剔除植被干扰信息的可行性。

     

  • 遥感技术在基础地质工作中应用越来越广泛,但是由于植被等因素的干扰,难以将植被覆盖下的信息从图像中分离出来,植被是遥感对地表探查的一道天然屏障。在基岩被植被完全覆盖的地区,遥感所获得的信息主要是植被信息,给遥感地质勘查工作带来极大不便。因此,找到一种可行的植被剔除方法,还原植被覆盖下的基岩信息,可更好地发挥遥感技术在基础地质工作中的应用能力。

    过去研究人员主要基于像元的植被指数对植被等干扰信息进行提取,但这种基于像元的模式不能有效地提高提取的精度。因为遥感图像中每个像元一般是多个地物的混合体,一个像元里记录了多类不同性质的地面目标,图像的光谱特征也是多个地物光谱特征的混合反映。混合像元的存在,产生了同物异谱、同谱异物现象,从而给地物的分类造成许多困难。因此,为了提高干扰信息提取的精度,就必须找到一种有效的进行混合像元分解的方法。线性光谱混合模型(Linear Spectrum Mixture Model,LSMM)是目前应用较多的混合像元分解方法[1~2],但LSMM因为要对整幅图像的每个像元都进行处理,所以处理速度较慢。

    蚁群算法是由意大利学者Dorigo等人,通过模拟自然界蚂蚁寻径的行为提出的一种全新仿生进化算法[3~4],是具有离散性、并行性、鲁棒性、正反馈性等特点的一种随机搜索方法。由于其概念简明、实现方便,迅速得到认可,并在优化问题求解、电力系统、计算机、冶金自动化等领域都有成功的应用[5]。但目前很少有人将蚁群算法引入遥感地质领域,由于遥感地质领域的信息提取也可以看作是一个组合优化问题,蚁群算法的全局性、离散性和基于概率选择路径等特点对于遥感图像非常适用。

    本文为了剔除植被干扰信息,综合考虑LSMM处理速度慢而蚁群算法识别目标速度快的特点,结合蚁群算法和线性光谱混合模型,建立基于蚁群搜索的光谱分解模型,剔除植被干扰信息,重构不含有植被信息的新的多波段图像,以期为后续基础地质工作提供基础影像。

    对于采用整数编码求解组合优化问题的蚁群系统,设求解问题因素有N个,蚁群中共有M只蚂蚁,τij(t)表示在t时刻ij组合之间信息素的数量。蚂蚁m在运动过程中根据各个路径上信息素的数量决定下一步的路径。用pijm(t)表示在t时刻蚂蚁m由城市i转移到城市j的概率,则:

    pmij(t)={ταij(t)ηβij(t)rTmταir(t)ηβir(t)jTm0otherwise
    (1)

    其中:Tm表示蚂蚁m下尝试组合情况的集合,该集合随蚂蚁m的行进过程而动态改变。

    信息量τij(t)随时间的推移会逐步衰减,用1-ρ表示它的衰减程度。经过n个时刻,要根据下式对各路径上的信息量作更新:

    τij(t+1)=ρτij(t)+Δτij
    (2)
    Δτij=Ms=1Δτkij
    (3)

    Δτijk表示蚂蚁m在本次循环中在组合ij之间留下的信息量,其计算方法根据蚁群系统的计算模型而定。

    正如Chialvo[6]和Millonas[7]所描述的,单只蚂蚁的状态可以用位置参数r和方向参数θ来表示。蚂蚁的状态转移规则可用一个加权函数表示:

    w(σ)=[1+σ1+δσ]β
    (4)

    这个方程描述了移动到信息素浓度为σ(r)的像素r处的相对概率。参数β表示一种随机度。β大,则w(σ)值较大,蚂蚁以较大的权系数跟随外激素浓度大的路径;反之,对蚂蚁路径选择影响不大。1/δ表示了蚂蚁感知外激素的能力。

    蚂蚁从像元k运动到像元i的归一化转移概率定义为:

    pik=w(σi)w(Δi)j/kw(σj)w(Δj)
    (5)

    w(Δ)是一加权因子,Δj表示蚂蚁在t-1时刻运动时的方向改变量,它的取值为8个离散的w值,一般蚂蚁的来向w定为1/20,原路径方向定为1,邻近的8个像元顺时针依次为1/20,1/12,1/4,1/2,1,1/2,1/4,1/12。

    蚂蚁爬行时,先计算8邻域每个像元的转移概率pik,然后根据轮赌算法选择移动到下一个像元。为了有利于全局搜索最优解,避免陷入局部最优解,设定了蚂蚁的兴奋阈值,当某支蚂蚁达到该值时,采用随机算法从8邻域中选择一个像元进行移动。

    线性光谱混合模型(Linear Spectrum Mixture Model,LSMM)主要的目的是绘制特定像元内地面物质的相对丰度图。因此,它必须要满足几个假设:

    ① 在一个像元里端元的混合必须是线性的。这种线性关系下光子在从太阳到传感器的路径上仅仅与一种端元成分交互,这要求端元组分分布在足够大的面积上。

    ② 在图像上所有的地物类型必须要有足够多的对照光谱,以便对它们进行分离。

    ③ 在图像里存在的端元组分必须在最少一个像元里保持纯状态,这个像元必须能正确识别且可空间定位。

    LSMM假定传感器测量的光谱是像元内所有成分光谱的线性组合[1]。LSMM模型描述如下:

    Ri=nk=1fkRik+εi
    (6)

    在这里i为光谱波段的数目;k是端元的数目;Ri是像元波段i的光谱反射率;fk是像元内端元k的面积比;Rik是波段i端元k的光谱反射率;εi是波段i的模拟误差。

    fk受下式约束:

    nk=1fk=10fk1
    (7)

    残差RMS用下式计算:

    RMS=(mi=1ε2i)/m
    (8)

    RMS残差要对所有图像像元都计算。RMS越大,模型的匹配就越差。所以残差图可以用来估计选择的端元是否适当,选择的端元数目是否足够。

    因为在LSMM模型中,端元反射率和像元反射率都是已知的,只有端元在像元里的面积比未知,只要波段的数目大于端元数目,就可以求出端元的面积比,从而得到各分量图像和残差图像,进而得到地物的分类图像。

    LSMM的优点在于它相对简单,并且是一种有物理意义的丰度测量方法;理论上有较好的科学性,对于解决像元内的混合现象有一定的效果。

    本文中基于蚁群算法的光谱分解算法采用的是混合像元分解蚁群算法,输入经过反射率转换预处理后的影像,输出不含有植被信息的新的多波段图像。

    【算法开始】

    ① 初始化

     1) 蚁群及相关参数;

     2) 信息素矩阵;

     3) 访问标志矩阵P

     4) 残差矩阵RMS

     5) 丰度矩阵F

    ② 迭代,直到达到停止条件

     1) 评价每只蚂蚁,对每只蚂蚁做

      如果没有访问过当前像元则完成下述工作

       (1) 建立端元系数矩阵x;(矩阵x主要存储像元内各子端元的面积)

       (2) 建立当前像元波段矩阵y;(矩阵y主要存储当前像元各波段的光谱反射率值)

       (3) 调用最小二乘法计算各端元的比例系数矩阵a;(重新计算各子端元的面积)

       (4) 对矩阵a进行负数漂移处理;

       (5) 对矩阵a进行归一化处理;

       (6) 计算该点残差存放到RMS矩阵;

       (7) 记录丰度信息到F

       (8) 根据矩阵a调整该像元的各波段值;

       (9) 设置当前像元访问标志;

       (10) 如果剔除端元所占比例超过某阈值,调整(增加)信息素;

     2) 每只蚂蚁爬行,做

      如果 蚁群迭代到某兴奋阈值 则

       蚂蚁采用随机算法从8邻域中选择一个像元进行移动

      否则

       (1) 根据公式(4)计算8邻域每个像元的w(σ)

       (2) 根据进入方向分配8邻域每个像元的权值

       (3) 根据公式(5)计算8邻域每个像元的转移概率pik

       (4) 按照轮赌算法选择移动到下一个像元

    ③ 按照BSQ格式保存图像数据;

    ④ 按照BSQ格式保存残差矩阵;

    ⑤ 按照BSQ格式保存剔除专题的丰度矩阵;

    【算法结束】

    选取青海黄南州吉地地区为试验区,该区为植被高度覆盖区。试验区遥感数据为Landsa7号卫星的ETM,景号为132-36,2000年7月18日获取的影像,太阳高度角为64.2°。试验区图像大小为408×402像素,经过反射率转换预处理后的影像作为输入图像。基于ENVI3.6及VC++6.0开发的蚁群光谱分解软件进行试验。

    通过求取植被指数,可以发现该区植被覆盖度为85.7157%(见图 1,绿色部分为植被),高植被覆盖率严重影响了矿化弱信息的提取。选用ETM1,2,3,4,5,7等6个波段进行主成分变换,其中PC1和PC2占95.3%。通过主成分分析,选用PC1和PC2制作二维散点图(见图 2)。通过散点图分析,确定图像上主要有植被及3种未知类别共4个端元。运行蚁群光谱分解软件,最后得到剔除植被后影像(见图 3)及残差图(见图 4)。

    图  1  预处理后影像(ETM7+ETM 4+ETM 1)
    Figure  1.  The image after pre-processing
    图  2  散点图
    Figure  2.  Scatter diagram
    图  3  剔除植被后影像(band 7+band 4+band 1)
    Figure  3.  The image after eliminating vegetation
    图  4  残差图(最大残差为0.103)
    蓝色:0~0.025;绿色:0.025~0.050;黄色:0.050~0.075;红色:0.075~0.103
    Figure  4.  The residual map

    从残差图上可以看出,最大残差为0.103,绝大部分残差在0.05以下。说明光谱分解的效果较好。对比图 2图 3可以发现,进行植被剔除后,原来绿色部分的植被覆盖区被还原为该区本来的面目。

    基于蚁群算法的光谱分解方法剔除植被信息,首次将蚁群这种全新的算法引入到遥感地质领域,综合考虑传统的光谱分解植被剔除方法处理速度慢和蚁群算法识别目标速度快的特点,通过残差图分析以及原图与剔除植被后影像对比分析,初步验证了基于蚁群算法的光谱分解方法来剔除植被信息的可行性。

    由于蚁群算法的理论体系还有待进一步完善,所以用这种全新的算法来处理遥感地质领域的一些复杂问题,还有许多工作要做。如蚁群的初始参数怎样优化,才能取得理想的效果,有待于进一步深入研究。

  • 图  1  预处理后影像(ETM7+ETM 4+ETM 1)

    Figure  1.  The image after pre-processing

    图  2  散点图

    Figure  2.  Scatter diagram

    图  3  剔除植被后影像(band 7+band 4+band 1)

    Figure  3.  The image after eliminating vegetation

    图  4  残差图(最大残差为0.103)

    蓝色:0~0.025;绿色:0.025~0.050;黄色:0.050~0.075;红色:0.075~0.103

    Figure  4.  The residual map

  • [1] Adams J B, Smith M O, ohnson P E.Pectral mixture modeling:A new analysis of rock and soil types at the Viking Lander 1 site[J].Journal of Geophysical Research, 1986, 91(B8):8098~8112. doi: 10.1029/JB091iB08p08098
    [2] Smith M O, Ustin S L, Adams J B, et al.Vegetation in deserts:A regional measure of abundance from multispectral images[J].Remote Sensing of Environment, 1990, 31(1):1~26. doi: 10.1016/0034-4257(90)90074-V
    [3] Dorigo M, Caro G D, Gambardella L M.Ant algorithms for discrete optimization[J].Artificial Life, 1999, 5(3):137~172.
    [4] Dorigo M, Bonabeau E, Theraulaz G.Ant algorithms and stigmergy[J].Future Generation Computer System, 2000, 16(6):851~871. http://www.sciencedirect.com/science/article/pii/S0167739X0000042X
    [5] 李士勇.蚁群优化算法及其应用研究进展[J].计算机测量与控制, 2003, 11(12):911~913. doi: 10.3321/j.issn:1671-4598.2003.12.001

    LI Shi-yong.Progresses in ant colony optimization algorithm with applications[J].Computer Measurement & Control, 2003, 11(12):911~913. doi: 10.3321/j.issn:1671-4598.2003.12.001
    [6] Chialvo D R, Millonas M M.How swarms build cognitive maps[C]//Steels L.The biology and technology of intelligent autonomous agents.NATO ASI Series, 1995:439~450.
    [7] 王树根, 杨耘, 林颖, 等.基于人工蚁群优化算法的遥感图像自动分类[J].计算机工程与应用, 2005, 29:77~80. doi: 10.3321/j.issn:1002-8331.2005.01.024

    WANG Shu-gen, YANG Yun, LIN Ying, et al.Automatic classification of remotely sensed images based on artificial ant colony algorithm[J].Computer Engineering and Application, 2005, 29:77~80. doi: 10.3321/j.issn:1002-8331.2005.01.024
  • 加载中
图(4)
计量
  • 文章访问数:  174
  • HTML全文浏览量:  87
  • PDF下载量:  4
  • 被引次数: 0
出版历程
  • 收稿日期:  2011-10-09
  • 刊出日期:  2012-03-01

目录

/

返回文章
返回