python - 如何从 Census FIPS 转换为 Lat Lon？

Question

我正在尝试通过绘制/操作公共数据集来学习一些编程（python）技能。目前，我正在尝试制作北极的交互式人口地图。我从阿拉斯加州劳工部下载了阿拉斯加人口普查数据，该数据给了我按 FIPS 中的人口普查区划分的人口（即 0201301090 是阿留申斯东区的阿库坦市）。有没有办法将 FIPS 转换为 Lat Lon 或在 python 中通过 FIPS 绘制？

Answer 1

您示例中的地点的 FIPS 代码（01090 部分）采用人口普查局很久以前停用的样式（FIPS 55.3），因此找到一个好的数据源有点棘手。我从人口普查局尝试了一堆花哨的 GIS 库和大量大型 shp 文件，但尽管您可以非常接近人口普查块和区域，但它们都没有足够的数据来追踪 FIPS 的位置。

但是，我最终发现美国地理名称委员会在这里有一个有用的数据库，它仍然使用这种格式。

如果您从该页面底部的“带有联邦代码的州文件”部分下载阿拉斯加文件，则只需使用 python 的内置 csv 阅读器即可轻松解析必要的信息。

import csv
datapath = r"C:\folder\where\you\saved\the\file\AK_FedCodes_20170201.txt"
data = []

with open(datapath) as f:
    reader = csv.reader(f, delimiter="|")
    rows = [row for row in reader]

    for row in rows[1:]:

        name = row[1]
        fips = row[7]+row[10]+row[3]
        lat = row[12]
        lon = row[13]

        data.append({"name": name,
                     "fips": fips,
                     "lat": lat,
                     "lon": lon})

        if name == "Akutan":
            print name, fips, lat, lon # just because we can

输出：

Akutan 0201301090 54.1355556 -165.7730556