基于卷积的气象传真图等值线交叉点检测方法与流程

本发明属于计算机图形学领域，特别是一种基于卷积的气象传真图等值线交叉点检测方法。

背景技术：

气象传真图(meteorologicalfacsimilechart)即传真天气图是一种通过气象传真机接收的简单直观的天气图，它通过在一张特制的底图上填出一定区域内各观测点同一时刻的气象数据，从而反映出该区域的天气状况及变化趋势。作为海上航行的重要工具，气象传真图可用于观察和监视航行海域内海况和天气系统的发生、演变和移动等情况，对船舶远离恶劣天气、选择合理航线、安排有效作业、保证航行安全等有重要意义。

在对气象传真图中等值线进行一系列的图像处理技术中，常常包括等值线的矢量化、等值线插值、等值线的断线修复、等值线的涡度线检测等，而在实现这些技术前需要得到完整平滑的等值线，从而需要对气象传真图进行预处理，包括底图去除、特殊文字去除、图像去噪等，针对于等值线上的毛刺噪声信息，通过检测出这些毛刺噪声信息与等值线形成的交叉点，再通过已有的方向链码的方法去除毛刺可以达到很好的效果，因此检测出准确的等值线交叉点具有重要意义。当前交叉点的检测方法中，有的基于模板匹配的方法进行检测，有的基于角点检测的方法进行检测，有的通过神经网络进行训练的方式进行检测，但是这些方法存在大量的误检和漏检，同时通过训练集的方法需要针对每类图像建立庞大的数据集而带来巨大的工作量，从而导致检测效率低、计算较为复杂。

技术实现要素：

本发明的目的在于提供一种基于卷积的气象传真图交叉点检测方法，能够在复杂的等值线中保持较高的检测精度，具有较高的鲁棒性，易于实现。

实现本发明目的的技术解决方案为：一种基于卷积的气象传真图交叉点检测方法，包括以下步骤：

步骤1，对气象传真图进行预处理，得到细化的等值线图；

步骤2，将细化的等值线与卷积核进行卷积运算，通过设定阈值以及判断图像的八邻接节点个数得到分类的候选交叉点；

步骤3，计算候选交叉点与每类中心点的最小欧式距离得到最终的交叉点；

步骤4，计算畸变交叉点所在等值线的宽度与最大角度，将满足阈值的点进行合并；

步骤5，在等值线上将交叉点进行标记。

本发明与现有技术相比，其显著优点为：(1)通过将细化的气象传真图与卷积核进行卷积计算方法简单，运算速度快，对于多种复杂情形的等值线，不会出现误检测以及漏检，通过设定阈值可以快速的检测到交叉点的位置坐标；(2)对于细化过程中存在畸变的交叉点，通过判断原等值线上的宽度以及细化等值线的最大角度来合并具有较高的严谨性，同时计算速度快，具有较高的识别率，将两种方法结合进行判定，有较高的检测效率，为后来的等值线的处理带来了极大的方便。

附图说明

图1是本发明基于卷积的气象传真图等值线交叉点检测方法流程图。

图2是本发明实施例的原始气象传真图。

图3是本发明实施例的预处理后细化等值线图。

图4是本发明实施例的等值线交叉点样例图。

图5是本发明实施例的等值线卷积样例图。

图6是本发明实施例的交叉点畸变示意图。

图7是本发明实施例的等值线的宽度以及最大角度示意图。

图8是本发明实施例的等值线交叉点检测结果图。

具体实施方式

结合图1，一种基于卷积的气象传真图等值线交叉点检测方法，包括以下步骤：

步骤1：对气象传真图进行预处理，得到细化的等值线图；

步骤2：将细化的等值线与卷积核进行卷积运算，通过设定阈值以及判断图像的八邻接节点个数得到分类的候选交叉点；

步骤3：计算候选交叉点与每类中心点的最小欧式距离得到最终的交叉点；

步骤4：计算畸变交叉点所在等值线的宽度与最大角度，将满足阈值的点进行合并；

步骤5：在等值线上将交叉点进行标记。

进一步的，步骤1进行的预处理包括底图去除、经纬海岸线去除、特殊字符去除、小面积噪声去除，图像二值化、细化，其中所述小面积噪声去除是指利用连通域的面积去除部分噪声。

进一步的，步骤2中首先对细化的等值线采用卷积核为进行卷积运算，则原细化图像中每个像素点的卷积值其中ai(j)为像素点ai的第j个邻接点的像素值，再将得到的卷积值与细化图像做点积，得到对应细化部分的卷积值。然后设定卷积值阈值，将满足卷积值阈值的像素点作为细化等值线的候选交叉点，对这些候选交叉点进行分类，若是在一个八邻域内，则把这些候选交叉点标记为一类。

进一步的，步骤3中对同一类的中心点选取到这些点中心距离最小的点(xl，yl)＝min{dist((xi，yi)，(xm，ym))}，其中(xi，yi)表示每一类中像素点的坐标，删去同一类中的其他点，得到最终的交叉点。

