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

我找不到某个城市或州的邮政编码边界文件。我从美国下载了一个大的邮政编码文件，并试图用以下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)

库（sp）
图书馆（地图工具）
gpclibPermit（）
zip您的代码是正确的；您需要考虑的唯一一点是，
zip$GEOID10
是因子，因此您需要首先将其转换为字符向量：

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

zip$GEOID10您也可以使用
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)

库（sp）
图书馆（地图工具）
库（zipcode）
数据（zipcode）
谢谢你给我一个很好的例子