专利名称：基于免疫克隆和投影寻踪的sar图像分割系统及分割方法
技术领域：
本发明属于图像处理技术领域，涉及SAR图像分割，可用于军事侦察手段以及民 用和工业领域。
背景技术：
SAR作为一种不可或缺的军事侦察手段，在民用和工业上也有广泛的用途。随着 SAR成像技术的日益发展，所获得的SAR图像数据越来越多，机器解译已经逐步代 替了人工解译，从而SAR图像处理也成为研究热点。SAR图像分割作为SAR图像处 理的基本问题之一，它在对SAR图像的研究和应用中，通过图像分割可以发现感兴 趣的目标区域，从而为SAR图像后期的分类和识别奠定基础，而且目标识别的准确 性很大程度上依赖于图像分割的质量。为了辨识和分析SAR图像中的目标，需要将 它们从图像中分离提取出来，在此基础上才有可能进一步对目标进行测量和对SAR 图像进行利用。对于SAR图像来讲，由于斑点噪声的影响同类地物场景的非平稳程 度较强，很难有效地选取有代表性的样本用于参数训练，聚类算法作为无监督的学习 方法，其自适应和实用性的优点，使其成为图像分割的研究热点。现阶段基于聚类算 法的SAR图像分割方法主要有K均值聚类算法以及谱聚类算法，其中K均值聚类SAR 图像分割方法结果不稳定且分割精度不高，基于谱聚类的SAR图像分割方法不能有 效处理大规模图像分割问题。
发明内容
本发明目的在于克服上述已有技术的缺点，提出了一种基于免疫克隆和投影寻踪 的SAR图像分割系统及分割方法，以实现直接对图像特征进行操作，提高图像分割 精度。
为实现上述目的，本发明的SAR图像分割系统包括
图像特征提取模块，对输入的图像完成灰度共生矩阵、小波、brushlet和contourlet 特征提取，并将其提取结果传送至投影寻踪聚类模块；
投影寻踪聚类模块，对输入的图像特征采用免疫克隆算法进行自适应投影寻踪聚类，将聚类标签与图像像素对应，获得图像分割结果。 所述的投影寻踪聚类模块包括-
初始标签选择子模块，用于对所提图像特征，进行聚类获得初始的标签，并将标 签传输给投影方向选择子模块；
投影方向选择子模块，根据输入标签计算投影寻踪算法中的线性判别分析投影指 标，并通过免疫克隆算法优化该投影指标，获得最优的投影方向传至子空间聚类子模 块；
子空间聚类子模块，将图像特征映射到最优投影方向，获得最优子空间，对最优 子空间进行聚类，获得聚类标签，将聚类标签返回至初始标签选择子模块，迭代数次， 获得最终的聚类标签，将最终聚类标签与图像像素对应，获得最终图像分割结果。
为实现上述目的，本发明的SAR图像分割方法，包括如下步骤
(1) 对输入的图像提取灰度共生矩阵、小波、brushlet和contourlet特征，将提 取的图像特征定义为J^，其中，凍示样本的个数，在图像中表示为第/个像素点，_/ 表示样本的维数，在图像中表示为像素点的第/维特征；
(2) 对提取的图像特征J^进行归一化，设定标签迭代次数g;
(3) 将归一化后的图像特征Xg，随机投影到一个低维子空间t/,并对该子空间 C/进行K均值聚类，获得一个初始标签i/;
(4) 根据获得的初始标签F，计算投影寻踪的线性判别分析投影指标，并通过 免疫克隆算法优化该投影指标，获得最优的投影方向；
(5) 将归一化后的图像特征Z"，投影到最优投影方向，获得最优的子空间f/;
(6) 对子空间t/进行K均值聚类，获得新的标签if';
(7) 将新的标签i/'作为初始标签H，返回步骤(4)进行迭代，经过数次迭代 后，若获得新标签/f'在迭代过程中达到收敛，输出收敛的标签/T，若在迭代g次后 没有收敛，则输出最后获得的标签//';
(8) 将收敛的标签/T或最后获得的标签i/'与输入图像的像素对应，输出图像分 割结果。
本发明与现有的技术相比具有以下优点-
l.本发明由于采用免疫克隆算法优化投影寻踪的线性判别分析投影指标，因而可 获得最优的投影方向以及最优子空间；
2.本发明由于采用投影寻踪的线性判别分析投影指标与K均值聚类算法进行自适
5应的迭代，提高了图像分割的精度；
仿真结果表明，本发明较K均值算法精度有一定的提高。
图1本发明的SAR图像分割系统示意图； 图2是本发明的SAR图像分割方法流程图； 图3是本发明采用的SAR图像原图； 图4是采用K均值聚类算法对图3进行分割的结果图； 图5是采用本发明对图3进行分割的结果图。
具体实施例方式
参照图1,本发明基于免疫克隆和投影寻踪的SAR图像分割系统包括图像特征
提取模块和投影寻踪聚类模块。其中投影寻踪聚类模块包括初始标签选择子模块、投 影方向选择子模块和子空间聚类子模块。
图像特征提取模块，对输入的图像完成灰度共生矩阵、小波、brushlet和contourlet 特征提取，并将其提取结果传至投影寻踪聚类模块中的初始标签选择子模块；初始标 签选择子模块，对所提图像特征，进行简单聚类获得初始的标签，并将标签传输给投 影方向选择子模块；投影方向选择子模块，根据输入标签计算投影寻踪算法中的线性 判别分析投影指标，并通过免疫克隆算法优化投影指标，获得最优的投影方向，并将 所得投影方向传至子空间聚类子模块；子空间聚类子模块，将图像特征映射到最优投 影方向，获得最优子空间，对最优子空间进行聚类，获得聚类标签，将聚类标签返回 至初始标签选择子模块，迭代数次，获得最终的聚类标签，将最终聚类标签与图像像 素对应，获得最终图像分割结果。
参照图2，本发明的基于免疫克隆和投影寻踪的SAR图像分割系统及分割方法， 包括如下步骤
步骤l:提取输入图像的灰度共生矩阵、小波、brushlet和contourlet特征。
la)输入原始的SAR图像，其大小为MxM'，如图3所示。该输入图像为四幅
SAR图像，其大小分别为256X256;
lb)对输入的原始SAR图像提取灰度共生矩阵特征8维、两层小波特征7维、
两层brushlet特征8维和三层contourlet特征17维。将提取的图像特征定义为；，
其中，/表示样本的个数，在图像中表示为第/个像素点，/表示样本的维数，在图像
中表示为像素点的第/维特征。步骤2:对提取的图像特征^进行归一化，设定标签迭代次数g。
2a)对提取的图像特征J^采用下式进行归一化
其中A咖和x—分别为第y维的最大最小值，《为归一化后的图像特征；
2b)设定标签迭代次数g为10次。
步骤3:将归一化后的图像特征J^，随机投影到一个低维子空间f/，并对该子 空间t/进行K均值聚类，获得一个初始标签//。
3a)对归一化后的图像特征《采用下式进行随机投影
其中",是随机产生的投影方向；
3b)对随机投影后获得的低维子空间t/，采用下式进行K均值聚类-
min/,钩x广附』2
其中，^是K类随机投影后获得的低维子i间C/中像素特征的类内散度，夂是类别 的个数，、是第A类的像素，附i是第/t类像素的均值，",.每一类中像素的个数；
3c)分别计算每个像素到w的欧式距离，在保证《4最小的前提小，将具有最小欧 式距离的类别数赋予像素的标签，获得一个初始标签^。
步骤4:设定参数并按照参数随机生成种群，将投影方向作为免疫克隆算法优化 的目标。 '
4a)设定进化代数gw为500，种群大小W为10，变异概率;^为0.9，克隆规模 为2;
4b)根据种群大小iV，随机产生初始种群^大小为wxiV，其中m为图像特征维数。
步骤5:根据所获得初始标签//，计算初始种群^中的个体适应度。
5a)根据所获得初始标签/f，采用下式计算类间散度&和类内散度^:
<formula>formula see original document page 7</formula>
其中《是第z'类图像特征x&的均值，即《-iJ^/巧，^是归一化后的图像特征;^的总均值，即1=|:￡;^/"，其中"=|>;， w,.为第/类图像特征《中像素点的个数， "为归一化后的图像特征^;.中总的像素点的个数，c为类别数。
5b)根据获得的类间散度&和类内散度&，采用下式计算初始种群^中的每
体的适应度
/(")=
1-
0 for
其中&=&+^为总散度，"为初始种群」中任意个体。
步骤6:根据克隆规模 ，对初始种群^进行克隆操作，获得克隆后的种群4。 6a)根据设定的克隆规模"e，采用下式对初始种群进行克隆操作
"(a) ="2 ，…,a"c ]
其中，7f表示克隆算子，fl为初始种群J中任意一个个体，a,，《，…，气分别为
对fl进行规模为"e的克隆操作后的克隆个体；
6b)对初始种群J中的每个个体进行克隆操作，获得克隆后的种群4。
步骤7:根据变异概率，对克隆后种群^进行变异操作，获得变异后的种群《。
7a)根据设定的变异概率&，采用下式对克隆后种群4中变异的个体进行变异
操作
]S' = P.邵(t'iV (0,1) + tW (O，l))
"+ y3'iV(0,1)
其中，G表示变异算子，s为概率/^选择出的变异个体，w为选择出的变异个体对
应的变异后个体，t和t'分别为(^广和(V^)—'， )3和)8'为设定参数，i^0,^是满 足均值为0,方差为1正态随机变量；
7b)对克隆后种群4中所有变异个体完成变异操作，没有变异的个体保持不变，
获得变异后的种群4 。
步骤8:计算变异后的种群4的适应度函数值，将种群4 中个体适应度与初始种
群J中个体适应度进行对比，依次选出适应度较大的个体，获得新的种群乂。
8a)比较变异后的种群《中个体适应度与初始种群^中对应原始的个体适应度，
8若/(C("))>/(")，则选择d")为新的种群乂，否则选择"为新的种群乂，其中， /(口)为适应度函数，7f表示选择算子，r,(a)为种群《中个体，a为对应的初始种 群j中个体。
8b)比较种群A与初始种群」中对应所有个体适应度，选出适应度较大的个体为. 新的种群乂。
步骤9:将新的种群乂作为初始种群J，返回步骤5，按照设定的进化代数gw进 行迭代。
步骤10:在进化代数^中，若每一代的最大适应度函数值收敛，则终止迭代，
输出最大适应度函数值。
步骤ll:将归一化后的图像特征Zy，釆用下式投影到最优的投影方向，获得最 优的子空间C/':
其中fl'是通过免疫克隆算法选择的最优的投影方向。
步骤I2:对投影后获得的低维子空间C/进行K均值聚类，获得新的标签H'。
12a)对投影后获得的低维子空间f/'，采用下式进行K均值聚类
<formula>formula see original document page 9</formula>
其中，^是尺类低维子空间C/中像素特征的类内散度，《是类别的个数，x,是第A类 的像素，m,.是第ifc类像素的均值，Wi每一类中像素的个数；
12b)分别计算每个像素到mi的欧式距离，在保证^最小的前提小，将具有最小 欧式距离的类别数赋予像素的标签，获得新的标签/Z'。
步骤13:将新的标签i/'作为初始标签//，返回步骤5进行迭代，经过多次迭代
后，若获得新标签i/'在迭代过程中达到收敛，输出收敛的标签/T，若在设定的迭代 次数g后没有收敛，则输出最后获得的标签W;
步骤14:将收敛的标签^或最后获得的标签//'与输入图像的像素对应，输出图
像分割结果。
本发明的效果可以通过以下对SAR图像分割进一步说明1、 仿真条件
本发明的仿真在windows XP， SPI， CPU Pentium (R) 4，基本频率2.4GHZ， 软件平台为Matlab7.0.1运行。仿真选用四幅SAR图像原图，如图3，其中图3(a)是 树木、桥梁、田野和河流，图3(b)是陆地和湖泊，图3(c)是陆地和河流，图3(d)是树 木和田野。
2、 仿真内容
(1) 用现有的K均值方法分别对图3(a)、图3(b)、图3(c)和图3(d)进行图像分害iJ。 仿真结果如图4所示。其中，图4(a)是对图3(a)的分割结果图，图4(b)是对图3(b)的 分割结果图，图4(c)是对图3(c)的分割结果图，图4(d)是对图3(d)的分割结果(2) 用本发明的方法分别对图3(a)、图3(b)、图3(c)和图3(d)进行图像分害U。仿真 结果如图5所示其中，图5(a)是对图3(a)的分割结果图，图5(b)是对图3(b)的分割结 果图，图5(c)是对图3(c)的分割结果图，图5(d)是对图3(d)的分割结果图。
.3、仿真结果分析
从图4(a)与图5(a)的分割结果图比较可以看出，图5(a)能够分割出桥梁下方树 木的阴影处，而图4(a)没有分割出。
从图4(b)与图5(b)的分割结果图比较可以看出，图5(b)图像陆地左下角的错分 点要少于图4(b)的分割结果。
从图4(c)与图5(c)的分割结果图比较可以看出，图5(c)在陆地部分没有杂点且 河流宽度与原图相近，而图4(c)的图像顶部有少量的杂点且分割结果的河流宽度大于 原始图3(c)。
从图4(d)与图5(d)的分割结果图比较可以看出，图5(d)左上角以及右下角田野 处错分点少于图4 (d)的分割结果，且图5 (d)的树木的下边缘较图4 (d)光滑。 由图4与图5比较可以看出本发明具有分割精度高的优点。
该实例是在以本发明技术方案为前提下进行实施的，给出了详细的实施方式和具 体的操作过程，但本发明的保护范围不限于上述实施例。
10
权利要求
1、一种基于免疫克隆和投影寻踪的SAR图像分割系统，包括图像特征提取模块，对输入的图像完成灰度共生矩阵、小波、brushlet和contourlet特征提取，并将其提取结果传至投影寻踪聚类模块；投影寻踪聚类模块，对输入的图像特征采用免疫克隆算法进行自适应投影寻踪聚类，将聚类标签与图像像素对应，获得图像分割结果。
2. 根据权利要求书1所述的SAR图像分割系统，其中投影寻踪聚类模块包括 初始标签选择子模块，用于对所提图像特征，进行聚类获得初始的标签，并将标签传输给投影方向选择子模块；投影方向选择子模块，根据输入标签计算投影寻踪算法中的线性判别分析投影 指标，并通过免疫克隆算法优化该投影指标，获得最优的投影方向传至子空间聚类 子模块；子空间聚类子模块，将图像特征映射到最优投影方向，获得最优子空间，对最 优子空间进行聚类，获得聚类标签，将聚类标签返回至初始标签选择子模块，迭代 数次，获得最终的聚类标签，将最终聚类标签与图像像素对应，获得最终图像分割 结果。
3. —种基于免疫克隆和投影寻踪的SAR图像分割方法，包括如下步骤(1) 对输入的图像提取灰度共生矩阵、小波、brushlet和contourlet特征，将提 取的图像特征定义为J^，其中，！'表示样本的个数，在图像中表示为第/个像素点， y表示样本的维数，在图像中表示为像素点的第y维特征；(2) 对提取的图像特征J^进行归一化，设定标签迭代次数g;(3) 将归一化后的图像特征《.，随机投影到一个低维子空间f/，并对该子空 间t/进行K均值聚类，获得一个初始标签/^(4) 根据获得的初始标签/f，计算投影寻踪的线性判别分析投影指标，并通过 免疫克隆算法优化该投影指标，获得最优的投影方向；(5) 将归一化后的图像特征《，投影到最优的投影方向，获得最优的子空间(6) 对子空间i/'进行K均值聚类，获得新的标签H';(7) 将新的标签/f'作为初始标签//,返回步骤(4)进行迭代，经过数次迭代 后，若获得新标签/f'在迭代过程中达到收敛，输出收敛的标签i/'，若在迭代g次后 没有收敛，则输出最后获得的标签//';(8) 将收敛的标签i/'或最后获得的标签/Z'与输入图像的像素对应，输出图像 分割结果。
4.根据权利要求书3所述的SAR图像分割方法，其中步骤(4)所述的通过免疫克 隆算法优化该投影指标，按如下过程进行4a)设定种群大小W，进化代数w,克隆规模 ，变异概率;^; 4b)根据种群大小iV,随机产生初始种群J; 4c)计算初始种群^的个体适应度；4d)根据克隆规模 ，对初始种群进行克隆，获得克隆后的种群^;4e)对克隆后种群4，根据变异概率&，进行变异操作，获得变异后的种群4f)计算克隆变异后的种群4的适应度函数值，将种群4中个体适应度与初始 种群4中个体适应度进行对比，依次选出适应度较大的个体，获得新的种群J'并输出最大适应度函数值，将新的种群乂作为初始种群」，返回步骤4c);4g)在进化gw代中，若每一代输出的最大适应度函数值收敛，则终止迭代。
全文摘要
本发明公开了一种基于免疫克隆和投影寻踪的SAR图像分割系统及分割方法。该系统包括图像特征提取模块、初始标签选择子模块、投影方向选择子模块和子空间聚类子模块。图像特征提取模块，提取输入图像的灰度共生矩阵、小波、brushlet和contourlet特征；初始标签选择子模块，对图像特征进行聚类获得初始标签，传输给投影方向选择子模块计算线性判别分析投影指标，并获得最优投影方向；子空间聚类子模块，将图像特征投影到最优投影方向，获得最优子空间，对最优子空间进行聚类，获得聚类标签，并返回至初始标签选择子模块迭代获得最终聚类标签，并与图像像素对应，获得最终图像分割结果。本发明具有分割精度高的优点，可用于军事侦察手段以及民用和工业领域。
文档编号G06T7/00GK101667292SQ20091002405
公开日2010年3月10日 申请日期2009年9月25日 优先权日2009年9月25日
发明者静 冯, 吴建设, 慕彩红, 朱虎明, 杨淑媛, 焦李成, 王宇琴, 缑水平, 桦 钟 申请人:西安电子科技大学