A partir de mon jeu de données (points de mesures effectuées au large des côtes marocaines) je souhaiterais :
1- Cartographier les points avec des valeurs associées en ajoutant le polygone du continent, en distinguant les mesures de jour et de nuit.
2- Calculer la distance entre chaque point et le trait de côte. Ce qui induit d'avoir un polygone à haute résolution.
Voici un exemple reproductible:
Code : Tout sélectionner
library(ggplot2)
library(maptools)
library(mapdata)
library(rworldmap)
lon = c(-16.80,-16.81,-16.82,-16.83,-16.84,-16.85)
lat = c(22.93,22.91,22.90,22.89,22.87,22.86)
values = c(1.84,2.06,1.54,1.47,1.48,1.17)
day_night = c(1,1,1,3,3,3)
dataset = data.frame(lon,lat,values,day_night)
maroc.map <- data.frame(map('worldHires', c('Morocco','Western Sahara'))[c('x', 'y')])
ggplot(maroc.map,aes(x=x,y=y))+ geom_path() +
geom_point(aes(x=lon, y=lat, color=values),inherit.aes=FALSE,data=dataset,alpha=1, size=2)+
facet_grid(~day_night)+
scale_x_continuous(limits = c(min(lon)-0.5,max(lon)+0.5),"Longitude") + scale_y_continuous(limits = c(min(lat)-0.5, max(lat)+0.5),"Latitude")+
theme_bw()
Pour le moment, j'utilise la librairie ggplot et worldHires pour le polygone continent. Le résultat n'est pas très satisfaisant car le trait de côte est d'une faible résolution et il est tronqué par geom_path:
Warning messages:
1: Removed 141 rows containing missing values (geom_path).
Y a t-il des librairies / fonctions plus adaptées à mes besoins ? Et je n'ai pas encore abordé le calcul de la distance entre chaque point et le trait de côte...
Merci d'avance,