Dear all,
A very basic terrain calculated as a matrix from Spatial Points Patterns:
#interpolate using the akima package
library(akima)
terrain=interp(ppoints$x,ppoints$y,ppoints$marks,xo=x0,yo=y0, linear=F)
> class(terrain)
[1] "list"
> class(terrain$x) #these are the x-coord i.e: [1...1000]
[1] "numeric"
> class(terrain$y)#these are the y-coord i.e: [1...1000]
[1] "numeric"
> class(terrain$z)#these are the height
[1] "matrix"
I would like to export this "terrain" (list object) to GRASS.
Reading the documentation that I found for interfacing between GRASS 6 and R
(library (spgrass6)). I understood the following (I already managed to
configured properly the mapset and so on):
1. I have to create a SpatialGridDataFrame object before to write into a Raster
file.
2. In order to create a SpatialGridDataFrame I need a GridTopology.
3. After doing it, I am allow to write the file into GRASS using the function
writeRAST6, like: writeRAST6 (spterrain, "maps.dem").
How to create the SpatialGridDataFrame and/or the GridTopology? I am doing
something wrong.
Thanks,
Javier Hidalgo.
_________________________________________________________________
[[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.
