[发明专利]一种海面溢油图像分割方法

说明书

技术领域

本发明涉及海洋污染处理领域，特别涉及一种海面溢油图像分割方法。

背景技术

随着海洋石油工业和海上石油运输业的迅速发展，海上石油勘探、开发、海底管线铺设规模的不断扩大，石油设施已成为海洋环境污染的重大危险源。海上溢油给海洋环境和海洋生态带来了严重的污染。保护海洋环境，防止海洋表面污染已成为目前迫切需要解决的问题。为防止海面溢油污染，多种卫星被部署用来执行跟踪监测海上溢油的日常任务。卫星传回的海面实时图像用于检测海面溢油情况。为有效的检测海面图像中的溢油区域，各式各样的溢油分割方法相继涌现。目前海面溢油分割方法大多是传统的图像分割方法，如阈值分割，区域分割，边缘检测等。尤其是阈值分割作为一种传统的图像分割方法，因其实现简单、计算量小、性能稳定而成为海面溢油图像分割常用的方法。

阈值分割利用图像中要提取的目标区域与其背景区域的差异，把图像看作具有不同灰度级的两类区域(目标区域和背景区域)的组合，选取一个比较合理的阈值，以确定图像中每个像素点所属的类别区域。为了得到海面溢油图像中溢油区域与背景区域的阈值，首先得到图像的灰度直方图。对于含有溢油现象的海面溢油图像，其灰度直方图都表现为存在两个峰值，一般选取两峰之间的极小值作为阈值来区分油膜和背景海水区域。基于阈值分割的海面溢油图像分割虽然能分割出溢油区域，但它只考虑了图像单个像素的灰度值，没有考虑像素之间的相关性，包含的信息少，分割准确性太差。假如图像包含溢油区域的大小不合适，灰度直方图只能出现一个峰值，阈值法不再适用，并且噪声的存在也使得检测出的溢油区域不尽如意。

为克服传统图像分割方法应用于海面溢油图像分割中的问题，Graph Cuts被用于海面溢油图像分割中。将海面溢油图像中的每个像素看作一个布尔变量，按照一定的规则构建二次布尔能量函数，然后根据能量函数构建相应的s/t图，二次布尔能量函数中每一布尔变量对应s/t图的一个节点，再利用最小割/最大流算法切割s/t图得到s/t图的最小割，此时s/t图各节点对应变量的取值使得二次布尔能量最小化。将二次布尔能量最小化时取值为1的布尔变量相对应的像素标为溢油区域，取值为0的布尔能量对应的像素标为非溢油区域，便得到了最优的溢油区域。

以上基于Graph Cuts的溢油分割算法能有效的分割出大片完整的油膜区域，对细长条的油膜带分割效果不佳。为有效分割细长条的溢油区域，高阶能量函数被用于海面溢油图像分割中，在一定程度上可以较为精确的分割长条状溢油带。但基于高阶能量函数的溢油分割方法会导致边缘扩散现象，分割出的溢油区域的边缘有向非溢油区域扩散的趋势，产生带有毛刺的溢油分割结果，影响溢油分割精度。

发明内容

本发明针对现有技术中的不足，本发明针对以上缺点，将保边滤波与高阶能量函数的图像分割算法相结合，发明了一种新的海面溢油图像分割方法，有效分割细长条的溢油区域，可以一定程度上抑制分割过程中的边缘扩散现象，从而提高了分割准确性。

本发明的技术方案是：

具体为构建高阶布尔能量函数，其包含一阶判定项与高阶耦合项。一阶判定项由单像素能量函数项与保边滤波器函数卷积求得，具有较好的保护溢油区域边缘的效果；高阶耦合项能够计算像素子集中的同类成对近邻像素的耦合性，有益于精确地分割出长条状溢油带。构造出的布尔能量函数由最大流/最小割方法求取最小值，能量最小时各布尔变量的取值对应海面溢油图像中溢油区域最优的分割方式。

一种海面溢油图像分割方法，包括以下步骤：

步骤1：首先设定初始一阶判定项

其中，p为像素的标号，I_p为像素p的灰度值，x_p是与像素p对应的布尔变量，η为阈值；