[发明专利]一种对海面SAR图像的滤波方法

摘要

本发明公开了一种对海面SAR图像的滤波方法，结合最大类间方差法和形态学滤波方法对海洋溢油的图像进行滤波；具体讲，通过最大类间方差法得到图像的二值图，再对二值图分别进行两次膨胀和两次腐蚀；利用腐蚀后二值图灰度值减去膨胀后图像的灰度值得到包含边界的局部区域二值图；最后，局部区域和原图像灰度进行内积，就可以把原图像进行简化，只剩下包含了溢油与疑似溢油边界的滤波区域，最后再通过滤波模板对滤波区域进行滤波。

主权项

1.一种对海面SAR图像的滤波方法，其特征在于，包括以下步骤：(1)、图像预处理提取海面SAR图像，并处理成灰度图，得到原始灰度图像；(2)、图像分割通过最大类间方差法对原始灰度图像进行二值分割，得到二值图；(3)、提取滤波区域利用形态学滤波对二值图分别进行两次膨胀和两次腐蚀，再利用两次膨胀后的二值图的灰度值减去对应坐标下两次腐蚀后的二值图的灰度值，得到由滤波区域组成的二值图，最后，将滤波区域组成的二值图与原始灰度图像的对应像素点的灰度值进行内积，得到滤波区域；(4)、对滤波区域进行滤波(4.1)、利用3x3模板对滤波区域进行初次滤波；设3x3模板的中心像素点m为目标像素点，3x3模板表示为如果点m为滤波区域的边界像素点时，则点m的灰度值保持不变；如果点m不为滤波区域的边界像素点时，3x3模板中各像素点的灰度值满足：满足条件1：则m₁＝(a1+a2)/2；或满足条件2：则m₂＝(b1+b2)/2；或满足条件3：则m₃＝(c1+c2)/2；或满足条件4：则m₄＝(d1+d2)/2；或满足条件5：则m₅＝(a1+a2)/2；或满足条件6：则m₆＝(b1+b2)/2；或满足条件7：则m₇＝(c1+c2)/2；或满足条件8：则m₈＝(d1+d2)/2；其中，k和n为常数，且满足k＜n；最终，点m的灰度值为同时满足上述条件所有m_i(i＝1,2,…,8)之和的平均值，即其中，h₁为3x3模板中满足各像素点灰度值条件的个数；(4.2)、利用5x5的模板对初次滤波后的滤波区域进行第二次滤波；设5x5模板的中心像素点m为目标像素点，5x5模板表示为：如果点m与滤波区域的边界像素点的距离小于等于3个像素点，则点m的灰度值保持不变；如果点m与滤波区域的边界距离大于3个像素点时，则5x5模板中各像素点的灰度值满足：满足条件1：则m₁＝a2＝(a1+a22)/2；或满足条件2：则m₂＝b2＝(b1+b22)/2；或满足条件3：则m₃＝c2＝(c1+c22)/2；或满足条件4：则m₄＝d2＝(d1+d22)/2；或满足条件5：则m₅＝a1＝(a11+a2)/2；或满足条件6：则m₆＝b1＝(b11+b2)/2；或满足条件7：则m₇＝c1＝(c11+c2)/2；或满足条件8：则m₈＝d1＝(d11+d2)/2；或满足条件9：则m₉＝a2＝(a1+a22)/2；或满足条件10：则m₁₀＝b2＝(b1+b22)/2；或满足条件11：则m₁₁＝c2＝(c1+c22)/2；或满足条件12：则m₁₂＝d2＝(d1+d22)/2；或满足条件13：则m₁₃＝a1＝(a11+a2)/2；或满足条件14：则m₁₄＝b1＝(b11+b2)/2；或满足条件15：则m₁₅＝c1＝(c11+c2)/2；或满足条件16：则m₁₆＝d1＝(d11+d2)/2；其中，k和n为常数，且满足k＜n；最终，点m的灰度值为同时满足上述条件所有m_j(j＝1,2,…,16)之和的平均值，即其中，h₂为5x5模板中满足各像素点灰度值条件的个数。