Now I understood. Thank you!



________________________________
De: Edzer Pebesma <edzer.pebesma@uni-muenster.de>
Para: r-sig-geo@r-project.org 
Enviadas: Terça-feira, 24 de Julho de 2012 14:01
Assunto: Re: [R-sig-Geo] 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@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@wwu.de

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
	[[alternative HTML version deleted]]

