緯度、経度を国勢調査地区にジオコードするAPI

米国に緯度と経度の座標（主に都市/町の重心）がたくさんあります。そのために、含まれている国勢調査地区を検索します。

例えば、

input => nn.nnn, -nn.nnn 
output => 010300

役立つ無料または低コストのAPIはありますか？（本質的にはYahoo Placefinder、国勢調査地区の追加の応答を伴う。）

（多くの都市には複数の国勢調査区が含まれていることを理解しています。一部の小さな町には1つ以上の国勢調査区の一部しか含まれていないと思います。今のところその問題を無視しても問題ありません。）

FCC国勢調査ブロック変換 APIは、あなたが探しているものを正確にです。

PostGISでは、ポイントのテーブルがあり、@ Sminbamanguが参照する国勢調査境界情報がある場合、次のアプローチを使用してこれを計算できます。

SELECT 
       c.census_tract, 
       p.point_id
FROM
       census_boundary as c,
       table_of_points as p
WHERE
       ST_Contains(c.geom, p.geom);

ST_Containsの説明はここにあります。これにより、ポイントのテーブルと、それらが配置されている国勢調査地区が返されます。

オンラインAPIの代わりにローカルで行うことができます。Rに関する1つのソリューション：米国の国勢調査境界データは、TIGERのcensus.govから入手できます。米国について考えている場合は、米国の州をダウンロードし、関数を使用して地域を取得できます。たとえば、カリフォルニアとランダム（またはランダムっぽい）ポイントを使用すると、次のようになります。

library(maps)
library(maptools)


tractLookup <- function(x, y, state) {
  pt <- SpatialPoints(data.frame(x = x, y = y))
  overlay.pt <- overlay(pt, state) # what index number does pt fall inside?
  return(census$TRACT[overlay.pt]) # give the Tract number from the census layer
}

やってみて：

california <- readShapePoly("~/Downloads/US_2000_Census_Tracts/tr06_d00_shp/tr06_d00.shp")
tractLookup(-123.123, 40.789, california)

0004を与えます、これは正しいです。

# Look at the map
plot(census)
map('state', c('California'), lwd = 2, col = 'green', add = F) # optional
points(-123.123, 40.789, col = 'red', lwd = 2)

これは機能しますが、私のMacで5秒かかります。これをたくさん行うとしたら、多すぎるでしょう。私は誰かが100万倍速いPostGISソリューションをすぐに検討するのではないかと思います...

@Simbamanguが述べたように、これは市場に出回っているほぼすべてのデスクトップGISソフトウェアによってローカルで簡単に実行できます。

無料のオープンソースについては、QGISの「場所による属性の結合」をご覧ください。

これを繰り返し行う必要がある場合は、データセットをPostGISに入れて空間結合を作成することを検討できます（オンラインヘルプのこの例を参照）。その後、結果をQGISで視覚化できます。