一种基于非均匀栅格划分进行台风轨迹路径模拟的方法技术

本发明专利技术公开了一种基于非均匀栅格划分进行台风轨迹路径模拟的方法，包括步骤：步骤S1、获取热带气旋数据，并对热带气旋数据进行预处理；步骤S2、将研究区域划分成非均匀栅格，其中，非均匀栅格是指大小不完全相等的栅格；步骤S3、基于步骤S2划分的非均匀栅格，通过经验全路径模型模拟台风轨迹。本发明专利技术采用了非均匀栅格划分方法，对台风密集的区域进行了精细化划分，提高了划分栅格的分辨率，同时能保证每个栅格有足够的历史数据用于参数拟合，以保证模拟的台风轨迹路径更加准确，接近实际的台风轨迹。轨迹。轨迹。

【技术实现步骤摘要】
一种基于非均匀栅格划分进行台风轨迹路径模拟的方法


[0001]本专利技术涉及台风路径模拟
，具体涉及一种基于非均匀栅格划分进行台风轨迹路径模拟的方法。

技术介绍

[0002]台风灾害作为常见的自然灾害之一，严重威胁着沿海地区的生命财产安全。强风能够毁坏房屋建筑和公共设施，伴随台风的强降雨能够引起洪水以及城市内涝。了解台风的轨迹，对减少台风的灾害损失非常重要。
[0003]Vickery等提出经验全路径是一种常用的台风轨迹模拟方法，该方法将研究区域划分为均匀的栅格，其中，均匀的栅格是指将研究区域划分为若干个大小完全相同的栅格，对每个栅格里面的历史轨迹数据，进行线性拟合，获得线性方程参数。然后代入方程参数，采用线性方程对台风行进速度、台风行进方向和台风中心气压进行模拟。但是该方法的栅格划分分辨率较低，而且不能保证每个栅格能为线性拟合提供足够的历史数，导致模拟的台风轨迹路径精度不高。

技术实现思路

