[R-sig-Geo] Uncondistional nonstationary RandomFields simulation

[Babak Naimi]

Dear list,

I'm trying to generate a 2-dimensional nonstationary field using GaussRF in RandomFields package. I have problem on how to define the covariance model. I tried to use the examples provided in the manual, but unfortunately I couldn't figure out how it works for my case. I was wondering if anyone can help by providing some examples.

Using the following command it seems that two models ("nonstWM" and "stp") are available for nonstationary fields:

> GetModelNames("nonstationary")
[1] "nonstWM" "stp"

x <- y <- seq(0, 100, 1) 
model <- ?

z <- GaussRF(x=x, y=y, model=model, grid=TRUE, method=?)

Best regards,
Babak