一种基于单点模糊初始化水平集演化的遥感图像溢油分割方法技术

本发明专利技术公开了一种基于单点模糊初始化水平集演化的遥感图像溢油分割方法，其步骤为：首先，在原始图像溢油区域内任选一像素作为种子像素，计算图像中各像素相对于种子像素的模糊连接度；然后，基于模糊连接度构建初始水平集函数；最后，基于区域可伸缩匹配(Region Scalable Fitting)水平集算法进行水平集演化，至收敛时，水平集函数大于零的区域对应于遥感图像中的溢油区域。本发明专利技术克服了区域可伸缩匹配水平集算法对于不同初始水平集函数会得到不同分割结果的缺点，能够将遥感溢油图像进行精确的分割。

Remote sensing image oil spill segmentation method based on single point fuzzy initialization level set evolution

The invention discloses a remote sensing image single point fuzzy initialization level set segmentation method based on the evolution of the oil spill, which comprises the following steps: firstly, the oil spill area to choose a pixel as seed pixels in the original image, each pixel in the image is calculated relative to the fuzzy seed pixel connectivity; then, fuzzy connectedness to construct the initial level set finally, based on function; based on region scalable matching (Region Scalable Fitting) level set algorithm for level set evolution, until convergence, oil spill domain function is greater than zero corresponds to remote sensing image level set. The invention overcomes the disadvantage that the region scalable matching level set algorithm can obtain different segmentation results for different initial level set functions, and can remotely segment the oil spill images.

