我正在尝试读取geojson数据以与leaflet
包一起使用。我尝试以多种方式读取数据,但无法绘制它。数据为here。
我尝试使用this StackOverflow post中的代码。第一块给了我一张空白的地图。第二个给了我一张没有添加新行的地图。
我还尝试了here中的代码。我对此并不陌生,因此不胜感激。
答案 0 :(得分:2)
这是您想要的吗?
library(leaflet)
library(jsonlite)
download.file("https://opendata.arcgis.com/datasets/772f3621fa354ec9abf3ba33f3ace59e_0.geojson", "RPD_Sections.geojson")
x = fromJSON("RPD_Sections.geojson", FALSE)
leaflet() %>% addTiles() %>% addGeoJSON(x) %>% setView(
lng = mean(vapply(x$features[[1]]$geometry$coordinates[[1]], "[[", 1, 1)),
lat = mean(vapply(x$features[[1]]$geometry$coordinates[[1]], "[[", 1, 2)),
zoom = 12)
答案 1 :(得分:0)
简洁一些,但结果相同
library(sf)
library(leaflet)
sf <- sf::st_read("https://opendata.arcgis.com/datasets/772f3621fa354ec9abf3ba33f3ace59e_0.geojson")
leaflet() %>%
addTiles() %>%
addPolygons(data = sf)
这也假设您的GeoJSON仅包含多边形