--- title: Explore and Download keywords: fastai sidebar: home_sidebar summary: "In this tutorial, the basics of Colabs are introduced and an American Community Survey (ACS) dataset is downloaded." description: "In this tutorial, the basics of Colabs are introduced and an American Community Survey (ACS) dataset is downloaded." ---
{% raw %}
/content/drive/My Drive/colabs/dataplay/dataplay/acsDownload.py:26: FutureWarning: Passing a negative integer is deprecated in version 1.0 and will not be supported in future version. Instead, use None to not limit the column width.
  pd.set_option('display.max_colwidth', -1)
/usr/local/lib/python3.6/dist-packages/psycopg2/__init__.py:144: UserWarning: The psycopg2 wheel package will be renamed from release 2.8; in order to keep installing from binary please use "pip install psycopg2-binary" instead. For details see: <http://initd.org/psycopg/docs/install.html#binary-install-from-pypi>.
  """)
{% endraw %}

This Coding Notebook is the first in a series.

An Interactive version can be found here Open In Colab.

This colab and more can be found at https://github.com/BNIA/colabs

  • Content covered in previous tutorials will be used in later tutorials.

  • New code and or information should have explanations and or descriptions attached.

  • Concepts or code covered in previous tutorials will be used without being explaining in entirety.

  • If content can not be found in the current tutorial and is not covered in previous tutorials, please let me know.

  • This notebook has been optimized for Google Colabs ran on a Chrome Browser.

  • Statements found in the index page on view expressed, responsibility, errors and ommissions, use at risk, and licensing extend throughout the tutorial.

About this Tutorial:

Whats inside?

The Tutorial

In this notebook, the basics of Colabs are introduced.

  • We will explore ACS data catalogs to locate data we like
  • We will programmatically download data from the American Community Survey (ACS)
    • Examples will use ACS table B19001, Baltimore City 2017 estimates
  • We will rework the datasets to be human friendly

Objectives

By the end of this tutorial users should have an understanding of:

  • Google Colabs
  • Census Data
  • Exploring, Retrieveing, and cleaning ACS Data programmatically
  • the 'retrieve_acs_data()' function, and how to use it in the future

Using Colabs:

Instructions: Read all text and execute all code in order.

How to execute code:

  • Locate labels taking the form: 'Run: (A Short Description)'
  • Left of this text you will see an open bracket [ ], possibly with a number inside it.
    • Hovering over the brackets will reveal a play button.
    • Click the button to execute code.
  • Alternately: Click on a box with code and hit 'Shift' + 'Enter'

If you would like to see the code you are executing, double click the label 'Run: '. Code is accompanied with brief descriptions inlined.

Try It! Go ahead and try running the cell below. What you will be shown as a result is a flow chart of how this current tutorial may be used.

{% raw %}
#@title Run: View User Path

%%html
<img src="https://charleskarpati.com/images/viewuserpath_short.png">
{% endraw %}

Background

About The Census Data

Census or ACS Data?

Census data comes in 2 flavors:


1) American Community Survey (ACS)

  • First released to the public in 2006
  • Derived using 5 Year Estimates (Past 5 years of data)
  • Delievered Annually
  • Delievered at Tract Level. (READ -> 'Geographic Granularity')
  • Margin of error generally prohibit a more granular view
  • ACS Data is accessible programmatically (READ -> 'ACS Programmatic Retrieval')(READ -> 'Geographic Reference Code')

2) Decienial Census

  • Estimates usings 10 years of ACS data
  • Decenial Census data are created, in part, by ACS data.
  • Delievered once every 10 years.
  • Delivered at Block Level. Most Accurate

Geographic Granularity

Census data can come in a variety of levels.

These levels define the specificity of the data.

Ie. Weather a data is reporing on individual communities, or entire cities is contingent on the data granularity.

The data we will be downloading in this tutorial, ACS Data, can be found at the Tract level and no closer.

Aggregating Tracts is the way BNIA calculates some of their yearly community indicators!

Each of the bolded words in the content below are levels that are identifiable through a (READ -> 'Geographic Reference Code') .

  • A census block is the smallest unit of measurement used by the Census
  • Information by census block is only available decenially (i.e. not ACS data)
  • Block groups are the next smallest unit of measurement used by the census and are composed of aggregate census blocks
  • Census tracts are composed of block groups and are the next largest unit of measurement used by the ACS
  • County, city and census designated places are composed of Tracts

For more information on Geographic Reference Codes, refer to the table of contents for the section on that matter.

Run the following code to see how these different levels nest into eachother!

{% raw %}
#@title Run: Census Granularities

%%html
<img src="https://charleskarpati.com/images/census_granularities.png">
{% endraw %}

Geographic Reference Codes

State, County, and Tract ID's are called Geographic Reference Codes.

This information is crucial to know when accessing data.

In order to successfully pull data, Census State and County Codes must be provided.

The code herin is configured by default to pull data on Baltimore City, MD and its constituent Tracts.

In order to find your State and County code:


Either

A) Click the link: https://geocoding.geo.census.gov/geocoder/geographies/address where upon entering a unique address you can locate state and county codes under the associated values 'Counties' and 'State'

OR

B) Conversly, click https://www.census.gov/geographies/reference-files/time-series/geo/tallies.html

  • The Geographies mainpage contains a lot of data assets.
  • The link to these tallies was located by accessing the geographical references subdirectory of the geographies mainpage and then filtered for publications made on the year 2010 (We are using the 2010 census boundaries)
  • Once clicked, simply scroll down to where you find the header 'Tallies of Geographic Entities By State'
  • Enter your state and press enter.
  • You will be redirected to a plain text file that will contain all the information on your state and its counties.

Working with the ACS Data

Searching for a dataset is the first step in the data processing pipeline.

In this tutorial we plan on processing ACS data in a programmatic fashion.

This tutorial will not just allow you to search/ explore ACS tables and inspect their contents (attributes), but also to download, format, and clean it!

Search Advice

Despite a table explorer section being provided, it is not suggested you use this approach, but rather, explore available data tables and retrieve their ID's using the dedicated websites provided below:

American Fact Finder may assist you in your data locating and download needs: https://factfinder.census.gov/faces/nav/jsf/pages/index.xhtml Fact Finder provides a nice interface to explore available datasets. From Fact Finder you can grab a Table's ID and continue the tutorial. Alternately, from Fact Finder, You can download the data for your community directly via an interface. From there, you may continue the tutorial by loading the downloaded dataset as an external resource, instructions on how to do this are provided further below in this tutorial.

Update : 12/18/2019 " American FactFinder (AFF) will remain as an "archive" system for accessing historical data until spring 2020. " - American Fact Finder Website

The New American Fact Finder : https://data.census.gov/cedsci/

This new website is provided by the Census Org. Within its 'Advanced Search' feature exist all the filtering abilities of the older, depricated, (soon discontinued) American Fact Finder Website. It is still a bit buggy to date and may not apply all filters. Filters include years(can only pick on year at a time), geography(state county tract), topic, surveys and Table ID. The filters you apply are shown at the bottom of the query and submitting the search will yield data tables ready for download as well as table ID's that you may snag for use in this tutorial.

ACS Programmatic Retrieval

Developer Resources

Notes on the Census API

Tutorial Notes:

  • Details and Subject tables are derived using the 5 year ACS data.
    • As a reminder, estimates using (ACS) 5-year estimates arrive at the tract level
  • These tables are created by the census and are pre-compiled views of the data.

  • The Detail Tables contain all possible ACS Data.
  • The Subjects Table contains ACS data in convenient groups
  • BNIA create their data mostly using Details table, but sometimes pulling the data from a Subject Table is more convenient (the data would otherwise be found along multiple details tables).

ACS Website Notes:

  • Detailed Tables contain the most detailed cross-tabulations, many of which are published down to block groups. The data are population counts. There are over 20,000 variables in this dataset.

  • Subject Tables provide an overview of the estimates available in a particular topic. The data are presented as population counts and percentages. There are over 18,000 variables in this dataset.

For more Information (via API) Please Visit

Guided Walkthrough

SETUP

You will need to this all installed and imported in order for anything following it to work

pip install ipywidgets geopandas
import ipywidgets as widgets
from IPython.core.interactiveshell import InteractiveShell
import ipywidgets as widgets
from ipywidgets import interact, interact_manual
import urllib.request as urllib 
from urllib.parse import urlencode 
import socket 
import pandas as pd
import json 
import numpy as np
from pandas.io.json import json_normalize
import csv
import geopandas as gpd
import psycopg2,pandas,numpy
from shapely import wkb
from shapely.wkt import loads
import os
import sys
import fiona
import matplotlib.pyplot as plt
import glob
import imageio
{% raw %}
/usr/local/lib/python3.6/dist-packages/ipykernel_launcher.py:22: FutureWarning: Passing a negative integer is deprecated in version 1.0 and will not be supported in future version. Instead, use None to not limit the column width.
/usr/local/lib/python3.6/dist-packages/psycopg2/__init__.py:144: UserWarning: The psycopg2 wheel package will be renamed from release 2.8; in order to keep installing from binary please use "pip install psycopg2-binary" instead. For details see: <http://initd.org/psycopg/docs/install.html#binary-install-from-pypi>.
  """)
{% endraw %}

