台风登陆预测方法技术

本发明专利技术提出一种台风登陆预测方法，包括：获取特定区域的台风活动记录；从台风活动记录中筛选出登陆台风和未登陆台风；筛选出登陆台风在登陆特定区域之前的预设时间段，登陆台风所在的区域，并利用聚类算法将区域化为多个子区域；针对多个子区域中的每个区域，筛选出进入每个区域的所有登陆台风观测点和未登陆台风观测点，并分别获取登陆台风观测点和未登陆台风观测点的特征因子；利用每个区域获取的登陆台风观测点和未登陆台风观测点的特征因子，利用分类算法生成每个区域的台风登陆判断标准；利用每个区域获取的登陆台风观测点和未登陆台风观测点，对每个区域的台风登陆判断标准进行验证，以获取台风登陆预测的概率。本方法准确性高，实时性好。

【技术实现步骤摘要】
台风登陆预测方法
本专利技术涉及台风预警
，尤其涉及一种台风登陆预测方法。
技术介绍
台风预报评估工作在国内外一直被广泛关注，由于台风具有随机性强、影响范围广、爆发能量巨大的特点，以往的研究主要集中于台风强度及路径预测方面。美国国家飓风中心、日本东京台风中心和中国中央气象台等多个机构一直对不同区域的热带气旋的预报情况进行统计分析、总结和改进。目前台风预测的主要方法有：变维分形模型预测，基于地理信息系统(GIS)的空间分析预测、气候持续(CLIPER)模式预测，基于遗传算法的人工神经元网络预测，数值预报模式并行化预测，卫星导风技术预测和分形分布模型预测等。以2004年为例，美国飓风中心分别预测预测大西洋飓风24/48/72小时位置平均距离误差分别为106/187/280公里。日本东京台风中心发布的台风位置预测平均误差分别为125/243/355公里。中国中央气象台24/48/72小时台风位置预报平均误差分别为120/215/326公里。由此可见台风路径预报误差较大的缺点虽有所改善，但仍旧无法满足各行业对台风预报精度的需求。而且以往的台风预测方法不能判断台风是否登陆，只能预测台风路径及相应特征因子。
技术实现思路
本专利技术旨在至少在一定程度上解决相关技术中的技术问题之一。为此，本专利技术的目的在于提出一种准确性高、实时性强的台风登陆预测方法。为了实现上述目的，本专利技术的实施例提出一种台风登陆预测方法，包括以下步骤：获取特定区域的台风活动记录；从所述台风活动记录中筛选出登陆台风和未登陆台风；筛选出所述登陆台风在登陆所述特定区域之前的预设时间段，所述登陆台风所在的区域，并利用聚类算法将所述区域化为多个子区域；针对所述多个子区域中的每个区域，筛选出进入所述每个区域的所有登陆台风观测点和未登陆台风观测点，并分别获取所述登陆台风观测点和所述未登陆台风观测点的特征因子；利用所述每个区域获取的所述登陆台风观测点和所述未登陆台风观测点的特征因子，利用分类算法生成所述每个区域的台风登陆判断标准；利用所述每个区域获取的所述登陆台风观测点和所述未登陆台风观测点，对所述每个区域的台风登陆判断标准进行验证，以获取所述台风登陆预测的概率。根据本专利技术实施例的台风登陆预测方法，依据特定区域的台风登陆记录，结合聚类分析和分类分析等的数据挖掘方法，计算出特定区域的台风的登陆的判断标准和预测概率，目标性强。准确性高，实时性好。在一些示例中，所述登陆台风定义为：台风中心与特定区域边界的最小距离小于或等于台风的影响半径。在一些示例中，所述最小距离的计算公式为：假设已知两点的经纬度坐标分别为d1＝(E1,N1),d2＝(E2,N2)，则所述已知两点之间的最小距离为：|d1d2|＝R·arccos[cosN1·cosN2·cos(E1-E2)+sinN1·sinN2]其中，E1,E2代表经度，N1,N2代表纬度，R为地球的平均半径。在一些示例中，所述聚类算法为K-means聚类算法。在一些示例中，所述特征因子包括：强度等级、纬度、经度、中心最低气压、2分钟平均近中心最大风速、纬度迁移速度和经度迁移速度。在一些示例中，所述纬度迁移速度和经度迁移速度的计算公式为：假设当前观测时刻和所述当前观测时刻之前的两个观测时刻的台风中心经纬度分别为(Et,Nt),(Et-1,Nt-1)和(Et-2,Nt-2)，则所述当前时刻的所述经度迁移速度和所述纬度迁移速度分别如下所示：其中，R为地球的平均半径，所述经度迁移速度和所述纬度迁移速度的单位为：km/h。在一些示例中，所述分类算法为分类回归树算法。在一些示例中，所述台风登陆预测的概率包括：正判概率、虚警概率和漏检概率。在一些示例中，所述正判概率、虚警概率和漏检概率的计算公式如下：假设某区域中的登陆台风和未登陆台风的观测点的数目分别为Z1和Z2，将登陆台风观测点的台风登陆的数目设为M1,将未登陆台风观测点未登陆台风的数目为M2，则所述正判概率、虚警概率和漏检概率分别为：本专利技术附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本专利技术的实践了解到。附图说明图1是根据本专利技术一个实施例的台风登陆预测方法的流程图；图2是本专利技术一个实施例的台风登陆前48小时所处区域划分示意图；图3是本专利技术一个实施例的一个区域的台风登陆判断标准流程图；和图4是本专利技术一个实施例的另一个区域的台风登陆判断标准流程图。具体实施方式下面详细描述本专利技术的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本专利技术，而不能理解为对本专利技术的限制。本专利技术的实施例提出一种台风登陆预测方法，如图1所示的根据本专利技术一个实施例的台风登陆预测方法的流程图，该方法包括：获取特定区域的台风活动记录；从台风活动记录中筛选出登陆台风和未登陆台风；筛选出登陆台风在登陆特定区域之前的预设时间段，登陆台风所在的区域，并利用聚类算法将区域化为多个子区域；针对多个子区域中的每个区域，筛选出进入每个区域的所有登陆台风观测点和未登陆台风观测点，并分别获取登陆台风观测点和未登陆台风观测点的特征因子；利用每个区域获取的登陆台风观测点和未登陆台风观测点的特征因子，利用分类算法生成每个区域的台风登陆判断标准；利用每个区域获取的登陆台风观测点和未登陆台风观测点，对每个区域的台风登陆判断标准进行验证，以获取台风登陆预测的概率。具体过程如下：步骤S101，获取特定区域的台风活动记录。步骤S102，从台风活动记录中筛选出登陆台风和未登陆台风。具体地，在本专利技术的一个实施例中，登陆台风定义为：台风中心与特定区域边界的最小距离小于或等于台风的影响半径。其中，最小距离的计算公式为：假设已知两点的经纬度(弧度制)坐标分别为d1＝(E1,N1),d2＝(E2,N2)，则已知两点之间的最小距离为：|d1d2|＝R·arccos[cosN1·cosN2·cos(E1-E2)+sinN1·sinN2]其中，E1,E2代表经度，N1,N2代表纬度，R为地球的平均半径。步骤S103，筛选出登陆台风在登陆特定区域之前的预设时间段，登陆台风所在的区域，并利用聚类算法将区域化为多个子区域。具体地，在本专利技术的一个实施例中，筛选出登陆台风在登陆前48小时所在的区域，然后利用K-means聚类算法，将该区域化为若干个子区域。子区域的数目根据实际情况而定，目的是缩小研究范围，提高登陆判断准则的正判概率。步骤S104，针对多个子区域中的每个区域，筛选出进入每个区域的所有登陆台风观测点和未登陆台风观测点，并分别获取登陆台风观测点和未登陆台风观测点的特征因子。具体地，在本专利技术的一个实施例中，针对每个子区域，筛选出进入该子区域的所有登陆台风观测点和未登陆台风观测点，分别得到这些观测点的7个特征因子，包括强度等级(SG)、纬度(LAT)、经度(LON)、中心最低气压(AP)、2分钟平均近中心最大风速(WS)、纬度迁移速度(LATMV)和经度迁移速度(LONMV)。其中，经度迁移速度和纬度迁移速度的计算方法如下：假设当前观测时刻和当前观测时刻之前的两个观测时刻的台风中心经纬度(弧度制)分别为(Et,Nt),(Et-1,本文档来自技高网...

【技术保护点】
一种台风登陆预测方法，其特征在于，包括以下步骤：获取特定区域的台风活动记录；从所述台风活动记录中筛选出登陆台风和未登陆台风；筛选出所述登陆台风在登陆所述特定区域之前的预设时间段，所述登陆台风所在的区域，并利用聚类算法将所述区域化为多个子区域；针对所述多个子区域中的每个区域，筛选出进入所述每个区域的所有登陆台风观测点和未登陆台风观测点，并分别获取所述登陆台风观测点和所述未登陆台风观测点的特征因子；利用所述每个区域获取的所述登陆台风观测点和所述未登陆台风观测点的特征因子，利用分类算法生成所述每个区域的台风登陆判断标准；利用所述每个区域获取的所述登陆台风观测点和所述未登陆台风观测点，对所述每个区域的台风登陆判断标准进行验证，以获取所述台风登陆预测的概率。

【技术特征摘要】
1.一种台风登陆预测方法，其特征在于，包括以下步骤：获取特定区域的台风活动记录；从所述台风活动记录中筛选出登陆台风和未登陆台风；筛选出所述登陆台风在登陆所述特定区域之前的预设时间段，所述登陆台风所在的区域，并利用聚类算法将所述区域化为多个子区域；针对所述多个子区域中的每个区域，筛选出进入所述每个区域的所有登陆台风观测点和未登陆台风观测点，并分别获取所述登陆台风观测点和所述未登陆台风观测点的特征因子；利用所述每个区域获取的所述登陆台风观测点和所述未登陆台风观测点的特征因子，利用分类算法生成所述每个区域的台风登陆判断标准；利用所述每个区域获取的所述登陆台风观测点和所述未登陆台风观测点，对所述每个区域的台风登陆判断标准进行验证，以获取所述台风登陆预测的概率。2.如权利要求1所述的方法，其特征在于，所述登陆台风定义为：台风中心与特定区域边界的最小距离小于或等于台风的影响半径。3.如权利要求2所述的方法，其特征在于，所述最小距离的计算公式为：假设已知两点的经纬度坐标分别为d1＝(E1,N1),d2＝(E2,N2)，则所述已知两点之间的最小距离为：|d1d2|＝R·arccos[cosN1·cosN2·cos(E1-E2)+sinN1·sinN2]其中，E1,E2代表经度，N1,N2代表纬度，R为地球的平均半径。4.如权利要求1所述的方法，其特征在于，所述聚类算法为K-means聚类算法。5.如权利要求1所述的方法，其特征在于，所述特征因子包括：强度等级、纬度、经度、中心最低气压、2分钟平均近中心最大风速、纬度迁移速度和经度迁移速度。6.如权利要求5所述的方法，其特征在于，所述纬度迁移速度和经度迁移速度的计算公式...

【专利技术属性】
技术研发人员：高文胜，周瑞旭，张博文，符祥干，陈钦柱，黄松，梁亚峰，
申请(专利权)人：清华大学，海南电网公司，
类型：发明
国别省市：北京;11