Bonjour,

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 : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
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:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
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,