我正在尝试使用r中的纬度和经度找到矩阵距离。我的数据框架有三列,第一列是地区(城市)的名称,第二列是经度,第三列是纬度。如果我只选择经度和纬度列,我可以使用下面的代码获得矩阵距离,但我想要的是在矩阵上也有城市的名称。请注意,下面代码中的数据只有两列,因为我不能让下面的代码在城市列上也工作。数据集中有1152行。非常感谢您的帮助。
dist_mat <- distm(data, fun = distGeo)
dist_mat
我制作了一个包含10个城市的数据集d
。此代码也可以扩展到1152个城市。
d <- data.frame(cities = letters[1:10], longitude = rnorm(10), latitude = rnorm(10))
library(geosphere)
dist_mat <- distm(d[, 2:3], fun = distGeo)
rownames(dist_mat) <- d$cities
colnames(dist_mat) <- d$cities
您只需将dist_mat
的rownames
和colnames
设置为您的citie
名称。