如何使用R在传单上绘制多边形？

Question

在尝试使用R在Leaflet顶部绘制SpatialPolygonsDataFrame时遇到一个简单的问题。我的代码如下：

leaflet() %>%
    addProviderTiles("CartoDB.Positron") %>%
    setView(lng = -80.8858673, lat = 41.1450276, zoom = 5) %>%
    addPolygons(data = SPDF, weight = 2, color = ~colorQuantile("red", SPDF$id)(id))

其中SPDF是我的SpatialPolygonsDataFrame。

当我执行此代码时，它“ PLOTS NOTHING”，但只有底图。 我一直在搜索，这个问题很相似，但是没有这个问题。

为了绘制多边形，我一直关注此链接。

这个问题似乎很简单，但是却为我浪费了时间。 期待提出建议。 感谢您的时间。

注意：SPDF包含从OSM导出的数据，这意味着（多边形）的坐标没有OSM数据中的小数点。

Answer 1

最后，我自己能够解决问题。 问题出在投影和CRS（坐标参考系统）中 。

默认的proj4string设置不正确，因此坐标不真实（无小数点）。 因此，首先，我设置了SpatialPolygonsDataFrame（SPDF）的默认proj4string：

SPDF@proj4string <-CRS("+init=epsg:3857")

设置此设置后，我提供了如下投影：

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

现在，当我执行以下代码传单行时，它可以完美运行。

leaflet() %>%
    addProviderTiles("CartoDB.Positron") %>%
    setView(lng = -80.8858673, lat = 41.1450276, zoom = 5) %>%
    addPolygons(data = SPDF, weight = 2, color = ~colorQuantile("red", SPDF$osm_id)(osm_id))

为了解决这个问题我也跟着上讨论此页。

我希望这对面临相同问题的其他人也有用。 尽管我仍然不是投影和制图专家，但是，如果有人可以推荐一些必要的信息来理解此类问题，那将是很好的。 感谢所有的时间。

PS：请确保已包含必要的软件包，例如传单，sp，magrittr等。