2

以此为参考，我正在尝试绘制下四十八的地图并添加图层以可视化状态之间的流动。

library(ggplot2)
library(maps)
library(geosphere) # to inter-polate a given pair of (lat,long) on the globe

# load map data for the US
all_states <- map_data("state")

# plot state map
p <- ggplot() + geom_polygon( data=all_states, 
                       aes(x=long, y=lat, group = group),
                       colour="white", fill="grey10" )

# sample origin - destination lat,long pairs
geo <- structure(list(orig_lat = c(36.17, 36.17, 36.17), 
orig_lon = c(-119.7462, -119.7462, -119.7462), dest_lat = c(33.7712, 36.17, 39.0646), 
    dest_lon = c(-111.3877, -119.7462, -105.3272)), .Names = c("orig_lat", 
"orig_lon", "dest_lat", "dest_lon"), row.names = c(NA, 3L), class = "data.frame")

#> geo
#  orig_lat  orig_lon dest_lat  dest_lon
#1    36.17 -119.7462  33.7712 -111.3877
#2    36.17 -119.7462  36.1700 -119.7462
#3    36.17 -119.7462  39.0646 -105.3272

# list to hold a dataframe of interpolated points for each origin-destination pair
list_lines <- list()

# use the geosphere package's gcIntermediate function to generate 50 interpolated  
# points for each origin-destination pair
for (i in 1:3) {
  inter <- as.data.frame(gcIntermediate(c(geo[i,]$orig_lon, geo[i,]$orig_lat), 
                                        c(geo[i,]$dest_lon, geo[i,]$dest_lat), 
                                        n=50, addStartEnd=TRUE))
  list_lines[i] <- list(inter)
  p <- p + geom_line( data = list_lines[[i]], aes(x = lon, y = lat), color = '#FFFFFF')
}
p

这是我尝试打印情节时得到的

p
Error in eval(expr, envir, enclos) : object 'lon' not found

我试图调试它，发现这有效

p + geom_line( data = list_lines[[1]], aes(x = lon, y = lat), color = '#FFFFFF')

但是为第二个列表元素添加另一层会破坏它，但这是我对 R 和 ggplot 的有限知识所能得到的！

4

3 回答 3

4

gcIntermediate返回不同的列名（因为 i=2 的起点和终点相同）：

for (i in 1:3) {
  inter <- as.data.frame(gcIntermediate(c(geo[i,]$orig_lon, geo[i,]$orig_lat), 
                                        c(geo[i,]$dest_lon, geo[i,]$dest_lat), 
                                        n=50, addStartEnd=TRUE))
  print(head(inter, n=2))
}
     lon   lat
1 -119.7 36.17
2 -119.6 36.13
      V1    V2
1 -119.7 36.17
2 -119.7 36.17
     lon   lat
1 -119.7 36.17
2 -119.5 36.24

以下几行应该有效：

for (i in 1:3) {
  inter <- as.data.frame(gcIntermediate(c(geo[i,]$orig_lon, geo[i,]$orig_lat), 
                                        c(geo[i,]$dest_lon, geo[i,]$dest_lat), 
                                        n=50, addStartEnd=TRUE))
  names(inter) <- c("lon", "lat")
  p <- p + geom_line(data=inter, aes(x=lon, y=lat), color='#FFFFFF')
}
于 2012-08-05T09:04:23.893 回答
1

令我感到奇怪的是，您以两种不同的方式引用经度：long在脚本的开头和lon结尾。如果您希望多个 geom 一起工作，则需要使这些名称保持一致。

此外，几乎不需要使用 for 循环添加相同的 geom。只需添加一个geom_line并使用color美学来绘制多条线。

于 2012-08-05T08:32:49.163 回答
0

有一个非常简单的解决方案，使用ggplot2. 这里有关于如何在R使用ggplot2中绘制流图的简单教程。

p + 
  geom_segment(data = geo, aes(x = orig_lon,    y = orig_lat, 
                               xend = dest_lon, yend = dest_lat,
                               color="#FFFFFF")) +  coord_equal()

在此处输入图像描述

于 2016-06-15T17:54:49.790 回答