r - 如何使用传单绘制 SpatialPolygonsDataframe 中的特定列？

Question

我想使用 leafletR::leaflet 创建分区统计图。我的数据位于 SpatialPolygonsDataFrame 中，我想选择要绘制的特定列。

使用 sp::spplot，这很容易，因为参数 zcol 允许我指定要绘制的图层/列:

library("maptools");library("sp");library("leafletR")
SP <- readShapePoly(system.file("shapes/sids.shp",  
                    package="maptools")[1],
                    proj4string=CRS("+proj=longlat +datum=WGS84 
                                     +no_defs +ellps=WGS84 
                                     +towgs84=0,0,0"))
spplot(SP, zcol="BIR79")

但是，使用leafletR，我不知道如何指定图层，它只是绘制普通 map 边框:

SP4leaflet <- toGeoJSON(data=SP, dest=tempdir(), name="BIR79")
SPleaflet  <- leaflet(data=SP4leaflet, dest=tempdir(), 
                  title="Trying to plot BIR79",
                  base.map="osm", popup="*")
SPleaflet

关于如何选择要使用 leafletR 绘制的所需图层/列有什么想法吗？

Answer 1

为了完整起见，我想提一下，现在可以使用library(mapview)轻松实现这一点。

library("maptools")
library("sp")
library("mapview")

SP <- readShapePoly(system.file("shapes/sids.shp",  
                package="maptools")[1],
                proj4string=CRS("+proj=longlat +datum=WGS84 
                                 +no_defs +ellps=WGS84 
                                 +towgs84=0,0,0"))

spplot(SP, zcol="BIR79")
mapview(SP, zcol="BIR79")

另请参阅http://edzer.github.io/sp/#interactive-maps-leaflet-mapview对于类似的例子和 http://environmentalinformatics-marburg.github.io/web-presentations/20150723_mapView.html了解mapview的更多功能。