主要對立面的配置檔案進行簡單的說明
其中乙個***.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...