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

【问题标题】：Selecting the range of zip code boundries in R from big data set从大数据集中选择 R 中的邮政编码边界范围
【发布时间】：2015-01-30 19:34:22
【问题描述】：

我找不到某个城市或州的邮政编码边界文件。我从我们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 to integer 时，相应的索引值会随机更改，我会得到随机映射。您知道如何将特定范围的邮政编码数据与大数据区分开来吗？

【问题讨论】：

标签： r plot zipcode

【解决方案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)

【讨论】：

感谢您提供的非常好的示例

【解决方案2】：

您的代码是正确的；您唯一需要考虑的是zip$GEOID10 是因子，因此您需要先将其转换为字符向量：

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

【讨论】：