如何在 R 中过滤大型 geojson 文件

Question

我在使用 R 过滤大型 geojson 文件时遇到问题。如果我只想显示一个国家/地区的地图，我不想加载整个欧洲地图，因为它很大。所以我想过滤这个数据集，例如保加利亚 -- CNTR_CODE == "BG" 但我无法管理。请在下面找到要下载的代码以及未产生预期结果的初步努力

link <- 'https://ec.europa.eu/eurostat/cache/GISCO/distribution/v2/nuts/download/ref-nuts-2013-01m.geojson.zip'
temp <- tempfile()
download.file(link,temp)
mapdata <- readLines(unzip(temp, "NUTS_RG_01M_2013_4326_LEVL_3.geojson"))
mapdata <- jsonlite::fromJSON(mapdata, simplifyVector = FALSE)
#glimpse(mapdata)
mapdata$features[[100]]$properties$CNTR_CODE
[1] "BG"


  library(sf)
  mapdata2 <- copy(mapdata)
  mapdata2 %>% 
    filter(CNTR_CODE %in% c('BG'))

谢谢。

Answer 1

如果你想要一个列表，那么 Filter 可以：

path = "NUTS_RG_01M_2013_4326_LEVL_3.geojson"
x <- jsonlite::fromJSON(path, simplifyVector = FALSE)
x$features <- Filter(function(z) z$properties$CNTR_CODE == "BG", x$features)
vapply(x$features, function(x) x$properties$CNTR_CODE, "")

如果你想以 geojson 字符格式保存数据，你可以使用 jqr

path = "NUTS_RG_01M_2013_4326_LEVL_3.geojson"
x <- paste0(readLines(path), collapse = "")
xx <- jqr::jq(x, '.features |= map(select(.properties.CNTR_CODE == "BG"))')
jqr::jq(xx, '.features[].properties.CNTR_CODE')

Answer 2

如果使用正确的工具，这个文件不会很大。

library(geojsonsf) 可以直接将 geojson 读取到 sf 对象
library(mapdeck)可以画出所有的多边形

link <- 'https://ec.europa.eu/eurostat/cache/GISCO/distribution/v2/nuts/download/ref-nuts-2013-01m.geojson.zip'
temp <- tempfile()
download.file(link,temp)

library(geojsonsf)
library(sf)

sf <- geojsonsf::geojson_sf(unzip(temp, "NUTS_RG_01M_2013_4326_LEVL_3.geojson"))

然后您可以过滤 sf 对象并绘制

library(mapdeck)
set_token( "YOUR_MAPBOX_API_TOKEN" )

mapdeck(
    style = mapdeck_style("dark")
) %>%
    add_polygon(
        data = sf[ sf$CNTR_CODE %in% c("BG"), ]
        , fill_colour = "NUTS_NAME"
        , legend = T
    )

或者绘制整个地段

mapdeck(
    style = mapdeck_style("dark")
) %>%
    add_polygon(
        data = sf
        , fill_colour = "NUTS_NAME"
        , legend = T
    )

如何在 R 中过滤大型 geojson 文件

How to filter a large geojson file in R

json

r

geojson

r-highcharter