将多边形添加到地图R传单

Question

我正在尝试使用rCharts包显示带有R的地图。 我从简单开始，所以我想向地图添加一个多边形。 但是我不知道如何。 有任何想法吗？ addPolygon不起作用。

map <- Leaflet$new()


map$tileLayer(provider = 'Stamen.TonerLite')

map$setView(c(48.1, 16.7), zoom = 10)
map$addPolygon(
  c(48.99831, 49.08815, 49.08815, 48.99831, 48.99831),
  c(13.42666, 13.42666, 13.56383, 13.56358, 13.42666),
  layerId=c("1"),
  options=opts,
  defaultOptions=opts)
map

Answer 1

通过转换为geoJSON格式将多边形添加到地图中，如rCharts源代码中的示例10： https ： //github.com/ramnathv/rCharts/blob/master/inst/libraries/leaflet/examples/example10.R

请注意geoJSON和setView中的xy坐标之间的纬度和经度有何不同。 此处的代码给了我一个靠近德国的捷克共和国的蓝框。

xy = cbind(
  c(13.42666, 13.42666, 13.56383, 13.56358, 13.42666),
    c(48.99831, 49.08815, 49.08815, 48.99831, 48.99831)
    )

xyjson = RJSONIO::toJSON(xy)

jsonX = paste(
    '{"type":"FeatureCollection","features":[
        {"type":"Feature",
         "properties":{"region_id":1, "region_name":"My Region"},
         "geometry":{"type":"Polygon","coordinates": [ ',xyjson,' ]}}
       ]
      }')

polys = RJSONIO::fromJSON(jsonX)
map = Leaflet$new()
map$tileLayer(provider = 'Stamen.TonerLite')
map$setView(c(49.1,13.5), zoom = 8)
map$geoJson(polys)
map
# or print(map) from a script probably.

如果您有多个多边形，则需要创建{"type": "Feature",多个结构{"type": "Feature",并在"FeatureCollection"的"features"的方括号内用逗号分隔。 为了使结构更好，我重新缩进了一点。 到brew包装等模板系统将帮助您的地步...