r - 从大数据集中选择R中邮政编码边界的范围

Question

我找不到特定城市或州的邮政编码边界文件。我从我们这里下载了一个大的邮政编码文件 census 并尝试使用以下 R 代码分隔佛罗里达州的邮政编码边界:

library(sp)
library(maptools)

gpclibPermit()

zip <- readShapeSpatial("~/data/cb_2013_us_zcta510_500k.shp")

zip$GEOID10 <- as.integer(zip$GEOID10)

bool_fl <- 32000 < zip$GEOID10 & zip$GEOID10 <34999

zip1 <-subset(zip, bool_fl)

plot(zip1)

我遇到的问题是，当我将zip$GEOID10更改为整数时，相应的索引值会随机更改，并且我会得到随机 map 。您知道如何将特定范围的邮政编码数据与大数据分开吗？

Answer 1

您还可以使用 zipcode 包来提供帮助:

library(sp)
library(maptools)
library(zipcode)

data(zipcode)

zip <- readShapeSpatial("cb_2013_us_zcta510_500k/cb_2013_us_zcta510_500k.shp")
fl <- zip[zip$GEOID10 %in% zipcode[zipcode$state=="FL",]$zip,]
plot(fl)