进一步的，步骤4中对于细化过程中造成的交叉点的畸变，按照两点间的欧式距离进行分组，将满足最小距离阈值的两个点划分成一组，同时这两个点需要同时满足两个条件：1.两点属于同一连通域，即两像素点之间存在一条等值线；2.两点间的阈值大小满足：d(fi，fj)＜k，其中fi、fj为畸变交叉点的坐标，k为距离阈值，对于满足以上两个约束条件的点，由于在原气象传真图等值线上产生交叉的部分的等值线的宽度以及畸变后的等值线的角度是一致的，利用细化等值线上各像素点的信息计算出交叉点各分叉线的方向，设pk为当前像素点，pk+s为细化等值线上距离pk步长为s的像素点，则当前像素点pk的方向为：其中通过细化图像上像素点的方向找到每条等值线的边缘，利用左右边界的像素空间坐标信息可以计算出等值线的宽度。像素点pk宽度计算公式为：其中为左边界的行列坐标；为右边界的行列坐标。再通过在两个分叉线方向上选取k1、k2个像素点长度的坐标(xk1，yk1)，(xk2，yk2)；计算分叉线选取的像素点的长度以及两点与交叉点的长度；

根据余弦定理求出每两个分支间的角度

θ＝arccos((a²+b²-c²)/2ab)

当两分支宽度以及最大角度满足的点进行合并，其中角度和宽度阈值wk，θk均为经验值。最终将交叉点在等值线上进行标记。

下面结合实施例和附图对本发明进行详细说明。

实施例

本实施例提出一种基于卷积的气象传真图等值线交叉点检测方法，通过将细化图像与卷积核进行卷积运算，在阈值条件的约束下可以快速的检测到交叉点的位置坐标，然后对于畸变的交叉点利用宽度和最大角度的判别方法进行合并，最终得到准确的气象传真图等值线交叉点。如图1所示，基于卷积的气象传真图等值线交叉点检测方法，包括以下步骤：

第一步，对气象传真图进行预处理，保留等值线关键信息。图2是日本jmh在2018年10月份发布的原始的一张气象传真图，从图中可见，除等值线外还包含经纬海岸线、底图、题图例、锋线、特殊字符等冗余信息。因此，要对其进行底图去除、图像二值化、细化、小面积噪声去除，其中所述小面积噪声去除是指利用连通域的面积去除部分噪声，最后得到细化的等值线图像，如图3所示。

第二步，对预处理后的细化等值线与卷积核为进行卷积运算，则原细化图像中每个像素点的卷积值其中ai(j)为像素点ai的第j个邻接点的像素值，得到的卷积值与细化图像做点积，得到对应细化部分的卷积值，然后设定卷积值阈值，将满足卷积值阈值的像素点作为细化等值线的候选交叉点，对这些候选交叉点进行分类，若是在一个八邻域内，则把这些候选交叉点标记为一类，其中等值线交叉点的样例图如图4所示，以及交叉点与卷积核卷积过程如图5所示。

第三步，对同一类的中心点选取到这些点中心距离最小的点(xl，yl)＝min{dist((xi，yi)，(xm，ym))}，其中(xi，yi)表示每一类中像素点的坐标，删去同一类中的其他点，得到最终的交叉点。同时对于细化过程中造成的交叉点的畸变，按照两点间的欧式距离进行分组，将满足最小距离阈值的两个点划分成一组，同时这两个点需要同时满足两个条件：1.两点属于同一连通域，即两像素点之间存在一条等值线2.两点间的阈值大小满足：d(fi，fj)＜k，其中fi、fj为畸变交叉点的坐标，k为距离阈值，对于满足以上两个约束条件的点，由于在原气象传真图等值线上产生交叉的部分的等值线的宽度以及畸变后的等值线的角度是一致的，利用细化等值线上各像素点的信息计算出交叉点各分叉线的方向，设pk为当前像素点，pk+s为细化等值线上距离pk步长为s的像素点，则当前像素点pk的方向为：其中通过细化图像上像素点的方向找到每条等值线的边缘，利用左右边界的像素空间坐标信息可以计算出等值线的宽度。像素点pk宽度计算公式为：其中为左边界的行列坐标；为右边界的行列坐标。再通过在两个分叉线方向上选取k1、k2个像素点长度的坐标(xk1，yk1)，(xk2，yk2)；计算分叉线选取的像素点的长度以及两点与交叉点的长度；

根据余弦定理求出每两个分支间的角度

θ＝arccos((a²+b²-c²)/2ab)

当两分支宽度以及最大角度满足的点进行合并，其中角度和宽度阈值wk，θk均为经验值，交叉点畸变示意图如图6所示，等值线的宽度以及最大角度示意图如图7所示，图中畸变交叉点两两分叉线构成的角度分别为α1、α2、β1、β2，分叉线的宽度分别为w1、w2、w3、w4。最终检测的等值线交叉点结果如图8所示，图中检测出的交叉点已经用星号进行标记，本实施例中一共检测出445个交叉点，合并了70个畸变交叉点，能够准确定位到原始等值线上的交叉区域，为进一步的检测研究奠定了良好的基础。