THA5 - Mapping

Useful CSE 163 Resources¶

• CSE 163 Code Quality Guide
• Common Python Errors
• Common THA Errors
• flake8 References
• Testing Tips

Learning objective: Analyze and plot geospatial data to investigate food deserts in Washington.

You can find the starter code here. Make sure to extract (unzip) the contents anywhere on your computer. If you are working in VS Code, navigate to File | Open or Open Folder, then select the hw5 folder.

• hw5.py is the file for you to put your implementations of the plotting functions and to call them in main.

• hw5_writeup.md is the file for your writeup. Instead of testing, this assessment emphasizes reflection on our data analysis.

• expected is a folder containing the expected output for the 5 plotting functions you’ll be writing. Don’t modify the contents or name of this folder.

• food_access.zip is a zip containing all relevant datasets to this homework. You will need to unzip the files and store them in the same folder as your other files.

• cse163_imgd.py is a helper file that checks your plot outputs against expected output, and creates an image showing any pixel differences.

• Playground.ipynb is a Jupyter Notebook playground. Feel free to edit this Jupyter Notebook to prototype ideas and explore the data.

Context¶

“Food deserts” are neighborhoods where residents do not have nearby access to grocery stores offering affordable and nutritious food. In a June 2009 report to the US Congress, the US Department of Agriculture reports:

According to data from the latest census (2000), about 23.5 million people, or 8.4 percent of the U.S. population, live in low-income neighborhoods that are more than a mile from a supermarket. Low-income neighborhoods are areas where more than 40 percent of the population has income less than or equal to 200 percent of the Federal poverty threshold ($44,000 per year for a family of four in 2008).

In this assessment, we’ll join 2010 US census data with food access data to investigate food deserts in Washington state. A census tract is defined as a food desert if enough people in the tract do not have nearby access to food sources.

• In urban areas, “low access” is defined as 0.5 miles, which we’ve stored in the column lapophalf.

• In rural areas, “low access” is defined as 10 miles, which we’ve stored in the column lapop10.

The 2010 US census dataset is geospatial data in shapefile format. The only columns you need to understand are CTIDFP00, the census tract identifier, and geometry, the geometric shape of the tract.

The food access dataset is tabular data in CSV format. Each row in the dataset corresponds to a census tract for every state in the country. The data has many columns, but you only need to understand the following:

• CensusTract is the census tract identifier.

• State is the state name for the census tract.

• County is the county name for the census tract.

• Urban is a flag (0 or 1) that indicates if this census tract is an urban environment.

• Rural is a flag that indicates if this census tract is a rural environment.

• LATracts_half is a flag that indicates if this census tract is “low access” in a half mile radius.

• LATracts10 is a flag that indicates if this census tract is “low access” in a 10 mile radius.

• POP2010 is the number of people in this census tract according to the 2010 census.

• lapophalf is the number of people in this census tract considered having “low access” in a half mile radius.

• lapop10 is the number of people in this census tract considered having “low access” in a 10 mile radius.

• lalowihalf is similar to lapophalf but only counts people considered low access and low income.

• lalowi10 is similar to lapop10 but only counts people considered low access and low income.

Technical Component¶

load_in_data¶

Task: Write a function load_in_data that takes two parameters, the filename for the census dataset and the filename for the food access dataset. load_in_data should merge the two datasets on CTIDFP00 / CensusTract and return the result as a GeoDataFrame. Assume the census identifier column names exist, but don’t assume any other columns in the datasets. The resulting merged dataset might have missing data. For the provided datasets, the shape of the resulting GeoDataFrame is (1318, 30) rows by columns.

percentage_food_data¶

Task: Write a function percentage_food_data that takes the merged data and returns the percentage of census tracts in Washington for which we have food access data. The percentage should be a float between 0 and 100. Do not round the result.

plot_map¶

Task: Write a function plot_map that takes the merged data and plots the shapes of all the census tracts in Washington in a file map.png. Give the plot a title of “Washington State” with plt.title(). Do not customize this plot or overlay any data.

plot_population_map¶

Task: Write a function plot_population_map that takes the merged data and plots the shapes of all the census tracts in Washington in a file population_map.png where each census tract is colored according to population. There will be some missing census tracts. Under the census tracts, plot the map of Washington in the background color #EEEEEE. Include a legend to indicate the meaning of each census tract color and give the plot a title of “Washington Census Tract Populations” with plt.title().

plot_population_county_map¶

Task: Write a function plot_population_county_map that takes the merged data and plots the shapes of all the census tracts in Washington in a file county_population_map.png where each county is colored according to population. This will involve aggregating all the census tract data in each county, and there will be some missing counties. Under the census tracts, plot the map of Washington in the background color #EEEEEE. Include a legend to indicate the meaning of each census tract color and give the plot a title of “Washington County Populations” with plt.title().

plot_food_access_by_county¶

Task: Write a function plot_food_access_by_county that takes the merged data and produces 4 plots on the same figure showing information about food access across income level.

First, compute the ratio of people in each category.

