Я пытаюсь найти несколько центроидов шейп-файла для датских муниципалитетов, а затем найти время в пути между ними. Я использую R
«s readShapePoly
функцию из в maptools
сочетании с gCentroid
функцией от rgeos
, и все работает. Тем не менее, я получаю пространственные точки, такие как
SpatialPoints:
x y
1 571860.7 6225016
Coordinate Reference System (CRS) arguments: NA
Это явно не то, что я могу использовать в Google, чтобы захватить время в пути. Я ищу способ преобразования этих чисел в долготу-широту, но понятия не имею, как.
Когда я читаю данные, используя readOGR
из rgdal
библиотеки, я получаю те же координаты, но они говорят мне следующее о том, что я предполагаю, является проекцией (но координаты одинаковы)
Slot "proj4string":
CRS arguments:
+proj=utm +zone=32 +ellps=intl +units=m +no_defs
Воспроизводимый пример: я поместил данные для примера здесь: https://github.com/sebastianbarfort/shapefiles
Это должно воспроизвести проблему:
library(maptools)
library(rgdal)
library(rgeos)
map = readShapePoly("~/Downloads/shapefiles-master/kommuner1983.shp")
centroid = gCentroid(map)
centroid
источник