us_diocese_mapper_V9.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",
    "I plan to add in more documentation/comments in the near future."
   ]
  },
  {
   "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": "markdown",
   "metadata": {},
   "source": [
    "Note: I used the diocese names on this list as the official names of almost all dioceses: https://www.usccb.org/about/bishops-and-dioceses/all-dioceses\n",
    "There are two exceptions: I use \"Springfield in Massachusetts\" instead of \"Springfield, Massachusetts\" and \"St. Thomas\" instead of \"St. Thomas, VI.\"\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",
    "from selenium.webdriver.common.keys import Keys \n",
    "# Source: https://selenium-python.readthedocs.io/getting-started.html\n",
    "import PIL.Image\n",
    "import IPython"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "generate_county_list_and_boundaries = False\n",
    "extract_FIPS_codes_from_website = False # Useful if working offline (or if \n",
    "# the website's layout changes in the future)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "if generate_county_list_and_boundaries == True:\n",
    "\n",
    "    shapefile_path = \\\n",
    "        r'C:\\Users\\kburc\\Downloads\\tl_2021_us_county_download\\tl_2021_us_county.shp'\n",
    "    if extract_FIPS_codes_from_website == True:\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', index = False)\n",
    "    else:\n",
    "        state_fips_codes = pd.read_csv('state_fips_codes.csv')\n",
    "\n",
    "    county_shapes = geopandas.read_file(shapefile_path)\n",
    "\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",
    "    # 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",
    "\n",
    "    # merged_county_shapes.to_file('original_county_outlines.geojson')\n",
    "    # This file is over 200 MB in size.\n",
    "\n",
    "    df_dioceses = merged_county_shapes.copy()\n",
    "\n",
    "else:\n",
    "    df_dioceses = geopandas.read_file('counties_by_diocese.geojson')\n",
    "    diocese_boundaries = geopandas.read_file('diocese_boundaries.geojson')\n",
    "    province_boundaries = geopandas.read_file('province_boundaries.geojson')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "regions_not_in_us_episcopal_conference = ['Puerto Rico', 'Guam', 'Northern Mariana Islands', 'American Samoa']\n",
    "df_dioceses = df_dioceses.query(\n",
    "    \"State_Name not in @regions_not_in_us_episcopal_conference\").copy() \n",
    "    # Puerto Rico and American Samoa belong to different episcopal conference, so they won't be included in this project (at least for now)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses['Diocese'] = 'Not_Set'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "condlist = [\n",
    "df_dioceses['State_Name'] == 'Alabama',\n",
    "df_dioceses['State_Name'] == 'Alaska',\n",
    "df_dioceses['State_Name'] == 'Arizona',\n",
    "df_dioceses['State_Name'] == 'Arkansas',\n",
    "df_dioceses['State_Name'] == 'California',\n",
    "df_dioceses['State_Name'] == 'Colorado',\n",
    "df_dioceses['State_Name'] == 'Connecticut',\n",
    "df_dioceses['State_Name'] == 'Delaware',\n",
    "df_dioceses['State_Name'] == 'Florida',\n",
    "df_dioceses['State_Name'] == 'Georgia',\n",
    "df_dioceses['State_Name'] == 'Hawaii',\n",
    "df_dioceses['State_Name'] == 'Idaho',\n",
    "df_dioceses['State_Name'] == 'Illinois',\n",
    "df_dioceses['State_Name'] == 'Indiana',\n",
    "df_dioceses['State_Name'] == 'Iowa',\n",
    "df_dioceses['State_Name'] == 'Kansas',\n",
    "df_dioceses['State_Name'] == 'Kentucky',\n",
    "df_dioceses['State_Name'] == 'Louisiana',\n",
    "df_dioceses['State_Name'] == 'Maine',\n",
    "df_dioceses['State_Name'] == 'Maryland',\n",
    "df_dioceses['State_Name'] == 'Massachusetts',\n",
    "df_dioceses['State_Name'] == 'Michigan',\n",
    "df_dioceses['State_Name'] == 'Minnesota',\n",
    "df_dioceses['State_Name'] == 'Mississippi',\n",
    "df_dioceses['State_Name'] == 'Missouri',\n",
    "df_dioceses['State_Name'] == 'Montana',\n",
    "df_dioceses['State_Name'] == 'Nebraska',\n",
    "df_dioceses['State_Name'] == 'Nevada',\n",
    "df_dioceses['State_Name'] == 'New Hampshire',\n",
    "df_dioceses['State_Name'] == 'New Jersey',\n",
    "df_dioceses['State_Name'] == 'New Mexico',\n",
    "df_dioceses['State_Name'] == 'New York',\n",
    "df_dioceses['State_Name'] == 'North Carolina',\n",
    "df_dioceses['State_Name'] == 'North Dakota',\n",
    "df_dioceses['State_Name'] == 'Ohio',\n",
    "df_dioceses['State_Name'] == 'Oklahoma',\n",
    "df_dioceses['State_Name'] == 'Oregon',\n",
    "df_dioceses['State_Name'] == 'Pennsylvania',\n",
    "df_dioceses['State_Name'] == 'Rhode Island',\n",
    "df_dioceses['State_Name'] == 'South Carolina',\n",
    "df_dioceses['State_Name'] == 'South Dakota',\n",
    "df_dioceses['State_Name'] == 'Tennessee',\n",
    "df_dioceses['State_Name'] == 'Texas',\n",
    "df_dioceses['State_Name'] == 'Utah',\n",
    "df_dioceses['State_Name'] == 'Vermont',\n",
    "df_dioceses['State_Name'] == 'Virginia',\n",
    "df_dioceses['State_Name'] == 'Washington',\n",
    "df_dioceses['State_Name'] == 'West Virginia',\n",
    "df_dioceses['State_Name'] == 'Wisconsin',\n",
    "df_dioceses['State_Name'] == 'Wyoming',\n",
    "df_dioceses['State_Name'] == 'District of Columbia',\n",
    "df_dioceses['State_Name'] == 'Virgin Islands'\n",
    "]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "choicelist = [\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Alabama',\n",
    "'Mobile',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Alaska',\n",
    "'Fairbanks',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Arizona',\n",
    "'Phoenix',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Arkansas',\n",
    "'Little Rock', # Statewide diocese\n",
    "\n",
    "# df_dioceses['State_Name'] == 'California',\n",
    "'Los Angeles', \n",
    "\n",
    "# df_dioceses['State_Name'] == 'Colorado',\n",
    "'Pueblo',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Connecticut',\n",
    "'Norwich', \n",
    "\n",
    "# df_dioceses['State_Name'] == 'Delaware',\n",
    "'Wilmington', # Statewide diocese\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Florida',\n",
    "'Orlando', \n",
    "\n",
    "# df_dioceses['State_Name'] == 'Georgia',\n",
    "'Savannah',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Hawaii',\n",
    "'Honolulu', # Statewide diocese\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Idaho',\n",
    "'Boise', #Statewide diocese\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Illinois',\n",
    "'Peoria',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Indiana',\n",
    "'Indianapolis', \n",
    "\n",
    "# df_dioceses['State_Name'] == 'Iowa',\n",
    "'Dubuque',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Kansas',\n",
    "'Salina',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Kentucky',\n",
    "'Lexington',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Louisiana',\n",
    "'Shreveport',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Maine',\n",
    "'Portland in Maine', # Statewide diocese\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Maryland',\n",
    "'Wilmington', \n",
    "\n",
    "# df_dioceses['State_Name'] == 'Massachusetts',\n",
    "'Boston',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Michigan',\n",
    "'Marquette',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Minnesota',\n",
    "'Duluth',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Mississippi',\n",
    "'Jackson',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Missouri',\n",
    "'Springfield-Cape Girardeau',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Montana',\n",
    "'Helena',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Nebraska',\n",
    "'Grand Island',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Nevada',\n",
    "'Reno',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'New Hampshire',\n",
    "'Manchester', # Statewide diocese\n",
    "\n",
    "# df_dioceses['State_Name'] == 'New Jersey',\n",
    "'Camden',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'New Mexico',\n",
    "'Santa Fe',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'New York',\n",
    "'Albany',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'North Carolina',\n",
    "'Raleigh',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'North Dakota',\n",
    "'Fargo',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Ohio',\n",
    "'Columbus',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Oklahoma',\n",
    "'Oklahoma City',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Oregon',\n",
    "'Portland in Oregon',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Pennsylvania',\n",
    "'Erie',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Rhode Island',\n",
    "'Providence', #Statewide diocese\n",
    "\n",
    "# df_dioceses['State_Name'] == 'South Carolina',\n",
    "'Charleston', # Statewide diocese\n",
    "\n",
    "# df_dioceses['State_Name'] == 'South Dakota',\n",
    "'Sioux Falls',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Tennessee',\n",
    "'Nashville',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Texas',\n",
    "'Lubbock',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Utah',\n",
    "'Salt Lake City', # Statewide Diocese\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Vermont',\n",
    "'Burlington', # Statewide diocese\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Virginia',\n",
    "'Richmond',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Washington',\n",
    "'Seattle',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'West Virginia',\n",
    "'Wheeling-Charleston', # Statewide diocese\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Wisconsin',\n",
    "'La Crosse',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Wyoming'],\n",
    "'Cheyenne', # Statewide diocese\n",
    "\n",
    "# df_dioceses['State_Name'] == 'District of Columbia',\n",
    "'Washington',\n",
    "\n",
    "# df_dioceses['State_Name'] == 'Virgin Islands'\n",
    "'St. Thomas'\n",
    "]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses['Diocese'] = np.select(condlist, choicelist, df_dioceses['Diocese'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "len(df_dioceses)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "initial_version_of_df_dioceses = \\\n",
    "df_dioceses.copy()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Drop counties that need to be divided here, and then add in custom county shapefiles to df_dioceses. (Use geojson.io to create these.)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses.set_index('county_state', drop = False, inplace = True)\n",
    "\n",
    "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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses.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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses.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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "fresno_county_list, state = 'California', diocese = 'Fresno')\n",
    "\n",
    "# Diocese of Monterey\n",
    "monterey_in_california_county_list = [\n",
    "'Monterey', 'San Benito County', 'San Luis Obispo County', 'Santa Cruz County'\n",
    "]\n",
    "\n",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "monterey_in_california_county_list, state = 'California', \n",
    "diocese = 'Monterey')\n",
    "\n",
    "# Diocese of Orange:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Orange\n",
    "df_dioceses.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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "sacramento_county_list, state = 'California', \n",
    "diocese = 'Sacramento')\n",
    "\n",
    "# Diocese of San Jose:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_San_Jose_in_California\n",
    "\n",
    "df_dioceses.at[\n",
    "\"Santa Clara County, California\", 'Diocese'] = 'San Jose' \n",
    "\n",
    "# Diocese of Santa Rosa:\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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "santa_rosa_in_california_county_list, state = 'California', \n",
    "diocese = 'Santa Rosa')\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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "palm_beach_county_list, state = 'Florida', \n",
    "diocese = 'Palm Beach')\n",
    "\n",
    "# Diocese of Venice:\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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "pensacola_tallahassee_county_list, state = 'Florida', \n",
    "diocese = 'Pensacola-Tallahassee')\n",
    "\n",
    "\n",
    "# Diocese of St. 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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "saint_petersburg_county_list, state = 'Florida', \n",
    "diocese = 'St. 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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "lake_charles_county_list, state = 'Louisiana', \n",
    "diocese = 'Lake Charles')\n",
    "\n",
    "# Diocese of Alexandria:\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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "alexandria_in_louisiana_county_list, state = 'Louisiana', \n",
    "diocese = 'Alexandria')\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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "baker_county_list, state = 'Oregon', \n",
    "diocese = 'Baker')\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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "jefferson_city_county_list, state = 'Missouri', \n",
    "diocese = 'Jefferson City')\n",
    "\n",
    "# Archdiocese of St. 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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "saint_paul_and_minneapolis_county_list, state = 'Minnesota', \n",
    "diocese = 'St. 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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "new_ulm_county_list, state = 'Minnesota', \n",
    "diocese = 'New Ulm')\n",
    "\n",
    "# Diocese of St. 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",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "saint_cloud_county_list, state = 'Minnesota', \n",
    "diocese = 'St. 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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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",
    "df_dioceses = add_dioceses(df_dioceses, \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\n",
    "\n",
    "\n",
    "df_dioceses.reset_index(drop=True, inplace = True)"
   ]
  },
  {
   "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(df_dioceses)) \n",
    "if df_dioceses.iloc[i, df_dioceses.columns.get_loc(\n",
    "    'Diocese')] != initial_version_of_df_dioceses.iloc[\n",
    "        i, initial_version_of_df_dioceses.columns.get_loc(\n",
    "            'Diocese')]]\n",
    "len(counties_changed) # The number of counties that were modified above"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "if 'Province' in df_dioceses.columns:\n",
    "    df_dioceses.drop('Province', axis = 1, inplace = True)\n",
    "    # This prevents two 'Province' columns from existing within df_dioceses.\n",
    "\n",
    "df_provinces = pd.read_csv('diocese_province_list.csv')\n",
    "df_dioceses = df_dioceses.merge(df_provinces, on = 'Diocese', how = 'left')\n",
    "df_dioceses.drop_duplicates('Diocese').sort_values('Province').tail(30)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "cathedral_list = pd.read_html('https://en.wikipedia.org/wiki/List_of_Catholic_cathedrals_in_the_United_States')[0][['Diocese', 'Cathedral']]\n",
    "cathedral_list['Diocese'] = cathedral_list['Diocese'].str.replace('Diocese of', '').str.replace('Archdiocese of', '').str.strip()\n",
    "# Adding str.strip() ended up being necessary for the merge to work\n",
    "cathedral_list['Cathedral_List'] = [[cathedral_list['Cathedral'][i]] for i in range(len(cathedral_list))]\n",
    "cathedral_list.drop('Cathedral', axis = 1, inplace = True)\n",
    "cathedral_list['Cathedral_coord_list'] = [[] for i in range(len(cathedral_list))]\n",
    "df_provinces_cathedrals = df_provinces.merge(cathedral_list, on = 'Diocese', how = 'left')\n",
    "df_provinces_cathedrals.to_csv('diocese_province_cathedral_list.csv', index = False)\n",
    "df_provinces_cathedrals"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "I then went to openstreetmap.org to determine the coordinates of each cathedral. I found these coordinates by searching for each cathedral; right clicking on the cathedral that appeared on the map and selecting 'show address'; and then copying in and pasting the geographic coordinates that appeared within the 'Results from Internal' section. (Alternately, I could have selected 'Centre Map Here' and then copied and pasted in the coordinates from the URL box.)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_provinces"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses.groupby('Diocese').count().sort_values(['STATEFP'], ascending = False).head(15)['COUNTYFP']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses.groupby('Diocese').count().sort_values(['STATEFP']).head(5)['COUNTYFP']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses.set_index('county_state', drop = False, inplace = True)\n",
    "df_dioceses['Diocese_Detail'] = df_dioceses['Diocese']\n",
    "df_dioceses['Province_Detail'] = df_dioceses['Province']\n",
    "\n",
    "df_dioceses.at['Plymouth County, Massachusetts', 'Diocese_Detail'] = 'Boston, Fall River'\n",
    "df_dioceses.at['Fayette County, Texas', 'Diocese_Detail'] = 'Austin, Victoria'\n",
    "df_dioceses.at['McMullen County, Texas', 'Diocese_Detail'] = 'Corpus Christi, San Antonio'\n",
    "df_dioceses.at['McMullen County, Texas', 'Province_Detail'] = 'Galveston-Houston, San Antonio'\n",
    "\n",
    "df_dioceses.at['Spencer County, Indiana', 'Diocese_Detail'] = 'Evansville, Indianapolis'\n",
    "df_dioceses.at['Rio Arriba County, New Mexico', 'Diocese_Detail'] = 'Gallup, Santa Fe'\n",
    "df_dioceses.at['Sandoval County, New Mexico', 'Diocese_Detail'] = 'Gallup, Santa Fe'\n",
    "df_dioceses.at['Bernalillo County, New Mexico', 'Diocese_Detail'] = 'Gallup, Santa Fe'\n",
    "df_dioceses.at['Valencia County, New Mexico', 'Diocese_Detail'] = 'Gallup, Santa Fe'\n",
    "df_dioceses.at['Pinal County, Arizona', 'Diocese_Detail'] = 'Phoenix, Tucson'\n",
    "df_dioceses.at['Coconino County, Arizona', 'Diocese_Detail'] = 'Gallup, Phoenix'\n",
    "df_dioceses.at['Jefferson Parish, Louisiana', 'Diocese_Detail'] = 'Houma-Thibodaux, New Orleans'\n",
    "df_dioceses.at['St. Mary Parish, Louisiana', 'Diocese_Detail'] = 'Houma-Thibodaux, Lafayette in Louisiana'\n",
    "df_dioceses.at['Suffolk County, New York', 'Diocese_Detail'] = 'Norwich, Rockville Centre'\n",
    "df_dioceses.at['Herkimer County, New York', 'Diocese_Detail'] = 'Albany, Ogdensburg'\n",
    "df_dioceses.at['Hamilton County, New York', 'Diocese_Detail'] = 'Albany, Ogdensburg'\n",
    "df_dioceses.at['Deuel County, Nebraska', 'Diocese_Detail'] = 'Grand Island, Lincoln'\n",
    "df_dioceses.at['Keith County, Nebraska', 'Diocese_Detail'] = 'Grand Island, Lincoln'\n",
    "df_dioceses.at['Lincoln County, Nebraska', 'Diocese_Detail'] = 'Grand Island, Lincoln'\n",
    "df_dioceses.at['Dawson County, Nebraska', 'Diocese_Detail'] = 'Grand Island, Lincoln'\n",
    "df_dioceses.at['Hall County, Nebraska', 'Diocese_Detail'] = 'Grand Island, Lincoln'\n",
    "df_dioceses.at['Toole County, Montana', 'Diocese_Detail'] = 'Great Falls-Billings, Helena'\n",
    "df_dioceses.at['Teton County, Montana', 'Diocese_Detail'] = 'Great Falls-Billings, Helena'\n",
    "df_dioceses.at['Bottineau County, North Dakota', 'Diocese_Detail'] = 'Bismark, Fargo'\n",
    "\n",
    "df_dioceses.reset_index(drop=True, inplace = True)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "if generate_county_list_and_boundaries == True:\n",
    "    # See https://geopandas.org/en/stable/docs/user_guide/aggregation_with_dissolve.html\n",
    "    diocese_boundaries = df_dioceses.dissolve(by = 'Diocese')\n",
    "    diocese_boundaries['geometry'] = diocese_boundaries.simplify(tolerance = 0.005)\n",
    "    diocese_boundaries['Diocese'] = diocese_boundaries.index \n",
    "    diocese_boundaries.reset_index(drop=True,inplace=True)\n",
    "    diocese_boundaries.to_file('diocese_boundaries.geojson', driver = 'GeoJSON')\n",
    "    diocese_boundaries.head(5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "if generate_county_list_and_boundaries == True:\n",
    "    # See https://geopandas.org/en/stable/docs/user_guide/aggregation_with_dissolve.html\n",
    "    province_boundaries = df_dioceses.dissolve(by = 'Province')\n",
    "    province_boundaries['geometry'] = province_boundaries.simplify(tolerance = 0.005) \n",
    "    province_boundaries.to_file('province_boundaries.geojson', driver = 'GeoJSON')\n",
    "    province_boundaries.head(5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses['geometry'] = df_dioceses.simplify(tolerance = 0.005) "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses.to_file('counties_by_diocese.geojson', driver = 'GeoJSON')\n",
    "# This raises a ValueError if the index is set to county_state, though I'm not sure why."
   ]
  },
  {
   "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. [Revised diocese maps]\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",
    "[Revised diocese maps]\n",
    "\n",
    "3. San Antonio/Corpus Christi: The Archdiocese of San Antonio includes only that part of McMullen County north of the Nueces River. The rest belongs to the Diocese of Corpus Christi. See: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_San_Antonio\n",
    "This also means that McMullen County is part of two different provinces (Galveston-Houston and San Antonio). I believe it is the only county to have this distinction. [Revised diocese and province maps]\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",
    "[Revised diocese and province maps]\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 [-106.87806 decimal degrees longitude] in New Mexico.\" The rest of those counties belong to Santa Fe. (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Gallup)\n",
    "(The Wikipedia page references https://web.archive.org/web/20120624033053/https://dioceseofgallup.org/about_history.php)\n",
    "(Perhaps it's related to the New Mexico meridian? https://en.wikipedia.org/wiki/New_Mexico_meridian That one is -106.8944 decimal degrees longitude).\n",
    "[Revised diocese maps]\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",
    "(This map of the Gila River Indian Community was also helpful: https://en.wikipedia.org/wiki/Gila_River_Indian_Community#/media/File:Pinal_County_Arizona_Incorporated_and_Unincorporated_areas_GRIC_highlighted.svg )\n",
    "[Revised diocese maps]\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",
    "[Revised diocese maps]\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",
    "Based on parish listings for each diocese (https://htdiocese.org/church-parishes-by-area and https://diolaf.org/all-parishes ), I allocated Berwick and Patterson to the Diocese of Lafayette. I ended up making the Atchafalaya River the dividing point between these two dioceses within St. Mary's. \n",
    "[Revised diocese maps]\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",
    "[Revised diocese maps]\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",
    "(https://www.rcdony.org/about/history.html and https://evangelist.org/Images/Images/2470.jpg were helpful in revising the diocesan borders.) [Revised diocese maps]\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) [Revised diocese maps]\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). [Revised diocese maps]\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",
    "[Revised diocese maps]\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. This change is not reflected on the current Bismark and Fargo diocese maps on Wikipedia. \n",
    "[Revised diocese maps]\n",
    "\n",
    "16. Gallup/Phoenix: The Navajo Nation part of Coconino County belongs to the Diocese of Gallup; the rest of Coconino County belongs to the Diocese of Phoenix. (This is reflected on the Phoenix and Santa Fe Province maps on Wikipedia but not on the Gallup map).\n",
    "[Revised diocese maps]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses.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",
    "    df_dioceses['Diocese']))))\n",
    "df_diocese_numbers.columns = ['Diocese_Color_Code', 'Diocese']\n",
    "df_diocese_numbers"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses.head(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_diocese_numbers"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "if 'Diocese_Color_Code' in df_dioceses:\n",
    "    df_dioceses.drop('Diocese_Color_Code', axis = 1, inplace = True)\n",
    "df_dioceses = df_dioceses.merge(df_diocese_numbers, on = 'Diocese', \n",
    "how = 'left') \n",
    "df_dioceses"
   ]
  },
  {
   "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_dioceses['Diocese_Color_Code'] = df_dioceses.apply(lambda row: 48 if row['Diocese'] in \n",
    "['Seattle', 'Dubuque', 'Charleston'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Light blue:\n",
    "df_dioceses['Diocese_Color_Code'] = df_dioceses.apply(lambda row: 0 if row['Diocese'] in [\n",
    "    'Santa Rosa', 'Crookston', 'Kansas City-St. Joseph', 'St. Petersburg', \n",
    "    'Gaylord'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "#Pink:\n",
    "df_dioceses['Diocese_Color_Code'] = df_dioceses.apply(lambda row: 64 if row['Diocese'] in [\n",
    "    'Reno', 'San Antonio', 'Cleveland'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "#Red:\n",
    "df_dioceses['Diocese_Color_Code'] = df_dioceses.apply(lambda row: 80 if row['Diocese'] in [\n",
    "    'Cheyenne', 'Sioux Falls', 'Baltimore', 'San Francisco', 'Scranton']\n",
    "    else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Blue:\n",
    "df_dioceses['Diocese_Color_Code'] = df_dioceses.apply(lambda row: 16 if row['Diocese'] in [\n",
    "    'Davenport', 'Tulsa', 'San Diego', \n",
    "    'Biloxi', 'St. Paul and Minneapolis'\n",
    "    ] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Light green:\n",
    "df_dioceses['Diocese_Color_Code'] = df_dioceses.apply(lambda row: 32 if row['Diocese'] in [\n",
    "    'Sacramento', 'Washington', 'Palm Beach', 'Birmingham', 'Paterson',\n",
    "    'Rockville Centre'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Yellow (changed to dark red within the mapping code to make the labels\n",
    "# easier to read):\n",
    "df_dioceses['Diocese_Color_Code'] = df_dioceses.apply(lambda row: 160 if row['Diocese'] in [\n",
    "    'Lubbock'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Light orange:\n",
    "df_dioceses['Diocese_Color_Code'] = df_dioceses.apply(lambda row: 96 if row['Diocese'] in [\n",
    "    'Youngstown', 'Denver'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Orange:\n",
    "df_dioceses['Diocese_Color_Code'] = df_dioceses.apply(lambda row: 112 if row['Diocese'] in [\n",
    "    'Alexandria', 'Covington', 'Helena', 'Raleigh'] \n",
    "    else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Purple:\n",
    "df_dioceses['Diocese_Color_Code'] = df_dioceses.apply(lambda row: 144 if row['Diocese'] in [\n",
    "    'Lafayette in Louisiana', 'Lexington', 'Gary', 'New York', \n",
    "    'St. Louis'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Brown:\n",
    "df_dioceses['Diocese_Color_Code'] = df_dioceses.apply(lambda row: 174 if row['Diocese'] in [\n",
    "    'Louisville', 'Bridgeport', 'Milwaukee', \n",
    "    'Colorado Springs'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Light purple:\n",
    "df_dioceses['Diocese_Color_Code'] = df_dioceses.apply(lambda row: 128 if row['Diocese'] in [\n",
    "    'St. Cloud', 'Memphis'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses.to_csv('counties_by_diocese.csv', index = False)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Finding mean latitude and longitude points for each diocese (in order to determine where to plot the diocese name)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "diocese_boundaries['center'] = diocese_boundaries['geometry'].representative_point()\n",
    "# See https://geopandas.org/en/stable/docs/reference/api/geopandas.GeoSeries.representative_point.html"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['Diocese']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "diocese_boundaries"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Portland in Oregon', 'Portland OR')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('San Francisco', 'SF')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('San Jose', 'SJ')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Oakland', 'Oakl.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Orange', 'Or.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Galveston-Houston', 'Galv.-Hou.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Victoria in Texas', 'Victoria')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('St. Paul and Minneapolis', 'St. Paul-Minn.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Chicago', 'Chi.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Lafayette in Indiana', 'Lafayette IN')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Lafayette in Louisiana', 'Lafyt. LA')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Shreveport', 'Shrvpt.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('St. Petersburg', 'St. Ptrsbg.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Springfield in Illinois', 'Springfield IL')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Washington', 'Was.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Baltimore', 'Balt.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Pittsburgh', 'Pitt.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Greensburg', 'Grnsbg.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Altoona-Johnstown', 'Alt-Jnstn.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Harrisburg', 'Hrsbg.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Allentown', 'Allntn.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Springfield in Massachusetts', 'Springfield MA.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Worcester', 'Wrcstr.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Bridgeport', 'Brgpt.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Providence', 'Prv.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Norwich', 'Nwch.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Brooklyn', 'Brkln.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Rockville Centre', 'Rockvl_Ctr.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Metuchen', 'Mtchn.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Newark', 'Nwrk.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Steubenville', 'Steub.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Philadelphia', 'Phil.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Paterson', 'Ptrsn.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Hartford', 'Htfd.')\n",
    "diocese_boundaries['name_for_plotting'] = diocese_boundaries['name_for_plotting'].str.replace('Manchester', 'Mnchstr.')\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "diocese_boundaries"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_provinces.query(\"Diocese == 'Newark'\")"
   ]
  },
  {
   "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_dioceses['Diocese_Color_Code'])\n",
    "colormap\n",
    "\n",
    "tooltip = folium.features.GeoJsonTooltip(fields = \n",
    "['county_state', 'Diocese', 'Province'], aliases = \n",
    "['County', 'Diocese', 'Province'])\n",
    "# Based on https://python-visualization.github.io/folium/modules.html#folium.features.GeoJsonTooltip\n",
    "\n",
    "\n",
    "style_function = lambda x: {\n",
    "        \"weight\": 0.5,\n",
    "        'color': '#777777',\n",
    "        \"fillColor\": '#770000' if\n",
    "        colormap(x['properties']['Diocese_Color_Code']) == '#ffff99ff'\n",
    "        else colormap(x['properties']['Diocese_Color_Code']),\n",
    "        # Note that 'properties' has to be added before the column name.\n",
    "        # This can be inferred from the GeoJson definition within Folium:\n",
    "        # https://python-visualization.github.io/folium/modules.html#folium.features.GeoJson\n",
    "        # The if/else statement checks whether the color corresponding to the\n",
    "        # colormap is yellow. (I found this value by printing out each color\n",
    "        # within the colormap, as shown later in the code.)\n",
    "        # If so, it changes the fill color to a dark red \n",
    "        # so that white text labels can be read atop it.\n",
    "        'fillOpacity': 0.8}\n",
    "\n",
    "folium.features.GeoJson(df_dioceses, \n",
    "style_function=style_function).add_to(diocese_map)\n",
    "\n",
    "\n",
    "\n",
    "diocese_style_function = lambda x: {\n",
    "        \"fillOpacity\": 0,\n",
    "        \"weight\": 2,\n",
    "        'color': '#FFFFFF',\n",
    "        'Fill': False}\n",
    "\n",
    "folium.features.GeoJson(diocese_boundaries['geometry'], \n",
    "style_function=diocese_style_function).add_to(diocese_map)\n",
    "# diocese_boundaries includes both geometry and point values. The point values\n",
    "# confuse the interpreter, so in order to avoid an error message, I \n",
    "# entered 'geometry' here to specify that those were the points to be plotted.\n",
    "\n",
    "\n",
    "province_style_function = lambda x: {\n",
    "        'weight': '2', \n",
    "        'fillOpacity': 0,\n",
    "        'color': \n",
    "        '#000000', \n",
    "        'Fill': False} \n",
    "# See https://python-visualization.github.io/folium/modules.html#folium.vector_layers.path_options\n",
    "# for style_function arguments.\n",
    "\n",
    "folium.GeoJson(province_boundaries, name = \"geojson\", style_function = province_style_function).add_to(diocese_map)\n",
    "\n",
    "tooltip_style_function = lambda x: {\n",
    "        \"fillOpacity\": 0,\n",
    "        \"weight\": 0,\n",
    "        'color': '#000000',\n",
    "        'Fill': False}\n",
    "\n",
    "folium.features.GeoJson(df_dioceses, tooltip = tooltip,\n",
    "style_function=tooltip_style_function).add_to(diocese_map)\n",
    "\n",
    "for i in range(len(diocese_boundaries)):\n",
    "    center_point = [float(diocese_boundaries.iloc[i, diocese_boundaries.columns.get_loc('center')].y),\n",
    "    float(diocese_boundaries.iloc[i, diocese_boundaries.columns.get_loc('center')].x)] \n",
    "    # The use of .y and .x here is based on Francois M.'s response at:\n",
    "    # https://gis.stackexchange.com/a/255867/195615\n",
    "    diocese_name = diocese_boundaries.iloc[i, diocese_boundaries.columns.get_loc('Diocese')]\n",
    "    diocese_name_for_plotting = diocese_boundaries.iloc[i, diocese_boundaries.columns.get_loc('name_for_plotting')]\n",
    "    split_name = diocese_name_for_plotting.replace('-', ' ').split(' ')\n",
    "    horizontal_offset = max([len(word) for word in split_name])*3\n",
    "    vertical_offset = len(split_name)*10\n",
    "#     print(split_name, horizontal_offset) # For debugging\n",
    "    style_tag = '<style> p{color:white; line-height: 1} </style>'\n",
    "#     style_tag = '<style> p{text-align:center} </style>'\n",
    "#     align_string = 'style = \"text-align:center\"'\n",
    "#     align_string = 'style = \"color:green\"'\n",
    "#     folium.CircleMarker(center_point, radius = 1).add_to(diocese_map) \n",
    "# The above line is very useful for debugging what offsets to use when placing\n",
    "# text on the map.\n",
    "    if diocese_name == diocese_boundaries.iloc[i, diocese_boundaries.columns.get_loc('Province')]: # E.g. if this is an archdiocese (as provinces are named after archdioceses)\n",
    "        folium.Marker(center_point, icon = folium.DivIcon(icon_anchor=(horizontal_offset, vertical_offset), html=f\"{style_tag}<p><b><i>{diocese_name_for_plotting}</i><b></p>\")).add_to(diocese_map)\n",
    "    else:\n",
    "        folium.Marker(center_point, icon = folium.DivIcon(icon_anchor=(horizontal_offset, vertical_offset), html=f\"{style_tag}<p><b>{diocese_name_for_plotting}</b></p>\")).add_to(diocese_map)\n",
    "        # folium.Marker(center_point, icon = folium.DivIcon(icon_anchor=(horizontal_offset, 10), html=f\"{style_tag}<b><p>{diocese_name_for_plotting.replace(' ', '_').replace('-', '_')}</p></b>\")).add_to(diocese_map)\n",
    "\n",
    "# The strategy of using an invisible folium.Marker object to display text on \n",
    "# the map via DivIcon elements comes from from Bob Haffner's response at:\n",
    "# https://stackoverflow.com/a/56878525/13097194 or:\n",
    "# and https://stackoverflow.com/a/46408144/13097194\n",
    "\n",
    "\n",
    "# folium.LatLngPopup().add_to(diocese_map) # Useful for debugging\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_list = list(np.arange(0, 177, 176/12))\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, added_text = '', \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','')+added_text+'.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','')+added_text+'.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')\n",
    "\n",
    "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": "markdown",
   "metadata": {},
   "source": [
    "# Generating a map with revised diocesan/provincial boundaries"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "revised_diocese_boundaries = geopandas.read_file('diocese_boundaries_revised.geojson')\n",
    "revised_diocese_boundaries.sort_values('Diocese', inplace = True) \n",
    "revised_diocese_boundaries = revised_diocese_boundaries.merge(df_dioceses[['Diocese', 'Diocese_Color_Code']].drop_duplicates(), on = 'Diocese')\n",
    "revised_diocese_boundaries = revised_diocese_boundaries.merge(diocese_boundaries[['Diocese', 'name_for_plotting']], on = 'Diocese')\n",
    "revised_diocese_boundaries\n",
    "revised_diocese_boundaries['center'] = revised_diocese_boundaries['geometry'].representative_point()\n",
    "revised_diocese_boundaries.head(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "revised_province_boundaries = geopandas.read_file('province_boundaries_revised.geojson')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_cathedral_locations = pd.read_csv('diocese_province_cathedral_list_revised.csv')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To obtain the geographic locations of each cathedral, I used multiple tools, including OpenStreetMap, WikiPedia, and GeoHack (which I believe uses the GPL license: https://bitbucket.org/magnusmanske/geohack/src/master/LICENSE). GeoHack proved to be the fastest option. I accessed each cathedral's GeoHack page by clicking on the coordinates listed on that cathedral's Wikipedia page. (I used GeoHack instead of the coordinates shown on Wikipedia because the former offered coordinates in decimal degree format.)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "def create_diocese_map(tile_option = 'Stamen Toner', starting_zoom = 6, starting_location = [38.7, -95], save_path = 'revised_diocese_map', include_cathedrals = True, include_names = True):    \n",
    "    revised_diocese_map = folium.Map(location=starting_location, zoom_start=starting_zoom, tiles = tile_option)\n",
    "    # tiles can be set to 'None' in order to produce a map without any background \n",
    "    # details. This produces a cleaner static image.\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 = revised_diocese_boundaries['Diocese_Color_Code'])\n",
    "    colormap\n",
    "\n",
    "    tooltip = folium.features.GeoJsonTooltip(fields = \n",
    "    ['county_state', 'Diocese_Detail', 'Province_Detail'], aliases = \n",
    "    ['County', 'Diocese', 'Province'])\n",
    "    # Based on https://python-visualization.github.io/folium/modules.html#folium.features.GeoJsonTooltip\n",
    "\n",
    "\n",
    "    style_function = lambda x: {\n",
    "            \"weight\": 0.5,\n",
    "            'color': '#000000',\n",
    "            # Note that 'properties' has to be added before the column name.\n",
    "            # This can be inferred from the GeoJson definition within Folium:\n",
    "            # https://python-visualization.github.io/folium/modules.html#folium.features.GeoJson\n",
    "            'fillOpacity': 0}\n",
    "\n",
    "    folium.features.GeoJson(df_dioceses, name = 'Counties',\n",
    "    style_function=style_function).add_to(revised_diocese_map)\n",
    "\n",
    "    diocese_style_function = lambda x: {\n",
    "            \"fillOpacity\": 0.8,\n",
    "            \"weight\": 2,\n",
    "            'color': '#FFFFFF',\n",
    "            \"fillColor\": '#770000' if\n",
    "            colormap(x['properties']['Diocese_Color_Code']) == '#ffff99ff'\n",
    "            else colormap(x['properties']['Diocese_Color_Code']),\n",
    "            'Fill': False}\n",
    "\n",
    "    folium.features.GeoJson(revised_diocese_boundaries[['geometry', 'Diocese_Color_Code']], \n",
    "    name = 'Dioceses', style_function=diocese_style_function).add_to(revised_diocese_map)\n",
    "    # revised_diocese_boundaries includes both geometry and point values. The point values\n",
    "    # confuse the interpreter, so in order to avoid an error message, I \n",
    "    # entered 'geometry' here to specify that those were the points to be plotted.\n",
    "\n",
    "\n",
    "    province_style_function = lambda x: {\n",
    "            'weight': '2', \n",
    "            'fillOpacity': 0,\n",
    "            'color': \n",
    "            '#000000', \n",
    "            'Fill': False} \n",
    "    # See https://python-visualization.github.io/folium/modules.html#folium.vector_layers.path_options\n",
    "    # for style_function arguments.\n",
    "\n",
    "\n",
    "    folium.GeoJson(revised_province_boundaries, name = \"Provinces\", style_function = province_style_function).add_to(revised_diocese_map)\n",
    "\n",
    "    # The tooltips are added after the diocese and province boundaries\n",
    "    # so that they will appear when the user\n",
    "    # hovers over them.\n",
    "\n",
    "    tooltip_style_function = lambda x: {\n",
    "            \"fillOpacity\": 0,\n",
    "            \"weight\": 0,\n",
    "            'color': '#000000',\n",
    "            'Fill': False}\n",
    "\n",
    "    # The following GeoJson function adds tooltips to the map so that users\n",
    "    # can see county, diocese, and province information when they hover over\n",
    "    # a given county.\n",
    "    folium.features.GeoJson(df_dioceses, tooltip = tooltip, name = 'Tooltips',\n",
    "    style_function=tooltip_style_function).add_to(revised_diocese_map)\n",
    "\n",
    "\n",
    "    if include_cathedrals == True:\n",
    "        icon_width = 10\n",
    "        icon_height = icon_width/(120/200)\n",
    "        # The original markers are 120 by 200 pixels, so the icon's width is divided by this ratio to retrieve the height.\n",
    "        cathedral_counter = 0\n",
    "        cathedral_feature_group = folium.FeatureGroup(name = 'Cathedrals')\n",
    "        for i in range(len(df_cathedral_locations)):\n",
    "            # print(i)\n",
    "            loc_list = df_cathedral_locations.iloc[i, df_cathedral_locations.columns.get_loc('Cathedral_coord_list')].split('],')\n",
    "            loc_list = [loc.replace('[','').replace(']','') for loc in loc_list]\n",
    "            loc_list = [loc.split(', ') for loc in loc_list]\n",
    "            loc_list = [[float(loc[0]), float(loc[1])] for loc in loc_list]\n",
    "\n",
    "            cathedral_list = df_cathedral_locations.iloc[i, df_cathedral_locations.columns.get_loc('Cathedral_List')].split(',')\n",
    "            cathedral_list = [cathedral.replace('[','').replace(']','').replace(\"'\", '') for cathedral in cathedral_list]\n",
    "            cathedral_list\n",
    "\n",
    "            diocese_name = df_cathedral_locations.iloc[i, df_cathedral_locations.columns.get_loc('Diocese')]\n",
    "            if diocese_name == df_cathedral_locations.iloc[i, df_cathedral_locations.columns.get_loc('Province')]:\n",
    "                # print(\"True\")\n",
    "                marker = 'globus_cruciger_archbishop.png'\n",
    "            else:\n",
    "                marker = 'globus_cruciger.png'\n",
    "            for j in range(len(loc_list)):\n",
    "                folium.Marker(loc_list[j], icon = folium.CustomIcon(marker, icon_size = (icon_width, icon_height)), tooltip = cathedral_list[j]).add_to(cathedral_feature_group)\n",
    "                cathedral_counter += 1\n",
    "        print(\"Added\",cathedral_counter,\"cathedrals to the map.\")\n",
    "\n",
    "        cathedral_feature_group.add_to(revised_diocese_map)\n",
    "\n",
    "\n",
    "\n",
    "    if include_names == True:\n",
    "        name_feature_group = folium.FeatureGroup(name = 'Diocese Names')\n",
    "        # This method of assigning markers to a FeatureGroup that can then be toggled\n",
    "        # on and off comes from Bob Haffner at\n",
    "        # https://stackoverflow.com/a/61283334/13097194 .\n",
    "\n",
    "        for i in range(len(revised_diocese_boundaries)):\n",
    "            center_point = [float(revised_diocese_boundaries.iloc[i, revised_diocese_boundaries.columns.get_loc('center')].y),\n",
    "            float(revised_diocese_boundaries.iloc[i, revised_diocese_boundaries.columns.get_loc('center')].x)] \n",
    "            # The use of .y and .x here is based on Francois M.'s response at:\n",
    "            # https://gis.stackexchange.com/a/255867/195615\n",
    "            diocese_name = revised_diocese_boundaries.iloc[i, revised_diocese_boundaries.columns.get_loc('Diocese')]\n",
    "            diocese_name_for_plotting = revised_diocese_boundaries.iloc[i, revised_diocese_boundaries.columns.get_loc('name_for_plotting')]\n",
    "            split_name = diocese_name_for_plotting.replace('-', ' ').split(' ')\n",
    "            horizontal_offset = max([len(word) for word in split_name])*3\n",
    "            vertical_offset = len(split_name)*10\n",
    "        #     print(split_name, horizontal_offset) # For debugging\n",
    "            style_tag = '<style> p{color:white; line-height: 1} </style>'\n",
    "        #     style_tag = '<style> p{text-align:center} </style>'\n",
    "        #     align_string = 'style = \"text-align:center\"'\n",
    "        #     align_string = 'style = \"color:green\"'\n",
    "        #     folium.CircleMarker(center_point, radius = 1).add_to(revised_diocese_map) \n",
    "        # The above line is very useful for debugging what offsets to use when placing\n",
    "        # text on the map.\n",
    "            if diocese_name == revised_diocese_boundaries.iloc[i, revised_diocese_boundaries.columns.get_loc('Province')]: # E.g. if this is an archdiocese (as provinces are named after archdioceses)\n",
    "                folium.Marker(center_point, icon = folium.DivIcon(icon_anchor=(horizontal_offset, vertical_offset), html=f\"{style_tag}<p><b><i>{diocese_name_for_plotting}</i><b></p>\")).add_to(name_feature_group)\n",
    "            else:\n",
    "                folium.Marker(center_point, icon = folium.DivIcon(icon_anchor=(horizontal_offset, vertical_offset), html=f\"{style_tag}<p><b>{diocese_name_for_plotting}</b></p>\")).add_to(name_feature_group)\n",
    "                # folium.Marker(center_point, icon = folium.DivIcon(icon_anchor=(horizontal_offset, 10), html=f\"{style_tag}<b><p>{diocese_name_for_plotting.replace(' ', '_').replace('-', '_')}</p></b>\")).add_to(revised_diocese_map)\n",
    "\n",
    "        name_feature_group.add_to(revised_diocese_map)\n",
    "\n",
    "    # The strategy of using an invisible folium.Marker object to display text on \n",
    "    # the map via DivIcon elements comes from from Bob Haffner's response at:\n",
    "    # https://stackoverflow.com/a/56878525/13097194 or:\n",
    "    # and https://stackoverflow.com/a/46408144/13097194\n",
    "\n",
    "    # for i in range(2):\n",
    "    #     loc_list = df_cathedral_locations.iloc[i, df_cathedral_locations.columns.get_loc('Cathedral_coord_list')].split('],')\n",
    "    #     loc_list = [loc.replace('[','').replace(']','') for loc in loc_list]\n",
    "    #     loc_list = [loc.split(', ') for loc in loc_list]\n",
    "    #     loc_list = [[float(loc[0]), float(loc[1])] for loc in loc_list]\n",
    "    #     for loc in loc_list:\n",
    "    #         folium.Marker(loc, icon = folium.Icon(icon = 'plane')).add_to(revised_diocese_map)\n",
    "\n",
    "\n",
    "            # folium.CircleMarker(loc).add_to(test_map)\n",
    "            \n",
    "        # Learn how to use Folium's FeatureGroup feature to toggle name and cathedral markers\n",
    "\n",
    "    folium.LatLngPopup().add_to(revised_diocese_map) # Useful for editing diocese\n",
    "    # boundaries\n",
    "\n",
    "    folium.LayerControl().add_to(revised_diocese_map)\n",
    "\n",
    "    revised_diocese_map.save(save_path+'.html')\n",
    "\n",
    "    # colormap.add_to(m)\n",
    "\n",
    "    return revised_diocese_map"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "create_diocese_map(tile_option = 'Stamen Toner')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "create_map_screenshot(absolute_path_to_map_folder = \n",
    "    r'C:\\Users\\kburc\\D1V1\\Documents\\!Dell64docs\\Programming\\py\\kjb3_programs\\us_diocese_mapper', \n",
    "    map_name = 'revised_diocese_map.html', screenshot_save_path = 'screenshots')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "def add_alaska_and_hawaii(starting_map_name, tile_option = 'Stamen Toner', \n",
    "include_cathedrals = True, include_names = True):\n",
    "    # This function creates a static version of the map defined by \n",
    "    # starting_map_name.html that includes Alaska and Hawaii at the bottom left.\n",
    "    # starting_map_name should not include the .html at the end. E.g. pass\n",
    "    # revised_diocese_map as an argument rather than revised_diocese_map.html.\n",
    "    create_diocese_map(tile_option = tile_option, starting_zoom = 4, starting_location = [65, -150], save_path = starting_map_name+'_alaska', include_cathedrals = include_cathedrals, include_names = include_names)\n",
    "    create_map_screenshot(absolute_path_to_map_folder = \n",
    "        r'C:\\Users\\kburc\\D1V1\\Documents\\!Dell64docs\\Programming\\py\\kjb3_programs\\us_diocese_mapper', \n",
    "        map_name = starting_map_name+'_alaska.html', screenshot_save_path = 'screenshots')\n",
    "    with PIL.Image.open('screenshots/'+starting_map_name+'_alaska.png') as map_image:\n",
    "        print(\"Opened\")\n",
    "        cropped_image = map_image.crop((2200, 1100, 3500, 2300)) # Left, upper, right,\n",
    "        # and lower points: https://pillow.readthedocs.io/en/stable/reference/Image.html#PIL.Image.Image.crop\n",
    "        # display(cropped_image) # https://stackoverflow.com/a/26649884/13097194\n",
    "        cropped_image.save('screenshots/'+starting_map_name+'_alaska_cropped.png')\n",
    "\n",
    "    create_diocese_map(tile_option = tile_option, starting_zoom = 6, starting_location = [21, -156], save_path = starting_map_name+'_hawaii', include_cathedrals = include_cathedrals, include_names = include_names)\n",
    "    create_map_screenshot(absolute_path_to_map_folder = \n",
    "        r'C:\\Users\\kburc\\D1V1\\Documents\\!Dell64docs\\Programming\\py\\kjb3_programs\\us_diocese_mapper', \n",
    "        map_name = starting_map_name+'_hawaii.html', screenshot_save_path = 'screenshots')\n",
    "    with PIL.Image.open('screenshots/'+starting_map_name+'_hawaii.png') as map_image:\n",
    "        print(\"Opened\")\n",
    "        cropped_image = map_image.crop((2500, 1400, 3300, 1900)) # Left, upper, right,\n",
    "        # # and lower points: https://pillow.readthedocs.io/en/stable/reference/Image.html#PIL.Image.Image.crop\n",
    "        cropped_image.save('screenshots/'+starting_map_name+'_hawaii_cropped.png')\n",
    "        \n",
    "    with PIL.Image.open('screenshots/'+starting_map_name+'.png') as map_image:\n",
    "        alaska = PIL.Image.open('screenshots/'+starting_map_name+'_alaska_cropped.png')\n",
    "        alaska = alaska.resize((int(alaska.width/1.25), int(alaska.height/1.25)))\n",
    "        hawaii = PIL.Image.open('screenshots/'+starting_map_name+'_hawaii_cropped.png')\n",
    "        map_image.paste(im = alaska, box = (0, 2300))\n",
    "        map_image.paste(im = hawaii, box = (1200, 2600))\n",
    "        map_image.save('screenshots/'+starting_map_name+'_50_states.png')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "add_alaska_and_hawaii('revised_diocese_map')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "create_diocese_map(tile_option = None, save_path = 'revised_diocese_map_tileless')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "add_alaska_and_hawaii(tile_option = None, starting_map_name = 'revised_diocese_map_tileless')"
   ]
  },
  {
   "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 = 'revised_diocese_map.html', screenshot_save_path = 'screenshots')\n",
    "\n",
    "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 = 'revised_diocese_map_tileless.html', screenshot_save_path = 'screenshots')\n",
    "\n",
    "convert_png_to_smaller_jpg(png_folder = 'screenshots', \n",
    "png_image_name = 'revised_diocese_map.png', jpg_folder = 'smaller_screenshots', \n",
    "quality_factor = 80)\n",
    "\n",
    "convert_png_to_smaller_jpg(png_folder = 'screenshots', \n",
    "png_image_name = 'revised_diocese_map_tileless.png', jpg_folder = 'smaller_screenshots', \n",
    "quality_factor = 80)\n",
    "\n",
    "IPython.display.Image(filename='smaller_screenshots/'+'revised_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": [
    "IPython.display.Image(filename='smaller_screenshots/'+'revised_diocese_map_tileless.jpg') \n",
    "# Based on https://ipython.readthedocs.io/en/stable/api/generated/IPython.display.html"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "create_diocese_map(tile_option = None, include_cathedrals = False, save_path = 'revised_diocese_map_no_cathedrals_tileless')\n",
    "create_map_screenshot(absolute_path_to_map_folder = \n",
    "    r'C:\\Users\\kburc\\D1V1\\Documents\\!Dell64docs\\Programming\\py\\kjb3_programs\\us_diocese_mapper', \n",
    "    map_name = 'revised_diocese_map_no_cathedrals_tileless.html', screenshot_save_path = 'screenshots')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "add_alaska_and_hawaii(tile_option = None, starting_map_name = 'revised_diocese_map_no_cathedrals_tileless', include_cathedrals = False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "create_diocese_map(tile_option = None, include_names = False, save_path = 'revised_diocese_map_no_names_tileless')\n",
    "create_map_screenshot(absolute_path_to_map_folder = \n",
    "    r'C:\\Users\\kburc\\D1V1\\Documents\\!Dell64docs\\Programming\\py\\kjb3_programs\\us_diocese_mapper', \n",
    "    map_name = 'revised_diocese_map_no_names_tileless.html', screenshot_save_path = 'screenshots')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "add_alaska_and_hawaii(tile_option = None, starting_map_name = 'revised_diocese_map_no_names_tileless', include_names = False)"
   ]
  },
  {
   "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"
   ]
  }
 ],
 "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
}