Question

我试图使用R“地图”包制作美国各地的热图/等值线/主题图。不幸的是，数据库中似乎缺少一些县。

library(maps)
data(county.fips)
which(county.fips[,1] == 35006)
which(county.fips[,2] == 'new mexico,cibola')

导致未找到任何条目

> library(maps)
> data(county.fips)
> which(county.fips[,1] == 35006)
integer(0)
> which(county.fips[,2] == 'new mexico,cibola')
integer(0)

我知道这个fips代码存在！（http://en.wikipedia.org/wiki/Cibola_County,_New_Mexico http://en.wikipedia.org/wiki/List_of_counties_in_New_Mexico）

任何人都知道如何制作完整的情节？我在这里找到的所有内容都列在这里：http://grokbase.com/t/r/r-sig-geo/0964515wmd/missing-county-in-maps

感谢您的帮助!!

Answer 1

啊哈，问题不是地图，而是country.fips。 Cibola County没有出现在该数据集中，因此当我将我的数据{variable-to-plot，FIPS code}与列表匹配时，Cibola总是丢失。

显示Cibola县的解决方法......为不涉及county.fips数据集的县编写明确的行。

countynames <- match(county.fips[,"fips"], mydatavector$fips)
cibola_value <- mydatavector[which(mydatavector$fips == 35006)]
m <- map("county")
map("county", countynames, col=rainbow(100)[mydatavector$var_to_plot], fill=TRUE, add=TRUE)
map("county", "new mexico,cibola", col=rainbow(100)[cibola_value], fill=TRUE, add=TRUE)

注意：我还发现地图功能无法在美国绘制所有县，但如果我将通话分成两组，则可以使用。

Answer 2

Cibola县在“县”数据库中没有丢失。应该足以绘制地图。例如见：

map.text('county', 'new mexico')

R映射包丢失的县

2 个答案: