ggplot2地图上的多个站点

Question

我正在尝试将所有网站绘制在ggplot地图上。我有5个位置，每个位置有多个站点，但是当我尝试绘制我的所有站点时，它每个站点只显示一个点而不是每个站点一个点。我有34个站点和5个位置，所以我需要一个34分而不是5分的地图！任何帮助将非常感激！非常感谢

这些是y和x作为我的long和lat的数据。

  

Location PAR y x Guam GFA 201.09952 144.8786111 13.495 Guam GFA2 179.04171 144.6597222 13.41638889 Guam GFA3 67.66379 144.2761111 13.47333333 Guam GFB1 201.09952 144.7105556 13.31416667 Guam GFB2 179.04171 144.655 13.50194444 Guam GFB3 67.66379 144.8697222 13.37472222

我需要为关岛绘制我的所有6个网站。这是我正在使用的代码和最终输出

map.world <- map_data(map="world")
gg <- ggplot()
gg <- gg + theme(legend.position="none")
gg <- gg + geom_map(data=map.world, map=map.world, aes(map_id=region, x=long, y=lat), fill="white", colour="black", size=0.25) + theme_bw()
gg

par<-read.csv("parmap.csv", header=T)
head(par)
g<-gg+ geom_polygon() + 
  geom_point(data=par, aes(x = y, y = x, color=PAR)) +theme_minimal()
g

Answer 1

您可能需要在感兴趣的点附近缩放地图，例如在这种情况下为“关岛”。

一种方法可能如下

library(ggplot2)
library(ggmap)

#get the map
center = paste(min(df$lat)+(max(df$lat)-min(df$lat))/2, 
               min(df$lon)+(max(df$lon)-min(df$lon))/2, sep=" ")
map <- get_map(location = center, maptype = "roadmap", zoom = 8, source = "google", color="bw")
p <- ggmap(map)

#plot points on the map
plt <- p +
  geom_point(data=df, aes(x = lon, y = lat, color = PAR)) +
  scale_color_gradientn(colours = rainbow(10)) +
  theme_minimal()
plt

更新：我稍微修改了您的示例数据。你可以注意到我在样本数据中又添加了一个位置来绘制地图上的所有7个点（即一个位置的6个站点和另一个位置的1个站点）。

输出图是：

注意：如果您收到geocode failed with status OVER_QUERY_LIMIT错误，则this链接可能会有所帮助。

示例数据：

df <- structure(list(Location = c("Guam", "Guam", "Guam", "Guam", "Guam", 
"Guam", "N_Mariana_Islands"), sites = c("GFA", "GFA2", "GFA3", 
"GFB1", "GFB2", "GFB3", "NMI1"), PAR = c(201.09952, 179.04171, 
67.66379, 201.09952, 179.04171, 67.66379, 100), lon = c(144.8786111, 
144.6597222, 144.2761111, 144.7105556, 144.655, 144.8697222, 
145.192522), lat = c(13.495, 13.41638889, 13.47333333, 13.31416667, 
13.50194444, 13.37472222, 14.150817)), .Names = c("Location", 
"sites", "PAR", "lon", "lat"), class = "data.frame", row.names = c(NA, 
-7L))
#           Location sites       PAR      lon      lat
#1              Guam   GFA 201.09952 144.8786 13.49500
#2              Guam  GFA2 179.04171 144.6597 13.41639
#3              Guam  GFA3  67.66379 144.2761 13.47333
#4              Guam  GFB1 201.09952 144.7106 13.31417
#5              Guam  GFB2 179.04171 144.6550 13.50194
#6              Guam  GFB3  67.66379 144.8697 13.37472
#7 N_Mariana_Islands  NMI1 100.00000 145.1925 14.15082