发明名称 基于分解的多目标遥感图像分割方法
摘要 本发明公开了一种基于分解的多目标遥感图像分割方法,主要解决现有图像分割技术中评价指标单一、计算复杂度高、分割效果不好的问题。实现步骤主要包括:输入待分割遥感图像;提取待分割图像特征;产生聚类数据;初始化初始种群;计算个体适应度值;初始化子问题;进化每一个子问题中的个体;判断是否满足终止条件分配类别标号;产生最优个体;输出分割图像。本发明提取图像每个像素的融合特征并结合分水岭粗分割产生超像素特征,通过分解的多目标方法,将多目标问题分解成一系列子问题实现遥感图像分割。本发明具有评价指标多样化、计算复杂度低、细节保持性能好等优点,对于图像的分割精度高、边缘定位准确,可用于对复杂图像的分割。
申请公布号 CN102800093B 申请公布日期 2015.12.09
申请号 CN201210241423.3 申请日期 2012.07.12
申请人 西安电子科技大学 发明人 李阳阳;焦李成;魏莹;刘若辰;缑水平;尚荣华;马文萍;于昕
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 程晓霞;王品华
主权项 一种基于分解的多目标遥感图像分割方法,包括以下步骤:步骤1:输入待分割遥感图像;步骤2:提取待分割图像特征;2.1利用灰度共生矩阵方法提取纹理特征向量;2.2利用非下采样小波分解方法获取小波特征向量;2.3用小波特征向量和纹理特征向量表示待分割图像的每一个像素点;步骤3:产生待聚类数据:用分水岭算法对待分割图像进行分水岭粗分割,得到不同的图像块;对每一个图像块的所有像素点特征取平均值,获得代表初始聚类数据的每一个超像素的特征向量,用所有超像素的特征向量作为待聚类数据;步骤4:利用待聚类数据初始化大小为N的初始种群X={x<sub>1</sub>,x<sub>2</sub>,...,x<sub>N</sub>},每个个体x<sub>n</sub>都代表一个分割结果,n=1,2,...,N,N为初始种群大小,取N=50;步骤5:根据类内方差F1和类间连接F2,计算每个个体的目标函数值F<sub>n</sub>F<sub>n</sub>=[F<sub>n1</sub>,F<sub>n2</sub>]其中,F<sub>n1</sub>=F1,F<sub>n2</sub>=F2;步骤6:初始化理想点Z<sup>*</sup>;其中<img file="FSB0000144242280000011.GIF" wi="305" he="83" /><img file="FSB0000144242280000012.GIF" wi="71" he="82" />是第1个目标函数F1到目前为止找到的最小值,<img file="FSB0000144242280000013.GIF" wi="64" he="72" />是第2个目标函数F2到目前为止找到的最小值;步骤7:将多目标问题F(x)=min(F1(x),F2(x))分解成N个子问题,每一个子问题的目标函数如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>min</mi><mi>imize</mi></mtd><mtd><msubsup><mi>g</mi><mi>j</mi><mi>te</mi></msubsup><mrow><mo>(</mo><mi>x</mi><mo>|</mo><msup><mi>&lambda;</mi><mi>j</mi></msup><mo>,</mo><msup><mi>z</mi><mo>*</mo></msup><mo>)</mo></mrow><mo>=</mo><mi>max</mi><mo>{</mo><msubsup><mi>&lambda;</mi><mi>i</mi><mi>j</mi></msubsup><mo>|</mo><mi>g</mi><mo>_</mo><msub><mi>F</mi><mi>ji</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>z</mi><mi>i</mi><mo>*</mo></msubsup><mo>|</mo><mo>}</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FSB0000144242280000014.GIF" wi="975" he="86" /></maths>其中,<img file="FSB0000144242280000015.GIF" wi="256" he="92" />表示第j个子问题的目标函数;<img file="FSB0000144242280000016.GIF" wi="267" he="95" />是第j个子问题的权值;<img file="FSB0000144242280000017.GIF" wi="187" he="121" /><img file="FSB0000144242280000018.GIF" wi="262" he="129" />j=1,2,...,N;x表示一个个体;g_F<sub>ji</sub>(x)表示第j个子问题的个体对应的第i个目标函数值;|是算术运算中的条件符号;||是算术运算中的绝对值符号;1≤i≤2;步骤8:根据每一个子问题<img file="FSB0000144242280000021.GIF" wi="278" he="104" />的权值λ<sup>j</sup>,计算每一个子问题的T个邻居子问题B(j)=(B<sub>j1</sub>,B<sub>j2</sub>,...,B<sub>jT</sub>),B<sub>ji</sub>表示第j个子问题的第i个邻居子问题的索引;取T=10;i=1,2,...,T;步骤9:根据初始种群X对每一个子问题<img file="FSB0000144242280000022.GIF" wi="260" he="99" />的父代个体P<sub>i</sub>(t)进行初始化,P<sub>i</sub>(t)=x<sub>i</sub>;其中t为迭代次数;t=0;步骤10:根据每个初始个体x<sub>j</sub>的目标值F<sub>j</sub>,对每个子问题<img file="FSB0000144242280000023.GIF" wi="268" he="100" />对应的两个目标函数值进行初始化,g_F<sub>j</sub>(t)=F<sub>j</sub>;其中t为迭代次数;t=0;步骤11:对每一个子问题<img file="FSB0000144242280000024.GIF" wi="257" he="89" />的个体进行进化操作11.1在第j个子问题的T个邻居子问题B(j)中随机选择3个邻居子问题s,k,l,对第s,k,l个邻居子问题的父代个体P<sub>s</sub>(t),P<sub>k</sub>(t),P<sub>l</sub>(t)进行交叉操作,得到一个新的临时子代个体Q<sub>j</sub>(t+1);11.2计算新的临时子代个体Q<sub>j</sub>(t+1)的两个目标函数值newF<sub>j1</sub>,newF<sub>j2</sub>,并根据newF<sub>j1</sub>,newF<sub>j2</sub>更新理想点Z<sup>*</sup>;通过新的临时子代个体Q<sub>j</sub>(t+1)和其目标值newF<sub>j1</sub>或newF<sub>j2</sub>来更新第j个子问题的所有T个邻居子问题的父代个体以及对应的目标函数值;步骤12:判断当前迭代次数t是否满足t<Tmax,如满足,则执行步骤13;否则,令迭代次数t加一次t=t+1,返回步骤11,其中Tmax为最大迭代次数,取Tmax=100;步骤13:分配类别标号:将每一个子问题<img file="FSB0000144242280000025.GIF" wi="258" he="92" />的父代个体P<sub>j</sub>(t)取出,作为最终的输出解集;对获得的解集中的每一个父代个体进行解码获得类别数及类别标号,作为可供选择的分割结果集合;步骤14:在可供选择的分割结果集合中产生最优分割结果;步骤15:输出分割图像。
地址 710071 陕西省西安市太白南路2号