(Optional) Local File Access

You can access Google Drive directories:

You can also import file directly into a temporary folder in the virutal colab enviornment

By default you are positioned in the ./content/ folder.

Explore Table Directories

Please Note: The following section details a programmatic way to access and explore the census data catalogs. It is advised that rather than use this portion of the section of the tutorial, you read the section 'Searching For Data' --> 'Search Advice' above and which provide links to dedicated websites hosted by the census bureaue explicitly for your data exploration needs!

Explore the Detailed Table Directory

Retrieve and search available ACS datasets through the ACS's table directory.

The table directory contains TableId's and Descriptions for each datatable the ACS provides.

By running the next cell, an interactive searchbox will filter the directory for keywords within the description.

Be sure to grab the TableId once you find a table with a description of interest.

response = urllib.urlopen('https://api.census.gov/data/2017/acs/acs5/groups/')
metaDataTable = json_normalize( json.loads(response.read())['groups'] )
metaDataTable.set_index('name', drop=True, inplace=True)
description = input("Search ACS Table Directory by Keyword: ")
metaDataTable[ metaDataTable['description'].str.contains(description.upper()) ]

Once you a table from the explorer has been picked, you can inspect its column names in the next part.

This will help ensure it has the data you need!

tableId = input("Please enter a Table ID to inspect: ")
url = f'https://api.census.gov/data/2017/acs/acs5/groups/{tableId}.json'
metaDataTable = pd.read_json(url).reset_index(inplace = True, drop=False) 
metaDataTable = pd.merge(
    json_normalize(data=metaDataTable['variables']), 
    metaDataTable['index'] , left_index=True, right_index=True )