[0004]本专利技术为了克服以上技术的不足，提供了一种基于非均匀栅格划分进行台风轨迹路径模拟的方法。
[0005]本专利技术克服其技术问题所采用的技术方案是：一种基于非均匀栅格划分进行台风轨迹路径模拟的方法，包括步骤：步骤S1、获取热带气旋数据，并对热带气旋数据进行预处理；步骤S2、将研究区域划分成非均匀栅格，其中，非均匀栅格是指大小不完全相等的栅格；步骤S3、基于步骤S2划分的非均匀栅格，通过经验全路径模型模拟台风轨迹。
[0006]进一步地，步骤S1具体包括：步骤S1.1、从公共网站下载预设年份的热带气旋数据；步骤S1.2、删除步骤S1.1中下载的热带气旋数据中时段数小于等于10的气旋，其中，时段间隔为6小时。
[0007]进一步地，步骤S2具体包括：步骤S2.1、设研究区域为；步骤S2.2、假设将划分为个大小相等的区域、、，其中为整数且的取值范围为；步骤S2.3、判断、、内的数据数是否均小于预设阈值：若是，则不会被划分；若否，则被划分为个大小相等的区域、、；其中，预设阈值用于保证经验全路径模型模拟台风轨迹时有足够的数据对经验全路径模型公式中的参数进行拟
合；步骤S2.4、针对、、中数据数小于预设阈值的区域，区域内的线性方程的参数将采用中的数据进行拟合；针对、、中数据数大于等于预设阈值的区域，对数据大于等于预设阈值的每个区域，分别重复步骤S2.2~步骤S2.3，直至不能继续划分或划分后的区域尺寸大小小于预设的最小尺寸。
[0008]进一步地，对于台风发生多受关注的区域，样本点多，栅格划分密集；对于台风发生少受关注少的区域，样本点少，栅格划分稀疏。
[0009]进一步地，步骤S3中，经验全路径模型采用如下公式：（1）（2）其中，和分别为台风中心所在位置的经度和纬度；和分别是公式（1）和公式（2）的残差；和分别是台风在时刻和时刻的行进角度，角度值按照顺时针方向从
‑
180
°
到180
°
，设正北方向角度值为0
°
；为时刻到时刻台风行进角度的变化量；为台风在时刻的行进速度；为时刻到时刻行进速度自然对数的变化量；、、、、为公式（1）的系数；、、、、、为公式（2）的系数。
[0010]步骤S3具体包括：步骤S3.1、针对步骤S2中划分的栅格，对每个栅格分别通过公式（1）和公式（2）对系数、、、、和、、、、、进行拟合，拟合方法采用稳健线性回归中的随机采样一致性算法RANSAC；步骤S3.2、对公式（1）的残差和公式（2）的残差，采用分布进行拟合；步骤S3.3、选取模拟台风的初始状态，初始状态包括台风的初始位置、台风初始行进角度和台风的初始行进速度，代入公式（1）和公式（2），获得模拟的台风轨迹，其中，台风的初始位置指初始时台风中心所在位置的经度和纬度。
[0011]进一步地，步骤S3.1中的随机采样一致性算法RANSAC具体如下：1）从栅格内所有样本点中随机选取若干个样本点，对选取的若干个样本点进行线性回归；2）判断栅格内所有样本点中非步骤1）中的样本点是否为经验全路径模型内的点；3）如果经验全路径模型内的样本点比例达到设定的阈值，则采用最小二乘法重新对这些满足经验全路径模型的点进行线性拟合，否则，重复步骤1）和2）；其中，最小二乘法的计算公式如下：（3）式（3）中，、、均为矩阵，为矩阵的转置，为矩阵的逆矩阵，为经验全路径模型公式系数的矩阵形式，为栅格内的所有样本点的自变量的矩阵形式，为栅格内的所有样本点的因变量的矩阵形式。
[0012]本专利技术的有益效果是：1.本专利技术采用了非均匀栅格划分方法，对台风密集的区域进行了精细化划分，提
高了划分栅格的分辨率，同时能保证每个栅格有足够的历史数据用于参数拟合，以保证模拟的台风轨迹路径更加准确，接近实际的台风轨迹。
[0013]2.在参数拟合中，部分栅格进行了数据共享，如步骤S2.4中区域，当的子区域内的数据数小于预设阈值时，区域的数据共享给了该子区域，这样能保证在区域划分子区域中，即使区域中大量样本数据聚集在某一子区域，此子区域仍能继续划分为更小的子区域；如若数据不能共享，当的子区域内的数据数小于阈值时，区域不能划分。
[0014]3.栅格内的参数拟合采用了稳健线性回归，减少了线性回归过程中离群值对拟合结果的影响。
附图说明
[0015]图1为本专利技术实施例所述的台风数据的栅格划分的示意图。
[0016]图2为本专利技术实施例所述的向东方向台风数据的栅格划分的示意图。
[0017]图3为本专利技术实施例所述的向西方向台风数据的栅格划分的示意图。
[0018]图4 为本专利技术实施例所述的部分历史台风轨迹的示意图。
[0019]图5为本专利技术实施例所述的部分模拟的台风轨迹的示意图。
具体实施方式
[0020]为了便于本领域人员更好的理解本专利技术，下面结合附图和具体实施例对本专利技术做进一步详细说明，下述仅是示例性的不限定本专利技术的保护范围。
[0021]本实施例所述的一种基于非均匀栅格划分进行台风轨迹路径模拟的方法，包括步骤：步骤S1、获取热带气旋数据，并对热带气旋数据进行预处理；步骤S2、将研究区域划分成非均匀栅格，其中，非均匀栅格是指大小不完全相等的栅格；步骤S3、基于步骤S2划分的非均匀栅格，通过经验全路径模型模拟台风轨迹。
[0022]为了更好地理解上述的技术方案，下面将参照附图详细地描述本专利技术的示例性实施例，虽然附图中显示了本专利技术的示例性实施例，然而应当理解，还可以以各种形式实现本专利技术而并不局限于这里所述的实施例。这些实施例是为了使本领域技术人员更清楚和更透彻地理解本专利技术。
[0023]具体地，本实施例所述的一种基于非均匀栅格划分进行台风轨迹路径模拟的方法，包括步骤如下：步骤S1、获取热带气旋数据，并对热带气旋数据进行预处理。
[0024]进一步地，步骤S1具体包括：步骤S1.1、从公共网站下载预设年份的热带气旋数据，本实施例以从某国气象局热带气旋资料中心下载194本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于非均匀栅格划分进行台风轨迹路径模拟的方法，其特征在于，包括步骤：步骤S1、获取热带气旋数据，并对热带气旋数据进行预处理；步骤S2、将研究区域划分成非均匀栅格，其中，非均匀栅格是指大小不完全相等的栅格；步骤S3、基于步骤S2划分的非均匀栅格，通过经验全路径模型模拟台风轨迹。2.根据权利要求1所述的基于非均匀栅格划分进行台风轨迹路径模拟的方法，其特征在于，步骤S1具体包括：步骤S1.1、从公共网站下载预设年份的热带气旋数据；步骤S1.2、删除步骤S1.1中下载的热带气旋数据中时段数小于等于10的气旋，其中，时段间隔为6小时。3.根据权利要求1所述的基于非均匀栅格划分进行台风轨迹路径模拟的方法，其特征在于，步骤S2具体包括：步骤S2.1、设研究区域为；步骤S2.2、假设将划分为个大小相等的区域、、，其中为整数且的取值范围为；步骤S2.3、判断、、内的数据数是否均小于预设阈值：若是，则不会被划分；若否，则被划分为个大小相等的区域、、；其中，预设阈值用于保证经验全路径模型模拟台风轨迹时有足够的数据对经验全路径模型公式中的参数进行拟合；步骤S2.4、针对、、中数据数小于预设阈值的区域，区域内的线性方程的参数将采用中的数据进行拟合；针对、、中数据数大于等于预设阈值的区域，对数据大于等于预设阈值的每个区域，分别重复步骤S2.2~步骤S2.3，直至不能继续划分或划分后的区域尺寸大小小于预设的最小尺寸。4.根据权利要求3所述的基于非均匀栅格划分进行台风轨迹路径模拟的方法，其特征在于，对于台风发生多受关注的区域，样本点多，栅格划分密集；对于台风发生少受关注少的区域，样本点少，栅格划分稀疏。5.根据权利要求1所述的基于非均匀栅格划分进行台风轨迹路径模拟的方法，其特征在于，步骤S3中，经验全路径模型采用如下公式：（1）（2）其中，和分别为台风中心所在位置的经度和纬度；和分别是公式（1）和公式（...

【专利技术属性】
技术研发人员：王乃玉，黄秀兵，
申请(专利权)人：浙江大学，
类型：发明
国别省市：