热带气旋路径相似程度确定方法及存储介质与流程

1.本发明涉及海洋气象预报技术领域，具体涉及一种热带气旋路径相似程度确定方法、装置及存储介质。


背景技术：

2.热带气旋，是发生在热带或副热带洋面上的低压涡旋，可以被描述为一个天气尺度的非锋面低压系统，在温暖水域上发展，有组织的对流，最大平均风速为34节或以上，在中心附近延伸超过一半，并持续至少6小时。热带气旋发生在热带北大西洋(包括加勒比海和墨西哥湾)、墨西哥西海岸外的北太平洋，偶尔向西延伸至夏威夷、北太平洋西部(包括菲律宾群岛和中国海)，孟加拉湾和阿拉伯海，马达加斯加海岸和澳大利亚西北海岸外的南印度洋，以及从澳大利亚东海岸到西140
°
的南太平洋。热带气旋是十分危险，可产生极端的风、暴雨和洪水以及破坏性的风暴潮，导致低洼沿海地区淹没，造成巨大人员和经济损失。
3.热带气旋的移动路径被认为是各种影响因素综合作用的结果。两条路径相似的热带气旋，在一定程度上能反映热带气旋的环境条件、演变过程的相似，反应了诸多影响因子综合作用的等效。热带气旋路径相似方法广泛应用于热带气旋预报、台风风暴潮预报。在台风风暴潮预报中，寻找与样本热带气旋路径相似的历史热带气旋，调阅历史热带气旋发展及引起的沿岸潮位及风暴增水实况，对开展当前样本热带气旋可能引发的风暴增水预测有十分重要的参考价值。
4.目前对于路径相似程度的评估或者路径相似热带气旋的筛选主要包括两类，一类是基于gis(地理信息系统)的空间分析功能，通过缓冲区或者叠加分析对热带气旋路径进行筛选；另一类是从路径相似出发，通过设置指标如欧氏距离、弗雷歇距离、豪斯多夫距离等方式，计算样本间的差异程度，进行相似性评估。但是上述两类方法也存在不足，主要体现在，在应用gis空间分析中的缓冲区时，缓冲区半径的取值一般只能依靠经验，但缓冲区大小则会直接影响到计算结果；在通过差异性指标计算相似性中如欧式距离、弗雷歇距离等，是通过度量路径上关键点的距离，来确定路径的相似(接近)程度。对于路径上的关键点的采样结果将直接影响计算结果，但如何合理的选取关键点采样，也没有合适的解决方案，一般人工依靠经验。此外，对于弗雷歇、动态时间归整等计算方法，在求解过程中需要递归计算，花费的计算时间也会相对偏长。


技术实现要素：

