Question

我有几千名成员的纬度和长度列表，并希望附加他们的人口普查区块组代码（即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'功能。我正在寻求一些想法或方向，以确定最佳方法。

谢谢

Answer 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)

这段代码将为您提供新泽西州的人口普查区形状文件。我知道这个晚会，但我希望这个过程可以帮助别人。

人口普查区块组形状文件

1 个答案: