用R子集geojson数据

Question

我有一个从这里获得的状态边界的geojson文件。特别是我正在使用2000万美国数据

我正在尝试对数据进行子集处理，以便可以使用传单仅映射某些状态。我可以使用以下方法将一个状态子集化：

states <- geojsonio::geojson_read("gz_2010_us_040_00_20m.json", what = "sp")
az <- subset(states, NAME == "Arizona")

但是，此方法似乎不适用于选择多个状态：

swStates <- subset(states, NAME == c("Arizona", "Nevada", "New Mexico"))

选择多个状态通常只会导致选择一个或两个状态，或者根本不会选择一个状态。我可以使用其他方法来对数据进行子集化吗？

Answer 1

您需要使用%in%而不是==

cities <- geojsonio::us_cities
TwoCities <- subset(cities, name %in% c("Seattle WA", "San Francisco CA"))