После пересечения двух шейп-файлов: "Ошибка в vapply(g2, st_is_empty, logic (1))"
Я запускаю следующее:
library(dplyr)
library(sf)
library(tigris)
library(tmap)
options(tigris_class = 'sf')
options(tigris_use_cache = TRUE)
nj = tigris::states(cb = T, year = 2015) %>%
filter(STUSPS == 'NJ')
nj_msas = tigris::core_based_statistical_areas(cb = T, year = 2015) %>%
filter(grepl('NJ', NAME)) %>%
sf::st_intersection(nj)
tmap_mode('plot')
nj_msas %>%
tm_shape() +
tm_polygons()
Последний кусок дает ошибку:
Ошибка в vapply(g2, st_is_empty, logic (1)): значения должны быть длиной 1, а результат FUN(X[[3]]) - длиной 148
Когда я удаляю st_intersectionЯ не получаю ошибки от последнего куска. Я не могу найти это сообщение об ошибке нигде в Google. Кто-нибудь знает, что происходит?
Кроме того, если я запускаю все вышеперечисленное, кроме последнего чанка, и я использую ggplot2::geom_sf создать карту, а не tmap функции, я получаю карту, которую я хочу. Нет ошибок
Я использую Ubuntu 18.04.1. RStudio v1.1.463. R 3.5.2. Тигр 0,7. Tmap 2.2. СФ 0,7-2.
2 ответа
Кажется, из-за того, что столбец геометрии, полученный в результате пересечения, имеет тип "ГЕОМЕТРИЯ", а не тип "ПОЛИГОН", вероятно, из-за незначительных смещений между tigris::states а также tigris::core_based_statistical_areasчто приводит к ошибке в tm_shape():
> nj_msas
Simple feature collection with 7 features and 17 fields
geometry type: GEOMETRY
dimension: XY
bbox: xmin: -75.55961 ymin: 38.92852 xmax: -73.89398 ymax: 41.35742
epsg (SRID): 4269
proj4string: +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
CSAFP CBSAFP AFFGEOID GEOID NAME LSAD ALAND AWATER STATEFP
1 428 36140 310M200US36140 36140 Ocean City, NJ M1 651209688 955666476 34
2 408 10900 310M200US10900 10900 Allentown-Bethlehem-Easton, PA-NJ M1 3763879943 58581593 34
3 428 47220 310M200US47220 47220 Vineland-Bridgeton, NJ M1 1252937239 502091171 34
4 408 45940 310M200US45940 45940 Trenton, NJ M1 581629671 11189320 34
5 428 37980 310M200US37980 37980 Philadelphia-Camden-Wilmington, PA-NJ-DE-MD M1 11920145588 693341961 34
6 408 35620 310M200US35620 35620 New York-Newark-Jersey City, NY-NJ-PA M1 21478408908 6689374328 34
7 428 12100 310M200US12100 12100 Atlantic City-Hammonton, NJ M1 1439256827 300767732 34
STATENS AFFGEOID.1 GEOID.1 STUSPS NAME.1 LSAD.1 ALAND.1 AWATER.1 geometry
1 01779795 0400000US34 34 NJ New Jersey 00 19048075783 3543447118 MULTIPOLYGON (((-74.55255 3...
2 01779795 0400000US34 34 NJ New Jersey 00 19048075783 3543447118 POLYGON ((-75.12051 40.9683...
3 01779795 0400000US34 34 NJ New Jersey 00 19048075783 3543447118 POLYGON ((-75.41956 39.4132...
4 01779795 0400000US34 34 NJ New Jersey 00 19048075783 3543447118 POLYGON ((-74.94228 40.3408...
5 01779795 0400000US34 34 NJ New Jersey 00 19048075783 3543447118 GEOMETRYCOLLECTION (LINESTR...
6 01779795 0400000US34 34 NJ New Jersey 00 19048075783 3543447118 POLYGON ((-74.02454 40.7094...
7 01779795 0400000US34 34 NJ New Jersey 00 19048075783 3543447118 POLYGON ((-74.98522 39.5148...
Вы можете решить эту проблему, извлекая только полигоны из результата пересечения:
library(dplyr)
library(sf)
library(tigris)
library(tmap)
options(tigris_class = 'sf')
options(tigris_use_cache = TRUE)
nj = tigris::states(cb = T, year = 2015) %>%
filter(STUSPS == 'NJ')
nj_msas = tigris::core_based_statistical_areas(cb = T, year = 2015) %>%
filter(grepl('NJ', NAME)) %>%
sf::st_intersection(nj) %>%
sf::st_collection_extract("POLYGON")
tmap_mode('plot')
#> tmap mode set to plotting
tm_shape(nj_msas) + tm_polygons()
или установив точность на nj Набор данных:
nj_msas = tigris::core_based_statistical_areas(cb = T, year = 2015) %>%
filter(grepl('NJ', NAME)) %>%
sf::st_intersection(sf::st_set_precision(nj ,1000))
tm_shape(nj_msas) + tm_polygons()
Создано в 2019-01-10 пакетом представлением (v0.2.1)
Ответ Ibusett работал для меня, добавив
%>% sf::st_collection_extract("POLYGON")
до конца моей трансформации изменить, пытаясь определить границы школьного округа только в пределах границ округа.
Это также все еще работает с ggplot.
Спасибо