一种海底沙波特征识别方法

技术领域

本发明涉及近岸海洋工程动力地质地貌环境技术领域，尤其涉及一种海底沙波地貌特征的识别方法。

背景技术

海底沙波是一种韵律形且活动性比较强的海底地貌形态，沙波波峰线通常与潮流主方向垂直。作为一种特殊形态的床面形式，沙波广泛发育在河口、海湾、海峡、陆架陆坡等区域。海底沙波在波浪、潮流、风暴潮、内波等海洋动力的强烈影响下会发生演化、迁移。沙波迁移能快速改变海床形态，导致海底管线悬跨或掩埋，造成海洋平台失稳等，对海洋油气管道、海底电缆、石油平台等海洋工程设施造成严重威胁。

海底沙波作为特殊的地貌形态具有两个显著的形态特征，即沙波波高和波长。随着海洋测深技术的迅速发展，多波束测深系统、数字测深技术、侧扫声呐、高分辨率浅地层剖面仪以及卫星遥感等手段的综合使用，使获取的地形数据精度越来越高。然而由于海底地形十分复杂，从实测水深数据中直接获取较准确的沙波波高及波长较为困难。

目前针对海底沙波地形特征识别的主要手段为剖面法。然而剖面法需要人为确定所要分析的剖面，不同剖面的选择会导致计算的沙波波高、波长出现差异。在此背景下，本发明公开一种能够对沙波区域测深数据自动进行波高及波长等特征识别的方法。

发明内容

本发明所要解决的技术问题是：通过图像识别和坐标转换的方式，提供一种能够自动识别沙波波长和波高等特征的技术，充分利用实测数据，获取更多的沙波特征样本，为分析海洋工程施工区域的沙波特征提供依据。

本发明目的是由以下技术解决方案实现的：

一种海底沙波特征识别方法，利用实测水深数据，采用图像识别的方法获取沙波波峰线、波谷线，并通过线性拟合及坐标转换将沙波波峰线、波谷线进行分组，进一步利用沙波波峰线、波谷线的信息自动获取大区域的沙波波长、波高等特征信息。

进一步的，一种海底沙波特征识别方法，包括如下步骤：

1）基于大区域的实测水深数据矩阵，估算沙波波长，将大区域分割成若干个小区域矩阵；

2）对每一个小区域矩阵利用图像识别获取沙波波峰线、波谷线，将波峰线及波谷线的位置及高程信息保存到变量中；

3）将步骤2）所获取的沙波波峰线及波谷线近似平行，对子区域中的每条波峰线、波谷线进行线性拟合，获取拟合后每条直线的斜率，根据斜率范围对波峰线及波谷线进行第一次分组；

4）对步骤3）中每一组沙波直线进行坐标转换，并对波峰线和波谷线分别进行编号和分组；

5）计算步骤4）中每一个分组内沙波的波长和波高。

优选的，所述步骤2）包括如下步骤：

a.对每一个小型矩阵进行平面拟合以消除整体坡度对识别的影响；

b.利用中值滤波技术消除孤立的噪声点；

c.利用二值化技术获取沙波波峰线的位置；

d.对水深取相反数，重复步骤a、b、c可以获取沙波波谷线的位置；

e.将波峰线及波谷线的位置及高程信息保存到变量中。

优选的，所述步骤b消除孤立的噪声点在Matlab中可采用im2bw函数。

优选的，所述步骤c二值化技术采用阈值变换法将灰度图转换为二值图，在Matlab中可采用im2bw函数。

进一步的，一种海底沙波特征识别方法，所述步骤4）包括如下步骤：首先，对沙波的波峰线、波谷线进行线性拟合，根据斜率差异，将所有波峰线、波谷线分组，其中，每组波峰线、波谷线间的斜率误差均小于0.25；其次，对每组波峰线、波谷线进行坐标变换；再次，获取转换后每条线中点的y值，根据y值由小到大，将直线进行编号，并将信息分别保存在ridge、valley两个变量中；最后根据编号将沙波域分为不同的组。

优选的，所述坐标转换方法包括：根据沙波波峰线、波谷线的斜率，将原坐标系进行旋转，使斜率与变换后的坐标系的x轴平行。

进一步的，一种海底沙波特征识别方法，所述步骤5）的计算方法为：将步骤3）中每一个沙波域中相邻波峰线y值的差和相邻波谷线y值的差，即为沙波波长；同样，相邻波峰线与波谷线上的平均高程差，即为沙波波高。

本发明的有益效果

本申请通过充分利用实测数据，自动获取研究区域大量的沙波波长、波高等特征数据，进而获取大量的沙波特征数据，获取更多的沙波特征样本，可以为进一步研究沙波波长、波高、水深、水动并对沙波波峰线、波谷线自动进行分类统计，力等之间的关系提供可靠、大量的基础数据。为进一步分析沙波波长、波高、水深及水动力等之间的关系提供了基础数据，节省了时间和成本，提高了效率。为分析海洋工程施工区域的沙波特征提供依据。

