[發(fā)明專利]一種基于克里金插值方程的水土流失空間監(jiān)測(cè)方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201410007114.9
• 申請(qǐng)日: 2014-01-08
• 公開（公告）號(hào): CN103699809B
• 公開（公告）日: 2017-02-08
• 發(fā)明（設(shè)計(jì)）人: 殷水清;劉寶元;劉峰;章文波;劉素紅;趙瑩;王友勝
• 申請(qǐng)（專利權(quán)）人: 北京師范大學(xué)
• 主分類號(hào): G06F19/00
• 分類號(hào): G06F19/00
• 代理公司: 北京孚睿灣知識(shí)產(chǎn)權(quán)代理事務(wù)所(普通合伙)11474
• 代理人: 舒麗亞,劉巖

1.一種基于克里金插值方程的水土流失空間監(jiān)測(cè)方法，其特征在于，其實(shí)施步驟如下：

1)獲取待監(jiān)測(cè)區(qū)域的各個(gè)調(diào)查單元的中心經(jīng)緯度，將各個(gè)所述調(diào)查單元分別記為各個(gè)相應(yīng)的點(diǎn)x_i，各個(gè)所述調(diào)查單元的中心經(jīng)緯度作為各個(gè)相應(yīng)所述點(diǎn)x_i的地理坐標(biāo)；

2)獲取待監(jiān)測(cè)區(qū)域的各個(gè)所述相應(yīng)點(diǎn)x_i的所有水土流失指標(biāo)，對(duì)各個(gè)所述調(diào)查單元的相應(yīng)點(diǎn)x_i的所有水土流失指標(biāo)根據(jù)克里金插值方程分別進(jìn)行插值，得到被監(jiān)測(cè)區(qū)域內(nèi)待插值點(diǎn)x₀的屬性估值，其中，i＝1，2，…，n，x₀為待插值點(diǎn)，Z^*(x₀)為待插值點(diǎn)x₀的屬性估值，Z(x_i)為點(diǎn)x_i的屬性值，λ_i為i的待定權(quán)重系數(shù)；

3)分別對(duì)所有待監(jiān)測(cè)區(qū)域的各個(gè)點(diǎn)x_i的所有水土流失指標(biāo)進(jìn)行空間插值，得到所述各待監(jiān)測(cè)區(qū)域的水土流失強(qiáng)度比例和面積比例；

4)根據(jù)各待監(jiān)測(cè)區(qū)域水土流失的強(qiáng)度和面積，繪制水土流失各指標(biāo)柵格圖，并根據(jù)所述水土流失柵格圖繪制水力侵蝕強(qiáng)度圖。

2.根據(jù)權(quán)利要求1所述的基于克里金插值方程的水土流失空間監(jiān)測(cè)方法，其特征在于，所述步驟2)中i的待定權(quán)重系數(shù)具體是指根據(jù)和計(jì)算得到i的待定權(quán)重系數(shù)λ_i，其中，γ(x_i，x_j)為被監(jiān)測(cè)區(qū)域內(nèi)的兩點(diǎn)x_i，x_j的變異函數(shù)，μ為拉格朗日乘子。

3.根據(jù)權(quán)利要求2所述的基于克里金插值方程的水土流失空間監(jiān)測(cè)方法，其特征在于，所述被監(jiān)測(cè)區(qū)域內(nèi)的兩點(diǎn)x_i，x_j的變異函數(shù)具體是根據(jù)計(jì)算得到x_i，x_j的變異函數(shù)，其中，Var表示方差，i，j＝1，2，…，n。

4.根據(jù)權(quán)利要求1所述的基于克里金插值方程的水土流失空間監(jiān)測(cè)方法，其特征在于，所述所有水土流失指標(biāo)包括4個(gè)類別67個(gè)指標(biāo)，第一類別包括7個(gè)指標(biāo)，分別為水蝕比例、侵蝕比例以及水蝕模數(shù)；第二類別包括20個(gè)指標(biāo)，分別為不同土地利用類型面積比例；第三類別包括20個(gè)指標(biāo)，分別為不同土地利用類型水蝕比例；第四類別包括20個(gè)指標(biāo)，分別為不同土地利用類型平均模數(shù)。

5.根據(jù)權(quán)利要求4所述的基于克里金插值方程的水土流失空間監(jiān)測(cè)方法，其特征在于，所述水蝕比例和所述侵蝕比例的插值以各調(diào)查單元對(duì)應(yīng)的各指標(biāo)分別進(jìn)行插值，所述水蝕模數(shù)的插值以各調(diào)查單元流失量為指標(biāo)進(jìn)行插值，所述不同土地利用類型面積比例進(jìn)行插值時(shí)，遇到某個(gè)調(diào)查單元沒有某種土地利用類型時(shí)，所述面積比例賦值為0，所述不同土地利用類型水蝕比例進(jìn)行插值時(shí)，遇到某個(gè)調(diào)查單元沒有某種土地利用類型時(shí)，所述水蝕比例賦值為0，所述不同土地利用類型平均模數(shù)進(jìn)行插值時(shí)，先對(duì)模比進(jìn)行插值，遇到某個(gè)調(diào)查單元沒有某種土地利用類型時(shí)，所述模比賦值為0。

6.根據(jù)權(quán)利要求5所述的基于克里金插值方程的水土流失空間監(jiān)測(cè)方法，其特征在于，所述模比為不同土地利用類型模數(shù)乘以不同土地利用類型面積比例。

7.根據(jù)權(quán)利要求1所述的基于克里金插值方程的水土流失空間監(jiān)測(cè)方法，其特征在于，還包括對(duì)水土流失面積比例進(jìn)行插值，所述水土流失面積比例包括輕度、中度、強(qiáng)烈、極強(qiáng)烈和劇烈5個(gè)等級(jí)，對(duì)所述水土流失面積的上述各個(gè)等級(jí)分別進(jìn)行插值并進(jìn)行相應(yīng)的調(diào)整：調(diào)整后的某侵蝕強(qiáng)度比例＝插值后某侵蝕強(qiáng)度比例*插值后水土流失面積比例/插值后5級(jí)侵蝕強(qiáng)度比例之和。