Hello everyone. 
I would kindly request your help concerning how R converts data between 
different structrures.

In the following example please keep attention on the following two
1) 
I create 
f <- GaussRF(x=x, y=y, model=model, grid=TRUE,param=c(mean, variance, nugget, 
scale, alpha))
with 
image(x,y,f)
   and
image(f)
I get exactly the same image.

then
2)I try to store f inside a raster layer using:
r <- setValues(r,as.matrix(f))

then comes the tricky part. I am trying to use image again and get the same 
output as the ouput I got by image(f).
I tried
image(as.matrix(getValues(r)))
which give a completely different output. 
image function expects a matrix... but for some reason as.matrix(getValues(r)) 
returns a huge vector :(
also tried
as.matrix(getValues(r),ncol=ncol(r),nrow=nrow(r)) to force as.matrix to return 
the appropriate matrix but this also failed and I only got back a vector again

Could you please let me understand why this might happening?
Best Regards
Alex

P.S Below you will find part of my code.



...................(lines omitted, declarations)
x <- seq(1, dimx, step)
y <- seq(1, dimy, step)
f <- GaussRF(x=x, y=y, model=model, grid=TRUE,param=c(mean, variance, nugget, 
scale, alpha))

# image(x, y, f) Displays the matrix
# f is a 2x2 matrix of dimension (x,y)

r <- raster(nrow=dimx, ncol=dimy)
r <- setValues(r,as.matrix(f))
# getValues(r)



      
	[[alternative HTML version deleted]]

