Próbuję dodać plik kształtu ESRI (.shp) do wykresu ggmap stanu Karolina Północna, który ma następujący kod:

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)

Shapefile, do którego załadowałem i wzmocniłemcitylim83.df . Oto kod używany do załadowania pliku shapefile do ggplot:

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

Jednostki projekcji granic miejskich, po przeprowadzeniu wyszukiwania w EPSG, są ft-US dla systemu płaszczyzny stanu. Mimo że ten plik .shp ma układ współrzędnych geograficznych NAD83, chcę również rzutować go na NAD83, aby pozbyć się systemu płaszczyzny stanu (kod EPSG dla NAD83 (UTM-17N) to 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)

Ta ramka danych została następnie użyta w kodzie ggmap pokazanym powyżej.

Mimo że jest teraz zdecydowanie rzutowany na NAD83, nadal nie pojawi się na podstawowej ggmap. Jaka jest projekcja bazowa importowanego obiektu get_map? Czy istnieje polecenie, aby to sprawdzić, więc mogę dopasować mapę do pliku kształtu, który chcę wyświetlić na wierzchu? Czy muszę „wyrenderować” mój obiekt miejski? FYI shapefile to granica miasta dla każdej gminy w Północnej Karolinie, jeśli nie była jasna. Każda pomoc byłaby bardzo mile widziana, ponieważ jestem bardzo nowa w społeczności ggplot2 / ggmap.