geostatsp (version 0.5.8)

anisoMatrix: Calculates a 2-D matrix for geometric anisotropy

Description

Calculates a 2-D matrix for geometric anisotropy, for use with GaussRF

Usage

anisoMatrix(angle, range=NULL,scale=1/range)

Arguments

angle

angle of rotation, in radians

range

range for x and y axes after rotation, ignored if scale is specified

scale

scale for x and y axes after rotation

Value

• A 2 by 2 matrix

Details

Produces matrix A such that (xnew, ynew)' = A (x,y)'

See Also

GaussRF

Examples

myAnisoMat = anisoMatrix(-pi/8, c(2, 0.5))

mymodel=list("$", var=2,  
		A=myAnisoMat,
		list("matern", nu=2))
		
myraster = raster(nrows=41,ncols=41,xmn=0,xmx=10,ymn=0,ymx=10)

myGRF=GaussRF(myraster, model=mymodel)

plot(myGRF)