使用带有geom_point的for循环向现有的ggplot对象添加点

Question

我是R的新手，正在尝试创建一张我的国家的地图，上面覆盖着点(代表城市，点按人口大小)。我已经用geom_sf()创建了地图，并将其存储为地图。

这是我试图复制的代码的一个示例，其中cities是一个包含列state、pop、lat和long (没有缺失值)的df：

map <- map +
  geom_point(aes(x = cities$long[1], y = cities$lat[1], size = cities$pop[1])) +
  geom_point(aes(x = cities$long[2], y = cities$lat[2], size = cities$pop[2])) +
  geom_point(aes(x = cities$long[3], y = cities$lat[3], size = cities$pop[3])) +
  geom_point(aes(x = cities$long[4], y = cities$lat[4], size = cities$pop[4])) +
  geom_point(aes(x = cities$long[5], y = cities$lat[5], size = cities$pop[5]))

这是我认为可以工作的循环，但它似乎会在每次迭代中覆盖geom_point，因为我的地图只是以最后一个点结束。当我在每次迭代中输出一张地图时，它只显示最近的点。

or(i in 1:nrow(cities)){
  if(!is.na(cities$lat[i]) && !is.na(cities$long[i])){
    map <- map +
      geom_point(aes(x = cities$long[i], y = cities$lat[i], size = cities$pop[i]))
  }
}

我希望我已经提供了足够的信息-如果没有，请让我知道。谢谢！

Answer 1

没有必要在循环中执行此操作。

没有你的数据就不能重现，但是假设你的城市对象看起来像：

cities <- data.frame(name = .., 
                     pop = .., 
                     lat = ..., 
                     long = ...)

假设您的底图是一个geom_sf对象，那么这就是您需要的所有代码。

cities_spatial <- cities %>% 
  st_as_sf(coords = c("lat","long"), crs = 4326)

map <- map +
  geom_sf(data = cities_spatial, aes(size = pop))

选项crs = 4326是为了确保您的数据以经纬度坐标的标准WGS84格式进行投影。如果底图map使用不同的投影进行投影，则可能需要使用sf::st_transform(...)对其进行调整以与经纬度数据相匹配