Plot Station Data from Map.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Objective\n",
    "\n",
    "Provide a map of current and historic gaging stations and plot data for a selected station. The individual sites will be shown as blue markers once you zoom into the area of intest. Hovering over a marker provides the station name and USGS site number. Clicking on the marker will plot available discharge and sediment concentration data. Note that for now the data is being pulled from the sediment portal, which doesn't have all the station data. Thus, some stations may not plot any information even though the USGS site information shows the site as having data available."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "4f7373dbbc2b4450b43bb12eb51acb78",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "Map(center=[30.661, -88.645], controls=(ZoomControl(options=['position', 'zoom_in_text', 'zoom_in_title', 'zoo…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "  agency_cd  site_no   datetime  DAILY_FLOW DAILY_FLOW_QUAL  DAILY_SSC  \\\n",
      "0      USGS  7010000 1861-01-01       25200               A        NaN   \n",
      "1      USGS  7010000 1861-01-02       27400               A        NaN   \n",
      "2      USGS  7010000 1861-01-03       27400               A        NaN   \n",
      "3      USGS  7010000 1861-01-04       27400               A        NaN   \n",
      "4      USGS  7010000 1861-01-05       29000               A        NaN   \n",
      "\n",
      "  DAILY_SSC_QUAL  DAILY_SSL DAILY_SSL_QUAL  \n",
      "0            NaN        NaN            NaN  \n",
      "1            NaN        NaN            NaN  \n",
      "2            NaN        NaN            NaN  \n",
      "3            NaN        NaN            NaN  \n",
      "4            NaN        NaN            NaN  \n",
      "  agency_cd  site_no   datetime  DAILY_FLOW DAILY_FLOW_QUAL  DAILY_SSC  \\\n",
      "0      USGS  7010000 1861-01-01       25200               A        NaN   \n",
      "1      USGS  7010000 1861-01-02       27400               A        NaN   \n",
      "2      USGS  7010000 1861-01-03       27400               A        NaN   \n",
      "3      USGS  7010000 1861-01-04       27400               A        NaN   \n",
      "4      USGS  7010000 1861-01-05       29000               A        NaN   \n",
      "\n",
      "  DAILY_SSC_QUAL  DAILY_SSL DAILY_SSL_QUAL  \n",
      "0            NaN        NaN            NaN  \n",
      "1            NaN        NaN            NaN  \n",
      "2            NaN        NaN            NaN  \n",
      "3            NaN        NaN            NaN  \n",
      "4            NaN        NaN            NaN  \n"
     ]
    }
   ],
   "source": [
    "import matplotlib\n",
    "%matplotlib ipympl\n",
    "from matplotlib import pyplot as plt\n",
    "import matplotlib.dates as mdates\n",
    "from ipyleaflet import Map, Marker, MarkerCluster, WMSLayer, Unicode, basemaps\n",
    "import numpy as np\n",
    "\n",
    "wms = WMSLayer(\n",
    "    url='https://basemap.nationalmap.gov/arcgis/rest/services/USGSImageryTopo/MapServer/tile/{z}/{y}/{x}',\n",
    "    attribution='US Geological Survey'\n",
    ")\n",
    "\n",
    "m = Map(basemap=basemaps.CartoDB.Positron, center=(30.661, -88.645), zoom=5)\n",
    "m.add_layer(wms)\n",
    "m"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "54b1a4eb79804677bb38e7de431a11e0",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig0,ax0 = plt.subplots()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "92186dbe3f414640902ec5673d5610c0",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig1,ax1 = plt.subplots()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "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>agency_cd</th>\n",
       "      <th>site_no</th>\n",
       "      <th>station_nm</th>\n",
       "      <th>site_tp_cd</th>\n",
       "      <th>dec_lat_va</th>\n",
       "      <th>dec_long_va</th>\n",
       "      <th>coord_acy_cd</th>\n",
       "      <th>dec_coord_datum_cd</th>\n",
       "      <th>alt_va</th>\n",
       "      <th>alt_acy_va</th>\n",
       "      <th>alt_datum_cd</th>\n",
       "      <th>huc_cd</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>USGS</td>\n",
       "      <td>2489500</td>\n",
       "      <td>Pearl River near Bogalusa, LA</td>\n",
       "      <td>ST</td>\n",
       "      <td>30.793243</td>\n",
       "      <td>-89.820907</td>\n",
       "      <td>U</td>\n",
       "      <td>NAD83</td>\n",
       "      <td>54.64</td>\n",
       "      <td>1.00</td>\n",
       "      <td>NAVD88</td>\n",
       "      <td>3180004.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>USGS</td>\n",
       "      <td>2490193</td>\n",
       "      <td>Pearl River at Pools Bluff near Bogalusa, LA</td>\n",
       "      <td>ST</td>\n",
       "      <td>30.705833</td>\n",
       "      <td>-89.846389</td>\n",
       "      <td>5</td>\n",
       "      <td>NAD83</td>\n",
       "      <td>0.32</td>\n",
       "      <td>0.16</td>\n",
       "      <td>NAVD88</td>\n",
       "      <td>3180004.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>USGS</td>\n",
       "      <td>2490200</td>\n",
       "      <td>Pearl River Navigation Canal Lock and Dam No. 3</td>\n",
       "      <td>ST</td>\n",
       "      <td>30.638056</td>\n",
       "      <td>-89.868056</td>\n",
       "      <td>5</td>\n",
       "      <td>NAD83</td>\n",
       "      <td>33.98</td>\n",
       "      <td>0.20</td>\n",
       "      <td>NAVD88</td>\n",
       "      <td>3180004.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>USGS</td>\n",
       "      <td>2491500</td>\n",
       "      <td>Bogue Chitto River at Franklinton, LA</td>\n",
       "      <td>ST</td>\n",
       "      <td>30.842962</td>\n",
       "      <td>-90.162030</td>\n",
       "      <td>U</td>\n",
       "      <td>NAD83</td>\n",
       "      <td>123.81</td>\n",
       "      <td>0.01</td>\n",
       "      <td>NGVD29</td>\n",
       "      <td>3180005.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>USGS</td>\n",
       "      <td>2492000</td>\n",
       "      <td>Bogue Chitto River near Bush, LA</td>\n",
       "      <td>ST</td>\n",
       "      <td>30.629356</td>\n",
       "      <td>-89.897296</td>\n",
       "      <td>U</td>\n",
       "      <td>NAD83</td>\n",
       "      <td>43.31</td>\n",
       "      <td>0.22</td>\n",
       "      <td>NAVD88</td>\n",
       "      <td>3180005.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  agency_cd  site_no                                       station_nm  \\\n",
       "0      USGS  2489500                    Pearl River near Bogalusa, LA   \n",
       "1      USGS  2490193     Pearl River at Pools Bluff near Bogalusa, LA   \n",
       "2      USGS  2490200  Pearl River Navigation Canal Lock and Dam No. 3   \n",
       "3      USGS  2491500            Bogue Chitto River at Franklinton, LA   \n",
       "4      USGS  2492000                 Bogue Chitto River near Bush, LA   \n",
       "\n",
       "  site_tp_cd  dec_lat_va  dec_long_va coord_acy_cd dec_coord_datum_cd  alt_va  \\\n",
       "0         ST   30.793243   -89.820907            U              NAD83   54.64   \n",
       "1         ST   30.705833   -89.846389            5              NAD83    0.32   \n",
       "2         ST   30.638056   -89.868056            5              NAD83   33.98   \n",
       "3         ST   30.842962   -90.162030            U              NAD83  123.81   \n",
       "4         ST   30.629356   -89.897296            U              NAD83   43.31   \n",
       "\n",
       "   alt_acy_va alt_datum_cd     huc_cd  \n",
       "0        1.00       NAVD88  3180004.0  \n",
       "1        0.16       NAVD88  3180004.0  \n",
       "2        0.20       NAVD88  3180004.0  \n",
       "3        0.01       NGVD29  3180005.0  \n",
       "4        0.22       NAVD88  3180005.0  "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import requests,pandas,io\n",
    "state_data=[]\n",
    "for state in ['la','ms','ar','al','tn','mo']:\n",
    "    resp = requests.post('https://waterservices.usgs.gov/nwis/site/?format=rdb&stateCd={0}&siteType=ST&siteStatus=active&hasDataTypeCd=iv,dv'.format(state))\n",
    "    ncomments=0\n",
    "    for line in resp.text.split('\\n'):\n",
    "        if '#' in line:\n",
    "            ncomments+=1\n",
    "        else:\n",
    "            break\n",
    "    state_data.append(pandas.read_table(io.BytesIO(resp._content),comment='#',skiprows=[ncomments+1]))\n",
    "f = pandas.concat(state_data, axis=0)\n",
    "f.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "def plot_station_daily(site_no):\n",
    "    import requests \n",
    "    import xml.etree.ElementTree as ET\n",
    "    from zipfile import ZipFile\n",
    "    from io import BytesIO\n",
    "    import numpy as np\n",
    "    # defining the api-endpoint  \n",
    "    URL = \"https://cida.usgs.gov/\"\n",
    "    API_ENDPOINT = \"https://cida.usgs.gov/sediment/data\"\n",
    "    ogcfilter = ET.Element(\"ogc:Filter\",attrib={\"xmlns:ogc\":\"http://www.opengis.net/ogc\"})\n",
    "    property_is = ET.SubElement(ogcfilter, \"ogc:PropertyIsEqualTo\", attrib={\"matchCase\":\"true\"})\n",
    "    ET.SubElement(property_is, \"ogc:PropertyName\").text = \"SITE_NO\"\n",
    "    ET.SubElement(property_is, \"ogc:Literal\").text = '{0:08d}'.format(site_no)\n",
    "    ogcfilter_string = ET.tostring(ogcfilter).decode()  \n",
    "    data = {'email' : 'cekees@gmail.com',\n",
    "            'format': 'tsv',\n",
    "            'directDownload':'true',\n",
    "            'dataTypes': 'sites_data_daily_discrete_',\n",
    "            'dailyFilter' : ogcfilter_string,\n",
    "            'discreteFilter' : ogcfilter_string\n",
    "           }\n",
    "    # sending post request and saving response as response object \n",
    "    search = requests.post(url = API_ENDPOINT, data = data)  \n",
    "    download = requests.get(url=URL+search.text)\n",
    "    zipfile = ZipFile(BytesIO(download._content))\n",
    "    def get_comments(filename):\n",
    "        ncomments=0\n",
    "        comments=''\n",
    "        with zipfile.open(filename) as f:\n",
    "            commentline = f.readline().decode()\n",
    "            while commentline:\n",
    "                if commentline[0] == '#' or commentline[:2] == r'\"#':\n",
    "                    ncomments +=1\n",
    "                    comments += commentline.strip(r'\"').strip('#')\n",
    "                    commentline=f.readline().decode()\n",
    "                else:\n",
    "                    commentline=None\n",
    "        return (ncomments,comments)   \n",
    "    if zipfile.filelist:\n",
    "        ax0.clear()\n",
    "        ax1.clear()\n",
    "        try:\n",
    "            (ncomments, comments) = get_comments('daily_data.tsv')\n",
    "            f = pandas.read_table(zipfile.open('daily_data.tsv'),header=ncomments).dropna(axis=1,how='all')\n",
    "            f.datetime = pandas.to_datetime(f.datetime)\n",
    "            print(f.head())\n",
    "            ax0.plot(f.datetime,f.DAILY_FLOW)\n",
    "            ax1.scatter(f.DAILY_FLOW,f.DAILY_SSC)\n",
    "            ax0.set_ylabel('Discharge [cfs]')\n",
    "            ax0.set_xlabel('Date')\n",
    "            ax0.set_title(\"Daily Discharge Data\")\n",
    "            ax1.set_xlabel('Discharge [cfs]')\n",
    "            ax1.set_ylabel('Sediment(SSC) [mg/l]')\n",
    "            ax1.set_title(\"Daily Sediment Data\")\n",
    "        except:\n",
    "            #print(\"daily flow failed\")\n",
    "            (ncomments, comments) = get_comments('discrete_data.tsv')\n",
    "            f = pandas.read_table(zipfile.open('discrete_data.tsv'),header=ncomments).dropna(axis=1,how='all')\n",
    "            f.DATETIME = pandas.to_datetime(f.DATETIME)\n",
    "            print(f.head())\n",
    "            try:\n",
    "                ax0.plot(f.DATETIME,f.INSTFLOW)\n",
    "            except:\n",
    "                #print(\"inst flow failed\")\n",
    "                pass\n",
    "            try:\n",
    "                ax0.plot(f.DATETIME,f.DAILYFLOW)\n",
    "                ssc = np.zeros(f.SSC.shape,'d')\n",
    "                for i,ssci in enumerate(f.SSC):\n",
    "                    if 'A' in ssci:\n",
    "                        ssc[i] = float(ssci[2:])\n",
    "                    else:\n",
    "                        ssc[i] = float(ssci)\n",
    "                ax1.scatter(f.DAILYFLOW,ssc)\n",
    "                if 'P63' in f.columns:    \n",
    "                    p63 = np.zeros(f.P63.shape,'d')\n",
    "                    for i,p63i in enumerate(f.P63):\n",
    "                        if type(p63i) is not float:\n",
    "                            if 'A' in p63i:\n",
    "                                p63[i] = float(p63i[2:])\n",
    "                            else:\n",
    "                                p63[i] = float(p63i)\n",
    "                        else:\n",
    "                            p63[i] = p63i\n",
    "                    ax1.scatter(f.DAILYFLOW,p63)               \n",
    "            except:\n",
    "                #print(\"daily flow data failed\")\n",
    "                pass\n",
    "            ax0.set_title(\"Discrete Discharge\")\n",
    "            ax0.set_ylabel('Discharge [cfs]')\n",
    "            times = np.array(f.DATETIME)\n",
    "            ax1.set_xlabel('Discharge [cfs]')\n",
    "            ax1.set_ylabel('Sediment(SSC) [mg/l]')\n",
    "            ax1.set_title(\"Discrete Sediment Data\")\n",
    "    else:\n",
    "        ax0.set_title(\"No Data\")\n",
    "        ax1.set_title(\"No Data\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "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>agency_cd</th>\n",
       "      <th>site_no</th>\n",
       "      <th>station_nm</th>\n",
       "      <th>site_tp_cd</th>\n",
       "      <th>dec_lat_va</th>\n",
       "      <th>dec_long_va</th>\n",
       "      <th>coord_acy_cd</th>\n",
       "      <th>dec_coord_datum_cd</th>\n",
       "      <th>alt_va</th>\n",
       "      <th>alt_acy_va</th>\n",
       "      <th>alt_datum_cd</th>\n",
       "      <th>huc_cd</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>USGS</td>\n",
       "      <td>2489500</td>\n",
       "      <td>Pearl River near Bogalusa, LA</td>\n",
       "      <td>ST</td>\n",
       "      <td>30.793243</td>\n",
       "      <td>-89.820907</td>\n",
       "      <td>U</td>\n",
       "      <td>NAD83</td>\n",
       "      <td>54.64</td>\n",
       "      <td>1.00</td>\n",
       "      <td>NAVD88</td>\n",
       "      <td>3180004.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>USGS</td>\n",
       "      <td>2490193</td>\n",
       "      <td>Pearl River at Pools Bluff near Bogalusa, LA</td>\n",
       "      <td>ST</td>\n",
       "      <td>30.705833</td>\n",
       "      <td>-89.846389</td>\n",
       "      <td>5</td>\n",
       "      <td>NAD83</td>\n",
       "      <td>0.32</td>\n",
       "      <td>0.16</td>\n",
       "      <td>NAVD88</td>\n",
       "      <td>3180004.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>USGS</td>\n",
       "      <td>2490200</td>\n",
       "      <td>Pearl River Navigation Canal Lock and Dam No. 3</td>\n",
       "      <td>ST</td>\n",
       "      <td>30.638056</td>\n",
       "      <td>-89.868056</td>\n",
       "      <td>5</td>\n",
       "      <td>NAD83</td>\n",
       "      <td>33.98</td>\n",
       "      <td>0.20</td>\n",
       "      <td>NAVD88</td>\n",
       "      <td>3180004.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>USGS</td>\n",
       "      <td>2491500</td>\n",
       "      <td>Bogue Chitto River at Franklinton, LA</td>\n",
       "      <td>ST</td>\n",
       "      <td>30.842962</td>\n",
       "      <td>-90.162030</td>\n",
       "      <td>U</td>\n",
       "      <td>NAD83</td>\n",
       "      <td>123.81</td>\n",
       "      <td>0.01</td>\n",
       "      <td>NGVD29</td>\n",
       "      <td>3180005.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>USGS</td>\n",
       "      <td>2492000</td>\n",
       "      <td>Bogue Chitto River near Bush, LA</td>\n",
       "      <td>ST</td>\n",
       "      <td>30.629356</td>\n",
       "      <td>-89.897296</td>\n",
       "      <td>U</td>\n",
       "      <td>NAD83</td>\n",
       "      <td>43.31</td>\n",
       "      <td>0.22</td>\n",
       "      <td>NAVD88</td>\n",
       "      <td>3180005.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  agency_cd  site_no                                       station_nm  \\\n",
       "0      USGS  2489500                    Pearl River near Bogalusa, LA   \n",
       "1      USGS  2490193     Pearl River at Pools Bluff near Bogalusa, LA   \n",
       "2      USGS  2490200  Pearl River Navigation Canal Lock and Dam No. 3   \n",
       "3      USGS  2491500            Bogue Chitto River at Franklinton, LA   \n",
       "4      USGS  2492000                 Bogue Chitto River near Bush, LA   \n",
       "\n",
       "  site_tp_cd  dec_lat_va  dec_long_va coord_acy_cd dec_coord_datum_cd  alt_va  \\\n",
       "0         ST   30.793243   -89.820907            U              NAD83   54.64   \n",
       "1         ST   30.705833   -89.846389            5              NAD83    0.32   \n",
       "2         ST   30.638056   -89.868056            5              NAD83   33.98   \n",
       "3         ST   30.842962   -90.162030            U              NAD83  123.81   \n",
       "4         ST   30.629356   -89.897296            U              NAD83   43.31   \n",
       "\n",
       "   alt_acy_va alt_datum_cd     huc_cd  \n",
       "0        1.00       NAVD88  3180004.0  \n",
       "1        0.16       NAVD88  3180004.0  \n",
       "2        0.20       NAVD88  3180004.0  \n",
       "3        0.01       NGVD29  3180005.0  \n",
       "4        0.22       NAVD88  3180005.0  "
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "f.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "stations = []\n",
    "site_dict={}\n",
    "\n",
    "def station_on_click_handler(event=None, id=None, properties=None, **args):\n",
    "    plot_station_daily(site_dict[tuple(args['coordinates'])])\n",
    "\n",
    "for site_no, station_nm, lat, long in zip(f.site_no,f.station_nm,f.dec_lat_va, f.dec_long_va):\n",
    "\n",
    "    stations.append(Marker(location=(float(lat),float(long)), \n",
    "                           draggable=False,\n",
    "                           title=\"Site: {0}\\nStation: {1}\".format(site_no,station_nm)))\n",
    "    site_dict[tuple(stations[-1].location)] = site_no\n",
    "    stations[-1].on_click(station_on_click_handler)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "marker_cluster = MarkerCluster(\n",
    "    markers=tuple(stations)\n",
    ")\n",
    "\n",
    "m.add_layer(marker_cluster);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.8.1"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": false,
   "sideBar": false,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": false,
   "toc_window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}