基于T样条的气象传真图等值线插值方法技术

本发明专利技术公开了一种基于T样条的气象传真图等值线插值方法，首先对气象传真图进行预处理，得到仅含等值线的图像；然后采用逐点筛选法对处理后的等值线进行采样，提取数据点；根据数据点建立T网格；最后使用T样条局部节点插值算法对图像进行插值，并将插值后的点拟合成曲线，形成插值后的等值线。本发明专利技术提高了运算速度和局部插值的能力，减少了误差的堆积和数据冗余。

【技术实现步骤摘要】
基于T样条的气象传真图等值线插值方法
本专利技术属于计算机图形学领域，特别是一种基于T样条的气象传真图等值线插值方法。
技术介绍
气象传真图(MeteorologicalFacsimileChart)即传真天气图是一种通过气象传真机接收的简单直观的天气图，它通过在一张特制的底图上填出一定区域内各观测点同一时刻的气象数据，从而反映出该区域的天气状况及变化趋势。作为海上航行的重要工具，气象传真图可用于观察和监视航行海域内海况和天气系统的发生、演变和移动等情况，对船舶远离恶劣天气、选择合理航线、安排有效作业、保证航行安全等有重要意义。现阶段气象传真图的数字化程度不高，无法显示图中所有地方的天气状况，因此需要进行插值，在现有数据的基础上预测未知数据。常用的空间插值方法有：Delaunay三角剖分法、最小曲率法、自然邻点插值法等，然而这些插值算法应用范围较小，只有在较为均匀的点分布下才能取得好的插值效果，并且计算过程复杂，对硬件要求高。Norman提出一种将不同等值线片段排列成连续线的算法，该算法适用于2D与3D物体的插值，但是需要大量的线段集。Gill提出一种在两个平行切片之间进行插值一个分段线性曲面的方法，该方法可以应用到医学成像、地理信息系统，但是在输入多边形的尖角附近产生不需要的长峰值，鲁棒性较差。Chai采用梯度控制部分微分方程同时实现了等值线的插值与有理曲面的重建，但该算法的计算较为复杂。
技术实现思路
本专利技术的目的在于提供一种基于T样条的气象传真图等值线的插值方法，提高了计算精度，降低了计算成本和操作难度，易于实现。实现本专利技术目的的技术解决方案为：一种基于T样条的气象传真图等值线插值方法，包括以下步骤：步骤1、对气象传真图进行预处理，得到仅含等值线的图像；步骤2、采用逐点筛选法对处理后的等值线进行采样，提取数据点；步骤3、根据数据点建立T网格；步骤4、基于T样条插值算法对图像进行插值拟合形成插值后的等值线；本专利技术与现有的技术相比，其显著的优点为：(1)控制网格的生成比较简单；(2)能够轻易实现控制网格的局部细化，多张样条曲面之间的高效拼接；(3)T样条还有一个特别的性质就是除去奇异点之外的每一个控制顶点都有一个基函数与之对应，从而T样条曲面在除奇异点外有一个统一的表达式，这为曲面的重建等相应的分析和处理带来了极大的便利。附图说明图1是本专利技术的方法流程图。图2是实施例的原始气象传真图。图3是实施例预处理后的仅含等值线的图像。图4是本专利技术提取原始数据点的示意图。图5是实施例采样后的等值线图。图6是本专利技术建立的T样条。图7是实施例插值后的等值线示意图。具体实施方法结合图1，一种基于T样条的气象传真图等值线插值方法，包括以下步骤：步骤1、对气象传真图进行预处理，得到仅含等值线的图像；步骤2、采用逐点筛选法对处理后的等值线进行采样，提取数据点；步骤3、根据数据点建立T网格；步骤4、基于T样条插值算法对图像进行插值拟合形成插值后的等值线；进一步的，步骤1进行的预处理包括二值化、底图去除、去噪、断点补全和线段细化操作，所述线段细化操作是指对线段进行单像素细化。进一步的，步骤2针对非封闭的等值线以线的端点为起点进行采样，针对封闭的等值线，任选一点为起点进行采样。提取数据点的方法为：设p(i)为筛选的起点，坐标为(xi,yi)，p(j)为待筛选点，坐标为(xj,yj)，p(k)为p(i)与p(j)之间的点，坐标为(xk,yk)，若于p(i)与p(j)之间的每一个点p(k)到p(i)与p(j)连线的距离小于阈值i，则将p(j)作为采样点，并将其作为下一轮待筛的选点。进一步的，步骤3根据数据点建立T网格的方法为：T网格的生成过程是通过数据建立平面网格，再在二维参数平面上利用四叉树细分为带T-交点的矩形T网格，由此确定T网格所需的边权值(节点间距)，从而确定T网格的拓扑结构。步骤3.1、以提取的数据为样本，建立平面网格；步骤3.2、在已经建立好的平面网格上定义四叉树细分的最大次数n以及最大容量值T，计算每个四边块里面包含的点数，与最大容量值进行比较之后再次进行细分直到所有小四方块所含点数都小于nT；步骤3.3、由四叉树细分的结果，记录下T网格的控制顶点数目，根据节点间距的定义给T网格的每条边都赋予一个非负值(节点间距)。进一步的，步骤4进行插值的方法为：在插值之前先输入三角网格曲面ST＝[V,E,T]，其中V＝{p1,p2...,pm}为顶点集，E＝{e1,e2...,em}为边集，T＝{t1,t2...,tm}为三角面片集。为了确定T网格控制点集合，设Q＝{q1,q2...,qm}为重建三角网格曲面ST＝[V,E,T]的T网格中待求的控制顶点集,V＝{p1,p2...,pn}为采样点集，m<n。由此可以定义T网格曲面为：其中Bi(ui,vi)是由T网格的节点间距矢量确定的基函数，wi是第i个控制点出的系数。为了求Q＝{q1,q2...,qm}需要建立目标方程：使其最小。其中(ui,vi)是三角网格曲面顶点{pi}在T样条中的参数值。通常地，为了曲面的光滑，在曲面重建过程中加一个附加光顺能量函数：其中Suu2(ui,vi)、Suv2(ui,vi)和Svv2(ui,vi)分别是和从而得到一个新的目标函数：E(q1,q2,...,qm)＝Efit+σEfair根据实验经验σ一般取值为0.01到0.1之间，重建曲面的目标就是在使E最小化时，解得控制顶点集Q，为求目标方程等价于求解方程组ATAQ＝ATP,矩阵A＝(aki)的元素为：ak,i＝Bi(uk,vk)再加上光顺项有：(ATA+σ(AuuTAuu+2AuvTAuv+AvvTAvv))Q＝ATP其中Auu,Auv,Avv的元素分别为基函数的二阶偏导数，解此线性方程可得到相应的T样条曲面控制点集Q，进而得到T样条曲面。步骤4.1、在T样条网格中插入一个新的控制点；步骤4.2、如果存在一个基函数不能从当前T网格中得出其相应的节点值，则对此其函数进行相应的节点插入步骤4.3、如果存在一个基函数，它有一个节点在T网格中找不到相应的控制顶点与之对应，则对此网格进行相应的控制顶点插入。步骤4.4、重复步骤4.2和4.3两步的操作直到没有新的节点与控制顶点的插入。最后根据新插入的控制点和原有的数据进行拟合形成新的等值线。下面结合实施例和附图对本专利技术进行详细说明。实施例本专利技术采用了一种基于T样条的等值线插值模型，该模型利用T样条的特别性质能够在插值过程中除奇异点外都有一个统一的表达式，能为曲面重建带来极大的便利，并用T样条局部插值算法细分等值线对其进行更好的插值。如图1所示，基于T样条气象传真图等值线的插值方法，包括以下步骤：第一步，对气象传真图进行预处理，仅保留等值线。图2表示的是日本JMH在2018年4月发布的原始的一张气象传真图，从图中可见，除所需等值线外还包含经纬线、底图、题图例、锋线等冗余信息。因此需要对其进行二值化、底图去除、断点补全、细化等处理得到仅包含等值线的图像，如图3所示。这些等值线是关于同一气压的海拔高度，同一根线上的高度是一样的。本专利技术以像素点在图中的坐标为采样点的坐标信息，即x,y值，以等值线上的点的实际海拔高度值作为采样点的高程信息，即z值，图中每一本文档来自技高网...

