基于局部Gamma拟合的活动轮廓SAR图像分割方法.pdf

(19)中华人民共和国国家知识产权局*CN102024260A*(12)发明专利申请(10)申请公布号CN102024260A(43)申请公布日2011.04.20(21)申请号201010614127.4(22)申请日2010.12.30(71)申请人浙江工商大学地址310018浙江省杭州市下沙高教园区学正街18号(72)发明人王勋孔丁科宋超吴海燕(74)专利代理机构杭州求是专利事务所有限公司33200代理人杜军(51)Int.Cl.G06T7/00(2006.01)权利要求书2页说明书5页附图2页(54)发明名称基于局部Gamma拟合的活动轮廓SAR图像分割方法(57)摘要本发明公开了一种基于局部Gamma拟合的活动轮廓SAR图像分割方法。现有的技术SAR图像分割精度不高，本发明首先利用Gamma分布表征原始SAR图像中每个像素一定邻域内的局部结构，以极大似然作为该邻域的区域分离准则，并定义局部拟合能量项；然后在定义域内对所有像素点进行积分，从而得到全局优化的能量泛函；最后采用基于曲线演化理论和水平集方法的几何活动轮廓模型描述和求解能量最小化过程，从而实现SAR图像的有效分割。本发明具有较强的斑点噪声抑制能力，能够实现SAR图像的有效分割，特别是针对深度凹陷等复杂边界仍有较为精确的分割结果。CN10246ACCNN102024260102024274AA权利要求书1/2页1.基于局部Gamma拟合的活动轮廓SAR图像分割方法，其特征在于该方法包括以下步骤：步骤1将定义域为的SAR图像划分为个全剖分的同质子区域，该同质区域集满足，；对定义域上的每个像素点的邻域内执行步骤1.1～步骤1.3，实现局部划分:步骤1.1根据相干斑噪声模型,目标区域均匀的SAR图像满足Gamma分布，利用Gamma分布拟合邻域内的局部结构，定义区域的概率模型，对于任意：其中是SAR图像的视数，是区域的强度均值，是Gamma函数；步骤1.2利用极大似然作为邻域的区域分离准则：，最大化等价于最小化，以定义能量泛函即可表征局部区域分离；步骤1.3结合用以控制邻域采样范围的权函数，得到邻域的局部拟合能量项，从而实现局部划分：；用水平集函数表示局部二分拟合能量项：，其中，，是单位阶跃函数；步骤2根据步骤1所得的点局部能量项，在定义域内对积分得到全局的能量泛函，其中用以纠正水平集函数与符号距离函数的偏差，避免演化过程中的周期性初始化过程，是长度约束项,是梯度算子；步骤3利用Euler-Lagrange变分极小化步骤2的能量泛函，得到曲线演化的控制方程:，2CCNN102024260102024274AA权利要求书2/2页是时间变量，式中是狄拉克函数，是演化曲线的曲率，，其中是局部区域均值，满足:；步骤4利用有限差分方法对步骤3所得的曲线演化控制方程进行迭代求解，最终得到闭合分割曲线就是由所组成，即是原始SAR图像的分割结果：,是时间步长，式中为控制方程右边的差分离散形式。3CCNN102024260102024274AA说明书1/5页基于局部Gamma拟合的活动轮廓SAR图像分割方法技术领域[0001]本发明属于图像处理领域，涉及一种图像分割方法，特别涉及一种SAR图像的活动轮廓图像分割方法。背景技术[0002]合成孔径雷达（SAR）以其全天时、全天候的高分辨率数据获取能力成为军事及民用领域普遍采用的对地观测技术。图像分割是SAR图像自动解译的基础性工作,受到了广泛的关注。由于SAR系统采用相干成像处理,致使SAR图像中存在大量相干斑噪声,这也是SAR图像精确分割的重要障碍。立足于消除斑点噪声的影响,传统SAR分割方法需求通过相干斑预处理,但由于预处理在抑制相干斑的同时,不可避免地损失了潜在边界等细节信息,致使分割效果不佳。[0003]目前,基于曲线演化理论和水平集方法的几何活动轮廓模型（ActiveContourModel,ACM）受到了众多研究者的关注,SAR图像的ACM分割研究也方兴未艾。ACM实现图像分割的基本思想是:将演化曲线（活动轮廓）隐含地表示为更高维曲面的零水平集,水平集在演化方程控制下进行演化,直至图像的目标边界。ACM有许多优点:曲线演化过程始终处于固定格网,数值实现简单；能够自动处理边界的拓扑变化；提供光滑、闭合的高精度分割曲线,方便后续的图像解译工作；特别是区域ACM能够利用图像全局信息定义能量泛函,对斑点噪声具有较好的抑制作用,从而避免了复杂的相干斑预处理过程。[0004]然而,现有SAR图像的ACM分割方法均基于图像区域内目标地物的强度统计,其分割精度取决于统计分布对图像数据的整体拟合程度,而忽略了边界局部结构的有效表征,从而导致分割结果不理想,特别是在处理深度凹陷等复杂边界情况时存在较大局限。发明内容[0005]本发明针对现有技术的