分割图像的方法及装置的制造方法

分割图像的方法及装置的制造方法
【专利摘要】本发明公开了一种图像分割的方法及装置。该方法包括：初始化水平集函数，根据水平集函数值的正负，将输入的待分割SAR图像分割为两个区域Ω1和Ω2，使用固定大小的局部窗口利用对数累积量方法对图像每个像素点估计其广义Gamma分布的局部参数；利用广义Gamma分布的累积分布函数设计能量函数；使用Kolmogorov?Smirnov检验计算两个区域内所有像素点对应经验累积分布函数的K?S距离，根据K?S距离对应的z来确定zm值；根据估计得到的广义Gamma分布的局部参数和zm值来计算能量函数；通过使代价函数最小化得到水平集函数，对水平集函数离散化得到水平集演化方程，且在演化过程中，根据K?S距离重新确定zm并计算能量函数ε。最终实现快速、准确和有效的SAR图像分割。
【专利说明】
分割图像的方法及装置
技术领域
[0001] 本发明设及图像处理领域，具体地，设及一种分割图像的方法及装置。
【背景技术】
[0002] 合成孔径雷达（英文：Synthetic Aperture Radar;缩写：SAR)是一种主动成像系 统，具有不受光照、天气条件影响，可W全天候、全天时对地观测，及透过地表和植被获取信 息等优点，使得SAR图像在农业、林业、地质、环境、水文、海洋、灾害、测绘与军事领域得到了 广泛地应用，尤其是在传统光学传感器成像困难的领域发挥了不可替代的作用。
[0003] SAR图像分割作为SAR图像处理与解译的基础和关键技术之一，其目的是按照一定 准则把SAR图像划分成若干个满足同质性、惟一性且互不交叠的区域，用分割区域表征图像 感兴趣的地方，实现图像描述方式的转换，为图像后续处理和解译提供重要的依据。然而， 相干成像的特性使SAR图像不可避免地受到相干斑噪声的影响，从而导致SAR图像分割更为 复杂和困难。
[0004] 至今为止，研究人员已经提出了多种不同的SAR图像分割方法，其中基于表面演化 的主动轮廓方法是一种十分简单且高效的方法。对于主动轮廓模型，最早由Kass提出了参 数主动轮廓方法，尽管结合似然比边缘检测已成功用于SAR图像分割，但是其不能处理轮廓 演化过程中的拓扑结构变化且算法稳定性不足。相较而言，水平集方法可W适应拓扑结构 的变化而且算法稳定性较高，特别是几何主动轮廓方法逐渐成为该领域的研究热点。几何 主动轮廓方法大致可分为基于边缘和基于区域的两种方法，前者主要利用图像梯度来终止 目标边界上的轮廓演化，而后者则利用了不同区域内的统计情况来引导轮廓运动的方向。 与基于边缘的几何主动轮廓方法相比，基于区域的主动轮廓方法在图像存在噪声和目标边 缘轮廓不明显的情况下具有更好的分割能力，且对初始轮廓位置具有较好的鲁棒性。
[0005] 然而，运些传统的水平集方法都是在W假设噪声为加性噪声为前提的情况下所提 出的，并不适用于SAR图像分割处理。为了消除相干斑噪声的影响，研究人员开始对基于统 计模型的水平集方法进行研究，用W实现SAR图像的分割。
[0006] 对于基于统计模型的SAR图像水平集分割方法，选取适当的统计模型具有十分重 要的意义，统计模型对SAR图像的拟合性能直接决定了 SAR图像水平集分割方法的性能。现 有方法中选取的统计模型有如Ayed I B[Ayed I B,Mitiche A,Belhadj Z.Multiregion Level-Set Partitioning of Synthetic Aperture Radar Images[J].Pattern Analysis and Machine Intelligence,I趾E IYansactions on,2005,27(5) :793-800.]所采用的 Gamma分布。虽然Gamma分布能较好地描述中低分辨率情况下的SAR图像分布均匀区域的统 计特性，但是随着分辨率的提高，SAR图像相干斑发育不充分，均匀度降低，使得Gamma分布 不能够很好地描述其统计特性。因此，现有基于Gamma分布的分割方法不适用于高分辨率的 SAR图像分割处理。
[0007] Marques R C P[Marques R C P,Medeiros F N,Nobre J S. SAR Image Segmentation Based on Level Set Approach and:蹲Model[J].Pattern Analysis and Machine Intelligence, I邸E IYansactions on,2012,:34(10) :2046-2057.]采用的G日分布 虽然能够描述高分辨率情况下SAR图像的分布均匀和分布不均匀区域的统计特性，但是它 只能用于描述幅度SAR图像，具有一定的局限性。
[000引此外，Ayed I B使用对数似然函数来构造能量函数，而Ma巧Ues R C P直接利用统 计模型的累计分布函数来设计能量函数，但其根据总能量变化来计算区域边界的能量波 动。定义SAR图像中区域边界的能量波动达到最大时对应的像素值为Zm(Zm的取值范围为SAR 图像的最小像素值与最大像素值之间），根据捜寻到的Zm通过能量函数计算各个区域的能 量值，利用计算出的能量值进行水平集演化，最终得到分割结果。运种方法得到的Zm并不是 最优的，从而影响水平集分割结果。

【发明内容】

[0009] 本发明的目的是提供一种分割图像的方法及装置，该方法用于实现快速、准确和 有效的SAR图像分割。该方法采用了广义Gamma分布来描述SAR图像的统计特性，该统计模型 具有简洁的表示形式，且为一个函数簇，瑞利分布、指数分布、化kagami分布、Gamma分布和 韦伯分布等都是其在不同参数下的特例。本发明通过对数累积量方法估计得到SAR图像中 的每个像素点的广义Gamma分布的幕参数V、尺度参数O和形状参数K，其中V声〇,〇>〇, K>〇。 然后利用KolmogoroV-SmirnoV检验化-S检验)确定Zm，结合每个像素点的V、O和KW及确定 的Zm来计算能量函数e。利用上述能量函数e计算参数闭合曲线的代价函数，通过使代价函 数取得最小，得到水平集函数。最后，对水平集函数离散化得到水平集演化方程，进而进行 有限次数的水平集函数的演化。在演化过程中，利用K-S检验重新计算Zm,从而重新计算能 量函数e，使两个区域的能量函数具有最佳判决(或区分)能力。经过水平集函数的迭代演化 最终实现快速、准确和有效的SAR图像分割。
[0010] 为了实现上述目的，本发明第一方面提供一种分割图像的方法，所述方法包括：
[0011] 通过K-S检验确定待分割图像的Zm,所述Zm是通过K-S检验获得的所述待分割图像 中的第一类区域内的像素点和第二类区域内的像素点对应的经验累积分布函数的K-S距 离，进而通过K-S距离对应的Z来确定的，所述第一类区域与所述第二类区域W初始分割曲 线为界线；
[0012] 根据所述Zm,确定所述第一类区域的能量函数为第一函数，并确定所述第二类区 域的能量函数为第二函数；
[0013] 根据所述第一函数和所述第二函数，通过水平集演化确定所述第一类区域与所述 第二类区域W最终分割曲线为界线，所述最终分割曲线与所述初始分割曲线不同。
[0014] 可选地，在所述通过K-S检验确定待分割图像的Zm之前，所述方法还包括：
[0015] 将所述待分割图像中所有像素点按图像坐标X对应到初始水平集函数护X的函数 值，所述初始水平集函数护X为：
[0016]
[0017]其中，说)/化为所述初始分割曲线，Q 1和Q 2分别为W所述初始分割曲线为界线的 所述第一类区域和所述第二类区域。
[001 8]可选地，根据所述Zm，确定所述待分割图像的能量函数，包括：
[0019] 确定广义Gamma分布的累积分布函数为所述待分割图像的能量函数，其中，所述广 义Gamma分布的累积分布函数为：
[0020]
[002。 ma函数，式中X为一个随机变量， 「?为Gamma函数，v、o和K分别为所述待分割图像中每个像素点的广义Gamma分布的幕参 数、尺度参数和形状参数。
[0022] 可选地，所述广义Gamma分布的累积分布函数是根据广义Gamma分布的概率密度函 数得到的，所述广义Gamma分布的概率密度函数用于描述所述待分割图像的统计特性，所述 广义Gamma分布的概率密度函数为：
[0023]
[0024] 其中，Z表示待分割图像的像素点的像素值，取值范围为[0,2S-1]，S为图像量化级 另Ij，且「?为Gamma函数。
[00巧]可选地，所述方法还包括：
[0026] 按照W下步骤确定所述待分割图像中的每个像素点的广义Gamma分布的幕参数V、 尺度参数O和形状参数的K估计量、S和《：
[0027] 步骤一:根据指定像素点所在局部窗口内的所有像素的像素值，确定所述指定像 素点的前=阶对数累积量、^2和％，所述指定像素点位于所述待分割图像上的预定大小 的局部窗口中屯、：
[002引
[0029] 其中，Z为所述待分割图像上的局部窗口内的像素点的像素值，N为所述待分割图 像上的局部窗口内的像素点的总数；
[0030] 步骤二:根据6;和續的大小关系，确定所述房：
[0031] 若比猫 < .句/巧< O.3巧根据f K的单调性，通过二分法求解f /璋=骑1:K /辞SK ,确 定所述沒，:其中，巫〇1，k>1/k+1/k1/2k2，巫
[0032] 若巧/片 > 化?>巧，未Il巧闭合古准曲I瑜巧所沐：
[0033]
[0034]
[0035]
[00；3 引
[0036] 若切/巧 <化2日，对所述局部窗口的大小进行扩展，直至满足巧;/巧' > O.25后，按照所述步骤二确定所述揉;[0037] 步骤S:按照W下公式确定所述#、& ;
[0039] ^ 化Iygamma函数。
[0040] 可选地，所述根据所述待分割图像的能量函数，通过水平集演化确定所述第一类 区域与所述第二类区域W最终分割曲线为界线，包括：
[0041 ]确定所述第一类区域与所述第二类区域在所述水平集演化每完成M次时的界线；
[0042] 根据确定出的界线，通过K-S检验重新确定所述待分割图像的Zm;
[0043] 根据重新确定出的Zm,重新计算所述待分割图像的能量函数；
[0044] 根据重新确定出的所述待分割图像的能量函数，通过水平集演化确定所述第一类 区域与所述第二类区域W所述最终分割曲线为界线。
[0045] 本发明第二方面提供一种分割图像的装置，所述装置包括：
[0046] 第一确定模块，用于通过K-S检验确定待分割图像的Zm,所述Zm是通过K-S检验获得 的所述待分割图像中的第一类区域内的像素点和第二类区域内的像素点对应的经验累积 分布函数的K-S距离，进而通过K-S距离对应的Z来确定的，所述第一类区域与所述第二类区 域W初始分割曲线为界线；
[0047] 第二确定模块，用于根据所述Zm,确定所述第一类区域的能量函数为第一函数，并 确定所述第二类区域的能量函数为第二函数；
[004引第=确定模块，根据所述第一函数和所述第二函数，通过水平集演化确定所述第 一类区域与所述第二类区域W最终分割曲线为界线，所述最终分割曲线与所述初始分割曲 线不同。
[0049] 可选地，所述装置还包括：
[0050] 对应模块，用于在所述通过K-S检验确定待分割图像的Zm之前，将所述待分割图像 中所有像素点按图像坐标X对应到初始水平集函数护X的函数值，所述初始水平集函数护X 为：
[0化1 ]
[005^ 其中，OQ/0X为所述初始分割曲线，Q 1和Q 2分别为W所述初始分割曲线为界线的 所述第一类区域和所述第二类区域。
[0053] 可选地，所述第二确定模块用于：
[0054] 确定广义Gamma分布的累积分布函数为所述待分割图像的能量函数，其中，所述广 义Gamma分布的累积分布函数为：
[0化5]
[0056] ma函数，式中X为一个随机变量，
r ?为Gamma函数，u、〇和K分别为所述待分割图像中每个像素点的广义Gamma分布的幕参 数、尺度参数和形状参数。
[0化7] 可选地，所述装置还包括：
[005引参数确定模块，用于按照W下步骤确定所述待分割图像中的每个像素点的广义 Gamma分布的幕参数V、尺度参数O和形状参数K的估计量矿、.S和《；
[0059] 步骤一:根据指定像素点所在局部窗口内的所有像素的像素值，确定所述指定像 素点的前=阶对数累积量n、^和戶3，所述指定像素点位于所述待分割图像上的预定大小 的局部窗口中屯、：
[0060]
[0061] 其中，Z为所述待分割图像上的局部窗口内的像素点的像素值，N为所述待分割图 像上的局部窗口内的像素点的总数；
[0062] 步骤二:根据诗和诗的大小关系，确定所述為:
[00创若0迎空哼/璋<0.疏，根据沪6的单调性，通过二分法求解K :。皆/璋-却1:,:k/*Uk , 确定所述《，其中，巫01，k>1/k+1/化2,巫02，k>-1/k2-1/k3;
[0064] 若巧/巧> O.3巧，利用闭合式准则确定所述A :
[00 化]
[0066] , 異二或
[0067] 若巧/巧 < 化巧，对所述局部窗口的大小进行扩展，直至满足巧/巧S化巧后，按 照所述步骤二确定所述店；
[0068] 步骤按照W下公式确定所述度、度;
[0069]
[0070] 其中，巫0 t = d log r t/dt为Digamma函数，巫0 n,t = (T+i log r t/dtn+i为 化Iygamma函数。
[0071 ]本发明与现有的技术相比具有W下优点：
[0072] 1、本发明所采用的广义Gamma分布模型具有拟合能力强和灵活性高的特点，能够 有效的描述SAR图像数据，能同时应用于强度图像和幅度图像，并且直接利用模型描述构造 能量泛函，不需要对SAR图像进行预处理，提高了精确性；
[0073] 2、本发明采用对数累积量方法，联合二分法和闭合式准则两种方法在不同条件下 完成SAR图像中的每个像素点的广义Gamma分布的幕参数V、尺度参数O和形状参数K的估计， 在保证参数准确估计的同时也提高了运算效率；
[0074] 3、本发明使用K-S检验，W确定使不同区域内图像像素的经验累积分布函数的差 异达到最大时所对应的像素值Zm,使得两个区域的能量函数具有最佳判决(或区分)能力。 与通过选取当区域边界的能量波动最大时确定的对应的像素值Zm运种方法相比较，本发明 中的方法可W直接根据不同区域内的所有像素值来确定Zm,结果更为准确；
[0075] 经过与现有分割方法的仿真比较，证明了本发明能够完成更加精细准确的分割， 且更具灵活性。
[0076] 本发明的其他特征和优点将在随后的【具体实施方式】部分予W详细说明。
【附图说明】
[0077] 附图是用来提供对本发明的进一步理解，并且构成说明书的一部分，与下面的具 体实施方式一起用于解释本发明，但并不构成对本发明的限制。在附图中：
[0078] 图1是本发明的流程图；
[0079] 图2是作为本发明具体实施例的真实SAR图像；
[0080] 图3和图4分别是作为本发明具体实施例的合成SAR图像；
[0081] 图5和图6分别是本发明与现有基于护分布的方法、现有基于Gamma分布的方法对 图2曰、图化的分割效果对比图；
[0082] 图7和图8分别是本发明与现有基于护分布的方法、现有基于Gamma分布的方法对 图3、图4的分割效果对比图；
[0083] 图9是不同参数下合成SAR图像分割效果评价指标的示意图； 图10为本发明提供的一种分割图像的装置的示意图。
【具体实施方式】
[0084] W下结合附图对本发明的【具体实施方式】进行详细说明。应当理解的是，此处所描 述的【具体实施方式】仅用于说明和解释本发明，并不用于限制本发明。
[0085] 参照附图1，本发明的实现步骤如下：
[0086] 步骤一，根据初始水平集函数护X数值符号的正负，将输入的待分割SAR图像的 整幅图像区域Q分割成两个区域Q1和Q 2;
[0087] 该步骤的具体实现过程为:初始化一个与待分割图像同样大小的矩阵，在矩阵内 部选择一定大小的区域，将所选区域部分的值设置为-1，矩阵剩余区域部分的值设置为1， 根据初始的水平集函数护X的正负，将SAR图像分割成两个区域Q1和Q2。
[0088] 其中，初始水平集函数护X表示为：
[0089]
[0090] 其中，X为所述待分割图像中的任一像素点的图像坐标，且SAR图像的每一个像素 点在矩阵中按坐标X对应都有一个相应值(对应值为-I或I)。Q I和Q 2分别为分割得到的两 个区域，说)/游为区域边界。
[0091] 步骤二，采用广义Gamma分布来描述SAR图像的统计特性，利用对数累积量方法计 算SAR图像中的每个像素点的广义Gamma分布的幕参数V、尺度参数O和形状参数K的估计量， 具体实现步骤如下：
[0092] 2a)广义Gamma分布的概率密度函数表示为：
[0093]
[0094] 其中，Z在式中为R+内的一个随机变量。当具体采用上述概率密度函数来描述SAR 图像时，Z表示SAR图像的像素值，取值范围为[0,2S-1]，S为图像量化级别，且r .为Gamma 函数。
[00M] 2b)考虑具体处理的SAR图像，选择一个适当的固定大小（如3X3、5X5等，单位为 像素)的局部窗口对指定像素点提取得到窗口内SAR图像的像素值，计算像素值的前S阶对 数累积量6、
[0096]
[0097] 其中，Z为当前局部窗口内SAR图像的像素值，N为当前局部窗口内SAR图像的像素 点的个数；
[0098] 2c)根据*-23和^的大小关系，分成如下S种情况计算参数K的估计量房：
[0099] 若〇批; '、O.3巧,根据K的单调性，通过二分法求解.r K二辦=和,U /吗2,,、'得 到估计量淀，其中，巫0 1，K>1/k+1/化2,巫〇2，k>-1/k2-1/k3;
t ,
[0101] 若?<(口。，将当前局部窗口的边长按2个像素单位逐步增加(如初始窗口为 3 X 3，扩展后为5 X 5，若再不满足继续扩展为7 X 7，单位均为像素，W此类推），按步骤化)重 新计算fl、气和弓，直至满足吗/巧含〇.朗，按照步骤扣)计算参数K的估计量在；
[0102] 2d)根据下式来计算参数￥和〇的估计量tK片：
[0103]
[0104] 其中，巫0 t = d Io 邑厂 t/dt 为 Digamma 函数，巫 0 n,t = d"i Io 邑厂 t/dt。"为 化Iygamma函数。
[01 05]步骤=，定义能量函数为广义Gamma分布的累计分布函数：
[0106] 3a)对于SAR图像中同一区域Qi的两个不同像素点P和Q，运两个像素点的累积分 布函数 P(Z《Zm)=F V,〇,K;Zm 之间的关系为 F Vp,〇P,Kp;Zm>F VQ,〇Q,KQ;Zm，其中 P(Z《Zm)中 的Z表示SAR图像的像素值，而Zm为Z取值范围内的某一个值。若像素点P和Q属于不同的区 域，则它们描述了不同区域的特性。基于累积分布函数的单一性，可W用不同的累积分布函 数表述不同的区域。因此，定义能量函数如下：
[0107] e(v,〇,K;Zm)=F(v,〇,K;Zm);
[0108] 3b)广义Gamma分布的累积分布函数表达式如下：
[0109]
[0110] 其中，Zm为采用广义Gamma分布描述的SAR图像像素值取值范围内的某一个值，运 一个值是SAR图像的不同区域的区分度达到最大时对应的像素值。巧
为 不完全Gamma函数，式中X为一个随机变量，「?为Gamma函数。
[0111 ] 步骤四，对两个区域Q 1和Q 2内所有像素值，矛Il用Kolmogorov-Smirnov检验计算两 个区域所有像素点对应经验累积分布函数的K-S距离D，根据K-S距离D对应的Z来确定Zm的 值。
[0112] 步骤五，根据步骤S所述的能量函数ev，〇，K;Zm，结合步骤二得到的广义Gamma分 布的局部参数V、0和K，W及步骤四得到的Zm值，计算两个区域Ql和Q 2的能量函数和 气。
[0113] 步骤六，结合步骤五所述的不同区域对应的能量函数％,和当前时刻的水平集函 数4，W及分割得到的两个区域Q 1和Q 2进行水平集辉化.且化爐讨讨賴々n下：
[0114] 6a)计算线性代价函数fc= I化寸2|，其中的
为区域Qi的平均 能量值
乐别为区域Qi的总能量和面积，而式中的Z是积分中的一个随机变 量，具体设及SAR图像时，表示SAR图像中区域Q i包含坐标位置所对应的像素点；
[0115] 6b)采用水平集模型而/流二一雖./试;^，使得参数闭合曲线則勺代价函数最小化， 得到水平集演化模型d式/航：
[0116]
[0117] g合曲线#正交的单位矢量；
[011引 也函数恥，其中的5e4= G A ? E2+护为近似狄 拉克S函数，4表示SAR图像像素所对应的水平集函数，而G为一个在r内的正则化常数，取 值为1;
[0119] 6d)设定演化过程的迭代步长At,根据当前(演化)时刻水平集函数护和能量函数 ek，通过离散化得到的演化方程批+1=护+A t .恥，计算下一(演化)时刻水平集函数批+1;
[0120] 6e)利用高斯滤波器对步骤6d)求得的下一(演化)时刻水平集函数护+1进行正则化 处理；
[0121] 6f)在上述步骤的演化过程中，当水平集演化方程每演化M次时(演化M次指V^i = 护+At ?恥迭代计算M次，即水平集函数演化M次），根据当前(演化)时刻水平集函数分割的 两个区域，按步骤四重新计算Zm,按步骤五计算能量函数S心来替代ek。
[0122] 步骤屯，判断是否收敛。
[0123] 本步骤的主要目的是开始进行水平集演化后，计算当前时刻化)和上一时刻化-1) 的前L次演化线性代价函数的均值的差^
ft为某一时刻下前L次演化的 线性代价函数的均值，L为一个有限大的整数(例如5、10或50)，此处的当前时刻化)和上一 时刻化-1)表示水平集某两个相邻的演化(如k表示第100次演化，而k-1表示第99次演化）。 根据Al来判断当前演化是否已达到迭代终止条件，如果当前演化未达来判到迭代终止条 件，则转到步骤五，用鸣和確替代哨-1和邱-1继续迭代;如果当前演化达到迭代终止条件， 则停止迭代，得到的两个区域Qf和即为最终的分割结果。运里的迭代终止条件包括达到 最大迭代次数(如最大迭代次数为1000次，根据具体处理的SAR图像设置实际的最大迭代次 数）和满田下古
[0124]
[01巧]其中，J为某一（演化)时刻前L次演化过程中的线性代价函数均值，A C为收敛阔 值。
[0126] 本发明的效果可W通过下面的仿真进一步说明：
[0127] 1.仿真平台
[012引 硬件平台为：Intel(R)Core(TM)i7-4700MQ CPU@2.40GHz，8.0 GB RAM。
[0129] 软件平台为:Mathworks Matlab R2014b(8.4.0.150421)。
[0130] 2.仿真内容与结果分析
[0131] 仿真中具体实施例的图像如图2、图3和图4所示，用于本发明的测试。其中，图2为 两幅真实SAR图像；图3和图4为不同参数下的合成SAR图像，并且图3a和图4a为幅度图像，图 3b和图4b为强度图像。本仿真中的合成图像数据均由K分布(本领域常用概率分布模型)产 生。
[0132] 应用本发明和现有的基于护分布和Gamma分布的分割方法分别对图2、图3和图4进 行分割实验，其分割结果如图5、图6、图7和图8所示，其中：图5曰、图化和图5c分别为本发明、 基于GD分布和Gamma分布的分割方法对图2a的分割结果；图6a、图化和图6C分别为本发明、 基于沪分布和Gamma分布的分割方法对图化的分割结果；图7a和图7b分别为本发明和基于沪 分布的分割方法对图3a的分割结果，图7c和图7d分别为本发明和基于Gamma分布的分割方 法对图3b的分割结果；图8a和图8b分别为本发明和基于护分布的分割方法对图4a的分割结 果，图8c和图8d分别为本发明和基于Gamma分布的分割方法对图4b的分割结果。
[0133]由图加、图5c、图化和图6c可W看出，现有的基于gD分布和Gamma分布的分割方法 检测出的边缘准确率不高。由图化可W看出，现有的基于护分布的分割方法检测出了较多 的虚假边缘，而由图化可知该方法检测不出一些边缘信息。由图5c可W看出，现有的基于 Gamma分布的分割方法检测出的边缘不够准确，存在一些无用的边缘，而由图6c同样可W看 出该方法检测不出一些边缘信息。相对于图加、图5c、图6b和图6c，由于本发明采用的广义 Gamma分布具有更强的拟合能力，从图5a和图6a可W看出，本发明得到的分割结果更为准 确，且分割曲线较为平滑。
[0134] 由图7和图8也可W看出，对于根据K分布产生的单视和多视合成SAR图像，本发明 具有更高的灵活性，对幅度图像和强度图像都能得到较好的分割结果。从图7a与图7b、图7c 与图7d、图8a与图8bW及图8c与图8d的对比可W看到，现有基于沪分布和Gamma分布的分割 方法得到的分割结果准确率不足。
[0135] 对于根据K分布产生的合成SAR图像，仿真中在不同等效视数和不同K分布形状参 数和尺度参数下进行了多组实验。其中，等效视数是多视SAR图像的一个属性参数，用来衡 量一幅图像斑点噪声相对强度的一种指标。本仿真中WKappa系数作为评价指标，Kappa系 数是一个衡量评估者间一致性的常用统计方法，取值范围为[-1，1 ]。该系数越高表示相应 的分割结果越准确，反之亦然。表1给出了不同参数下合成SAR图像分割效果评价指标。表中 "背景/目标"表示不同区域相应的K分布形状参数和尺度参数，其中横向为尺度参数，纵向 为形状参数;采用方法中"GO"和％amma"分别表示现有基于6<^分布和6曰111111曰分布的分割方 法，而"GGD-A"和"GGD-r分别表示本发明方法用于幅度图像和强度图像分割的情况。
[0136] 通过表1的分割效果评价指标可W看出，对于不同等效视数、形状和尺度参数下的 合成SAR图像，本发明用于幅度图像和强度图像的分割都能得到较好且稳定的结果。从表中 结果可W看出，现有基于沪分布和Gamma分布的分割方法对不同参数下合成SAR图像的分割 准确率不稳定，且均不如本发明的分割效果好。
[0137] 基于同一发明构思，本发明还提供一种分割图像的装置。请参考图10,图10为本发 明提供的一种分割图像的装置的示意图。该装置100包括：
[013引第一确定模块101，用于通过K-S检验确定待分割图像的Zm,所述Zm是通过K-S检验 获得的所述待分割图像中的第一类区域内的像素点和第二类区域内的像素点对应的经验 累积分布函数的K-S距离，进而通过K-S距离对应的Z来确定的，所述第一类区域与所述第二 类区域W初始分割曲线为界线；
[0139] 第二确定模块102,用于根据所述Zm,确定所述第一类区域的能量函数为第一函 数，并确定所述第二类区域的能量函数为第二函数；
[0140] 第=确定模块103,根据所述第一函数和所述第二函数，通过水平集演化确定所述 第一类区域与所述第二类区域W最终分割曲线为界线，所述最终分割曲线与所述初始分割 曲线不同。
[0141] 可选地，所述装置还包括：
[0142] 对应模块，用于在所述通过K-S检验确定待分割图像的Zm之前，将所述待分割图像 中所有像素点按图像坐标X对应到初始水平集函数护X的函数值，所述初始水平集函数护X 为：
[0143]
[0144] 其中，饼)/'淑为所述初始分割曲线，Qi和分别为W所述初始分割曲线为界线的 所述第一类区域和所述第二类区域。可选地，所述第二确定模块102用于：
[0146]
[0145] 确定广义Gamma分布的累积分布函数为所述待分割图像的能量函数，其中，所述广 义Gamma分布的累积分布函数为：
[0147] ma函数，式中X为一个随机变量， r .为Gamma函数，u、〇和K分别为所述待分割图像中每个像素点的广义Gamma分布的幕参 数、尺度参数和形状参数。
[0148] 可选地，所述装置还包括：
[0149] 参数确定模块，用于按照W下步骤确定所述待分割图像中的每个像素点的广义 Gamma分布的幕参数V、尺度参数O和形状参数K的估计量6、和该；
[0150] 步骤一:根据指定像素点所在局部窗口内的所有像素的像素值，确定所述指定像 素点的前=阶对数累积量r'2和rl?，所述指定像素点位于所述待分割图像上的预定大小 的局部窗口中屯、：
[0151]
[0152] 其中，Z为所述待分割图像上的局部窗口内的像素点的像素值，N为所述待分割图 像上的局部窗口内的像素点的总数；
[0153] 步骤二:根据璋和/"'32的大小关系，确定所述A :
[0154] 若0'进空弓轉<0巧5 ,根据沪K的单调性，通过二分法求解f K ;;;巧/S 丕I U/本旨2非，确 定所述这;，其中，巫0 1，K>1/k+1/化2,巫〇2,k^j-1/k2-1/k3;
[0155] 若巧/巧 > 化:仍，利用闭合式准则确定所述沒；
[0156]
[0157]
[015 引
[0159] 若巧V弓<〇.2日，对所述局部窗口的大小进行扩展，直至满足殘/每 >化踞后，按 照所述步骤二确定所述店；
[0160] 步骤S:按照W下公式确定所述令、S ;
[0161]
[0162] 其中，巫 〇t = d log r t/dt 为 Digamma 函数，巫 0 n,t = d"i log r t/dt。"为 化Iygamma函数。
[0163] 关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法 的实施例中进行了详细描述，此处将不做详细阐述说明。
[0164] W上结合附图详细描述了本发明的优选实施方式，但是，本发明并不限于上述实 施方式中的具体细节，在本发明的技术构思范围内，可W对本发明的技术方案进行多种简 单变型，运些简单变型均属于本发明的保护范围。
[0165] 另外需要说明的是，在上述【具体实施方式】中所描述的各个具体技术特征，在不矛 盾的情况下，可W通过任何合适的方式进行组合，为了避免不必要的重复，本发明对各种可 能的组合方式不再另行说明。
[0166] 此外，本发明的各种不同的实施方式之间也可W进行任意组合，只要其不违背本 发明的思想，其同样应当视为本发明所公开的内容。
【主权项】
1. 一种分割图像的方法，其特征在于，所述方法包括： 通过Kolmogorov-Smirnov检验确定待分割图像的Zm，所述Zm是通过κ-s检验获得的所述 待分割图像中的第一类区域内的像素点和第二类区域内的像素点对应经验累积分布函数 的Κ-S距离，进而通过Κ-S距离对应的z来确定的，所述第一类区域与所述第二类区域以初始 分割曲线为界线； 根据所述^，确定所述第一类区域的能量函数为第一函数，并确定所述第二类区域的能 量函数为第二函数； 根据所述第一函数和所述第二函数，通过水平集演化确定所述第一类区域与所述第二 类区域以最终分割曲线为界线，所述最终分割曲线与所述初始分割曲线不同。2. 根据权利要求1所述的方法，其特征在于，在所述通过Kolmogorov-Smirnov检验确定 待分割图像的之前，所述方法还包括： 将所述待分割图像中所有像素点按图像坐标X对应到初始水平集函数的函数值，所 述初始水平集函数fx为：其中，8Ω/3Χ为所述初始分割曲线，Ω潮Ω 2分别为以所述初始分割曲线为界线的所述 第一类区域和所述第二类区域。3. 根据权利要求1所述的方法，其特征在于，所述根据所述^，确定所述待分割图像的能 量函数，包括： 确定广义Gamma分布的累积分布函数为所述待分割图像的能量函数，其中，所述广义 Ga_a*布的累积分布函数为：其中，P f 为不完全Gamma函数，式中X为一个随机变量，Γ ·为 { a 0 Gamma函数，v、σ和κ分别为所述待分割图像中每个像素点的广义Gamma分布的幂参数、尺度 参数和形状参数。4. 根据权利要求3所述的方法，其特征在于，所述广义Gamma分布的累积分布函数是根 据广义Gamma分布的概率密度函数得到的，所述广义Gamma分布的概率密度函数用于描述所 述待分割图像的统计特性，所述广义Gamma分布的概率密度函数为：其中，z表示待分割图像的像素点的像素值，取值范围为[0,2s-l]，S为图像量化级别，且 Γ ·为Gamma函数。5. 根据权利要求3或4所述的方法，其特征在于，所述方法还包括： 按照以下步骤确定所述待分割图像中的每个像素点的广义Gamma分布的幂参数v、尺度 参数σ和形状参数κ的估计量￡>_、々和泠 步骤一:根据指定像素点所在局部窗口内的所有像素的像素值，确定所述指定像素点 的前三阶对数累积量A、4和&，所述指定像素点位于所述待分割图像上的预定大小的局 部窗口中心：其中，z为所述待分割图像上的局部窗口内的像素点的像素值，N为所述待分割图像上 的局部窗口内的像素点的总数； 步骤二:根据If和續的大小关系，确定所述总： 5^0.25 < r2!/f5- < 0.375 /φο * 确定所述，?..，其中，Φ〇 1，κ~1/κ+1/2κ2，Φ〇 2，κ~-1/κ2-1/κ3; 若2 0.375，利用闭合式准则确定所述応：若g/ff <0.25,对所述局部窗口的大小进行扩展，直至满足2 0.25后，按照所 述步骤二确定所述/? ; 步骤三:按照以下公式确定所述沒;其中，Φ〇 t = d logrt/dt 为Digamma 函数，Φ〇 n，t = dn+1logrt/dtn+1 为Polygamma 函 数。6. 根据权利要求1所述的方法，其特征在于，根据所述待分割图像的能量函数，通过水 平集演化确定所述第一类区域与所述第二类区域以最终分割曲线为界线，包括： 确定所述第一类区域与所述第二类区域在所述水平集演化每完成Μ次时的界线； 根据确定出的界线，通过Kolmogorov-Smirnov检验重新确定所述待分割图像的zm; 根据重新确定出的zm，重新确定所述待分割图像的能量函数； 根据重新确定出的所述待分割图像的能量函数，通过水平集演化确定所述第一类区域 与所述第二类区域以所述最终分割曲线为界线。7. -种分割图像的装置，其特征在于，所述装置包括： 第一确定模块，用于通过Kolmogorov-Smirnov检验确定待分割图像的zm，所述zm是通过 K-S检验获得的所述待分割图像中的第一类区域内的像素点和第二类区域内的像素点对应 的经验累积分布函数的K-S距离，进而通过K-S距离对应的z来确定的，所述第一类区域与所 述第二类区域以初始分割曲线为界线； 第二确定模块，用于根据所述^，确定所述第一类区域的能量函数为第一函数，并确定 所述第二类区域的能量函数为第二函数； 第三确定模块，根据所述第一函数和所述第二函数，通过水平集演化确定所述第一类 区域与所述第二类区域以最终分割曲线为界线，所述最终分割曲线与所述初始分割曲线不 同。8. 根据权利要求7所述的装置，其特征在于，所述装置还包括： 对应模块，用于在所述通过Ko lmogorov-Smirnov检验确定待分割图像的zm之前，将所述 待分割图像中所有像素点按图像坐标X对应到初始水平集函数的函数值，所述初始水平 集函数fx为：其中肩为所述初始分割曲线，ω#ρω2分别为以所述初始分割曲线为界线的所述 第一类区域和所述第二类区域。9. 根据权利要求7所述的装置，其特征在于，所述第二确定模块用于： 确定广义Gamma分布的累积分布函数为所述待分割图像的能量函数，其中，所述广义 Ga_a*布的累积分布函数为：其中％不完全Gamma函数，式中X为一个随机变量，Γ ·为 Gamma函数，v、σ和κ分别为所述待分割图像中每个像素点的广义Gamma分布的幂参数、尺度 参数和形状参数。10. 根据权利要求9所述的装置，其特征在于，所述装置还包括： 参数确定模块，用于按照以下步骤确定所述待分割图像中的每个像素点的广义Gamma 分布的幂参数V、尺度参数σ和形状参数κ的估计量、$和澆 步骤一:根据指定像素点所在局部窗口内的所有像素的像素值，确定所述指定像素点 的前三阶对数累积量4、匕和&，所述指定像素点位于所述待分割图像上的预定大小的局 部窗口中心：其中，ζ为所述待分割图像上的局部窗口内的像素点的像素值，Ν为所述待分割图像上 的局部窗口内的像素点的总数； 步骤二:根据巧和續的大小关系，确定所述总： 若0.25?;@Α':卜;0.375，根据JF ?的单调性，通过二分法求解f κ =亡"?=封1，《 /辞2,《，确 定所述总，其中，Φ〇 1，κ~1/κ+1/2κ2，Φ0 2，κ~-1/κ2-1/κ3; 若巧斤I 2 0.373，利用闭合式准则确定所述応：若?t/校< 0.25 ,对所述局部窗口的大小进行扩展，直至满足2 0.25后，按照所 述步骤二确定所述总I 步骤三:按照以下公式确定所述#、士:其中，Φ〇 t = d log rt/dt为Digamma函数，Φ〇 n，t = dn+1logr t/dtn+1 为Polygamma函 数。
【文档编号】G06T7/00GK105957050SQ201610237599
【公开日】2016年9月21日
【申请日】2016年4月15日
【发明人】李恒超, 邱允亮, 李静静
【申请人】西南交通大学