我正在尝试使用文化距离矩阵中包含的数据来描绘地理地图上的网络。 E.g:
AT BE CH CZ
AT 0 0.00276 0.148 0.109
BE 0.00276 0 0.145 0.112
CH 0.148 0.145 0 0.257
CZ 0.109 0.112 0.257 0
网络线路的起点和终点应位于不同的国家(即AT,BE,CH和CZ)。
当矩阵的相应条目低于某个阈值(例如,所有矩阵条目的平均值)时,应在国家之间描绘线。 (我想 dplyr 包可用于过滤数据,如示例http://www.gis-blog.com/flight-connection-map-with-r/中所示)
该地图包括欧亚大陆国家。我使用Trimble Data Marketplace获取Shapefile并在R中绘制地理地图,如下所示:
使用以下代码获取此地图:
> shapefile <- readOGR("directory_with_file", "name_of_file")
> shapefile_df <- fortify(shapefile)
> map <- ggplot() + geom_path(data = shapefile_df, aes(x = long, y = lat,
group = group),color = ‘black', size = .2)
> print(map)
如何使用矩阵'数据?
在此地理地图上绘制网络(网络将代表各国的文化接近度及其随时间的演变)
答案 0 :(得分:1)
您需要为您正在寻找的内容建立一个网络;我从ggmap :: geocode包中获得的坐标; 拥有网络后,您可以设置&#34; Edge&#34;用来表示你所拥有的文化距离值的参数,因为它从0.002变为0.2你必须放大它,否则你会得到非常小的线条,这段代码会帮助你上轨,你只需要添加剩余的文化距离
library(maps)
library(igraph)
df<-data.frame(from = c("at", "be", "ch", "cz"), to= c("be", "ch", "cz", "at"),
weight=c(0.02,0.145,0.257,.109))
meta <- data.frame("name"=c("at", "be", "ch", "cz"),
"lon"=c(14.55,4.46,8.227,14.4738),
"lat"=c(47.51,50.5,46.818,50.0755))
g <- graph.data.frame(df, directed=F, vertices=meta)
E(g)$color <- "brown"
E(g)$width <- E(g)$weight*10
lo <- as.matrix(meta[,2:3])
map("world", xlim = c(-8, 30),
ylim = c(35, 55), asp=1)
plot(g, layout=lo, add = TRUE, rescale = FALSE)
其他建议:在ggplot2
中绘制地图时使用+ coord_equal()