Как скачать данные геойсона и прочитать их в R
# From http://leafletjs.com/examples/choropleth/us-states.js
states <- geojsonio::geojson_read("json/us-states.geojson", what = "sp")
bins <- c(0, 10, 20, 50, 100, 200, 500, 1000, Inf)
pal <- colorBin("YlOrRd", domain = states$density, bins = bins)
labels <- sprintf(
"<strong>%s</strong><br/>%g people / mi<sup>2</sup>",
states$name, states$density
) %>% lapply(htmltools::HTML)
leaflet(states) %>%
setView(-96, 37.8, 4) %>%
addProviderTiles("MapBox", options = providerTileOptions(
id = "mapbox.light",
accessToken = Sys.getenv('MAPBOX_ACCESS_TOKEN'))) %>%
addPolygons(
fillColor = ~pal(density),
weight = 2,
opacity = 1,
color = "white",
dashArray = "3",
fillOpacity = 0.7,
highlight = highlightOptions(
weight = 5,
color = "#666",
dashArray = "",
fillOpacity = 0.7,
bringToFront = TRUE),
label = labels,
labelOptions = labelOptions(
style = list("font-weight" = "normal", padding = "3px 8px"),
textsize = "15px",
direction = "auto")) %>%
addLegend(pal = pal, values = ~density, opacity = 0.7, title = NULL,
position = "bottomright")
Приведенный выше код скопирован с https://rstudio.github.io/leaflet/choropleths.html.
Я пытаюсь воспроизвести вывод. Однако я застрял на первом этапе - загрузке файла геойсона. Я использовал ссылку, показанную в первой строке, и сохранил ее в виде текстового файла, а затем переименовал в файл геоджон. Но я не смог прочитать этот файл. Очевидно, что-то не так с загрузкой файлов или загрузкой в R, но я понятия не имею, где это.
Может кто-нибудь дать какие-либо инструкции? Я никогда не имел дело с данными геойсон раньше. Мне просто нужна помощь с первыми двумя строками кода, и я могу справиться со всеми остальными самостоятельно.
1 ответ
Решение
Загружаемый файл имеет назначение javascript во главе. Удаление, кажется, решает проблему,
library(geojson)
library(geojsonio)
url <- "http://leafletjs.com/examples/choropleth/us-states.js"
# read as text file
doc <- readLines(url)
# remove the javascript assignment at the front
doc2 <- gsub("var statesData = ", "", doc)
# write out as a temp file and read
write(doc2, file = "tempgeo.json")
states <- geojson_read("tempgeo.json", what = "sp")