Re: [R] GRASS raster data processing
Perhaps reading ?readRAST6 would help? It does say: vname: A vector of GRASS 6.0 raster file names so you could write: Kar - readRAST6(c(Incis_Kar, DEM_Kar), plugin=FALSE) and plot(Incis_Kar ~ DEM_Kar, data=Kar) should get you there - in spearfish: sp - readRAST6(c(erosion1, elevation.dem), plugin=FALSE) plot(erosion1 ~ elevation.dem, data=sp) works for me. (Aside - using @ in S4 classes should be avoided unless you both know what you are doing and need to be doing it, use the provided access mechanisms, here treating a Spatial*DataFrame as an ordinary data frame.) Roger Nikos Alexandris wrote: Maayt: I just imported two raster maps into R using the SPGRASS6 package, one containing elevation data and the other containing an erosion index: Kar_inc -readRAST6(Incis_Kar, plugin=FALSE) Kar_dem - readRAST6(DEM_Kar, plugin=FALSE) I just wanted to make a xy plot of erosion parameter vs elevation. How does this work? I don't get how to handle SpatialGridDataFrames... Maarten, you can check some web-pages with respect to GRASS R [1][2]. This is a question for grass-stats actually [3]. A quick answer: check the structure of the newly created object and you will find that the numbers are to be found in (e.g. for Incis_Kar ) kar_...@data$incis_kar. So use @data and the $ to access a slot. Kind regards, Nikos --- [1] http://grass.osgeo.org/wiki/R [2] http://grass.ibiblio.org/statsgrass/index.php#grassR [3] http://grass.osgeo.org/statsgrass/index.php -- View this message in context: http://www.nabble.com/GRASS-raster-data-processing-tp23981740p23996501.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org 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.
[R] GRASS raster data processing
Hi, I just imported two raster maps into R using the SPGRASS6 package, one containing elevation data and the other containing an erosion index: Kar_inc -readRAST6(Incis_Kar, plugin=FALSE) Kar_dem - readRAST6(DEM_Kar, plugin=FALSE) I just wanted to make a xy plot of erosion parameter vs elevation. How does this work? I don't get how to handle SpatialGridDataFrames... Thanks a lot Maarten -- View this message in context: http://www.nabble.com/GRASS-raster-data-processing-tp23981740p23981740.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org 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.
Re: [R] GRASS raster data processing
Maayt: I just imported two raster maps into R using the SPGRASS6 package, one containing elevation data and the other containing an erosion index: Kar_inc -readRAST6(Incis_Kar, plugin=FALSE) Kar_dem - readRAST6(DEM_Kar, plugin=FALSE) I just wanted to make a xy plot of erosion parameter vs elevation. How does this work? I don't get how to handle SpatialGridDataFrames... Maarten, you can check some web-pages with respect to GRASS R [1][2]. This is a question for grass-stats actually [3]. A quick answer: check the structure of the newly created object and you will find that the numbers are to be found in (e.g. for Incis_Kar ) kar_...@data$incis_kar. So use @data and the $ to access a slot. Kind regards, Nikos --- [1] http://grass.osgeo.org/wiki/R [2] http://grass.ibiblio.org/statsgrass/index.php#grassR [3] http://grass.osgeo.org/statsgrass/index.php __ R-help@r-project.org 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.