m203-ejscreen.ipynb

{
  "cells": [
    {
      "cell_type": "raw",
      "metadata": {},
      "source": [
        "---\n",
        "title: \"EJScreen: Environmental Justice Screening and Mapping Tool\"\n",
        "format: html\n",
        "author: \"Elaine Famutimi and Camilla Green\"\n",
        "date: \"November 3, 2024\"\n",
        "---"
      ],
      "id": "bcdc4e99"
    },
    {
      "cell_type": "markdown",
      "metadata": {},
      "source": [
        "## Overview\n",
        "\n",
        "This lesson introduces EJSCREEN, the EPA’s mapping tool that integrates environmental and demographic data into a consistent national framework. It enables users to visualize and analyze environmental burdens alongside socioeconomic indicators. The lesson demonstrates features of the tool, including mapping features, generating community reports, and comparing EJ indexes. A school district shapefile from Wayne County is then uploaded to EJScreen, which has boundaries based on school county.\n",
        "\n",
        "The lesson also includes a Python component. Specifically, students will learn how to create bar graphs comparing separate regions using EJScreen-generated data in Python. While the lesson primarily focuses on using the web-based EJScreen tool, the Python portion enhances the analysis capabilities by allowing users to create custom visualizations of the data obtained from EJScreen reports. Participants learn to extract meaningful insights by leveraging EJScreen’s datasets and integrating Python for tasks like querying data and generating custom charts. This approach bridges technical and environmental sciences, fostering practical skills in analyzing equity and environmental health issues using programming.\n",
        "\n",
        "## Learning Objectives\n",
        "\n",
        "After completing this lesson, you should be able to:\n",
        "\n",
        "-   Select and visualize Maps and Places available in EJScreen.\n",
        "-   Use the tools tab to create a report for the area of interest by drawing a polygon.\n",
        "-   Review and interpret the report outputs.\n",
        "-   Use the Tools section to upload census tract shapefile.\n",
        "-   Create bar graphs comparing separate regions from EJScreen generated graphs on Python.\n",
        "-   Interpret and compare different types of indexes.\n",
        "\n",
        "## Introduction\n",
        "\n",
        "### What is EJScreen?\n",
        "\n",
        "EJScreen, the Environmental Justice Screening and Mapping Tool, was developed by the U.S. Environmental Protection Agency (EPA) to help identify and address environmental justice issues. The tool is designed to provide communities, NGOs, and policymakers with vital data on environmental and demographic factors that affect vulnerable populations.\n",
        "\n",
        "### Background on the development of EJScreen\n",
        "\n",
        "The concept of environmental justice gained prominence in the 1980s and 1990s when it became clear that certain communities, especially low-income and minority populations, were disproportionately affected by environmental hazards. This concern led to the signing of Executive Order 12898 in 1994, which directed federal agencies, including the EPA, to address environmental justice in minority and low-income populations \\[https://www.epa.gov/laws-regulations/summary-executive-order-12898-federal-actions-address-environmental-justice\\] (Zahra et al., 2009) (USAID, 2024)\n",
        "\n",
        "### Importance of EJScreen\n",
        "\n",
        "EJScreen is important in empowering communities. Access to EJScreen helps communities understand the specific environmental challenges they face, such as pollution levels or exposure to hazardous substances. By providing this data, the tool empowers residents to advocate for their health, safety, and well-being. According to the EPA, EJScreen allows communities to participate more effectively in public discussions, regulatory processes, and decision-making, backed by concrete data to support their concerns and needs. Environmental justice focuses on ensuring that no group of people, particularly minority and low-income communities, bears a disproportionate share of environmental burdens.\\[https://www.epa.gov/ejscreen/what-ejscreen\\] Environmental and social justice organizations often rely on data to build and support their campaigns. EJScreen provides this kind of data, which can highlight environmental injustices and help NGOs propose solutions and mobilize public support. Many NGOs also use these tools to strengthen their grant applications, demonstrating the specific needs of the communities they serve. This is especially important for organizations advocating for low-income and minority populations, who are often disproportionately impacted by environmental hazards. Policymakers require accurate data to create effective and equitable regulations. EJScreen helps identify areas of concern, prioritize actions, and design policies that better protect vulnerable populations. For example, the EPA uses this tool to inform regulatory actions, compliance monitoring, and enforcement activities.\\[https://www.epa.gov/environmentaljustice\\]. By making environmental data accessible to the public, EJScreen also promotes transparency and helps hold governments and industries accountable for their environmental impact. This public accessibility allows for trust and collaboration between communities, regulatory agencies, and policymakers, which is crucial for successfully implementing environmental policies \\[https://www.epa.gov/environmentaljustice/national-environmental-justice-advisory-council\\].\n",
        "\n",
        "Environmental issues often require collaborative solutions involving multiple stakeholders, such as government agencies, NGOs, businesses, and community groups. EJScreen provides a common platform for discussing and addressing these challenges, which fosters more effective partnerships and solutions. EJScreen’s development reflects the EPA's commitment to enhancing the accessibility of environmental data, promoting transparency, and addressing environmental justice concerns by equipping all stakeholders with the tools needed to assess and respond to environmental health disparities.\n",
        "\n",
        "::: callout-note\n",
        "## Knowledge Check\n",
        "\n",
        "What is the importance of EJScreen?\n",
        "\n",
        "-   Helping environmental and social justice organizations support their campaigns\n",
        "\n",
        "-   Providing policymakers with data to create effective regulations\n",
        "\n",
        "-   Empowering residents to advocate for their health and safety\n",
        "\n",
        "-   Providing a common platform for collaborative solutions between multiple stakeholders\n",
        "\n",
        "-   All of the above\n",
        ":::\n",
        "\n",
        "## How to access EJScreen\n",
        "\n",
        "To begin, click on the EJScreen link here: <https://www.epa.gov/ejscreen>. On that page, along with resources on how to use EJScreen, you will see a blue, bold text labeled “Launch the EJScreen Tool.” You can also go straight to the mapper using this link:[https://ejscreen.epa.gov/mapper/.](https://ejscreen.epa.gov/mapper/)\n",
        "\n",
        "![](docs/m203-ej-images/image9.png){fig-align=\"center\"}\n",
        "\n",
        "Once you click the link, you'll be directed to a Welcome screen. Here, you'll find an introduction, resource links, and a video overview of the tool. We recommend exploring these resources before proceeding. After that, close the Welcome dialog to continue.\n",
        "\n",
        "![](docs/m203-ej-images/image13.png){fig-align=\"center\"}\n",
        "\n",
        "You can use the search bar to explore a region. For this lesson, we will be focusing on Detroit, Michigan.\n",
        "\n",
        "![]()![](docs/m203-ej-images/image15.png){width=\"689\"}\n",
        "\n",
        "Type 'Metro Detroit' in the search bar to focus in on the metropolitan region of Detroit. Alternatively, you can use latitude and longitude coordinates. You can also zoom in or out of a region by clicking on the “+” or “-” icons on the bottom right.\n",
        "\n",
        "![](docs/m203-ej-images/image11.png){fig-align=\"center\"}\n",
        "\n",
        "We can start exploring maps by noticing the Widget toolbar. Hovering your cursor over each icon will show the name of each tab. From right to left to right, the icons are ‘Maps’, ‘Places’, ‘Reports’, and ‘Tools.’ Clicking each icon will display a dropdown list of indicators to choose from. For now, we will focus on the ‘Maps’ icon.\n",
        "\n",
        "![](docs/m203-ej-images/image28.gif){fig-align=\"center\"}\n",
        "\n",
        "Now, click on the Socioeconomic Indicators icon. Clicking that will expose seven types of socioeconomic indicators, and two indexes that can be explored. Click on People of Color, and your mapper should look like this. We will be mapping People of Color compared to the US. One can also compare it to just the state. To better view the map, one can click on the ‘\\<’ on the Map widget.\n",
        "\n",
        "![](docs/m203-ej-images/image2.png){fig-align=\"center\"}\n",
        "\n",
        "Notice the ‘Map Contents’ on the upper right. In statistics, a percentile is a score that shows how a particular score compares to other scores in a dataset. For example, people of color in the 95-100 percentile in a block group means that 95-100% of people living in that area are non-white. You can learn more about how EJScreen uses percentiles here <https://www.epa.gov/ejscreen/how-interpret-ejscreen-data>.\n",
        "\n",
        "![](docs/m203-ej-images/image5.png){fig-align=\"center\"}\n",
        "\n",
        "You can hide the Map Contents by clicking on the sign. Noice the symbols in the Map Contents. To the right of that is the sign has the layer turned on; clicking on it turns it off. The allows you to view the metadata. On the opposite side, is an sign, which allows you to set the transparency of the map. A pop up with a sliding scale will appear when that is clicked. Clicking on removes the layer. Selecting provides a description of the index. Clicking minimizes the Map Contents pane, to expand it one can simply click on it again (this time the arrows will be pointing downwards).\n",
        "\n",
        "::: callout-note\n",
        "## FAIR Standards and Metadata\n",
        "\n",
        "Metadata is structured information that describes a dataset. It is often thought of as \"data about data,\" allowing users to understand the format, content, quality, temporality, and many more characteristics of data. Providing metadata is important for the FAIR (Findability, Accessibility, Interoperability, and Reusability) principles, helping users find and understand data.\n",
        ":::\n",
        "\n",
        "![](docs/m203-ej-images/image26.png){fig-align=\"center\"}\n",
        "\n",
        "One can generate reports and charts given an area of interest. Go to the ‘Reports’ tab. One option to select an area of interest is by drawing a polygon. Click on ‘Draw an Area’.\n",
        "\n",
        "![](docs/m203-ej-images/image12.png){fig-align=\"center\"}\n",
        "\n",
        "Your cursor should turn into a crosshair. You can now click on the map and start drawing your polygon. Move your cursor over to where you want the next point of the polygon to be. To form the final polygon, just connect to the first point you made, and the polygon will turn into a transparent green.\n",
        "\n",
        "![](docs/m203-ej-images/image33.gif){fig-align=\"center\"}\n",
        "\n",
        "A pop-up will appear where you can name your study area, add a buffer, download reports, and explore charts. Name it 'Detroit Study Area' to proceed.\n",
        "\n",
        "![](docs/m203-ej-images/image31.png){fig-align=\"center\"}\n",
        "\n",
        "Clicking on EJScreen Community Report will open up a PDF file with detailed information about the study area. Let us explore this PDF. On the right side, you can find information about other characteristics of the community of interest, such as low income residents, the education level of residents, and the gender breakdown of the community. Right below Community Information is Breakdown By Race, which shows the percent of each race in this area. You will also find the percentage of age brackets, and limited English speaking breakdown.\n",
        "\n",
        "![](docs/m203-ej-images/image3.png){fig-align=\"center\"}\n",
        "\n",
        "On the left side of the document, there is be a map showing the study area, including the legend from the Map Contents tab. Centering the map or zooming in or out on the Mapper will result in a different view on the report. This view will also highlight the languages spoken at home.\n",
        "\n",
        "![](docs/m203-ej-images/image25.png){fig-align=\"center\"}\n",
        "\n",
        "Page 2 of the report presents bar graphs comparing EJ indexes at state and national levels. Review these graphs for differences in environmental factors like drinking water non-compliance and wastewater discharge. Note: Although the graph does not include wastewater discharge, the table below shows values for each index from left to right.\n",
        "\n",
        "![](docs/m203-ej-images/image8.png){fig-align=\"center\"}\n",
        "\n",
        "![](docs/m203-ej-images/image27.png){fig-align=\"center\"}\n",
        "\n",
        "![](docs/m203-ej-images/image4.png){fig-align=\"center\"}\n",
        "\n",
        "Let's compare this report to one generated outside of Detroit. Try to pick an area that shows a lower percentile of people of color, or less red, like Sterling Heights or Westland. What differences do you see? Are the EJ indexes higher or lower?\n",
        "\n",
        "## Importing Image Layers from ESRI Living Atlas into EJScreen\n",
        "\n",
        "You can also access and visualize demographic and environmental data from other sources, including Map Services. The first step to incorporating external data into EJScreen is to access the ESRI Living Atlas, which contains a variety of publicly available datasets. Users can search for relevant layers by entering keywords related to their area of interest, such as air pollution, heat vulnerability, or land use. Each layer in the Living Atlas comes with detailed metadata and a \"Service URL\" (often referred to as the REST URL). This URL is essential for linking the data to other mapping tools like EJScreen.\n",
        "\n",
        "Navigate to <https://livingatlas.arcgis.com/en/home/>\n",
        "\n",
        "Search for 'TEMPO NO2' <https://gis.earthdata.nasa.gov/portal/home/item.html?id=0252904123a74e74a7cff652d52a5b19>[![](images/image27-01.png)](https://gis.earthdata.nasa.gov/portal/home/item.html?id=0252904123a74e74a7cff652d52a5b19)\n",
        "\n",
        "Once you have identified a relevant layer and copied its Service URL in the bottom right of the page, you can add it to EJScreen The EJScreen interface includes an option to \"Add Maps & Data.\" Then, select \"Add Map Services\" where you will paste the copied URL. This action imports the chosen geospatial layer into your EJScreen map. Depending on the type of layer, you may be able to customize its display, such as adjusting transparency, color, or symbology, to better integrate it with the existing EJScreen data.\n",
        "\n",
        "![](images/ezgif-4acdcb73a594b.gif){fig-align=\"center\"}\n",
        "\n",
        "After the Living Atlas layer is successfully added, it can be used in combination with EJScreen's built-in indicators. This allows for deeper analysis by overlaying environmental justice data with external geospatial information, such as pollution levels or land use patterns. By integrating both datasets, users can gain a more comprehensive view of the factors influencing environmental justice in a particular area. Another resource for EJI by CDC: <https://www.atsdr.cdc.gov/place-health/php/eji/eji-explorer.html> To learn more about TEMPO: Air Quality Monitoring from Space by NASA using TEMPO, check out this StoryMap: <https://storymaps.arcgis.com/stories/01e82aefbc8b4d7a951fe089c818bc0c>\n",
        "\n",
        "## Using the \"Explore Charts\" Function to Generate Data for Analysis\n",
        "\n",
        "One can also view other reports, and download excel data. Lets compare the bar graphs in a red population of Detroit and a region in Sterling Heights. Delete the original polygon and create new ones within those regions.\n",
        "\n",
        "![](docs/m203-ej-images/image24.png){fig-align=\"center\"}\n",
        "\n",
        "Click on a polygon and select ‘Explore Charts.’ A pop up will open, showing tabs for Environmental Justice Indexes, Environmental Burden Indicators, Socioeconomic Indicators, and Supplemental Indexes. Placing the cursor over each bar will show the exact percentile for each index. You can also unselect and select indexes of your choosing. For now, let's keep all of them selected. We will also only focus on the USA percentile. Environmental Justice Indexes, Environmental Burden Indicators, and Supplemental Indexes focus on the same environmental indexes. Selecting the Socioeconomic indicators will list different types of indexes focusing on socioeconomic factors.\n",
        "\n",
        "![](docs/m203-ej-images/image7.png){fig-align=\"center\"}\n",
        "\n",
        "One can also explore the Socioeconomic (ACS) Report, which includes data from the US Census Bureau American Community Survey (ACS). Select Get Data Table, and a tabular view of the data should pop up. This will include all the categories, selected variables, percentile in state and percentile in USA.\n",
        "\n",
        "![](docs/m203-ej-images/image16.png){fig-align=\"center\"}\n",
        "\n",
        "Click on the save icon![](docs/m203-ej-images/image22.png) to export the data, and an excel file in .csv format will automatically download.\n",
        "\n",
        "Explore the charts for Sterling Heights.\n",
        "\n",
        "![](docs/m203-ej-images/image14.png){fig-align=\"center\"}\n",
        "\n",
        "Notice the lower percentiles for the EJ Indexes. Is this data in line with what was seen in the generated report? Download the data for this as well.\n",
        "\n",
        "Let's perform some basic data analysis!\n",
        "\n",
        "We can now use this data to create graphs on our own and better visually compare the indexes between separate regions. Let us compare the data for Sterling Heights and the region we selected in the metro area of Detroit. Save the files that were downloaded in a working directory, and rename them so that it makes sense. One file can be named ejscreen_detroit and the other ejscreen_sterlingheights. Open one of the csv files and take note of how the data is formatted. Under the main header that describes the EPA region, there should be 8 columns, each one with its own header (Category, Selected Variables, etc.) The Category column has the category of each index and indicator. The Selected Variables are the indexes for each category. For the EJ Indexes category, you should see 13 indexes in the Selected Variable column. To work with these files, we will be using the pandas library, which allows us to handle and manipulate excel files as a dataframe.\n",
        "\n",
        "Import the pandas and matplotlib libraries. Documentation for pandas can be found here https://pandas.pydata.org/docs/, and for matplotlib, here https://matplotlib.org/stable/users/getting_started/\n"
      ],
      "id": "40c9d0ef"
    },
    {
      "cell_type": "code",
      "metadata": {},
      "source": [
        "import pandas as pd\n",
        "import matplotlib.pyplot as plt"
      ],
      "id": "ab9627cf",
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {},
      "source": [
        "In order to use the data that you just downloaded as a CSV file, you have to tell the computer to read the CSV. The pandas read_csv() function reads CSV files. You can specify if the data has a header by adding header = 1, otherwise it treats the column names as part of the data. By specifying the header, it is extracting the column names we noticed in the CSV file and treating them as the header. Reading these files will now turn the csv files into a dataframe.\n"
      ],
      "id": "bf0b3913"
    },
    {
      "cell_type": "code",
      "metadata": {},
      "source": [
        "df_detroit = pd.read_csv(r'data/ejscreen_detroit.csv', header = 1)"
      ],
      "id": "f65fe0f4",
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {},
      "source": [
        "Printing out the shape of the dataframe lets one see the number of rows and columns the data has. It is important to know the dimensions of your data before performing any analysis on it, to ensure that you will be correctly manipulating the data. You can also print just the dataframe variable for Sterling Heights and/or Detroit. The output looks like the contents in the CSV file.\n"
      ],
      "id": "a6a68957"
    },
    {
      "cell_type": "code",
      "metadata": {},
      "source": [
        "print(df_detroit.shape)\n",
        "\n",
        "#Read in Sterling Heights csv file.\n",
        "df_sterlingheights = pd.read_csv(r'data/ejscreen_sterlingheights.csv', header = 1)\n",
        "\n",
        "#Print the dataframe to examine the contents.\n",
        "print(df_sterlingheights)"
      ],
      "id": "9be155b8",
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {},
      "source": [
        "We are going to merge the dataframes using the pandas merge() function. The merge function combines dataframes using a chosen parameter. For this case, we will be merging using the ‘#’, ‘Category’, and ‘Selected Variables’ column. The ‘on’ parameter chooses which column both dataframes will join on; the columns must be found in both dataframes. Since the ‘#’, ‘Category’, and ‘Selected Variables’ columns are the same in both dataframes, we will merge on those columns, meaning that the final dataframe will only have one of each of those columns. The rest of the columns will double, and the numbers of rows will remain the same. The ‘suffixes’ parameter allows you to add a suffix to the columns to help identify which column belongs to which dataset. Since the ‘Value’, ‘State Avg.’, %ile in State, USA Avg., and %ile in USA differ between the two datasets, those column headers will carry the suffixes.\n"
      ],
      "id": "1b8f7802"
    },
    {
      "cell_type": "code",
      "metadata": {},
      "source": [
        "#Suffix ‘_dt’ for Detroit; ‘_sh’ for Sterling Heights.\n",
        "df_merge = pd.merge(df_detroit, df_sterlingheights, on = ('#', 'Category', 'Selected Variables'), suffixes = ('_dt', '_sh'))\n",
        "print(df_merge)"
      ],
      "id": "e691dab5",
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {},
      "source": [
        "Extract the 'EJ Index' rows from the ‘Category’ columns, the corresponding '%ile in State' values and the ‘Selected Variables’ for the Detroit and Sterling Heights.\n"
      ],
      "id": "6b48843e"
    },
    {
      "cell_type": "code",
      "metadata": {},
      "source": [
        "#Extracting only the rows that have ‘EJ Index’ as the category from the merged dataframe, and turning it into a new dataframe called ‘ej_index_df’ that only has that category. The rest of the variables will be extracted from this dataframe.\n",
        "ej_index_df = df_merge[(df_merge['Category'] == 'EJ Index')]\n",
        "\n",
        "#Extracting the '%ile in State' for Detroit and Sterling Heights.\n",
        "ej_index_detroit = ej_index_df['%ile in State_dt']\n",
        "ej_index_sterlingheights = ej_index_df['%ile in State_sh']\n",
        "\n",
        "#Extracting the ‘Selected Variables.’\n",
        "selected_variables = ej_index_df['Selected Variables']"
      ],
      "id": "f874c8ec",
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {},
      "source": [
        "Notice how each of the indexes in the selected variables dataframe starts with ‘EJ Index for’ before the name of the index. When we plot this, the x axis will have the label of ‘EJ Index’, so having it before each index seems redundant. We can remove the first three words by using the .split() function. We will create a function that removes the first three words of each of the indexes in the selected_variables dataframe, and then returns the remaining words. Each selected variable is a string, which in this case, is a sequence of words. The .split() function, without specifying any parameters, splits the string into individual words. The .join() function then combines the separate words back into a string.\n",
        "\n",
        "::: callout-note\n",
        "A function is a unit of code that performs a specific task. Its main advantage is that it can be reused, making code more efficient. def() is the keyword used to define the function, and is placed before the name of the function. What goes in the parentheses is the optional parameters, which are the input values. Sometimes, a function ends with ‘return’, which outputs a specific value from the function.\n",
        ":::\n"
      ],
      "id": "a6450057"
    },
    {
      "cell_type": "code",
      "metadata": {},
      "source": [
        "def remove_first_three_words(text):\n",
        "\t#Splitting the string in between spaces.\n",
        "\twords = text.split()\n",
        "\t#Returning a new string that excludes the first three words.\n",
        "\treturn ' '.join(words[3:])\n",
        "\n",
        "#Applying the function to the selected_variables dataframe.\n",
        "selected_variables = selected_variables.apply(remove_first_three_words)"
      ],
      "id": "af23cd29",
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {},
      "source": [
        "A bar graph can now be created using the information that was extracted. The ‘index’ variable is storing a sequence of indices from 0 to the length of selected_variables minus 1. The len() function calculates the size of the selected_variables dataframe, which is 13 since there are 13 indexes. The range() function then creates a sequence of numbers starting from 0 to minus 1 of the length of selected variables.\n"
      ],
      "id": "1ebb851a"
    },
    {
      "cell_type": "code",
      "metadata": {},
      "source": [
        "#Configuring size of bar graph.\n",
        "plt.figure(figsize=(12, 8))\n",
        "bar_width = 0.35\n",
        "\n",
        "#Creating sequence of numbers from 0 to 12 for plotting the 13 indexes as labels.\n",
        "index = range(len(selected_variables))"
      ],
      "id": "31402b03",
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {},
      "source": [
        "plt.bar(index, ej_index_detroit, bar_width, label='Detroit') plots a series of bars at the x-positions specified by index and sets the height according to the ej_index_detroit values, with each bar having a width defined by the bar_width variable. The second line, plt.bar(\\[i + bar_width for i in index\\], ej_index_sterlingheights, bar_width, label='Sterling Heights'), plots another series of bars that are shifted to the right by bar_width per each index value to position them next to the Detroit bars. Using only ‘index’ without the adjustment overlaps the bars instead. A similar thing is done to plot the x ticks in the plt.xticks(\\[i + bar_width/2 for i in index\\], selected_variables, rotation=45, ha='right') line, where selected_variables is used as the labels at a 45 degree angle, and ensures that the ticks and labels are positioned between the Detroit and Sterling Heights bars.\n"
      ],
      "id": "5e18a9a8"
    },
    {
      "cell_type": "code",
      "metadata": {},
      "source": [
        "plt.bar(index, ej_index_detroit, bar_width, label='Detroit')\n",
        "\n",
        "plt.bar([i + bar_width for i in index], ej_index_sterlingheights, bar_width, label='Sterling Heights')\n",
        "\n",
        "#Labeling the x-axis and y-axis.\n",
        "plt.xlabel('EJ Index')\n",
        "plt.ylabel('%ile in State')\n",
        "plt.title('Comparison of EJ Index (%ile in State) between Detroit and Sterling Heights')\n",
        "plt.xticks([i + bar_width/2 for i in index], selected_variables, rotation=45, ha='right')\n",
        "\n",
        "#Adding a legend\n",
        "plt.legend()\n",
        "\n",
        "plt.tight_layout()\n",
        "plt.show()"
      ],
      "id": "3ef7ccdf",
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {},
      "source": [
        "What conclusions can you draw from looking at the resulting bar graph?\n",
        "\n",
        "Next, we can start mapping other indexes. Let's increase our study area. Go back to the mapper and zoom out so that you can see the metro area of Detroit, and the surrounding area, up to around Sterling Heights. Delete the polygons that were created by selecting them and choosing “Delete this site.”\n",
        "\n",
        "We are going to do some side by side map comparisons of other indexes. Go to tools and select ‘Side by Side Map Comparisons.’ That will take you to a different web page which should show a left and right hand side of the map. At the top, on the right hand side, click on “Map Data”. A tab similar to the mapper should appear. Select Threshold map for the theme of the map, choose EJ Indexes as type of Index, and Compare to US for location. Keep the lower and upper bounds the same. Select ‘Particulate Matter 2.5’. Then click on ‘Update Map’.\n",
        "\n",
        "![](docs/m203-ej-images/image34.png){fig-align=\"center\"}\n",
        "\n",
        "For the right hand side, select ‘More’ for the theme of the map, choose Health Disparities as the service, and Asthma as the Layer.\n",
        "\n",
        "![](docs/m203-ej-images/image32.png){fig-align=\"center\"}\n",
        "\n",
        "The result should look similar to this:\n",
        "\n",
        "![](docs/m203-ej-images/image29.png){fig-align=\"center\"}\n",
        "\n",
        "Expand the Legend to further understand the data. Do you see a correlation between the data? Play around with the comparisons. Try choosing People of Color while keeping asthma.\n",
        "\n",
        "You can also create reports based on boundaries. Under Tools, select Boundaries. The Map Contents will show you several types of boundaries you can add to your maps, like zip codes, counties, etc. You can generate reports based on these boundaries, and boundaries you want to import, as well.\n",
        "\n",
        "We can add a shapefile to the Mapper. Head on to the Wayne County datapage \\[https://www.waynecounty.com/departments/technology/gis-data.aspx\\]. The page has different datasets under the label “GIS Data”, including “Commission Districts”, “Municipal Boundaries”, etc. Please click on the one labeled “School Districts”. A zip file will immediately download. Once it is done downloading, unzip the contents into a folder.\n",
        "\n",
        "Next, go to Tools, select Add Shapefile, and then select Add File.\n",
        "\n",
        "![](docs/m203-ej-images/image6.png){fig-align=\"center\"}\n",
        "\n",
        "Now, you can generate reports for each school district. What conclusions do you think you can come to by comparing school districts that overlay high percentiles of certain indexes or indicators, in Detroit and other counties?\n",
        "\n",
        "![](docs/m203-ej-images/image21.png){fig-align=\"center\"}\n",
        "\n",
        "Feel free to explore EJScreen more. There are many additional analyses you can perform with this helpful tool. If you ever want to come back to your work, you can save it as a session. In Tools, select Save Session, and add a name to it. Select save, and then under Actions, select ‘Save to File.’ Your session will be downloaded as a .json, which can then be added back in a later time by selecting ‘Load from file’ and choosing the .json file from your folder. For more information, please check out the EJScreen User Guide: https://ejscreen.epa.gov/mapper/help/ejscreen_help.pdf\n",
        "\n",
        "## Conclusion\n",
        "\n",
        "Congratulations! Now you should be able to:\n",
        "\n",
        "-   Map different types of indices and understand percentiles on the Map Contents.\n",
        "-   Interpret the results of the Community Report and understand demographic and indices information for a specific area.\n",
        "-   Generate charts and download the charts as an excel file to use for simple data analysis.\n",
        "-   Create a side by side comparison of health disparities and socioeconomic indexes.\n",
        "-   Upload a shapefile to the EJScreen Mapper.\n",
        "\n",
        "## Continue to Lesson 4\n",
        "\n",
        "[Lesson 3](https://ciesin-geospatial.github.io/TOPSTSCHOOL-air-quality/m204-grdiv1-pm25.html){.btn .btn-primary .btn role=\"button\"}"
      ],
      "id": "13189553"
    }
  ],
  "metadata": {
    "kernelspec": {
      "name": "python3",
      "language": "python",
      "display_name": "Python 3 (ipykernel)"
    }
  },
  "nbformat": 4,
  "nbformat_minor": 5
}