5.为至少实现路径相似程度更为准确的判断，本发明提出了一种热带气旋路径相似程度确定方法及存储介质，结合相近、重叠、连续和强度4个方面，确定热带气旋之间的相似程度。
6.本发明第一方面，提供了一种热带气旋路径相似程度确定方法，包括：
7.历史数据获取步骤，获取包括路径点数据集的热带气旋历史数据，所述路径点数据集包括热带气旋序号、时间、经度、纬度、强度、中心最低气压和热带气旋尺度，所述热带
气旋历史数据包括一条或多条历史热带气旋的数据；空间相近程度确定步骤，确定样本热带气旋路径与历史热带气旋路径的空间相近程度；重叠程度确定步骤，确定样本热带气旋路径与历史热带气旋路径的重叠程度；连续程度确定步骤，确定样本热带气旋路径与历史热带气旋路径的连续程度；强度相似程度确定步骤，确定样本热带气旋路径与历史热带气旋路径的强度相似程度；路径相似程度确定步骤，综合样本热带气旋路径与历史热带气旋路径的空间相近程度、重叠程度、连续程度和强度相似程度，确定样本热带气旋和历史热带气旋的路径相似程度。
8.可选地，所述热带气旋历史数据的数据源为来自中国气象局热带气旋资料中心的cma热带气旋最佳路径数据集和卫星分析热带气旋尺度资料，所述历史数据获取步骤包括将cma热带气旋最佳路径与卫星分析热带气旋尺度资料按热带气旋序号、时间进行拼接处理，形成数据并集，删除副中心情况，形成包括热带气旋序号、时间、经度、纬度、强度、中心最低气压、热带气旋尺度的路径点数据集。
9.可选地，所述历史数据获取步骤还包括查找历史热带气旋路径中是否存在前后相邻的两个路径点位置相同的情况，若是，则删除其中一个路径点。
10.可选地，所述历史数据获取步骤还包括为热带气旋路径设定空间约束，获取第一个和最后一个满足空间约束的路径点后，再依次取出中间路径点，包括中间不满足空间约束的路径点，按时间序列排序，作为满足空间约束的路径点集合。
11.可选地，所述空间相近程度确定步骤通过构造热带气旋椭圆轨迹，分别确定历史热带气旋路径中的路径点与样本热带气旋路径所对应的椭圆轨迹之间的匹配程度和样本热带气旋路径中的路径点与历史热带气旋路径所对应的椭圆轨迹之间的匹配程度，获得历史热带气旋路径与样本热带气旋路径的空间相近程度。
12.可选地，所述重叠程度确定步骤通过构造热带气旋椭圆轨迹，分别计算历史热带气旋路径中的路径点与样本热带气旋路径所对应的椭圆轨迹之间的归一化距离和样本热带气旋路径中的路径点与历史热带气旋路径所对应的椭圆轨迹之间的归一化距离，获得历史热带气旋路径与样本热带气旋路径的重叠程度。
13.可选地，所述连续程度确定步骤通过构造热带气旋椭圆轨迹，分别计算历史热带气旋路径中的路径点与样本热带气旋路径所对应的椭圆轨迹的第一匹配位置的序列和样本热带气旋路径中的路径点与历史热带气旋路径所对应的椭圆轨迹第一匹配位置的序列，获得历史热带气旋路径与样本热带气旋路径的连续程度。
14.可选地，所述强度相似程度确定步骤通过分别计算历史热带气旋路径中的路径点和样本热带气旋路径中的路径点的强度集合，分别构建历史热带气旋路径与样本热带气旋路径的强度矢量，通过计算两个强度矢量的余弦相似度值，确定历史热带气旋路径与样本热带气旋路径的强度相似程度。
15.可选地，还包括外接矩形相近程度确定步骤，用于确定样本热带气旋路径与历史热带气旋路径的外接矩形空间重叠程度，并根据所述外接矩形空间重叠程度的大小选择如下处理方式：
16.当外接矩形的重叠程度大于或等于第一设定值时，循环执行空间相近程度确定步骤、重叠程度确定步骤和连续程度确定步骤三次，其中路径点数据的使用，第一次执行每m个点取一个路径点，第二次执行每n个点取一个路径点，第三次执行取全部路径点，综合三
次执行所确定的空间相近程度、重叠程度和连续程度，作为最终的空间相近程度、重叠程度和连续程度，其中m大于n，n大于1；
17.当外接矩形的重叠程度小于第一设定值且大于或等于低于第一设定值的第二设定值时，循环执行两次空间相近程度确定步骤、重叠程度确定步骤和连续程度确定步骤，其中路径点数据的使用，第一次执行每a个点取一个路径点，第二次执行取全部路径点，综合两次执行所确定的空间相近程度、重叠程度和连续程度求取其均值，作为最终的相近程度、重叠程度和连续程度，其中a大于1；
18.当外接矩形的重叠程度大于0且小于第二设定值时，取全部路径点通过空间相近程度确定步骤、重叠程度确定步骤和连续程度确定步骤确定空间相近程度、重叠程度和连续程度；
19.当外接矩形不重叠的情况下，判别为历史热带气旋路径和样本热带气旋路径不相似。
20.本发明的第二方面，提供一种存储介质，所述存储介质包括存储的程序，其中，在所述程序运行时由处理器执行前述任意一种方法。
21.本发明的第三方面，提供一种热带气旋路径相似程度确定装置，包括：
22.历史数据获取模块，用于获取包括路径点数据集的热带气旋历史数据，所述路径点数据集包括热带气旋序号、时间、经度、纬度、强度、中心最低气压和热带气旋尺度，所述热带气旋历史数据包括一条或多条历史热带气旋的数据；空间相近程度确定模块，用于确定样本热带气旋路径与历史热带气旋路径的空间相近程度；重叠程度确定模块，用于确定样本热带气旋路径与历史热带气旋路径的重叠程度；连续程度确定模块，用于确定样本热带气旋路径与历史热带气旋路径的连续程度；强度相似程度确定模块，用于确定样本热带气旋路径与历史热带气旋路径的强度相似程度；路径相似程度确定模块，用于综合样本热带气旋路径与历史热带气旋路径的空间相近程度、重叠程度、连续程度和强度相似程度，确定样本热带气旋和历史热带气旋的路径相似程度。
23.本发明的第四方面，提供一种热带气旋路径相似程度确定装置，包括：处理器；以及存储器，与所述处理器连接，存储器存储有程序，其中，在所述程序运行时由处理器执行前述任意一种方法。
24.本发明提供的上述技术方案，结合相近、重叠、连续和强度4个方面，能够更为准确地确定历史热带气旋数据库中历史热带气旋个例与样本热带气旋的路径相似程度，以此为基础可进一步给出与样本热带气旋的相似排序结果，筛选排序出最为相似的热带气旋，进而再筛选背景环境场，选择与样本热带气旋相近的热带气旋，查询这些历史热带气旋的发展变化的情况以及沿岸验潮站的潮位及增水情况，结合数值预报模式输出结果、人工经验分析判断，给出对样本热带气旋可能造成风暴增水的预测结果。
25.根据下文结合附图对本发明的具体实施例的详细描述，本领域技术人员将会更加明了本发明的上述以及其他目的、优点和特征。
附图说明
26.后文将参照附图以示例性而非限制性的方式详细描述本发明的一些具体实施例。附图中相同的附图标记标示了相同或类似的部件或部分。本领域技术人员应该理解，这些
附图未必是按比例绘制的。附图中：
27.图1为本发明实施例1中热带气旋路径相似程度确定方法的基本步骤示意图；
28.图2为本发明实施例4中热带气旋路径相似程度确定装置的结构示意图；
29.图3为本发明实施例5中热带气旋路径相似程度确定装置的结构示意图。
具体实施方式
30.为了使本技术领域的人员更好地理解本发明，下面将结合附图，对本发明实施例中的技术方案进行清楚、完整地描述。需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
31.实施例1
32.本实施例提供了一种热带气旋路径相似程度确定方法，下面以图1作为参考，说明该方法的基本步骤及各步骤示例性的具体实现方式。
33.历史数据获取步骤101，获取包括路径点数据集的热带气旋历史数据，路径点数据集包括热带气旋序号、时间(utc)、经度、纬度、强度、中心最低气压和热带气旋尺度。该热带气旋历史数据包括一条或多条历史热带气旋的数据。
34.在一个具体示例中，热带气旋历史数据的数据源为来自中国气象局热带气旋资料中心的cma热带气旋最佳路径数据集和卫星分析热带气旋尺度资料。cma热带气旋最佳路径数据集包括1949年以来的西北太平洋(赤道以北，东经180
°
以西)海域热带气旋每6小时的位置和强度，包括时间(utc)、强度、纬度、经度、中心最低气压等。卫星分析热带气旋尺度资料包括西北太平洋所有被卫星捕捉到的热带气旋，空间范围涵盖赤道以北，180
°
以西。自1980年起的热带气旋逐6小时位置、强度和尺度信息，包括时间(utc)、纬度、经度、中心最低气压、中心最大风速、热带气旋尺度等。
35.在这种情况下，为获得包括热带气旋序号、时间(utc)、经度、纬度、强度、中心最低气压和热带气旋尺度在内的全部路径点数据，需先对cma热带气旋最佳路径数据集和卫星分析热带气旋尺度资料进行数据整合处理，即将cma热带气旋最佳路径与卫星分析热带气旋尺度资料按热带气旋序号、时间进行拼接处理，形成数据并集，删除副中心情况，形成包括热带气旋序号、时间(utc)、经度、纬度、强度、中心最低气压、热带气旋尺度的路径点数据集。
36.为提高后续分析计算准确性，还可对路径点数据集进行以下处理：
37.查找每条热带气旋路径中是否存在前后相邻的两个路径点位置(经纬度坐标)相同的情况，若是，则删除其中一个路径点。存在多条历史热带气旋路径时，可按热带气旋序号排序，遍历每条路径。
38.在台风风暴潮预报中对热带气旋路径相似性的分析中，重点关注的是热带气旋登陆前后的路径相似性，而非自热带气旋的生成源地开始的全部路径，可被理解为局部相似性。因而，在实际预报工作中，热带气旋路径相似程度可仅针对给定空间范围(一定空间约束)内的局部路径进行确定。在这种情况下，需首先遍历热带气旋路径，获取其中满足空间约束的路径点，形成满足空间约束的路径点集合。为保证不出现在时间序列上割裂热带气旋路径的情况，即在针对给出的空间约束，某条热带气旋路径出现路径中时间排序靠前的路径点满足空间约束要求，时间序列中间的某个或者某几个路径点不满足，但时间序列靠
后的路径点又满足空间约束的情况，还可进行以下处理：获取第一个和最后一个满足空间约束的路径点后，再依次取出中间路径点，包括中间不满足需求的路径点，按时间序列排序，作为满足空间约束的路径点集合。
39.空间相近程度确定步骤102，确定样本热带气旋路径与历史热带气旋路径的空间相近程度。
40.该步骤在具体执行中，可通过构造热带气旋椭圆轨迹，分别确定历史热带气旋路径中的路径点与样本热带气旋路径所对应的椭圆轨迹之间的匹配程度和样本热带气旋路径中的路径点与历史热带气旋路径所对应的椭圆轨迹之间的匹配程度，获得历史热带气旋路径与样本热带气旋路径的空间相近程度。每条历史热带气旋路径均可采用相同的方法处理。举例说明如下：
41.设样本热带气旋路径为cyclone
p
＝《p1,p2,...pn》，pi为热带气旋中的路径点(已按时间序列升序排列)；
42.设历史热带气旋路径为cycloneq＝《q1,q2,...qn》，qi为热带气旋中的路径点(已按时间序列升序排列)；
43.以样本热带气旋路径中的前后两个路径点为焦点，构造一个椭圆，依此类推，生成一系列的椭圆，形成椭圆序列e
p
＝《e1,e2,...en》，e为椭圆序列的椭圆；同样方式构造历史热带气旋路径的一系列椭圆，形成椭圆序列sq＝《s1,s2,...sn》，s为椭圆序列的椭圆；
44.根据样本热带气旋路径中的路径点落入历史热带气旋路径的椭圆序列的数量，和历史热带气旋路径中的路径点落入样本热带气旋路径的椭圆序列的数量，确定两条热带气旋路径的空间相近程度。
45.在一个具体示例中，两条热带气旋路径的空间相近程度的确定可通过以下具体步骤实现：
46.子步骤1：椭圆参数设置。
47.设样本热带气旋路径上依次选取相邻的两个路径点，如p1和p2，为其椭圆e1的焦点，焦距为e1，路径点p1和p2的尺度参数为scale
p1
和scale
p2
，其10级风圈半径分别为r10
p1
和r10
p2
；
48.设样本热带气旋路径上依次选取相邻的两个路径点，如q1和q2，为其椭圆s1的焦点，焦距为s1，路径点q1和q2的尺度参数为scale
q1
和scale
q2
，其10级风圈半径分别r10
q1
和r10
q2
；
49.通过公式确定椭圆长半轴a，其中dist为两个焦点间的球面距离，计算公式如下：
[0050][0051]
其中r为地球平均半径；
[0052]
k为尺度系数，根据以下情况分别计算：
[0053][0054]
情况1：路径点尺度参数存在且平均尺度参数＞dist值；
[0055]
情况2：路径点尺度参数不存在且10级风圈半径平均值＞dist值；
[0056]
情况3：除去情况1和情况2的其它情况；
[0057]
椭圆短半轴b通过公式确定，其中，a为椭圆长半轴，e为椭圆焦距。
[0058]
子步骤2：分别计算样本热带气旋路径中的路径点与历史热带气旋路径的椭圆轨迹以及历史热带气旋路径中的路径点与样本热带气旋路径的椭圆轨迹的匹配程度。
[0059]
以样本热带气旋路径中的路径点与历史热带气旋路径的椭圆轨迹为例，路径点与椭圆轨迹的匹配程度用m(pi,sq)表达，计算方法如下：
[0060][0061]
情况1：当pi与sq中任意一个椭圆处于包含关系时；
[0062]
情况2：除去情况1的其它情况。
[0063]
同样的方法计算历史热带气旋路径中的路径点与样本热带气旋路径的椭圆轨迹的匹配程度。
[0064]
子步骤3：计算样本热带气旋路径与历史热带气旋路径的空间相近程度。
[0065]
累加子步骤2求得的样本热带气旋路径中的路径点与历史热带气旋路径的椭圆轨迹的匹配程度除以样本热带气旋路径中的路径点数量；累加子步骤2求得的历史热带气旋路径中的路径点与样本热带气旋路径的椭圆轨迹的匹配程度除以历史热带气旋路径中的路径点数量，计算公式如下：
[0066][0067]
m(p,q)取值范围在0-1，越接近1表明两条路径空间相近程度更高。
[0068]
存在多条历史热带气旋路径的情况下，则循环计算样本热带气旋路径与每条历史热带气旋路径的空间相近程度，并存储计算结果。
[0069]
重叠程度确定步骤103，确定样本热带气旋路径与历史热带气旋路径的重叠程度。
[0070]
该步骤在具体执行中，可通过构造热带气旋椭圆轨迹，分别计算历史热带气旋路径中的路径点与样本热带气旋路径所对应的椭圆轨迹之间的归一化距离和样本热带气旋路径中的路径点与历史热带气旋路径所对应的椭圆轨迹之间的归一化距离，获得历史热带气旋路径与样本热带气旋路径的重叠程度。每条历史热带气旋路径均可采用相同的方法处理。举例说明如下：
[0071]
设样本热带气旋路径为cyclone
p
＝《p1,p2,...pn》，pi为样本热带气旋路径中的路径点(已按时间序列升序排列)；
[0072]
设历史热带气旋路径为cycloneq＝《q1,q2,...qn》，qi为历史热带气旋路径中的路径点(已按时间序列升序排列)；
[0073]
以样本热带气旋路径中的前后两个路径点为焦点，生成一个椭圆，依此类推，生成一系列的椭圆，形成椭圆序列e
p
＝《e1,e2,...en》，e为椭圆序列的椭圆；同样方式生成历史热带气旋路径的一系列椭圆，形成椭圆序列sq＝《s1,s2,...sn》，s为椭圆序列的椭圆；
[0074]
椭圆长半轴、短半轴以及焦距的参数分别为a、b、e。
[0075]
具体计算步骤如下：
[0076]
子步骤1：计算路径点与椭圆轨迹的归一化距离。
[0077]
计算公式如下：
[0078]
设p为样本热带气旋路径中的路径点，p∈cyclone
p
，s为历史热带气旋路径的椭圆序列中的椭圆，s∈sq；
[0079][0080]
其中，k可取值为椭圆s的焦点s
f1
,s
f2
或者椭圆原点s
fc
。dist
euc
(p,k)为p点到k点的欧式距离；
[0081]
α取值公式为其中b为椭圆短半轴、a为椭圆长半轴；
[0082]
β取值公式为其中e为所对应的椭圆焦距。
[0083]
子步骤2：计算样本热带气旋路径与历史热带气旋路径的重叠程度。
[0084]
计算公式如下：
[0085][0086]
ovp(p,q)取值范围在0-1，越接近1表明两条路径重叠程度更高。
[0087]
存在多条历史热带气旋路径的情况下，则循环计算样本热带气旋路径与每条历史热带气旋路径的重叠程度，并存储计算结果。
[0088]
连续程度确定步骤104，确定样本热带气旋路径与历史热带气旋路径的连续程度。
[0089]
该步骤在具体执行中，可通过构造热带气旋椭圆轨迹，分别计算历史热带气旋路径中的路径点与样本热带气旋路径所对应的椭圆轨迹的第一匹配位置的序列和样本热带气旋路径中的路径点与历史热带气旋路径所对应的椭圆轨迹第一匹配位置的序列，确定获得历史热带气旋路径与样本热带气旋路径的连续程度。每条历史热带气旋路径均可采用相同的方法处理。举例说明如下：
[0090]
设样本热带气旋路径为cyclone
p
＝《p1,p2,...pn》，pi为样本热带气旋路径中的路
径点(已按时间序列升序排列)；
[0091]
设历史热带气旋路径为cycloneq＝《q1,q2,...qn》，qi为历史热带气旋路径中的路径点(已按时间序列升序排列)；
[0092]
以样本热带气旋路径中的前后两个路径点为焦点，生成一个椭圆，依此类推，生成一系列的椭圆，形成椭圆序列e
p
＝《e1,e2,...en》，e为椭圆序列的椭圆；同样方式生成历史热带气旋路径的一系列椭圆，形成椭圆序列sq＝《s1,s2,...sn》，s为椭圆序列的椭圆；
[0093]
具体计算步骤如下：
[0094]
子步骤1：计算路径点与椭圆轨迹的第一匹配位置。
[0095]
将热带气旋路径中的路径点与椭圆轨迹匹配的顺序序列定义为seque(pi,sq)，如seque(p1,sq)＝《0,2,3|。将第一匹配位置定义为f(pi,sq)，其按以下情况取值：
[0096]
1)若i＝1，且seque(pi,sq)不为空的情况下，f(p1,sq)＝minseque(p1,sq)；
[0097]
2)若i＞1，且seque(pi,sq)不为空的情况下，f(pi,sq)＝min(seque(pi,sq),f(p
i-1
,sq))；
[0098]
3)除去情况1)、2)的其它情况时，f(pi,sq)＝-1。
[0099]
子步骤2：计算样本热带气旋路径与历史热带气旋路径的连续程度。
[0100]
样本热带气旋路径可以表示为其每个路径点的第一匹配位置序列，表达为seque
p
＝《f(p1,sq),f(p2,sq),...f(pn,sq)》；同样的，历史热带气旋路径可以表示为其每个路径点的第一匹配位置序列，表达为sequeq＝《f(q1,e
p
),f(q2,e
p
),...f(qn,e
p
)》；
[0101]
二者的连续程度用如下公式表达：
[0102][0103]
compute(parameter)函数计算规则如下：parameter是第一匹配位置序列；
[0104]
compute为分段函数，当遍历seque
p
时，设i为顺序数字自0开始，seque
p
[i]＝f(pi,sq)，
[0105]
1)当i＝0且seque
p
[i]的值≠-1时，compute(seque
p
[i])＝1；
[0106]
2)当i＞0且seque
p
[i]的值≠-1且seque
p
[i]≥seque
p
[i-1]时，compute(seque
p
[i])＝1；
[0107]
3)除去情况1)、2)的其它情况时，compute(seque
p
[i])＝0；
[0108]
con(p,q)取值范围在0-1，越接近1表明两条路径连续性更高。
[0109]
存在多条历史热带气旋路径的情况下，则循环计算样本热带气旋路径与每条历史热带气旋路径的连续程度，并存储计算结果。
[0110]
强度相似程度确定步骤105，确定样本热带气旋路径与历史热带气旋路径的强度相似程度。
[0111]
该步骤在具体执行中，可通过分别计算历史热带气旋路径中的路径点和样本热带气旋路径中的路径点的强度集合，分别构建历史热带气旋路径与样本热带气旋路径的强度矢量，通过计算两个强度矢量的余弦相似度值，确定历史热带气旋路径与样本热带气旋路径的强度相似程度。每条历史热带气旋路径均可采用相同的方法处理。举例说明如下：
[0112]
设样本热带气旋路径为cyclone
p
＝《p1,p2,...pn》，pi为样本热带气旋路径中的路
径点(已按时间序列升序排列)；
[0113]
设历史热带气旋路径为cycloneq＝《q1,q2,...qn》，qi为历史热带气旋路径中的路径点(已按时间序列升序排列)；
[0114]
具体计算步骤如下：
[0115]
热带气旋的强度分为弱于热带低压、热带低压、热带风暴、强热带风暴、台风、强台风、超强台风和变性8种；
[0116]
遍历样本热带气旋路径中的路径点，读取路径点上的热带气旋强度信息，记录每个强度出现的频次，生成强度矢量，并生成样本热带气旋路径的强度矢量，设为vector
p
＝《pd,p
td
,...p
de
》；
[0117]
遍历历史热带气旋路径中的路径点，读取路径点上的热带气旋强度信息，记录每个强度出现的频次，生成强度矢量，并生成历史热带气旋路径的强度矢量，设为vectorq＝《qd,q
td
,...q
de
》；
[0118]
通过计算两个矢量的余弦值来确定强度相似程度，公式如下：
[0119][0120]
cos{p,q)取值范围在0-1，越接近1表明历史热带气旋路径与样本热带气旋路径的强度更为接近。
[0121]
存在多条历史热带气旋路径的情况下，则循环计算样本热带气旋路径与每条历史热带气旋路径的强度相似程度，并存储计算结果。
[0122]
路径相似程度确定步骤106，综合样本热带气旋路径与历史热带气旋路径的空间相近程度、重叠程度、连续程度和强度相似程度，确定二者的路径相似程度。
[0123]
该步骤在具体执行中，综合空间相近程度确定步骤到强度相似程度确定步骤所确定的结果，确定样本热带气旋路径与历史热带气旋路径的路径相似程度。举例说明如下：
[0124]
按照热带气旋编号，提取空间相近程度确定步骤到强度相似程度确定步骤的计算结果，按如下公式计算：total(p,q)＝wmm(p,q)+wo
vp
ovp(p,q)+wconcon(p,q)+wcoscos(p,q)
[0125]
wm、w
ovp
、w
con
和w
cos
分别为空间相近程度、重叠程度、连续程度和强度相似程度的权重系数，wm+w
ovp
+w
con
+w
cos
的权重和为1，wm、w
ovp
、w
con
和w
cos
在实际应用中可例如取值为0.45、0.2、0.25和0.10，也就是说，可以空间相近程度的权重为最大，重叠程度和连续程度的权重小于空间相近程度的权重，以强度相似程度的权重最小。
[0126]
实施例2
[0127]
与实施例1相比，实施例2的不同在于在历史数据获取步骤之后增加外接矩形相近程度确定步骤，该步骤用于确定样本热带气旋路径与历史热带气旋路径的外接矩形空间重叠程度，并根据该外接矩形空间重叠程度的大小选择后续步骤的计算方式，具体如下：
[0128]
设样本热带气旋路径为cyclone
p
＝《p1,p2,...pn》，pi为样本热带气旋路径中的路径点(已按时间序列升序排列)；
[0129]
设历史热带气旋路径为cycloneq＝《q1,q2,...qn》，qi为历史热带气旋路径中的路径点(已按时间序列升序排列)；
[0130]
热带气旋路径中的路径点由经纬度对组成，如p1表示为(lat,lon)，因此一条热带
气旋的外接矩形可以表示为[lat
min
,lat
max
,lon
min
,lon
max
]。通过样本热带气旋路径和历史热带气旋路径的外接矩形，判别两个矩形在地理空间的重叠程度，并根据所述外接矩形空间重叠程度的大小选择如下处理方式：
[0131]
当外接矩形的重叠程度大于或等于第一设定值时，循环执行空间相近程度确定步骤、重叠程度确定步骤和连续程度确定步骤三次，对于路径点数据的选取，第一次执行每m个点取一个路径点，第二次执行每n个点取一个路径点，第三次执行取全部路径点，综合三次执行所确定的空间相近程度、重叠程度和连续程度，作为最终的空间相近程度、重叠程度和连续程度，其中m大于n，n大于1；
[0132]
当外接矩形的重叠程度小于第一设定值且大于或等于低于第一设定值的第二设定值时，循环执行两次空间相近程度确定步骤、重叠程度确定步骤和连续程度确定步骤，对于路径点数据的选取，第一次执行每a个点取一个路径点，第二次执行取全部路径点，综合两次执行所确定的空间相近程度、重叠程度和连续程度求取其均值，作为最终的相近程度、重叠程度和连续程度，其中a大于1；
[0133]
当外接矩形的重叠程度大于0且小于第二设定值时，取全部路径点通过空间相近程度确定步骤、重叠程度确定步骤和连续程度确定步骤确定空间相近程度、重叠程度和连续程度；
[0134]
当外接矩形不重叠的情况下，判别为历史热带气旋路径和样本热带气旋路径不相似。
[0135]
举例说明如下：
[0136]
(1)当外接矩形的重叠程度大于或等于70％的情况下，做如下处理：
[0137]
通过空间相近程度确定步骤、重叠程度确定步骤和连续程度确定步骤开展三次计算，分别求得三次的空间相近程度、重叠程度和连续程度，其中，第一次计算的数据源是样本热带气旋路径为cyclone
p
＝《p1,p2,...pn》的每4个点取一个点组成的序列，即cyclone
p3
＝《p1,p5,p9...》，对于历史热带气旋路径采用相同取样方法，第二次计算的数据源是样本热带气旋路径为cyclone
p
＝《p1,p2,...pn》的每2个点取一个点组成的序列，即cyclone
p3
＝《p1,p3,p5...》，对于历史热带气旋路径采用相同取样方法，第三次计算采用原始序列，三次计算分别求得的空间相近程度、重叠程度和连续程度求取其均值，作为最终的空间相近程度、重叠程度和连续程度的数值。
[0138]
(2)当外接矩形的重叠程度在大于或等于34％且小于70％的情况下，做如下处理：
[0139]
通过空间相近程度确定步骤、重叠程度确定步骤和连续程度确定步骤开展两次计算，分别求得两次的空间相近程度、重叠程度和连续程度，其中，第一次计算的数据源是样本热带气旋路径为cyclone
p
＝《p1,p2,...pn》的每2个点取一个点组成的序列，即cyclone
p3
＝《p1,p3,p5...》，对于历史热带气旋路径采用相同取样方法，第二次计算采用原始序列，二次计算分别求得的空间相近程度、重叠程度和连续程度求取其均值，作为最终的相近程度、重叠程度和连续程度的数值。
[0140]
(3)当外接矩形的重叠程度在大于0且小于34％的情况下，做如下处理：
[0141]
开展1次计算，采用原始序列通过空间相近程度确定步骤、重叠程度确定步骤和连续程度确定步骤求得1次的空间相近程度、重叠程度和连续程度。
[0142]
(4)当外接矩形不重叠的情况下，做如下处理：
[0143]
直接判别为不相似。
[0144]
通过增加该步骤，能够更好地计算热带气旋路径的相似性，避免出现在局部相似度高，但整体相似度偏低的情况。
[0145]
本发明的上述方法可以通过在移动终端、计算机终端、服务器或者类似的计算设备中执行，因而对应上述方法，本发明还提供了存储介质、热带气旋路径相似程度确定装置的三个实施例。其中：
[0146]
实施例3
[0147]
本实施例提供了一种存储介质，所述存储介质包括存储的程序，其中，在所述程序运行时由处理器执行实施例1或2提供的方法。
[0148]
实施例4
[0149]
本实施例提供了一种热带气旋路径相似程度确定装置，该装置与根据实施例1和实施例2的方法相对应。参考图2所示，该装置包括：
[0150]
历史数据获取模块201，用于获取包括路径点数据集的热带气旋历史数据，所述路径点数据集包括热带气旋序号、时间、经度、纬度、强度、中心最低气压和热带气旋尺度，所述热带气旋历史数据包括一条或多条历史热带气旋的数据；
[0151]
空间相近程度确定模块202，用于确定样本热带气旋路径与历史热带气旋路径的空间相近程度；
[0152]
重叠程度确定模块203，用于确定样本热带气旋路径与历史热带气旋路径的重叠程度；
[0153]
连续程度确定模块204，用于确定样本热带气旋路径与历史热带气旋路径的连续程度；
[0154]
强度相似程度确定模块205，用于确定样本热带气旋路径与历史热带气旋路径的强度相似程度；
[0155]
路径相似程度确定模块206，用于综合样本热带气旋路径与历史热带气旋路径的空间相近程度、重叠程度、连续程度和强度相似程度，确定样本热带气旋和历史热带气旋的路径相似程度。
[0156]
在此基础上，还可增加外接矩形相近程度确定模块，用于确定样本热带气旋路径与历史热带气旋路径的外接矩形空间重叠程度，并根据所述外接矩形空间重叠程度的大小选择如下处理方式：
[0157]
当外接矩形的重叠程度超过第一设定值时，循环执行三次空间相近程度确定、重叠程度确定和连续程度确定，对于路径点数据的选取，第一次执行每m个点取一个路径点，第二次执行每n个点取一个路径点，第三次执行取全部路径点，综合三次执行所确定的空间相近程度、重叠程度和连续程度，作为最终的空间相近程度、重叠程度和连续程度，其中m大于n，n大于1；
[0158]
当外接矩形的重叠程度处于第一设定值和低于第一设定值的第二设定值之间时，循环执行两次空间相近程度确定、重叠程度确定和连续程度确定，对于路径点数据的选取，第一次执行每a个点取一个路径点，第二次执行取全部路径点，综合两次执行所确定的空间相近程度、重叠程度和连续程度求取其均值，作为最终的相近程度、重叠程度和连续程度，其中a大于1；
[0159]
当外接矩形的重叠程度处于0和第二设定值之间时，取全部路径点通过空间相近程度确定步骤、重叠程度确定步骤和连续程度确定步骤确定空间相近程度、重叠程度和连续程度；
[0160]
当外接矩形不重叠的情况下，判别为历史热带气旋路径和样本热带气旋路径不相似。
[0161]
实施例5
[0162]
本实施例提供了一种带气旋路径相似程度确定装置，该装置与根据实施例1和实施例2的方法相对应。参考图3所示，该装置包括：
[0163]
处理器301；
[0164]
存储器302，与所述处理器301连接，存储器302存储有程序，其中，在所述程序运行时由处理器301执行实施例1和实施例2提供的方法。
[0165]
尽管本文对于实施例3-5的描述未做全面展开，但本领域普通技术人员可以明了其具体实现与实施例1和实施例2中的方法是相对应，且可以获得与之一致的技术优势和效果。
[0166]
本发明上述实施例，结合相近、重叠、连续和强度4个方面，能够更为准确地确定历史热带气旋数据库中历史热带气旋个例与样本热带气旋的路径相似程度，以此为基础可进一步给出与样本热带气旋的相似排序结果，筛选排序出最为相似的热带气旋，进而再筛选背景环境场，选择与样本热带气旋相近的热带气旋，查询这些历史热带气旋的发展变化的情况以及沿岸验潮站的潮位及增水情况，结合数值预报模式输出结果、人工经验分析判断，给出对样本热带气旋可能造成风暴增水的预测结果。
[0167]
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。