[发明专利]一种基于泰森多边形动态调整的面雨量计算方法

权利要求书

1.一种基于泰森多边形动态调整的面雨量计算方法，其特征在于，包括以下步骤:

S1，获取某一时段各雨量站降水观测值P＝{P₁，P₂，P₃……P_n},其中n为有数据的雨量站个数；

S2，获取有数据的雨量站序列Sta＝{Sta₁，Sta₂…….Sta_n},其中n为有数据的雨量站个数；

S3，判断雨量站序列Sta是否存在于历史数据库中的雨量站序列集合Stas；若存在，则直接获取历史数据中该Sta对应的泰森多边形划分方案Vor，并跳转至步骤S6；否则跳转步骤S4；

S4，重新基于待测的雨量站序列Sta进行泰森多边形的划分，形成泰森多边形划分方案Vor’；

S5，将新设计的泰森多边形划分方案添加至历史数据库中的泰森多边形划分方案集Vors；

S6，利用Vor’或者Vor计算待测区域的面雨量。

2.根据权利要求1所述的机遇泰森多边形动态调整的面雨量计算方法，其特征在于，步骤S2和步骤S3之间还包括：

获取历史上出现过的泰森多边形划分方案Vors所涉及的面雨量序列集合Stas，这里的Vors是用于计算当前区域面雨量的，历史上出现过的划分方案的集合，Stas是指Vors中的每一种划分方案所分别包含的雨量站序列的集合。

3.根据权利要求1所述的基于泰森多边形动态调整的面雨量计算方法，其特征在于，步骤S4中具体包括：

S41，针对雨量站的离散点构建Delaunay三角网，对离散点和形成的三角形编号，记录每个三角形是由哪三个离散点构成的；

S42，找出与每个离散点相邻的所有三角形的编号，并记录下来；

S43，对与每个离散点相邻的三角形按顺时针或逆时针方向排序，计算每个三角形的外接圆圆心，并记录之；

S44，根据每个离散点的相邻三角形，连接这些相邻三角形的外接圆圆心，即得到泰森多边形。

4.根据权利要求3所述的基于泰森多边形动态调整的面雨量计算方法，其特征在于，步骤S43中的对离散点相邻的三角形按顺时针或逆时针方向排序的具体步骤是：

设某个离散点为o，找出以o为顶点的一个三角形，设为A；取三角形A除o以外的另一顶点，设为a，则另一个顶点也可找出，即为f；则下一个三角形必然是以of为边的，即为三角形F；三角形F的另一顶点为e，则下一三角形是以oe为边的；如此重复进行，直到回到oa边。

5.根据权利要求1所述的基于泰森多边形动态调整的面雨量计算方法，其特征在于，步骤S6中计算待测区域的面雨量的公式具体为：

其中，S_i为第i个泰森多边形的面积，P_i为第i个雨量站某一时段监测的雨量，P为该时段流域上的面雨量，n为流域上雨量站的数量，i∈[1,6]。