J'aimerais savoir si quelqu'un saurait extraire une pente d'un modèle numérique de terrain s'il vous plait (MNT).
Par exemple :
J'ai mon data :
Code : Tout sélectionner
data <- unique(data.frame(INSEE = data$INSEE, LONGITUDE=long,
LATITUDE=lat))
J'ai mon fichier MNT dispo sur le net :
Code : Tout sélectionner
uu <- raster("bdalti_france_250m_L93.tif")
Ici je transforme mes coordonnées :
Code : Tout sélectionner
coordinates(data) <- c("LONGITUDE", "LATITUDE")
je spécifie la projection :
Code : Tout sélectionner
proj4string(data) = CRS("+proj=longlat +datum=WGS84")
data2 <- spTransform(data, "+init=epsg:2154")
et j'extrais l'altitude :
Code : Tout sélectionner
altitude <- raster::extract(uu,pccosp2)
J'ai cru comprendre que que le package
Code : Tout sélectionner
SDMTools
Code : Tout sélectionner
slope
Merci à tous !