基于广义gamma分布的SAR图像超像素生成方法
【技术领域】
[0001] 本发明涉及成像雷达图像分割技术领域，更具体地说，涉及一种SAR(Synthetic Aperture Radar，合成孔径雷达）图像超像素生成方法。
【背景技术】
[0002] SAR是一种主动微波成像传感器，具有全天时、全天候成像能力，这使得SAR图像 在很多军事和民用领域中得到了广泛应用。SAR图像分割是SAR图像自动解译的重要内容， 是大量其他SAR图像解译算法包括目标检测与识别、地物分类等的重要环节。近年来，随着 SAR成像技术的发展，SAR图像的分辨率不断提高。与传统的中低分辨率SAR图像相比，高 分辨率SAR图像能够获取地物和目标更为丰富的信息，如何有效地挖掘这些有效信息是当 前SAR图像解译的重要研究内容。SAR图像超像素生成是一种通过对图像像素进行聚类，得 到一系列大小相近的具有特定意义的图像区域的方法。SAR图像超像素生成本质上是一种 图像分割方法，伴随着高分辨率SAR图像的出现而兴起，在基于区域的SAR图像解译方法研 究和应用系统构建中具有广泛的应用前景。
[0003] 当前，典型的图像超像素生成方法有基于图论的Normalized-Cut方法、凝聚式聚 类（agglomerative clustering)方法、Turbopixel方法和SLIC(Simple Linear Iterative Clustering，简单线性迭代聚类）方法等，其研究主要针对光学图像开展。与其他方法相 比，SLIC方法在边界拟合以及效率上都具有更优的性能。
[0004] SLIC方法主要包括两个步骤：（1)利用局部k-均值算法对图像像素聚类；（2)利 用连通分量算法消除孤立的图像小区域。SLIC方法在聚类过程中，利用了各像素到聚类 中心（即图像区域的均值）的欧氏距离进行聚类。然而，受SAR图像固有相干斑噪声以 及SAR图像幅度值潜在的大动态变化特征的影响，SLIC方法难以生成有效的SAR图像超 像素。另外，连通分量算法仅通过最近距离原则消除孤立小区域，没有充分利用图像本身 的特征，难以处理步骤（1)生成的初始超像素的边界与实际地物边界偏离的问题。由于 F1DF(Probability Density Function,概率密度函数）比图像区域的均值能够更全面地描 述SAR图像区域的统计特性，而广义ga_a分布是一种用于描述SAR图像区域统计特性的 重要TOF，因此，研究基于广义gamma分布的SAR图像超像素生成方法具有重要的理论意义 和应用价值。
【发明内容】
[0005] 本发明为了有效解决SAR图像超像素生成问题，提供了一种基于广义gamma分布 的SAR图像超像素生成方法。本方法能够生成与实际地物的边界吻合较好的超像素，可以 为基于区域的SAR图像解译方法研究和应用系统构建提供重要支撑。
[0006] 本发明的技术方案是：一种基于广义gamma分布的SAR图像超像素生成方法，将 SAR图像划分成若干区域，其特征在于，包括下述步骤：
[0007] 利用广义gamma分布描述SAR图像的区域IW ;结合SAR图像的区域PDF和空间距 离信息构造局部聚类准则，对SAR图像的每一个像素进行聚类，得到SAR图像的若干区域；
[0008] 利用基于局部贝叶斯准则的边界演化方法消除SAR图像若干区域中的孤立小区 域，则剩余的每个区域对应一个SAR图像超像素。
[0009] 特别地，局部聚类准则是基于待分类像素和待聚类中心的幅度相似度和距离相似 度的加权和进行判定。
[0010] 本发明的有益效果是：结合SAR图像的区域PDF和空间距离信息构造局部聚类 准则能够更全面地利用区域的统计信息，生成均匀性更好的初始图像区域。利用基于局部 贝叶斯准则的边界演化方法能够更好地利用SAR图像的上下文信息，不仅能消除孤立小区 域，还能使得生成的超像素的边界与实际地物的边界更加吻合。
【附图说明】
[0011] 图1为本发明的SAR图像超像素生成方法流程图；
[0012] 图2是利用实测图像进行实验的结果，其中（a)原始SAR图像；(b)采用标准SLIC 的结果图；（c)是利用本发明的结果图；
[0013] 图3是对应实测图像中区域1的实验结果，其中（a)原始SAR图像的区域l;(b) 采用标准SLIC的结果图；（c)是利用本发明的结果图；
[0014] 图4是对应实测图像中区域2的实验结果，其中（a)原始SAR图像的区域2;(b) 采用标准SLIC的结果图；（c)是利用本发明的结果图；
【具体实施方式】
[0015] 下面结合附图对本发明提供的SAR图像超像素生成方法进行详细说明。
[0016] 图1是本发明的SAR图像超像素生成方法流程图。该流程的第一步为算法初始 化，包括SAR图像标准化、SAR图像区域初始划分、SAR图像像素标签的初始化、SAR图像聚 类中心初始化以及相似度矩阵初始化。第二步为基于SAR图像的区域PDF和空间距离信息 的SAR图像局部聚类。对应每一个聚类中心，计算其搜索范围内的各像素与其之间的综合 相似度，利用它对相似度矩阵进行更新，进而更新相应像素的标签。第三步为基于局部贝叶 斯准则的SAR图像区域边界演化。由于受到相干斑噪声等因素的影响，第二步得到的SAR 图像区域包含很多孤立的小区域，并且SAR图像区域的边界可能与实际地物的边界存在一 定的偏差，通过采用基于局部贝叶斯准则的边界演化处理，不仅可以消除孤立的小区域，还 使得生成的超像素的边界与地物的边界更加吻合。
[0017] 下面具体说明上述各个步骤：
[0018] 第一步：算法初始化。
[0019] 步骤①，SAR图像标准化。设给定一幅原始SAR图像I。的幅度值为 I。(X，y)，1彡X彡M, 1彡y彡N，其中（X，y)表示该SAR图像的像素坐标值。对原 始SAR图像I。进行标准化可得标准化SAR图像I i,其幅度值为I i (X，y) = I。(X，y) / μ0，1彡X彡M，1彡y彡N，其中内为原始SAR图像I。的均值。标准 M * 7V ν_[· 化的目的是为了压缩图像的动态范围，便于构造合理的像素聚类准则。
[0020] 步骤②，SAR图像区域初始划分。将标准化SAR图像I1划分为 一系列相邻的大小均为LsXLs像素的图像区域R D R2,. . .，Rk，. . .，Rns，即 I1 Ui?2 U...U & U …= ，其中 JVs =M _ /V /€ 为区域总数，为计算方便， 假设M和N均为1^的整数倍，并且L洒取值根据需要确定。
[0021] 步骤③，SAR图像像素标签初始化。对由步骤②划分得到的每个图像区域中的像 素赋予一个与图像区域序号相同的标签，得到标签图像L，即：
[0022] L： = |L.(x, = k,￥(x, _y) 6.?.,. ?<χ< Μ,Ι < j? <
[0023] 步骤④，SAR图像聚类中心初始化。计算标准化SAR图像I1的梯度模图像G，其每 一个像素 G (X，y)由下式进行计算：
[0024]
[0025] 然后将图像区域&中具有最小梯度模的像素作为图像区域1的初始聚类中心Ck， 即：若ivyirgminGho')，则令〇 j ()。
[0026] 步骤⑤，相似度矩阵初始化。初始化相似度矩阵S为M行N列的全零矩阵，即令S ={S (X，y) I S (X，y) = 0, 1彡X彡M, 1彡y彡N}。相似度矩阵记录标准化SAR图像I1中的 每一个像素点与某个聚类中心的相似度。
[0027] 第二步：基于图像区域的PDF和空间距离信息的SAR图像局部聚类。
[0028] 步骤①，估计标准化SAR图像^中各图像区域的TOF。设图像各区域的PDF为广
义gamma分布，第k个图像区域&的PDF的形状参数、能量参数和尺度参数为(6,成、Pi, 卜 按照下述卞'式彳+管.
[0029]
[0030]
[0031]
α = I. 5 ; DiP为％：!:的指示矩阵 的第W行第列元素。
[0055] 步骤④，重复执行步骤③，直到执行次数达到设定的最大次数AC，然后将原始图 像Ic中标签L(x，y) = k的像素点划归到图像区域R k, k = 1，2,. . .，Ns,则图像区域Rk即为 原始图像Ic的第k个超像素。
[0056] 图2(a)给出了一幅大小为300X450像素的原始SAR图像，图2(b)和图2(c)分 别给出了采用标准SLIC和本发明方法生成的SAR图像超像素结果，其中本发明方法中参数 设置为 Ls= 15、W = 0.6、iV￡|x =10、=10 和C =10，α =1.5。从图中可以看出， 标准SLIC方法生成的超像素的边界能够较好吻合对比度较大的地物的边界，但是对于对 比度较小的地物的边界吻合较差。相比之下，本发明方法生成的超像素的边界不仅与对比 度大的地物的边界吻合很好，还与对比度相对较小的地物的边界吻合较好。
[0057] 为了更好地观察标准SLIC方法和本发明方法生成的SAR图像超像素结果，图3和 图4分别给出了利用两种方法对包含不同地物的典型区域（即分别对应图2(a)中左下角 矩形框所示的区域1和右下角矩形框所示的区域2)的处理结果。图3(a)给出了与区域1 对应的原始SAR图像，图3 (b)和图3 (c)分别给出了采用标准SLIC和本发明方法生成的与 区域1对应的超像素结果，其中曲线表示超像素的边界；图4(a)给出了与区域2对应的原 始SAR图像，图4 (b)和图4 (c)分别给出了采用标准SLIC和本发明方法生成的与区域2对 应的原始SAR图像超像素结果，其中曲线表示超像素的边界。
[0058] 从图3和图4可以清楚地看出，与标准SLIC方法相比，本发明方法能够生成与实 际地物边界更加吻合的超像素，典型的如对应图3(a)中左侧中部椭圆所示的区域A和右下 角椭圆所示区域B以及图4(a)中底部中间椭圆所示的区域C的超像素生成结果，从而验证 了本发明方法的有效性。
【主权项】
1. 一种基于广义gamma分布的SAR图像超像素生成方法，SAR是指合成孔径雷达，将 SAR图像划分成若干区域，其特征在于，还包括下述步骤： 利用广义gamma分布描述SAR图像每个区域的区域H)F，PDF是指概率密度函数；结合 SAR图像的区域PDF和空间距离信息构造局部聚类准则，对SAR图像的每一个像素进行聚 类，更新SAR图像的区域划分；利用基于局部贝叶斯准则的边界演化方法消除SAR图像若干 区域中的孤立小区域，则剩余的每个区域对应一个SAR图像超像素。2. 根据权利要求1所述的基于广义gamma分布的SAR图像超像素生成方法，其特征在 于，局部聚类准则是基于待分类像素和待聚类中心的幅度相似度和距离相似度的加权和进 行判定。
【专利摘要】本发明提供一种基于广义gamma分布的SAR图像超像素生成方法。技术方案是：将SAR图像划分成若干区域，利用广义gamma分布描述SAR图像每个区域的区域PDF，PDF是指概率密度函数；结合SAR图像的区域PDF和空间距离信息构造局部聚类准则，对SAR图像的每一个像素进行聚类，更新SAR图像的区域划分；利用基于局部贝叶斯准则的边界演化方法消除SAR图像若干区域中的孤立小区域，则剩余的每个区域对应一个SAR图像超像素。本发明能够更好地利用SAR图像的上下文信息，不仅能消除孤立小区域，还能使得生成的超像素的边界与实际地物的边界更加吻合。
【IPC分类】G06T7/00
【公开号】CN105005987
【申请号】CN201510348973
【发明人】邹焕新, 秦先祥, 周石琳, 计科峰, 雷琳
【申请人】中国人民解放军国防科学技术大学
【公开日】2015年10月28日
【申请日】2015年6月23日