经度和纬度在ggplot2中无法正确显示

Question

我现在正试图绘制加拿大的地图,ggplot2我发现经度和纬度没有在地图中正确显示.有什么解决方案吗？非常感谢.

arcgis shapfile是从https://www.arcgis.com/home/item.html?id=dcbcdf86939548af81efbd2d732336db下载的

library(ggplot2)
library(rgdal)
countries<-readOGR("Canada.shp", layer="Canada")
ggplot()+geom_polygon(data=countries,aes(x=long,y=lat,group=group),fill='white',color = "black")

地图中的经度应该是110W,100W,90W.并且地图中的纬度应该像50N,60N,70N.但是,现在情况并非如此.

Answer 1

坐标不是纬度:

> summary(countries)
Object of class SpatialPolygonsDataFrame
Coordinates:
         min     max
x -2314694.5 3093025
y   321591.9 4811137
Is projected: TRUE 
proj4string :
[+proj=aea +lat_1=50 +lat_2=70 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0
+datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0]

它们是"aea",它是具有给定参数的Albers Equal Area.

要转换为lat-long,请使用spTransform投影"epsg:4326"转换为WGS84 lat-long,如GPS系统中所使用的那样.

> ca = spTransform(countries, "+init=epsg:4326")
> summary(ca)
Object of class SpatialPolygonsDataFrame
Coordinates:
         min       max
x -141.00301 -52.62028
y   41.91332  83.10832