【技术实现步骤摘要】
一种基于单点模糊初始化水平集演化的遥感图像溢油分割方法
本专利技术涉及遥感图像处理
，具体涉及一种基于单点模糊初始化水平集演化的遥感图像溢油分割方法。
技术介绍
海上石油开采时常会有泄露的情况发生。为减小溢油造成的影响，需及时检测溢油区域。基于海面遥感图像的溢油区域检测任务通常是将遥感图像分割为溢油区域和非溢油区域的过程。区域可伸缩匹配(RegionScalableFitting)水平集演化算法(参见文献：LiC,KaoCY,GoreJC,etal.Minimizationofregion-scalablefittingenergyforimagesegmentation.[J].IEEETransactionsonImageProcessing,2008,17(10):1940-1949)是一类广泛应用的图像分割方法。该方法需要设定待分割区域的初始轮廓。在水平集演化算法中，初始轮廓对应于初始水平集中的零水平集函数。实验证明，该方法对初始轮廓具有依赖性：对于不同的初始轮廓，其分割结果不同。将该方法应用到溢油分割时，对于与溢油区域相差较大的初始轮廓，其分割结果不准确。如图2、图3、图4和图5中(b)和(c)中白实线分别是手动设置的不同初始轮廓以及对应的分割结果。由上可知，区域可伸缩匹配水平集算法对初始轮廓的依赖，提升了初始轮廓设定的难度，给溢油检测的实际操作造成了困难，为溢油自动检测造成了困扰。
技术实现思路
为了解决上述问题，本专利技术根据模糊连接度(参见文献：UdupaJK,SahaPK.Fuzzyconnectednessandimagesegmentation[J].ProceedingsoftheIEEE,2003,91(10):1649-1669)提出了一种基于单点模糊初始化水平集演化的遥感图像溢油分割方法。该方法有效解决了
技术介绍
中提到的水平集演化对初始化轮廓的依赖问题，能够实现对遥感溢油图像的精确分割。本专利技术的技术方案是：一种基于单点模糊初始化水平集演化的遥感图像溢油分割方法，包含以下步骤：步骤一：计算图像中各像素相对于种子像素的模糊连接度：定义I表示一张分辨率为X×Y的遥感溢油灰度图像，I是X×Y维矩阵，I(x,y)是矩阵(x,y)位置的元素，表示溢油图像(x,y)位置像素的灰度值。在原始图像溢油区域内任选一个(x0,y0)位置的像素作为种子像素。计算图像中所有位置中任意两位置所对应像素间的模糊邻近度，位于任意(xi,yi)和(xj,yj)位置两像素间模糊邻近度ρi,j如下计算：其中，k1是距离递减因子，为非负常数，用来调整图像中两像素间的距离权重；计算图像中所有位置中任意两位置所对应像素间的模糊相似度，位于任意(xi,yi)与(xj,yj)位置两像素间模糊相似度ai,j如下计算：其中，k2是梯度权重，为非负常数，是平衡像素强度差异的权重；计算图像中所有位置的像素与(x0,y0)位置的种子像素之间的模糊连接度，位于任意(x,y)位置的像素与位于(x0,y0)位置的种子像素之间的模糊连接度如下计算：其中，P是位置(x,y)到(x0,y0)的所有路径，p是P中的一条路径。步骤二：基于模糊连接度构建初始水平集函数：设定阈值θ，基于θ和构建初始水平集函数φ0(x,y)：步骤三：以φ0(x,y)为初始水平集函数，基于区域可伸缩匹配进行水平集演化，具体执行过程为：定义平滑阶跃函数为：H∈(φ)导数为：定义基于上述定义，构建区域可伸缩匹配函数：其中，Kσ为高斯核函数σ是它的尺度参数。基于Kσ和qi(x,y)，定义数据拟合函数：ei(x,y)＝∫∫Kσ(x-u,y-υ)|I(u,υ)-qi(x,y)|2dudυ(8)基于公式(6)和公式(8)，梯度流方程如下：其中，λi是正数，用来平衡各个数据拟合函数的权重；v、η均是非负参数，保证水平集演化的正则性；以公式(4)为t＝0时初始水平集函数，构建水平集演化方程为：式中，α是演化速率。经公式(10)迭代演化至收敛，得到最终的水平集函数，其大于零的区域分割为遥感图像中的溢油区域。本专利技术的有益效果为：图2、图3、图4、图5中，(b)、(c)中白实线是采用区域可伸缩匹配水平集算法设定的不同初始轮廓，(g)、(h)是对应的分割结果，其中黑色区域是分割出的溢油区域。由分割结果可知，初始轮廓不同，分割结果不同，且分割结果不准确。(d)、(e)中溢油区域内的白点是对于同一幅遥感溢油图像采用本专利技术所提单点模糊初始化方法分别设定的不同初始种子像素，图像中的白虚线是由设定的不同初始种子像素所得的初始轮廓(即公式(4)中初始水平集φ0(x,y)的零水平集函数φ0(x,y)＝0曲线)，(i)、(j)是对应的分割结果，其中黑色区域是分割出的溢油区域。由分割结果可知，初始种子像素不同，所得分割结果相同，且分割结果准确。附图说明图1为本专利技术一种基于单点模糊初始化水平集演化的遥感图像溢油分割方法的原理框图。图2、图3、图4和图5是四组实验结果图。其中，(a)、(f)分别是原始遥感图像和真实的分割结果；(b)、(c)中的白实线是采用区域可伸缩匹配初始化水平集算法分别设定的不同初始轮廓，(g)、(h)是对应的分割结果，其中黑色区域是分割出的溢油区域；(d)、(e)中，溢油区域内的白点是采用本专利技术所提方法设定的不同初始种子像素，图像中的白虚线是由设定的不同初始种子像素所得的初始轮廓(即公式(4)中初始水平集φ0(x,y)的零水平集函数φ0(x,y)＝0曲线)，(i)、(j)是对应的分割结果，其中黑色区域是分割出的溢油区域。具体实施方式为了更好地说明本专利技术的技术方案，以下结合附图对本专利技术的实施方式作进一步描述。一种基于单点模糊初始化水平集演化的遥感图像溢油分割方法，包含以下步骤：步骤一：参见图1，计算图像中各像素相对于种子像素的模糊连接度：定义I表示一张分辨率为X×Y的遥感溢油灰度图像，I是X×Y维矩阵，I(x,y)是矩阵(x,y)位置的元素，表示溢油图像(x,y)位置像素的灰度值。在原始图像溢油区域内任选一个(x0,y0)位置的像素作为种子像素。计算图像中所有位置中任意两位置所对应像素间的模糊邻近度，位于任意(xi,yi)和(xj,yj)位置两像素间模糊邻近度ρi,j如下计算：其中，k1是距离递减因子，为非负常数，用来调整图像中两像素间的距离权重；计算图像中所有位置中任意两位置所对应像素间的模糊相似度，位于任意(xi,yi)与(xj,yj)位置两像素间模糊相似度ai,j如下计算：其中，k2是梯度权重，为非负常数，是平衡像素强度差异的权重；计算图像中所有位置的像素与(x0,y0)位置的种子像素之间的模糊连接度，位于任意(x,y)位置的像素与位于(x0,y0)位置的种子像素之间的模糊连接度如下计算：其中，P是位置(x,y)到(x0,y0)的所有路径，p是P中的一条路径。具体的，运用动态规划方法实现公式(13)的计算。步骤二：参见图1，基于模糊连接度构建初始水平集函数：设定阈值θ，基于θ和构建初始水平集函数φ0(x,y)：步骤三：参加图1，以φ0(x,y)为初始水平集函数，基于区域可伸缩匹配进行水平集演化，具体执行过程为：定义平滑阶跃函数为：H∈(φ)导数为：定义基于上述定义，构造区域可伸缩本文档来自技高网...

