Carte d'une région du monde (exemple : UK)

install.packages("raster");library(raster)

getData('ISO3') # pour voir le code de l'United Kingdom

monde <- getData('GADM', country='GBR', level=2) # Charger la carte

# afficher le fond de carte (vérification)

plot(monde, col="#FFB765")

# TRACER L'ECHELLE (attention de régler sa position en longitude, latitude)

# d : distance en km

# xy = c(longitude, latitude)  ==> position de la barre d'échelle

scalebar(d = 100, xy = c(-11, 54.8), type = "bar", below = "km",lwd = 4, divs = 2, col = "black", cex = 0.75, lonlat = T)

############# FONCTION northarrow qui tracera la rosace

northarrow <- function(loc, size, bearing = 0, cex = 1) {

cols <- rep(c("white", "black"), 8)

# Coordonnees des polygones de la rose des vents

radii <- rep(size/c(1, 4, 2, 4), 4) ; x <- radii[(0:15) + 1] * cos((0:15) * pi/8 + bearing) + loc[1] ; y <- radii[(0:15) + 1] * sin((0:15) * pi/8 + bearing) + loc[2]

# Trace des polygones

for (i in 1:15) {x1 <- c(x[i], x[i + 1], loc[1]);y1 <- c(y[i], y[i + 1], loc[2]);polygon(x1, y1, col = cols[i])}

# Trace du dernier polygone

x1 <- c(x[16], x[1], loc[1]);y1 <- c(y[16], y[1], loc[2]);polygon(x1, y1, col = cols[16])

# Rajout des lettres

b <- c("E", "N", "O", "S")

for (i in 0:3) {text((size + par("cxy")[1]) * cos(bearing + i * pi/2) + loc[1],(size + par("cxy")[2]) * sin(bearing + i * pi/2) + loc[2],b[i + 1], cex = cex)}

}

# =========================================

# A MODIFIER EN RÉGLANT LA POSITION DE LA ROSACE

# loc = c(longitude, latitude)

# size = la taille de l'indication

northarrow(loc = c(-6, 54.7), size = 0.7, cex = 0.7)