geostatsp (version 0.4.4)

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 rotaiton, 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(0.5, 0.1))

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

myGRF=GaussRF(myraster, model=mymodel)

plot(myGRF)