1. Slice the dataframe to keep only the columns County, geometry, POP2010, lapophalf, lapop10, lalowihalf, lalowi10.

2. Aggregate this dataset by County, summing up all the numeric columns.

3. For each County, compute the ratio of people in each category: lapophalf_ratio, lapop10_ratio, lalowihalf_ratio, lalowi10_ratio. Add these columns to the sliced copy of the dataset to simplify access in later steps.

Then, set up the figure using subplots.

fig, [[ax1, ax2], [ax3, ax4]] = plt.subplots(2, 2, figsize=(20, 10))

Finally, plot the data on each ax1 ax2 ax3 ax4 subplot axis. For each subplot:

1. Plot the map of Washington in the background color #EEEEEE.

2. Call the plot function on the appropriate column with the arguments ax (to specify the subplot axis) and vmin=0 and vmax=1 (so they all share the same scale). Include a legend.

3. Set the titles for each axis to match the expected output picture, such as ax1.set_title('Low Access: Half').

Save the figure to a file county_food_access.png.

plot_low_access_tracts¶

Task: Write a function plot_low_access_tracts that takes the merged data and plots all census tracts considered “low access” in a file low_access.png.

First, compute the above statistics for each census tract depending on its classification as Urban or Rural. Don’t use the LATracts_half or LATracts10 columns in the GeoDataFrame since this computation is slightly more subtle. The thresholds for urban tracts are considered differently from rural tracts.

• Urban: If the census tract is “urban”, the distance of interest is half a mile from a food source. The threshold for low access in an urban census tract is at least 500 people or at least 33% of the people in the census tract being more than half a mile from a food source. An urban census tract that satisfies either of these conditions is considered low access.

• Rural (i.e. non-urban): If the census tract is “rural”, the distance of interest is 10 miles from a food source. The threshold for low access in a rural census tract is at least 500 people or at least 33% of the people in the census tract being more than 10 miles from a food source. A rural census tract that satisfies either of these conditions is considered low access.

Then, produce a layered plot (all on the same axes) to highlight low access census tracts. Each plot will draw on top of the previous plot.

1. Plot the map of Washington in the background with color #EEEEEE.

2. Plot all the census tracts for which we have food access data in the color #AAAAAA.

3. Plot all the census tracts considered low access in the default color, blue.

Give it a title of “Low Access Census Tracts” with plt.title().

Writeup¶

Task: In hw5_writeup.md, apply critical thinking to address the following questions about data collection and analysis. You could spend an entire course talking about any of these topics, but we’re looking for at least 2 to 4 sentences on each question.

1. Is the plot produced by plot_food_access_by_county an effective visualization? Why or why not?

2. Is the plot produced by plot_low_access_tracts an effective visualization? Why or why not?

3. What is one way that government officials could use plot_food_access_by_county or plot_low_access_tracts to shape public policy on how to improve food access?

4. What is a limitation or concern with using the plots in this way? What information might be missing from these plots?

Creative Component: Accessible Maps¶

An important consideration when making any kind of visualization is ensuring that our visualizations are accessible to all users, including those with visual impairments or color vision deficiencies.

Task: Choose 2 plots from your Technical Component implementations and conduct an accessibility audit on each. After identifying accessibility issues, modify your code to improve the accessibility of these visualizations.

Requirements¶

For each of your 2 chosen plots:

1. Audit: Use the provided accessibility tools to evaluate your plot for common accessibility issues, such as:

   • Color contrast and distinguishability
   • Colorblind-friendly palettes
   • Font size and readability
   • Alternative text descriptions

2. Modify: Update your code to address the accessibility issues you identified. Your modifications might include:

   • Adjusting color schemes to improve contrast or accommodate color vision deficiencies
   • Increasing font sizes for better readability
   • Adding or improving alt text descriptions that convey the key insights of the visualization (note: matplotlib does not currently support adding alt text directly to a plot, but you can add alt text to images in Markdown!)
   • Modifying plotting elements (borders, patterns, labels) to improve visual clarity

3. Explain: In the provided markdown cells in hw5_creative.ipynb, for each plot, write 3-5 sentences explaining:

   • What accessibility issues you identified
   • What specific changes you made to your code
   • How these changes improve accessibility for users with different needs

You do not need to save your plots to new files, but they should be visible in the output of your code!

Quality¶

Assessment submissions should pass these checks: flake8 and code quality guidelines. The code quality guidelines are very thorough. For this assessment, the most relevant rules can be found in these sections (new sections bolded):

• Naming

• Documentation

• Efficiency and Redundancy

Submission¶

Submit your work by uploading the following files to Gradescope:

• hw5.py

• hw5_writeup.md

• hw5_creative.ipynb

• map.png

• population_map.png

• county_population_map.png

• county_food_access.png

• low_access.png

Please make sure you are familiar with the resources and policies outlined in the syllabus and the take-home assessments page.

THA 5 - Mapping

Initial Submission by Thursday 02/26 at 11:59 pm.

Submit on Gradescope