2

我有几千名成员的经纬度列表，并希望附加他们的人口普查区块组代码（即 36001038141345）。我的第一个想法是将这些成员绘制在具有美国人口普查区块组图层的地图上。然后，例如使用其中一个映射包中的叠加函数，或者不使用形状文件将文件转换为数据框，并使用例如“ggplot”和“ggmap”中的连接函数来匹配 lat 和 long 到 14数字人口普查块组值。

## Example of results: ##
GEOID10              LAT        LONG  
1. 1005950300123    +31.7908933 -085.5670514  
2. 01005950900134   +31.8467221 -085.1462332  
3. 01005950800145   +31.9206930 -085.1760317  
4. 01005950700156   +31.9440855 -085.2620842  
5. 01005950600167   +31.8783526 -085.2729215  
6. 01005950100178   +32.0353533 -085.2477678  
7. 01005950400189   +31.6794432 -085.5577196  
8. 01005950200200   +31.9222753 -085.4498309  
9. 01005950500211   +31.7634183 -085.2718649  
10. 01027959000222  +33.3093759 -085.8820883  
11. 01027959200233  +33.1781831 -086.0215486  
12. 01027958900244  +33.3872783 -085.7690615

我正在搜索要在 R 中使用的美国人口普查块组形状文件。对于美国所有州使用美国人口普查网站上的 Tiger 文件以及必要的支持文件（“.shp, .dbf , .prj 等”）使使用这个文件有问题，因为它会导致几百个需要的文件。

还检查了 UScensus2010 包，但它还没有设置“install.blkgrp”功能。我正在寻求一些关于什么是最好的方法的想法或方向。

谢谢

4

1 回答 1

1

我发现这是使用 R 下载人口普查区形状文件的最简单方法：

if (("tigris" %in% rownames(installed.packages()))==FALSE)install.packages("tigris"); library(tigris)
if (("ggplot2" %in% rownames(installed.packages()))==FALSE)install.packages("ggplot2"); library(ggplot2)

census.tracts <- tracts(state="NJ")
plot(census.tracts)
census.tracts.fort <- fortify(census.tracts)

这段代码将为您提供新泽西州的人口普查区形状文件。我知道派对迟到了，但我希望这个过程可以帮助其他人。

于 2016-04-05T18:28:12.510 回答