В R существует пакет, называемый «потоками» (подробности здесь) , который позволяет построить карту взаимодействия потоков (с помощью функции plotMapDomFlows () из квадратной матрицы потоков. Вот пример с данными, содержащимися в пакете, как выглядит вывод:
# Import data
data(nav)
myflows <- prepflows(mat = nav, i = "i", j = "j", fij = "fij")
# Remove the matrix diagonal
diag(myflows) <- 0
# Select flows that represent at least 20% of the sum of outgoing flows for
# each urban area.
flowSel1 <- firstflows(mat = myflows/rowSums(myflows)*100, method = "xfirst",
k = 20)
# Select the dominant flows (incoming flows criterion)
flowSel2 <- domflows(mat = myflows, w = colSums(myflows), k = 1)
# Combine selections
flowSel <- myflows * flowSel1 * flowSel2
# Node weights
inflows <- data.frame(id = colnames(myflows), w = colSums(myflows))
# Plot dominant flows map
opar <- par(mar = c(0,0,2,0))
sp::plot(GE, col = "#cceae7", border = NA)
plotMapDomFlows(mat = flowSel, spdf = UA, spdfid = "ID", w = inflows, wid = "id",
wvar = "w", wcex = 0.05, add = TRUE)
![Output of the code above](https://i.stack.imgur.com/FUOWT.png)
Я хотел бы, чтобы эти линии и многоугольники были на карте листовки с использованием R.
Есть идеи?