[發明專利]一種反距離平方加權空間插值計算方法

在審

• 申請號: 202010123012.9
• 申請日: 2020-02-27
• 公開（公告）號: CN111339476A
• 公開（公告）日: 2020-06-26
• 發明（設計）人: 陳皓銳;米博宇;劉靜;陶園;王少麗;高黎輝;戴瑋
• 申請（專利權）人: 中國水利水電科學研究院
• 主分類號: G06F17/10
• 分類號: G06F17/10
• 代理公司: 成都正華專利代理事務所(普通合伙) 51229
• 代理人: 李蕊

1.一種反距離平方加權空間插值計算方法，定義x-y直角坐標系內的樣本點個數N≥3，第i個樣本點為Ai，Ai的樣本值為z_i，待插值點為P，i＝1,2,…N；

其特征在于，所述方法包括以下步驟：

S1、計算樣本點Ai與待插值點P之間的距離d_i；

S2、根據距離d_i，計算樣本點Ai相對于待插值點P的距離加權系數wd_i；

S3、判斷樣本點Ai與待插值點P形成的向量在x-y直角坐標系中的位置；

S4、根據向量在x-y直角坐標系中的位置，計算向量的方位角α_i；

S5、根據方位角α_i，計算樣本點Ai的方位占角θ_i；

S6、根據方位占角θ_i，計算樣本點Ai的方位加權系數wr_i；

S7、根據方位加權系數wr_i、距離加權系數wd_i計算樣本點Ai的綜合加權系數w_i；

S8、根據綜合加權系數w_i、樣本值z_i計算待插值點P的插值z。

2.根據權利要求1所述的反距離平方加權空間插值計算方法，其特征在于，設樣本點Ai的空間坐標為(x_i,y_i)，待插值點P的空間坐標為(x，y)，則所述步驟S1中，

3.根據權利要求2所述的反距離平方加權空間插值計算方法，其特征在于，所述步驟S2中，其中d_k為第k個樣本點與待插值點P之間的距離。

4.根據權利要求3所述的反距離平方加權空間插值計算方法，其特征在于，所述步驟S3中，向量在x-y直角坐標系中的位置判斷方法具體包括：

若x_i-x0，且y_i-y0，則位于第一象限；

若x_i-x0，且y_i-y0，則位于第二象限；

若x_i-x0，且y_i-y0，則位于第三象限；

若x_i-x0，且y_i-y0，則位于第四象限；

若(x_i-x)×(y_i-y)＝0，則位于坐標軸上。

5.根據權利要求4所述的反距離平方加權空間插值計算方法，其特征在于，所述步驟S4中，方位角α_i的計算方法包括：

若位于第一象限，則

若位于第二象限，則

若位于第三象限，則

若位于第四象限，則

若位于坐標軸上，且y_i-y＝0，x_ix，則

若位于坐標軸上，且y_i-y＝0，x_ix，則

若位于坐標軸上，且x_i-x＝0，y_iy，則α_i＝π；

若位于坐標軸上，且x_i-x＝0，y_iy，則α_i＝0。

6.根據權利要求5所述的反距離平方加權空間插值計算方法，其特征在于，所述步驟S5中，方位占角θ_i的計算方法包括：

當i＝1時，

當2≤i≤N-1時，

當i＝N時，

7.根據權利要求6所述的反距離平方加權空間插值計算方法，其特征在于，所述步驟S6中，方位加權系數

8.根據權利要求7所述的反距離平方加權空間插值計算方法，其特征在于，所述步驟S7中，綜合加權系數

9.根據權利要求8所述的反距離平方加權空間插值計算方法，其特征在于，所述步驟S8中，待插值點P的插值