us_diocese_mapper_V6.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# US Diocese Mapper\n",
    "\n",
    "By Kenneth Burchfiel\n",
    "\n",
    "Released under the MIT License\n",
    "\n",
    "(The county_shapes.csv, county_shapes.geojson, state_fips_codes.csv, and counties_by_diocese.csv files are released into in the public domain.)\n",
    "\n",
    "## Note:\n",
    "This is still a work in progress, as I need to update the diocesan boundaries to account for situations where one county lies in multiple dioceses. (Until then, the map will contain a number of inaccuracies.) I also plan to add in cathedral locations and diocese names.\n",
    "\n",
    "I also plan to add in more documentation/comments."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Preliminary steps\n",
    "\n",
    "First, I downloaded county shapefiles from the US Census bureau here: https://www.census.gov/cgi-bin/geo/shapefiles/index.php?year=2021&layergroup=Counties+%28and+equivalent%29\n",
    "\n",
    "See this note regarding use of the shapefiles: https://www2.census.gov/geo/pdfs/maps-data/data/tiger/tgrshp2021/TGRSHP2021_TechDoc_Ch1.pdf\n",
    "\n",
    "Once you download the shapefiles to your computer, extract them using an unzipping utility. The shapefile document within this unzipped folder ends in .shp; for the 2021 county data, the file name is tl_2021_us_county.shp. \n",
    "\n",
    "It appears that the other files within this folder may also be used by Geopandas in the creation of GeoDataFrames, so I recommend accessing the file within the unzipped folder rather than copying it into your project folder.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import time\n",
    "start_time = time.time()\n",
    "import geopandas\n",
    "import pandas as pd\n",
    "import folium\n",
    "import branca\n",
    "import numpy as np\n",
    "from selenium import webdriver\n",
    "import PIL.Image\n",
    "import IPython"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "generate_county_list = True"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "if generate_county_list == True:\n",
    "\n",
    "    shapefile_path = \\\n",
    "        r'C:\\Users\\kburc\\Downloads\\tl_2021_us_county_download\\tl_2021_us_county.shp'\n",
    "    state_fips_codes = pd.read_html(\n",
    "        'https://www.nrcs.usda.gov/wps/portal/nrcs/detail/national/technical/nra/nri/results/?cid=nrcs143_013696')[0]\n",
    "    state_fips_codes = state_fips_codes.iloc[0:55].copy()\n",
    "    state_fips_codes['FIPS'] = state_fips_codes['FIPS'].astype('int')\n",
    "    state_fips_codes.rename(columns={'Name':'State_Name',\n",
    "    'Postal Code':'State_Code'},inplace=True) \n",
    "    state_fips_codes.loc[55] = ['District of Columbia', 'DC', 11] # DC wasn't \n",
    "    # present in the original webpage, so I added it in. I located its FIPS code\n",
    "    # within the county shapefile dataset.\n",
    "    state_fips_codes.to_csv('state_fips_codes.csv')\n",
    "\n",
    "    county_shapes = geopandas.read_file(shapefile_path)\n",
    "\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"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses = merged_county_shapes.copy()\n",
    "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": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>STATEFP</th>\n",
       "      <th>COUNTYFP</th>\n",
       "      <th>COUNTYNS</th>\n",
       "      <th>GEOID</th>\n",
       "      <th>NAME</th>\n",
       "      <th>NAMELSAD</th>\n",
       "      <th>LSAD</th>\n",
       "      <th>geometry</th>\n",
       "      <th>State_Name</th>\n",
       "      <th>State_Code</th>\n",
       "      <th>FIPS</th>\n",
       "      <th>county_state</th>\n",
       "      <th>Diocese</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>31</td>\n",
       "      <td>039</td>\n",
       "      <td>00835841</td>\n",
       "      <td>31039</td>\n",
       "      <td>Cuming</td>\n",
       "      <td>Cuming County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-96.55515 41.91587, -96.55515 41.914...</td>\n",
       "      <td>Nebraska</td>\n",
       "      <td>NE</td>\n",
       "      <td>31</td>\n",
       "      <td>Cuming County, Nebraska</td>\n",
       "      <td>Not_Set</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>53</td>\n",
       "      <td>069</td>\n",
       "      <td>01513275</td>\n",
       "      <td>53069</td>\n",
       "      <td>Wahkiakum</td>\n",
       "      <td>Wahkiakum County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-123.49077 46.38358, -123.48813 46.3...</td>\n",
       "      <td>Washington</td>\n",
       "      <td>WA</td>\n",
       "      <td>53</td>\n",
       "      <td>Wahkiakum County, Washington</td>\n",
       "      <td>Not_Set</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>35</td>\n",
       "      <td>011</td>\n",
       "      <td>00933054</td>\n",
       "      <td>35011</td>\n",
       "      <td>De Baca</td>\n",
       "      <td>De Baca County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-104.38368 34.69213, -104.37658 34.6...</td>\n",
       "      <td>New Mexico</td>\n",
       "      <td>NM</td>\n",
       "      <td>35</td>\n",
       "      <td>De Baca County, New Mexico</td>\n",
       "      <td>Not_Set</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>31</td>\n",
       "      <td>109</td>\n",
       "      <td>00835876</td>\n",
       "      <td>31109</td>\n",
       "      <td>Lancaster</td>\n",
       "      <td>Lancaster County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-96.68140 41.04566, -96.68139 41.045...</td>\n",
       "      <td>Nebraska</td>\n",
       "      <td>NE</td>\n",
       "      <td>31</td>\n",
       "      <td>Lancaster County, Nebraska</td>\n",
       "      <td>Not_Set</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>31</td>\n",
       "      <td>129</td>\n",
       "      <td>00835886</td>\n",
       "      <td>31129</td>\n",
       "      <td>Nuckolls</td>\n",
       "      <td>Nuckolls County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-98.04802 40.35066, -98.04674 40.350...</td>\n",
       "      <td>Nebraska</td>\n",
       "      <td>NE</td>\n",
       "      <td>31</td>\n",
       "      <td>Nuckolls County, Nebraska</td>\n",
       "      <td>Not_Set</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3229</th>\n",
       "      <td>13</td>\n",
       "      <td>123</td>\n",
       "      <td>00351260</td>\n",
       "      <td>13123</td>\n",
       "      <td>Gilmer</td>\n",
       "      <td>Gilmer County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-84.30237 34.57832, -84.30329 34.577...</td>\n",
       "      <td>Georgia</td>\n",
       "      <td>GA</td>\n",
       "      <td>13</td>\n",
       "      <td>Gilmer County, Georgia</td>\n",
       "      <td>Not_Set</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3230</th>\n",
       "      <td>27</td>\n",
       "      <td>135</td>\n",
       "      <td>00659513</td>\n",
       "      <td>27135</td>\n",
       "      <td>Roseau</td>\n",
       "      <td>Roseau County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-95.25857 48.88666, -95.25707 48.885...</td>\n",
       "      <td>Minnesota</td>\n",
       "      <td>MN</td>\n",
       "      <td>27</td>\n",
       "      <td>Roseau County, Minnesota</td>\n",
       "      <td>Not_Set</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3231</th>\n",
       "      <td>28</td>\n",
       "      <td>089</td>\n",
       "      <td>00695768</td>\n",
       "      <td>28089</td>\n",
       "      <td>Madison</td>\n",
       "      <td>Madison County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-90.14883 32.40026, -90.14890 32.400...</td>\n",
       "      <td>Mississippi</td>\n",
       "      <td>MS</td>\n",
       "      <td>28</td>\n",
       "      <td>Madison County, Mississippi</td>\n",
       "      <td>Not_Set</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3232</th>\n",
       "      <td>48</td>\n",
       "      <td>227</td>\n",
       "      <td>01383899</td>\n",
       "      <td>48227</td>\n",
       "      <td>Howard</td>\n",
       "      <td>Howard County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-101.18138 32.21252, -101.18138 32.2...</td>\n",
       "      <td>Texas</td>\n",
       "      <td>TX</td>\n",
       "      <td>48</td>\n",
       "      <td>Howard County, Texas</td>\n",
       "      <td>Not_Set</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3233</th>\n",
       "      <td>54</td>\n",
       "      <td>099</td>\n",
       "      <td>01550056</td>\n",
       "      <td>54099</td>\n",
       "      <td>Wayne</td>\n",
       "      <td>Wayne County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-82.30872 38.28106, -82.30874 38.280...</td>\n",
       "      <td>West Virginia</td>\n",
       "      <td>WV</td>\n",
       "      <td>54</td>\n",
       "      <td>Wayne County, West Virginia</td>\n",
       "      <td>Not_Set</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3146 rows × 13 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      STATEFP COUNTYFP  COUNTYNS  GEOID       NAME          NAMELSAD LSAD  \\\n",
       "0          31      039  00835841  31039     Cuming     Cuming County   06   \n",
       "1          53      069  01513275  53069  Wahkiakum  Wahkiakum County   06   \n",
       "2          35      011  00933054  35011    De Baca    De Baca County   06   \n",
       "3          31      109  00835876  31109  Lancaster  Lancaster County   06   \n",
       "4          31      129  00835886  31129   Nuckolls   Nuckolls County   06   \n",
       "...       ...      ...       ...    ...        ...               ...  ...   \n",
       "3229       13      123  00351260  13123     Gilmer     Gilmer County   06   \n",
       "3230       27      135  00659513  27135     Roseau     Roseau County   06   \n",
       "3231       28      089  00695768  28089    Madison    Madison County   06   \n",
       "3232       48      227  01383899  48227     Howard     Howard County   06   \n",
       "3233       54      099  01550056  54099      Wayne      Wayne County   06   \n",
       "\n",
       "                                               geometry     State_Name  \\\n",
       "0     POLYGON ((-96.55515 41.91587, -96.55515 41.914...       Nebraska   \n",
       "1     POLYGON ((-123.49077 46.38358, -123.48813 46.3...     Washington   \n",
       "2     POLYGON ((-104.38368 34.69213, -104.37658 34.6...     New Mexico   \n",
       "3     POLYGON ((-96.68140 41.04566, -96.68139 41.045...       Nebraska   \n",
       "4     POLYGON ((-98.04802 40.35066, -98.04674 40.350...       Nebraska   \n",
       "...                                                 ...            ...   \n",
       "3229  POLYGON ((-84.30237 34.57832, -84.30329 34.577...        Georgia   \n",
       "3230  POLYGON ((-95.25857 48.88666, -95.25707 48.885...      Minnesota   \n",
       "3231  POLYGON ((-90.14883 32.40026, -90.14890 32.400...    Mississippi   \n",
       "3232  POLYGON ((-101.18138 32.21252, -101.18138 32.2...          Texas   \n",
       "3233  POLYGON ((-82.30872 38.28106, -82.30874 38.280...  West Virginia   \n",
       "\n",
       "     State_Code  FIPS                  county_state  Diocese  \n",
       "0            NE    31       Cuming County, Nebraska  Not_Set  \n",
       "1            WA    53  Wahkiakum County, Washington  Not_Set  \n",
       "2            NM    35    De Baca County, New Mexico  Not_Set  \n",
       "3            NE    31    Lancaster County, Nebraska  Not_Set  \n",
       "4            NE    31     Nuckolls County, Nebraska  Not_Set  \n",
       "...         ...   ...                           ...      ...  \n",
       "3229         GA    13        Gilmer County, Georgia  Not_Set  \n",
       "3230         MN    27      Roseau County, Minnesota  Not_Set  \n",
       "3231         MS    28   Madison County, Mississippi  Not_Set  \n",
       "3232         TX    48          Howard County, Texas  Not_Set  \n",
       "3233         WV    54   Wayne County, West Virginia  Not_Set  \n",
       "\n",
       "[3146 rows x 13 columns]"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_dioceses['Diocese'] = 'Not_Set'\n",
    "df_dioceses"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "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": 7,
   "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', # 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',\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",
    "'Saint Thomas'\n",
    "]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses['Diocese'] = np.select(condlist, choicelist, df_dioceses['Diocese'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "3146"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(df_dioceses)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>STATEFP</th>\n",
       "      <th>COUNTYFP</th>\n",
       "      <th>COUNTYNS</th>\n",
       "      <th>GEOID</th>\n",
       "      <th>NAME</th>\n",
       "      <th>NAMELSAD</th>\n",
       "      <th>LSAD</th>\n",
       "      <th>geometry</th>\n",
       "      <th>State_Name</th>\n",
       "      <th>State_Code</th>\n",
       "      <th>FIPS</th>\n",
       "      <th>county_state</th>\n",
       "      <th>Diocese</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>31</td>\n",
       "      <td>039</td>\n",
       "      <td>00835841</td>\n",
       "      <td>31039</td>\n",
       "      <td>Cuming</td>\n",
       "      <td>Cuming County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-96.55515 41.91587, -96.55515 41.914...</td>\n",
       "      <td>Nebraska</td>\n",
       "      <td>NE</td>\n",
       "      <td>31</td>\n",
       "      <td>Cuming County, Nebraska</td>\n",
       "      <td>Grand Island</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>53</td>\n",
       "      <td>069</td>\n",
       "      <td>01513275</td>\n",
       "      <td>53069</td>\n",
       "      <td>Wahkiakum</td>\n",
       "      <td>Wahkiakum County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-123.49077 46.38358, -123.48813 46.3...</td>\n",
       "      <td>Washington</td>\n",
       "      <td>WA</td>\n",
       "      <td>53</td>\n",
       "      <td>Wahkiakum County, Washington</td>\n",
       "      <td>Seattle</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>35</td>\n",
       "      <td>011</td>\n",
       "      <td>00933054</td>\n",
       "      <td>35011</td>\n",
       "      <td>De Baca</td>\n",
       "      <td>De Baca County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-104.38368 34.69213, -104.37658 34.6...</td>\n",
       "      <td>New Mexico</td>\n",
       "      <td>NM</td>\n",
       "      <td>35</td>\n",
       "      <td>De Baca County, New Mexico</td>\n",
       "      <td>Santa Fe</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>31</td>\n",
       "      <td>109</td>\n",
       "      <td>00835876</td>\n",
       "      <td>31109</td>\n",
       "      <td>Lancaster</td>\n",
       "      <td>Lancaster County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-96.68140 41.04566, -96.68139 41.045...</td>\n",
       "      <td>Nebraska</td>\n",
       "      <td>NE</td>\n",
       "      <td>31</td>\n",
       "      <td>Lancaster County, Nebraska</td>\n",
       "      <td>Grand Island</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>31</td>\n",
       "      <td>129</td>\n",
       "      <td>00835886</td>\n",
       "      <td>31129</td>\n",
       "      <td>Nuckolls</td>\n",
       "      <td>Nuckolls County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-98.04802 40.35066, -98.04674 40.350...</td>\n",
       "      <td>Nebraska</td>\n",
       "      <td>NE</td>\n",
       "      <td>31</td>\n",
       "      <td>Nuckolls County, Nebraska</td>\n",
       "      <td>Grand Island</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3229</th>\n",
       "      <td>13</td>\n",
       "      <td>123</td>\n",
       "      <td>00351260</td>\n",
       "      <td>13123</td>\n",
       "      <td>Gilmer</td>\n",
       "      <td>Gilmer County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-84.30237 34.57832, -84.30329 34.577...</td>\n",
       "      <td>Georgia</td>\n",
       "      <td>GA</td>\n",
       "      <td>13</td>\n",
       "      <td>Gilmer County, Georgia</td>\n",
       "      <td>Savannah</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3230</th>\n",
       "      <td>27</td>\n",
       "      <td>135</td>\n",
       "      <td>00659513</td>\n",
       "      <td>27135</td>\n",
       "      <td>Roseau</td>\n",
       "      <td>Roseau County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-95.25857 48.88666, -95.25707 48.885...</td>\n",
       "      <td>Minnesota</td>\n",
       "      <td>MN</td>\n",
       "      <td>27</td>\n",
       "      <td>Roseau County, Minnesota</td>\n",
       "      <td>Duluth</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3231</th>\n",
       "      <td>28</td>\n",
       "      <td>089</td>\n",
       "      <td>00695768</td>\n",
       "      <td>28089</td>\n",
       "      <td>Madison</td>\n",
       "      <td>Madison County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-90.14883 32.40026, -90.14890 32.400...</td>\n",
       "      <td>Mississippi</td>\n",
       "      <td>MS</td>\n",
       "      <td>28</td>\n",
       "      <td>Madison County, Mississippi</td>\n",
       "      <td>Jackson</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3232</th>\n",
       "      <td>48</td>\n",
       "      <td>227</td>\n",
       "      <td>01383899</td>\n",
       "      <td>48227</td>\n",
       "      <td>Howard</td>\n",
       "      <td>Howard County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-101.18138 32.21252, -101.18138 32.2...</td>\n",
       "      <td>Texas</td>\n",
       "      <td>TX</td>\n",
       "      <td>48</td>\n",
       "      <td>Howard County, Texas</td>\n",
       "      <td>Lubbock</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3233</th>\n",
       "      <td>54</td>\n",
       "      <td>099</td>\n",
       "      <td>01550056</td>\n",
       "      <td>54099</td>\n",
       "      <td>Wayne</td>\n",
       "      <td>Wayne County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-82.30872 38.28106, -82.30874 38.280...</td>\n",
       "      <td>West Virginia</td>\n",
       "      <td>WV</td>\n",
       "      <td>54</td>\n",
       "      <td>Wayne County, West Virginia</td>\n",
       "      <td>Wheeling-Charleston</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3146 rows × 13 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      STATEFP COUNTYFP  COUNTYNS  GEOID       NAME          NAMELSAD LSAD  \\\n",
       "0          31      039  00835841  31039     Cuming     Cuming County   06   \n",
       "1          53      069  01513275  53069  Wahkiakum  Wahkiakum County   06   \n",
       "2          35      011  00933054  35011    De Baca    De Baca County   06   \n",
       "3          31      109  00835876  31109  Lancaster  Lancaster County   06   \n",
       "4          31      129  00835886  31129   Nuckolls   Nuckolls County   06   \n",
       "...       ...      ...       ...    ...        ...               ...  ...   \n",
       "3229       13      123  00351260  13123     Gilmer     Gilmer County   06   \n",
       "3230       27      135  00659513  27135     Roseau     Roseau County   06   \n",
       "3231       28      089  00695768  28089    Madison    Madison County   06   \n",
       "3232       48      227  01383899  48227     Howard     Howard County   06   \n",
       "3233       54      099  01550056  54099      Wayne      Wayne County   06   \n",
       "\n",
       "                                               geometry     State_Name  \\\n",
       "0     POLYGON ((-96.55515 41.91587, -96.55515 41.914...       Nebraska   \n",
       "1     POLYGON ((-123.49077 46.38358, -123.48813 46.3...     Washington   \n",
       "2     POLYGON ((-104.38368 34.69213, -104.37658 34.6...     New Mexico   \n",
       "3     POLYGON ((-96.68140 41.04566, -96.68139 41.045...       Nebraska   \n",
       "4     POLYGON ((-98.04802 40.35066, -98.04674 40.350...       Nebraska   \n",
       "...                                                 ...            ...   \n",
       "3229  POLYGON ((-84.30237 34.57832, -84.30329 34.577...        Georgia   \n",
       "3230  POLYGON ((-95.25857 48.88666, -95.25707 48.885...      Minnesota   \n",
       "3231  POLYGON ((-90.14883 32.40026, -90.14890 32.400...    Mississippi   \n",
       "3232  POLYGON ((-101.18138 32.21252, -101.18138 32.2...          Texas   \n",
       "3233  POLYGON ((-82.30872 38.28106, -82.30874 38.280...  West Virginia   \n",
       "\n",
       "     State_Code  FIPS                  county_state              Diocese  \n",
       "0            NE    31       Cuming County, Nebraska         Grand Island  \n",
       "1            WA    53  Wahkiakum County, Washington              Seattle  \n",
       "2            NM    35    De Baca County, New Mexico             Santa Fe  \n",
       "3            NE    31    Lancaster County, Nebraska         Grand Island  \n",
       "4            NE    31     Nuckolls County, Nebraska         Grand Island  \n",
       "...         ...   ...                           ...                  ...  \n",
       "3229         GA    13        Gilmer County, Georgia             Savannah  \n",
       "3230         MN    27      Roseau County, Minnesota               Duluth  \n",
       "3231         MS    28   Madison County, Mississippi              Jackson  \n",
       "3232         TX    48          Howard County, Texas              Lubbock  \n",
       "3233         WV    54   Wayne County, West Virginia  Wheeling-Charleston  \n",
       "\n",
       "[3146 rows x 13 columns]"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_dioceses"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "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": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "def add_dioceses(df, county_list, state, diocese):\n",
    "    # Adds the diocese for all counties in a given state. If a diocese spans\n",
    "    # multiple stats, this function will need to be called once for each\n",
    "    # state.\n",
    "    # If a given 'county' variable does not have any spaces, \n",
    "    # the loop below will add ' County' to that string before \n",
    "    # adding in the state. Otherwise, the code assumes that the 'county' is\n",
    "    # either a non-county (e.g. a municipality, city, census area, parish,\n",
    "    # or borough), and will only add the state to this string.\n",
    "    # Therefore, in order to run this function successfully, enter the \n",
    "    # area designation (county, city, etc.) after all non-counties and\n",
    "    # counties with more than one word when building the county list.\n",
    "    # However, you don't need to add in 'County' for counties with only\n",
    "    # one word.\n",
    "    for county in county_list:\n",
    "        if ' ' not in county:\n",
    "            county_string = county + ' County, ' + state\n",
    "        else:\n",
    "            county_string = county + ', ' + state\n",
    "        if county_string not in df.index:\n",
    "            raise ValueError(county_string, \"not found in index\")\n",
    "        df.at[county_string, 'Diocese'] = diocese\n",
    "    return df\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "List of dioceses:\n",
    "https://en.wikipedia.org/wiki/List_of_Catholic_dioceses_in_the_United_States\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_dioceses.set_index('county_state', drop = False, inplace = True)\n",
    "\n",
    "# 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 in California:\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 in California')\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 in California:\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 in California' \n",
    "\n",
    "# Diocese of Santa Rosa in California:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Santa_Rosa_in_California\n",
    "\n",
    "santa_rosa_in_california_county_list = [\n",
    "'Del Norte County', 'Humboldt', 'Lake', 'Mendocino', 'Napa', 'Sonoma'\n",
    "]\n",
    "\n",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "santa_rosa_in_california_county_list, state = 'California', \n",
    "diocese = 'Santa Rosa in California')\n",
    "\n",
    "\n",
    "# Diocese of Stockton:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Stockton\n",
    "\n",
    "stockton_county_list = [\n",
    "'Alpine', 'Calaveras', 'Mono', 'San Joaquin County', 'Stanislaus', 'Tuolumne'\n",
    "]\n",
    "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 in Florida:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Venice_in_Florida\n",
    "\n",
    "venice_in_florida_county_list = [\n",
    "'Charlotte', 'Collier', 'DeSoto', 'Glades', 'Hardee', 'Hendry', 'Highlands', \n",
    "'Lee', 'Manatee', 'Sarasota'\n",
    "]\n",
    "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 Saint Petersburg:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Saint_Petersburg\n",
    "saint_petersburg_county_list = [\n",
    "'Pinellas', 'Hillsborough', 'Pasco', 'Hernando', 'Citrus'\n",
    "]\n",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "saint_petersburg_county_list, state = 'Florida', \n",
    "diocese = 'Saint Petersburg')\n",
    "\n",
    "# Diocese of St. Augustine:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_St._Augustine\n",
    "st_augustine_county_list = [\n",
    "'Alachua', 'Baker', 'Bradford', 'Clay', 'Columbia', 'Dixie', 'Duval', 'Flagler',\n",
    "'Gilchrist', 'Hamilton', 'Lafayette', 'Levy', 'Nassau', 'Putnam', \n",
    "'St. Johns County', 'Suwannee', 'Union'\n",
    "]\n",
    "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 in Louisiana:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Alexandria_in_Louisiana\n",
    "\n",
    "alexandria_in_louisiana_county_list = [\n",
    "'Avoyelles Parish', 'Rapides Parish', 'Vernon Parish', 'Natchitoches Parish',\n",
    "'Winn Parish', 'Caldwell Parish', 'Madison Parish', 'Franklin Parish', \n",
    "'Tensas Parish', 'Concordia Parish', 'Catahoula Parish', 'LaSalle Parish',\n",
    "'Grant Parish'\n",
    "]\n",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "alexandria_in_louisiana_county_list, state = 'Louisiana', \n",
    "diocese = 'Alexandria in Louisiana')\n",
    "\n",
    "# Archdiocese of New York:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_New_York\n",
    "new_york_county_list = ['New York County', 'Richmond', 'Bronx', 'Dutchess',\n",
    "'Orange', 'Putnam', 'Rockland', 'Sullivan', 'Ulster', 'Westchester'\n",
    "]\n",
    "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 Saint Paul and Minneapolis:\n",
    "# Source: https://commons.wikimedia.org/wiki/File:Map_of_the_Catholic_archdiocese_of_Saint_Paul_%26_Minneapolis.svg\n",
    "# See also: https://s3.amazonaws.com/archspmmainsite/Resources/ArchMap--Parishes-Deaneries-COLORRamped_061516.pdf\n",
    "\n",
    "saint_paul_and_minneapolis_county_list = [\n",
    "'Goodhue', 'Rice', 'Le Sueur County', 'Scott', 'Dakota', 'Washington', 'Ramsey',\n",
    "'Hennepin', 'Wright', 'Anoka', 'Chisago', 'Carver'\n",
    "]\n",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "saint_paul_and_minneapolis_county_list, state = 'Minnesota', \n",
    "diocese = 'Saint Paul and Minneapolis')\n",
    "\n",
    "# Diocese of Crookston:\n",
    "# https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Crookston\n",
    "\n",
    "crookston_county_list = [\n",
    "'Kittson', 'Roseau', 'Lake of the Woods County', 'Marshall', 'Polk', \n",
    "'Red Lake County', 'Pennington', 'Clearwater', 'Beltrami', 'Norman', 'Mahnomen', \n",
    "'Hubbard', 'Clay', 'Becker'\n",
    "]\n",
    "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 Saint Cloud:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Saint_Cloud\n",
    "\n",
    "saint_cloud_county_list = [\n",
    "'Benton', 'Douglas', 'Grant', 'Isanti', 'Kanabec', 'Mille Lacs County', \n",
    "'Morrison', 'Otter Tail County', 'Pope', 'Sherburne', 'Stearns', 'Stevens',\n",
    "'Todd', 'Traverse', 'Wadena', 'Wilkin'\n",
    "]\n",
    "df_dioceses = add_dioceses(df_dioceses, \n",
    "saint_cloud_county_list, state = 'Minnesota', \n",
    "diocese = 'Saint Cloud')\n",
    "\n",
    "\n",
    "# Diocese of Winona-Rochester:\n",
    "# Source: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Winona%E2%80%93Rochester\n",
    "\n",
    "winona_rochester_county_list = [\n",
    "'Blue Earth County', 'Cottonwood', 'Dodge', 'Faribault', 'Fillmore', 'Freeborn',\n",
    "'Houston', 'Jackson', 'Martin', 'Mower', 'Murray', 'Nobles', 'Olmsted',\n",
    "'Pipestone', 'Rock', 'Steele', 'Wabasha', 'Waseca', 'Watonwan', 'Winona'\n",
    "]\n",
    "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",
    "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": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1765"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "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": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>STATEFP</th>\n",
       "      <th>COUNTYFP</th>\n",
       "      <th>COUNTYNS</th>\n",
       "      <th>GEOID</th>\n",
       "      <th>NAME</th>\n",
       "      <th>NAMELSAD</th>\n",
       "      <th>LSAD</th>\n",
       "      <th>geometry</th>\n",
       "      <th>State_Name</th>\n",
       "      <th>State_Code</th>\n",
       "      <th>FIPS</th>\n",
       "      <th>county_state</th>\n",
       "      <th>Diocese</th>\n",
       "      <th>Province</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>48</td>\n",
       "      <td>011</td>\n",
       "      <td>01383791</td>\n",
       "      <td>48011</td>\n",
       "      <td>Armstrong</td>\n",
       "      <td>Armstrong County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-101.08728 35.02708, -101.08729 35.0...</td>\n",
       "      <td>Texas</td>\n",
       "      <td>TX</td>\n",
       "      <td>48</td>\n",
       "      <td>Armstrong County, Texas</td>\n",
       "      <td>Amarillo</td>\n",
       "      <td>San Antonio</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>37</th>\n",
       "      <td>48</td>\n",
       "      <td>077</td>\n",
       "      <td>01383824</td>\n",
       "      <td>48077</td>\n",
       "      <td>Clay</td>\n",
       "      <td>Clay County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-97.97784 33.81417, -97.97784 33.813...</td>\n",
       "      <td>Texas</td>\n",
       "      <td>TX</td>\n",
       "      <td>48</td>\n",
       "      <td>Clay County, Texas</td>\n",
       "      <td>Fort Worth</td>\n",
       "      <td>San Antonio</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>48</td>\n",
       "      <td>189</td>\n",
       "      <td>01383880</td>\n",
       "      <td>48189</td>\n",
       "      <td>Hale</td>\n",
       "      <td>Hale County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-101.73504 33.82864, -101.73507 33.8...</td>\n",
       "      <td>Texas</td>\n",
       "      <td>TX</td>\n",
       "      <td>48</td>\n",
       "      <td>Hale County, Texas</td>\n",
       "      <td>Lubbock</td>\n",
       "      <td>San Antonio</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>48</td>\n",
       "      <td>327</td>\n",
       "      <td>01383949</td>\n",
       "      <td>48327</td>\n",
       "      <td>Menard</td>\n",
       "      <td>Menard County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-99.77120 30.71079, -99.78240 30.710...</td>\n",
       "      <td>Texas</td>\n",
       "      <td>TX</td>\n",
       "      <td>48</td>\n",
       "      <td>Menard County, Texas</td>\n",
       "      <td>San Angelo</td>\n",
       "      <td>San Antonio</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1206</th>\n",
       "      <td>48</td>\n",
       "      <td>377</td>\n",
       "      <td>01383974</td>\n",
       "      <td>48377</td>\n",
       "      <td>Presidio</td>\n",
       "      <td>Presidio County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-103.79577 29.72453, -103.79577 29.7...</td>\n",
       "      <td>Texas</td>\n",
       "      <td>TX</td>\n",
       "      <td>48</td>\n",
       "      <td>Presidio County, Texas</td>\n",
       "      <td>El Paso</td>\n",
       "      <td>San Antonio</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>120</th>\n",
       "      <td>48</td>\n",
       "      <td>247</td>\n",
       "      <td>01383909</td>\n",
       "      <td>48247</td>\n",
       "      <td>Jim Hogg</td>\n",
       "      <td>Jim Hogg County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-98.49291 27.17469, -98.49292 27.174...</td>\n",
       "      <td>Texas</td>\n",
       "      <td>TX</td>\n",
       "      <td>48</td>\n",
       "      <td>Jim Hogg County, Texas</td>\n",
       "      <td>Laredo</td>\n",
       "      <td>San Antonio</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>359</th>\n",
       "      <td>15</td>\n",
       "      <td>003</td>\n",
       "      <td>00365281</td>\n",
       "      <td>15003</td>\n",
       "      <td>Honolulu</td>\n",
       "      <td>Honolulu County</td>\n",
       "      <td>06</td>\n",
       "      <td>MULTIPOLYGON (((-166.38653 23.84636, -166.3842...</td>\n",
       "      <td>Hawaii</td>\n",
       "      <td>HI</td>\n",
       "      <td>15</td>\n",
       "      <td>Honolulu County, Hawaii</td>\n",
       "      <td>Honolulu</td>\n",
       "      <td>San Francisco</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>618</th>\n",
       "      <td>6</td>\n",
       "      <td>075</td>\n",
       "      <td>00277302</td>\n",
       "      <td>06075</td>\n",
       "      <td>San Francisco</td>\n",
       "      <td>San Francisco County</td>\n",
       "      <td>06</td>\n",
       "      <td>MULTIPOLYGON (((-122.40537 37.70826, -122.4054...</td>\n",
       "      <td>California</td>\n",
       "      <td>CA</td>\n",
       "      <td>6</td>\n",
       "      <td>San Francisco County, California</td>\n",
       "      <td>San Francisco</td>\n",
       "      <td>San Francisco</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1885</th>\n",
       "      <td>6</td>\n",
       "      <td>085</td>\n",
       "      <td>00277307</td>\n",
       "      <td>06085</td>\n",
       "      <td>Santa Clara</td>\n",
       "      <td>Santa Clara County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-121.47115 36.97797, -121.47121 36.9...</td>\n",
       "      <td>California</td>\n",
       "      <td>CA</td>\n",
       "      <td>6</td>\n",
       "      <td>Santa Clara County, California</td>\n",
       "      <td>San Jose in California</td>\n",
       "      <td>San Francisco</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>340</th>\n",
       "      <td>6</td>\n",
       "      <td>009</td>\n",
       "      <td>01675885</td>\n",
       "      <td>06009</td>\n",
       "      <td>Calaveras</td>\n",
       "      <td>Calaveras County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-120.37620 38.14265, -120.37616 38.1...</td>\n",
       "      <td>California</td>\n",
       "      <td>CA</td>\n",
       "      <td>6</td>\n",
       "      <td>Calaveras County, California</td>\n",
       "      <td>Stockton</td>\n",
       "      <td>San Francisco</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>149</th>\n",
       "      <td>32</td>\n",
       "      <td>003</td>\n",
       "      <td>00863600</td>\n",
       "      <td>32003</td>\n",
       "      <td>Clark</td>\n",
       "      <td>Clark County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-115.89692 36.84208, -115.89159 36.8...</td>\n",
       "      <td>Nevada</td>\n",
       "      <td>NV</td>\n",
       "      <td>32</td>\n",
       "      <td>Clark County, Nevada</td>\n",
       "      <td>Las Vegas</td>\n",
       "      <td>San Francisco</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>444</th>\n",
       "      <td>6</td>\n",
       "      <td>097</td>\n",
       "      <td>01657246</td>\n",
       "      <td>06097</td>\n",
       "      <td>Sonoma</td>\n",
       "      <td>Sonoma County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-122.63498 38.58301, -122.63496 38.5...</td>\n",
       "      <td>California</td>\n",
       "      <td>CA</td>\n",
       "      <td>6</td>\n",
       "      <td>Sonoma County, California</td>\n",
       "      <td>Santa Rosa in California</td>\n",
       "      <td>San Francisco</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>49</td>\n",
       "      <td>033</td>\n",
       "      <td>01448030</td>\n",
       "      <td>49033</td>\n",
       "      <td>Rich</td>\n",
       "      <td>Rich County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-111.04653 41.77964, -111.04651 41.7...</td>\n",
       "      <td>Utah</td>\n",
       "      <td>UT</td>\n",
       "      <td>49</td>\n",
       "      <td>Rich County, Utah</td>\n",
       "      <td>Salt Lake City</td>\n",
       "      <td>San Francisco</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1098</th>\n",
       "      <td>6</td>\n",
       "      <td>013</td>\n",
       "      <td>01675903</td>\n",
       "      <td>06013</td>\n",
       "      <td>Contra Costa</td>\n",
       "      <td>Contra Costa County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-121.91046 37.73039, -121.91057 37.7...</td>\n",
       "      <td>California</td>\n",
       "      <td>CA</td>\n",
       "      <td>6</td>\n",
       "      <td>Contra Costa County, California</td>\n",
       "      <td>Oakland</td>\n",
       "      <td>San Francisco</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>6</td>\n",
       "      <td>091</td>\n",
       "      <td>00277310</td>\n",
       "      <td>06091</td>\n",
       "      <td>Sierra</td>\n",
       "      <td>Sierra County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-120.55587 39.50874, -120.55614 39.5...</td>\n",
       "      <td>California</td>\n",
       "      <td>CA</td>\n",
       "      <td>6</td>\n",
       "      <td>Sierra County, California</td>\n",
       "      <td>Sacramento</td>\n",
       "      <td>San Francisco</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>81</th>\n",
       "      <td>32</td>\n",
       "      <td>001</td>\n",
       "      <td>00858638</td>\n",
       "      <td>32001</td>\n",
       "      <td>Churchill</td>\n",
       "      <td>Churchill County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-118.18700 40.00209, -118.17322 40.0...</td>\n",
       "      <td>Nevada</td>\n",
       "      <td>NV</td>\n",
       "      <td>32</td>\n",
       "      <td>Churchill County, Nevada</td>\n",
       "      <td>Reno</td>\n",
       "      <td>San Francisco</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>35</td>\n",
       "      <td>011</td>\n",
       "      <td>00933054</td>\n",
       "      <td>35011</td>\n",
       "      <td>De Baca</td>\n",
       "      <td>De Baca County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-104.38368 34.69213, -104.37658 34.6...</td>\n",
       "      <td>New Mexico</td>\n",
       "      <td>NM</td>\n",
       "      <td>35</td>\n",
       "      <td>De Baca County, New Mexico</td>\n",
       "      <td>Santa Fe</td>\n",
       "      <td>Santa Fe</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>35</td>\n",
       "      <td>035</td>\n",
       "      <td>00929104</td>\n",
       "      <td>35035</td>\n",
       "      <td>Otero</td>\n",
       "      <td>Otero County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-105.93164 33.38990, -105.93089 33.3...</td>\n",
       "      <td>New Mexico</td>\n",
       "      <td>NM</td>\n",
       "      <td>35</td>\n",
       "      <td>Otero County, New Mexico</td>\n",
       "      <td>Las Cruces</td>\n",
       "      <td>Santa Fe</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>391</th>\n",
       "      <td>4</td>\n",
       "      <td>027</td>\n",
       "      <td>00023901</td>\n",
       "      <td>04027</td>\n",
       "      <td>Yuma</td>\n",
       "      <td>Yuma County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-114.47325 33.02788, -114.45989 33.0...</td>\n",
       "      <td>Arizona</td>\n",
       "      <td>AZ</td>\n",
       "      <td>4</td>\n",
       "      <td>Yuma County, Arizona</td>\n",
       "      <td>Tucson</td>\n",
       "      <td>Santa Fe</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>35</td>\n",
       "      <td>003</td>\n",
       "      <td>00929108</td>\n",
       "      <td>35003</td>\n",
       "      <td>Catron</td>\n",
       "      <td>Catron County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-108.20839 33.20097, -108.21384 33.2...</td>\n",
       "      <td>New Mexico</td>\n",
       "      <td>NM</td>\n",
       "      <td>35</td>\n",
       "      <td>Catron County, New Mexico</td>\n",
       "      <td>Gallup</td>\n",
       "      <td>Santa Fe</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>832</th>\n",
       "      <td>4</td>\n",
       "      <td>013</td>\n",
       "      <td>00037026</td>\n",
       "      <td>04013</td>\n",
       "      <td>Maricopa</td>\n",
       "      <td>Maricopa County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-111.89308 33.20471, -111.90205 33.2...</td>\n",
       "      <td>Arizona</td>\n",
       "      <td>AZ</td>\n",
       "      <td>4</td>\n",
       "      <td>Maricopa County, Arizona</td>\n",
       "      <td>Phoenix</td>\n",
       "      <td>Santa Fe</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1129</th>\n",
       "      <td>53</td>\n",
       "      <td>013</td>\n",
       "      <td>01513273</td>\n",
       "      <td>53013</td>\n",
       "      <td>Columbia</td>\n",
       "      <td>Columbia County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-117.73737 46.44933, -117.73742 46.4...</td>\n",
       "      <td>Washington</td>\n",
       "      <td>WA</td>\n",
       "      <td>53</td>\n",
       "      <td>Columbia County, Washington</td>\n",
       "      <td>Spokane</td>\n",
       "      <td>Seattle</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1168</th>\n",
       "      <td>53</td>\n",
       "      <td>037</td>\n",
       "      <td>01531926</td>\n",
       "      <td>53037</td>\n",
       "      <td>Kittitas</td>\n",
       "      <td>Kittitas County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-120.00214 46.99732, -120.00193 46.9...</td>\n",
       "      <td>Washington</td>\n",
       "      <td>WA</td>\n",
       "      <td>53</td>\n",
       "      <td>Kittitas County, Washington</td>\n",
       "      <td>Yakima</td>\n",
       "      <td>Seattle</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>53</td>\n",
       "      <td>069</td>\n",
       "      <td>01513275</td>\n",
       "      <td>53069</td>\n",
       "      <td>Wahkiakum</td>\n",
       "      <td>Wahkiakum County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-123.49077 46.38358, -123.48813 46.3...</td>\n",
       "      <td>Washington</td>\n",
       "      <td>WA</td>\n",
       "      <td>53</td>\n",
       "      <td>Wahkiakum County, Washington</td>\n",
       "      <td>Seattle</td>\n",
       "      <td>Seattle</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>727</th>\n",
       "      <td>29</td>\n",
       "      <td>221</td>\n",
       "      <td>00758563</td>\n",
       "      <td>29221</td>\n",
       "      <td>Washington</td>\n",
       "      <td>Washington County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-90.81877 37.73553, -90.82091 37.735...</td>\n",
       "      <td>Missouri</td>\n",
       "      <td>MO</td>\n",
       "      <td>29</td>\n",
       "      <td>Washington County, Missouri</td>\n",
       "      <td>St. Louis</td>\n",
       "      <td>St. Louis</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33</th>\n",
       "      <td>29</td>\n",
       "      <td>011</td>\n",
       "      <td>00758460</td>\n",
       "      <td>29011</td>\n",
       "      <td>Barton</td>\n",
       "      <td>Barton County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-94.38276 37.64914, -94.37935 37.649...</td>\n",
       "      <td>Missouri</td>\n",
       "      <td>MO</td>\n",
       "      <td>29</td>\n",
       "      <td>Barton County, Missouri</td>\n",
       "      <td>Springfield-Cape Girardeau</td>\n",
       "      <td>St. Louis</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>29</td>\n",
       "      <td>083</td>\n",
       "      <td>00758496</td>\n",
       "      <td>29083</td>\n",
       "      <td>Henry</td>\n",
       "      <td>Henry County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-93.51565 38.37395, -93.51565 38.373...</td>\n",
       "      <td>Missouri</td>\n",
       "      <td>MO</td>\n",
       "      <td>29</td>\n",
       "      <td>Henry County, Missouri</td>\n",
       "      <td>Kansas City-St. Joseph</td>\n",
       "      <td>St. Louis</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>91</th>\n",
       "      <td>29</td>\n",
       "      <td>103</td>\n",
       "      <td>00758506</td>\n",
       "      <td>29103</td>\n",
       "      <td>Knox</td>\n",
       "      <td>Knox County</td>\n",
       "      <td>06</td>\n",
       "      <td>POLYGON ((-91.95166 40.21395, -91.94724 40.213...</td>\n",
       "      <td>Missouri</td>\n",
       "      <td>MO</td>\n",
       "      <td>29</td>\n",
       "      <td>Knox County, Missouri</td>\n",
       "      <td>Jefferson City</td>\n",
       "      <td>St. Louis</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>79</th>\n",
       "      <td>78</td>\n",
       "      <td>030</td>\n",
       "      <td>02378250</td>\n",
       "      <td>78030</td>\n",
       "      <td>St. Thomas</td>\n",
       "      <td>St. Thomas Island</td>\n",
       "      <td>10</td>\n",
       "      <td>POLYGON ((-65.15388 18.28114, -65.15364 18.293...</td>\n",
       "      <td>Virgin Islands</td>\n",
       "      <td>VI</td>\n",
       "      <td>78</td>\n",
       "      <td>St. Thomas Island, Virgin Islands</td>\n",
       "      <td>Saint Thomas</td>\n",
       "      <td>Washington</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>613</th>\n",
       "      <td>11</td>\n",
       "      <td>001</td>\n",
       "      <td>01702382</td>\n",
       "      <td>11001</td>\n",
       "      <td>District of Columbia</td>\n",
       "      <td>District of Columbia</td>\n",
       "      <td>00</td>\n",
       "      <td>POLYGON ((-77.11975 38.93435, -77.11886 38.935...</td>\n",
       "      <td>District of Columbia</td>\n",
       "      <td>DC</td>\n",
       "      <td>11</td>\n",
       "      <td>District of Columbia, District of Columbia</td>\n",
       "      <td>Washington</td>\n",
       "      <td>Washington</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      STATEFP COUNTYFP  COUNTYNS  GEOID                  NAME  \\\n",
       "12         48      011  01383791  48011             Armstrong   \n",
       "37         48      077  01383824  48077                  Clay   \n",
       "10         48      189  01383880  48189                  Hale   \n",
       "6          48      327  01383949  48327                Menard   \n",
       "1206       48      377  01383974  48377              Presidio   \n",
       "120        48      247  01383909  48247              Jim Hogg   \n",
       "359        15      003  00365281  15003              Honolulu   \n",
       "618         6      075  00277302  06075         San Francisco   \n",
       "1885        6      085  00277307  06085           Santa Clara   \n",
       "340         6      009  01675885  06009             Calaveras   \n",
       "149        32      003  00863600  32003                 Clark   \n",
       "444         6      097  01657246  06097                Sonoma   \n",
       "23         49      033  01448030  49033                  Rich   \n",
       "1098        6      013  01675903  06013          Contra Costa   \n",
       "7           6      091  00277310  06091                Sierra   \n",
       "81         32      001  00858638  32001             Churchill   \n",
       "2          35      011  00933054  35011               De Baca   \n",
       "29         35      035  00929104  35035                 Otero   \n",
       "391         4      027  00023901  04027                  Yuma   \n",
       "30         35      003  00929108  35003                Catron   \n",
       "832         4      013  00037026  04013              Maricopa   \n",
       "1129       53      013  01513273  53013              Columbia   \n",
       "1168       53      037  01531926  53037              Kittitas   \n",
       "1          53      069  01513275  53069             Wahkiakum   \n",
       "727        29      221  00758563  29221            Washington   \n",
       "33         29      011  00758460  29011                Barton   \n",
       "20         29      083  00758496  29083                 Henry   \n",
       "91         29      103  00758506  29103                  Knox   \n",
       "79         78      030  02378250  78030            St. Thomas   \n",
       "613        11      001  01702382  11001  District of Columbia   \n",
       "\n",
       "                  NAMELSAD LSAD  \\\n",
       "12        Armstrong County   06   \n",
       "37             Clay County   06   \n",
       "10             Hale County   06   \n",
       "6            Menard County   06   \n",
       "1206       Presidio County   06   \n",
       "120        Jim Hogg County   06   \n",
       "359        Honolulu County   06   \n",
       "618   San Francisco County   06   \n",
       "1885    Santa Clara County   06   \n",
       "340       Calaveras County   06   \n",
       "149           Clark County   06   \n",
       "444          Sonoma County   06   \n",
       "23             Rich County   06   \n",
       "1098   Contra Costa County   06   \n",
       "7            Sierra County   06   \n",
       "81        Churchill County   06   \n",
       "2           De Baca County   06   \n",
       "29            Otero County   06   \n",
       "391            Yuma County   06   \n",
       "30           Catron County   06   \n",
       "832        Maricopa County   06   \n",
       "1129       Columbia County   06   \n",
       "1168       Kittitas County   06   \n",
       "1         Wahkiakum County   06   \n",
       "727      Washington County   06   \n",
       "33           Barton County   06   \n",
       "20            Henry County   06   \n",
       "91             Knox County   06   \n",
       "79       St. Thomas Island   10   \n",
       "613   District of Columbia   00   \n",
       "\n",
       "                                               geometry            State_Name  \\\n",
       "12    POLYGON ((-101.08728 35.02708, -101.08729 35.0...                 Texas   \n",
       "37    POLYGON ((-97.97784 33.81417, -97.97784 33.813...                 Texas   \n",
       "10    POLYGON ((-101.73504 33.82864, -101.73507 33.8...                 Texas   \n",
       "6     POLYGON ((-99.77120 30.71079, -99.78240 30.710...                 Texas   \n",
       "1206  POLYGON ((-103.79577 29.72453, -103.79577 29.7...                 Texas   \n",
       "120   POLYGON ((-98.49291 27.17469, -98.49292 27.174...                 Texas   \n",
       "359   MULTIPOLYGON (((-166.38653 23.84636, -166.3842...                Hawaii   \n",
       "618   MULTIPOLYGON (((-122.40537 37.70826, -122.4054...            California   \n",
       "1885  POLYGON ((-121.47115 36.97797, -121.47121 36.9...            California   \n",
       "340   POLYGON ((-120.37620 38.14265, -120.37616 38.1...            California   \n",
       "149   POLYGON ((-115.89692 36.84208, -115.89159 36.8...                Nevada   \n",
       "444   POLYGON ((-122.63498 38.58301, -122.63496 38.5...            California   \n",
       "23    POLYGON ((-111.04653 41.77964, -111.04651 41.7...                  Utah   \n",
       "1098  POLYGON ((-121.91046 37.73039, -121.91057 37.7...            California   \n",
       "7     POLYGON ((-120.55587 39.50874, -120.55614 39.5...            California   \n",
       "81    POLYGON ((-118.18700 40.00209, -118.17322 40.0...                Nevada   \n",
       "2     POLYGON ((-104.38368 34.69213, -104.37658 34.6...            New Mexico   \n",
       "29    POLYGON ((-105.93164 33.38990, -105.93089 33.3...            New Mexico   \n",
       "391   POLYGON ((-114.47325 33.02788, -114.45989 33.0...               Arizona   \n",
       "30    POLYGON ((-108.20839 33.20097, -108.21384 33.2...            New Mexico   \n",
       "832   POLYGON ((-111.89308 33.20471, -111.90205 33.2...               Arizona   \n",
       "1129  POLYGON ((-117.73737 46.44933, -117.73742 46.4...            Washington   \n",
       "1168  POLYGON ((-120.00214 46.99732, -120.00193 46.9...            Washington   \n",
       "1     POLYGON ((-123.49077 46.38358, -123.48813 46.3...            Washington   \n",
       "727   POLYGON ((-90.81877 37.73553, -90.82091 37.735...              Missouri   \n",
       "33    POLYGON ((-94.38276 37.64914, -94.37935 37.649...              Missouri   \n",
       "20    POLYGON ((-93.51565 38.37395, -93.51565 38.373...              Missouri   \n",
       "91    POLYGON ((-91.95166 40.21395, -91.94724 40.213...              Missouri   \n",
       "79    POLYGON ((-65.15388 18.28114, -65.15364 18.293...        Virgin Islands   \n",
       "613   POLYGON ((-77.11975 38.93435, -77.11886 38.935...  District of Columbia   \n",
       "\n",
       "     State_Code  FIPS                                county_state  \\\n",
       "12           TX    48                     Armstrong County, Texas   \n",
       "37           TX    48                          Clay County, Texas   \n",
       "10           TX    48                          Hale County, Texas   \n",
       "6            TX    48                        Menard County, Texas   \n",
       "1206         TX    48                      Presidio County, Texas   \n",
       "120          TX    48                      Jim Hogg County, Texas   \n",
       "359          HI    15                     Honolulu County, Hawaii   \n",
       "618          CA     6            San Francisco County, California   \n",
       "1885         CA     6              Santa Clara County, California   \n",
       "340          CA     6                Calaveras County, California   \n",
       "149          NV    32                        Clark County, Nevada   \n",
       "444          CA     6                   Sonoma County, California   \n",
       "23           UT    49                           Rich County, Utah   \n",
       "1098         CA     6             Contra Costa County, California   \n",
       "7            CA     6                   Sierra County, California   \n",
       "81           NV    32                    Churchill County, Nevada   \n",
       "2            NM    35                  De Baca County, New Mexico   \n",
       "29           NM    35                    Otero County, New Mexico   \n",
       "391          AZ     4                        Yuma County, Arizona   \n",
       "30           NM    35                   Catron County, New Mexico   \n",
       "832          AZ     4                    Maricopa County, Arizona   \n",
       "1129         WA    53                 Columbia County, Washington   \n",
       "1168         WA    53                 Kittitas County, Washington   \n",
       "1            WA    53                Wahkiakum County, Washington   \n",
       "727          MO    29                 Washington County, Missouri   \n",
       "33           MO    29                     Barton County, Missouri   \n",
       "20           MO    29                      Henry County, Missouri   \n",
       "91           MO    29                       Knox County, Missouri   \n",
       "79           VI    78           St. Thomas Island, Virgin Islands   \n",
       "613          DC    11  District of Columbia, District of Columbia   \n",
       "\n",
       "                         Diocese       Province  \n",
       "12                      Amarillo    San Antonio  \n",
       "37                    Fort Worth    San Antonio  \n",
       "10                       Lubbock    San Antonio  \n",
       "6                     San Angelo    San Antonio  \n",
       "1206                     El Paso    San Antonio  \n",
       "120                       Laredo    San Antonio  \n",
       "359                     Honolulu  San Francisco  \n",
       "618                San Francisco  San Francisco  \n",
       "1885      San Jose in California  San Francisco  \n",
       "340                     Stockton  San Francisco  \n",
       "149                    Las Vegas  San Francisco  \n",
       "444     Santa Rosa in California  San Francisco  \n",
       "23                Salt Lake City  San Francisco  \n",
       "1098                     Oakland  San Francisco  \n",
       "7                     Sacramento  San Francisco  \n",
       "81                          Reno  San Francisco  \n",
       "2                       Santa Fe       Santa Fe  \n",
       "29                    Las Cruces       Santa Fe  \n",
       "391                       Tucson       Santa Fe  \n",
       "30                        Gallup       Santa Fe  \n",
       "832                      Phoenix       Santa Fe  \n",
       "1129                     Spokane        Seattle  \n",
       "1168                      Yakima        Seattle  \n",
       "1                        Seattle        Seattle  \n",
       "727                    St. Louis      St. Louis  \n",
       "33    Springfield-Cape Girardeau      St. Louis  \n",
       "20        Kansas City-St. Joseph      St. Louis  \n",
       "91                Jefferson City      St. Louis  \n",
       "79                  Saint Thomas     Washington  \n",
       "613                   Washington     Washington  "
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "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": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 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.to_file('diocese_boundaries.geojson', driver = 'GeoJSON')\n",
    "diocese_boundaries.head(5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 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.to_file('counties_by_diocese.geojson', driver = 'GeoJSON')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Inaccurate parts of the map at this point:\n",
    "\n",
    "1. Fall River/Boston: The Diocese of Fall River includes \"the towns of Marion, Mattapoisett, and Wareham along the south coast of Plymouth County in Massachusetts, in the New England region of the United States.\" (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Fall_River) The rest of Plymouth County belongs to the Archdiocese of Boston. \n",
    "\n",
    "2. Austin/Victoria in Texas: Only part of Fayette County belongs to the Diocese of Victoria in Texas. \"The part of Fayette County north of the Colorado River\" is part of the Diocese of Austin. (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Austin)\n",
    "\n",
    "3. San Antonio/Corpus Christi: The Archdiocese of San Antonio includes only part of McMullen County. The rest belongs to the Diocese of Corpus Christi. See: https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Corpus_Christi \n",
    "This also means that McMullen County is part of two different provinces (Galveston-Houston and San Antonio).\n",
    "\n",
    "4. Evansville/Indianapolis: Wikipedia notes that \"Harrison Township in Spencer County, the location of St. Meinrad Archabbey, is part of the Archdiocese of Indianapolis.\" (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Evansville)\n",
    "\n",
    "5. Gallup/Santa Fe: Wikipedia notes that the Diocsee of Gallup also includes \"parts of Rio Arriba, Sandoval, Bernalillo, and Valencia Counties west of 106,52',41\" meridian in New Mexico.\" The rest of those counties belong to Santa Fe. (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Gallup)\n",
    "\n",
    "6. Tucson/Phoenix: The parts of Pinal County that lie within \"the territorial boundaries of the Gila River Indian Community\" are part of the Diocese of Phoenix, not the Diocese of Tucson. (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Tucson)\n",
    "\n",
    "7. New Orleans/Houma-Thibodaux: The Archdiocese of New Orleans does not include Grand Isle (within Jefferson Parish); that instead belongs to the Diocese of Houma-Thibodaux. (https://en.wikipedia.org/wiki/Roman_Catholic_Archdiocese_of_New_Orleans)\n",
    "\n",
    "8. Lafayette in Louisiana/Houma-Thibodaux: The eastern portion of St. Mary Parish (including Morgan City) belongs to the Diocese of Houma-Thibodaux; the rest belongs to the Diocese of Lafayette in Louisiana.\n",
    "\n",
    "9. Rockville Centre/Norwich: Fisher Island (within Suffolk County, NY) belongs to the Diocese of Norwich, not the Diocese of Rockville Centre.\n",
    "\n",
    "10. Albany/Ogdensburg: The northern part of Herkimer County belongs to the Diocese of Ogdensburg; the rest belongs to the Diocese of Albany. (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Ogdensburg\n",
    "\n",
    "11. Albany/Ogdensburg: The southern part of Hamilton County belongs to the Diocese of Albany; the rest belongs to the Diocese of Ogdensburg. (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Ogdensburg)\n",
    "\n",
    "13. Grand Island/Lincoln: the parts of Deuel, Keith, Lincoln, Dawson, and Hall counties that lie south of the Platte River belong to the Diocese of Lincoln. The northern parts of these counties belong to the Diocese of Grand Island.\n",
    "Although the Archdiocese of Omaha reports that \"parts of Hall County\" belong to its diocese (https://archomaha.org/about-us/), I believe this is incorrect, as the Diocese of Lincoln's parish map shows the southern portion of Hall County as belonging instead instead to the Diocsee of Lincoln: https://www.lincolndiocese.org/directory/parish-map\n",
    "For now, I have assigned Lincoln County to the Diocese of Lincoln and Deuel, Keith, Dawson, and Hall counties to the Diocese of Grand Island (based on which diocese has the majority of each's county's area). \n",
    "It appears that a small portion of Buffalo County may lie south of the Platte, but I'm keeping all of this county within the Diocese of Grand Island for now, since it's possible the Platte's course changed following the diocese's establishment (and the bottom edge of the county does appear to follow the course of a river).\n",
    "\n",
    "14. Helena/Great Falls-Billings: The Diocese of Great Falls-Billings includes parts of 2 counties, given that the Diocese of Helena does also: https://diocesehelena.org/about/ Based on the Diocese of Helena's Wikipedia map (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Helena#/media/File:Diocese_of_Helena.PNG), it looks like the Diocese of Great Falls-Billings includes portions of Toole and Teton counties, the majority of which belong to the Diocese of Helena.\n",
    "However, the Diocese of Great Falls-Billings map doesn't appear to reflect these modifications:\n",
    "https://diocesegfb.org/home-page/about/\n",
    "Meanwhile, the Diocese of Helena's wikipedia page reports that the diocese includes \"Parts of Meagher, Musselshell, and Toole counties in Montana.\" However, Musselshell county is landlocked within the Diocese of Great Falls-Billings, so this is likely a typo. I also don't see any evidence on the Wikipedia map (https://en.wikipedia.org/wiki/Roman_Catholic_Diocese_of_Helena#/media/File:Diocese_of_Helena.PNG) of part of Meagher county belonging to the Diocsee of Great Falls-Billings.\n",
    "\n",
    "15. Bismark/Fargo: The western section of Bottineau County belongs to the Diocese of Bismark; the rest belongs to the Diocese of Fargo. See https://bismarckdiocese.com/about and page 3 of https://d2y1pz2y630308.cloudfront.net/2950/documents/2022/3/DiocesanDirectory.pdf."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "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_cd = df_dioceses.merge(df_diocese_numbers, on = 'Diocese', \n",
    "how = 'left') # the 'cd' in df_cd stands for 'counties' and 'dioceses'\n",
    "\n",
    "# Change df_cd to df_dioceses once you have all the colors right"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Updating Color Codes:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# # This use of df.apply() with if/else statements is based on an example\n",
    "# by Professor Hardeep Johar.\n",
    "\n",
    "# Green:\n",
    "df_cd['Diocese_Color_Code'] = df_cd.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_cd['Diocese_Color_Code'] = df_cd.apply(lambda row: 0 if row['Diocese'] in [\n",
    "    'Santa Rosa in California', 'Crookston', 'Kansas City-St. Joseph', \n",
    "    'Gaylord'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "#Pink:\n",
    "df_cd['Diocese_Color_Code'] = df_cd.apply(lambda row: 64 if row['Diocese'] in [\n",
    "    'Reno', 'San Antonio'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "#Red:\n",
    "df_cd['Diocese_Color_Code'] = df_cd.apply(lambda row: 80 if row['Diocese'] in [\n",
    "    'Cheyenne', 'Sioux Falls', 'Baltimore', 'San Francisco']\n",
    "    else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Blue:\n",
    "df_cd['Diocese_Color_Code'] = df_cd.apply(lambda row: 16 if row['Diocese'] in [\n",
    "    'Davenport', 'Tulsa', 'San Diego', \n",
    "    'Biloxi'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Light green:\n",
    "df_cd['Diocese_Color_Code'] = df_cd.apply(lambda row: 32 if row['Diocese'] in [\n",
    "    'Sacramento', 'Washington', 'Palm Beach', 'Birmingham',\n",
    "    'Rockville Centre'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Yellow:\n",
    "df_cd['Diocese_Color_Code'] = df_cd.apply(lambda row: 160 if row['Diocese'] in [\n",
    "    'Lubbock'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Light orange:\n",
    "df_cd['Diocese_Color_Code'] = df_cd.apply(lambda row: 96 if row['Diocese'] in [\n",
    "    'Youngstown', 'Denver'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Orange:\n",
    "df_cd['Diocese_Color_Code'] = df_cd.apply(lambda row: 112 if row['Diocese'] in [\n",
    "    'Alexandria in Louisiana', 'Covington', 'Helena', 'Raleigh'] \n",
    "    else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Purple:\n",
    "df_cd['Diocese_Color_Code'] = df_cd.apply(lambda row: 144 if row['Diocese'] in [\n",
    "    'Lafayette in Louisiana', 'Lexington', 'Gary', 'New York', \n",
    "    'St. Louis'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Brown:\n",
    "df_cd['Diocese_Color_Code'] = df_cd.apply(lambda row: 174 if row['Diocese'] in [\n",
    "    'Louisville', 'Bridgeport', 'Milwaukee'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n",
    "# Light purple:\n",
    "df_cd['Diocese_Color_Code'] = df_cd.apply(lambda row: 128 if row['Diocese'] in [\n",
    "    'Saint Cloud', 'Memphis'] else row['Diocese_Color_Code'], axis = 1)\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_cd.to_csv('counties_by_diocese.csv')"
   ]
  },
  {
   "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_lat'] = ''\n",
    "diocese_boundaries['center_lon'] = ''\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "for i in range(len(diocese_boundaries)):\n",
    "    geometry = str(diocese_boundaries.iloc[i,diocese_boundaries.columns.get_loc('geometry')])\n",
    "    # print(geometry)\n",
    "    if ('MULTIPOLYGON') not in geometry:\n",
    "        coord_list = geometry.replace('POLYGON ((', '').replace('))', '').replace(')', '').replace('(', '').split(', ')\n",
    "        lon_list = [float(entry.split(' ')[0]) for entry in coord_list]\n",
    "        lat_list = [float(entry.split(' ')[1]) for entry in coord_list]\n",
    "        \n",
    "        max_lat = np.max(lat_list)\n",
    "        min_lat = np.min(lat_list)\n",
    "\n",
    "        max_lon = np.max(lon_list)\n",
    "        min_lon = np.min(lon_list)\n",
    "        \n",
    "        center_lat = (max_lat + min_lat) / 2\n",
    "        center_lon = (max_lon + min_lon) /2\n",
    "        \n",
    "        diocese_boundaries.iloc[i,diocese_boundaries.columns.get_loc('center_lat')] = center_lat\n",
    "        diocese_boundaries.iloc[i,diocese_boundaries.columns.get_loc('center_lon')] = center_lon"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "diocese_boundaries['center'] = diocese_boundaries['geometry'].representative_point()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "center_point = [float(diocese_boundaries.iloc[i, diocese_boundaries.columns.get_loc('center')].x),\n",
    "float(diocese_boundaries.iloc[i, diocese_boundaries.columns.get_loc('center')].y)] \n",
    "center_point"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "diocese_map = folium.Map(location=[38.7, -95], zoom_start=6, tiles = 'Stamen Toner')\n",
    "\n",
    "# Although Folium has a choropleth library, I wasn't able to find a way\n",
    "# to disable the default legend. Therefore, I am instead using a custom\n",
    "# choropleth mapping function. Much of this function is based on \n",
    "# Amodiovalerio Verde's code at:\n",
    "# https://vverde.github.io/blob/interactivechoropleth.html . \n",
    "\n",
    "colormap = branca.colormap.linear.Paired_12.to_step(\n",
    "    12, data = df_cd['Diocese_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",
    "        \"fillOpacity\": 0.5,\n",
    "        \"weight\": 0.5,\n",
    "        'color': '#777777',\n",
    "        \"fillColor\": 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",
    "        'fillOpacity': 0.8}\n",
    "\n",
    "folium.features.GeoJson(df_cd, \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",
    "        '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_cd, 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",
    "    try:\n",
    "        folium.Marker(center_point, icon = folium.DivIcon(html=f\"<b><center>{diocese_boundaries.index[i]}</center></b>\")).add_to(diocese_map)\n",
    "# The strategy of using an invisible folium.Marker object to display text on \n",
    "# the map 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",
    "    except:\n",
    "        continue\n",
    "\n",
    "folium.LatLngPopup().add_to(diocese_map)\n",
    "\n",
    "diocese_map.save('diocese_map.html')\n",
    "\n",
    "# colormap.add_to(m)\n",
    "\n",
    "diocese_map"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "colormap"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "color_code_list = list(np.arange(0, 177, 176/11))\n",
    "color_code_list"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "color_list = list(np.arange(0, 177, 176/11))\n",
    "[colormap(color_list[i]) for i in range(len(color_list))]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "def create_map_screenshot(absolute_path_to_map_folder, map_name, \n",
    "screenshot_save_path = None):\n",
    "\n",
    "    '''\n",
    "    This function uses the Selenium library to create a screenshot \n",
    "    of a map so that it can be shared as a .png file.\n",
    "    See https://www.selenium.dev/documentation/ for more information on \n",
    "    Selenium. \n",
    "    \n",
    "    absolute_path_to_map_folder designates the absolute path where\n",
    "    the map is stored. (I wasn't able to get this code to work using\n",
    "    just relative paths.) \n",
    "\n",
    "    map_name specifies the name of the map, including its extension.\n",
    "\n",
    "    screenshot_save_path designates the folder where you wish to save\n",
    "    the map screenshot. This can be a relative path.\n",
    "\n",
    "    Note that some setup work is required for the Selenium code\n",
    "    to run correctly; if you don't have time right now to complete this \n",
    "    setup, you can comment out any code that calls this function.\n",
    "    '''\n",
    "\n",
    "    ff_driver = webdriver.Firefox() \n",
    "    # See https://www.selenium.dev/documentation/webdriver/getting_started/open_browser/\n",
    "    # For more information on using Selenium to get screenshots of .html \n",
    "    # files, see my get_screenshots.ipynb file within my route_maps_builder\n",
    "    # program, available here:\n",
    "    # https://github.com/kburchfiel/route_maps_builder/blob/master/get_screenshots.ipynb\n",
    "    window_width = 3000 # This produces a large window that can better\n",
    "    # capture small details (such as zip code shapefiles).\n",
    "    ff_driver.set_window_size(window_width,window_width*(9/16)) # Creates\n",
    "    # a window with an HD/4K/8K aspect ratio\n",
    "    ff_driver.get(f'{absolute_path_to_map_folder}\\\\{map_name}') \n",
    "    # See https://www.selenium.dev/documentation/webdriver/browser/navigation/\n",
    "    time.sleep(2) # This gives the page sufficient\n",
    "    # time to load the map tiles before the screenshot is taken. \n",
    "    # You can also experiment with longer sleep times.\n",
    "\n",
    "    if screenshot_save_path != None:\n",
    "        # If specifying a screenshot save path, you must create this path\n",
    "        # within your directory before the function is run; otherwise,\n",
    "        # it won't return an image. \n",
    "        ff_driver.get_screenshot_as_file(\n",
    "            screenshot_save_path+'/'+map_name.replace('.html','')+'.png') \n",
    "    else: # If no save path was specified for the screenshot, the image\n",
    "        # will be saved within the project's root folder.\n",
    "        ff_driver.get_screenshot_as_file(\n",
    "            map_name.replace('.html','')+'.png') \n",
    "    # Based on:\n",
    "    # https://www.selenium.dev/selenium/docs/api/java/org/openqa/selenium/TakesScreenshot.html\n",
    "\n",
    "    ff_driver.quit()\n",
    "    # Based on: https://www.selenium.dev/documentation/webdriver/browser/windows/\n",
    "\n",
    "def convert_png_to_smaller_jpg(png_folder, png_image_name, jpg_folder, \n",
    "reduction_factor = 1, quality_factor = 50):\n",
    "    ''' This function converts a .png image into a smaller .jpg image, which\n",
    "    helps reduce file sizes and load times when displaying a series of images\n",
    "    within a notebook or online.\n",
    "    png_folder and png_image_name specify the location of the original .png\n",
    "    image.\n",
    "    jpg folder specifies the location where the .jpg screenshot should be\n",
    "    saved.\n",
    "    reduction_factor specifies the amount by which you would like to reduce\n",
    "    the image's dimensions. For instance, to convert a 4K (3840*2160) image\n",
    "    to a full HD (1920*1080) one, use a reduction factor of 2. If you do not\n",
    "    wish to reduce the image's size, use the default reduction factor of 1.\n",
    "    '''\n",
    "    with PIL.Image.open(f'{png_folder}/{png_image_name}') as map_image:\n",
    "        (width, height) = (map_image.width // reduction_factor, \n",
    "        map_image.height // reduction_factor)\n",
    "        jpg_image = map_image.resize((width, height))\n",
    "        # The above code is based on:\n",
    "        # https://pillow.readthedocs.io/en/stable/reference/Image.html\n",
    "        jpg_image = jpg_image.convert('RGB')\n",
    "        # The above conversion is necessary in order to save .png files as \n",
    "        # .jpg files. It's based on Patrick Artner's answer at: \n",
    "        # https://stackoverflow.com/a/48248432/13097194\n",
    "        jpg_image_name = png_image_name.replace('png', 'jpg') \n",
    "        jpg_image.save(f'{jpg_folder}/{jpg_image_name}',\n",
    "        format = 'JPEG', quality = quality_factor, optimize = True)\n",
    "        # See https://pillow.readthedocs.io/en/stable/handbook/image-file-formats.html#jpeg"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "create_map_screenshot(\n",
    "    absolute_path_to_map_folder = \n",
    "    r'C:\\Users\\kburc\\D1V1\\Documents\\!Dell64docs\\Programming\\py\\kjb3_programs\\us_diocese_mapper', \n",
    "    map_name = 'diocese_map.html', screenshot_save_path = 'screenshots')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "convert_png_to_smaller_jpg(png_folder = 'screenshots', \n",
    "png_image_name = 'diocese_map.png', jpg_folder = 'smaller_screenshots', \n",
    "quality_factor = 80)\n",
    "\n",
    "IPython.display.Image(filename='smaller_screenshots/'+'diocese_map.jpg') \n",
    "# Based on https://ipython.readthedocs.io/en/stable/api/generated/IPython.display.html"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# m = folium.Map(location=[38.7, -95], zoom_start=6, tiles = 'Stamen Toner')\n",
    "\n",
    "# style_function = lambda x: {'weight': '1', 'color': '#000000', 'fillOpacity': 0} \n",
    "# # See https://python-visualization.github.io/folium/modules.html#folium.vector_layers.path_options\n",
    "# # for style_function arguments.\n",
    "\n",
    "# folium.GeoJson(merged_county_shapes, name = \"geojson\", style_function = style_function).add_to(m)\n",
    "\n",
    "# m"
   ]
  },
  {
   "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
}