技术摘要：
本发明公开一种基于二维海浪谱的海浪成分分区方法，该分区方法能够区分二维海浪谱中不同的海浪成分，并且采用了“边查找边分区，边查找边合并”的方法，在一定程度上解决了现有的分区方法中过度分割的问题，同时解决了现有分区方法过分依赖于图像分割技术、缺乏对实际  全部
背景技术：
： 海洋无时无刻不在运动，其运动形式多种多样，纷繁复杂，而海浪是海水运动最常 见的形式之一。海浪就是海水质点在它的平衡位置附近产生一种周期性的震动运动和能量 的传播，其对海上运输、沿海经济、海洋开发以及军事方面有着重要的影响。船舶在海上航 行时，需要避开大浪区；而沿海地带会经常遭受大浪带来的经济损失；对于海洋开发来说， 海浪又是一种取之不尽用之不竭的能源；对于军事来说，海浪对水面舰艇航行具有重要影 响。因此，研究海浪有很重要的意义。 由于海洋有复杂的气象条件，会产生各种不同波长不同振幅的波浪，而这些波浪 往往混合在一起，因此仅仅通过观测很难区分不同的海浪成分。二维海浪谱可以描述海浪 的总能量在频率和方向上的分布，其包含了很多海浪能量成分的不同信息和特征。目前，卫 星和各种现场观测仪器测量的二维海浪谱数据日益增多，通过二维海浪谱分区算法，我们 可以将这些不同成分的波浪能量区分开，从而为海浪模式提供准确的数据支撑，进一步优 化海浪预报精度。现有的海浪谱分区方法主要存在过度分割问题，并且其主要依赖于图像 分割技术，缺乏对实际二维海浪谱物理意义的考虑、对滤波较为依赖以及过滤噪声的同时 会失去部分数据的问题。
技术实现要素：
： 为解决上述问题，本发明提供一种基于二维海浪谱的海浪成分分区方法。 本发明提出的技术方案为： 一种基于二维海浪谱的海浪成分分区方法，包括如下具体步骤： 步骤一、获取目标海域的二维海浪谱，并确定所述二维海浪谱的限制关系； 步骤二、剔除所述二维海浪谱中异常部分，并判断所述二维海浪谱是否存在异常 缺测值，如存在，则该二维海浪谱不完整，不对其进行分区；如不存在，则进入步骤三； 步骤三、对所述二维海浪谱去噪； 步骤四、在所述二维海浪谱没有分区的像素点中找出能量最大的像素点，如果该 像素点的能量为0，则该二维海浪谱已分区完毕，进入步骤九；如果该像素点的能量不为0， 则根据该像素点是否有相邻分区进行如下判断： ①如果该像素点没有相邻分区，则将该像素点划分为一个新分区； ②如果该像素点至少有1个相邻分区，则根据该像素点的相邻分区数目进行如下 判断： (a)如果该像素点仅有1个相邻分区，则判断如果该像素点合并入该相邻分区后， 该相邻分区是否满足步骤1)中限制关系：如果满足，则将该像素点合并入该相邻分区；如果 不满足，则将该像素点划分为一个新分区； 4 CN 111611975 A 说　明　书 2/5 页 (b)如果该像素点至少有2个相邻分区，则判断该像素点分别合并入每个相邻分区 后，是否满足步骤1)中限制关系： (I)如果相邻分区均不满足，则将该像素点划分为一个新分区； (II)如果仅一个相邻分区满足，则将该像素点合并入该相邻分区； (III)如果至少两个相邻分区满足，则先将该像素点与任意一个相邻分区合并，得 到一个新分区；顺次判断该新分区与下一个相邻分区是否满足步骤1)中限制关系，满足则 合并成一个新分区，并以该新分区继续顺次判断下一相邻分区，不满足则直接判断下一相 邻分区；直至判断完所有相邻分区； 步骤五、不断重复步骤四，直到所述二维海浪谱中有一个分区的总能量达到能量 阈值，将该分区记为显著分区； 步骤六、遍历所述二维海浪谱，将所有满足与显著分区合并条件的像素点合并入 该显著分区，所述满足与显著分区合并条件的像素点为与显著分区相邻、且该像素点的能 量不大于与显著分区中与其相邻的像素点的能量、且没有分区的像素点； 步骤七、不断重复步骤六，直到所述二维海浪谱中不存在满足与显著分区合并条 件的像素点，进入步骤八； 步骤八、将步骤七中显著分区中的所有像素点和所有与显著分区相邻的像素点中 的孤立部分并入显著分区； 步骤九、存储步骤八中得到的显著分区，对所述二维海浪谱其余部分不断重复步 骤四～步骤八，直到所述二维海浪谱中所有能量大于0的像素点都有分区； 步骤十、将所述二维海浪谱中所有分区按分区能量的大小进行排序，取其中能量 最大的N个分区作为该二维海浪谱海浪成分分区结果。 优选地，所述异常部分为二维海浪谱中波长小于30m和波长大于600m的部分； 优选地，述限制关系包括波长波向限制关系。 优选地，步骤三中所述去噪包括如下具体步骤： 3.1)去除观测仪器导致二维海浪谱中的热噪声； 3.2)采用高斯滤波去除二维海浪谱中的随机噪声。 优选地，所述高斯滤波的高斯内核为(3,3)，高斯核函数在X和Y方向的标准偏差均 为0.8。 优选地，步骤五中所述能量阈值为所述二维海浪谱总能量的5％-10％。 优选地，步骤五中所述能量阈值为所述二维海浪谱总能量的6.25％。 优选地，步骤八中所述孤立部分为相邻的三个及其以上像素点都属于同一分区的 像素点以及上、下相邻或者左、右相邻的像素点都属于同一分区的像素点。 优选地，步骤十中N取3-6。 本发明相比于现有技术具有如下有益效果 本发明提出的基于二维海浪谱的海浪成分分区方法，能够区分二维海浪谱中不同 的海浪成分，并且采用了“边查找边分区，边查找边合并”的方法，在一定程度上解决了现有 的分区方法中过度分割的问题，同时解决了现有分区方法过分依赖于图像分割技术、缺乏 对实际情况的考虑、对滤波较为依赖以及过滤噪声的同时会失去部分数据的问题。 本发明分区方法结果更多的受实际海浪谱影响，而现有技术对海浪谱使用的分水 5 CN 111611975 A 说　明　书 3/5 页 岭算法则侧重图像识别，本发明更能考虑海浪谱的实际物理意义。 本发明分区方法对于滤波的依赖是轻微的，这也就决定了本发明可以采用较小的 滤波参数，从而既能过滤较明显的噪声，又不会丢失太多海浪谱数据，同时轻微的噪声也有 合并算法进行处理。 附图说明： 图1为实施例中分区方法的流程示意图； 图2为实施例中步骤四的流程示意图； 图3为实施例中海浪谱波长波向限制关系示意图； 图4为实施例中采用本发明分区方法对目标海域的分区结果示意图；其中，图4(a) 为第76个box左侧二维海浪谱分区结果示意图；图4(b)为第79个box左侧二维海浪谱分区结 果示意图； 图5为实施例中采用超强滤波参数的方法对目标海域的分区结果示意图；其中，图 5(a)为第76个box左侧二维海浪谱分区结果示意图；图5(b)为第79个box左侧二维海浪谱分 区结果示意图； 图6为实施例中采用去除滤波参数的方法对目标海域的分区结果示意图；其中，图 6(a)为第76个box左侧二维海浪谱分区结果示意图；图6(b)为第79个box左侧二维海浪谱分 区结果示意图；