【技术保护点】
1.一种基于T样条的气象传真图等值线插值方法，其特征在于，包括以下步骤：步骤1、对气象传真图进行预处理，得到仅含等值线的图像；步骤2、采用逐点筛选法对处理后的等值线进行采样，提取数据点；步骤3、根据数据点建立T网格；步骤4、基于T样条插值算法对图像进行插值和拟合，形成插值后的等值线。

【技术特征摘要】
1.一种基于T样条的气象传真图等值线插值方法，其特征在于，包括以下步骤：步骤1、对气象传真图进行预处理，得到仅含等值线的图像；步骤2、采用逐点筛选法对处理后的等值线进行采样，提取数据点；步骤3、根据数据点建立T网格；步骤4、基于T样条插值算法对图像进行插值和拟合，形成插值后的等值线。2.根据权利要求1所述的基于T样条的气象传真图等值线插值方法，其特征在于，步骤1进行的预处理包括二值化、底图去除、去噪、断点补全和线段细化操作，所述线段细化操作是指对线段进行单像素细化。3.根据权利要求1所述的基于T样条的气象传真图等值线插值方法，其特征在于，步骤2针对非封闭的等值线以线的端点为起点进行采样，针对封闭的等值线，任选一点为起点进行采样。4.根据权利要求1所述的基于T样条的气象传真图等值线插值方法，其特征在于，步骤2提取数据点的方法为：设p(i)为筛选的起点，坐标为(xi，yi)，p(j)为待筛选点，坐标为(xj，yj)，p(k)为p(i)与p(j)之间的点，坐标为(xk，yk)，若p(i)与p(j)之间的每一个点p(k)到p(i)与p(j)连线的距离小于阈值i，则将p(j)作为采样点，并将其作为下一轮待筛的选点。5.根据权利要求1所述的基于T样条的气象传真图等值线插值方法，其特征在于，步骤3根据数据点建立T网格的方法为：T网格的生成过程是通过数据点建立平面网格，再在二维参数平面上利用四叉树细分为带T-交点的矩形T网格，由此确定T网格所需的边权值，即为节点间距，从而确定T网格的拓扑结构；具体为：步骤3.1、以提取的数据点为样本，建立平面网格；步骤3.2、在已经建立好的平面网格上定义四叉树细分的最大次数n以及最大容量值T，计算每个四边块里面包含的点数，与最大容量值进行比较之后再次进行细分直到所有小四方块所含点数都小于n*T；步骤3.3、由四叉树细分的结果，记录下T网格的控制顶点数目，根据节点间距的定义给T网格的每条边都赋予一个非负值，该非负值即为节点间距。6.根据权利要...

【专利技术属性】
技术研发人员：张芮，张雪禹，潘海朗，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏,32