r - 使用 R 在传单 map 上投影我的 shapefile 数据

Question

我(在 Win7 上)试图使用 leaflet 显示我的 shapefile 数据(Here 是数据文件。)包裹。但没有任何成功。我只得到背景 map 图块，而不是我的数据。不过，我可以在 QGIS 中看到我的数据。我从 here 得到以下代码和描述.

library(rgdal)
shapeData <- readOGR(".",'myGIS')
ogrInfo(".",'myGIS')
leaflet()  %>% addTiles() %>% setView(lng = -106.363590, lat=31.968483,zoom=11) %>% addPolygons(data=shapeData,weight=5,col = 'red') %>% addMarkers(lng = -106.363590,lat=31.968483,popup="Hi there") 

这是我在互联网浏览器中看到的。我没有看到 shapeData 被投影到它上面:

Answer 1

Victorp 的回答有效。但是我建议使用:

shapeData <- spTransform(shapeData, CRS("+proj=longlat +datum=WGS84 +no_defs"))

此 CRS 对应于 EPSG:4326。
Leaflet 实际上负责将 EPSG:4326 转换为 EPSG:3857(称为“Google Mercator”)。

现在，差异(GRS80 与 WGS84)可以忽略不计(在一个轴上为 0.1 毫米)。似乎传单认为两者都是同一个椭圆体。但是，出于测试目的，建议使用 EPSG:4326。