us_diocese_mapper_V5.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# US Diocese Mapper\n",
    "\n",
    "By Kenneth Burchfiel\n",
    "\n",
    "Released under the MIT License\n",
    "\n",
    "(The county_shapes.csv, county_shapes.geojson, state_fips_codes.csv, and counties_by_diocese.csv files are released into in the public domain.)\n",
    "\n",
    "## Note:\n",
    "This is still a work in progress, as I need to update the diocesan boundaries to account for situations where one county lies in multiple dioceses. (Until then, the map will contain a number of inaccuracies.) I also plan to add in cathedral locations and diocese names.\n",
    "\n",
    "I also plan to add in more documentation/comments."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Preliminary steps\n",
    "\n",
    "First, I downloaded county shapefiles from the US Census bureau here: https://www.census.gov/cgi-bin/geo/shapefiles/index.php?year=2021&layergroup=Counties+%28and+equivalent%29\n",
    "\n",
    "See this note regarding use of the shapefiles: https://www2.census.gov/geo/pdfs/maps-data/data/tiger/tgrshp2021/TGRSHP2021_TechDoc_Ch1.pdf\n",
    "\n",
    "Once you download the shapefiles to your computer, extract them using an unzipping utility. The shapefile document within this unzipped folder ends in .shp; for the 2021 county data, the file name is tl_2021_us_county.shp. \n",
    "\n",
    "It appears that the other files within this folder may also be used by Geopandas in the creation of GeoDataFrames, so I recommend accessing the file within the unzipped folder rather than copying it into your project folder.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import time\n",
    "start_time = time.time()\n",
    "import geopandas\n",
    "import pandas as pd\n",
    "import folium\n",
    "import branca\n",
    "import numpy as np\n",
    "from selenium import webdriver\n",
    "import PIL.Image\n",
    "import IPython"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "generate_county_list = False"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "if generate_county_list == True:\n",
    "\n",
    "    shapefile_path = \\\n",
    "        r'C:\\Users\\kburc\\Downloads\\tl_2021_us_county_download\\tl_2021_us_county.shp'\n",
    "    state_fips_codes = pd.read_html(\n",
    "        'https://www.nrcs.usda.gov/wps/portal/nrcs/detail/national/technical/nra/nri/results/?cid=nrcs143_013696')[0]\n",
    "    state_fips_codes = state_fips_codes.iloc[0:55].copy()\n",
    "    state_fips_codes['FIPS'] = state_fips_codes['FIPS'].astype('int')\n",
    "    state_fips_codes.rename(columns={'Name':'State_Name',\n",
    "    'Postal Code':'State_Code'},inplace=True) \n",
    "    state_fips_codes.loc[55] = ['District of Columbia', 'DC', 11] # DC wasn't \n",
    "    # present in the original webpage, so I added it in. I located its FIPS code\n",
    "    # within the county shapefile dataset.\n",
    "    state_fips_codes.to_csv('state_fips_codes.csv')\n",
    "\n",
    "    county_shapes = geopandas.read_file(shapefile_path)\n",
    "    county_shapes['geometry'] = county_shapes.simplify(tolerance = 0.005) \n",
    "    county_shapes.drop(['CLASSFP', 'MTFCC', 'CSAFP', 'CBSAFP', 'METDIVFP', \n",
    "    'FUNCSTAT', 'ALAND', 'AWATER', 'INTPTLAT', 'INTPTLON'], \n",
    "    axis = 1, inplace = True)\n",
    "    county_shapes['STATEFP'] = county_shapes['STATEFP'].astype('int')\n",
    "\n",
    "    merged_county_shapes = county_shapes.merge(state_fips_codes, \n",
    "    left_on = 'STATEFP', right_on = 'FIPS', how = 'left')\n",
    "    merged_county_shapes['county_state'] = merged_county_shapes[\n",
    "        'NAMELSAD'] + ', ' + merged_county_shapes['State_Name']\n",
    "\n",
    "    merged_county_shapes.to_csv('county_shapes.csv')\n",
    "    merged_county_shapes.to_file('county_shapes.geojson',\n",
    "        driver = 'GeoJSON') \n",
    "    # The above line exports the GeoDataFrame created by prepare_zip_table \n",
    "    # so that it can be imported back into the program, which takes less time\n",
    "    # than does recreating the GeoDataFrame.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "merged_county_shapes_test = geopandas.read_file('county_shapes.geojson')\n",
    "merged_county_shapes_test"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "merged_county_shapes_test.set_index('county_state', inplace = True)\n",
    "merged_county_shapes_test"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "merged_county_shapes_test[\n",
    "    'county_state_column'] = merged_county_shapes_test.index\n",
    "merged_county_shapes_test.to_file('county_shapes_test.geojson',\n",
    "    driver = 'GeoJSON') \n",
    "merged_county_shapes_test\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "merged_county_shapes_test['Diocese'] = 'Not_Set'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "condlist = [\n",
    "merged_county_shapes_test['State_Name'] == 'Alabama',\n",
    "merged_county_shapes_test['State_Name'] == 'Alaska',\n",
    "merged_county_shapes_test['State_Name'] == 'Arizona',\n",
    "merged_county_shapes_test['State_Name'] == 'Arkansas',\n",
    "merged_county_shapes_test['State_Name'] == 'California',\n",
    "merged_county_shapes_test['State_Name'] == 'Colorado',\n",
    "merged_county_shapes_test['State_Name'] == 'Connecticut',\n",
    "merged_county_shapes_test['State_Name'] == 'Delaware',\n",
    "merged_county_shapes_test['State_Name'] == 'Florida',\n",
    "merged_county_shapes_test['State_Name'] == 'Georgia',\n",
    "merged_county_shapes_test['State_Name'] == 'Hawaii',\n",
    "merged_county_shapes_test['State_Name'] == 'Idaho',\n",
    "merged_county_shapes_test['State_Name'] == 'Illinois',\n",
    "merged_county_shapes_test['State_Name'] == 'Indiana',\n",
    "merged_county_shapes_test['State_Name'] == 'Iowa',\n",
    "merged_county_shapes_test['State_Name'] == 'Kansas',\n",
    "merged_county_shapes_test['State_Name'] == 'Kentucky',\n",
    "merged_county_shapes_test['State_Name'] == 'Louisiana',\n",
    "merged_county_shapes_test['State_Name'] == 'Maine',\n",
    "merged_county_shapes_test['State_Name'] == 'Maryland',\n",
    "merged_county_shapes_test['State_Name'] == 'Massachusetts',\n",
    "merged_county_shapes_test['State_Name'] == 'Michigan',\n",
    "merged_county_shapes_test['State_Name'] == 'Minnesota',\n",
    "merged_county_shapes_test['State_Name'] == 'Mississippi',\n",
    "merged_county_shapes_test['State_Name'] == 'Missouri',\n",
    "merged_county_shapes_test['State_Name'] == 'Montana',\n",
    "merged_county_shapes_test['State_Name'] == 'Nebraska',\n",
    "merged_county_shapes_test['State_Name'] == 'Nevada',\n",
    "merged_county_shapes_test['State_Name'] == 'New Hampshire',\n",
    "merged_county_shapes_test['State_Name'] == 'New Jersey',\n",
    "merged_county_shapes_test['State_Name'] == 'New Mexico',\n",
    "merged_county_shapes_test['State_Name'] == 'New York',\n",
    "merged_county_shapes_test['State_Name'] == 'North Carolina',\n",
    "merged_county_shapes_test['State_Name'] == 'North Dakota',\n",
    "merged_county_shapes_test['State_Name'] == 'Ohio',\n",
    "merged_county_shapes_test['State_Name'] == 'Oklahoma',\n",
    "merged_county_shapes_test['State_Name'] == 'Oregon',\n",
    "merged_county_shapes_test['State_Name'] == 'Pennsylvania',\n",
    "merged_county_shapes_test['State_Name'] == 'Rhode Island',\n",
    "merged_county_shapes_test['State_Name'] == 'South Carolina',\n",
    "merged_county_shapes_test['State_Name'] == 'South Dakota',\n",
    "merged_county_shapes_test['State_Name'] == 'Tennessee',\n",
    "merged_county_shapes_test['State_Name'] == 'Texas',\n",
    "merged_county_shapes_test['State_Name'] == 'Utah',\n",
    "merged_county_shapes_test['State_Name'] == 'Vermont',\n",
    "merged_county_shapes_test['State_Name'] == 'Virginia',\n",
    "merged_county_shapes_test['State_Name'] == 'Washington',\n",
    "merged_county_shapes_test['State_Name'] == 'West Virginia',\n",
    "merged_county_shapes_test['State_Name'] == 'Wisconsin',\n",
    "merged_county_shapes_test['State_Name'] == 'Wyoming',\n",
    "merged_county_shapes_test['State_Name'] == 'District of Columbia',\n",
    "merged_county_shapes_test['State_Name'] == 'Virgin Islands'\n",
    "]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "choicelist = [\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Alabama',\n",
    "'Mobile',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Alaska',\n",
    "'Fairbanks',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Arizona',\n",
    "'Phoenix',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Arkansas',\n",
    "'Little Rock', # Statewide diocese\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'California',\n",
    "'Los Angeles', \n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Colorado',\n",
    "'Pueblo',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Connecticut',\n",
    "'Norwich', \n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Delaware',\n",
    "'Wilmington', # Statewide diocese\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Florida',\n",
    "'Orlando', \n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Georgia',\n",
    "'Savannah',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Hawaii',\n",
    "'Honolulu', # Statewide diocese\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Idaho',\n",
    "'Boise', #Statewide diocese\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Illinois',\n",
    "'Peoria',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Indiana',\n",
    "'Indianapolis', \n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Iowa',\n",
    "'Dubuque',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Kansas',\n",
    "'Salina',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Kentucky',\n",
    "'Lexington',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Louisiana',\n",
    "'Shreveport',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Maine',\n",
    "'Portland', # Statewide diocese\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Maryland',\n",
    "'Wilmington', \n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Massachusetts',\n",
    "'Boston',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Michigan',\n",
    "'Marquette',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Minnesota',\n",
    "'Duluth',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Mississippi',\n",
    "'Jackson',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Missouri',\n",
    "'Springfield-Cape Girardeau',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Montana',\n",
    "'Helena',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Nebraska',\n",
    "'Grand Island',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Nevada',\n",
    "'Reno',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'New Hampshire',\n",
    "'Manchester', # Statewide diocese\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'New Jersey',\n",
    "'Camden',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'New Mexico',\n",
    "'Santa Fe',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'New York',\n",
    "'Albany',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'North Carolina',\n",
    "'Raleigh',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'North Dakota',\n",
    "'Fargo',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Ohio',\n",
    "'Columbus',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Oklahoma',\n",
    "'Oklahoma City',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Oregon',\n",
    "'Portland',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Pennsylvania',\n",
    "'Erie',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Rhode Island',\n",
    "'Providence', #Statewide diocese\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'South Carolina',\n",
    "'Charleston', # Statewide diocese\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'South Dakota',\n",
    "'Sioux Falls',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Tennessee',\n",
    "'Nashville',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Texas',\n",
    "'Lubbock',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Utah',\n",
    "'Salt Lake City', # Statewide Diocese\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Vermont',\n",
    "'Burlington', # Statewide diocese\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Virginia',\n",
    "'Richmond',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Washington',\n",
    "'Seattle',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'West Virginia',\n",
    "'Wheeling-Charleston', # Statewide diocese\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Wisconsin',\n",
    "'La Crosse',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Wyoming'],\n",
    "'Cheyenne', # Statewide diocese\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'District of Columbia',\n",
    "'Washington',\n",
    "\n",
    "# merged_county_shapes_test['State_Name'] == 'Virgin Islands'\n",
    "'Saint Thomas'\n",
    "]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "merged_county_shapes_test['Diocese'] = np.select(condlist, choicelist, merged_county_shapes_test['Diocese'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "merged_county_shapes_test = merged_county_shapes_test.query(\n",
    "    \"State_Name != 'Puerto Rico'\").copy() # Puerto Rico belongs to a different\n",
    "    # episcopal conference, so it won't be included in this dataset."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "initial_version_of_merged_county_shapes_test = \\\n",
    "merged_county_shapes_test.copy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "def add_dioceses(df, county_list, state, diocese):\n",
    "    # Adds the diocese for all counties in a given state. If a diocese spans\n",
    "    # multiple stats, this function will need to be called once for each\n",
    "    # state.\n",
    "    # If a given 'county' variable does not have any spaces, \n",
    "    # the loop below will add ' County' to that string before \n",
    "    # adding in the state. Otherwise, the code assumes that the 'county' is\n",
    "    # either a non-county (e.g. a municipality, city, census area, parish,\n",
    "    # or borough), and will only add the state to this string.\n",
    "    # Therefore, in order to run this function successfully, enter the \n",
    "    # area designation (county, city, etc.) after all non-counties and\n",
    "    # counties with more than one word when building the county list.\n",
    "    # However, you don't need to add in 'County' for counties with only\n",
    "    # one word.\n",
    "    for county in county_list:\n",
    "        if ' ' not in county:\n",
    "            county_string = county + ' County, ' + state\n",
    "        else:\n",
    "            county_string = county + ', ' + state\n",
    "        if county_string not in df.index:\n",
    "            raise ValueError(county_string, \"not found in index\")\n",
    "        df.at[county_string, 'Diocese'] = diocese\n",
    "    return df\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "List of dioceses:\n",
    "https://en.wikipedia.org/wiki/List_of_Catholic_dioceses_in_the_United_States\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Diocese of Anchorage-Juneau:\n",
    "# Source: Wikipedia province map\n",
    "\n",
    "anchorage_juneau_county_list = [\n",
    "'Aleutians West Census Area',\n",
    "'Aleutians East Borough', 'Lake and Peninsula Borough',\n",
    "'Bristol Bay Borough', 'Dillingham Census Area', 'Kodiak Island Borough',\n",
    "'Kenai Peninsula Borough', 'Chugach Census Area', 'Yakutat City and Borough',\n",
    "'Aleutians West Census Area', 'Hoonah-Angoon Census Area',\n",
    "'Petersburg Borough', 'Wrangell City and Borough', 'Ketchikan Gateway Borough',\n",
    "'Prince of Wales-Hyder Census Area', 'Sitka City and Borough',\n",
    "'Copper River Census Area', 'Matanuska-Susitna Borough',\n",
    "'Anchorage Municipality', 'Juneau City and Borough', 'Skagway Municipality',\n",
    "'Haines Borough'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "anchorage_juneau_county_list, state = 'Alaska', diocese = 'Anchorage-Juneau')\n",
    "\n",
    "# Archdiocese of Atlanta (Source: Diocesan website:\n",
    "# https://archatl.com/about/about-the-archdiocese-of-atlanta/)\n",
    "# This site noted the counties within the Archdiocese of Atlanta that border\n",
    "# the Diocese of Savannah. I entered in those counties first, then added in\n",
    "# all other counties to the north of those.\n",
    "\n",
    "atlanta_county_list = [\n",
    "'Gilmer', 'Carroll', 'Douglas', 'Clayton', 'Spalding', 'Henry', 'Rockdale',\n",
    "'Oconee', 'Barrow', 'Jackson', 'Madison', 'Elbert', 'Rabun', 'Towns',\n",
    "'Union', 'Fannin', 'Murray', 'Catoosa', 'Whitfield', 'Dade', 'Walker',\n",
    "'Chattooga', 'Gordon', 'Pickens', 'Dawson', 'Lumpkin', 'White', 'Habersham',\n",
    "'Stephens', 'Hart', 'Franklin', 'Banks', 'Hall', 'Forsyth', 'Fulton',\n",
    "'Cherokee', 'Bartow', 'Floyd', 'Polk', 'Haralson', 'Paulding', 'Cobb',\n",
    "'DeKalb', 'Gwinnett', 'Walton', 'Clarke', 'Oglethorpe', 'Coweta', 'Wilkes',\n",
    "'Coweta', 'Taliaferro', 'Greene', 'Morgan', 'Coweta', 'Newton', 'Butts', \n",
    "'Lamar', 'Pike', 'Coweta', 'Fayette', 'Coweta', 'Heard', 'Troup', 'Meriwether',\n",
    "'Upson', 'Monroe', 'Jasper', 'Putnam', 'Baldwin', 'Hancock', 'Warren',\n",
    "'McDuffie', 'Lincoln'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "atlanta_county_list, state = 'Georgia', diocese = 'Atlanta')\n",
    "\n",
    "# # Diocese of Charlotte:\n",
    "# Source: https://web.archive.org/web/20110504140939/http://www.charlottediocese.org/atlasofthediocese.html (See also:\n",
    "# https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Charlotte#/media/File:Diocese_of_Charlotte.jpg)\n",
    "\n",
    "charlotte_county_list = [\n",
    "'Richmond', 'Montgomery', 'Randolph', 'Guilford', 'Rockingham', 'Mecklenburg', \n",
    "'Davie', 'Yadkin', 'Surry', 'Alleghany', 'Wilkes', 'Iredell', 'Catawba', \n",
    "'Lincoln', 'Gaston', 'Cleveland', 'Stokes', 'Forsyth', 'Davidson', 'Rowan', \n",
    "'Stanly', 'Anson', 'Union', 'Cabarrus', 'Cherokee', 'Graham', 'Swain', \n",
    "'Haywood', 'Madison', 'Ashe', 'Watauga', 'Caldwell', 'Alexander', 'Burke', \n",
    "'McDowell', 'Rutherford', 'Polk', 'Henderson', 'Buncombe', 'Yancey', 'Mitchell', \n",
    "'Clay', 'Macon', 'Jackson', 'Transylvania', 'Avery'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "charlotte_county_list, state = 'North Carolina', diocese = 'Charlotte')\n",
    "\n",
    "# Archdiocese of Baltimore:\n",
    "# https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_Baltimore\n",
    "\n",
    "baltimore_county_list = [\n",
    "'Anne Arundel County', 'Baltimore', 'Carroll', 'Frederick', 'Garrett', \n",
    "'Harford', 'Howard', 'Washington', 'Baltimore city','Allegany'\n",
    "] \n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "baltimore_county_list, state = 'Maryland', diocese = 'Baltimore')\n",
    "\n",
    "# Diocese of Arlington:\n",
    "# https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Arlington\n",
    "\n",
    "arlington_county_list = [\n",
    "'Arlington', 'Clarke', 'Culpeper', 'Fairfax', 'Fauquier', 'Frederick', \n",
    "'King George County', 'Lancaster', 'Loudoun', 'Madison', 'Northumberland', \n",
    "'Orange', 'Page', 'Prince William County', 'Rappahannock', 'Richmond',\n",
    "'Shenandoah', 'Spotsylvania', 'Stafford', 'Warren', 'Westmoreland', \n",
    "'Alexandria city', 'Falls Church city', 'Fairfax city', 'Fredericksburg city',\n",
    "'Manassas city', 'Manassas Park city', 'Winchester city'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "arlington_county_list, state = 'Virginia', diocese = 'Arlington')\n",
    "\n",
    "# Archdiocese of Washington:\n",
    "# https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_Washington\n",
    "# Note that the Archdiocese of Washington includes a number of counties\n",
    "# in Maryland, which will be added here. The District of Columbia is already\n",
    "# listed as belonging to the Diocese of Washington.\n",
    "washington_county_list = [\n",
    "'Calvert', 'Charles', 'Montgomery', \n",
    "\"Prince George's County\", \"St. Mary's County\"\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "washington_county_list, state = 'Maryland', diocese = 'Washington')\n",
    "\n",
    "\n",
    "# # Diocese of Fall River:\n",
    "# # https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Fall_River\n",
    "\n",
    "fall_river_county_list = [\n",
    "'Barnstable', 'Bristol', 'Dukes', 'Nantucket'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "fall_river_county_list, state = 'Massachusetts', diocese = 'Fall River')\n",
    "\n",
    "# Diocese of Springfield in Massachusetts:\n",
    "# https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Springfield_in_Massachusetts\n",
    "\n",
    "springfield_in_massachusetts_county_list = [\n",
    "'Berkshire', 'Franklin', 'Hampshire', 'Hampden'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "springfield_in_massachusetts_county_list, state = 'Massachusetts', diocese = \n",
    "'Springfield in Massachusetts')\n",
    "\n",
    "# Diocese of Worcester: \n",
    "# https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Worcester\n",
    "\n",
    "merged_county_shapes_test.at[\n",
    "\"Worcester County, Massachusetts\", 'Diocese'] = 'Worcester' \n",
    "# The Diocese of Worcester includes only one county, so it can\n",
    "# be added directly to the DataFrame.\n",
    "\n",
    "# Archdiocese of Chicago: \n",
    "# https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_Chicago\n",
    "\n",
    "chicago_county_list = [\n",
    "'Cook', 'Lake'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "chicago_county_list, state = 'Illinois', diocese = \n",
    "'Chicago')\n",
    "\n",
    "# # Diocese of Belleville:\n",
    "# # I couldn't find a list of counties on the diocese's website\n",
    "# # or on the Diocese's Wikipedia page, so I based the county\n",
    "# # list off this map:\n",
    "# # https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Belleville#/media/File:Diocese_of_Belleville_map_1.png\n",
    "\n",
    "belleville_county_list = [\n",
    "\"St. Clair County\", \"Clinton\", \"Marion\", \"Clay\", \"Richland\", \"Lawrence\", \n",
    "\"Wabash\", \"Edwards\", \"Wayne\", \"Monroe\", \"Washington\", \"Jefferson\", \"Hamilton\", \n",
    "\"White\", \"Gallatin\", \"Saline\", \"Williamson\", \"Jackson\", \"Union\", \"Johnson\", \n",
    "\"Pope\", \"Massac\", \"Pulaski\", \"Alexander\", \"Randolph\", \"Perry\", \"Franklin\", \n",
    "\"Hardin\"\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "belleville_county_list, state = 'Illinois', diocese = \n",
    "'Belleville')\n",
    "\n",
    "# # Diocese of Joliet:\n",
    "# # https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Joliet_in_Illinois\n",
    "\n",
    "joliet_county_list = [\n",
    "\"DuPage\", \"Ford\", \"Grundy\", \"Iroquois\", \"Kankakee\", \"Kendall\", \"Will\"\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "joliet_county_list, state = 'Illinois', diocese = \n",
    "'Joliet')\n",
    "\n",
    "# Diocese of Rockford:\n",
    "# https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Rockford\n",
    "\n",
    "rockford_county_list = [\n",
    "\"Boone\", \"Carroll\", \"DeKalb\", \"Jo Daviess County\", \"Kane\", \"Lee\", \"McHenry\",\n",
    "\"Ogle\", \"Stephenson\", \"Whiteside\", \"Winnebago\"\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "rockford_county_list, state = 'Illinois', diocese = \n",
    "'Rockford')\n",
    "\n",
    "springfield_in_illinois_county_list = [\n",
    "\"Adams\", \"Bond\", \"Brown\", \"Calhoun\", \"Cass\", \"Christian\", \"Clark\", \"Coles\", \n",
    "\"Crawford\", \"Cumberland\", \"Douglas\", \"Edgar\", \"Effingham\", \"Fayette\", \"Greene\", \n",
    "\"Jasper\", \"Jersey\", \"Macon\", \"Macoupin\", \"Madison\", \"Menard\", \"Moultrie\", \n",
    "\"Montgomery\", \"Morgan\", \"Pike\", \"Sangamon\", \"Scott\", \"Shelby\"\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "springfield_in_illinois_county_list, state = 'Illinois', diocese = \n",
    "'Springfield in Illinois')\n",
    "\n",
    "# Archdiocese of Cincinnati:\n",
    "# Source: https://catholicaoc.org/about/geography\n",
    "\n",
    "cincinnati_county_list = [\n",
    "'Mercer', 'Auglaize', 'Logan', 'Darke', 'Shelby', 'Miami', 'Clark', 'Preble', \n",
    "'Montgomery', 'Greene', 'Butler', 'Warren', 'Clinton', 'Hamilton', 'Clermont', \n",
    "'Brown', 'Highland', 'Adams', 'Champaign'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "cincinnati_county_list, state = 'Ohio', diocese = \n",
    "'Cincinnati')\n",
    "\n",
    "# Diocese of Cleveland: \n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Cleveland\n",
    "cleveland_county_list = [\n",
    "'Ashland', 'Cuyahoga', 'Geauga', 'Lake', 'Lorain',\n",
    "'Medina', 'Summit', 'Wayne'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "cleveland_county_list, state = 'Ohio', diocese = \n",
    "'Cleveland')\n",
    "\n",
    "# Diocese of Steubenville:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Steubenville\n",
    "steubenville_county_list = [\n",
    "'Athens', 'Belmont', 'Carroll', 'Gallia', 'Guernsey', 'Harrison', 'Jefferson', \n",
    "'Lawrence', 'Meigs', 'Morgan', 'Monroe', 'Noble', 'Washington'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "steubenville_county_list, state = 'Ohio', diocese = \n",
    "'Steubenville')\n",
    "\n",
    "# Diocese of Toledo:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Toledo\n",
    "\n",
    "toledo_county_list = [\n",
    "'Allen', 'Crawford', 'Defiance', 'Erie', 'Fulton',\n",
    "'Hancock', 'Henry', 'Huron', 'Lucas', 'Ottawa', 'Paulding', 'Putnam', \n",
    "'Richland', 'Sandusky', 'Seneca', 'Van Wert County', 'Williams', 'Wood',\n",
    "'Wyandot'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "toledo_county_list, state = 'Ohio', diocese = \n",
    "'Toledo')\n",
    "\n",
    "# Diocese of Youngstown:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Youngstown\n",
    "\n",
    "youngstown_county_list = [\n",
    "'Mahoning', 'Trumbull', 'Columbiana', 'Stark', 'Portage', 'Ashtabula'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "youngstown_county_list, state = 'Ohio', diocese = \n",
    "'Youngstown')\n",
    "\n",
    "# Archdiocese of Denver:\n",
    "# Source: https://archden.org/wp-content/uploads/2021/11/2021_media_kit_updated_11-3-21.pdf\n",
    "# This source incorrectly spells Broomfield County as \n",
    "# 'Bromfield.' \n",
    "\n",
    "denver_county_list = [\n",
    "'Adams', 'Arapahoe', 'Boulder', 'Broomfield', \n",
    "'Clear Creek County', 'Denver', 'Eagle', 'Garfield',\n",
    "'Gilpin', 'Grand', 'Jackson', 'Jefferson', 'Larimer',\n",
    "'Logan', 'Moffat', 'Morgan', 'Phillips', 'Pitkin',\n",
    "'Rio Blanco County', 'Routt', 'Sedgwick', 'Summit',\n",
    "'Washington', 'Weld', 'Yuma'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "denver_county_list, state = 'Colorado', diocese = \n",
    "'Denver')\n",
    "\n",
    "\n",
    "# # Diocese of Colorado Springs:\n",
    "# # Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Colorado_Springs\n",
    "\n",
    "colorado_springs_county_list = [\n",
    "'Chaffee', 'Lake', 'Park', 'Teller', 'Douglas',\n",
    "'El Paso County', 'Elbert', 'Lincoln', \n",
    "'Kit Carson County', 'Cheyenne'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "colorado_springs_county_list, state = 'Colorado', diocese = \n",
    "'Colorado Springs')\n",
    "\n",
    "# Archdiocese of Detroit:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_Detroit\n",
    "\n",
    "detroit_county_list = [\n",
    "'Lapeer', 'Macomb', 'Monroe', 'Oakland', 'St. Clair County', 'Wayne'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "detroit_county_list, state = 'Michigan', diocese = \n",
    "'Detroit')\n",
    "\n",
    "\n",
    "# Diocese of Gaylord:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Gaylord\n",
    "\n",
    "gaylord_county_list = [\n",
    "'Charlevoix', 'Emmet', 'Cheboygan', 'Presque Isle County', \n",
    "'Leelanau', 'Antrim', 'Otsego', 'Montmorency', 'Alpena', 'Benzie', \n",
    "'Grand Traverse County', 'Kalkaska', 'Crawford', 'Oscoda', 'Alcona',\n",
    "'Manistee', 'Wexford', 'Missaukee', 'Roscommon', 'Ogemaw', 'Iosco'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "gaylord_county_list, state = 'Michigan', diocese = \n",
    "'Gaylord')\n",
    "\n",
    "# Diocese of Grand Rapids:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Grand_Rapids\n",
    "\n",
    "grand_rapids_county_list = [\n",
    "'Ottawa', 'Kent', 'Ionia', 'Muskegon', 'Newaygo', 'Oceana', \n",
    "'Montcalm', 'Mecosta', 'Lake', 'Mason', 'Osceola'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "grand_rapids_county_list, state = 'Michigan', diocese = \n",
    "'Grand Rapids')\n",
    "\n",
    "\n",
    "# Diocese of Kalamazoo:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Kalamazoo\n",
    "kalamazoo_county_list = [\n",
    "'Allegan', 'Barry', 'Van Buren County', 'Kalamazoo', 'Calhoun', \n",
    "'Berrien', 'Cass', 'St. Joseph County', 'Branch'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "kalamazoo_county_list, state = 'Michigan', diocese = \n",
    "'Kalamazoo')\n",
    "\n",
    "# Diocese of Lansing:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Lansing\n",
    "\n",
    "lansing_county_list = [\n",
    "'Clinton', 'Eaton', 'Genesee', 'Hillsdale', 'Ingham', 'Jackson', \n",
    "'Lenawee', 'Livingston', 'Shiawassee', 'Washtenaw'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "lansing_county_list, state = 'Michigan', diocese = \n",
    "'Lansing')\n",
    "\n",
    "\n",
    "# Diocese of Saginaw:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Saginaw\n",
    "saginaw_county_list = [\n",
    "'Arenac', 'Bay', 'Clare', 'Gladwin', 'Gratiot', 'Huron', 'Isabella',\n",
    "'Midland', 'Saginaw', 'Sanilac', 'Tuscola']\n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "saginaw_county_list, state = 'Michigan', diocese = \n",
    "'Saginaw')\n",
    "\n",
    "# Diocese of Davenport:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Davenport\n",
    "# See also:\n",
    "# https://d2y1pz2y630308.cloudfront.net/13543/documents/2020/8/Parishes%20Map%207-1-20.pdf\n",
    "davenport_county_list = [\n",
    "'Jasper', 'Poweshiek', 'Iowa', 'Johnson', 'Cedar', 'Clinton', 'Marion',\n",
    "'Monroe', 'Appanoose', 'Mahaska', 'Wapello', 'Davis', 'Keokuk', \n",
    "'Washington', 'Muscatine', 'Scott', 'Louisa', 'Des Moines County', \n",
    "'Henry', 'Lee', 'Van Buren County', 'Jefferson']\n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "davenport_county_list, state = 'Iowa', diocese = \n",
    "'Davenport')\n",
    "\n",
    "# Diocese of Sioux CIty:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Sioux_City#/media/File:Diocese_of_Sioux_City.jpg\n",
    "\n",
    "sioux_city_county_list = [\n",
    "'Boone', 'Greene', 'Carroll', 'Crawford', 'Monona', 'Webster', 'Humboldt', \n",
    "'Kossuth', 'Dickinson', 'Lyon', 'Sioux', 'Plymouth', 'Woodbury', 'Osceola',\n",
    "'Emmet', \"O'Brien\", 'Clay', 'Palo Alto County', 'Cherokee', \n",
    "'Buena Vista County', 'Pocahontas', 'Ida', 'Sac', 'Calhoun'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "sioux_city_county_list, state = 'Iowa', diocese = \n",
    "'Sioux City')\n",
    "\n",
    "# Diocese of Des Moines:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Des_Moines#/media/File:Diocese_of_Des_Moines.jpg\n",
    "\n",
    "des_moines_county_list = [\n",
    "'Harrison', 'Shelby', 'Audubon', 'Guthrie', 'Dallas', 'Polk',\n",
    "'Pottawattamie', 'Cass', 'Adair', 'Madison', 'Warren', 'Mills',\n",
    "'Montgomery', 'Adams', 'Union', 'Clarke', 'Lucas', 'Fremont',\n",
    "'Page', 'Taylor', 'Ringgold', 'Decatur', 'Wayne'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "des_moines_county_list, state = 'Iowa', diocese = \n",
    "'Des Moines')\n",
    "\n",
    "# Archdiocese of Galveston-Houston:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_Galveston%E2%80%93Houston\n",
    "galveston_houston_county_list = [\n",
    "'Galveston', 'Harris', 'Austin', 'Brazoria', 'Fort Bend County', 'Grimes',\n",
    "'Montgomery', 'San Jacinto County', 'Walker', 'Waller'] \n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "galveston_houston_county_list, state = 'Texas', diocese = \n",
    "'Galveston-Houston')\n",
    "\n",
    "\n",
    "# Diocese of Austin:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Austin\n",
    "austin_county_list = [\n",
    "'Bastrop', 'Bell', 'Blanco', 'Brazos', 'Burleson', 'Burnet', 'Caldwell',\n",
    "'Coryell', 'Falls', 'Hamilton', 'Hays', 'Lampasas', 'Lee', 'Limestone',\n",
    "'Llano', 'Mason', 'McLennan', 'Milam', 'Mills', 'Robertson', 'San Saba County',\n",
    "'Travis', 'Washington', 'Williamson'\n",
    "]\n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "austin_county_list, state = 'Texas', diocese = \n",
    "'Austin')\n",
    "\n",
    "# Diocese of Beaumont:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Beaumont\n",
    "\n",
    "beaumont_county_list = [\n",
    "'Chambers', 'Hardin', 'Jasper', 'Jefferson', 'Liberty', 'Newton', 'Orange',\n",
    "'Polk', 'Tyler'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "beaumont_county_list, state = 'Texas', diocese = \n",
    "'Beaumont')\n",
    "\n",
    "# Diocese of Brownsville:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Brownsville\n",
    "brownsville_county_list = [\n",
    "'Starr', 'Willacy', 'Hidalgo', 'Cameron'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "brownsville_county_list, state = 'Texas', diocese = \n",
    "'Brownsville')\n",
    "\n",
    "# Diocese of Corpus Christi:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Corpus_Christi\n",
    "\n",
    "corpus_christi_county_list = [\n",
    "'Aransas', 'Bee', 'Brooks', 'Duval', 'Jim Wells County', 'Kenedy', 'Kleberg', \n",
    "'Live Oak County', 'Nueces', 'Refugio', 'San Patricio County'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "corpus_christi_county_list, state = 'Texas', diocese = \n",
    "'Corpus Christi')\n",
    "\n",
    "# Diocese of Tyler:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Tyler\n",
    "\n",
    "tyler_county_list = ['Anderson', 'Angelina', 'Bowie', 'Camp', 'Cass',\n",
    "'Cherokee', 'Delta', 'Franklin', 'Freestone', 'Gregg', 'Harrison',\n",
    "'Henderson', 'Hopkins', 'Houston', 'Lamar', 'Leon', 'Madison', 'Marion', \n",
    "'Morris', 'Nacogdoches', 'Panola', 'Rains', 'Red River County', 'Rusk', \n",
    "'Sabine', 'San Augustine County', 'Shelby', 'Smith', 'Titus', 'Trinity', \n",
    "'Upshur', 'Van Zandt County', 'Wood']\n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "tyler_county_list, state = 'Texas', diocese = \n",
    "'Tyler')\n",
    "\n",
    "# Diocese of Victoria:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Victoria_in_Texas\n",
    "\n",
    "victoria_in_texas_county_list = [\n",
    "'Calhoun', 'Colorado', 'DeWitt', 'Fayette', 'Goliad', 'Jackson', 'Lavaca', \n",
    "'Matagorda', 'Victoria', 'Wharton'\n",
    "] \n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "victoria_in_texas_county_list, state = 'Texas', diocese = \n",
    "'Victoria in Texas')\n",
    "\n",
    "# Archddiocese of San Antonio:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_San_Antonio\n",
    "\n",
    "san_antonio_county_list = [\n",
    "'Val Verde County', 'Edwards', 'Kerr', 'Gillespie', 'Kendall', 'Comal',\n",
    "'Guadalupe', 'Gonzales', 'Uvalde', 'Bandera', 'Real', 'Kinney', 'Medina', 'Bexar', 'Wilson', 'Karnes', 'Frio', 'Atascosa', 'McMullen'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "san_antonio_county_list, state = 'Texas', diocese = \n",
    "'San Antonio')\n",
    "\n",
    "# Diocese of Amarillo:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Amarillo\n",
    "\n",
    "amarillo_county_list = [\n",
    "'Armstrong', 'Briscoe', 'Carson', 'Castro', 'Childress', 'Collingsworth', \n",
    "'Dallam', 'Deaf Smith County', 'Donley', 'Gray', 'Hall', 'Hansford', 'Hartley',\n",
    "'Hemphill', 'Hutchinson', 'Lipscomb', 'Moore', 'Ochiltree', 'Oldham', 'Parmer',\n",
    "'Potter', 'Randall', 'Roberts', 'Sherman', 'Swisher', 'Wheeler'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "amarillo_county_list, state = 'Texas', diocese = \n",
    "'Amarillo')\n",
    "\n",
    "# Diocese of Dallas:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Dallas\n",
    "\n",
    "dallas_county_list = [\n",
    "'Collin', 'Dallas', 'Ellis', 'Fannin', 'Grayson', 'Hunt', 'Kaufman', 'Navarro',\n",
    "'Rockwall'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "dallas_county_list, state = 'Texas', diocese = \n",
    "'Dallas')\n",
    "\n",
    "# Diocese of El Paso:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_El_Paso\n",
    "el_paso_county_list = [\n",
    "'El Paso County', 'Brewster', 'Culberson', 'Hudspeth', 'Jeff Davis County',\n",
    "'Loving', 'Presidio', 'Reeves', 'Ward', 'Winkler'\n",
    "] \n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "el_paso_county_list, state = 'Texas', diocese = \n",
    "'El Paso')\n",
    "\n",
    "# Diocese of Fort Worth:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Fort_Worth\n",
    "\n",
    "fort_worth_county_list = [\n",
    "'Archer', 'Baylor', 'Bosque', 'Clay', 'Comanche', 'Cooke', 'Denton', 'Eastland',\n",
    "'Erath', 'Foard', 'Hardeman', 'Hill', 'Hood', 'Jack', 'Johnson', 'Knox', \n",
    "'Montague', 'Palo Pinto County', 'Parker', 'Shackelford', 'Somervell', \n",
    "'Stephens', 'Tarrant', 'Throckmorton', 'Wichita', 'Wilbarger', 'Wise', 'Young'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "fort_worth_county_list, state = 'Texas', diocese = \n",
    "'Fort Worth')\n",
    "\n",
    "# Diocese of Laredo:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Laredo\n",
    "\n",
    "laredo_county_list = [\n",
    "'Dimmit', 'Jim Hogg County', 'La Salle County', 'Maverick', 'Webb', 'Zapata', \n",
    "'Zavala']\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "laredo_county_list, state = 'Texas', diocese = \n",
    "'Laredo')\n",
    "\n",
    "# Diocese of San Angelo:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_San_Angelo\n",
    "\n",
    "san_angelo_county_list = [\n",
    "'Andrews', 'Brown', 'Callahan', 'Coke', 'Coleman', 'Concho', 'Crane', \n",
    "'Crockett', 'Ector', 'Glasscock', 'Howard', 'Irion', 'Kimble', 'Martin', \n",
    "'McCulloch', 'Menard', 'Midland', 'Mitchell', 'Nolan', 'Pecos', 'Reagan',\n",
    "'Runnels', 'Schleicher', 'Sterling', 'Sutton', 'Taylor', 'Terrell', \n",
    "'Tom Green County', 'Upton'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "san_angelo_county_list, state = 'Texas', diocese = \n",
    "'San Angelo')\n",
    "\n",
    "# Archdiocese of Hartford:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_Hartford\n",
    "\n",
    "hartford_county_list = [\n",
    "'Hartford', 'Litchfield', 'New Haven County'\n",
    "]\n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "hartford_county_list, state = 'Connecticut', diocese = \n",
    "'Hartford')\n",
    "\n",
    "# Diocese of Bridgeport:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Bridgeport\n",
    "\n",
    "merged_county_shapes_test.at[\n",
    "\"Fairfield County, Connecticut\", 'Diocese'] = 'Bridgeport' \n",
    "# The Diocese of Bridgeport includes only one county, so it can\n",
    "# be added directly to the DataFrame.\n",
    "\n",
    "# Diocese of Evansville:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Evansville\n",
    "evansville_county_list = ['Warrick', 'Greene', 'Knox', 'Sullivan', 'Daviess', \n",
    "'Martin', 'Dubois', 'Spencer', 'Gibson', 'Pike', 'Posey', 'Vanderburgh']\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "evansville_county_list, state = 'Indiana', diocese = \n",
    "'Evansville')\n",
    "\n",
    "\n",
    "# Diocese of Fort Wayne-South Bend:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Fort_Wayne%E2%80%93South_Bend\n",
    "\n",
    "fort_wayne_south_bend_county_list = [\n",
    "'Adams', 'Allen', 'DeKalb', 'Elkhart', 'Huntington', 'Kosciusko', 'LaGrange',\n",
    "'Marshall', 'Noble', 'Steuben', 'St. Joseph County', 'Wabash', 'Wells', \n",
    "'Whitley'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "fort_wayne_south_bend_county_list, state = 'Indiana', diocese = \n",
    "'Fort Wayne-South Bend')\n",
    "\n",
    "\n",
    "# Diocese of Gary:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Gary\n",
    "\n",
    "gary_county_list = [\n",
    "'Lake', 'Porter', 'LaPorte', 'Starke'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "gary_county_list, state = 'Indiana', diocese = \n",
    "'Gary')\n",
    "\n",
    "# Diocese of Lafayette in Indiana:\n",
    "# Source: https://dol-in.org/uih-maps\n",
    "\n",
    "lafayette_in_indiana_county_list = [\n",
    "'Newton', 'Jasper', 'Pulaski', 'Fulton', 'Benton', 'White', 'Carroll', 'Cass',\n",
    "'Miami', 'Warren', 'Tippecanoe', 'Clinton', 'Tipton', 'Grant', 'Blackford', \n",
    "'Jay', 'Fountain', 'Montgomery', 'Boone', 'Madison', 'Hamilton', 'Delaware',\n",
    "'Randolph', 'Howard'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "lafayette_in_indiana_county_list, state = 'Indiana', diocese = \n",
    "'Lafayette in Indiana')\n",
    "\n",
    "\n",
    "# Diocese of Kansas City in Kansas:\n",
    "# Source: \n",
    "\n",
    "kansas_city_in_kansas_county_list = [\n",
    "'Anderson', 'Atchison', 'Brown', 'Coffey', 'Doniphan', 'Douglas', 'Franklin',\n",
    "'Jackson', 'Jefferson', 'Johnson', 'Leavenworth', 'Linn', 'Lyon', 'Marshall',\n",
    "'Miami', 'Nemaha', 'Osage', 'Pottawatomie', 'Shawnee', 'Wabaunsee', 'Wyandotte'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "kansas_city_in_kansas_county_list, state = 'Kansas', diocese = \n",
    "'Kansas City in Kansas')\n",
    "\n",
    "# Diocese of Dodge City:\n",
    "# Source: https://www.dcdiocese.org/about/history\n",
    "# Note: This website misspells Comanche as Commanche and Kearny as Kearney.\n",
    "dodge_city_county_list = [\n",
    "'Barber', 'Barton', 'Clark', 'Comanche', 'Edwards', 'Finney', 'Ford', 'Grant',\n",
    "'Gray', 'Greeley', 'Hamilton', 'Haskell', 'Hodgeman', 'Kearny', 'Kiowa', \n",
    "'Lane', 'Meade', 'Morton', 'Ness', 'Pawnee', 'Pratt', 'Rush', 'Scott', 'Seward',\n",
    "'Stafford', 'Stanton', 'Stevens', 'Wichita'\n",
    "]\n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "dodge_city_county_list, state = 'Kansas', diocese = \n",
    "'Dodge City')\n",
    "\n",
    "# Diocese of Wichita:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Wichita\n",
    "# Note: The map on this page shows Ellsworth County as belonging to\n",
    "# the Diocese of Wichita. However, Ellsworth County actually belongs\n",
    "# to the Diocese of Salina: https://salinadiocese.org/wp-content/uploads/2021/08/2021-Salina-Diocese-Map-PARISH-1.pdf\n",
    "\n",
    "wichita_county_list = [\n",
    "'Allen', 'Bourbon', 'Butler', 'Chase', 'Chautauqua', 'Cherokee', 'Cowley',\n",
    "'Crawford', 'Elk', 'Greenwood', 'Harper', 'Harvey', 'Kingman', 'Labette',\n",
    "'Marion', 'McPherson', 'Montgomery', 'Morris', 'Neosho', 'Reno', 'Rice', \n",
    "'Sedgwick', 'Sumner', 'Wilson', 'Woodson'\n",
    "]\n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "wichita_county_list, state = 'Kansas', diocese = \n",
    "'Wichita')\n",
    "\n",
    "\n",
    "# Diocese of Fresno:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Fresno\n",
    "\n",
    "fresno_county_list = [\n",
    "'Fresno', 'Inyo', 'Kern', 'Kings', 'Madera', 'Mariposa', 'Merced', 'Tulare'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "fresno_county_list, state = 'California', diocese = 'Fresno')\n",
    "\n",
    "# Diocese of Monterey in California:\n",
    "monterey_in_california_county_list = [\n",
    "'Monterey', 'San Benito County', 'San Luis Obispo County', 'Santa Cruz County'\n",
    "]\n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "monterey_in_california_county_list, state = 'California', \n",
    "diocese = 'Monterey in California')\n",
    "\n",
    "# Diocese of Orange:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Orange\n",
    "merged_county_shapes_test.at[\n",
    "\"Orange County, California\", 'Diocese'] = 'Orange' \n",
    "\n",
    "# Diocese of San Bernardino:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_San_Bernardino\n",
    "san_bernardino_county_list = ['San Bernardino County', 'Riverside']\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "san_bernardino_county_list, state = 'California', \n",
    "diocese = 'San Bernardino')\n",
    "\n",
    "# Diocese of San Diego:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_San_Diego\n",
    "\n",
    "san_diego_county_list = [\n",
    "'San Diego County', 'Imperial'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "san_diego_county_list, state = 'California', \n",
    "diocese = 'San Diego')\n",
    "\n",
    "\n",
    "# Diocese of San Francisco:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_San_Francisco\n",
    "\n",
    "san_francisco_county_list = [\n",
    "'San Francisco County', 'Marin', 'San Mateo County'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "san_francisco_county_list, state = 'California', \n",
    "diocese = 'San Francisco')\n",
    "\n",
    "# Diocese of Oakland:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Oakland\n",
    "\n",
    "oakland_county_list = [\n",
    "'Alameda', 'Contra Costa County'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "oakland_county_list, state = 'California', \n",
    "diocese = 'Oakland')\n",
    "\n",
    "# Diocese of Sacramento:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Sacramento\n",
    "sacramento_county_list = [\n",
    "'Siskiyou', 'Modoc', 'Trinity', 'Shasta', 'Lassen', 'Tehama', 'Plumas', 'Glenn',\n",
    "'Butte', 'Sierra', 'Colusa', 'Sutter', 'Yuba', 'Nevada', 'Yolo', 'Placer', \n",
    "'Solano', 'Sacramento', 'El Dorado County', 'Amador'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "sacramento_county_list, state = 'California', \n",
    "diocese = 'Sacramento')\n",
    "\n",
    "# Diocese of San Jose in California:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_San_Jose_in_California\n",
    "\n",
    "merged_county_shapes_test.at[\n",
    "\"Santa Clara County, California\", 'Diocese'] = 'San Jose in California' \n",
    "\n",
    "# Diocese of Santa Rosa in California:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Santa_Rosa_in_California\n",
    "\n",
    "santa_rosa_in_california_county_list = [\n",
    "'Del Norte County', 'Humboldt', 'Lake', 'Mendocino', 'Napa', 'Sonoma'\n",
    "]\n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "santa_rosa_in_california_county_list, state = 'California', \n",
    "diocese = 'Santa Rosa in California')\n",
    "\n",
    "\n",
    "# Diocese of Stockton:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Stockton\n",
    "\n",
    "stockton_county_list = [\n",
    "'Alpine', 'Calaveras', 'Mono', 'San Joaquin County', 'Stanislaus', 'Tuolumne'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "stockton_county_list, state = 'California', \n",
    "diocese = 'Stockton')\n",
    "\n",
    "# Diocese of Las Vegas:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Las_Vegas\n",
    "\n",
    "las_vegas_county_list = [\n",
    "'Clark', 'Esmeralda', 'Lincoln', 'Nye', 'White Pine County'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "las_vegas_county_list, state = 'Nevada', \n",
    "diocese = 'Las Vegas')\n",
    "\n",
    "\n",
    "\n",
    "# Diocese of Gallup:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Gallup\n",
    "gallup_county_list_nm = [\n",
    "'San Juan County', 'McKinley', 'Cibola', 'Catron']\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "gallup_county_list_nm, state = 'New Mexico', \n",
    "diocese = 'Gallup')\n",
    "\n",
    "# Gallup also has territory in Arizona:\n",
    "gallup_county_list_az = [\n",
    "'Navajo', 'Apache'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "gallup_county_list_az, state = 'Arizona', \n",
    "diocese = 'Gallup')\n",
    "\n",
    "# Diocese of Las Cruces:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Las_Cruces\n",
    "\n",
    "las_cruces_county_list = [\n",
    "'Hidalgo', 'Grant', 'Luna', 'Sierra', 'Doña Ana County', 'Otero', 'Lincoln', \n",
    "'Chaves', 'Eddy', 'Lea'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "las_cruces_county_list, state = 'New Mexico', \n",
    "diocese = 'Las Cruces')\n",
    "\n",
    "# Diocese of Tucson:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Tucson\n",
    "tucson_county_list = [\n",
    "'Gila', 'Graham', 'Greenlee', 'Pinal', 'Cochise', 'Santa Cruz County', \n",
    "'Pima', 'Yuma', 'La Paz County' \n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "tucson_county_list, state = 'Arizona', \n",
    "diocese = 'Tucson')\n",
    "\n",
    "# Archdiocese of Louisville:\n",
    "# Source: https://www.archlou.org/about-the-archdiocese/history/statistics/\n",
    "# Note: This list misspells Metcalfe County as 'Metcalf' and Russell County\n",
    "# as 'Russel.'  \n",
    "louisville_county_list = [\n",
    "'Meade', 'Trimble', 'Oldham', 'Henry', 'Jefferson',\n",
    "'Shelby', 'Spencer', 'Bullitt', 'Hardin', 'Nelson', 'Washington', 'Larue', \n",
    "'Marion', 'Hart', 'Green', 'Taylor', 'Casey', 'Barren', 'Metcalfe', 'Adair', \n",
    "'Russell', 'Monroe', 'Cumberland', 'Clinton'\n",
    "] \n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "louisville_county_list, state = 'Kentucky', \n",
    "diocese = 'Louisville')\n",
    "\n",
    "# The Diocese of Lexington counties were added in earlier, but here's a map of\n",
    "# the diocese's counties anyway: \n",
    "# https://www.stjoanhershey.org/blog/special-collection-mission-cooperative-aug-5-6-2017\n",
    "\n",
    "# Diocese of Covington:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Covington\n",
    "covington_county_list = [\n",
    "'Boone', 'Kenton', 'Campbell', 'Gallatin', 'Carroll', 'Grant', 'Owen', \n",
    "'Pendleton', 'Harrison', 'Bracken', 'Robertson', 'Mason', 'Fleming', 'Lewis']\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "covington_county_list, state = 'Kentucky', \n",
    "diocese = 'Covington')\n",
    "\n",
    "# Diocese of Owensboro:\n",
    "# Source: https://owensborodiocese.org/our-diocese/\n",
    "owensboro_county_list = [\n",
    "'Fulton', 'Hickman', 'Carlisle', 'Ballard', 'McCracken', 'Graves', 'Marshall',\n",
    "'Calloway', 'Trigg', 'Lyon', 'Caldwell', 'Crittenden', 'Union', 'Henderson', \n",
    "'Webster', 'Hancock', 'Ohio', 'Butler', 'Warren', 'Edmonson', 'Grayson',\n",
    "'Breckinridge', 'Todd', 'Logan', 'Simpson', 'Allen', 'Livingston', \n",
    "'Christian', 'Hopkins', 'Muhlenberg', 'McLean', 'Daviess' \n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "owensboro_county_list, state = 'Kentucky', \n",
    "diocese = 'Owensboro')\n",
    "\n",
    "# Diocese of Knoxville:\n",
    "# Source: https://dioknox.org/about (Doesn't include county names, so I looked\n",
    "# on the map of US County shapefiles to identify them)\n",
    "\n",
    "knoxville_county_list = [\n",
    "'Polk', 'Bradley', 'Hamilton', 'Marion', 'Sequatchie',\n",
    "'Bledsoe', 'Cumberland', 'Fentress', 'Pickett', 'Scott', 'Morgan', 'Roane',\n",
    "'Rhea', 'Meigs', 'McMinn', 'Monroe', 'Blount', 'Knox', 'Anderson', 'Campbell',\n",
    "'Claiborne', 'Union', 'Grainger', 'Jefferson', 'Sevier', 'Cocke', 'Greene',\n",
    "'Hawkins', 'Sullivan', 'Washington', 'Unicoi', 'Carter', 'Johnson', 'Loudon',\n",
    "'Hamblen', 'Hancock'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "knoxville_county_list, state = 'Tennessee', \n",
    "diocese = 'Knoxville')\n",
    "\n",
    "# Diocese of Memphis:\n",
    "# Source: https://cdom.org/wp-content/uploads/2021/08/Aug-2021-Parish-Directory-Public.pdf\n",
    "# (Doesn't include county names, so I used the map of US County shapefiles\n",
    "# to identify them)\n",
    "memphis_county_list = [\n",
    "'Shelby', 'Tipton', 'Lauderdale', 'Dyer', 'Lake', 'Obion', 'Crockett', \n",
    "'Hardin', 'Decatur', 'Benton', 'Henry', 'Carroll', 'Henderson', 'Chester', \n",
    "'McNairy', 'Hardeman', 'Fayette', 'Haywood', 'Gibson', 'Weakley', 'Madison'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "memphis_county_list, state = 'Tennessee', \n",
    "diocese = 'Memphis')\n",
    "\n",
    "# Archdiocese of Miami:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_Miami\n",
    "miami_county_list = [\n",
    "'Broward', 'Miami-Dade', 'Monroe'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "miami_county_list, state = 'Florida', \n",
    "diocese = 'Miami')\n",
    "\n",
    "# Diocese of Palm Beach:\n",
    "# Source: https://www.diocesepb.org/about-us/\n",
    "\n",
    "palm_beach_county_list = [\n",
    "'Palm Beach County', 'Martin', 'St. Lucie County', 'Indian River County', \n",
    "'Okeechobee'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "palm_beach_county_list, state = 'Florida', \n",
    "diocese = 'Palm Beach')\n",
    "\n",
    "# Diocese of Venice in Florida:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Venice_in_Florida\n",
    "\n",
    "venice_in_florida_county_list = [\n",
    "'Charlotte', 'Collier', 'DeSoto', 'Glades', 'Hardee', 'Hendry', 'Highlands', \n",
    "'Lee', 'Manatee', 'Sarasota'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "venice_in_florida_county_list, state = 'Florida', \n",
    "diocese = 'Venice')\n",
    "\n",
    "# Diocese of Pensacola-Tallahassee:\n",
    "# Source: https://ptdiocese.org/pictures/2017/1/Standard%20Deanery%20Map%202017.png\n",
    "\n",
    "pensacola_tallahassee_county_list = [\n",
    "'Escambia', 'Santa Rosa County', 'Okaloosa', 'Walton', 'Holmes', 'Washington',\n",
    "'Jackson', 'Bay', 'Calhoun', 'Gulf', 'Liberty', 'Franklin', 'Wakulla', \n",
    "'Leon', 'Gadsden', 'Jefferson', 'Madison', 'Taylor'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "pensacola_tallahassee_county_list, state = 'Florida', \n",
    "diocese = 'Pensacola-Tallahassee')\n",
    "\n",
    "\n",
    "# Diocese of Saint Petersburg:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Saint_Petersburg\n",
    "saint_petersburg_county_list = [\n",
    "'Pinellas', 'Hillsborough', 'Pasco', 'Hernando', 'Citrus'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "saint_petersburg_county_list, state = 'Florida', \n",
    "diocese = 'Saint Petersburg')\n",
    "\n",
    "# Diocese of St. Augustine:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_St._Augustine\n",
    "st_augustine_county_list = [\n",
    "'Alachua', 'Baker', 'Bradford', 'Clay', 'Columbia', 'Dixie', 'Duval', 'Flagler',\n",
    "'Gilchrist', 'Hamilton', 'Lafayette', 'Levy', 'Nassau', 'Putnam', \n",
    "'St. Johns County', 'Suwannee', 'Union'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "st_augustine_county_list, state = 'Florida', \n",
    "diocese = 'St. Augustine')\n",
    "\n",
    "# Archdiocese of Milwaukee:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_Milwaukee\n",
    "\n",
    "milwaukee_county_list = [\n",
    "'Dodge', 'Fond du Lac County', 'Kenosha', 'Milwaukee', 'Ozaukee', 'Racine',\n",
    "'Sheboygan', 'Walworth', 'Washington', 'Waukesha'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "milwaukee_county_list, state = 'Wisconsin', \n",
    "diocese = 'Milwaukee')\n",
    "\n",
    "# Diocese of Green Bay:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Green_Bay\n",
    "\n",
    "green_bay_county_list = [\n",
    "'Brown', 'Calumet', 'Door', 'Florence', 'Forest', 'Kewaunee', 'Langlade', \n",
    "'Manitowoc', 'Marinette', 'Menominee', 'Oconto', 'Outagamie', 'Shawano',\n",
    "'Waupaca', 'Waushara', 'Winnebago'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "green_bay_county_list, state = 'Wisconsin', \n",
    "diocese = 'Green Bay')\n",
    "\n",
    "# Diocese of Madison:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Madison\n",
    "\n",
    "madison_county_list = ['Columbia', 'Dane', 'Grant', 'Green', \n",
    "'Green Lake County', 'Iowa', 'Jefferson', 'Lafayette', 'Marquette', \n",
    "'Rock', 'Sauk'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "madison_county_list, state = 'Wisconsin', \n",
    "diocese = 'Madison')\n",
    "\n",
    "# Diocese of Superior:\n",
    "# https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Superior\n",
    "\n",
    "superior_county_list = [\n",
    "'Ashland', 'Barron', 'Bayfield', 'Burnett', 'Douglas', 'Iron', 'Lincoln', \n",
    "'Oneida', 'Price', 'Polk', 'Rusk', 'Sawyer', 'St. Croix County', 'Taylor',\n",
    "'Vilas', 'Washburn'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "superior_county_list, state = 'Wisconsin', \n",
    "diocese = 'Superior')\n",
    "\n",
    "# Diocese of Biloxi:\n",
    "# Source: https://biloxidiocese.org/geographics\n",
    "\n",
    "biloxi_county_list = [\n",
    "'Covington', 'Forrest', 'George', 'Greene', 'Hancock', 'Harrison', 'Jackson', \n",
    "'Jefferson Davis County', 'Jones', 'Lamar', 'Lawrence', 'Marion', \n",
    "'Pearl River County', 'Perry', 'Stone', 'Walthall', 'Wayne'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "biloxi_county_list, state = 'Mississippi', \n",
    "diocese = 'Biloxi')\n",
    "\n",
    "\n",
    "# Diocese of Birmingham:\n",
    "# Source: https://bhmdiocese.org/documents/2021/12/Dio%20info%2012-13-21.pdf\n",
    "\n",
    "birmingham_county_list = [\n",
    "'Lauderdale', 'Limestone', 'Madison', 'Jackson', 'Colbert', 'Franklin', \n",
    "'Lawrence', 'Morgan', 'Marshall', 'DeKalb', 'Marion', 'Winston', 'Cullman', \n",
    "'Blount', 'Etowah', 'Cherokee', 'Lamar', 'Fayette', 'Walker', 'Jefferson',\n",
    "'St. Clair County', 'Calhoun', 'Cleburne', 'Pickens', 'Tuscaloosa', 'Bibb', \n",
    "'Shelby', 'Talladega', 'Clay', 'Randolph', 'Sumter', 'Greene', 'Hale', 'Perry', \n",
    "'Chilton', 'Coosa', 'Tallapoosa', 'Chambers', 'Marengo'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "birmingham_county_list, state = 'Alabama', \n",
    "diocese = 'Birmingham')\n",
    "\n",
    "# Archdiocese of New Orleans:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_New_Orleans\n",
    "new_orleans_county_list = [\n",
    "'Jefferson Parish', 'Orleans Parish', 'Plaquemines Parish', \n",
    "'St. Bernard Parish', 'St. Charles Parish', 'St. John the Baptist Parish',\n",
    "'St. Tammany Parish', 'Washington Parish'\n",
    "] \n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "new_orleans_county_list, state = 'Louisiana', \n",
    "diocese = 'New Orleans')\n",
    "\n",
    "# Diocese of Houma-Thibodaux:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Houma%E2%80%93Thibodaux\n",
    "\n",
    "houma_thibodaux_county_list = ['Terrebonne Parish', 'Lafourche Parish']\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "houma_thibodaux_county_list, state = 'Louisiana', \n",
    "diocese = 'Houma-Thibodaux')\n",
    "\n",
    "\n",
    "# Diocese of Lafayette in Louisiana:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Lafayette_in_Louisiana\n",
    " \n",
    "lafayette_in_louisiana_county_list = [\n",
    "'St. Landry Parish', 'Evangeline Parish', 'Lafayette Parish', \n",
    "'St. Martin Parish', 'Iberia Parish', 'St. Mary Parish', 'Acadia Parish',\n",
    "'Vermilion Parish'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "lafayette_in_louisiana_county_list, state = 'Louisiana', \n",
    "diocese = 'Lafayette in Louisiana')\n",
    "\n",
    "# Diocese of Baton Rouge:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Baton_Rouge\n",
    "\n",
    "baton_rouge_county_list = [\n",
    "'Ascension Parish', 'Assumption Parish', 'East Baton Rouge Parish', \n",
    "'East Feliciana Parish', 'Iberville Parish', 'Livingston Parish', \n",
    "'Pointe Coupee Parish', 'Tangipahoa Parish', 'St. Helena Parish', \n",
    "'St. James Parish', 'West Baton Rouge Parish', 'West Feliciana Parish'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "baton_rouge_county_list, state = 'Louisiana', \n",
    "diocese = 'Baton Rouge')\n",
    "\n",
    "# Diocese of Lake Charles:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Lake_Charles\n",
    "\n",
    "lake_charles_county_list = [\n",
    "'Allen Parish', 'Beauregard Parish', 'Calcasieu Parish', 'Cameron Parish',\n",
    "'Jefferson Davis Parish'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "lake_charles_county_list, state = 'Louisiana', \n",
    "diocese = 'Lake Charles')\n",
    "\n",
    "# Diocese of Alexandria in Louisiana:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Alexandria_in_Louisiana\n",
    "\n",
    "alexandria_in_louisiana_county_list = [\n",
    "'Avoyelles Parish', 'Rapides Parish', 'Vernon Parish', 'Natchitoches Parish',\n",
    "'Winn Parish', 'Caldwell Parish', 'Madison Parish', 'Franklin Parish', \n",
    "'Tensas Parish', 'Concordia Parish', 'Catahoula Parish', 'LaSalle Parish',\n",
    "'Grant Parish'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "alexandria_in_louisiana_county_list, state = 'Louisiana', \n",
    "diocese = 'Alexandria in Louisiana')\n",
    "\n",
    "# Archdiocese of New York:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_New_York\n",
    "new_york_county_list = ['New York County', 'Richmond', 'Bronx', 'Dutchess',\n",
    "'Orange', 'Putnam', 'Rockland', 'Sullivan', 'Ulster', 'Westchester'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "new_york_county_list, state = 'New York', \n",
    "diocese = 'New York')\n",
    "\n",
    "# Diocese of Brooklyn:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Brooklyn\n",
    "brooklyn_county_list = [\n",
    "'Kings', 'Queens'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "brooklyn_county_list, state = 'New York', \n",
    "diocese = 'Brooklyn')\n",
    "\n",
    "# Diocese of Buffalo:\n",
    "# https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Buffalo\n",
    "\n",
    "buffalo_county_list = [\n",
    "'Erie', 'Niagara', 'Genesee', 'Orleans', 'Chautauqua', 'Wyoming', \n",
    "'Cattaraugus', 'Allegany'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "buffalo_county_list, state = 'New York', \n",
    "diocese = 'Buffalo')\n",
    "\n",
    "# Diocese of Rockville Centre:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Rockville_Centre\n",
    "\n",
    "rockville_centre_county_list = [\n",
    "'Nassau', 'Suffolk'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "rockville_centre_county_list, state = 'New York', \n",
    "diocese = 'Rockville Centre')\n",
    "\n",
    "# Diocese of Ogdensburg:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Ogdensburg\n",
    "\n",
    "ogdensburg_county_list = [\n",
    "'Clinton', 'Essex', 'Franklin', 'Jefferson', 'Lewis', 'St. Lawrence County', \n",
    "'Hamilton'\n",
    "] \n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "ogdensburg_county_list, state = 'New York', \n",
    "diocese = 'Ogdensburg')\n",
    "\n",
    "# Diocese of Rochester:\n",
    "# Source: \n",
    "\n",
    "rochester_county_list = [\n",
    "'Monroe', 'Cayuga', 'Livingston', 'Wayne', 'Tioga', 'Tompkins', 'Ontario', \n",
    "'Seneca', 'Schuyler', 'Yates', 'Steuben', 'Chemung'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "rochester_county_list, state = 'New York', \n",
    "diocese = 'Rochester')\n",
    "\n",
    "# Diocese of Syracuse: \n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Syracuse\n",
    "\n",
    "syracuse_county_list = [\n",
    "'Broome', 'Chenango', 'Cortland', 'Madison', 'Oneida', 'Onondaga', 'Oswego']\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "syracuse_county_list, state = 'New York', \n",
    "diocese = 'Syracuse')\n",
    "\n",
    "# Archdiocese of Neawrk:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_Newark\n",
    "\n",
    "newark_county_list = ['Bergen', 'Union', 'Hudson', 'Essex']\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "newark_county_list, state = 'New Jersey', \n",
    "diocese = 'Newark')\n",
    "\n",
    "# Diocese of Metuchen:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Metuchen\n",
    "\n",
    "metuchen_county_list = [\n",
    "'Middlesex', 'Somerset', 'Hunterdon', 'Warren'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "metuchen_county_list, state = 'New Jersey', \n",
    "diocese = 'Metuchen')\n",
    "\n",
    "# Diocese of Paterson:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Paterson\n",
    "\n",
    "paterson_county_list = [\n",
    "'Passaic', 'Morris', 'Sussex'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "paterson_county_list, state = 'New Jersey', \n",
    "diocese = 'Paterson')\n",
    "\n",
    "# Diocese of Trenton:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Trenton\n",
    "\n",
    "trenton_county_list = [\n",
    "'Burlington', 'Mercer', 'Monmouth', 'Ocean'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "trenton_county_list, state = 'New Jersey', \n",
    "diocese = 'Trenton')\n",
    "\n",
    "\n",
    "# Diocese of Tulsa:\n",
    "# Source: https://d2y1pz2y630308.cloudfront.net/2682/documents/2022/3/revised%20vicariate%20full%20map.pdf\n",
    "# (This map does not include county names, so I used the county shapefile \n",
    "# map to identify counties within the diocese.)\n",
    "\n",
    "tulsa_county_list = [\n",
    "'McCurtain', 'Choctaw', 'Bryan', 'Atoka', 'Pushmataha', 'Latimer', 'Coal', \n",
    "'Hughes', 'Okfuskee', 'Creek', 'Payne', 'Pawnee', 'Osage', \n",
    "'Washington', 'Nowata', 'Craig', 'Ottawa', 'Delaware', 'Adair', 'Sequoyah',\n",
    "'Le Flore County', 'Pittsburg', 'McIntosh', 'Okmulgee',\n",
    "'Tulsa', 'Rogers', 'Wagoner', 'Cherokee', 'Muskogee', 'Haskell', 'Mayes']\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "tulsa_county_list, state = 'Oklahoma', \n",
    "diocese = 'Tulsa')\n",
    "\n",
    "# Archdiocese of Ohama:\n",
    "# Source: https://archomaha.org/about-us/\n",
    "# Note: Parish data on the Grand Island and Lincoln diocesan webpages \n",
    "# indicate that Wheeler, Greeley, Howard, and Hall Counties do not belong\n",
    "# to the Archdiocese of Omaha, so I've excluded them from the list below.\n",
    "# See https://gidiocese.org/parishfinder (note that one of the parishes is\n",
    "# in Wheeler County) and https://www.lincolndiocese.org/directory/parish-map .\n",
    "\n",
    "omaha_county_list = [\n",
    "'Boyd', 'Holt', 'Merrick', 'Nance', 'Boone', 'Antelope', 'Knox', 'Pierce', \n",
    "'Madison', 'Platte', 'Colfax', 'Stanton', 'Wayne', 'Cedar', 'Dixon', 'Dakota',\n",
    "'Thurston', 'Cuming', 'Dodge', 'Burt', 'Washington',\n",
    "'Douglas', 'Sarpy']\n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "omaha_county_list, state = 'Nebraska', \n",
    "diocese = 'Omaha')\n",
    "\n",
    "# Archdiocese of Lincoln:\n",
    "# Source: https://commons.wikimedia.org/wiki/File:Diocese_of_Lincoln_map_1.png\n",
    "# See also: https://www.lincolndiocese.org/diocese/about-the-diocese/history-of-the-diocese\n",
    "\n",
    "lincoln_county_list = [\n",
    "'Richardson', 'Nemaha', 'Otoe', 'Cass', 'Pawnee', 'Johnson', 'Lancaster', \n",
    "'Saunders', 'Butler', 'Polk', 'Hamilton', 'Clay', 'Nuckolls', 'Thayer', \n",
    "'Jefferson', 'Gage', 'Webster', 'Adams', 'Kearney', 'Franklin', 'Harlan', \n",
    "'Furnas', 'Gosper', 'Lincoln', 'Frontier', 'Red Willow County', 'Hitchcock',\n",
    "'Dundy', 'Hayes', 'Chase', 'Perkins', 'Phelps', 'York', 'Seward', 'Fillmore', \n",
    "'Saline'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "lincoln_county_list, state = 'Nebraska', \n",
    "diocese = 'Lincoln')\n",
    "\n",
    "# Archdiocese of Philadelphia:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_Philadelphia\n",
    "\n",
    "philadelphia_county_list = [\n",
    "'Bucks', 'Chester', 'Delaware', 'Montgomery', 'Philadelphia'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "philadelphia_county_list, state = 'Pennsylvania', \n",
    "diocese = 'Philadelphia')\n",
    "\n",
    "# Diocese of Allentown:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Allentown\n",
    "\n",
    "allentown_county_list = [\n",
    "'Berks', 'Carbon', 'Lehigh', 'Northampton', 'Schuylkill'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "allentown_county_list, state = 'Pennsylvania', \n",
    "diocese = 'Allentown')\n",
    "\n",
    "# Diocese of Altoona-Johnstown:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Altoona%E2%80%93Johnstown\n",
    "\n",
    "\n",
    "altoona_johnstown_county_list = [\n",
    "'Bedford', 'Blair', 'Cambria', 'Centre', 'Clinton', 'Fulton', 'Huntingdon',\n",
    "'Somerset'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "altoona_johnstown_county_list, state = 'Pennsylvania', \n",
    "diocese = 'Altoona-Johnstown')\n",
    "\n",
    "# Diocese of Harrisburg:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Harrisburg\n",
    "\n",
    "harrisburg_county_list = [\n",
    "'Adams', 'Columbia', 'Cumberland', 'Dauphin', 'Franklin', 'Juniata', \n",
    "'Lancaster', 'Lebanon', 'Mifflin', 'Montour', 'Northumberland', 'Perry', \n",
    "'Snyder', 'Union', 'York'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "harrisburg_county_list, state = 'Pennsylvania', \n",
    "diocese = 'Harrisburg')\n",
    "\n",
    "# Diocese of Greensburg:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Greensburg\n",
    "\n",
    "greensburg_county_list = [\n",
    "'Armstrong', 'Fayette', 'Indiana', 'Westmoreland'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "greensburg_county_list, state = 'Pennsylvania', \n",
    "diocese = 'Greensburg')\n",
    "\n",
    "# Diocese of Pittsburgh:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Pittsburgh\n",
    "\n",
    "pittsburgh_county_list = [\n",
    "'Allegheny', 'Beaver', 'Butler', 'Greene', 'Lawrence', 'Washington'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "pittsburgh_county_list, state = 'Pennsylvania', \n",
    "diocese = 'Pittsburgh')\n",
    "\n",
    "# Diocese of Scranton:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Scranton\n",
    "scranton_county_list = [\n",
    "'Lackawanna', 'Luzerne', 'Bradford', 'Susquehanna', 'Wayne', 'Tioga', \n",
    "'Sullivan', 'Wyoming', 'Lycoming', 'Pike', 'Monroe'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "scranton_county_list, state = 'Pennsylvania', \n",
    "diocese = 'Scranton')\n",
    "\n",
    "# Diocese of Baker:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Baker\n",
    "\n",
    "baker_county_list = [\n",
    "'Baker', 'Crook', 'Deschutes', 'Gilliam', 'Grant', 'Harney', \n",
    "'Hood River County', 'Jefferson', 'Klamath', 'Lake', 'Malheur', 'Morrow', \n",
    "'Sherman', 'Umatilla', 'Union', 'Wallowa', 'Wasco', 'Wheeler'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "baker_county_list, state = 'Oregon', \n",
    "diocese = 'Oregon')\n",
    "\n",
    "# Diocese of Great Falls-Billings:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Great_Falls%E2%80%93Billings\n",
    "# Note: I updated the source to include McCone county, which is encompassed\n",
    "# by other counties in the diocese.\n",
    "\n",
    "great_falls_billings_county_list = [\n",
    "'Big Horn County', 'Blaine', 'Carbon', 'Carter', 'Cascade', 'Chouteau', \n",
    "'Custer', 'Daniels', 'Dawson', 'Fallon', 'Fergus', 'Garfield', \n",
    "'Golden Valley County', 'Hill', 'Judith Basin County', 'Liberty', 'McCone',\n",
    "'Musselshell', 'Park', 'Petroleum', 'Phillips', 'Powder River County', \n",
    "'Prairie', 'Richland', 'Roosevelt', 'Rosebud', 'Sheridan', 'Stillwater', \n",
    "'Sweet Grass County', 'Treasure', 'Valley', 'Wibaux', 'Yellowstone'\n",
    "]\n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "great_falls_billings_county_list, state = 'Montana', \n",
    "diocese = 'Great Falls-Billings')\n",
    "\n",
    "# Archdiocese of St. Louis:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_St._Louis\n",
    "\n",
    "st_louis_county_list = [\n",
    "'Franklin', 'Jefferson', 'Lincoln', 'Perry', 'St. Charles County', \n",
    "'St. Francois County', 'Ste. Genevieve County', 'St. Louis County', 'Warren', \n",
    "'Washington', 'St. Louis city']\n",
    "\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "st_louis_county_list, state = 'Missouri', \n",
    "diocese = 'St. Louis')\n",
    "\n",
    "# Diocese of Kansas City-St. Joseph:\n",
    "# Source: https://www.lifeandjusticekcsj.org/diocese-map.html\n",
    "# (Website is affiliated with two diocesan offices)\n",
    "\n",
    "kansas_city_st_joseph_county_list = [\n",
    "'Andrew', 'Atchison', 'Bates', 'Buchanan', 'Caldwell', 'Carroll', 'Cass', \n",
    "'Clay', 'Clinton', 'Daviess', 'DeKalb', 'Gentry', 'Grundy', 'Harrison', 'Henry',\n",
    "'Holt', 'Jackson', 'Johnson', 'Lafayette', 'Livingston', 'Mercer', 'Nodaway',\n",
    "'Platte', 'Ray', 'St. Clair County', 'Vernon', 'Worth']\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "kansas_city_st_joseph_county_list, state = 'Missouri', \n",
    "diocese = 'Kansas City-St. Joseph')\n",
    "\n",
    "# Diocese of Jefferson City:\n",
    "# Source: https://diojeffcity.org/about-us/\n",
    "\n",
    "jefferson_city_county_list = [\n",
    "'Putnam', 'Schuyler', 'Scotland', 'Clark', 'Sullivan', 'Adair', 'Knox', 'Lewis',\n",
    "'Linn', 'Macon', 'Shelby', 'Marion', 'Chariton', 'Randolph', 'Monroe', 'Ralls',\n",
    "'Pike', 'Saline', 'Howard', 'Boone', 'Audrain', 'Pettis', 'Cooper', 'Moniteau',\n",
    "'Cole', 'Callaway', 'Osage', 'Montgomery', 'Gasconade', 'Benton', 'Morgan', \n",
    "'Hickory', 'Camden', 'Miller', 'Maries', 'Pulaski', 'Phelps', 'Crawford'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "jefferson_city_county_list, state = 'Missouri', \n",
    "diocese = 'Jefferson City')\n",
    "\n",
    "# Archdiocese of Saint Paul and Minneapolis:\n",
    "# Source: https://commons.wikimedia.org/wiki/File:Map_of_the_Catholic_archdiocese_of_Saint_Paul_%26_Minneapolis.svg\n",
    "# See also: https://s3.amazonaws.com/archspmmainsite/Resources/ArchMap--Parishes-Deaneries-COLORRamped_061516.pdf\n",
    "\n",
    "saint_paul_and_minneapolis_county_list = [\n",
    "'Goodhue', 'Rice', 'Le Sueur County', 'Scott', 'Dakota', 'Washington', 'Ramsey',\n",
    "'Hennepin', 'Wright', 'Anoka', 'Chisago', 'Carver'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "saint_paul_and_minneapolis_county_list, state = 'Minnesota', \n",
    "diocese = 'Saint Paul and Minneapolis')\n",
    "\n",
    "# Diocese of Crookston:\n",
    "# https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Crookston\n",
    "\n",
    "crookston_county_list = [\n",
    "'Kittson', 'Roseau', 'Lake of the Woods County', 'Marshall', 'Polk', \n",
    "'Red Lake County', 'Pennington', 'Clearwater', 'Beltrami', 'Norman', 'Mahnomen', \n",
    "'Hubbard', 'Clay', 'Becker'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "crookston_county_list, state = 'Minnesota', \n",
    "diocese = 'Crookston')\n",
    "\n",
    "# Diocese of New Ulm:\n",
    "# Source: https://www.dnu.org/about\n",
    "\n",
    "new_ulm_county_list = [\n",
    "'Big Stone County', 'Brown', 'Chippewa', 'Kandiyohi', 'Lac qui Parle County',\n",
    "'Lincoln', 'Lyon', 'McLeod', 'Meeker', 'Nicollet', 'Redwood', 'Renville',\n",
    "'Sibley', 'Swift', 'Yellow Medicine County'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "new_ulm_county_list, state = 'Minnesota', \n",
    "diocese = 'New Ulm')\n",
    "\n",
    "# Diocese of Saint Cloud:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Saint_Cloud\n",
    "\n",
    "saint_cloud_county_list = [\n",
    "'Benton', 'Douglas', 'Grant', 'Isanti', 'Kanabec', 'Mille Lacs County', \n",
    "'Morrison', 'Otter Tail County', 'Pope', 'Sherburne', 'Stearns', 'Stevens',\n",
    "'Todd', 'Traverse', 'Wadena', 'Wilkin'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "saint_cloud_county_list, state = 'Minnesota', \n",
    "diocese = 'Saint Cloud')\n",
    "\n",
    "\n",
    "# Diocese of Winona-Rochester:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Winona%E2%80%93Rochester\n",
    "\n",
    "winona_rochester_county_list = [\n",
    "'Blue Earth County', 'Cottonwood', 'Dodge', 'Faribault', 'Fillmore', 'Freeborn',\n",
    "'Houston', 'Jackson', 'Martin', 'Mower', 'Murray', 'Nobles', 'Olmsted',\n",
    "'Pipestone', 'Rock', 'Steele', 'Wabasha', 'Waseca', 'Watonwan', 'Winona'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "winona_rochester_county_list, state = 'Minnesota', \n",
    "diocese = 'Winona-Rochester')\n",
    "\n",
    "# Diocese of Bismark:\n",
    "# Source https://bismarckdiocese.com/about\n",
    "\n",
    "bismark_county_list = [\n",
    "'Adams', 'Billings', 'Bowman', 'Burke', 'Burleigh', 'Divide', 'Dunn', 'Emmons', \n",
    "'Golden Valley County', 'Grant', 'Hettinger', 'McKenzie', 'McLean', 'Mercer', \n",
    "'Morton', 'Mountrail', 'Oliver', 'Renville', 'Sioux', 'Slope', 'Stark', 'Ward', \n",
    "'Williams'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "bismark_county_list, state = 'North Dakota', \n",
    "diocese = 'Bismark')\n",
    "\n",
    "# Diocese of Rapid City:\n",
    "# Source: https://www.rapidcitydiocese.org/wp-content/uploads/2022/01/MapDeaneryDec21-1030x752-1.jpeg\n",
    "\n",
    "rapid_city_county_list = [\n",
    "'Corson', 'Dewey', 'Stanley', 'Lyman', 'Gregory', 'Tripp', 'Todd', 'Bennett',\n",
    "'Oglala Lakota County', 'Fall River County', 'Custer', 'Pennington', 'Lawrence',\n",
    "'Meade', 'Butte', 'Harding', 'Perkins', 'Ziebach', 'Haakon', 'Jackson',\n",
    "'Jones', 'Mellette']\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "rapid_city_county_list, state = 'South Dakota', \n",
    "diocese = 'Rapid City')\n",
    "\n",
    "# Diocese of Spokane:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Spokane\n",
    "\n",
    "spokane_county_list = [\n",
    "'Okanogan', 'Ferry', 'Stevens', 'Pend Oreille County', 'Lincoln', 'Spokane',\n",
    "'Adams', 'Whitman', 'Franklin', 'Walla Walla County', 'Columbia', 'Garfield',\n",
    "'Asotin'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "spokane_county_list, state = 'Washington', \n",
    "diocese = 'Spokane')\n",
    "\n",
    "# Diocese of Yakkima:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Yakima\n",
    "\n",
    "yakima_county_list = [\n",
    "'Benton', 'Chelan', 'Douglas', 'Grant', 'Kittitas', 'Klickitat', 'Yakima'\n",
    "]\n",
    "merged_county_shapes_test = add_dioceses(merged_county_shapes_test, \n",
    "yakima_county_list, state = 'Washington', \n",
    "diocese = 'Yakima')\n",
    "\n",
    "\n",
    "# county_list = [\n",
    "# '', '', '', '', '', '', '', '', '', '',\n",
    "# '', '', '', '', '', '', '', '', '', '',\n",
    "# '', '', '', '', '', '', '', '', '', '',\n",
    "# '', '', '', '', '', '', '', '', '', '',\n",
    "# '', '', '', '', '', '', '', '', '', '',\n",
    "# ]\n",
    "\n",
    "# List of Catholic Dioceses in the US:\n",
    "# https://en.wikipedia.org/wiki/List_of_Catholic_dioceses_in_the_United_States"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here's a look at how many counties' dioceses were added in via the above block:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "counties_changed = [1 for i in range(len(merged_county_shapes_test)) \n",
    "if merged_county_shapes_test.iloc[i, merged_county_shapes_test.columns.get_loc(\n",
    "    'Diocese')] != initial_version_of_merged_county_shapes_test.iloc[\n",
    "        i, initial_version_of_merged_county_shapes_test.columns.get_loc(\n",
    "            'Diocese')]]\n",
    "len(counties_changed) # The number of counties that were modified above"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Inaccurate parts of the map at this point:\n",
    "\n",
    "1. Fall River/Boston: The Diocese of Fall River includes \"the towns of Marion, Mattapoisett, and Wareham along the south coast of Plymouth County in Massachusetts, in the New England region of the United States.\" (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Fall_River) The rest of Plymouth County belongs to the Archdiocese of Boston. \n",
    "\n",
    "2. Austin/Victoria in Texas: Only part of Fayette County belongs to the Diocese of Victoria in Texas. \"The part of Fayette County north of the Colorado River\" is part of the Diocese of Austin. (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Austin)\n",
    "\n",
    "3. San Antonio/Corpus Christi: The Archdiocese of San Antonio includes only part of McMullen County. The rest belongs to the Diocese of Corpus Christi. See: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Corpus_Christi \n",
    "This also means that McMullen County is part of two different provinces (Galveston-Houston and San Antonio).\n",
    "\n",
    "4. Evansville/Indianapolis: Wikipedia notes that \"Harrison Township in Spencer County, the location of St. Meinrad Archabbey, is part of the Archdiocese of Indianapolis.\" (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Evansville)\n",
    "\n",
    "5. Gallup/Santa Fe: Wikipedia notes that the Diocsee of Gallup also includes \"parts of Rio Arriba, Sandoval, Bernalillo, and Valencia Counties west of 106,52',41\" meridian in New Mexico.\" The rest of those counties belong to Santa Fe. (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Gallup)\n",
    "\n",
    "6. Tucson/Phoenix: The parts of Pinal County that lie within \"the territorial boundaries of the Gila River Indian Community\" are part of the Diocese of Phoenix, not the Diocese of Tucson. (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Tucson)\n",
    "\n",
    "7. New Orleans/Houma-Thibodaux: The Archdiocese of New Orleans does not include Grand Isle (within Jefferson Parish); that instead belongs to the Diocese of Houma-Thibodaux. (https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_New_Orleans)\n",
    "\n",
    "8. Lafayette in Louisiana/Houma-Thibodaux: The eastern portion of St. Mary Parish (including Morgan City) belongs to the Diocese of Houma-Thibodaux; the rest belongs to the Diocese of Lafayette in Louisiana.\n",
    "\n",
    "9. Rockville Centre/Norwich: Fisher Island (within Suffolk County, NY) belongs to the Diocese of Norwich, not the Diocese of Rockville Centre.\n",
    "\n",
    "10. Albany/Ogdensburg: The northern part of Herkimer County belongs to the Diocese of Ogdensburg; the rest belongs to the Diocese of Albany. (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Ogdensburg\n",
    "\n",
    "11. Albany/Ogdensburg: The southern part of Hamilton County belongs to the Diocese of Albany; the rest belongs to the Diocese of Ogdensburg. (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Ogdensburg)\n",
    "\n",
    "13. Grand Island/Lincoln: the parts of Deuel, Keith, Lincoln, Dawson, and Hall counties that lie south of the Platte River belong to the Diocese of Lincoln. The northern parts of these counties belong to the Diocese of Grand Island.\n",
    "Although the Archdiocese of Omaha reports that \"parts of Hall County\" belong to its diocese (https://archomaha.org/about-us/), I believe this is incorrect, as the Diocese of Lincoln's parish map shows the southern portion of Hall County as belonging instead instead to the Diocsee of Lincoln: https://www.lincolndiocese.org/directory/parish-map\n",
    "For now, I have assigned Lincoln County to the Diocese of Lincoln and Deuel, Keith, Dawson, and Hall counties to the Diocese of Grand Island (based on which diocese has the majority of each's county's area). \n",
    "It appears that a small portion of Buffalo County may lie south of the Platte, but I'm keeping all of this county within the Diocese of Grand Island for now, since it's possible the Platte's course changed following the diocese's establishment (and the bottom edge of the county does appear to follow the course of a river).\n",
    "\n",
    "14. Helena/Great Falls-Billings: The Diocese of Great Falls-Billings includes parts of 2 counties, given that the Diocese of Helena does also: https://diocesehelena.org/about/ Based on the Diocese of Helena's Wikipedia map (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Helena#/media/File:Diocese_of_Helena.PNG), it looks like the Diocese of Great Falls-Billings includes portions of Toole and Teton counties, the majority of which belong to the Diocese of Helena.\n",
    "However, the Diocese of Great Falls-Billings map doesn't appear to reflect these modifications:\n",
    "https://diocesegfb.org/home-page/about/\n",
    "Meanwhile, the Diocese of Helena's wikipedia page reports that the diocese includes \"Parts of Meagher, Musselshell, and Toole counties in Montana.\" However, Musselshell county is landlocked within the Diocese of Great Falls-Billings, so this is likely a typo. I also don't see any evidence on the Wikipedia map (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Helena#/media/File:Diocese_of_Helena.PNG) of part of Meagher county belonging to the Diocsee of Great Falls-Billings.\n",
    "\n",
    "15. Bismark/Fargo: The western section of Bottineau County belongs to the Diocese of Bismark; the rest belongs to the Diocese of Fargo. See https://bismarckdiocese.com/about and page 3 of https://d2y1pz2y630308.cloudfront.net/2950/documents/2022/3/DiocesanDirectory.pdf."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "merged_county_shapes_test.sort_values('Diocese', inplace = True) \n",
    "# Sorts the dioceses alphabetically so that the colors are not determined by the order in which the dioceses were added\n",
    "df_diocese_numbers = pd.DataFrame(list(enumerate(pd.unique(\n",
    "    merged_county_shapes_test['Diocese']))))\n",
    "df_diocese_numbers.columns = ['Diocese_Code', 'Diocese']\n",
    "df_diocese_numbers"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "merged_county_shapes_test.head(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_cd = merged_county_shapes_test.merge(df_diocese_numbers, on = 'Diocese', \n",
    "how = 'left') # the 'cd' in df_cd stands for 'counties' and 'dioceses'\n",
    "df_cd.set_index('county_state_column', inplace = True, drop = False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_cd.index.rename('county_state', inplace = True)\n",
    "df_cd.head(2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Updating Color Codes:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# # This use of df.apply() with if/else statements is based on an example\n",
    "# by Professor Hardeep Johar.\n",
    "\n",
    "# Green:\n",
    "df_cd['Diocese_Code'] = df_cd.apply(lambda row: 48 if row['Diocese'] in \n",
    "['Seattle', 'Dubuque', 'Charleston'] else row['Diocese_Code'], axis = 1)\n",
    "\n",
    "# Light blue:\n",
    "df_cd['Diocese_Code'] = df_cd.apply(lambda row: 0 if row['Diocese'] in [\n",
    "    'Santa Rosa in California', 'Crookston', 'Kansas City-St. Joseph', \n",
    "    'Gaylord'] else row['Diocese_Code'], axis = 1)\n",
    "\n",
    "#Pink:\n",
    "df_cd['Diocese_Code'] = df_cd.apply(lambda row: 64 if row['Diocese'] in [\n",
    "    'Reno', 'San Antonio'] else row['Diocese_Code'], axis = 1)\n",
    "\n",
    "#Red:\n",
    "df_cd['Diocese_Code'] = df_cd.apply(lambda row: 80 if row['Diocese'] in [\n",
    "    'Cheyenne', 'Sioux Falls', 'Baltimore'] else row['Diocese_Code'], axis = 1)\n",
    "\n",
    "# Blue:\n",
    "df_cd['Diocese_Code'] = df_cd.apply(lambda row: 16 if row['Diocese'] in [\n",
    "    'Davenport', 'Tulsa', 'San Diego', \n",
    "    'Biloxi'] else row['Diocese_Code'], axis = 1)\n",
    "\n",
    "# Light green:\n",
    "df_cd['Diocese_Code'] = df_cd.apply(lambda row: 32 if row['Diocese'] in [\n",
    "    'Sacramento', 'Washington', 'Palm Beach', 'Birmingham',\n",
    "    'Rockville Centre'] else row['Diocese_Code'], axis = 1)\n",
    "\n",
    "# Yellow:\n",
    "df_cd['Diocese_Code'] = df_cd.apply(lambda row: 160 if row['Diocese'] in [\n",
    "    'Lubbock'] else row['Diocese_Code'], axis = 1)\n",
    "\n",
    "# Light orange:\n",
    "df_cd['Diocese_Code'] = df_cd.apply(lambda row: 96 if row['Diocese'] in [\n",
    "    'Youngstown', 'Denver'] else row['Diocese_Code'], axis = 1)\n",
    "\n",
    "# Orange:\n",
    "df_cd['Diocese_Code'] = df_cd.apply(lambda row: 112 if row['Diocese'] in [\n",
    "    'Alexandria in Louisiana', 'Covington', \n",
    "    'Raleigh'] else row['Diocese_Code'], axis = 1)\n",
    "\n",
    "# Purple:\n",
    "df_cd['Diocese_Code'] = df_cd.apply(lambda row: 144 if row['Diocese'] in [\n",
    "    'Lafayette in Louisiana', 'Lexington', 'Gary', 'New York', \n",
    "    'St. Louis'] else row['Diocese_Code'], axis = 1)\n",
    "\n",
    "# Brown:\n",
    "df_cd['Diocese_Code'] = df_cd.apply(lambda row: 174 if row['Diocese'] in [\n",
    "    'Louisville', 'Bridgeport', 'Milwaukee'] else row['Diocese_Code'], axis = 1)\n",
    "\n",
    "# Light purple:\n",
    "df_cd['Diocese_Code'] = df_cd.apply(lambda row: 128 if row['Diocese'] in [\n",
    "    'Saint Cloud', 'Memphis'] else row['Diocese_Code'], axis = 1)\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_cd.to_csv('counties_by_diocese.csv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "diocese_map = folium.Map(location=[38.7, -95], zoom_start=6, tiles = 'Stamen Toner')\n",
    "\n",
    "# Although Folium has a choropleth library, I wasn't able to find a way\n",
    "# to disable the default legend. Therefore, I am instead using a custom\n",
    "# choropleth mapping function. Much of this function is based on \n",
    "# Amodiovalerio Verde's code at:\n",
    "# https://vverde.github.io/blob/interactivechoropleth.html . \n",
    "\n",
    "colormap = branca.colormap.linear.Paired_12.to_step(\n",
    "    12, data = df_cd['Diocese_Code'])\n",
    "colormap\n",
    "\n",
    "tooltip = folium.features.GeoJsonTooltip(fields = \n",
    "['county_state_column', 'Diocese'], aliases = \n",
    "['County', 'Diocese'])\n",
    "# Based on https://python-visualization.github.io/folium/modules.html#folium.features.GeoJsonTooltip\n",
    "\n",
    "\n",
    "style_function = lambda x: {\n",
    "        \"fillOpacity\": 0.5,\n",
    "        \"weight\": 1,\n",
    "        'color': 'black',\n",
    "        \"fillColor\": colormap(x['properties']['Diocese_Code']),\n",
    "        # Note that 'properties' has to be added before the column name.\n",
    "        'fillOpacity': 0.8}\n",
    "\n",
    "folium.features.GeoJson(df_cd, tooltip = tooltip, \n",
    "style_function=style_function).add_to(diocese_map)\n",
    "\n",
    "diocese_map.save('diocese_map.html')\n",
    "\n",
    "# colormap.add_to(m)\n",
    "\n",
    "diocese_map"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "colormap"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "color_code_list = list(np.arange(0, 177, 176/11))\n",
    "color_code_list"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "color_list = list(np.arange(0, 177, 176/11))\n",
    "[colormap(color_list[i]) for i in range(len(color_list))]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "def create_map_screenshot(absolute_path_to_map_folder, map_name, \n",
    "screenshot_save_path = None):\n",
    "\n",
    "    '''\n",
    "    This function uses the Selenium library to create a screenshot \n",
    "    of a map so that it can be shared as a .png file.\n",
    "    See https://www.selenium.dev/documentation/ for more information on \n",
    "    Selenium. \n",
    "    \n",
    "    absolute_path_to_map_folder designates the absolute path where\n",
    "    the map is stored. (I wasn't able to get this code to work using\n",
    "    just relative paths.) \n",
    "\n",
    "    map_name specifies the name of the map, including its extension.\n",
    "\n",
    "    screenshot_save_path designates the folder where you wish to save\n",
    "    the map screenshot. This can be a relative path.\n",
    "\n",
    "    Note that some setup work is required for the Selenium code\n",
    "    to run correctly; if you don't have time right now to complete this \n",
    "    setup, you can comment out any code that calls this function.\n",
    "    '''\n",
    "\n",
    "    ff_driver = webdriver.Firefox() \n",
    "    # See https://www.selenium.dev/documentation/webdriver/getting_started/open_browser/\n",
    "    # For more information on using Selenium to get screenshots of .html \n",
    "    # files, see my get_screenshots.ipynb file within my route_maps_builder\n",
    "    # program, available here:\n",
    "    # https://github.com/kburchfiel/route_maps_builder/blob/master/get_screenshots.ipynb\n",
    "    window_width = 3000 # This produces a large window that can better\n",
    "    # capture small details (such as zip code shapefiles).\n",
    "    ff_driver.set_window_size(window_width,window_width*(9/16)) # Creates\n",
    "    # a window with an HD/4K/8K aspect ratio\n",
    "    ff_driver.get(f'{absolute_path_to_map_folder}\\\\{map_name}') \n",
    "    # See https://www.selenium.dev/documentation/webdriver/browser/navigation/\n",
    "    time.sleep(2) # This gives the page sufficient\n",
    "    # time to load the map tiles before the screenshot is taken. \n",
    "    # You can also experiment with longer sleep times.\n",
    "\n",
    "    if screenshot_save_path != None:\n",
    "        # If specifying a screenshot save path, you must create this path\n",
    "        # within your directory before the function is run; otherwise,\n",
    "        # it won't return an image. \n",
    "        ff_driver.get_screenshot_as_file(\n",
    "            screenshot_save_path+'/'+map_name.replace('.html','')+'.png') \n",
    "    else: # If no save path was specified for the screenshot, the image\n",
    "        # will be saved within the project's root folder.\n",
    "        ff_driver.get_screenshot_as_file(\n",
    "            map_name.replace('.html','')+'.png') \n",
    "    # Based on:\n",
    "    # https://www.selenium.dev/selenium/docs/api/java/org/openqa/selenium/TakesScreenshot.html\n",
    "\n",
    "    ff_driver.quit()\n",
    "    # Based on: https://www.selenium.dev/documentation/webdriver/browser/windows/\n",
    "\n",
    "def convert_png_to_smaller_jpg(png_folder, png_image_name, jpg_folder, \n",
    "reduction_factor = 1, quality_factor = 50):\n",
    "    ''' This function converts a .png image into a smaller .jpg image, which\n",
    "    helps reduce file sizes and load times when displaying a series of images\n",
    "    within a notebook or online.\n",
    "    png_folder and png_image_name specify the location of the original .png\n",
    "    image.\n",
    "    jpg folder specifies the location where the .jpg screenshot should be\n",
    "    saved.\n",
    "    reduction_factor specifies the amount by which you would like to reduce\n",
    "    the image's dimensions. For instance, to convert a 4K (3840*2160) image\n",
    "    to a full HD (1920*1080) one, use a reduction factor of 2. If you do not\n",
    "    wish to reduce the image's size, use the default reduction factor of 1.\n",
    "    '''\n",
    "    with PIL.Image.open(f'{png_folder}/{png_image_name}') as map_image:\n",
    "        (width, height) = (map_image.width // reduction_factor, \n",
    "        map_image.height // reduction_factor)\n",
    "        jpg_image = map_image.resize((width, height))\n",
    "        # The above code is based on:\n",
    "        # https://pillow.readthedocs.io/en/stable/reference/Image.html\n",
    "        jpg_image = jpg_image.convert('RGB')\n",
    "        # The above conversion is necessary in order to save .png files as \n",
    "        # .jpg files. It's based on Patrick Artner's answer at: \n",
    "        # https://stackoverflow.com/a/48248432/13097194\n",
    "        jpg_image_name = png_image_name.replace('png', 'jpg') \n",
    "        jpg_image.save(f'{jpg_folder}/{jpg_image_name}',\n",
    "        format = 'JPEG', quality = quality_factor, optimize = True)\n",
    "        # See https://pillow.readthedocs.io/en/stable/handbook/image-file-formats.html#jpeg"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "create_map_screenshot(\n",
    "    absolute_path_to_map_folder = \n",
    "    r'C:\\Users\\kburc\\D1V1\\Documents\\!Dell64docs\\Programming\\py\\kjb3_programs\\us_diocese_mapper', \n",
    "    map_name = 'diocese_map.html', screenshot_save_path = 'screenshots')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "convert_png_to_smaller_jpg(png_folder = 'screenshots', \n",
    "png_image_name = 'diocese_map.png', jpg_folder = 'smaller_screenshots', \n",
    "quality_factor = 80)\n",
    "\n",
    "IPython.display.Image(filename='smaller_screenshots/'+'diocese_map.jpg') \n",
    "# Based on https://ipython.readthedocs.io/en/stable/api/generated/IPython.display.html"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "end_time = time.time()\n",
    "run_time = end_time - start_time\n",
    "run_minutes = run_time // 60\n",
    "run_seconds = run_time % 60\n",
    "print(\"Completed run at\",time.ctime(end_time),\"(local time)\")\n",
    "print(\"Total run time:\",'{:.2f}'.format(run_time),\n",
    "\"second(s) (\"+str(run_minutes),\"minute(s) and\",'{:.2f}'.format(run_seconds),\n",
    "\"second(s))\") \n",
    "# Only meaningful when the program is run nonstop from start to finish"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "folium.Map(location=[38.7, -95], zoom_start=6, tiles = 'Stamen Toner')"
   ]
  }
 ],
 "metadata": {
  "interpreter": {
   "hash": "373fc5a02dc1a4f4d7122438e12fcda035e925b32ee1fe85dbff58fb3bdd1371"
  },
  "kernelspec": {
   "display_name": "Python 3.10.4 ('ga15pyd')",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.4"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}