[发明专利]平面散乱点集的快速网格化方法、网格化系统及存储介质

权利要求书

1.一种平面散乱点集的快速网格化方法，其特征在于，所述平面散乱点集的快速网格化方法包括：

获取平面散乱点集，确定散乱点集分布的最大矩形区域；采用分治法将平面散乱点集生成非结构化三角网并对所有三角形按顺序进行编号；通过设置网格间距，获得网格化的待插结点；逐一选择待插网格结点，快速搜索其所在的三角形并进行插值，即可。

2.如权利要求1所述平面散乱点集的快速网格化方法，其特征在于，所述平面散乱点集的快速网格化方法包括以下步骤：

步骤一，给定平面散乱点集(x,y,z)，按读取顺序编号[1,2,...,m]，并确定散乱点集分布的最大矩形区域:[X_min,X_max],[Y_min,Y_max]；采用分治法将平面散乱点集生成非结构化三角网；

步骤二，对所有三角形按顺序编号[1,2,...,n]，并记录每个三角形对应的3个顶点编号以及与其三边相邻的三角形编号；计算并一次性保存所有三角形的线性插值权系数；设置x,y方向的网格间距，获得网格化的待插结点；

步骤三，选择一个待插网格结点，并将三角网中的某个三角形作为种子三角形，利用待插点与种子三角形两顶点的向量积判别待插点与所述三角形边的相对位置关系；

步骤四，反复判别待插点与三角形边的关系，快速逼近并定位待插结点所在的三角形；若在三角网中能找到待插结点所在的三角形，则利用所述三角形的线性插值权系数对所述结点进行插值；若未找到，说明所述待插点不在三角网内部，则不作插值；

步骤五，重复步骤四，直至完成所有网格结点的插值，并将网格化结果按*.grd格式输出到文件；

步骤六，重复执行步骤三至步骤五，得到不同网格尺度的网格化结果。

3.如权利要求2所述平面散乱点集的快速网格化方法，其特征在于，步骤二中，所述对所有三角形按顺序编号[1,2,...,n]，并记录每个三角形对应的3个顶点编号以及与其三边相邻的三角形编号包括：

以平面散乱点集的坐标编号作为三角网中三角形的顶点编号，对于任一三角形三边相邻的三角形编号，则以三角网中三角形的顺序号作为编号，若三角形位于三角网的边部，部分边无相邻三角形时，则编号置为-1，所述编号的正、负符号表示待插网格结点位于三角网内部还是外部。

4.如权利要求2所述平面散乱点集的快速网格化方法，其特征在于，步骤二中，所述记录每个三角形对应的3个顶点编号以及与其三边相邻的三角形编号包括：

若三角形位于三角网内部则有3个相邻的三角形；若三角形位于三角网边部则有1或2个相邻的三角形。

5.如权利要求2所述平面散乱点集的快速网格化方法，其特征在于，步骤二中，所述计算并一次性保存所有三角形的线性插值权系数包括:

利用下式计算并一次性保存所有三角形的插值权系数a,b,c：

若y₁-y₃≠0，则

若y₁-y₃＝0，则

其中，A(x₁,y₁,z₁)，B(x₂,y₂,z₂)，C(x₃,y₃,z₃)表示三角形的三个顶点坐标；ax+by+c＝z表示平面方程；

对于三角形ABC内一点(x_p,y_p)，z_p＝ax_p+by_p+c_p。