GRASS 6.4.0RC5 (geol):~ > R
R version 2.9.2 (2009-08-24)
Copyright (C) 2009 The R Foundation for Statistical Computing
[....]
> library(spgrass6) # chargement du package spgrass6 qui permet d'interagir "en direct" avec Grass.
[....]
# lecture d'un raster et d'un MNT provenant de Grass dans la variable montrait
> G <- gmeta6() # paramètres de Grass, région etc.
> montrait <- readRAST6(c("geologie", "mnt"),cat=c(TRUE, FALSE), ignore.stderr=TRUE,plugin=NULL)
# puisque Grass utilise des rasters "intelligents" (catégories), il y a moyen d'obtenir le nombre de cellules dans chaque catégories (cat) qui sont ici les formations géologiques
> table(montrait$geologie)
ACO AIFA ANS BAU ....
20116 2067 4485 12733 ....
# on peut visualiser dans R les rasters de Grass
> image(montrait, attr = 2, col = terrain.colors(20))
# l'auteur s'amuse ensuite à créer un diagramme "en boîte" (box plot) de la répartition des formations géologiques en fonction de l'altitude
> boxplot(montrait$mnt ~ montrait$geologie, medlwd = 1)
#création d'un nouveau MNT
> montrait$nvmnt <- montrait$mnt * 2
# renvoi du nouveau MNT dans Grass
> writeRAST6(montrait, "mnt2", zcol="nvmnt", ignore.stderr=TRUE)
# on quitte R
> q()
GRASS 6.4.0RC5 (geol):~ > r.info map=mnt2@moi
+----------------------------------------------------------------------------+
| Layer: mnt2@moi Date: Mon Mar 22 22:12:59 2010 |
| Mapset: moi Login of Creator: yo |
| Location: geol |
|[....] |
| Type of Map: raster Number of Categories: 255 |
| Data Type: FCELL |
| Rows: 191 |
| Columns: 270 |
| Total Cells: 51570 |
| Projection: Lambert Conformal Conic |
| N: 118892.8197047 S: 107481.49881847 Res: 59.74513553 |
| E: 178100.24901493 W: 161969.06242183 Res: 59.74513553 |
[....]
GRASS 6.4.0RC5 (geol):~ >
