[英]Adding points to represent cities on a World Map using rworldmap
我想制作一个包含特定国家/地区的世界地图,但如果可能的话,我还想指出代表这些国家/地区中特定城市的精确点。
我一直在使用rworldmap
包并与其他人进行实验,但是我很难理解如何确定特定的城市位置以及在哪里包括mapbubbles
。
也许我可以添加一条代表圣地亚哥的线,但是我不知道它会去哪里:
mapBubbles(dF="CEAMap", nameX = "-117.16", nameY = "32.71", nameZSize = "CEAMap",)
我一直在阅读开发人员的PDF,老实说我真的不明白如何实现地图气泡以及像nameZsize这样的东西首先应该包含的内容。
我使用本指南来学习如何突出显示国家/地区,但就精确定位城市而言,我似乎是一个人,因为我无法真正理解开发人员的PDF。 如何在R中填写特定国家/地区的世界地图?
到目前为止,这是我的完整代码:
library(rworldmap)
library(ggmap)
library(maptools)
library(maps)
theCountries <- c("USA",
"CAN", "DEU", "FRA", "IND",
"GBR", "NLD", "ITA",
"CHN", "KOR", "JPN",
"ESP", "PRT", "RUS",
"NOR", "SGP", "AUS",
"CHL", "MEX", "PHL", "RWA",
"JOR", "HND", "PAN", "THA", "DOM",
"ZAF", "TUR", "CHE", "FIN",
"SEN", "BOL", "OMN", "PAK", "CMR", "MUS", "BEL", "MYS",
"UAE", "BRA", "MLI", "MOZ", "NAM", "EGY", "ARG", "UKR", "ZMB", "KEN",
"VNM", "NGA", "DNK", "IRN", "AFG")
# These are the ISO3 names of the countries you'd like to plot in red
CEAMap <- data.frame(country = c("USA",
"CAN", "DEU", "FRA", "IND",
"GBR", "NLD", "ITA",
"CHN", "KOR", "JPN",
"ESP", "PRT", "RUS",
"NOR", "SGP", "AUS",
"CHL", "MEX", "PHL", "RWA",
"JOR", "HND", "PAN", "THA", "DOM",
"ZAF", "TUR", "CHE", "FIN",
"SEN", "BOL", "OMN", "PAK", "CMR", "MUS", "BEL", "MYS",
"UAE", "BRA", "MLI", "MOZ", "NAM", "EGY", "ARG", "UKR", "ZMB", "KEN",
"VNM", "NGA", "DNK", "IRN", "AFG"),
involvement = c(1,
2, 2, 2, 2,
3, 3, 3,
4, 4, 4,
5, 5, 5,
6, 6, 6,
7, 7, 7, 7,
8, 8, 8, 8, 8,
9, 9, 9, 9,
10, 10, 10, 10, 10, 10, 10, 10,
11, 11, 11, 11, 11, 11, 11, 11, 11, 11,
12, 12, 12, 12, 12))
# CEAMap is a data.frame with the ISO3 country names plus a variable to
# merge to the map data
CEAcountries <- joinCountryData2Map(CEAMap, joinCode = "ISO3",
nameJoinColumn = "country")
# This will join your CEAMap data.frame to the country map data
mapCountryData(CEAcountries, nameColumnToPlot="country",
catMethod = "categorical",
mapTitle='CEA Locations',
missingCountryCol = gray(.8))
我想提出一些要点,以突出显示这些国家/地区中的城市。
我想提出一些要点,以突出显示这些国家/地区中的城市。
那么mapBubbles()
是错误的函数。 您可以像在常规R图中一样添加points()
。
在获取以下城市的位置数据之前:
library(maps)
data("world.cities")
plotcities <- subset(world.cities, capital == 1)
我不知道您对哪个城市感兴趣,所以我只带了首都
mapCountryData(CEAcountries, nameColumnToPlot="country",
catMethod = "categorical",
mapTitle='CEA Locations',
missingCountryCol = gray(.8), addLegend = FALSE)
points(plotcities$long, plotcities$lat, pch = 18, col = "black")
我建议尽管看一下tmap
软件包。 他们制作了更好的地图。 ggplot2
也可能是一个更好的选择。
编辑:要使用您提供的城市,只需选择相应的:
plotcities <- subset(world.cities,
name %in% c("Cologne", "Chennai", "Denver", "Madrid", "Manila", "San Diego", "Seattle", "Shanghai")
& country.etc %in% c("Germany", "USA", "Spain", "China", "Philippines", "India"))
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.