Sextante et R

Opérateurs sur données Raster | terrain

##Cartographie avec R=group
##Raster=raster
##Slope= output raster
##Aspect= output raster
##Hillshade= output raster
library(raster)
rMNT_RO=Raster[[1]]
proj_RGF93=CRS("+init=epsg:2154")
projection(rMNT_RO)=proj_RGF93
rSLOPE_RO=terrain(rMNT_RO,opt="slope")
rASP_RO=terrain(rMNT_RO,opt="aspect")
rHILL_RO=hillShade(terrain(rMNT_RO,opt="slope"),terrain(rMNT_RO,opt="aspect"))
Slope=rSLOPE_RO
Aspect=rASP_RO
Hillshade=rHILL_RO

Agrégation d’un raster | aggregate

##Cartographie avec R=group
##Raster=raster
##nombre=number 20
##Output= output raster
library(raster)
rMNT_RO=Raster[[1]]
proj_RGF93=CRS("+init=epsg:2154")
projection(rMNT_RO)=proj_RGF93
rMNT_RO_LowRes<-aggregate(rMNT_RO,fact=nombre) 
Output=rMNT_RO_LowRes

Opérateurs sur données Raster | focal

##Cartographie avec R=group
##Raster=raster
##nombre=number 9
##Raster_Lisse= output raster
library(raster)
rSLOPE_RO=Raster[[1]]
proj_RGF93=CRS("+init=epsg:2154")
FenetreLissage=matrix(1,nombre,nombre)
rSLOPE_RO_med<- focal(rSLOPE_RO,fun=median,w=FenetreLissage)
Raster_Lisse=rSLOPE_RO_med

Accès aux données slot des vecteurs | slot

##Cartographie avec R=group
##Layer=vector
##showplots
library(sp)
vOCSOL_BO=Layer
listPolygons=slot(vOCSOL_BO, "polygons")
vecArea=sapply(listPolygons, slot, "area")
hist(vecArea,breaks=10,xlab="Aire des polygones",ylab="Frequence",main="")

Manipulation de la table attributaire de données vectorielles | $

##Cartographie avec R=group
##Layer=vector
##Output= output vector
library(sp)
vOCSOL_BO=Layer
listPolygons=slot(vOCSOL_BO, "polygons")
vecArea=sapply(listPolygons, slot, "area")
vOCSOL_BO$Area=vecArea
vOCSOL_BO$AreaHa=vecArea/10000
vOCSOL_BO_lines=as(vOCSOL_BO,"SpatialLines")
vecPerimeters=SpatialLinesLengths(vOCSOL_BO_lines)
vOCSOL_BO$Perimeter=vecPerimeters
vOCSOL_BO$Ratio=vecPerimeters/vecArea
Output=vOCSOL_BO

Création d’une zone tampon | gBuffer

##Cartographie avec R=group
##Layer=vector
##nombre=number 30
##Un_buffer_par_segment=boolean
##Vecteur_Buffer= output vector
library(rgeos)
vHYDRO_RO=Layer
vHYDRO_RO_buffer=gBuffer(vHYDRO_RO,width=nombre,byid=Un_buffer_par_segment)
if(!Un_buffer_par_segment)
vHYDRO_RO_buffer=SpatialPolygonsDataFrame(vHYDRO_RO_buffer,
                                          data=data.frame(1),
                                          match.ID=F)
if(Un_buffer_par_segment)
vHYDRO_RO_buffer=SpatialPolygonsDataFrame(vHYDRO_RO_buffer,
                                          data=data.frame(vHYDRO_RO),
                                          match.ID=F)
Vecteur_Buffer=vHYDRO_RO_buffer

Calcul de l’aire de chaque cellule d’un MNT | surfaceArea

##Cartographie avec R=group
##Layer=vector
##nombre=number 30
##Un_buffer_par_segment=boolean
##Raster_SurfaceArea= output raster
library(sp)
library(raster)
rMNT_RO=Raster[[1]]
proj_RGF93=CRS("+init=epsg:2154")
projection(rMNT_RO)=proj_RGF93
rMNT_RO_sp=as(rMNT_RO,"SpatialGridDataFrame")
rSA_RO=surfaceArea(rMNT_RO_sp,byCell=T)
rSA_RO=raster(rSA_RO)
Raster_SurfaceArea=rSA_RO