Code : Tout sélectionner
require (raster)
require (SDMTools)
require(adehabitat)
library(ade4) # livre de recettes analyses multivariées
setwd("C:/Users/faye/......") # Définit répertoire de travail
##Clippage des couches##.........faye/Desktop/6K", pattern='.tif', full.names=TRUE)
Bio.LIG
rasterBio=list();
length(rasterBio)=length(Bio.LIG)
for (i in 1 : length(Bio.LIG)) {
rasterBio[[i]]=raster(Bio.LIG[[i]])
cat(i,"\n")
}
afrique<- shapefile("C:/Users/....../shp ocean afrique/ocean afrique2.shp")
ex<-extent(1,20,-10,12) # extension de l'aire d'interet
Maskocean<- crop(afrique,ex,snap= "out")
# masquer toutes les couches
maskBio<-list(); length(maskBio)=length(rasterBio)
cropBio<-list();length(cropBio)=length(rasterBio)
for (i in 1: length(rasterBio)){
maskBio[[i]]<-mask(rasterBio[[i]], Maskocean,inverse=T, snap='out')# en 1er
cropBio[[i]]<- crop(maskBio[[i]], ex, snap= "out")#en sd
cat(i,"\n")
plot (cropBio[[i]])
}
# exporter les rasters decoupes
ascBio=list();length(ascBio)=length(cropBio)
for (i in 1: length(cropBio)){
write.asc(asc.from.raster( cropBio[[i]]),paste("C:/Users/......./bio",i,".asc",sep="" ))
ascBio[[i]]<-read.asc(paste("C:/Users/......../bio",i,".asc",sep=""))
cat(i,"\n")
}