r - 使用 Leaflet(对于 R)加载 SpatialPolygonsDataFrame 不起作用

Question

首先，我是 R 的新手，所以请耐心等待。

我最终想要完成的是用 Leaflet 显示阿姆斯特丹的交互式 map 。为此，我使用 RGDAL 来读取 shapefile。

此 link包含阿姆斯特丹的 shapefile。

我正在使用以下代码读取 shapefile 并显示 map 。

amsterdam <- readOGR(".", layer = "sd2010zw_region", verbose = FALSE)

leaflet(amsterdam) %>%
  addProviderTiles("CartoDB.Positron", options= providerTileOptions(opacity = 0.99)) %>%
    addPolygons(
      stroke = FALSE, fillOpacity = 0.5, smoothFactor = 0.5
    )

我得到的是来自 CartoDB.Positron 的 map ，而不是作为第二层的“多边形 map ”。我得到的是一个包含各种数据的 SpatialPolygonsDataFrame。

另一方面，当我使用 plot 方法时，我得到了阿姆斯特丹的 map

plot(amsterdam, axes=TRUE, border="gray")

但我不想使用情节，我想使用传单:)

我在这里做错了什么？

Answer 1

问题是投影。您需要将您的数据投影到 longlat使用 spTransform来自 rgdal 或 sp .另外，请提供您的 SpatialPolygonsDataFrame在 addPolygons()称呼。

library(leaflet)
library(rgdal)

amsterdam <- readOGR(".", layer = "sd2010zw_region", verbose = FALSE)

ams_ll <- spTransform(amsterdam, CRS("+init=epsg:4326"))

leaflet() %>%
  addProviderTiles("CartoDB.Positron", options= providerTileOptions(opacity = 0.99)) %>%
  addPolygons(data = ams_ll,
    stroke = FALSE, fillOpacity = 0.5, smoothFactor = 0.5
  )