Converter coordenadas de readShapePoly em R em coordenadas long-lat

Estou tentando encontrar alguns centróides de um shapefile para os municípios dinamarqueses e depois encontrar o tempo de condução entre eles. Eu uso Ra readShapePolyfunção de maptoolscombinada com a gCentroidfunção de rgeose tudo funciona. No entanto, recebo pontos espaciais como

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

O que claramente não é algo que eu possa usar no Google para aproveitar os tempos de viagem. Estou procurando uma maneira de converter esses números em longitude-latitude, mas não tenho idéia de como.

Quando leio os dados usando readOGRda rgdalbiblioteca, obtenho as mesmas coordenadas, mas ele me diz o seguinte sobre o que presumo ser a projeção (mas as coordenadas são as mesmas)

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

Exemplo reproduzível: coloquei os dados para o exemplo aqui: https://github.com/sebastianbarfort/shapefiles

Isso deve reproduzir o problema:

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

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

Use spTransformpara transformar as coordenadas em 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é capaz de ler as informações da projeção automaticamente. maptoolsAs funções não leem nem gravam informações de projeção, cabendo a você gerenciar esses detalhes manualmente.