【技术保护点】
一种基于单点模糊初始化水平集演化的遥感图像溢油分割方法，其特征在于：包含以下步骤：步骤一：设定溢油区域中的种子像素，计算图像中各像素相对于种子像素的模糊连接度：定义I表示一张分辨率为X×Y的遥感溢油灰度图像，I是X×Y维矩阵，I(x,y)是矩阵中(x,y)位置的元素，表示溢油图像(x,y)位置像素的灰度值。在原始图像溢油区域内任选一个(x0,y0)位置的像素作为种子像素。计算图像中所有位置中任意两位置所对应像素间的模糊邻近度，位于任意(xi,yi)和(xj,yj)位置两像素间模糊邻近度ρi,j如下计算：

【技术特征摘要】
1.一种基于单点模糊初始化水平集演化的遥感图像溢油分割方法，其特征在于：包含以下步骤：步骤一：设定溢油区域中的种子像素，计算图像中各像素相对于种子像素的模糊连接度：定义I表示一张分辨率为X×Y的遥感溢油灰度图像，I是X×Y维矩阵，I(x,y)是矩阵中(x,y)位置的元素，表示溢油图像(x,y)位置像素的灰度值。在原始图像溢油区域内任选一个(x0,y0)位置的像素作为种子像素。计算图像中所有位置中任意两位置所对应像素间的模糊邻近度，位于任意(xi,yi)和(xj,yj)位置两像素间模糊邻近度ρi,j如下计算：其中，k1是距离递减因子，为非负常数，用来调整图像中两像素间的距离权重；计算图像中所有位置中任意两位置所对应像素间的模糊相似度，位于任意(xi,yi)与(xj,yj)位置两像素间模糊相似度ai,j如下计算：其中，k2是梯度权重，为非负常数，是平衡像素强度差异的权重；计算图像中所有位置的像素与(x0,y0)位置的种子像素之间的模糊连接度，位于任意(x,y)位置的像素与位于(x0,y0)位置的种子像素之间的模糊连接度如下计算：其中，P是位置(x,y)到(x0,y0)的所有路径，p是P中的一条路径。步骤二：基于模糊连接度构建初始水平集函数：设定阈值θ，基于θ和构建初始水平集函数φ0(x,y)：步骤三：以φ0(x,y)为初始水平集函数，基于区域可伸缩匹配进行水平集演化，具体执行过程为：定义平滑阶跃函...

【专利技术属性】
技术研发人员：任鹏，许敏，陈芳，余兴瑞，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东,37