Converti le coordinate da readShapePoly in R in coordinate long-lat

Sto cercando di trovare alcuni centroidi di uno shapefile per i comuni danesi e successivamente trovare il tempo di guida tra di loro. Uso Rla readShapePolyfunzione di maptoolscombinata con la gCentroidfunzione di rgeose tutto funziona. Tuttavia, ottengo punti spaziali come

SpatialPoints:
     x       y
1 571860.7 6225016
Coordinate Reference System (CRS) arguments: NA 

Che chiaramente non è qualcosa che posso usare in Google per prendere i tempi di viaggio. Sto cercando un modo per convertire questi numeri in longitudine-latitudine, ma non ho idea di come.

Quando leggo i dati utilizzando readOGRdalla rgdallibreria ottengo le stesse coordinate ma mi dice quanto segue su ciò che presumo sia la proiezione (ma le coordinate sono le stesse)

Slot "proj4string":
CRS arguments:
+proj=utm +zone=32 +ellps=intl +units=m +no_defs

Esempio riproducibile: ho inserito i dati per l'esempio qui: https://github.com/sebastianbarfort/shapefiles

Questo dovrebbe riprodurre il problema:

library(maptools)
library(rgdal)
library(rgeos)

map = readShapePoly("~/Downloads/shapefiles-master/kommuner1983.shp")
centroid = gCentroid(map)
centroid

Utilizzare spTransformper trasformare le coordinate in WGS84:

library("rgdal")
library("rgeos")

map <- readOGR(".", "kommuner1983")
map_wgs84 <- spTransform(map, CRS("+proj=longlat +datum=WGS84"))
plot(map_wgs84, axes=TRUE)

gCentroid(map_wgs84)
# SpatialPoints:
#       x     y
# 1 10.05 55.96
# Coordinate Reference System (CRS) arguments: +proj=longlat +datum=WGS84
# +ellps=WGS84 +towgs84=0,0,0 

rgdal::readOGRè in grado di leggere automaticamente le informazioni sulla proiezione. maptoolsle funzioni non leggono né scrivono informazioni sulla proiezione, lasciando a te la gestione manuale di questi dettagli.