基于sar图像的海洋溢油信息识别与提取方法
【专利摘要】本发明公开了一种基于SAR图像的海洋溢油信息识别与提取方法，通过预处理消除星载SAR图像的噪声和畸变，然后进行图像的普适性处理将非海面区域遮掩，接着使用二维最大类间方差法分割图像获取感兴趣的区域，在此基础上利用形态学滤波矫正误分割现象，进一步从图像分割区域提取选定的特征参量，最后结合先验知识利用经过同样处理的样本图像集训练的BP神经网络对图像分割的结果进行分类，区分海面溢油与疑似溢油，准确提取海洋溢油信息。
【专利说明】
基于SAR图像的海洋溢油信息识别与提取方法
技术领域
[0001] 本发明属于图像处理技术领域，更为具体地讲，设及一种基于SAR图像的海洋溢油 信息识别与提取方法。
【背景技术】
[0002] 随着人类海洋活动的逐渐增加，石油开采、船舶事故等海上溢油事件频繁发生，给 海洋生态环境造成了巨大的威胁，同时也给社会经济造成了重大的损失。海上溢油发生W 后，对溢油进行检测的实时性和准确性对保护海洋生态环境、减少社会经济损失具有重要 的意义。微波遥感技术是利用波长为1毫米到1米的电磁波远距离检测目标地物信息的技 术，通常采用合成孔径雷达（Synthetic Aperture Radar, SAR)运行平台。与光学遥感技术 相比，微波遥感技术具有明显优势:不依赖光照，能够穿透云、雨、雪、雾，因此具有全天时、 全天候工作的能力。运用微波遥感技术监测海面溢油开始于20世纪70年代，不同型号遥感 卫星的升空提供了多种监测海面溢油的途径，许多专家学者应用SAR平台监测海面溢油，取 得了良好效果。随着SAR技术与各种检测算法日益完善，微波遥感技术已经成为海面溢油监 测中不可或缺的方式，同时为海洋环保部口做出最佳的灾害方案提供了最有力的支持。
[0003] 由于星载SAR系统投入业务化运行，与海上溢油相关的SAR图像数据量急剧增加， 日常化、业务化地监测海上溢油已经远远超出传统的专家解译过程的工作量。在处理海量 的与海洋上溢油相关的SAR图像数据时，为了能够及时准确地进行监测，需要建立一套完善 的利用SAR图像数据自动监测海上溢油的系统。
[0004] 海面油膜阻尼了溢油海面的毛细波、短重力波，降低了其表面粗糖度，相较于清洁 海面，溢油海面的后向散射强度较低，在SAR图像上表现为阴影区域。同时，存在其他的一些 海洋现象在SAR图像上表现出相似的情形。例如海面低风速区域、背风区域、浅海区域、海洋 生物活动形成的生物油膜、海洋内波区域等与溢油区域具有相似的微波散射特性，在SAR图 像上也表现为阴影区域，对海上溢油信息的识别与提取造成了不可忽视的影响，产生极高 的虚警率。因此，快速而准确地区分溢油和疑似溢油区域已经成为海上溢油识别亟待解决 的热口问题。
【发明内容】
[0005] 本发明的目的在于克服现有技术的不足，提供一种基于SAR图像的海洋溢油信息 识别与提取方法，利用BP神经网络来实现高精度、快速获取海上溢油信息。
[0006] 为实现上述发明目的，本发明为一种基于SAR图像的海洋溢油信息识别与提取方 法，其特征在于，包括W下步骤：
[0007] (1)、获取原始SAR图像
[000引通过星载SAR传感器，在待监测区域获取原始SAR图像；
[0009] (2)、图像预处理
[0010] 使用ENVI软件对原始SAR图像依次进行福射定标、地理编码及滤波，得到初始SAR 图像；
[0011] (3)、对初始SAR图像进行普适性处理
[0012] (3.1)、使用掩膜遮蔽初始SAR图像中包含大片陆地的区域，得到海洋区域SAR图 像；
[OOK] (3.2)、对海洋区域SAR图像中舰船、岛蜗形成的高亮区域进行处理：引入正态分布 的30原理，设定上限灰度值：
[0014] limit=]i+3〇
[0015] 其中，ii为初始SAR图像灰度均值，O为初始SAR图像灰度标准差；
[0016] 捜索海洋区域SAR图像中灰度值超过该上限灰度值的像素点，再使用灰度值未超 过该上限灰度值的邻近像素点进行替换，得到经过普适性处理的标准SAR图像；
[0017] (4)、对标准SAR图像进行分割
[0018] (4.1)、将标准SAR图像中与后向散射系数相关的像素点的灰度值映射到0到255的 整数范围内，得到像素点的灰度阔值S，同时引入像素点的n*n邻域灰度均值阔值t，组成二 维阔值(s，t);
[0019] (4.2)、根据二维阔值(s，t)，建立二维灰度直方图；利用二维灰度直方图将灰度值 小于S且邻域灰度均值小于t的区域设置为目标，将灰度值大于S且邻域灰度均值大于t的区 设置为背景，其它区域设置为噪声；
[0020] (4.3)、利用二维灰度直方图计算类间方差d(s，t):
[0021] 设二维灰度直方图中的像素点总数为N，设灰度值为f同时n*n邻域灰度均值为/的 像素点个数为％巧^么运一类像素点出现的概率为：
[0022]
[0023] 那么当阔值为(s，t)时，目标区域和背景区域发生的概率分别为：
[0024]
[0025] 其中，下角标a代表目标区域，b代表背景区域；
[00%] 那么两区域对应的灰度均值分别为：
[0027]
[002引其中，A;，/V伯，/?分别表示目标区域和背景区域的灰度均值和邻域灰度 均值的平均值；
[0029] 二维灰度直方图的灰度均值和邻域灰度均值的平均值为：
[0030]
[0031] 由此可得，类间方差d(s，t)为：
[0032] d(S , t) = W a[ (yai-yi)2+(Jiai-yi)2]+ W b[(化广叫）2+(化广化）2];
[003引（4.4)、利用灰度阔值S和邻域灰度均值阔值t分别遍历巧I帖4,获得使d(S，T)在所 有的d(S，t)中取得最大值的最优阔值(S，T);
[0034] (4.5)、根据最优阔值(S，T)对标准SAR图像中的像素点按照步骤(4.2)所述方法进 行分类，分割出标准SAR图像中的目标、背景和噪声，得到目标SAR图像；
[0035] (5)、提取特征参量
[0036] 先利用形态学滤波器对目标SAR图像进行多次形态学滤波，再提取几何、灰度、纹 理=方面的9类特征参量值；
[0037] (6)、利用BP神经网络分类来区分溢油与疑似溢油
[0038] (6.1)、将经过形态学滤波后的目标SAR图像作为训练样本图像，使用训练样本图 样中提取的9个特征参量构成特征向量组作为未经过训练的BP神经网络的输入，根据先验 知识确定滤波后的目标SAR图像的阴影区域是否属于溢油，如果是溢油则设置输出"1 0"， 如果是疑似溢油则设置输出"0 r，并作为该BP神经网络的输出，通过调整内部的权值和阔 值，训练BP神经网络；
[0039] (6.2)、将待监测的原始SAR图像经过上述步骤（1巧Ij步骤(5)处理后提取的特征参 量构成特征向量组，输入经过训练后的BP神经网络，根据BP神经网络的输出结果确定出目 标区域属于溢油还是疑似溢油。
[0040] 本发明的发明目的是运样实现的：
[0041] 本发明基于SAR图像的海洋溢油信息识别与提取方法，通过预处理消除星载SAR图 像的噪声和崎变，然后进行图像的普适性处理将非海面区域遮掩，接着使用二维最大类间 方差法分割图像获取感兴趣的区域，在此基础上利用形态学滤波矫正误分割现象，进一步 从图像分割区域提取选定的特征参量，最后结合先验知识利用经过同样处理的样本图像集 训练的BP神经网络对图像分割的结果进行分类，区分海面溢油与疑似溢油，准确提取海洋 溢油信息。
[0042] 同时，本发明基于SAR图像的海洋溢油信息识别与提取方法还具有W下有益效果：
[0043] (1)、通过对SAR图像使用针对海洋场景区域的普适性处理技术进行处理，能够使 其适应提出的基于SAR图像的海洋溢油信息识别与提取方法，扩展了该方法的适用范围，一 定程度上避免了两种典型的误分割现象，提高了图像分割算法的精度。
[0044] (2)对通过2D-〇tsu算法分割获取的目标区域SAR图像进行多次形态学滤波处理， 在预处理的基础上进一步抑制了SAR图像斑点噪声对溢油信息提取与分类结果的影响，保 留了目标区域SAR图像的几何和纹理特征。
[0045] (3)结合大量的样本数据，多次训练调整内部权值和阔值，获取了稳定的BP人工神 经网络，对溢油和疑似溢油进行了有效地分类，经过验证，分类准确度达到80% W上。
【附图说明】
[0046] 图1是本发明基于SAR图像的海洋溢油信息识别与提取方法流程图；
[0047] 图2是原始SAR图像；
[004引图3是海洋区域SAR图像；
[0049] 图4是经过普适性处理的标准SAR图像；
[0050] 图5是二维灰度直方图的底座划分示意图；
[0051] 图6图像分割后的目标SAR图像；
[0052] 图7是多次形态学滤波后的效果图；
[0053] 图8是特征参量提取结果图；
[0054]图9是特征参量输入到BP神经网络时的结构示意图；
[0化5]图10是BP神经网络分类结果的示意图。
【具体实施方式】
[0056] 下面结合附图对本发明的【具体实施方式】进行描述，W便本领域的技术人员更好地 理解本发明。需要特别提醒注意的是，在W下的描述中，当已知功能和设计的详细描述也许 会淡化本发明的主要内容时，运些描述在运里将被忽略。
[0057] 实施例
[0058] 图1是本发明基于SAR图像的海洋溢油信息识别与提取方法流程图。
[0059] 在本实施例中，如图1所示，本发明提出一种基于SAR图像的海洋溢油信息识别与 提取方法，包括W下步骤：
[0060] Sl、获取原始SAR图像
[0061] 通过星载SAR传感器，在待监测区域获取原始SAR图像;在本实施例中，通过搭载C 波段SAR传感器的高分辨率商用卫星RADARSAT-2获取的高分辨率的原始SAR图像，如图2所 /J、- O
[0062] S2、图像预处理
[0063] 使用肥ST软件对原始SAR图像依次进行福射定标、地理编码及增强型Lee滤波，得 到初始SAR图像。
[0064] S3、对初始SAR图像进行普适性处理
[0065] S3.1、海面溢油现象有时发生在海岸线附近，运种情况下，与海面溢油相关的SAR 图像通常包含一部分陆地区域，图像分割算法一般会将陆地视为背景区域，导致无法有效 获取海面阴影区。因此，在分割包含大片陆地区域的初始SAR图像前，使用掩膜遮蔽，得到海 洋区域SAR图像，具体的掩膜可W通过感兴趣区域获取，如图3所示；
[0066] S3.2、截取感兴趣的海洋区域时，不可避免的存在包含由于大的舰船，小的岛蜗等 后向散射系数明显高于海水的典型地物的像素点，在海洋区域SAR图像中表现为高亮区域。 此类像素点达到一定规模时，针对海洋溢油的图像分割算法通常将高亮区域划分为背景， 其他区域划分为目标。
[0067] 解决该问题的方法是引入正态分布的30原理，设定上限灰度值：
[0068] limit=]i+3〇
[0069] 其中，y为初始SAR图像灰度均值，O为初始SAR图像灰度标准差；
[0070] 捜索海洋区域SAR图像中灰度值超过该上限灰度值的像素点，再使用灰度值未超 过该上限灰度值的邻近像素点进行替换，得到经过普适性处理的标准SAR图像，如图4所示；
[0071] S4、对标准SAR图像进行分割
[0072] S4.1、为了适应2D-〇tsu图像分割算法的适用条件，首先将标准SAR图像中与后向 散射系数相关的像素点的灰度值映射到0到255的整数范围内，得到像素点的灰度阔值S，其 映射关系为：
[0073]
[0074] 其中，Aij为标准SAR图像中像素点Q，j)的灰度值，Bij为该像素点映射到0~255时 的灰度值，Amin和Amax分别为标准SAR图像的所有像素点中的最小、最大灰度值，round ()为取 整符号；
[0075] 同时，引入像素点的n*n邻域灰度均值阔值t，进而组成二维阔值（S，t);在本实施 例中，n邻域的取值一般为3、5、7;
[0076] S4.2、根据二维阔值(s，t)，建立二维灰度直方图，再利用二维灰度直方图将灰度 值小于S且邻域灰度均值小于t的区域设置是目标，将灰度值大于S且邻域灰度均值大于t的 区设置为是背景，其它区域设置为噪声；
[0077] 其中，二维灰度直方图的底座是由灰度值与邻域灰度均值形成的二维平面，在二 维平面上的数据绝大部分集中在其对角线附近，因此可根据阔值将底座划分为4个区域，如 图所示5;
[0078] S4.3、利用二维灰度直方图计算类间方差d(s，t):
[0079] 设二维灰度直方图中的像素点总数为N，设灰度值为?"同时n袖邻域灰度均值为I的 像素点个数为么运一类像素点出现的概率为：
[0080] 巧=巧
[0081] 那么当阔值为(s，t)时，目标区域和背景区域发生的概率分别为：
[0082]
[0083]
[0084]
[0085]
[0086] 其中，，/V伯，分别表示目标区域和背景区域的灰度均值和邻域灰度 均值的平均值；
[0087] 二维灰度直方图的灰度均值和邻域灰度均值的平均值为：
[008引
[0089] 由此可得，类间方差d(s，t)为：
[0090] d(S , t) = W a[ (yai-yi)2+(Jiai-yi)2]+ W b[(化广叫）2+(化广化）2];
[0091] S4.4、利用灰度阔值S和邻域灰度均值t分别遍历巧Ij254,获得使d(S，T)在所有的d (S，t)中取得最大值的最优阔值(S，T);
[0092] S4.5、根据最优阔值(S，T)对标准SAR图像中的像素点按照步骤S4.2所述方法进行 分类，分割出标准SAR图像中的目标、背景和噪声，得到目标SAR图像，如图6所示；
[0093] S5、形态学滤波
[0094] 基于形态学理论，先利用抑制相干斑点噪声的滤波器对目标SAR图像进行多次的 形态学滤波，在目标SAR图像的滤波处理中，1次开运算加一次闭运算作为1次滤波过程;为 了取得更好的相干斑点噪声抑制效果，需要逐步加大选择的结构体尺寸，进行多次滤波，滤 波后的目标SAR图像如图7所示。
[00巧]S6、提取特征参量
[0096] 根据形态学滤波后的目标SAR图像，可W准确的分割出目标和背景两部分，再从中 提取几何、灰度、纹理=方面的9类特征参量值，其提取结果如图8所示；
[0097] 在滤波后的目标SAR图像中，疑似溢油区域与溢油区域呈现近似的灰度值，都表现 为阴影区域，因此需要提取该SAR图像中几何、灰度、纹理=方面的多种特征参量进行进一 步的判断，但是，过多的特征参量不仅增加运算成本，而且与海上溢油特征关系较小的参量 对分类会产生消极影响，降低判定准确性。因此，需要对各种的SAR图像特征参量进行取舍， 本实施例中，提取出9类特征参量值包括：
[0098] a、几何特征
[0099] 面积(Area):由于低风速区海面的后向散射系数非常小而面积非常大，因此，面积 在一定程度上反映了目标区域的大小特征。
[0100] 复杂度(complexity, COM):不同海面目标呈现各不相同的形状特征。油田溢油形 成粗条带或较大面积的溢油区域;船舶漏油在SAR图像中表现为白亮点后的由细到粗的条 带或较粗的均匀条带；自然界的生物油膜一般呈不规则的形状，周长较大而面积较小。因 此，复杂度在一定程度上反映了目标区域的形状特征。
[0101] b、灰度特征
[0102] 灰度对比度(RBIO):目标区、背景区后向散射系数均值的比值。
[0103] RBSD0:背景区后向散射系数均值与其方差的比值。
[0104] 边缘梯度(Gradient Of edge,G0E):目标区、背景区的边界处梯度的均值。
[01化]C、纹理特征
[0106]灰度共生矩阵(化CM)中的元素片7. 表示图像中灰度值为J的像素点，在指定 的像素空间d及方向0上出现灰度值为)的像素点的概率，取d = 1，0 = 0°时，简记为P(Jj)。 [0107] 角二阶矩(Angular Second Moment,ASM):图像灰度分布均匀性的度量。当图像纹 理粗糖时，ASM取值较大；图像纹理平滑时，ASM取值较小;GLCM中各元素值相等时，ASM取得 最小值，其计算公式为：
[010 引
[0109] 纹理对比度(Con化ast，C0N):图像明暗对比的度量值。图像的局部灰度值的差别 越大，纹理对比度越大，视觉效果也越清晰，其计算公式为：
[0110]
[0111] Ic有的信息量的度量，表征了图像纹理的复杂度。图 像纹理自 崎值越小，其计算公式为：
[0112]
[0113] Terence Moment,畑M):图像整体对比度的度量。图 像整体 > 巧比度较高，RDM较小其计算公式为：
[0114]
[0115] S7、利用BP神经网络分类来区分溢油与疑似溢油
[0116] S7.1、将滤波后的目标SAR图像作为训练样本图像，将提取的9个特征参量构成特 征向量组，再将特征向量组作为未经过训练的BP神经网络的输入，如图9所示，根据先验知 识确定滤波后的目标SAR图像的阴影区域是否属于溢油，如果是溢油输出为"1 0"，如果是 疑似溢油输出"0 1"，并作为该BP神经网络的输出，调整内部的权值和阔值，训练BP神经网 络；
[0117] S7.2、将待监测的原始SAR图像经过上述处理后构成特征向量组，再输入经过训练 后的BP神经网络，根据BP神经网络的输出结果确定出目标区域属于溢油还是疑似溢油，其 最终的分类结果如图10所示。
[0118] 尽管上面对本发明说明性的【具体实施方式】进行了描述，W便于本技术领域的技术 人员理解本发明，但应该清楚，本发明不限于【具体实施方式】的范围，对本技术领域的普通技 术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，运些 变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
【主权项】
1. 一种基于SAR图像的海洋溢油信息识别与提取方法，其特征在于，包括以下步骤： (1) 、获取原始SAR图像 通过星载SAR传感器，在待监测区域获取原始SAR图像； (2) 、图像预处理 使用ENVI软件对原始SAR图像依次进行辐射定标、地理编码及滤波，得到初始SAR图像； (3 )、对初始SAR图像进行普适性处理 (3.1) 、使用掩膜遮蔽初始SAR图像中包含大片陆地的区域，得到海洋区域SAR图像； (3.2) 、对海洋区域SAR图像中舰船、岛屿形成的高亮区域进行处理：引入正态分布的3〇 原理，设定上限灰度值： 1?πι?? = μ+3σ 其中，μ为初始SAR图像灰度均值，σ为初始SAR图像灰度标准差； 搜索海洋区域SAR图像中灰度值超过该上限灰度值的像素点，再使用灰度值未超过该 上限灰度值的邻近像素点进行替换，得到经过普适性处理的标准SAR图像； (4)、对标准SAR图像进行分割 (4.1) 、将标准SAR图像中与后向散射系数相关的像素点的灰度值映射到O到255的整数 范围内，得到像素点的灰度阈值s，同时入像素点的η*η邻域灰度均值阈值t，组成二维阈值 (s，t); (4.2) 、根据二维阈值(s，t)，建立二维灰度直方图；利用二维灰度直方图将灰度值小于 S且邻域灰度均值小于t的区域设置为目标，将灰度值大于s且邻域灰度均值大于t的区设置 为背景，其它区域设置为噪声； (4.3) 、利用二维灰度直方图计算类间方差d(s，t): 设二维灰度直方图中的像素点总数为N，设灰度值为『同时n*n邻域灰度均值为)的像素 点个数为;% >那么这一类像素点出现的概率为： Pr1=nTjl N 那么当阈倌为(S，t )时，目标K域和背景区域发生的概率分别为：其中，下角标a代表目标区域，b代表背景区域； 那么两区域对应的灰度均值分别为：其中，夂〃 分别表示目标区域和背景区域的灰度均值和邻域灰度均值 的平均值； 二维灰度直方图的灰度均值和邻域灰度均值的平均值为： d(4.4) 、利用灰度阈值s和邻域灰度均值阈值t分别遍历1到254，获得使d(S，T)在所有的 d (s，t)中取得最大值的最优阈值(S，T); (4.5) 、根据最优阈值(S，T)对标准SAR图像中的像素点按照步骤(4.2)所述方法进行分 类，分割出标准SAR图像中的目标、背景和噪声，得到目标SAR图像； (5) 、提取特征参量 先利用形态学滤波器对目标SAR图像进行多次形态学滤波，再提取几何、灰度、纹理三 方面的9类特征参量值； (6) 、利用BP神经网络分类来区分溢油与疑似溢油 (6.1) 、将经过形态学滤波后的目标SAR图像作为训练样本图像，使用训练样本图样中 提取的9个特征参量构成特征向量组作为未经过训练的BP神经网络的输入，根据先验知识 确定滤波后的目标SAR图像的阴影区域是否属于溢油，如果是溢油则设置输出"1 0"，如果 是疑似溢油则设置输出"1 0"，并作为该BP神经网络的输出，通过调整内部的权值和阈值， 训练BP神经网络； (6.2) 、将待监测的原始SAR图像经过上述步骤（1)到步骤(5)处理后提取的特征参量构 成特征向量组，输入经过训练后的BP神经网络，根据BP神经网络的输出结果确定出目标区 域属于溢油还是疑似溢油。2.根据权利要求1所述的基于SAR图像的海洋溢油信息识别与提取方法，其特征在于， 所述步骤(4.1中，将标准SAR图像中与后向散射系数相关的像素点的灰度值映射到0到255 @敕撒茁困*的古社乇.其中，Aij为标准SAR图像中像素点（i，j)的灰度值，Bij为该像素点映射到0~255时的灰 度值，Ami4PAmax分别为标准SAR图像的所有像素点中的最小、最大灰度值，roundO为取整符 号。
【文档编号】G06K9/36GK106022288SQ201610370277
【公开日】2016年10月12日
【申请日】2016年5月30日
【发明人】陈彦, 赵鹏, 赵讯, 童玲
【申请人】电子科技大学