[发明专利]SAR图像海陆分割方法、装置、计算机设备和可读存储介质

权利要求书

1.一种SAR图像海陆分割方法，其特征在于，所述方法包括：

(1)获取待分割SAR图像；

(2)获取待分割SAR图像的能量函数，在所述能量函数中引入所述待分割SAR图像的水平集函数，获取所述能量函数最小时的第一分割模板，包括：

获取所述待分割SAR图像的基于C-V水平集模型的能量函数；

在所述能量函数中引入水平集，求解关于水平集的偏微分方程，得到迭代函数；

对迭代函数进行优化，利用变分法得到水平集演化方程；

利用有限差分进行离散化得到所述第一分割模板；

(3)根据所述待分割SAR图像，对所述第一分割模板进行迭代更新预设次数，得到第二分割模板，包括：

获取待分割SAR图像的初始能量值和灰度值；

根据所述初始能量值、所述灰度值和所述迭代函数，对所述第一分割模板迭代预设次数，得到所述第二分割模板，预设次数为十次；

(4)根据预设算子检测边缘，提取所述第二分割模板中的疑似误分割孔洞，包括：

根据预设算子检测边缘和所述第二分割模板，得到二值化模板；

提取所述二值化模板中的连通分量，得到所述疑似误分割孔洞；

(5)获取所述疑似误分割孔洞的纹理特征，根据所述纹理特征填充所述疑似误分割孔洞中的误判孔洞，得到第三分割模板，包括：

根据所述疑似误分割孔洞的所述纹理特征，查询所述疑似误分割孔洞中的误判孔洞和非误判孔洞；

其中，利用一阶偏导的有限差分，计算二值孔洞I(x,y)的梯度幅度，E_x(i,j)和E_y(i,j)分别为沿x和y方向的偏导数：

E_x(i,j)＝[I(i,j)-I(i+1,j)+I(i,j+1)-I(i+1,j+1)]/2

E_y(i,j)＝[I(i,j)-I(i,j+1)+I(i+1,j)-I(i+1,j+1)]/2

孔洞中每个像素点(i,j)的梯度幅度值M(i,j)为：

计算每个孔洞的平均梯度幅度，作为孔洞的纹理特征，生成特征矩阵作为FCM的聚类对象，判断其属于毛刺较多型孔洞或平滑型孔洞，毛刺较多型孔洞为误判孔洞，平滑型孔洞为非误判孔洞；

将所述误判孔洞的灰度值设置为第一预设值，将所述非误判孔洞的灰度值设置为第二预设值，得到第三分割模板；

(6)根据所述第三分割模板和所述第二分割模板获取第四分割模板，将所述待分割SAR图像输入所述第四分割模板，得到分割图像。

2.根据权利要求1所述的方法，其特征在于，所述获取待分割SAR图像，包括：

获取初始待分割SAR图像；

对所述初始待分割SAR图像进行图像增强和高斯模糊平滑噪声，得到所述待分割SAR图像。

3.根据权利要求1所述的方法，其特征在于，所述根据所述第三分割模板和所述第二分割模板获取第四分割模板，包括：

将所述第三分割模板和所述第二分割模板进行异或运算，得到初始第四分割模板；

对所述初始第四分割模板进行形态学过滤，得到所述第四分割模板。