R объединить choropleth участок с координатами планировавших

Question

Допустим, у меня есть два набора данных, countries и cities:

set.seed(42) 

countries <- data.frame(country = c("DEU", "FRA", "ITA", "POL"), 
         intensity = abs(rnorm(4))) 

cities <- data.frame(city = c("Berlin", "Paris", "Roma", "Warszawa"), 
        longitude = c(13.408333, 2.351667, 12.483333, 21.033333), 
        latitude = c(52.518611, 48.856667, 41.883333, 52.216667)) 

countries описывает уровень странового набора данных интенсивностей. cities - это набор данных, содержащий координаты. До сих пор я использовал rworldmap для построения интенсивности:

library(rworldmap) 

map <- joinCountryData2Map(countries, joinCode = "ISO3", nameJoinColumn = "country") 

mapCountryData(map, nameColumnToPlot = "intensity", mapRegion = "Europe") 

Как я могу объединить choropleth карту данных на уровне стран с наложением точек в соответствии с набором данных городов?

Answer 1

Если я правильно понимаю, вы можете добавлять точки на карту с помощью points:

map <- joinCountryData2Map(countries, joinCode = "ISO3", nameJoinColumn = "country") 
mapCountryData(map, nameColumnToPlot = "intensity", mapRegion = "Europe") 
points(cities$longitude, cities$latitude, col = "black", cex = 1.5)