Hello everyone.
I am trying to build up understanding in R by trying to develop just some
simple
scenarios.
I would like to explain you what I am trying to do and what I did so far.
I would like to put inside a RasterLayer (raster package) a Gaussian field (for
given covariance) using grf function (geoR package)
1. First I created a Raster Layer object
r <- raster() # Default values are ok
2. Then I set some values to test how if setValues worked
r <- setValues(r,1:ncell(r)/100) # Every cell of the RasterLayer takes as data
its cell number/100.
3. Tested that 2 works with
getValuesBlock(r,1,nrow(r),1,ncol(r))
4. Then I tried to generate a Gaussian random field for given covariance
parameters using grf.
temp <- grf(1,nx=nrow(r),ny=ncol(r), cov.pars=c(1, .25)) # Note: The first
input
parameter in grf is locations. At first I tried nrow(r)*ncol(y) but I got the
ouput that this will produce 15,5 Gb of data :(. I just wanted to created a x*y
array of this data. Then I tried to pass the nrow(r) and ncol(r) to the nx and
ny arguments. Now this seems ok to me.
5. Final step is to set the values created in 4th step to the RasterLayer
object
that was created at first step.
r <- setValues(r,temp) # Unfortunately this ends with the message: Error in
setValues(r, temp) : values must be a vector
and at this point I am completely stuck :(
I would like to request some of your kindness to help me unstuck
Best Regards
Alex.
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
