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

坐标不是lat-long：

> 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，请使用投影“epsg：4326”的spTransform转换为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