API để đảo ngược vĩ độ mã địa lý, kinh độ cho đường điều tra dân số?

Tôi có một loạt các tọa độ kinh độ, vĩ độ (chủ yếu là thành phố / thị trấn) ở Hoa Kỳ mà tôi muốn tìm đường điều tra dân số có chứa.

Ví dụ,

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

Có một API miễn phí hoặc chi phí thấp ngoài kia có thể giúp đỡ không? (Về cơ bản là Yahoo Placefinder, với phản hồi bổ sung của điều tra dân số.)

(Tôi hiểu rằng nhiều thành phố bao gồm nhiều vùng điều tra dân số và tôi hy vọng rằng một số thị trấn nhỏ chỉ chứa một phần của một hoặc nhiều vùng điều tra dân số. Tôi hiện đang thoải mái bỏ qua vấn đề đó.)

Các FCC khối điều tra dân số chuyển đổi API là chính xác những gì bạn đang tìm kiếm.

Trong PostGIS, nếu bạn có bảng điểm và thông tin ranh giới điều tra dân số mà @Sminbamangu đề cập đến bạn có thể tính toán điều này bằng cách sử dụng phương pháp sau:

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

Bạn có thể xem mô tả ST_Contains tại đây . Điều này sẽ trả về một bảng các điểm và đường điều tra dân số nơi chúng được đặt.

Bạn có thể làm điều đó cục bộ thay vì API trực tuyến. Một giải pháp với dữ liệu ranh giới R: C điều tra cho Hoa Kỳ có sẵn thông qua TIGER tại c Điều tra địa chỉ ; nếu bạn đang nghĩ về Hoa Kỳ thì bạn có thể tải xuống các tiểu bang Hoa Kỳ và sử dụng một chức năng để có được các vùng. Ví dụ: sử dụng California và một điểm ngẫu nhiên (hoặc ngẫu nhiên):

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
}

Dùng thử:

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

Cung cấp 0004, đó là chính xác.

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

Điều này hoạt động, nhưng mất 5 giây trên máy Mac của tôi, có lẽ sẽ quá nhiều nếu bạn sẽ làm nhiều việc này; Tôi nghi ngờ ai đó sẽ cân nhắc trong một thời gian ngắn với giải pháp PostGIS sẽ nhanh hơn gấp triệu lần ...

Như được đề cập bởi @Simbamangu, điều này có thể dễ dàng thực hiện tại địa phương chỉ bằng bất kỳ phần mềm GIS trên máy tính để bàn nào trên thị trường.

Đối với nguồn mở và miễn phí, hãy xem " tham gia thuộc tính theo vị trí " trong QGIS .

Nếu bạn có kế hoạch cần thực hiện việc này nhiều lần theo thời gian, bạn có thể xem xét đưa các bộ dữ liệu vào PostGIS và tạo một liên kết không gian (xem ví dụ này trong trợ giúp trực tuyến). Sau đó, bạn có thể hình dung các kết quả trong QGIS.