克里金插值c程式 克里金插值方法的原理

2021-10-12 08:57:57 字數 826 閱讀 8290

克里金插值方法原理(步驟1-4用來說明半方差函式的構建,步驟5說明了**模型,即如何求取未知點的數值):

半方差函式訓練樣本的獲取公式:

2. 構建散點圖(x軸:距離,y軸:半方差值)

3. 根據已有的函式擬合經驗半方差圖

arcgis中提供了五種函式:circular 、spherical 、exponential 、gaussian 、linear

4. 根據擬合後得到的半方差函式,可知塊金(nugget)、基台(sill)和變程(range)。

塊金是距離為0時的半方差函式值。根據半方差計算公式可知,當距離為0時,半方差應該為0,因此塊金用來反映重複測量的偶然誤差和測量裝置的系統誤差。

在變程範圍內,資料的空間自相關明顯,當超過變程時,可認為資料之間不再存在相關性。

5. **未知點的數值

與idw插值方法相似,均是對未知點的特定鄰域範圍的測量點或者特定數量的相鄰測量點的數值進行加權相加,以求得未知點的數值,實現對未知點的**。其中,周圍測量點的權重根據半方差函式來確定。

均來自arcgis幫助文件。

克里金插值

由於用supermap objects 沒有解決插值範圍的問題 見本版帖子 求助!哪位大俠在用supermap objects,請教乙個插值區域的問題 改用arcgis engine來做,現在遇到同樣的問題。使用iinterpolationop的krige方法已經實現了插值,但範圍侷限於氣象站點的外...

普通克里金插值

最近因為專案需要,研究了下克里金插值演算法。在地質學中,克里金插值演算法是一種使用的空間屬性估計技術,克里金插值說到底是個回歸問題,且依據的因素只有兩個位置之間的距離。克里金插值演算法又分為很多中,比如普通克里金插值,簡單克里金插值等,不同的克里金插值演算法只是假設條件不同。下面以普通克里金為例來說...

克里金插值java實現

主要對立面的配置檔案進行簡單的說明 其中乙個 par檔案 內容如下 parameters for jk3d start of parameters e test.dat file with data 1 2 3 4 0 columns for x,y,z,var,sec var 1.0e21 1.0...