专利名称：基于局部Gamma拟合的活动轮廓SAR图像分割方法
技术领域：
本发明属于图像处理领域，涉及一种图像分割方法，特别涉及一种SAR图像的 活动轮廓图像分割方法。
背景技术：
合成孔径雷达(SAR)以其全天时、全天候的高分辨率数据获取能力成为军事 及民用领域普遍采用的对地观测技术。图像分割是SAR图像自动解译的基础性工作，受 到了广泛的关注。由于SAR系统采用相干成像处理，致使SAR图像中存在大量相干斑噪 声，这也是SAR图像精确分割的重要障碍。立足于消除斑点噪声的影响，传统SAR分割 方法需求通过相干斑预处理，但由于预处理在抑制相干斑的同时，不可避免地损失了潜在 边界等细节信息，致使分割效果不佳。目前，基于曲线演化理论和水平集方法的几何活动轮廓模型(Active Contour Model, ACM)受到了众多研究者的关注，SAR图像的ACM分割研究也方兴未艾。ACM 实现图像分割的基本思想是将演化曲线(活动轮廓)隐含地表示为更高维曲面的零水平 集，水平集在演化方程控制下进行演化，直至图像的目标边界。ACM有许多优点曲线演 化过程始终处于固定格网，数值实现简单；能够自动处理边界的拓扑变化；提供光滑、 闭合的高精度分割曲线，方便后续的图像解译工作；特别是区域ACM能够利用图像全局 信息定义能量泛函，对斑点噪声具有较好的抑制作用，从而避免了复杂的相干斑预处理过 程。然而,现有SAR图像的ACM分割方法均基于图像区域内目标地物的强度统计， 其分割精度取决于统计分布对图像数据的整体拟合程度，而忽略了边界局部结构的有效表 征，从而导致分割结果不理想，特别是在处理深度凹陷等复杂边界情况时存在较大局限。
发明内容
本发明针对现有技术的不足，提出一种基于局部Gamma拟合的活动轮廓SAR图 像分割方法。该方法首先利用Gamma分布表征SAR图像中每个像素一定邻域内的局部 结构，以极大似然作为该邻域的区域分离准则并以此定义局部拟合能量项；然后通过全 局积分，得到全局优化的能量泛函；最后利用几何活动轮廓模型描述和求解能量最小化 过程，更好地实现SAR图像的精确分割。本发明是通过以下技术方案实现，包含的步骤如下
步骤1将定义域为 eTR2的SAR图像I:Ω -IR划分为个全剖分的同质子区域
I,该同质区域集满足
权利要求
1.基于局部Gamma拟合的活动轮廓SAR图像分割方法，其特征在于该方法包括以下 步骤步骤1将定义域为Ω€ R2的SAR图像/ Ω ^R划分为个全剖分的同质子区域4N，该同质区域集满足UQi = Q, Ω<ηΩ , . .= 0 ；对定义域上的每个像素点的邻域内执行步骤1.1 步骤1.3，实现局部划分步骤1.1根据相干斑噪声模型，目标区域均勻的SAR图像满足Gamma分布，利用Gamma分布拟合邻域Os.内的局部结构，定义区域Ω; Π A的概率模型，对于任意 J^qnq / …、"“MMUβ其中￡是SAR图像的视数，、是区域ΑΠΟ,的强度均值，r是Gamma函数； 步骤1.2利用极大似然作为邻域4的区域分离准则MML{Q\I) = U Π U' ),最大化Ml等价于最小化-Iog(M1)，以-Iog(Afl)定义能量泛函即可表征局部区域 分离；步骤1.3结合用以控制邻域采样范围的权函数％，得到邻域4的局部拟合能量项 E，从而实现局部划分K = —ΣJrjilic, A^-y)loSP-^i(ι(.>,))# ；用水平集函数#表示局部二分拟合能量项E、= —Σviy))^ ，其中均⑷=//⑷，Μ2(φ) = Ι-Η(Φ), η是单位阶跃函数； 步骤2根据步骤1所得的Z点局部能量项，在定义域内对积分得到全局的能量泛函 5 = J Ex (φ, cu, C2^ )ck + aP (+ ν￡ (,其中卜if咖用以纠正水平集函数与符号距离函数的偏差，避免演化过程中的周期性初始化过程，丄⑷=|V#〔F)楸是长度约束项，V是梯度算子；步骤3利用Euler-Lagnmge变分极小化步骤2的能量泛函，得到曲线演化的控制方程ft2是时间变量，式中= 幻是狄拉克函数，r是演化曲线的曲率，
全文摘要
本发明公开了一种基于局部Gamma拟合的活动轮廓SAR图像分割方法。现有的技术SAR图像分割精度不高，本发明首先利用Gamma分布表征原始SAR图像中每个像素一定邻域内的局部结构，以极大似然作为该邻域的区域分离准则，并定义局部拟合能量项；然后在定义域内对所有像素点进行积分，从而得到全局优化的能量泛函；最后采用基于曲线演化理论和水平集方法的几何活动轮廓模型描述和求解能量最小化过程，从而实现SAR图像的有效分割。本发明具有较强的斑点噪声抑制能力，能够实现SAR图像的有效分割，特别是针对深度凹陷等复杂边界仍有较为精确的分割结果。
文档编号G06T7/00GK102024260SQ20101061412
公开日2011年4月20日 申请日期2010年12月30日 优先权日2010年12月30日
发明者吴海燕, 孔丁科, 宋超, 王勋 申请人:浙江工商大学