克里金插值java實現

2021-06-28 12:05:45 字數 2194 閱讀 4419

主要對立面的配置檔案進行簡單的說明

其中乙個***.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.0e21                 -trimming limits

0                                -option: 0=grid, 1=cross, 2=jackknife  **模式

xvk.dat                          -file with jackknife data

1   2   0    3    0              -columns for x,y,z,vr and sec var

3                                -debugging level: 0,1,2,3

testdata-iw3d-3d.dbg             -file for debugging output

real.out             -file for kriged output    計算後的檔案輸出路徑

10 0 1                         -nx,xmn,xsiz  次數,起始值,步長

10 0 1                         -ny,ymn,ysiz

1 0 1                         -nz,zmn,zsiz

1    1      1                    -x,y and z block discretization

1    1600                          -min, max data for kriging   最小、最大值

16    8                           -max per octant (0-> not used), blank if more than this many octants are empty

1  1  1.1                    -maximum search radii  搜尋半徑

0.0   0.0   0.0                  -angles for search ellipsoid

1     -1001                      -0=sk,1=ok,2=non-st sk,3=exdrift

0 0 0 0 0 0 0 0 0                -drift: x,y,z,xx,yy,zz,xy,xz,zy

0                                -0, variable; 1, estimate trend

extdrift.dat                     -gridded file with drift/mean

4                                -column number in gridded file

1    0.0                         -nst, nugget effect

3   0.04    0.0   0.0   0.0      \it,cc,ang1,ang2,ang3

0.994  0.994   0.84            \a_hmax, a_hmin, a_vert

3   0.0336  0.0   0.0  90.0      \it,cc,ang1,ang2,ang3

0.85 0.85 0.0                 \a_hmax, a_hmin, a_vert

***.dat  內容如下:

primary data4x

yzval1 1 0 200      --x  y  z 座標   值

2 2 0 300

3 3 0 1000

4 4 0 500

5 5 0 100

6 6 0 100

7 7 0 100

8 8 0 100

9 9 0 100

克里金插值

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

普通克里金插值

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

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

克里金插值方法原理 步驟1 4用來說明半方差函式的構建,步驟5說明了 模型,即如何求取未知點的數值 半方差函式訓練樣本的獲取公式 2.構建散點圖 x軸 距離,y軸 半方差值 3.根據已有的函式擬合經驗半方差圖 arcgis中提供了五種函式 circular spherical exponential...