Построение большого круга ar c (или соединяющих маршрутов) на объекте r maptools - PullRequest
Новая
       105

Построение большого круга ar c (или соединяющих маршрутов) на объекте r maptools

0 голосов
/ 25 марта 2020

Я использую пакет R maptools для визуализации доли сектора обслуживания для выбранных стран, но я также хочу добавить большие дуги окружности (или соединительные маршруты, похожие на маршруты полетов, отображаемые на экране авиакомпаний) между лон-лат столицы каждой страны. Мне интересно, как это сделать в среде maptools, поскольку я делал это только в map?

Ниже приведен мой код, и в качестве примера я привожу три города (показаны красными точками) будет очень полезно, если кто-то продемонстрирует добавление дуг между любыми парами городов на карте. 

rm(list=ls())
library(foreign)
library(maptools)
library(rworldmap)

# load the data
share <- read.csv("https://www.dropbox.com/s/t1hst2z1kfljmu3/share.csv?dl=1")

# rename column
names(share)[2] <- "Percentage GDP"

#create a map-shaped window
mapDevice('x11')
#join to a coarse resolution map
share <- joinCountryData2Map(share, joinCode="NAME", nameJoinColumn="country")
mymap <- mapCountryData(share, nameColumnToPlot="Percentage GDP", catMethod="fixedWidth", addLegend=FALSE, borderCol
 = "grey38")

do.call(addMapLegend, c(mymap
                       ,legendLabels="all"
                       ,legendWidth=0.3
                       ,legendShrink=0.75
                       ,legendIntervals="data"))

# add three cities as example: Bueno Aries, Paris, Melbourne

points(x=-58, y=-34, col="red", pch = 16, cex = 0.1)
points(x=2, y=49, col = "red", pch = 16, cex = 0.1)
points(x=145, y=-38, col = "red", pch = 16, cex = 0.1)
...