本方法通过程序对沙波特征进行识别，具有数据处理量大、速度快的优点。极大提高效率的同时，排除了人工标注的主观因素影响。利用该方法获取的大量的沙波特征数据可以用于如沙波波长波高之间的关系、水深与沙波特征的关系等方面的分析。

附图说明：

图1为南海北部某沙波区域实测水深图。

图2为图1所述某一小区域的沙波区域灰度图。

图3为经过二值化后获取沙波波峰线和波谷线图（叠加在沙波区域上）。

图4为坐标转换分析沙波波峰线及波谷线的原理图。

图5为计算得到的沙波波长统计图。

图6为计算得到的沙波波高统计图。

具体实施方式

实施例1

南海北部区域发育有大量的沙波地形地貌，下面以该区域为例，结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。

步骤一、南海北部存在大量的海底沙波，经过实际水深测量获取南海某沙波区域的水深图，如图1所示。为了分析该区域沙波特征，将该区域分割成7行17列的小区域，逐个区域进行沙波特征分析，其中第5行第6列（以左下角为第1行第1列）的沙波形态见图1右上角。

步骤二、以选出的第5行第6列（该位置下文以R5C6表示）的水深数据为例，首先对该区域进行平面拟合，用实测水深减去平面拟合值，消除该区域整体坡度的影响；然后采用中值滤波技术消除孤立的噪声点（在Matlab中可采用medfilt2函数）；再次采用阈值变换法（Matlab中可采用im2bw函数）将灰度图转换为二值图，所获取的二值图见图2。对水深取相反数，采用同样的步骤可以获取波谷线的二值图。

步骤三、对图像进行填充孔洞修正后（在Matlab中可采用imfill函数）标记（在Matlab可采用bwlabel函数），即将图中的各个连通区域（即沙波波峰线）在二值图像的标记矩阵中用连续的阿拉伯数字表示，每个数字即代表图中一条连续的沙波波峰线。通过波峰线的像素坐标可以绘制波峰线（在Matlab中可依次采用sort, find, plot函数），对波谷线二值图做同样处理可获取沙波的波谷线。R5C6位置的沙波波峰线与波谷线见图3（为直观反映识别的沙波波峰线和波谷线，将波峰线和峰谷线叠加在了沙波底图上），其中实线为波峰线，虚线为波谷线。

步骤四、步骤三中获取了沙波的波峰线和波谷线，并将每条线上的坐标和高程保存在了变量中。为获取沙波的波高、波长等特征值，采用坐标变换并对沙波波峰线、波谷线进行分组。具体方法为：首先，对沙波的每条波峰线、波谷线进行线性拟合，根据斜率差异，将所有特征线分为1~3组（每组特征线间的斜率误差均小于0.25）。其次，对每组波峰线、波谷线进行坐标变换，坐标转换原理如图4所示，实线坐标系为原坐标系，黑色实线为沙波波峰线，灰色实线为沙波波谷线。根据波峰线、波谷线的斜率，将原坐标系进行旋转，使斜率与变换后的坐标系的x轴平行，转换后的坐标系见图4中虚线坐标系。再次，获取每条线段中点的y值（对应于变换后的坐标系）；根据y值由小到大，将线段进行编号，并将信息分别保存在ridge、valley两个变量中。对应于图4，ridge中线段的编号为1、2、4、7、9；valley中线段编号为3、5、6、8。最后根据编号将沙波域分为不同的组（以ridge和valley变量中连续的数进行分界，本例中ridge和valley的连续数分别为1、2和5、6，即1为1组，2-5为1组，6-9为1组，由于第1组仅一条沙波波谷线故舍去），即编号为[2，3，4，5]为一个沙波域、编号为[6、7、8、9]为第二个沙波域。

步骤五、步骤四中每一个沙波域中相邻波峰线y值（对应于变换后的坐标系）的差及相邻波谷线y值的差，即为沙波波长。R5C6位置提取的波长分别为线段4与线段2的y值差、线段5与线段3的y值差、线段8与线段6的y值差、线段9与线段7的y值差。同样的，相邻波峰线及波谷线上的平均高程差，即为沙波波高。

步骤六、经计算，本区域的沙波波长共识别出沙波波高及波长各443个。其统计特征见图5。利用该方法计算出的该区域沙波平均波长为114.89m，平均波高为1.48m；沙波波长集中在75~125m范围内，沙波波高集中在0.75~2.0m范围内。

利用该方法获取的大量的沙波特征数据可以用于如沙波波长波高之间的关系、水深与沙波特征的关系等方面的分析。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。