海面溢油区域图像探测方法与流程

本发明属于图像处理领域，具体地说是通过选择适当的波段并进行相应的处理来将遥感数据转化成灰度图像，然后对灰度图像进行处理进而得到溢油区的轮廓。

背景技术：

溢油污染是严重的海洋环境灾害之一，溢油灾害一旦发生就会直接危害海洋生态环境，给沿岸的社会经济和人类的身体健康带来直接的影响。海洋石油污染有多种途径,既有天然来源如海底油气藏烃渗漏和沉积岩石的侵蚀，也有沿岸工业污水和生活废水的排放、海洋倾废,更有海上石油运输和生产所造成的石油泄漏。海洋溢油发生后，能否准确及时的监测溢油对于海洋环境保护具有重要意义。

过去检测油膜主要依靠直接测量，具有检测覆盖面积小、判断主观等缺点。航空遥感技术是目前世界上发达国家进行海洋检测、监视的重要手段，具有速度快、机动灵活、覆盖面积较大、视距范围较宽、成本低、光谱和空间分辨率高等特点。利用卫星数据不仅可以大面积监测海上溢油的面积、种类、厚度，及时引导海监船只和飞机进行执法监测，作为执法索赔依据，而且可以利用卫星连续遥感跟踪油污范围和溢油扩散方向，确定最佳溢油清除方法。美国发射了搭载modis(moderate-resolutionimagingspectroradiometer)中分辨率成像光谱仪的人造卫星—terra(1999年12月)和aqua(2002年4月)卫星。modis传感器是卫星上唯一能将实时观测数据通过x波段向全世界直接广播，并可以免费接收数据并无偿使用的星载仪器。几何分辨率为250米、500米和1000米。目前仍然是世界上使用最广泛的卫星传感器数据。另外，由于海面水体部分和溢油区域对不同波段的反射率存在差异，因此，选择合适的波段并结合图像处理的方法就可以更有效地对溢油区域进行突出。

本发明是通过波段选取与图像处理的一些方法从遥感数据中得到海面溢油区的轮廓。首先通过海面水体部分与溢油区域的反射率差异来选取能够更好区分水体与溢油区的波段；然后利用所选波段来对自定义的参数进行计算来突出两者之间的差异；最后将计算得到的参数视为图像的灰度值来得到灰度图像，并通过相应的处理方法提取出溢油区的轮廓。

技术实现要素：

本发明的目的是从modis遥感数据中选择合适的波段并通过相关参数的计算构建灰度图像，使用图像处理的方法得到溢油区的轮廓。首先通过海面水体部分与溢油区域的反射率差异来对波段进行选取，使其能够更好地反映这两部分的差异；然后利用所选波段来对自定义的参数进行计算；最后将计算得到的参数视为图像的灰度值来得到灰度图像，并通过相应的处理方法提取出溢油区的轮廓。

本发明采用的技术方案是：

第一步：下载modis数据，使用envi软件对读入的数据进行处理，通过以下步骤实现：

(a)使用envi软件读入modis遥感数据，并分别将第1、4、3波段的数据对应为rgb模型中的r、g和b分量的值，进而得到遥感图像irgb；

(b)根据各区域是否有溢油发生这一先验知识选取遥感图像irgb中一块典型的水体区域和溢油区，得到各波段水体反射率的平均值rsea和油区反射率的平均值roil，并通过式(1)计算各波段的c值；

(c)将modis遥感数据转化为tiff格式以便下一步的处理；

第二步：使用matlab读入tiff文件,将包含数据内容的三维矩阵记为m；

第三步：通过式(2)计算坐标(x,y)处ndoi的值ndoi(x,y)，其中chcmax(x,y)和chcmin(x,y)分别为使式(1)中c的值最大和最小的波段所对应的反射率在坐标(x,y)处的取值，然后将得到的各坐标处ndoi的取值构成一个double类型的矩阵ιndoi；

本步操作可以通过matlab中的ιndoi＝(m(:,:,ncmax)-m(:,:,ncmin))./(m(:,:,ncmax)+m(:,:,ncmin))语句来实现，其中ncmax和ncmin分别为使式(1)中c的值最大和最小的波段在三维矩阵m中坐标第三个分量的取值；

第四步：将矩阵ιndoi中各元素的取值视为归一化后的灰度值，得到图像i0；

第五步：采用式(3)对图像i0进行灰度拉伸，得到图像is，其中i0(x,y)和is(x,y)分别为图像i0和is在坐标(x,y)处的取值，a和b分别为图像i0中灰度的最小值和最大值；

第六步：将图像is进行模糊化，通过以下步骤实现：

(a)利用matlab中的g＝fspecial('gaussian',[33],5)语句构造一个大小为3×3，标准差为5的高斯模板g；

(b)将图像is与高斯模板g进行卷积，得到的模糊后的图像ib；

