[发明专利]多视SAR图像分割方法及装置

摘要

本发明提供了一种多视SAR图像分割方法及装置，其中，该方法包括：读取待分割多视SAR图像；初始化双权重w；重复执行下述步骤：计算Gamma分布尺度参数β；计算上述图像的或然率p(z|w)；计算上述双权重w的分布函数p(w)；计算品质函数L；根据梯度法更新上述双权重w；将更新后的双权重w代入品质函数L；直到|L^(t+1)‑L^(t)|小于预设的阈值ε时，停止执行上述步骤，根据当前的双权重w确定上述图像中各个像素所属的类别；按照各个像素所属的类别输出分割结果。本发明分割图像的抗噪性能较好，分割结果误分现象较少，分割边界拟合精确。

主权项

1.一种多视SAR图像分割方法，其特征在于，包括：读取待分割多视SAR图像；所述待分割多视SAR图像使用特征场z表示：z＝{zi(xi,yi):i＝1,…,n}其中，i为像素索引，n为总像素数，zi为像素i的强度，(xi，yi)∈D为像素i的格点位置，D为图像域；初始化双权重w：w＝(wil:l＝1,…,k)其中，w_il为类属l包含像素i的类属权重，满足k为类数；重复执行下述步骤：计算Gamma分布尺度参数βl：其中，α为Gamma分布形态参数；计算所述图像的或然率p(z|w)；其中p(zi|wil)为以双权重Gamma混合模型定义的zi的概率密度，如下：其中Ga(zi|βl)为以zi为自变量，βl为尺度参数的Gamma分布概率密度函数，Γ(α)为以α为参数的Gamma函数；计算所述双权重w的分布函数p(w)：其中，A为归一化系数，η为邻域作用系数，Ni为以(xi,yi)为中心的8邻域像素集合，且满足i’∈Ni,i’≠i；计算品质函数L，所述品质函数L为所述或然率p(z|w)与所述分布函数p(w)联合概率分布函数的对数函数；根据梯度法更新所述双权重w：w(t+1)＝w(t)+ξΔw(t)其中，t为迭代次数，ξ为步长，Δw(t)为梯度，如下：将更新后的所述双权重w代入所述品质函数L；直到|L(t+1)‑L(t)|小于预设的阈值ε时，停止执行上述步骤，根据当前的所述双权重w确定所述图像中各个像素所属的类别；按照所述各个像素所属的类别输出分割结果。