[发明专利]一种针对相交线状地理要素的空间剖分方法

权利要求书

1.一种针对相交线状地理要素的空间剖分方法，其特征在于，具体步骤如下：

步骤一：加载GIS线要素集合L＝{l_i|i＝1，2，…，n}，n为线要素的个数；

步骤二：求取线要素集合L中不同线要素l_i和l_j间的交叉点C_ij(x，y)，并进行内插处理；

步骤三：将线要素集合L离散为点集P＝{p_i(x_i，y_i)|i＝1，2，…，k}，k为点要素的个数，其中包含交叉辅助点的个数，并使用点集P的“线要素”属性来记录该点所归属的线要素；

步骤四：基于点集P，生成Delaunay三角网D；

步骤五：基于Delaunay三角网D，生成Voronoi图V；

步骤六：将点集P中各点的“线要素”属性值，赋给包含该点的相应Thicssen多边形t_i的“原线要素”的属性值；

步骤七：对Voronoi图V进行裁剪处理，通过求解给定矩形范围R与图V中各Thicssen多边形t_i的交集，生成裁剪后的Thicssen多边形集合T_i'；

步骤八：根据裁剪后的Thicssen多边形集合T_i'的“原线要素”属性，依次将Voronoi图V中具有相同“原线要素”属性值的Thicssen多边形t_i'合并为一个多边形，每个多边形区域便是一个线要素的影响范围，从而实现相交地理要素的空间剖分。

2.根据权利要求1所述的一种针对相交线状地理要素的空间剖分方法，其特征在于，所述步骤二的具体过程为：

(a)首先，针对线要素l_i，读取线要素l_i上的点集，查找线要素l_i上交叉点邻近的两个点，分别记为C_a(x_a，y_a)，C_b(x_b，y_b)；

(b)根据下式(1)，计算待插入点I_c(x_c，y_c)、I_d(x_d，y_d)坐标，其中(x，y)为交叉点C_ij的坐标，(x_a，y_a)和(x_b，y_b)为交叉点附近的两点C_a，C_b的坐标，(x_c，y_c)和(x_d，y_d)为待插入点I_c，I_d的坐标，s为步长因子，s可取10，100，1000，...10ⁿ，

$\left\{\begin{array}{c}\mathrm{\Δ}d=\frac{\sqrt{{(y_b-y_a)}^2+{(x_b-x_a)}^2}}{s}\\ k=\frac{y_b-y_a}{x_b-x_a}\\ x_c=x+\mathrm{\Δ}d\\ y_c=k\×x_c+y_a-k\×x_a\\ x_d=x-\mathrm{\Δ}d\\ y_d=k\×x_d+y_a-k\×x_a\end{array}\right.---\left(1\right)$

(c)将点I_c和I_d插入到线要素l_i中的点C_a、C_b之间；

(d)其次，针对该交叉点所在的另一线要素l_j，重复执行步骤(a)-(c)，将生成的相应插入点插入到线要素l_j中。