使用R在地图中绘制电子邮件流

使用R在地图中绘制电子邮件流,r,email,map,plot,R,Email,Map,Plot,我想在这个世界上策划电子邮件往来。例如,我从以下国家收到了以下数量的电子邮件,我住在美国 recievedcountry <- c("India", "China", "France", "Chile", "Australia", "Chad", "Nepal", "Burma") rfrequency <- c(12, 20, 5, 2, 12, 1, 3, 2) # frequency of emails sendcountry <

我想在这个世界上策划电子邮件往来。例如,我从以下国家收到了以下数量的电子邮件,我住在美国

recievedcountry <- c("India", "China", "France", "Chile", "Australia", "Chad",
                    "Nepal", "Burma")
rfrequency <- c(12, 20, 5, 2, 12, 1, 3, 2) # frequency of emails 
sendcountry <- c("Canda", "USA", "France", "India", "China", "Japan")
sfrequency <- c(14, 108, 12, 15, 18, 4)
以下是我的假设模型:


这是以前考虑过的。看看这篇博文,它非常仔细地描述了如何构建这样一幅地图


基于上面提供的地址(谢谢!)链接,我可以回答我自己的问题,除了与电子邮件数量成比例的行思考

下面是代码:我使用随机数据,没有搜索每个位置的确切经度和纬度

library(maps)
library(geosphere)
map("world", col = "green4", bg="#F5FFFA", lwd=0.05)
myposition <- c(-74, 40) # my position (where I am opening emails)

rlong <- c(75, 105, 135, - 10.2,  45.2, -30.4, 105, 35, -150, 
   10.2,  145.2, 30.4) # received lat
rlat <- c(30, 43, 23, 12, 68, 55.6, 30, 43, 23, 12, 68, 55.6) # received long
nrecived <- c(4, 10, 5, 2, 4, 10, 4, 10, 5, 2, 4, 10 )     # number of email received
slong <- c(85, 85, 55, -40.2,  45.2, -30.4,45, 95, 55, 40.2,  55.2, 60.4 ) # send lat
slat <- c(10, 43, 13, 12, 68, 55.6,10, 43, 13, 12, 68, 55.6 ) # send long
nsend <- c(4, 10, 5, 2, 4, 10, 4, 10, 5, 2, 4, 10 )     # number of email send

mydf <- data.frame (rlat, rlong, nrecived, slat, slong, nsend)

for (i in 1: length (mydf) ) {
            send <- gcIntermediate(c(mydf[i,]$slong, mydf[i,]$slat), c(-74, 40), 
                   n=100, addStartEnd=TRUE)
            lines (send, col = "blue", lwd = mydf[i, "nsend"]) # edited 
                                                                following suggestion 
            received <- gcIntermediate(c(mydf[i,]$rlong, mydf[i,]$rlat), c(-74, 40), 
             n=100,  addStartEnd=TRUE)
            lines (received , col = "red", lwd = mydf[i, "nrecived"])
            }
库(地图)
图书馆(地球圈)
地图(“世界”,col=“green4”,bg=“#F5FFFA”,lwd=0.05)
myposition这是一个解决方案。可以制作平面贴图或球面贴图:

如果您将
lwd=1.2
更改为
lwd=mydf[i,“nsend”]
,它是否解决了您的测线厚度问题?您可以将当前的lwd重新缩放到更具功能的位置。你可以除以0.8或者类似的数字。
library(maps)
library(geosphere)
map("world", col = "green4", bg="#F5FFFA", lwd=0.05)
myposition <- c(-74, 40) # my position (where I am opening emails)

rlong <- c(75, 105, 135, - 10.2,  45.2, -30.4, 105, 35, -150, 
   10.2,  145.2, 30.4) # received lat
rlat <- c(30, 43, 23, 12, 68, 55.6, 30, 43, 23, 12, 68, 55.6) # received long
nrecived <- c(4, 10, 5, 2, 4, 10, 4, 10, 5, 2, 4, 10 )     # number of email received
slong <- c(85, 85, 55, -40.2,  45.2, -30.4,45, 95, 55, 40.2,  55.2, 60.4 ) # send lat
slat <- c(10, 43, 13, 12, 68, 55.6,10, 43, 13, 12, 68, 55.6 ) # send long
nsend <- c(4, 10, 5, 2, 4, 10, 4, 10, 5, 2, 4, 10 )     # number of email send

mydf <- data.frame (rlat, rlong, nrecived, slat, slong, nsend)

for (i in 1: length (mydf) ) {
            send <- gcIntermediate(c(mydf[i,]$slong, mydf[i,]$slat), c(-74, 40), 
                   n=100, addStartEnd=TRUE)
            lines (send, col = "blue", lwd = mydf[i, "nsend"]) # edited 
                                                                following suggestion 
            received <- gcIntermediate(c(mydf[i,]$rlong, mydf[i,]$rlat), c(-74, 40), 
             n=100,  addStartEnd=TRUE)
            lines (received , col = "red", lwd = mydf[i, "nrecived"])
            }