Как сохранить идентификаторы с пакетом geojsonio?

Question

Я пытаюсь прочитать файл topoJSON, созданный с geojsonio пакет: https://raw.githubusercontent.com/pachamaltese/chilemaps/master/data-raw/r15_arica_y_parinacota.json

Вот предварительный просмотр соответствующей части первой записи в этом файле:

{"type":"Topology", ... , "id":295,"properties":{"comuna":"Arica","id":295}} ...

Для импорта я запускаю:

r15 <- "https://raw.githubusercontent.com/pachamaltese/chilemaps/master/data-raw/r15_arica_y_parinacota.json"
r15 <- geojsonio::topojson_read(r15)

И затем, чтобы преобразовать обратно в topoJSON, я запускаю:

geojsonio::topojson_json(r15)

Который возвращает:

{"type":"Topology", ... , "id":0,"properties":{"id":"295","comuna":"Arica"}} ...

Быть тем id=0 неправильно, если это не я, как и должно быть id=295 из JSON я читаю.

Как я могу быть уверен, что я правильно читаю и пишу в JSON?

Answer 1

Спасибо @sckott

Ваш пример дал мне простую идею, которая решает это. Будучи моим первоначальным намерением визуализировать карту, вот полностью воспроизводимый пример того, что сработало. Я решил создать новый id столбец в таблице, который я буду использовать с данными topoJSON:

if (!require("pacman")) { install.packages("pacman") }
pacman::p_load(geojsonio, dplyr)
pacman::p_load_gh("pachamaltese/d3plus", "pachamaltese/chilemaps")

# data to visualize on a map

data3 <- tibble(
  id = c(307, 295, 302, 331),
  comuna = c("General Lagos", "Arica", "Camarones", "Putre"),
  valor = c(400,300,200,100)
) %>%
  arrange(id) %>%
  mutate(id2 = row_number() - 1) # trick to match modified JSON "ids" after reading from GH

# visualize using SpatialPolygonsDataFrame

r15 <-  subset(chilemaps::comunas, region_id == 15)

d3plus() %>%
  d3p_data(data3) %>%
  d3p_map(coords = geojsonio::topojson_json(r15), text = "comuna", tooltip = "valor") %>%
  d3p_id(c("id")) %>%
  d3p_colour("valor")

# visualize using topoJSON read from GitHub

r15_2 <- "https://raw.githubusercontent.com/pachamaltese/chilemaps/master/data-json/r15_arica_y_parinacota.json"
r15_2_out <- geojsonio::topojson_read(r15_2)

d3plus() %>%
  d3p_data(data3) %>%
  d3p_map(coords = geojsonio::topojson_json(r15_2_out), text = "comuna", tooltip = "valor") %>%
  d3p_id(c("id2")) %>% # here I use id2 instead of id
  d3p_colour("valor")

Answer 2

Идентификаторы становятся искаженными sp пкг я думаю

library(geojsonio)
r15 <- "https://raw.githubusercontent.com/pachamaltese/chilemaps/master/data-raw/r15_arica_y_parinacota.json"
out <- topojson_read(r15)
x <- topojson_list(out)
vapply(x$objects$foo$geometries, "[[", 1, "id")
#> [1] 0 1 2 3
vapply(x$objects$foo$geometries, function(z) z$properties$id, "")
#> [1] "295" "302" "307" "331"

Мы используем SP внутри topojson_read,

Таким образом, "295", "302", "307", "331" вкладываются внутрь других идентификаторов.

Если мы делаем

s <- sf::st_read(r15)
topojson_list(zz)

Я думаю, что получает то, что вы хотите