[R-sig-Geo] Variogram model 3D

[Edzer Pebesma]

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