Eu estou tentando adicionar um shapefile ESRI (.shp) para um gráfico de ggmap do estado da Carolina do Norte, que tem o seguinte código:

x <- get_map(location="North Carolina", zoom=6, maptype="terrain")

ggmap(x) + geom_polygon(data=citylim83.df, aes(x=long, y=lat), fill="red", alpha=0.2)

O shapefile que eu carreguei e fortifiqueicitylim83.df . Aqui está o código usado para carregar o shapefile no ggplot:

    epsgs <- make_EPSG()
    citylim <- readOGR(dsn=".", layer="MunicipalBoundaries_polys")`

As unidades da projeção dos Limites Municipais, depois de fazer uma pesquisa no EPSG, são ft-US para o sistema de planos de estados. Mesmo que este .shp tenha um sistema de coordenadas geográficas de NAD83, eu também quero projetá-lo para o NAD83 para se livrar do sistema de planos de estado (o código EPSG para NAD83 (UTM-17N) é 26917):

    citylim83 <- spTransform(citylim, CRS("+init=epsg:26917"))
    summary(citylim83)

    Object of class SpatialPolygonsDataFrame
    Is projected: TRUE 
    [+init=epsg:26917 +proj=utm +zone=17 +ellps=GRS80 +datum=NAD83 +units=m

    citylim83.df <- fortify(citylim83)

Este quadro de dados foi então usado no código ggmap mostrado acima.

Mesmo que agora esteja definitivamente projetado no NAD83, ele ainda não aparecerá no ggmap base. Qual é a projeção base do objeto get_map que eu importei? Existe um comando para descobrir isso para que eu possa coincidir com o meu mapa com o shapefile que eu quero exibir em cima dele? Eu tenho que "desprojeto" meu objeto citylim? FYI o shapefile é o limite de cidade limite de cada município da Carolina do Norte, se não fosse claro. Qualquer ajuda seria muito apreciada, pois sou muito nova na comunidade ggplot2 / ggmap.