source("gensurf.R")
#install.packages("fields")
library(fields)

locs=expand.grid(1:24,1:24)
cga=covpow(locs,pow=2,scale=4)
z=rmultnorm(1,rep(0,24*24),cga)
persp(matrix(z,ncol=24))

cga=covpow(locs,pow=2,scale=8)
z=rmultnorm(1,rep(0,24*24),cga)
persp(matrix(z,ncol=24))

cga=covpow(locs,pow=2,scale=2)
z=rmultnorm(1,rep(0,24*24),cga)
persp(matrix(z,ncol=24))


h=seq(from=0,to=15,by=.1)
plot(h,exp(-(h/4)^2),type="l")
lines(h,exp(-(h/2)^2),lty=2)
lines(h,exp(-(h/8)^2),lty=3)


z.vg=vgram(locs,z,N=32)
plot(h,1-exp(-(h/4)^2),type="l")
points(z.vg$centers,z.vg$stats[2,])
var(as.numeric(z))