metaDataTable = metaDataTable[['index', 'concept']].dropna(subset=['concept'])

Explore the Subject Table Directory

The Data Structure we recieve is different than the prior table.

Intake and processing is different as a result.

Now lets explore what we got, just like before.

Only difference is that the column names are automatically included in this query.

url = 'https://api.census.gov/data/2017/acs/acs5/subject/variables.json'
data = json.loads(urllib.urlopen(url).read())['variables']
objArr = []
for key, value in data.items():
  value['name'] = key
  objArr.append(value)
metaDataTable = json_normalize(objArr).set_index('name', drop=True, inplace=True)
metaDataTable = metaDataTable[ ['attributes', 'concept', 'group', 'label', 'limit', 'predicateType' ] ]
concept = input("Search ACS Subject Table Directory by Keyword")
metaDataTable[ metaDataTable['concept'].str.contains(concept.upper(), na=False) ]

Get Table Data

Intro

Hopefully, by now you know which datatable you would like to download!

The following Python function will do that for you.

  • It can be imported and used in future projects or stand alone.
{% raw %}
{% endraw %} {% raw %}

retrieve_acs_data[source]

retrieve_acs_data(state, county, tract, tableId, year, saveOriginal, save)

{% endraw %}

Function Explanation

Description: This function returns ACS data given appropriate params.

