05 - Carte 3D et interactives en langage R

library(raster)  # vérifier aussi que rgdal est installé

# Charger les contours administratifs de la Grèce

adm_gr <- getData('GADM', country='GRC', level=1)

# Charger les données altimétriques de la grèce

srtm <- getData('SRTM', lon=23, lat=36) # srtm : meilleure résolution que alt

# Réduire les contours administratifs à une région

adm <- adm_gr[match(toupper(adm_gr$NAME_1[7]),toupper(adm_gr$NAME_1)),] # Péloponnèse

# Réduire les données altimétriques à une zone

vatikasrtm <- crop(srtm,extent(c(xlim=c(23,23.3),ylim=c(36.422,36.54))))

# Lisser les données altimétriques

vatikasrtm[is.na(vatikasrtm[])] <- minValue(vatikasrtm)

vatikasrtm_mask_brut <- disaggregate(vatikasrtm, fact=17)

FenetreLissage=matrix(1,17,17) # Création de la fenêtre de lissage

vatikasrtm_mask<- focal(vatikasrtm_mask_brut,fun=mean,w=FenetreLissage,pad=T,padValue=0)  # Lissage avec filtre mean ou median (fun)médian

#vatikasrtm_mask <- mask(vatikasrtm_mask ,mask= adm, inverse=F); # optionnel, supprimer les pixels hors des contours adm

# Définir une palette de couleurs

liste_colors <- c("#BCF5A9","#36a013","#E58544","#C6753F","#B6652F")

colfunc<-colorRampPalette(liste_colors,interpolate="spline")

colors <- colfunc(1000)

colors[1] <- "#B1FAF3"

# Visualiser

plot(vatikasrtm_mask,col=colors)