[发明专利]一种基于线性内插的规则格网地形建模方法

权利要求书

1.一种基于线性内插的规则格网地形建模方法，包括以下步骤：

步骤(1)：设定阈值T；

步骤(2)：选择格网单元，计算高程值的方差Var；

步骤(3)：判断阈值T和方差Var的大小，若阈值大于方差值转步骤（4），否则转步骤（5）；

步骤(4)：采用双线性多项式进行地形的插值计算，建立当前格网的地形结构；

步骤(5)：采用双三次多项式进行地形的插值计算，建立当前格网的地形结构；

步骤(6)：计算所有格网单元，完成基于规则格网的地形建模。

2.根据权利要求1所述的一种基于线性内插的规则格网地形建模方法，其特征在于：所述步骤（1）中的阈值T，是根据地形平坦程度和建模精度要求进行设定。

3.根据权利要求2所述的一种基于线性内插的规则格网地形建模方法，其特征在于：所述步骤(2)中高程值的方差Var计算如下：

将区域平面划分为规则的格网单元，各单元顶点对应一个高程值；选择其中的一个单元，计算高程值的方差Var，如式（1）和式（2）所示：

（1）

（2）

式中：z _{i, j}、z_i+1,j、z_i,j+1、z_i+1,j+1为规则格网中第i行第j列格网单元顶点的高程值；为该格网单元高程的平均值。

4.根据权利要求3所述的一种基于线性内插的规则格网地形建模方法，其特征在于：所述的步骤（4）中采用双线性多项式进行地形的插值计算，利用待插值点所在格网的四个顶点坐标计算多项式系数，待插值点高程z的表达式为：

（3）

式中：x、y为待插值点的平面坐标，a₁₁、a₁₂、a₂₁、a₂₂为双线性多项式系数。

5.根据权利要求4所述的一种基于线性内插的规则格网地形建模方法，其特征在于：所述的步骤（5）中采用双三次多项式进行地形的插值计算，多项式系数由待插值点所在格网和邻近格网的16个顶点坐标决定，高程z表达式为：

（4）

式中：x、y为待插值点的平面坐标，为双三次多项式系数。