Purpose: Retrieves ACS data from the web

Services

  • Download an ACS dataset from an Subject (S) table
  • Download an ACS dataset from a Details (B) table

Input:

  • state
  • county
  • tract
  • tableId
  • year
  • saveAcs

Output:

  • Acs Data.
  • Prints to ../../data/2_cleaned/acs/

How it works

  • Before our program retrieve the actual data, it will want the table's metadata.

    • This metadata will be used as a crosswalk to replace the awkward column names
    • If this is not done, only a column ID would denote each column. not human readable.
  • The Function changes the URL it requests data from depending on if it is an S or B type table the user has requested

  • Multiple calls for data must be made as a single table may have several hundred columns in them.

    • Constructing a table requires merging the data from multiple responces
  • Our program not just pulls tract level data but the aggregate for the county.

    • County totals are included automatically as 'tract 010000'.
      • The County total is not the sum of all other tracts but a seperate, indendent and unique query.
    • Tract and County Datatables must be merged to form a single dataset
  • Finally, we will download the data in two different formats if desired.

  • If we choose to save the data, we save it with the Table IDs + ColumnNames, and once without the TableIDs.

Function Diagrams

{% raw %}
#@title Run: Class Diagram retrieve_acs_data()

%%html
<img src="https://charleskarpati.com/images/class_diagram_retrieve_acs_data.png">
{% endraw %} {% raw %}
#@title Run: retrieve_acs_data Flow Chart

%%html
<img src="https://charleskarpati.com/images/flow_chart_retrieve_acs_data.png">
{% endraw %} {% raw %}
#@title Run: Gannt Chart  retrieve_acs_data()

%%html
<img src="https://charleskarpati.com/images/gannt_chart_retrieve_acs_data.png">
{% endraw %} {% raw %}
#@title Run: Sequence Diagram  retrieve_acs_data()

%%html
<img src="https://charleskarpati.com/images/sequence_diagram_retrieve_acs_data.png">
{% endraw %}

Function Examples

Now use this function to Download the Data!

{% raw %}
# Our download function will use Baltimore City's tract, county and state as internal paramters
# Change these values in the cell below using different geographic reference codes will change those parameters
tract = '*'
county = '153' # '059' # 153 '510'
state = '51'

# Specify the download parameters the function will receieve here
tableId = 'B19049' # 'B19001'
year = '17'
saveAcs = True
{% endraw %} {% raw %}
# state, county, tract, tableId, year, saveOriginal, save

df = retrieve_acs_data(state, county, tract, tableId, year, saveAcs, True)
df.head()
Number of Columns 5
B19049_001E_Median_household_income_in_the_past_12_months_(in_2017_inflation-adjusted_dollars)_Total B19049_002E_Median_household_income_in_the_past_12_months_(in_2017_inflation-adjusted_dollars)_Householder_under_25_years B19049_003E_Median_household_income_in_the_past_12_months_(in_2017_inflation-adjusted_dollars)_Householder_25_to_44_years B19049_004E_Median_household_income_in_the_past_12_months_(in_2017_inflation-adjusted_dollars)_Householder_45_to_64_years B19049_005E_Median_household_income_in_the_past_12_months_(in_2017_inflation-adjusted_dollars)_Householder_65_years_and_over state county tract
NAME
Census Tract 9002.01, Prince William County, Virginia 68250 -666666666 54900 86705 128077 51 153 900201
Census Tract 9002.03, Prince William County, Virginia 51329 -666666666 52469 58958 37125 51 153 900203
Census Tract 9014.13, Prince William County, Virginia 124801 -666666666 121402 127206 -666666666 51 153 901413
Census Tract 9014.11, Prince William County, Virginia 163086 -666666666 155556 190333 108542 51 153 901411
Census Tract 9007.01, Prince William County, Virginia 68304 34120 79868 66863 57885 51 153 900701
{% endraw %}