从GeoJson到DataFrame？

Question

我有一个秘鲁和其州（西班牙语中称为“Departamentos”）的GeoJson文件。

我可以使用leaflet绘制秘鲁的各个州，但由于GeoJson文件没有我需要的所有数据，因此我考虑将其转换为data.frame并添加所需的列数据，然后将其转换回GeoJson格式进行绘图。

数据：您可以从此处下载秘鲁的GeoJson数据：

这是我正在使用的数据，我需要向其添加一列销售额数据，并为每个州（“NOMBDEP”-共24个州）添加一行。

library(leaflet)
library(jsonlite)
library(dplyr)


states <- geojsonio::geojson_read("https://raw.githubusercontent.com/juaneladio/peru-geojson/master/peru_departamental_simple.geojson", what = "sp")

我考虑使用“jsonlite”包将“GeoJson”转换为数据框，但是出现了以下错误：

library(jsonlite)

states <- fromJSON(states)

Error: Argument 'txt' must be a JSON string, URL or file.

我原本期望在拥有数据框之后，我可以做类似下面的事情：

states$sales #sales is a vector with the sales for every department

states <- toJson(states)

- Omar Gonzales

Answer 1

您可以使用library(geojsonsf)将GeoJSON和sf相互转换。

library(geojsonsf)
library(sf)       ## for sf print methods

states <- geojsonsf::geojson_sf("https://raw.githubusercontent.com/juaneladio/peru-geojson/master/peru_departamental_simple.geojson")

states$myNewValue <- 1:nrow(states)

geo <- geojsonsf::sf_geojson(states)

substr(geo, 1, 200)
# [1] "{\"type\":\"FeatureCollection\",\"features\":[{\"type\":\"Feature\",\"properties\":{\"COUNT\":84,\"FIRST_IDDP\":\"01\",\"HECTARES\":3930646.567,\"NOMBDEP\":\"AMAZONAS\",\"myNewValue\":1},\"geometry\":{\"type\":\"Polygon\",\"coordinat"


.

您可以看到myNewValue在GeoJSON中。

- SymbolixAU

Answer 2

您不需要来回转换，只需向states SPDF添加另一列即可：

states <- geojsonio::geojson_read("https://raw.githubusercontent.com/juaneladio/peru-geojson/master/peru_departamental_simple.geojson", what = "sp")
states$sales <- abs(rnorm(nrow(states), sd=1000))
plot(states, col=states$sales)

产生以下图片：

- Stedy