第七步：以阈值0.036对图像ib进行roberts算子检测，将得到的图像记为ir,本步操作可以通过matlab中的

ir＝edge(ib,'roberts',0.036)语句来实现；

第八步：对溢油区域进行提取，通过以下步骤实现：

(a)对图像ir进行二值化处理，将图像ir中小于0.5的值置为0，其余值置为1，这样就可以得到一幅二值图像，记为ib；

(b)通过matlab中的se＝strel('disk',8)语句构造半径为8的圆形结构元素se，然后采用结构元素se对图像ib进行闭操作，得到图像ic，该操作可以通过matlab中的ic＝imclose(ib，se)语句来实现；

(c)通过matlab中的iclose＝imfill(ic,'holes')语句对图像ic进行孔洞填充，得到图像iclose；

第九步：提取溢油区域的轮廓，通过以下步骤实现：

(a)将图像iclose与图像iclose1进行差分，得到图像idifference,其中iclose1为已知无溢油时相同海域数据经前八步处理后得到的图像；

(b)获取图像idifference中的最大连通分量，并提取该连通区域的轮廓，该轮廓即为溢油区轮廓。

本发明优点在于：

1、可以直接对下载的modis遥感数据进行处理；

2、能够选择合适波段，并通过相关参数计算及图像处理方法有效突出海面溢油区，并提取出溢油区的轮廓。

附图说明

图1为本发明的流程图；

图2为第四步中得到的图像i0；

图3为第五步中经过灰度拉伸得到的图像is；

图4为第七步中得到的图像ir；

图5为第八步中经过闭操作及孔洞填充得到的图像iclose；

图6为最后得到的溢油区的轮廓。

具体实施方式：

下面结合具体实例对本发明做详细说明。

第一步：下载modis数据，使用envi软件对读入的数据进行处理，通过以下步骤实现：

(a)使用envi软件读入modis遥感数据，并分别将第1、4、3波段的数据对应为rgb模型中的r、g和b分量的值，进而得到遥感图像irgb；

(b)根据各区域是否有溢油发生这一先验知识选取遥感图像irgb中一块典型的水体区域和溢油区，得到各波段水体反射率的平均值rsea和油区反射率的平均值roil，并通过式(1)计算各波段的c值；

(c)将modis遥感数据转化为.giff格式以便下一步的处理；

第二步：使用matlab读入tiff文件,将包含数据内容的三维矩阵记为m；

第三步：通过式(2)计算坐标(x,y)处ndoi的值ndoi(x,y)，其中(x,y)和(x,y)分别为使式(1)中c的值最大和最小的波段所对应的反射率在坐标(x,y)处的取值，然后将得到的各坐标处ndoi的取值构成一个double类型的矩阵ιndoi；

本步操作可以通过matlab中的ιndoi＝(m(:,:,ncmax)-m(:,:,ncmin))./(m(:,:,ncmax)+m(:,:,ncmin))语句来实现，其中ncmax和ncmin分别为使式(1)中c的值最大和最小的波段在三维矩阵m中坐标第三个分量的取值；

第四步：将矩阵ιndoi中各元素的取值视为归一化后的灰度值，得到图像i0；

第五步：采用式(3)对图像i0进行灰度拉伸，得到图像is，其中i0(x,y)和is(x,y)分别为图像i0和is在坐标(x,y)处的取值，a和b分别为图像i0中灰度的最小值和最大值；

第六步：将图像is进行模糊化，通过以下步骤实现：

(a)利用matlab中的g＝fspecial('gaussian',[33],5)语句构造一个大小为3×3，标准差为5的高斯模板g；

(b)将图像is与高斯模板g进行卷积，得到的模糊后的图像ib；

第七步：以阈值0.036对图像ib进行roberts算子检测，将得到的图像记为ir,本步操作可以通过matlab中的

ir＝edge(ib,'roberts',0.036)语句来实现；

第八步：对溢油区域进行提取，通过以下步骤实现：

(a)对图像ir进行二值化处理，将图像ir中小于0.5的值置为0，其余值置为1，这样就可以得到一幅二值图像，记为ib；

(b)通过matlab中的se＝strel('disk',8)语句构造半径为8的圆形结构元素se，然后采用结构元素se对图像ib进行闭操作，得到图像ic，该操作可以通过matlab中的ic＝imclose(ib，se)语句来实现；

(c)通过matlab中的iclose＝imfill(ic,'holes')语句对图像ic进行孔洞填充，得到图像iclose；

第九步：提取溢油区域的轮廓，通过以下步骤实现：

(a)将图像iclose与图像iclose1进行差分，得到图像idifference,其中iclose1为已知无溢油时相同海域数据经前八步处理后得到的图像；

(b)获取图像idifference中的最大连通分量，并提取该连通区域的轮廓，该轮廓即为溢油区轮廓。