我正在使用人口普查数据中的邮政编码信息以及我对每个邮政编码的数据创建地图。
我从https://www.census.gov/geo/maps-data/data/tiger-cart-boundary.html获取了邮政编码的shapefile。
我使用readOGR使用以下代码成功读取了shapefile:
Zip = readOGR("C:\\...pathname....\\cb_2013_us_zcta510_500k.shp")
当我尝试使用fortify时遇到问题,以便我可以将shapefile中的信息与我的数据结合起来。
names(Zip)
[1] "ZCTA5CE10" "AFFGEOID10" "GEOID10" "ALAND10" "AWATER10"
test = ggplot2::fortify(Zip, region="AWATER10")
Error in fortify.SpatialPolygonsDataFrame(Zip, region = "AWATER10") :
`region' not supported
当我尝试从名称(Zip)中获取所有结果时,会复制错误。我不确定哪个列包含我感兴趣的邮政编码信息,所以我只是尝试了所有这些。
我的最终目标是使用ggmap的邮政编码信息根据我的数据中的信息填写不同颜色的邮政编码。