距離行列の地理的なポイント（首都）

Question

私は世界中の国々の首都の距離行列を作成しようとしています。私はこのような国の重心の距離行列を作成しました。シェープファイルはhttp://www.gadm.org/version2から派生さ：

library(rgeos) 
library(rgdal) 
shapefile <- readOGR("./Map/World Map", layer = "TM_WORLD_BORDERS-0.3") 

centroids <- gCentroid(shapefile, byid = TRUE, as.character(shapefile@data$UN) 
dist_matrix <- as.data.frame(geosphere::distm(centroids) 
colnames(dist_matrix) <- shapefile@data$NAME 
rownames(dist_matrix) <- shapefile@data$NAME 

今、私は同じことをしたいが、国の首都の代わりに、国の重心と。理想的には、資本ではない他の地理的なポイントがあっても機能するメソッドが必要です。これまで私は与えられたshapefileに大文字をプロットすることにしましたが、それらの中から距離行列を作成する方法を見つけることはできません。

library(rgeos) 
library(rgdal) 
library(maps) 
library(tidyverse) 
shapefile <- readOGR("./Map/World Map", layer = "TM_WORLD_BORDERS-0.3") 

data("world.cities) 
world..cities <- world.cities %>% 
filter(world.cities$capital == 1) 
plot(shapefile) 
points(world.cities$long, world.cities$lat, col ="red", cex = .6, pch = 22, add = TRUE) 

理想的な出力形式は、xがメートルまたはキロメートルで、それぞれの国の首都の間の距離を表し、この

  Algeria Albania Azerbaijan 
Algeria  0  x   x 
Albania  x  0   x 
Azerbaijan x  x   0 

ようになります。

Answer 1

あなたが必要とするのは、緯度/経度で与えられたポイント間の距離を計算する関数だと思います。これは地球圏パッケージで利用可能です。

library(geosphere) 
Caps = cbind(world.cities$long, world.cities$lat) 
CapDistMatrix = distGeo(Caps, Caps)