Variogram model 3D
Could you please tell which part of the help of vgm is unclear to you? As in GSLIB, X is the major direction, Y the first minor, and Z the second minor direction. Then, if anis = c(p,q,r,s,t) t is the anisotropy ration between the X and Z direction.
On 07/24/2012 03:49 AM, Tiago Vieira wrote:
hi,
Excuse my ignorance, but I've read the help of the function vgm in
package gstat and still do not understand how to change the vertical
range of the variogram model in 3D, as is possible in gslib (i.e. vmodel
and kt3d).
Thank you,
Tiago Vieira.
hypothetical examplein gslib:
Parameters for KT3D
*******************
START OF PARAMETERS:
data3D.dat \file with data
0 1 2 3 4 0 \ columns for DH, 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
kt3d.dbg \file for debugging output
teste3D.out \file for kriged output
51 0.0 2.0 \nx,xmn,xsiz
51 0.0 2.0 \ny,ymn,ysiz
51 0.0 2.0 \nz,zmn,zsiz
1 1 1 \x,y and z block discretization
4 8 \min, max data for kriging
0 \max per octant (0-> not used)
60.0 60.0 60.0 \maximum search radii
0.0 0.0 0.0 \angles for search ellipsoid
1 2.302 \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
1 3.0 0.0 0.0 0.0 \it,cc,ang1,ang2,ang3
40.0 40.0 0.0 \a_hmax, a_hmin, a_vert
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Edzer Pebesma Institute for Geoinformatics (ifgi), University of M?nster Weseler Stra?e 253, 48151 M?nster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de http://www.52north.org/geostatistics e.pebesma at wwu.de