3 # © Copyright 2022, Scott Gasch
5 """Wrapper around US Census address geocoder API described here:
6 https://www2.census.gov/geo/pdfs/maps-data/data/Census_Geocoder_User_Guide.pdf
7 https://geocoding.geo.census.gov/geocoder/Geocoding_Services_API.pdf
12 --form benchmark=2020 \
13 https://geocoding.geo.census.gov/geocoder/locations/addressbatch \
14 --output geocoderesult.csv
19 from typing import Any, Dict, List, Optional
22 from requests.utils import requote_uri
26 logger = logging.getLogger(__name__)
29 def geocode_address(address: str) -> Optional[Dict[str, Any]]:
30 """Send a single address to the US Census geocoding API. The response
31 is a parsed JSON chunk of data with N addressMatches in the result
32 section and the details of each match within it. Returns None on error.
34 >>> json = geocode_address('4600 Silver Hill Rd,, 20233')
35 >>> json['result']['addressMatches'][0]['matchedAddress']
36 '4600 SILVER HILL RD, WASHINGTON, DC, 20233'
38 >>> json['result']['addressMatches'][0]['coordinates']
39 {'x': -76.92743, 'y': 38.84599}
42 url = 'https://geocoding.geo.census.gov/geocoder/geographies/onelineaddress'
43 url += f'?address={address}'
44 url += '&returntype=geographies&layers=all&benchmark=4&vintage=4&format=json'
45 url = requote_uri(url)
46 logger.debug('GET: %s', url)
49 except Exception as e:
53 if r.status_code != 200:
55 logger.error('Unexpected response code %d, wanted 200. Fail.', r.status_code)
57 logger.debug('Response: %s', json.dumps(r.json(), indent=4, sort_keys=True))
61 def batch_geocode_addresses(addresses: List[str]):
62 """Send up to addresses for batch geocoding. Each line of the input
63 list should be a single address of the form: STREET ADDRESS, CITY,
64 STATE, ZIP. Components may be omitted but the commas may not be.
65 Result is an array of the same size as the input array with one
66 answer record per line. Returns None on error.
68 This code will deal with requests >10k addresses by chunking them
69 internally because the census website disallows requests > 10k lines.
71 >>> batch_geocode_addresses(
73 ... '4600 Silver Hill Rd, Washington, DC, 20233',
74 ... '935 Pennsylvania Avenue, NW, Washington, DC, 20535-0001',
75 ... '1600 Pennsylvania Avenue NW, Washington, DC, 20500',
76 ... '700 Pennsylvania Avenue NW, Washington, DC, 20408',
79 ['"1"," 4600 Silver Hill Rd, Washington, DC, 20233","Match","Exact","4600 SILVER HILL RD, WASHINGTON, DC, 20233","-76.92743,38.84599","76355984","L","24","033","802405","2004"', '"2"," 935 Pennsylvania Avenue, NW, Washington, DC","No_Match"', '"3"," 1600 Pennsylvania Avenue NW, Washington, DC, 20500","Match","Exact","1600 PENNSYLVANIA AVE NW, WASHINGTON, DC, 20500","-77.03534,38.898754","76225813","L","11","001","980000","1034"', '"4"," 700 Pennsylvania Avenue NW, Washington, DC, 20408","Match","Exact","700 PENNSYLVANIA AVE NW, WASHINGTON, DC, 20408","-77.02304,38.89362","76226346","L","11","001","980000","1025"']
83 url = 'https://geocoding.geo.census.gov/geocoder/geographies/addressbatch'
84 payload = {'benchmark': '4', 'vintage': '4'}
86 for chunk in list_utils.shard(addresses, 9999):
89 raw_file += f'{n}, {address}\n'
91 files = {'addressFile': ('input.csv', raw_file)}
92 logger.debug('POST: %s', url)
94 r = requests.post(url, files=files, data=payload)
95 except Exception as e:
99 if r.status_code != 200:
101 logger.error('Unexpected response code %d, wanted 200. Fail.', r.status_code)
103 logger.debug('Response: %s', r.text)
104 for line in r.text.split('\n'):
111 if __name__ == '__main__':