00-habitat-suitability-grasslands.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Import necessary packages\n",
    "import os\n",
    "import pathlib\n",
    "\n",
    "import geopandas as gpd # Working with vector data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'C:\\\\Users\\\\riede\\\\earth-analytics\\\\data\\\\habitat-suitability'"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Set up a project data directory\n",
    "data_dir = os.path.join(\n",
    "    pathlib.Path.home(),\n",
    "    'earth-analytics',\n",
    "    'data',\n",
    "    'habitat-suitability'\n",
    ")\n",
    "\n",
    "os.makedirs(data_dir, exist_ok=True)\n",
    "\n",
    "data_dir"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Define study area:\n",
    "\n",
    "* Download the USFS National Grassland Units\n",
    "\n",
    "* Select study sites: Caddo National Grasslands and Sheyenne National Grasslands\n",
    "    * Caddo National Grasslands is mainly in Texas, USA, although a small portion of it is in Oklahoma too. It is comprised of two units: Bois d' Arc Creek Unit (13,360 acres) and Ladonia Unit (2,780 acres). It is 17,785 acres and has 3 lakes: Lake Coffee Mill (651 acres), Lake Crockett (388 acres), and Lake Fannin (45 acres). In the Caddo Grasslands you can do several recreational activities such as camping, picnicking, and fishing. The Caddo is also home to wildlife such as coyotes, waterfowl, turkey, largemouth bass, and sunfish. \n",
    "    * Sheyenne National Grasslands is in the southeastern corner of North Dakota, USA and includes 70,180 acres of public land amid 64,769 acres of private land. Similar to Caddo National Grasslands, you can camp, hike, hunt, and backpack as well as other recreational activities. This tallgrass prairie is home to prairie chickens, the Dakota skipper (a small butterfly), the Regal Fritillary (another type of butterfly), ferns, and the western prairie fringed orchid. \n",
    "\n",
    "* Citations:\n",
    "    1. “Caddo-LBJ National Grasslands.” National Forests and Grasslands in Texas, U.S. Forest Service, U.S. Department of Agriculture, www.fs.usda.gov/detail/texas/about-forest/districts/?cid=fswdev3_008440. Accessed 2 Dec. 2024. \n",
    "    2. Caddo National Grasslands WMA, Texas Parks and Wildlife Department, tpwd.texas.gov/huntwild/hunt/wma/find_a_wma/list/?id=4. Accessed 2 Dec. 2024.\n",
    "    3. “Sheyenne National Grassland.” Dakota Prairie Grasslands, U.S. Forest Service, U.S. Department of Agriculture, www.fs.usda.gov/recarea/dpg/recarea/?recid=79470. Accessed 2 Dec. 2024. \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "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>OBJECTID</th>\n",
       "      <th>NATIONALGR</th>\n",
       "      <th>GRASSLANDN</th>\n",
       "      <th>GIS_ACRES</th>\n",
       "      <th>SHAPE.AREA</th>\n",
       "      <th>SHAPE.LEN</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>187245</td>\n",
       "      <td>281771010328</td>\n",
       "      <td>Fort Pierre National Grassland</td>\n",
       "      <td>209044.225</td>\n",
       "      <td>0.095149</td>\n",
       "      <td>1.455518</td>\n",
       "      <td>POLYGON ((-100.0841 44.28163, -100.0841 44.278...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>187246</td>\n",
       "      <td>295507010328</td>\n",
       "      <td>Butte Valley National Grassland</td>\n",
       "      <td>19489.170</td>\n",
       "      <td>0.008557</td>\n",
       "      <td>0.853736</td>\n",
       "      <td>MULTIPOLYGON (((-121.99601 41.84049, -121.9960...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>187247</td>\n",
       "      <td>295508010328</td>\n",
       "      <td>Kiowa National Grassland</td>\n",
       "      <td>144281.321</td>\n",
       "      <td>0.058543</td>\n",
       "      <td>9.858642</td>\n",
       "      <td>MULTIPOLYGON (((-104.30415 36.08063, -104.3041...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>187248</td>\n",
       "      <td>295509010328</td>\n",
       "      <td>Sheyenne National Grassland</td>\n",
       "      <td>70428.175</td>\n",
       "      <td>0.033356</td>\n",
       "      <td>4.097398</td>\n",
       "      <td>MULTIPOLYGON (((-97.31082 46.51457, -97.3056 4...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>187249</td>\n",
       "      <td>295510010328</td>\n",
       "      <td>Cedar River National Grassland</td>\n",
       "      <td>6717.517</td>\n",
       "      <td>0.003157</td>\n",
       "      <td>0.999947</td>\n",
       "      <td>MULTIPOLYGON (((-101.82222 45.95897, -101.8170...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>187250</td>\n",
       "      <td>295511010328</td>\n",
       "      <td>Black Kettle National Grassland</td>\n",
       "      <td>33103.349</td>\n",
       "      <td>0.013340</td>\n",
       "      <td>4.561060</td>\n",
       "      <td>MULTIPOLYGON (((-99.9166 35.71892, -99.91662 3...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>187251</td>\n",
       "      <td>295512010328</td>\n",
       "      <td>Rita Blanca National Grassland</td>\n",
       "      <td>94127.091</td>\n",
       "      <td>0.038271</td>\n",
       "      <td>6.652113</td>\n",
       "      <td>MULTIPOLYGON (((-102.62994 36.44073, -102.6298...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>187252</td>\n",
       "      <td>295513010328</td>\n",
       "      <td>Thunder Basin National Grassland</td>\n",
       "      <td>626249.208</td>\n",
       "      <td>0.282888</td>\n",
       "      <td>44.088050</td>\n",
       "      <td>MULTIPOLYGON (((-105.46006 43.31909, -105.4601...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>187253</td>\n",
       "      <td>295514010328</td>\n",
       "      <td>McClellan Creek National Grassland</td>\n",
       "      <td>1401.715</td>\n",
       "      <td>0.000562</td>\n",
       "      <td>0.115902</td>\n",
       "      <td>POLYGON ((-100.86004 35.20951, -100.86009 35.2...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>187254</td>\n",
       "      <td>295515010328</td>\n",
       "      <td>Caddo National Grassland</td>\n",
       "      <td>68479.549</td>\n",
       "      <td>0.026940</td>\n",
       "      <td>1.159342</td>\n",
       "      <td>MULTIPOLYGON (((-95.85493 33.79814, -95.85495 ...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>187255</td>\n",
       "      <td>295516010328</td>\n",
       "      <td>Lyndon B. Johnson National Grassland</td>\n",
       "      <td>115408.461</td>\n",
       "      <td>0.045238</td>\n",
       "      <td>1.177130</td>\n",
       "      <td>POLYGON ((-97.67965 33.43397, -97.67952 33.433...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>187256</td>\n",
       "      <td>295517010328</td>\n",
       "      <td>Grand River National Grassland</td>\n",
       "      <td>154665.301</td>\n",
       "      <td>0.072336</td>\n",
       "      <td>9.796405</td>\n",
       "      <td>MULTIPOLYGON (((-102.08143 45.68619, -102.0814...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>187257</td>\n",
       "      <td>295518010328</td>\n",
       "      <td>Buffalo Gap National Grassland</td>\n",
       "      <td>654877.177</td>\n",
       "      <td>0.295196</td>\n",
       "      <td>33.006274</td>\n",
       "      <td>MULTIPOLYGON (((-102.3384 43.93679, -102.3384 ...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>187258</td>\n",
       "      <td>295519010328</td>\n",
       "      <td>Crooked River National Grassland</td>\n",
       "      <td>173593.026</td>\n",
       "      <td>0.079499</td>\n",
       "      <td>2.815873</td>\n",
       "      <td>POLYGON ((-120.98026 44.72935, -120.98032 44.7...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>187259</td>\n",
       "      <td>295520010328</td>\n",
       "      <td>Little Missouri National Grassland</td>\n",
       "      <td>1025317.517</td>\n",
       "      <td>0.492748</td>\n",
       "      <td>60.421675</td>\n",
       "      <td>MULTIPOLYGON (((-104.01689 47.51794, -104.0115...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>187260</td>\n",
       "      <td>295521010328</td>\n",
       "      <td>Oglala National Grassland</td>\n",
       "      <td>215804.927</td>\n",
       "      <td>0.096279</td>\n",
       "      <td>1.970612</td>\n",
       "      <td>POLYGON ((-103.72478 43.001, -103.72008 43.000...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>187261</td>\n",
       "      <td>295522010328</td>\n",
       "      <td>Comanche National Grassland</td>\n",
       "      <td>444413.904</td>\n",
       "      <td>0.183064</td>\n",
       "      <td>26.658022</td>\n",
       "      <td>MULTIPOLYGON (((-104.02264 37.69225, -104.0225...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>187262</td>\n",
       "      <td>295523010328</td>\n",
       "      <td>Pawnee National Grassland</td>\n",
       "      <td>208424.885</td>\n",
       "      <td>0.089972</td>\n",
       "      <td>15.341594</td>\n",
       "      <td>MULTIPOLYGON (((-104.58107 40.82665, -104.5810...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>187263</td>\n",
       "      <td>295524010328</td>\n",
       "      <td>Cimarron National Grassland</td>\n",
       "      <td>109101.348</td>\n",
       "      <td>0.044765</td>\n",
       "      <td>5.539623</td>\n",
       "      <td>MULTIPOLYGON (((-101.98978 37.11934, -101.9852...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>187264</td>\n",
       "      <td>295525010328</td>\n",
       "      <td>Curlew National Grassland</td>\n",
       "      <td>74783.572</td>\n",
       "      <td>0.032983</td>\n",
       "      <td>1.771642</td>\n",
       "      <td>MULTIPOLYGON (((-112.69951 42.24751, -112.6995...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    OBJECTID    NATIONALGR                            GRASSLANDN    GIS_ACRES  \\\n",
       "0     187245  281771010328        Fort Pierre National Grassland   209044.225   \n",
       "1     187246  295507010328       Butte Valley National Grassland    19489.170   \n",
       "2     187247  295508010328              Kiowa National Grassland   144281.321   \n",
       "3     187248  295509010328           Sheyenne National Grassland    70428.175   \n",
       "4     187249  295510010328        Cedar River National Grassland     6717.517   \n",
       "5     187250  295511010328       Black Kettle National Grassland    33103.349   \n",
       "6     187251  295512010328        Rita Blanca National Grassland    94127.091   \n",
       "7     187252  295513010328      Thunder Basin National Grassland   626249.208   \n",
       "8     187253  295514010328    McClellan Creek National Grassland     1401.715   \n",
       "9     187254  295515010328              Caddo National Grassland    68479.549   \n",
       "10    187255  295516010328  Lyndon B. Johnson National Grassland   115408.461   \n",
       "11    187256  295517010328        Grand River National Grassland   154665.301   \n",
       "12    187257  295518010328        Buffalo Gap National Grassland   654877.177   \n",
       "13    187258  295519010328      Crooked River National Grassland   173593.026   \n",
       "14    187259  295520010328    Little Missouri National Grassland  1025317.517   \n",
       "15    187260  295521010328             Oglala National Grassland   215804.927   \n",
       "16    187261  295522010328           Comanche National Grassland   444413.904   \n",
       "17    187262  295523010328             Pawnee National Grassland   208424.885   \n",
       "18    187263  295524010328           Cimarron National Grassland   109101.348   \n",
       "19    187264  295525010328             Curlew National Grassland    74783.572   \n",
       "\n",
       "    SHAPE.AREA  SHAPE.LEN                                           geometry  \n",
       "0     0.095149   1.455518  POLYGON ((-100.0841 44.28163, -100.0841 44.278...  \n",
       "1     0.008557   0.853736  MULTIPOLYGON (((-121.99601 41.84049, -121.9960...  \n",
       "2     0.058543   9.858642  MULTIPOLYGON (((-104.30415 36.08063, -104.3041...  \n",
       "3     0.033356   4.097398  MULTIPOLYGON (((-97.31082 46.51457, -97.3056 4...  \n",
       "4     0.003157   0.999947  MULTIPOLYGON (((-101.82222 45.95897, -101.8170...  \n",
       "5     0.013340   4.561060  MULTIPOLYGON (((-99.9166 35.71892, -99.91662 3...  \n",
       "6     0.038271   6.652113  MULTIPOLYGON (((-102.62994 36.44073, -102.6298...  \n",
       "7     0.282888  44.088050  MULTIPOLYGON (((-105.46006 43.31909, -105.4601...  \n",
       "8     0.000562   0.115902  POLYGON ((-100.86004 35.20951, -100.86009 35.2...  \n",
       "9     0.026940   1.159342  MULTIPOLYGON (((-95.85493 33.79814, -95.85495 ...  \n",
       "10    0.045238   1.177130  POLYGON ((-97.67965 33.43397, -97.67952 33.433...  \n",
       "11    0.072336   9.796405  MULTIPOLYGON (((-102.08143 45.68619, -102.0814...  \n",
       "12    0.295196  33.006274  MULTIPOLYGON (((-102.3384 43.93679, -102.3384 ...  \n",
       "13    0.079499   2.815873  POLYGON ((-120.98026 44.72935, -120.98032 44.7...  \n",
       "14    0.492748  60.421675  MULTIPOLYGON (((-104.01689 47.51794, -104.0115...  \n",
       "15    0.096279   1.970612  POLYGON ((-103.72478 43.001, -103.72008 43.000...  \n",
       "16    0.183064  26.658022  MULTIPOLYGON (((-104.02264 37.69225, -104.0225...  \n",
       "17    0.089972  15.341594  MULTIPOLYGON (((-104.58107 40.82665, -104.5810...  \n",
       "18    0.044765   5.539623  MULTIPOLYGON (((-101.98978 37.11934, -101.9852...  \n",
       "19    0.032983   1.771642  MULTIPOLYGON (((-112.69951 42.24751, -112.6995...  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "<Axes: >"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGdCAYAAADAAnMpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABECUlEQVR4nO3deXxU5d3//9fsM0lmJntCFsK+G1FUwKKIgGi5lWIXrZZa26qtcru13t7aVVtFa+vP7Ve91S5a2tLF4o21UrQC1tstEDZZBNmyJ5BlJtvs5/sHNZqyZuMk4f18PObxYM45c85nrp4671znOtexGIZhICIiImISq9kFiIiIyKlNYURERERMpTAiIiIiplIYEREREVMpjIiIiIipFEZERETEVAojIiIiYiqFERERETGV3ewC/l0ikaCqqgqv14vFYjG7HBERETkKwzBobm4mLy8Pq7X7/Rv9LoxUVVVRWFhodhkiIiJygsrLyykoKOj25/tdGPF6vcChL+bz+UyuRkRERI4mGAxSWFjY8dvdXf0ujHx0acbn8ymMiIiIDAA9HVahAawiIiJiKoURERERMZXCiIiIiJhKYURERERMpTAiIiIiplIYEREREVMpjIiIiIipFEZERETEVAojIiIiYiqFERERETGVwoiIiIiYSmFERERETNWjMLJkyRIsFgu33nprx7KWlhYWL15MQUEBHo+H8ePH8+STT/a0ThERERmkuv3U3pKSEp5++mmKi4s7Lb/ttttYvXo1S5cuZdiwYaxatYobb7yRvLw8FixY0OOCRUSkd6zb18CwzGRawzGi8QT1LRGmFKVht6nTXE6ubp1xLS0tXH311TzzzDOkpaV1Wvf2229zzTXXcMEFFzBs2DCuv/56Tj/9dNatW9crBYuISM+1RWLc+cJm7n1pK/e+tI2/b63ljV0HePqfe3h5c7XZ5ckpplth5KabbmL+/PnMmTPnsHUzZsxgxYoVVFZWYhgGq1evZufOncybN++I+wqHwwSDwU4vERHpW+FYgqvOKeDsYen898VjuWnWKO6YN47PTylkdE4yJfsazC5RTiFdvkyzbNkySktLKSkpOeL6xx57jOuuu46CggLsdjtWq5Vnn32WGTNmHHH7JUuWcM8993S1DBER6YH2cIz2iMEV5xSQ4vr4pyDL6yLL6zKxMjkVdalnpLy8nFtuuYWlS5fidruPuM1jjz3GO++8w4oVK1i/fj0/+9nPuPHGG3nttdeOuP1dd91FIBDoeJWXl3f9W4iISJfUNoeoaGqnNhjqWNYSjplYkZzKLIZhGCe68YsvvsjChQux2Wwdy+LxOBaLBavVSiAQIC0tjeXLlzN//vyObb7+9a9TUVHBypUrj3uMYDCI3+8nEAjg8/m6+HVERORExeMJbP8arBqJJXjgb9v41OgsZo/PMbkyGSh66ze7S5dpZs+ezZYtWzotu/baaxk3bhx33nkn8XicaDSK1dq5w8Vms5FIJLpdpIiIHFkiYVDW0EZFYzv1rWFmj8sixe08oc/aPnHXjNNu5fyx2SRO/O9TkV7TpTDi9XqZNGlSp2XJyclkZGR0LJ85cyZ33HEHHo+HoqIi1q5dy/PPP8/DDz/ce1WLiAjReIJfvrmXX7y5l1jC4NsXjaE1kiDlyFfRj+uCsdm9W6DICer2PCNHs2zZMu666y6uvvpqGhoaKCoq4r777uMb3/hGbx9KROSU9vQbe3h9Rx3NoRhWC/g8DtKSTqxXJBKN8+GBVpKcVoZlpvRxpSLH1uMwsmbNmk7vc3Nz+dWvftXT3YqIyHFsKm+iJRRldE4KH9Q0U5SejNN+YvcllJY18t3/3crIzGS+cHYhmSkuTi9M7duCRY6i13tGRETk5Mj1u3HYrVQ1teO0W4nE48ChcSQAVqsFOHSXTE0ghGEYjM7xsu9gKz9fs5vJhamMzfGyoayJ9GQno7JTSHbpZ0FOPp11IiID1D+212GxHLoTZvqIDH788nbmnzaEsoY2tlUF+cFlExib4+P17bXkpXrYUNbEsIxkhmUm86lRmSyaVkSSwof0AzoLRUQGqAvHZbOxookvnTOUJ9fu5mBLhJ/X76ahNcK4XC9L3y5jWGYy10wvImYYpCY5eXdvPRWNbVQ1tdISjiqMSL+gs1BEZID6Q0k5F03M4XfvlXHPgomUNbTzzu56tlQG8LrtTCny0RJOsPdgCy+UVpKe7OSNnQcB2FjexKLpI8jWdE7SD3Rp0rOTQZOeiYgcn2EY/PLNvUwpSuOX/7ePffWtfGZyHoYBGSkuyutb2VffxrghXn791n6+cFYh8yblMCbby/INlTSHolx5zlDcDtvxDyZyFKZMeiYiIv2DYYABbKtuZlhGEhnJTnbUNPPHdRX84D/G8n5lI80Rg3Vljdz3mUnc+ZfNNIejXH/eCD47pcDs8kU6URgRERmg3tvbQDxhEIknuP78EWyrDrLwjHy8bhc//cIZLCspJ9Pr4kcvb+PGC0by2vY60pKcbCpvpKk9xrkjM3DYuvXwdpFepTAiIjIAxRIGTpuV0XleXttey+/eLWPvwVaC7VHOH53J7gMtrNl5gIrGdjL+NVZk/mm5fH/F+4zMPDQfycQhPjL1hF7pBzRmRERkAIonDCob29h9sBW/x8Gu2mYsFgvJTjt5fhf3v7KDkn2N5PndNLZFyfG5uGnWKLJ9brKSnSS5bOSnJalnRHpEY0ZERE5hNquFoRnJDM1IBuDMoWkd69bsqGV4etKhMJLq4WBLhKnDM3h9ex0zx2Zx9rA0nDYrdgUR6ScURkREBokdNUGWb6ikNRRlVLaXc4als+dgC+PzfFQHQ7yx8wA//XwxlU1tjNY9vdKPKIyIiAwS/9x5kDU7DmAYBhPy/RxsDTM620tZQxtzx2czbXg6pWVNnFGUdvydiZxECiMiIoPE3Ak5eN12GtoixOIGC07P49k39zI6O4Vo3GDuhCya2mJmlylyGIUREZFBYljmoefOrNhYxaptNRRlJJGW5CTJaWdCno/xuT5sGici/ZDOShGRQWTx70rZWN5EeUMb43K9fHl6EdleF5vKm2hoi5pdnsgRqWdERGQQCUUTbK5oYkRWCsUFqRRlJBNoi/BCaQX/3HmAyzX7qvRD6hkRERlEzhmeRrbPxcQ8H0nOj/7etFDRGCIcS/Dgyh2m1idyJOoZEREZRBIJGOJz897eBvJS3YzN9jIiO4XvXzqBnTXNhKJxs0sUOYzCiIjIIFLe2MZv3y3j5gtHkZHi5O/bagisjzEsM5mzh6Vx6eQ8s0sUOYzCiIjIIDJ1RAYHmsOcPSydIX4PKzZtw2m3Uh1oZ+rwdNKSnGaXKHIYjRkRERlEslKc3DJnNGkpDtbvb8RmteC0WXHabfx01Qc0tETMLlHkMAojIiKDyIcHWtlWFaS2KUTp/iYisQQ2q4XSskZmjc2moTVsdokih9FlGhGRQeTS04ZQ1tjKztoWRuck87ct1SyYnMekfD9pSU6GZSWbXaLIYRRGREQGkXv/uo1wLE5VIERZfRtfmFJAZVM726ub+cqnhuHUDKzSDymMiIgMQImEwUubqyhM91BW38bBlgjDM5OZPT6bSXl+0pIdvLKliupAiGhLgukjMyhKT8Iw4L199dS3hDl7WAZZXpfZX0VEY0ZERAaihGHQ3B5lV20L2Skusn1uYvEEI7KSsVigorGNcUP8VDaF+eyUAuw2uO0PG/njujJGZKbgdzvYVN5k9tcQIByLs2prjdllmMpiGIZhdhGfFAwG8fv9BAIBfD6f2eWIiPRLiUSCHTVBXtpcw6byJs4ZnsZnzsjH53bw+/fKSSQMJuR5qQ1GOC3PT0G6B6sVwILf46CxNUJasm7zNdsrW6r51KhM/riuDIwELoeDz59ViNthM7u0E9Jbv9m6TCMiMgAZBrRH4iw4fQiRaIx/7jxIRoqTYFuMgvQk5ozPIdllJxKNY7dZsVotnT6vIGK+vQdaWLWtlle31fLa9lpunj2aS0/PGzBBpDfpMo2IyABktVqwWizUBMOMyErmm7NGEmyPc9W0IkZmpfDIazt5Y9cB6tsivLixkkgsYXbJ8m+iCYPXttWS5XVw2+xRTMpLIfUUnZROYUREZIAakuohGovT1BrBarFwzfQi0pKcvPnhQb46YzjtkThVTSFGZ6cQ1jNp+p2GljD+JAdbKoOMzEqmrCFMyb4Gs8syhS7TiIgMUAdbIgzLTKGpLcLSd8o40Bxm9vhsvjStCLvVwuRCP1Pvf52N35+L6xTs+u/vpo3M5Jkvn4XHYePxf+wiFEtQFWjn7GHpZpd20mkAq4jIAPX27npSPXYONodIT3GR4/MQisZoj8XZVB5kcuGhic4yUnT7bn8XjSewWizY/m1sT3+nAawiIqe46SMz/vUv/yeWuojFE+SnJmG1WE7JwZADkeMUn4xOYUREZJCx26zYT/EfNxlYenS2LlmyBIvFwq233tpp+fbt27nsssvw+/14vV6mTZtGWVlZTw4lIiIig1S3w0hJSQlPP/00xcXFnZbv3r2bGTNmMG7cONasWcOmTZv43ve+h9vt7nGxIiIip4rNFU1sqWgyu4yToluXaVpaWrj66qt55pln+PGPf9xp3Xe+8x0+/elP85Of/KRj2YgRI3pWpYiIDEjNoShOuxWX/eOxK5FYAqddl5GOpyjj0NT+p4JunQ033XQT8+fPZ86cOZ2WJxIJXn75ZcaMGcO8efPIzs5m6tSpvPjii0fdVzgcJhgMdnqJiMjgkOy0dzwpOBZPEI0nsA+wO0bM4vc48LkdZpdxUnQ5jCxbtozS0lKWLFly2Lq6ujpaWlp44IEHuPjii1m1ahULFy7k8ssvZ+3atUfc35IlS/D7/R2vwsLCrn8LERHpN/YcaKG8sQWg0zT0dpsVx7+mpm+PxDUrrHTo0mWa8vJybrnlFlatWnXEMSCJxKETa8GCBdx2220ATJ48mbfeeounnnqKmTNnHvaZu+66i9tvv73jfTAYVCARETFBPGHw9601JDttTC5MpaEtSigaZ/yQE58/4kBziEg8QZLT0XE55mBLmCxv598Mj1O3HMvHuhRG1q9fT11dHVOmTOlYFo/HeeONN3jiiSdobW3FbrczYcKETp8bP348b7755hH36XK5cLk0IY+IiNliiQRrPqjj6qlFlDe2s7++jTd2HmD2+GzSU5wM8blpDscYl3v0cJLldR8WPP79vci/61IYmT17Nlu2bOm07Nprr2XcuHHceeeduFwuzj77bD744INO2+zcuZOioqKeVysiIn3GarEwOtvLik1VXDwxh7G5KdS3hnlr90EmDPHz4oZKhmemHDOMiHRHl8KI1+tl0qRJnZYlJyeTkZHRsfyOO+7giiuu4Pzzz2fWrFmsXLmSl156iTVr1vRa0SIi0vscNivXnX/o7sdEwqC0rJEDzWHOHJpGltdNVaCdy07P69h+z4EW/ndjFd+8YKRmepUe6fUZWBcuXMhTTz3FkiVLuPnmmxk7diwvvPACM2bM6O1DiYhIH7FaLZw1LJ0P61qYMTqL9GTnJ6afPyTJaeeCsVkdd8uIdJcelCciIkcUiyew26ys29dAZVM7E/P8jMpOYf3+Rl4orWBTeRPnjc5iVHYKn5tSYHa5YoLe+s1WnBURkSP6x446tlQEeOz1D8lKcTEqOwWAmkCItR8c4OvnDWdyoZ9ZY7NMrlQGOj0oT0REjmjO+Byu+J+3GZObgs1qIRSNAzC/eAizx2fjsluxWCy0R2L8fWsNU4rSyEzR3ZHSdQojIiJyRDarhYl5PiLxBCu31tAejZOR7OK0An+nAavhWIIn1+zGaoFvzxvLuSMzTaxaBiKFEREROaqrpxUxJsfLO3vqmTYi44jb2G1Wpo5IJ8/v4cyhaSe5QhkMFEZEROSoxuR4AY4aRAAcNgvlDW3sqm3B7bByxdlDT1Z5MkgojIiISI84bVZunTOasvp2po5IN7scGYAURkREpEcsFgtjcnyMydF0DNI9urVXRERETKUwIiIiIqZSGBERERFTKYyIiIiIqRRGRERExFQKIyIiImIqhRERERExlcKIiIiImEphREREREylMCIiIiKmUhgRERERUymMiIiIiKkURkRERMRUemqviEgfC0Xj1LdGsABtkTijslPMLkmkX1EY6Ybyhjb8SQ58bofZpYjIAPCX0gru/es2Jhem4rbbeOCzxeT63WaXJdJv6DJNFyUSBtf88j121baYXYqIDADhaJyVW2u54fyRNIdi3DBzBKlJ+kNG5JMURk5AZVM78YQBQEVjO6NzUpiY5zO5KhEZCOpbI3z2jHx21TUzZWgquT434VjC7LJE+hVdpjkBT675kDE5Xr48fRjf+9/32V/fitOmHCcix7errpkkp5U543NwO6zEEwYeh83sskT6FYWR4/jqr0tobItww/kj2V4dpLyhDZ/HgdVqMbs0ERkAZo7JNrsEkX5Pf94fQyga5+3d9RgGFKYnsb++jT0HW9lWFWRXbbPZ5YmIiAwK6hk5ippAiEufeJP2aJwvnlPIvP/vDeqaQwDEEgbBUMzkCkVERAYH9YwcxQulFRxoDgPw+Osfkul14rJ/fJ03HI2bVZqIiMigop6Ro0h22vjclAKMQzfRMHV4OsFQlJJ9DRgGpCU7zS1QRERkkLAYxkc/t/1DMBjE7/cTCATw+XT7rIiISH/VW7/ZukwjIiIiptJlGhERGTDK6ttY+u5+Xnm/utPyqcMz+OnnTzepKukphRERERkwXtpcxdNv7Dls+YjMsAnVSG/p0WWaJUuWYLFYuPXWW4+4/oYbbsBisfDII4/05DAiIoPartpmPqgJEGgLs706wGvbaghH47RHNIXAiXpr90GicU2zP1B1O4yUlJTw9NNPU1xcfMT1L774Iu+++y55eXndLk5E5FRQsq+BmkCYnXXNOOwWmtqj7KgJcv/fdhDSNAKd7DvYesTlJd+Zg10zYw9Y3QojLS0tXH311TzzzDOkpaUdtr6yspLFixfz29/+FodDT6cUETmWyYV+8tM8pCc5eGljNWcUpvH46x8yc2wW++vbzC6vXxmWmXzYskXTikhNcmKxKIwMVN0KIzfddBPz589nzpw5h61LJBIsWrSIO+64g4kTJx53X+FwmGAw2OklInKqiCcMagNh0pOdZHmTuOS0XPY3tDJ3Qg6//r99/Or/9qp35BjsVgufP6vA7DKkh7o8gHXZsmWUlpZSUlJyxPUPPvggdrudm2+++YT2t2TJEu65556uliEiMihYgLw0D7F4grZIjHG5fsZk+/iwrpkLx2Xzz10HMXs2qGAoSiJhkJr08WSPew604LRbKUhLOqm1TMzzMb94CH/bcuhuGrvNQnFB6kmtQXpfl8JIeXk5t9xyC6tWrcLtdh+2fv369Tz66KOUlpaecHfZXXfdxe23397xPhgMUlhY2JWyREQGLKvVQls0xhC7m5pgiETCQkGah/zUJJLddi4/s+//6m+PxGkJR0lPdmE7wriL8oY2/v5+DdNGZFATDOF12bn35W18b/4ESvY1MCbHy9u767FY4D9OG8Lft9UyOttLbTDEZ87I79VaLxibzQVj9STkwaZLM7C++OKLLFy4EJvt42e0xONxLBYLVquVBx98kDvuuAOr1dppvdVqpbCwkH379h33GJqBVURONYmEQTyeoDLQSm1zhDML09hzsIUhqUn43H0/7u6qZ95h/f5GPjulgC0VTcwvzuOy0/PIS/V0bLOxvInW8KG7e/66uZq/lFZQXODn/cogXz9vOK9uq+XL04oIxxKMyfGysaKJ1CQHV08t6vP6xTy99ZvdpZ6R2bNns2XLlk7Lrr32WsaNG8edd97JkCFDmDdvXqf18+bNY9GiRVx77bXdLlJEZDBqbouwtTrIWcPSiQPxhJUJQ/z8pbSSAy1hFl84mmfe2MOuumZ+8rm+m9DrazOGc97oLF55v5qm9ihv7a7nq58a1mmbyYWpHf8+e1g600akE4rG+bBuB3/fWkO2180/PzzIxRNzWLGpkrrmMMMzU/qsZhlcuhRGvF4vkyZN6rQsOTmZjIyMjuUZGRmd1jscDnJzcxk7dmwPSxURGbjqmtv566ZqUpMcFOf7sFisNLVHCMUSWC0WyhpaqQ2EqGhq40BLmHmTctlW3UQskSAcS7DvYOsR7yTpDbPH5zB7fA6fnZLP794t4+yiNOqCIUKxBB/WtXB6YSpD/B/3kjjtVmaMyiQjxYXf4yA/NYlN5Y24nTbWlzWxo6aZEVkp7Kxt5tl/7sHvcXDJpFxSTkIvjwxMmoFVROQkSHU7GZeTQkskQTwBT639kG3VzQBMW5zByKwUrFgItEcYlZWM3+3g1W0NfPq0IXx2SgHZ3sPH6fWW2mCI3XUtJAz4zwtH88zaXby+o464AQdbwrSG41w4LptYwuCP68pZsbGKhrYIX5sxnMvPzCeeMLhq6lDe2VPPt/+0mTy/m9MLU8nxuTjQEmZUdgpOux6FJkfX4zCyZs2aY64/kXEiIiKDRSQaJ2EYuJ126oIhDMOgOhBiRFYyDruNPRVBhmUmkeVz4zzQQnMoxvbqIM3tUUbnesHioDDdw4rNVUwp9FMTbGdsTt+Nn/vhiq1MHZ7OlKI02qNx9tY18/w7ZVQFwpxVlMYXzxnKa9treXlLFdedN4KfrfqA1CQn7ZE4m8qb2FEdZOGZBURiCVZtq+N788fREorym3creO7as0hPcZLicuC0245fjJyyujSA9WTQAFYRGQxW76ilJhjCgoUcn5v9DS1YgAPNUf5SWsHcCTkMzUgiyWln6Tv7+dLUoUzI83FafirWf93RUtHQSun+Js4ant5pMGlveb8ywA2/WY/FAuNyfaQlOThneDp7DrYyxO8m13eoh+PXb+3j2nOH0dAW4YOaZnbUNDNjVCbnDE+nvKGNEVkpROMJNpY3kZXiJN3joCkUY4jPjcOhEDKYmTKAVUREjs0wDAwDhqZ7mJTnZ/mmSrwROxnJLp5/ez/jclJYfOEodh9oJT/Vw9hcLxdPzGX5hgq8bjsVje3k+d3Y7VacDisHWsNke119Uuu2qiAOm4V5E3O5Zc5oWkIxsn2HXw668+JxAGT73IzL9bHgE+tGZB0apOqwWTl7WHrHcl9y39Qsg5PCiIhIL4nFExgGNLSGefbNfUwdkU4kmuDP6yuIJQyicYNsn4d99a0caA7hdtioCYTYe6CV+cW5HGyJUrq/nrkTssjxJYFhUJzvIxJPYLf17piLmkA72T4HN80czmVnFOJy2Ehy6idBzKEzT0SkFyQSBgeaw2QkOTnY2saPLptIXUuYqsYqZo7JYktlkKwUA3+Sg5Vba6gJhAALbZEY10wvwmm1MTEviYl5ftrCMT480MxLG6soLkylqqmdUdneXq33g9pmDgYj/MfpBbh0KUVMpjAiItJLcv1umtsiVDWGSXG14XE48Dht7K9vwzAMtlUHmTcxh/31bTjtVsrqW/jhggmMyPBh+UTHx5bKAEWZSYzN9TN+iJ/8PhgvMqUoHafNSjim596I+XSvlYhIL7BYoKqpjeZIjFljc/E4HGytDpKe7GRDeRMHmkOEYwl++145WV4XbZEY+enJJDmd+JIceD8xB8fUERlkp7iZNTaLIUcYw9EbUlx2XtpUybf+uKlP9i/SFeoZERHpBRaLhapAmNZwjO01LbjtNjAMMpJdGMD2mmZy/W4Mw+Cqc4ayv6GVC8ZmMSbnyJdfrFYLSa5j/ye6oTWM1+3A0c3xJE+s3s0RHkUjctIpjIiI9JLfvrOfZJedRdOGUpCezO66Fh5+dSe761q4ZFIudc1hNpQ10hyKsb++lT0HWhmb46MwvetPvq1sauf+l7cxKc/HtTNG4O7GuI+DLX13p45IVyiMiIj0khFZKXx2SgFuu5XG1jB/3VLFnAk5TB+RzmWT86lobOO9vQ28XxXkplkjSfXYuxVE4NBllvLGdmqCYWZPyD1qD8uRVAfaefS1XZw/Josd1cFuHV+kN2nMiIhIL0lxWVm9o46dtc14nHa+PL2IsTkpzJ2QQ20wRHUgxMIz8pkxKoO/balhzc6DJBLdm3eyviVMjs/N3oOt/HVTVZdCxRC/h3NHZZKiW3mln1AYERHpJV88aygFqW7SUxz4nHZ21rTw7T9tZveBVnYfaGXVtlpe3lLNn9ZVsPqDOl7fcYBgKNKtY9msFnbUBCnKSOKiiTnsPtBywp+NxBK8vbueHbXNDPH33TNvRE6UYrGISC/xuB1MGZbecWdMYVoSj145mdPy/RjA5WfkU98aIcVlJyPFxbTh6di6OYLUarHw1U8Nx261UN8aZdrw9MO2if1rsrRYPAHQMXFaSzjGf108lide38X6/U3dOr5Ib1IYERHpRZ+8RXdM7uHjOLK8LrK8LrwuOxgGKW5nt47jtts4d1QGja1RdlQHyU/1kPGJwwXao4SjcTJTXJQ1tOH3OMhIOTRY9e9bazgQDJHktHPxxNxuHV+kNymMiIicZHMn5JJIGB0PxOuOQCjCs2v3gAX+sqGKRdNa+dLUodQEI0zK8+J22kly2rBaLQzxuztu/00kDAr8bt7bU4/baWfE8OTe+loi3aYwIiJigp4EkUMsLN9YRSxhMHNMJpefUYDFaiXFZWNHbQs+t5X3q1pIJAw+c0Y+26uDnFaQyo6aIA2tEeaMz6E5EtMMrNIvKIyIiAxAFgvMHpfNZ6cUUB1oZ/nGSsKxBOUNbTz02dP4w7oKVm6t5cyhqVx5zlAc/7rdeOk7+xmemUJZQytN7TH2HWzlyrOLzP46coqzGIbRvfvK+kgwGMTv9xMIBPD5fGaXIyLSL31Y18ych9846nrLvzpePup/efbLZ3HB2CwO/SffAhYLVgskDLo9iFakt36z1TMiIjIA+TwO5hcPwWqxUFbfyqaKQKf1H/2Z+dFfmxaLBav18NkcbMoh0g8ojIiIDEDZXjf//1VnAvDatlocdiuRWILXd9RS1tDG/31Y32n7cCxhRpkiJ0RhRERkgJszIafj33Mn5BAMRbn08TfZX9/WsXz1jjounqTbeKV/UhgRERlkfG4HRRnJRGMJwrEEboeN1GTH8T8oYhKFERGRQej5r55DKBpn3b5GZozONLsckWPSs2lERAYpt8OmICIDgsKIiIiImEphREREREylMCIiIiKmUhgRERERUymMiIiIiKkURkRERMRUCiMiIiJiKoURERERMZXCiEgfu+l3pbyypdrsMkRE+i2FEZE+9qMFk/ivFzZz3fPrqA2G2H2gxeySRET6FYURkT6WnuzE53ZQmOZh/mP/ZO7Da1m9o87sskRE+o0ehZElS5ZgsVi49dZbAYhGo9x5552cdtppJCcnk5eXx5e//GWqqqp6o1aRAeueyyZy0cRczh+dxYLJ+aQlOTAMw+yyRET6hW6HkZKSEp5++mmKi4s7lrW1tVFaWsr3vvc9SktL+ctf/sLOnTu57LLLeqVYkYFqzoQchvjdXDY5j2/MHMnkoWlYLBazyxIR6Rfs3flQS0sLV199Nc888ww//vGPO5b7/X5effXVTts+/vjjnHPOOZSVlTF06NCeVSsygH311yXsPtDK56YU8NPPn252OSIi/Ua3ekZuuukm5s+fz5w5c467bSAQwGKxkJqa2p1DiQwaP/386RQX+Llv4SSzSxER6Ve63DOybNkySktLKSkpOe62oVCI//7v/+aqq67C5/MdcZtwOEw4HO54HwwGu1qSyIBwxtA0ViyeYXYZIiL9Tpd6RsrLy7nllltYunQpbrf7mNtGo1GuvPJKEokEP//5z4+63ZIlS/D7/R2vwsLCrpTUZe9XBnhvb0OfHkNEREROnMXowpD+F198kYULF2Kz2TqWxeNxLBYLVquVcDiMzWYjGo3yhS98gT179vD666+TkZFx1H0eqWeksLCQQCBw1N6U7orEEnxz6XoK05P470vG4XbYjv8hEREROaJgMIjf7+/xb3aXLtPMnj2bLVu2dFp27bXXMm7cOO68885OQWTXrl2sXr36mEEEwOVy4XK5ul55NzjtVgrTk0hNcrDnQCsT8no37IiIiEjXdSmMeL1eJk3qPPguOTmZjIwMJk2aRCwW43Of+xylpaX89a9/JR6PU1NTA0B6ejpOp7P3Ku+mz00pINfnJtN7cgKQiIiIHFu3bu09moqKClasWAHA5MmTO61bvXo1F1xwQW8erlsm5fvNLkFEREQ+ocdhZM2aNR3/HjZsmGaVFBERkS7Rs2lERETEVAojIiIiYiqFERERETGVwoiIiIiYSmFERERETKUwIiIiIqZSGBERERFTKYyIiIiIqRRGRERExFQKIyIiImIqhRERERExlcKIiIiImEphREREREylMCIiIiKmUhgRERERUymMiIiIiKnsZhcgIiJd1xqOEYrF2VIRIBJPkJ/qZlh6Cv/34QEuHJ+D3aa/NWXg0NkqItLPBNsjhKLxY27jdthobI1QmJ5EY2uEXJ+HxrYIQ1LdbK8O0haOnaRqRXpOYUREpB9oi8RoaovwwxVbWVZSzu4DzcfdflS2l0K/h6a2CO/ubSAr2cn4IX4a28JsrGjEMIyTVL1IzyiMiIj0A7vrWrn+N+tx2q0MTU/ijyUVnPbDv/PC+nL+UFLWadv6ljBet4NYPEE0nuCLZw/l06cNwW63UtbQxosbqzhraDoWi8WkbyPSNQojIiL9wPINlSQSBsMyPGyuCFDfGuG0fB9rPjiA3+PotG1Te5Q3dtbxgxVbCUaitMeiRKJxorEEqR4HXzl3GJF4wqRvItJ1GsAqItIPfGf+eNbvr+fpN/Zy3phMvjlzBL9+ez+pHgf5qUmdto0nDM4dkUlGsou9B9sYm5OMBbBYLWCxsLOmhUQCJub7cWggqwwACiMiIv2AzWphcmE6159vIRSNY7NZWbevkZpAiEjcYNwQb0ewsFstbKsJUNkYYkx2Cs+/VcZp+T6SXQ7OGJrKheOzSU92mfyNRE6cIrOISD/htFs5qyid7dXNbCxroj0aY+aYLOqaQ9itH4//yPa6eOafe/nDunI2VwUxgN0H25gyNI2YAenJLla+X23eFxHpIvWMiIj0I1arhRtmjmRbVYAlC4tx2q38bUs1hgEfjUdd/cEBZo3NZn7xENbta+Rgc5gcn5NwPMHanQcYN8RHQVrSsQ8k0o9YjH5271cwGMTv9xMIBPD5fGaXIyIiIkfRW7/ZukwjIiIiplIYEREREVMpjIiIiIipFEZERETEVAojIiIiYiqFERERETGVwoiIiIiYSmFERERETNWjMLJkyRIsFgu33nprxzLDMPjhD39IXl4eHo+HCy64gK1bt/a0ThERERmkuh1GSkpKePrppykuLu60/Cc/+QkPP/wwTzzxBCUlJeTm5jJ37lyam5t7XKyIiIgMPt0KIy0tLVx99dU888wzpKWldSw3DINHHnmE73znO1x++eVMmjSJ5557jra2Nn73u9/1WtEiIiIyeHQrjNx0003Mnz+fOXPmdFq+d+9eampquOiiizqWuVwuZs6cyVtvvdWzSkVERGRQ6vJTe5ctW0ZpaSklJSWHraupqQEgJyen0/KcnBz2799/xP2Fw2HC4XDH+2Aw2NWSREREZADrUs9IeXk5t9xyC0uXLsXtdh91O8tHz7n+F8MwDlv2kSVLluD3+ztehYWFXSlJREREBrguhZH169dTV1fHlClTsNvt2O121q5dy2OPPYbdbu/oEfmoh+QjdXV1h/WWfOSuu+4iEAh0vMrLy7v5VURERGQg6tJlmtmzZ7Nly5ZOy6699lrGjRvHnXfeyYgRI8jNzeXVV1/ljDPOACASibB27VoefPDBI+7T5XLhcrm6Wb6IiIgMdF0KI16vl0mTJnValpycTEZGRsfyW2+9lfvvv5/Ro0czevRo7r//fpKSkrjqqqt6r2oREREZNLo8gPV4/uu//ov29nZuvPFGGhsbmTp1KqtWrcLr9fb2oURERGQQsBiGYZhdxCcFg0H8fj+BQACfz2d2OSIiInIUvfWbrWfTiIiIiKkURkRERMRUCiMiIiJiKoURERERMZXCiIiIiJhKYURERERMpTAiIiIiplIYEREREVMpjIiIiIipFEZERETEVAojIiIiYiqFERERETGVwoiIiIiYSmFERERETKUwIiIiIqZSGBERERFTKYyIiIiIqRRGRERExFQKIyIiImIqhRERERExlcKIiIiImEphREREREylMCIiIiKmUhgRERERUymMiIiIiKkURkRERMRUCiMiIiJiKoURERERMZXCiIiIiJhKYURERERMpTAiIiIiplIYEREREVMpjIiIiIipFEZERETEVAojIiIiYqouhZEnn3yS4uJifD4fPp+P6dOn88orr3Ssb2lpYfHixRQUFODxeBg/fjxPPvlkrxctIiIig4e9KxsXFBTwwAMPMGrUKACee+45FixYwIYNG5g4cSK33XYbq1evZunSpQwbNoxVq1Zx4403kpeXx4IFC/rkC4iIiMjAZjEMw+jJDtLT03nooYf42te+xqRJk7jiiiv43ve+17F+ypQpfPrTn+ZHP/rRCe0vGAzi9/sJBAL4fL6elCYiIiJ9qLd+s7s9ZiQej7Ns2TJaW1uZPn06ADNmzGDFihVUVlZiGAarV69m586dzJs376j7CYfDBIPBTi8RERE5dXTpMg3Ali1bmD59OqFQiJSUFJYvX86ECRMAeOyxx7juuusoKCjAbrdjtVp59tlnmTFjxlH3t2TJEu65557ufwMREREZ0LrcMzJ27Fg2btzIO++8wze/+U2uueYatm3bBhwKI++88w4rVqxg/fr1/OxnP+PGG2/ktddeO+r+7rrrLgKBQMervLy8+99GROQU1RKOHbasORTlnT31bK0KmFCRyInr8ZiROXPmMHLkSB555BH8fj/Lly9n/vz5Heu//vWvU1FRwcqVK09ofxozIiLSNX/fWsP3//d9PjelgOvPG4k/yQHAs//cw5NrdjM6J4XffX0aVqvF5EplsDF9zMhHDMMgHA4TjUaJRqNYrZ13abPZSCQSPT2MiIgcwfuVAd7d08BnJucTjxu8tLmKQFuURMJgW1WQL00byrkjM80uU+SYujRm5O677+aSSy6hsLCQ5uZmli1bxpo1a1i5ciU+n4+ZM2dyxx134PF4KCoqYu3atTz//PM8/PDDfVW/iMgpq7Kxne+8+D5el50vnFVASzjOK+9Xc+bQNN7YdYBowuBgS4Tv/ccE9YpIv9alMFJbW8uiRYuorq7G7/dTXFzMypUrmTt3LgDLli3jrrvu4uqrr6ahoYGioiLuu+8+vvGNb/RJ8SIip6r99a28tKkKgPLGNqYUpbOjJkhBWhJ3Ld9CVoqLq6YWMjHPj9thM7lakWPr8ZiR3qYxIyIix/d+ZYB7/7qNoelJtIRizBybhc1i4S8bKvjiOUNJctqZOyHH7DJlkOs3Y0ZEROTkmzDEx2fPzKeuOczNc0YxJieFF0or8LodrNpWy7L3yviftbvNLlPkhKhnRERkAGsORUly2rF9YkxIXTDEH0rKWXzhKCwWjRWRvtNbv9ldnvRMRET6D6/b0el9cyjKik1VOO1WBREZMBRGREQGkdpgiCvOLjwspIj0ZwojIiKDyKhsr9kliHSZBrCKiAwCdc0h2o4wJbzIQKAwIiIyCCQ77SS5utbZXd7Qxs7aIK2hj0NMLJ6gvKGNVgUbOYl0mUZEZBBIPk4QaQ3HaI3EKG9oY3S2l6a2KM++uYdJeT5KjUZmjskm1+/mob9/wKR8P363nXNGZODSQFg5CRRGREROAVWNbWyuDPC798q59PQhnFnop7jAz7L3yll2/TTsNisHmsMEQ1F8bjtej4M/rjv0FPWLxudQGQgxpSjtqPsPx+IAuOya7VW6TvOMiIgMQturg/x1cxVzxuewtSrIiMxkyhvamTU+C4fFgtVisL2mGYfNjsdpJW4YXPr4//HmnbOwWiy8sqWaBZPz+Nv7NRgGzJuYS67ffcRjfVjTTEN7hHOGZ9DQGiE92dlpfU0ghNdtP27vjQw8vfWbrTAiIjII3fnnzUzI8xFoizBrXA7DMpOob4nQ1BbFabdQGwzTHIqS43OxtTrIzNHZeJw2vG47XreDeMLoNJFaoD1KOBqnqT3KmJzOd+xsrwoAkJfqwWKxkOKy8/etNew52MrU4emcVuBXj8kgpUnPRETkqG6YOQKA4ZnJHWM+Pjn3yIQ8SCQMNpY3EWyP8eLGStbvb+SKswupaGynviXCLbNH4/PYsVgsBNqiVDS1MTIr5fCDWSxkpThpbIvyj+11lOxrYGt1AIfNyqdPG6IgIselMCIiMgiNOFJo+DdWq4Wmtgg7a5tx2KwM8XtoaouSSBj8x+lDePQfO7ni7KH4PHb8SQ7y0zI69ZZ8JJEweOS1Xby9p56aQIhsn5uhaR4ev+pMYnGDgy1hMlNcffE1ZZBQGBEROUUlEgkm5ftx2SHZ5SLb6+TZN/dyx7xxvLevgeGZyby+o45V22rYUNbEPZdNJNvr4pLThnTaz0OrPuDNXQcZnpnMnReP47LJedzz0jZ+vno38ybm4PU4FEbkmBRGREROUbFYAr/LRo4vifZonLgBmysCVAfaeb8ywK/f2sdp+X4ALj8zn6ff2MOPF046bD+RWILLJufxpWlFFKZ6aAlFmTEqkynD/BSlp+jWYDkuhRERkVOUw3FoLEdashN3NI7f4+C680fQHIoyND2JZ798Fi6HldHZXlrDUdbvb2Jk5uGXf647bwRnDk3D5bDgdhz6WSlITz6p30UGNt1NIyIiIt3SW7/Zmg5eRERETKUwIiIiIqZSGBERERFTKYyIiIiIqRRGRERExFS6tVdE5BRSWtbIloompg7PID/Ng91qxePUdO1iLvWMiIicQs4cmobP7aCxLUzp/kberwyYXZKIekZERE41OT434ViC4ZnJFGV8PDlZOBbHZbcRiyew2/S3qpw8CiMiIqeYc0dlsqGskVe31XLxpFxsFgu7D7aQ5LCTkewgHEvgslnJTfPoibtyUiiMiIicgnL8LvJTPWyrCjIiK5lE3MDrteNy2NhSGaSqqY1pIzKw2mBYRgopLofZJcsgpjAiInIKammPk5/mprggDYCijGTCsQQOm4UJeT6mjcwg0BahuimEO+fj3pG2SAy71UpjW5gcn8es8mWQURgRETkFpXrs/PPDQEcYaQnF+NHL2xid7eUf22upbGrny9OLmDkmi1ji0BiSF9ZXcNHEHDwOC067jUTCwGrVE3ml5zRCSUTkFNIcigJQ3tTOrroWYvEEb+0+iNVicOvsMeyvb8XncdASipHssuOyW9laGQTAZrXQ1BalrjnMloomBRHpNeoZERE5hXjdDkr2NbC1MoDLbqU5FGNSnp/2aJyhGUk88Nli4gmDP68vZ/rIDIamJ/P2hwcpa2hlVHYKhelJhGNxzh6WbvZXkUFEYURE5BRT2dhOMBQjFIsRiydIS3bj83w8QNVmtXDF2UM73ie57JTub8LjtBFPJDi9MI1QNG5G6TJIKYyIiAwAq3fUUZSRhNthI9vrAiBhgNPe9avtnzkjH4C1Ow9gOYFLLacXpnJ6YWqnZW6HbvmV3qMxIyIiA8DQjCR2H2jh0sff5KXNVdz6h41sqWzq0T5njskiM8XVOwWK9ECXwsiTTz5JcXExPp8Pn8/H9OnTeeWVVzpts337di677DL8fj9er5dp06ZRVlbWq0WLiJwqYvEE8YRBUXoSWyqD/OeFo/jH9jqu/dQwphRp3IYMDl26TFNQUMADDzzAqFGjAHjuuedYsGABGzZsYOLEiezevZsZM2bwta99jXvuuQe/38/27dtxu919UryIyGBXHQjxzp56RmQmEY8nqGsOE4omsFnVsS2Dh8UwDKMnO0hPT+ehhx7ia1/7GldeeSUOh4Pf/OY33d5fMBjE7/cTCATw+Xw9KU1EZMBLJAzmP/4mt8weycGWMMMyUnh7dz1Wq4UvnFVIYXqS2SXKKay3frO7Ha3j8TjLli2jtbWV6dOnk0gkePnllxkzZgzz5s0jOzubqVOn8uKLLx5zP+FwmGAw2OklIiKHtEXjeF12QlGD88Zkc87wDL40vQivy64gIoNGl8PIli1bSElJweVy8Y1vfIPly5czYcIE6urqaGlp4YEHHuDiiy9m1apVLFy4kMsvv5y1a9cedX9LlizB7/d3vAoLC3v0hUREBpMUl53vXzqBaDyBx2HDabcyxO/h+pkjzS5NpNd0+TJNJBKhrKyMpqYmXnjhBZ599lnWrl1Lamoq+fn5fPGLX+R3v/tdx/aXXXYZycnJ/P73vz/i/sLhMOFwuON9MBiksLBQl2lERP4lnjA42BImK8WlWU+lX+mtyzRdnmfE6XR2DGA966yzKCkp4dFHH+Xxxx/HbrczYcKETtuPHz+eN99886j7c7lcuFy6tUxE5GhsVgs5Pt0IcCTBUBSnzap5Twa4Hk96ZhgG4XAYp9PJ2WefzQcffNBp/c6dOykqKurpYURERA7jczuOv5H0e10KI3fffTeXXHIJhYWFNDc3s2zZMtasWcPKlSsBuOOOO7jiiis4//zzmTVrFitXruSll15izZo1fVG7iIhIJ83tUVojMQwDqgPtnPmvuVjqW8NEowlyUz0mVyhH0qUwUltby6JFi6iursbv91NcXMzKlSuZO3cuAAsXLuSpp55iyZIl3HzzzYwdO5YXXniBGTNm9EnxIiIinxQMRalobOW0/FT2HIxT3tBGKBJjSKqH5vao2eXJUfR4npHepnlGRESku2qb2vmgNsj2mmYumpiD3+MiEouT6/dQ3xImQ9Pf9yrTBrCKiIj0V1HDIC/Vg9WSIM3jJDXZ2bHOYbNgGAYWi+5I6m8URkREZNAoSDs0EdxbuxtIdrdxRrKLukCIcCKO1+VQEOmnFEZERGTQOXdkBvWtYSKxBF6PHZ/Fodt/+zGFERERGXRG5XgZhfdf7/RQwf5O/wuJiIiIqRRGRERExFQKIyIiImIqhRERERExlcKIiIiImEphREREREylMCIiIv1GPGGwqbzJ7DLkJFMYERGRfsNmtXB6YarZZchJpjAiIiIiplIYEREREVMpjIiIiIipFEZERETEVAojIiIiYiqFERERETGVwoiIiIiYSmFERERETKUwIiIiIqZSGBERERFTKYyIiIiIqRRGRERExFQKIyIiImIqhRERERExlcKIiIiImEphREREREylMCIiIiKmUhgRERERUymMiIiIiKkURkRERMRUCiMiIiJiKoURERERMVWXwsiTTz5JcXExPp8Pn8/H9OnTeeWVV4647Q033IDFYuGRRx7pjTpFRERkkOpSGCkoKOCBBx5g3bp1rFu3jgsvvJAFCxawdevWTtu9+OKLvPvuu+Tl5fVqsSIiIjL4dCmMXHrppXz6059mzJgxjBkzhvvuu4+UlBTeeeedjm0qKytZvHgxv/3tb3E4HL1esIiIiAwu9u5+MB6P86c//YnW1lamT58OQCKRYNGiRdxxxx1MnDjxhPYTDocJh8Md74PBYHdLEhERkQGoywNYt2zZQkpKCi6Xi2984xssX76cCRMmAPDggw9it9u5+eabT3h/S5Yswe/3d7wKCwu7WpKIiIgMYF3uGRk7diwbN26kqamJF154gWuuuYa1a9fS3t7Oo48+SmlpKRaL5YT3d9ddd3H77bd3vA8GgwokIiIipxCLYRhGT3YwZ84cRo4cyfjx47n99tuxWj/ubInH41itVgoLC9m3b98J7S8YDOL3+wkEAvh8vp6UJiIiIn2ot36zuz1m5COGYRAOh1m0aBFz5szptG7evHksWrSIa6+9tqeHERERkUGqS2Hk7rvv5pJLLqGwsJDm5maWLVvGmjVrWLlyJRkZGWRkZHTa3uFwkJuby9ixY3u1aBERETmyqqZ2hvjdhGMJ3th5gPrWCNWBEP954Sgctv4512mXwkhtbS2LFi2iuroav99PcXExK1euZO7cuX1Vn4iIiHSB024l0B7l2X/u5YnVH3Ysv/GCkThsJhZ2DF0KI7/4xS+6tPMTHSciIiIivSMzxQVAerKTsTlegqEo00dkYLee+M0lJ1uPB7D2Ng1gFRERGRh66ze7f148EhERkVOGwoiIiIiYSmFERERETKUwIiIiIqbq8aRnIiIicnSt4Rjv7KknnjCwWixYLJCf5mFcbt/epPHHknL+XFrBjxZMYmyut0+P1VMKIyIiIn3IZrVw3fPrSHzi3lWPw8bym87t00DSFomRluTAYeu/t/R+RGFERESkD7ns1k5BBGBSvo/CtKQ+Pe5XPjWcr3xqeJ8eo7dozIiIiEgfMgz4zOQ8vK5Df///9POn8/xXp5LsUn/AR9QSIiIifSgSTzAk1UMknuC80ZksPCMfWz+eDdUMCiMiIiJ9yO2wMS7Xy8+vPpPMFJeCyBEojIiIiPSxBZPzzS6hX9OYERERETGVwoiIiIiYSmFERERETKUwIiIiIqZSGBERERFTKYyIiIiIqRRGRERExFQKIyIiImIqhRERERExlcKIiIiImEphREREREylMCIiIiKmUhgRERERUymMiIiIiKnsZhfw7wzDACAYDJpciYiIiBzLR7/VH/12d1e/CyPNzc0AFBYWmlyJiIiInIjm5mb8fn+3P28xehpnelkikaCqqgqv14vFYjG7nD4TDAYpLCykvLwcn89ndjmDhtq196lN+4batfepTfvGsdrVMAyam5vJy8vDau3+yI9+1zNitVopKCgwu4yTxufz6f80fUDt2vvUpn1D7dr71KZ942jt2pMekY9oAKuIiIiYSmFERERETKUwYhKXy8UPfvADXC6X2aUMKmrX3qc27Rtq196nNu0bJ6Nd+90AVhERETm1qGdERERETKUwIiIiIqZSGBERERFTKYyIiIiIqRRGToL77ruPc889l6SkJFJTUw9bv2nTJr74xS9SWFiIx+Nh/PjxPProo4dtt2XLFmbOnInH4yE/P5977723x88DGKiO16YAt9xyC1OmTMHlcjF58uTD1u/btw+LxXLYa+XKlX1bfD/WG+0KOlc/6UTatKysjEsvvZTk5GQyMzO5+eabiUQiHet1rh6uN9oVdK4eS2lpKXPnziU1NZWMjAyuv/56WlpaOm1zpPPyqaee6vKx+t0MrINRJBLh85//PNOnT+cXv/jFYevXr19PVlYWS5cupbCwkLfeeovrr78em83G4sWLgUPT8c6dO5dZs2ZRUlLCzp07+cpXvkJycjLf+ta3TvZXMt3x2hQOTVP81a9+lXfffZfNmzcfdV+vvfYaEydO7Hifnp7e6/UOFL3RrjpXOztem8bjcebPn09WVhZvvvkm9fX1XHPNNRiGweOPP95pW52rH+uNdtW5enRVVVXMmTOHK664gieeeIJgMMitt97KV77yFf785z932vZXv/oVF198ccf7bs3IashJ86tf/crw+/0ntO2NN95ozJo1q+P9z3/+c8Pv9xuhUKhj2ZIlS4y8vDwjkUj0dqkDxom06Q9+8APj9NNPP2z53r17DcDYsGFDn9Q2kPWkXXWuHtnR2vRvf/ubYbVajcrKyo5lv//97w2Xy2UEAgHDMHSuHktP2lXn6tH9z//8j5GdnW3E4/GOZRs2bDAAY9euXR3LAGP58uU9Pp4u0/RTgUCg0189b7/9NjNnzuw06cy8efOoqqpi3759JlQ4eFx22WVkZ2fzqU996rDEL12nc7Vr3n77bSZNmkReXl7Hsnnz5hEOh1m/fn2nbXWunrgTaVedq0cXDodxOp2dHn7n8XgAePPNNzttu3jxYjIzMzn77LN56qmnSCQSXT6ewkg/9Pbbb/PHP/6RG264oWNZTU0NOTk5nbb76H1NTc1JrW+wSElJ4eGHH+bPf/4zf/vb35g9ezZXXHEFS5cuNbu0AU3natccqb3S0tJwOp0d7aVztetOpF11rh7dhRdeSE1NDQ899BCRSITGxkbuvvtuAKqrqzu2+9GPfsSf/vQnXnvtNa688kq+9a1vcf/993f5eAoj3fTDH/7wiAN3Pvlat25dl/e7detWFixYwPe//33mzp3baZ3FYun03vjXIKt/Xz5Q9VWbHk1mZia33XYb55xzDmeddRb33nsvN954Iz/5yU967Rj9wcluV9C52tU2PVK7GIbRsVznat+065G2GWzn6r870TaeOHEizz33HD/72c9ISkoiNzeXESNGkJOTg81m69jfd7/7XaZPn87kyZP51re+xb333stDDz3U5bo0gLWbFi9ezJVXXnnMbYYNG9alfW7bto0LL7yQ6667ju9+97ud1uXm5h6W1Ovq6gAOS/YDVV+0aVdNmzaNZ599tk+PcbKd7HbVuXrIibZpbm4u7777bqdljY2NRKPRY7aXztVjO5F2PRXO1X/XlTa+6qqruOqqq6itrSU5ORmLxcLDDz/M8OHDj/rZadOmEQwGqa2t7VIbKox0U2ZmJpmZmb22v61bt3LhhRdyzTXXcN999x22fvr06dx9991EIhGcTicAq1atIi8vr89/oE+W3m7T7tiwYQNDhgwxtYbedrLbVedq10yfPp377ruP6urqjnNv1apVuFwupkyZctTP6Vw9thNp11PhXP133Wnjj0LFL3/5S9xu92G99p+0YcMG3G73UW+3PhqFkZOgrKyMhoYGysrKiMfjbNy4EYBRo0aRkpLC1q1bmTVrFhdddBG33357R1K32WxkZWUBhxLqPffcw1e+8hXuvvtudu3axf3338/3v//9QdudeCzHa1OADz/8kJaWFmpqamhvb+/YZsKECTidTp577jkcDgdnnHEGVquVl156iccee4wHH3zQpG9lvt5oV52rnR2vTS+66CImTJjAokWLeOihh2hoaODb3/421113HT6fD0Dn6hH0RrvqXD22J554gnPPPZeUlBReffVV7rjjDh544IGOoPHSSy9RU1PD9OnT8Xg8rF69mu985ztcf/31XX/Cb4/vx5HjuuaaawzgsNfq1asNwzh0i+SR1hcVFXXaz+bNm43zzjvPcLlcRm5urvHDH/7wlL397HhtahiGMXPmzCNus3fvXsMwDOPXv/61MX78eCMpKcnwer3GlClTjN/85jfmfKF+ojfa1TB0rn7SibTp/v37jfnz5xsej8dIT083Fi9e3Ol2U52rh+uNdjUMnavHsmjRIiM9Pd1wOp1GcXGx8fzzz3da/8orrxiTJ082UlJSjKSkJGPSpEnGI488YkSj0S4fy2IYmmpOREREzKO7aURERMRUCiMiIiJiKoURERERMZXCiIiIiJhKYURERERMpTAiIiIiplIYEREREVMpjIiIiIipFEZERETEVAojIiIiYiqFERERETGVwoiIiIiY6v8BU1vZL9TrCp4AAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Set up the grassland URL\n",
    "grasslands_url = (\"https://apps.fs.usda.gov/arcx/rest/services\"\n",
    "    \"/EDW/EDW_NationalGrassland_01/MapServer/0\"\n",
    "    \"/query?where=1%3D1&outFields=*&geometry=&geometry\"\n",
    "    \"Type=esriGeometryEnvelope&inSR=4326&spatialRel=esriSpatialRelIntersects&outSR=4326&f=json\")\n",
    "\n",
    "# Set up a path to save the data on your machine\n",
    "grasslands_dir = os.path.join(data_dir, 'grasslands_dir')\n",
    "\n",
    "# Make the grasslands directory\n",
    "os.makedirs(grasslands_dir, exist_ok=True)\n",
    "\n",
    "# Join grasslands shapefile path\n",
    "grasslands_path = os.path.join(grasslands_dir, 'S_USA.NationalGrassland.shp')\n",
    "\n",
    "# Only download once\n",
    "if not os.path.exists(grasslands_path):\n",
    "    grasslands_gdf = gpd.read_file(grasslands_url)\n",
    "    grasslands_gdf.to_file(grasslands_path)\n",
    "\n",
    "# Check that the URL is correct\n",
    "grasslands_gdf = gpd.read_file(grasslands_path)\n",
    "\n",
    "# Display grasslands_gdf and plot grasslands to check download\n",
    "display(grasslands_gdf)\n",
    "grasslands_gdf.plot()"
   ]
  },
  {
   "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>OBJECTID</th>\n",
       "      <th>NATIONALGR</th>\n",
       "      <th>GRASSLANDN</th>\n",
       "      <th>GIS_ACRES</th>\n",
       "      <th>SHAPE.AREA</th>\n",
       "      <th>SHAPE.LEN</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>187248</td>\n",
       "      <td>295509010328</td>\n",
       "      <td>Sheyenne National Grassland</td>\n",
       "      <td>70428.175</td>\n",
       "      <td>0.033356</td>\n",
       "      <td>4.097398</td>\n",
       "      <td>MULTIPOLYGON (((-97.31082 46.51457, -97.3056 4...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   OBJECTID    NATIONALGR                   GRASSLANDN  GIS_ACRES  SHAPE.AREA  \\\n",
       "3    187248  295509010328  Sheyenne National Grassland  70428.175    0.033356   \n",
       "\n",
       "   SHAPE.LEN                                           geometry  \n",
       "3   4.097398  MULTIPOLYGON (((-97.31082 46.51457, -97.3056 4...  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "<Axes: >"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVgAAAGdCAYAAABXSWtHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA7eklEQVR4nO3deXhU9aE38O+ZmWSyTSYLhEASlrAk7EuIoAYQQUGxRKiWKpUW8WpvucpyL6VofSyVt0G0b1HRSwErat82tRflFrERFwIiFkMCCBiIhCUhCwGyzCQhk2TmvH9MZmCSSebMzDkzk+T7eZ48kDPnnPmdQL755bcKoiiKICIi2an8XQAiop6KAUtEpBAGLBGRQhiwREQKYcASESmEAUtEpBAGLBGRQhiwREQK0fi7AHKxWCwoLy+HTqeDIAj+Lg4R9WCiKMJoNGLAgAFQqTqvp/aYgC0vL0dSUpK/i0FEvUhpaSkSExM7fb3HBKxOpwNgfeDIyEg/l4aIejKDwYCkpCR77nSmxwSsrVkgMjKSAUtEPuGqOZKdXERECmHAEhEphAFLRKQQBiwRkUIYsERECmHAEhEphAFLRKQQBiwRkUIYsERECmHAEhEphAFLRKQQBiwRkUIYsERECukxq2mR/N7+6gKu1ZugEgQIggCVAPxwUiKSYsIczvv0uyu4eK0BKpUAtQCoVAJUggC1yvqhaftzSJ9wjEuM8s/DEPkBA5Y69be8UpypNDocu21ITIeA/fDYZXx8stLl/RZPGciApV6FTQTkM9zJh3obBix5TRSlnSeACUu9CwOWfEbFfKVehm2w5Jb9Z6pQWt0IQbDWR1WCgLLaG5Ku5W6/1NswYMkt27+84PG1otS2BKIegk0E5DMW5iv1MgxY8hkza7DUyzBgyWfYREC9DQOWfMbMNgLqZRiw1Kkgtbz/PVrNDFjqXRiw1Klgjbz/PZrNFlnvRxToGLDUKa3cAdvKgKXehQFLnZI9YFmDpV6GAUud0mrUst6PNVjqbRiw1CltkMydXBxFQL0Mp8pSp4I9GEUwbXgfaDUqVNQ14XS5weE1DtOi3oYBS53ypAb7fx4ci4GxYSivvYGjl2ogiiLMFhGiCESHBylQSqLAxYClTnnSBmubDjsgKhTzo0LlLhJRt8KAldmu/Mu4dL0BZlGE2QJ7Dc5s+9MiwiKKsFisK/xbP24u/ScIN/8UICA6LAhPzxrus/LnXaxGlcEElQCU1UhbhlAfGoT4yBCYRRFqPy9JeKPZjANFVzF6QCQKSmoQFRaMGSP6+rVM1HsxYGX2wbHL+Orcddnu118f4tOA3fHleXxy+opb18wdHY+XHhqnUInc868L11F0xYhZI+Owfs93iA4Lwls/TcfgPuH+Lhr1QhxFIDOVzDU4CxdIcUtqvA4/vWMwqhuaMX/8AOhDg/DO1xf9XSzqpViDlZncq/Z3h573A0VX8fjOPKhVAtbPH40BLtpeLRYR9c2tKK1uxMCYMOhC5Ov86q+3vrc+NAi/mT8ap8rqZF9T4VyVEb/adRKtFhG/njcSkwfHyHp/6jkYsDKTe98pXwesJxsTVhqaUGloAgD8ck6Ky/OrjCasyD6G78oN+Pw/Z8gasO2NSdDLvkxivcmMo5dqAAC1jS2y3pt6FgaszOTu4vF1wI5PirKPBDhXVY8L1xpcXqML0aCvTguzRYRGQm0xXh+C7Cen4nS5QfZwfffri/i8sAoAYPvKrZ2TgtEJepfXNrWYcdVosnZCitbmGbHt77bOSVG0fl1suIg4dYUBKzP522BlvZ1L/37XUPvff7/vLF7/4pzLa+4bE49ND413630EQcAYCaHnruKqehwouupw7MXM0ZKuPXqxBj9564hb78d8pa6wk0tmvbENtq9O6+8i2Dlrb5X6NQxSu/9vx05I6gprsDKTexhoIP8KunL2cAyKDUOE1rcztC7XNOKv35TYa/cCbn7d80tqOpzf1CJtkZkgD1YP2/Hleew9WYHY8GD8NnOM29dTz8aAlZncnVz+3MfK1aM8PDkJZrMIQ5NvO3re+9cl/PHAeaevJfftON7145MVGDUg0uV9g1TuB2xBSS1QUosB+hAGLHXAgJWZJ73wNpEhGtw+NBZmi4gW881ZX/6SFBOG24bEoNVssQYJrMOfYsKD0dK2tuvA2DCflun81Xp8V27AtOF98OX315y83rFTrqS6UdK9NR40EdhwpTByhgErM28CcUifcPzxsckylsY7D09OwsOTkwAAw579GK0WEQ+lJeL5B0b5rUzJfSPw3rIpAKy1+5Tnc1yuM/vPUxUY+5sqqAQBKf10+LasFrHhWgyPi0BpTSNiw7W4cL0BV40mj8vVHdrKyfcYsDLz6vvMz/P4u2J7rEAqoSAIsEj4greYRbSYWwEAYcFq3DMqHiEaFf6efxkAUHy1AemDo2FsapHcXtsed2sgZxiwMvOmzTSQwqs9W81c6s+AS9cb8M9TlVALAlQqASoB6BOhxQ/GD/C4DFWGJsRFhjgcc7cTMLfdEC6bvIs1GJ8UhROltR6VrYUBS04wYGXmTQVW7g4yOdlyTOowtOKr9dj4zzMOx0b1j5QcsM2tFhRWGJAQHYo+EVpcqzeh3tSKuFvOyfq4UNZxqJ6GK2CtJRO1x4CV2dzR8RjaNxxmy82ZQGbbzKC2Tivba7ZlDG3Hh/SJ8Hfxnbq1Vu6LnwGiKOJkWS0mDYzG5ZobOFFai6iwICT3dfz6lNVKW07RF2xLUaoD+ack+RwDVmY/Sk/ydxFkd2stUWoN1lnN8vsqI+7c+AVEUbTX9EURECG2/Wn14S/uwLjEKAiCgKSYMCTFOB+pEGiD/FvMFqhV8m4U2d4P//swzlYaAdz8Ybfz8duQNiha0fclzzBgyaVbY8ybfrgWsyip1imKzmdktRdoPfctZgtCgpQN2KYWM+pNrQ7HAu3rQDcxYMklV00ERVeMOHKh2tqh1bYjw5m2WpYnpI4pDbR+JV9sSx7AA03ICQYsueSqBvvNhWo8v/uUbO8ntUbmz1luzvhisoGzxYQC7etAN3GxFwLQ9TepKALRYUEIC1ZDJQiKf0NLbVsNtHUafFODZRW2OxHEHvLjz2AwQK/Xo66uDpGRruedk6PKuibE60MgimKX38SbPyvCa59/j2CNClqNGlqNCs1mi6wLT6fG6xCsUdlHXNg2jrx1nVazRcRVowkmCaH2g/ED8G/Thjgc++rcdbyU4ziMLKWfDr//kXXZRdt3hbGpBY/ukLaE4bC4CGjaRhHYrn/1kQlIjZfv/+Pp8jo0mMxtm2Jaf6MY0U+n6KLl1JHUvGETAQGw7ib70bflSBsUjSenD+30POuQMusKVZ7OenLFm/ZbZ6JCgzAuMcrhWLmTzraQYHWHNWprG5slv8+tC3HbNDabJV8vxegB8q+hS8phwBIAYFfBZagEAf82LRnfXzHi8zNViNNpcarMgGCNqq3GqrJvlQIA88cPgFkUUVxVL3soyqnV0vEHgdOF0RX4ZU7KVF7quRiwBMDasfRl8TW8/MlZ/OyOwXj5k7P27VI68/LD46DVqPHO4Yt44R+nfVdYNzmbZeWrCQHM196NAUsAgNa25RHfzC3GqAGRCFILaGrpOh1stcBAb8ZvdTKey1nAKvEUHKPauzFgCQCwbUkaRABqQYBGLeDUb+bg6b8ewz9PVTqc94u7huLhyUloMVtuduhIfI/kvuFYNXsERNzsoNnwUaF9R1qlOKvBjk3QY8cS69KQtld1IfJ/OwTabDPyLa+GaWVlZUEQBKxcudLheGFhIebPnw+9Xg+dToepU6eipKSk0/vs3LkTgiB0+GhqUvYbj27ShQQhMiQI4VoNtBo1NGqV03GdMeHBGNInHCP66eyjDaRmSJ9w62pa88cPwA/GD8AD4wYoEmrtOVvpKjZCi9mj+mH2qH64p+1janKs7O/NgO3dPP7fnZeXh23btmHcuHEOx4uLi5GRkYFly5Zh/fr10Ov1KCwsREhISCd3soqMjMTZs2cdjrm6hpTl7NdbpwPdpd7QSbOn3LvwOuPPpQTZRNC7eRSw9fX1WLx4MbZv344NGzY4vPbcc8/h/vvvx6ZNm+zHkpOTXd5TEATEx8d7UhxSiLNgctp2KbGW5qxfyRfj5v25nQsrsL2bRwG7fPlyzJs3D7Nnz3YIWIvFgr179+KXv/wl5syZg2PHjmHIkCFYt24dHnzwwS7vWV9fj0GDBsFsNmPChAl48cUXMXHixE7PN5lMMJlubvFhMBg8eRTJzlQa0GAyd1wJqq2n3XZ0bIK+w6Dvoxercbnmhn2wvCgCoxMiA35MY6uTtkuVF73vvqitOuOLGVadYRNB7+Z2wGZnZ6OgoAB5eXkdXquqqkJ9fT02btyIDRs24KWXXkJOTg4WLlyI/fv3Y8aMGU7vmZqaip07d2Ls2LEwGAx49dVXceedd+LEiRMYPny402uysrKwfv16d4vvsTV//xYny+pcnvfhL+7AxIGOS8ftPHwRH31b4Xi/OSkBFbDL/1KAS9cbIIrWoUWiKOLS9Y6bBb791QXsO10JtUqAWhCgVgmSNxW8eK0BGz76DoJws2ZX5cU+WFJ5U4MNCVLjuftHdlhSURQBtcq6yaXQtsDNrX8KbQvfjOink+chqFtyK2BLS0uxYsUK7Nu3z2n7qKVtQHdmZiZWrVoFAJgwYQIOHz6MrVu3dhqwU6dOxdSpU+2f33nnnZg0aRJef/11vPbaa06vWbduHVavXm3/3GAwIClJubVYpe44KrXdMtB8f8WIoisdZyK1d/5qg9OdW6Uor2vCjkMXPLrWG9589UOC1Pi36a6buIiccStg8/PzUVVVhbS0NPsxs9mMgwcPYsuWLWhoaIBGo8GoUY67jo4cORKHDh2S/D4qlQrp6en4/vvvOz1Hq9VCq9W6U3yvBKmkDbhwNiSIi9z7Vzf4+UY9lFsBO2vWLJw8edLh2NKlS5Gamoq1a9dCq9UiPT29w2iAoqIiDBo0SPL7iKKI48ePY+zYse4UT1E9vQbbkwkBvZ0k9WRuBaxOp8OYMWMcjoWHhyM2NtZ+fM2aNVi0aBGmT5+OmTNnIicnB3v27EFubq79miVLliAhIQFZWVkAgPXr12Pq1KkYPnw4DAYDXnvtNRw/fhxvvPGGl48nH6lTK1uczXtnFZaoV5J9lPeCBQuwdetWZGVl4ZlnnkFKSgp27dqFjIwM+zklJSVQ3fIrd21tLZ588klUVlZCr9dj4sSJOHjwIG677Ta5i+cxKVuYAM4X9+gO+dqjO7u7wdefeiavA/bWmqnN448/jscff1zyNX/4wx/whz/8wduiKEojtQbrtA2W3+H+xK8++Qt3NJBIahuss6XxuAq9f/HLT/7CgJVII3EUgfNOLrlLQ+5gJxf5CwNWovTB0YjQut6S2Tb7SRRF+0dnNahAX+avp5D4s5FIdtyTSyJRFHGm0oCfvZ2HKwZ5Zx/ZAti6hJ9gX8pPgIAfjB9g3yfKJv9SNdb8/Vto1AI0KhWC1AKC1Cpo2v4MvuXv1g8BmrbjQWrh5qiGW/7l/3a0VNZ9taS6c1gsdNqgttlP1mPfXKjBtXrXX+OJA6MwqW3WnNg2Xdn2v9k2pVkUrcskLr1zSOc3InIT9+SSmSAIGNlfjw9+cSd++qdvnO6/5Cl7KNz6SdsRs5M23cZmM85f82w2VaD5zQ9GY3i76aSPvXUEX37vOmCnD++LVfeMUKpoRF7jL09uSogKxa6f34HbBsf45P2cdZD1jN85vMfOKwp0DFgP6MOC8O6y2zAzpa/i79XUYkaVoQnX6032rU96er5K/QHCzisKdGwi8FBIkBrx+hCMS9Tj28uuV9ny1D9PVdq3bclZOQ2p8ZE9qnPM6eauEn+ESJz7QeQ3DFgvXDWaUFih7Dq0t1qZfRwRWg2MTa0+e0+lFVYYYWhqdai1Gm50fL5R/SOREB1qXSJRLUCjEpASL39nJpGcGLBeMFtEpzO3lHKm0uiz9/KVp/96TNJ5T0wbgoWTEhUuDZG8+EuWF/y5FUlvww4t6o4YsF5wth3IiH4R0IcGOTmbvMEOLeqOGLBecDYttuhKPQZEhSJO57vFwHsD1mCpO2IbrERPvJOHi217VNlmWnW2F1VhhQHJfcMRrFHhcs0Nye/xUFoiHhjXHyKsyx6KIrD/bBX+35ESGZ6g++gTEYywYA3UKmtnlkatQlRYsL+LReQ2BqxE56+5txfV+asNSIwORXLfcMnXDe0bgbtS4hyOlddJD+ie4s3FabhtiG8mchApiU0EEjlbSNuVyzU3YGxqxcj+0oYTOfs12FkzBBF1DwxYicweDu6/ajShrKYREwdGuTxX7SRhe2O+9qSJFNS7sYmgnRvNZvzjRBkAxymb3gzuNzS14nS5AbcNicE3F6o7Pc/prCaGDVG3xYBtp6axGWt3nXR9opuaWy04erEaU4bEIP9SjdMxtB+frMDF6w1QCYL941S5ctNwA5VFBFrNFphFERaL9bcHiyjCYhFhtoiwiNYfPGFaDSK0/C9MgYv/O9txNrZVvnsDRy5UI31wNE6XG9DYbHZ4vaCkFgUltYq9f3fxyPZ/STrvP+8ZgadnDVe4NESeYxtsO06WX5Vd3sUaDI+LQHQYJyR4w9N2cSJfYcC246tv2hOX69AvMgR9OSHBY56M7CDyJQZsO74cFnWm0ojIEA0So0N99p49CWuwFOjYBttOv0gt3nh0EgDbvlhWK7KPo9ksf/tBsQcTEqSK0GqQNigaZouIVosFFot1W3GLaG1rNts7jax/2meQAW2dStbOJGsnU9vf2143W0Sg7T43Wsx+GU6mwD8HkawYsO3oQoIwb1z/DsdXZB9X7D0v19xAnE6Lkf0jZV1fdlBsGN55/DbZ7teZ5X8pwN5vKxR/n/aU7JAkkgObCCRqVbj3q8poQnntDYxPilL0fZTgr7G6HCNMgY41WAksbWMvlVZ3owVFlUakD45G3sUa5d9QJs7ardMHRyM1PhKtFgtazdZmhta2xclvPXv/mSrUmzybxMEmAgp0DFgJfNmZcqPFjIKSWtzWNiHBm043Xy3x56yIc8f0x7KMIS6vnf1/D3i8BbrUvbuI/IUBK4EA4JdzU+xTZ22/mlo7fqzf6NY/ra+JojWUxbaOoY++rUBZrfRVscwWEd9cqEbaoGgUlNR4vE13c6sFl2vallgUrEtWa1QC4iJDPLthJ5yVzxfZzhYCCnQMWAk0ahV+cdcwj68/WVbnVsDa5F+qQUJUKBqbW1HT2OL29UVX6pHx0n6HY311WuQ9N9vte3XFWVuo1NqzN+2o7OSiQMdOLh/wZu+ustobiNOFIF6mWqdKgaqls6BT+aB9gvlKgY4B6wPezjg6e8WIkCAVBsaEeV0WJfa2cvZ4kmuwXr0vE5YCGwPWB+TYffbi9Ua0mi0YFhfh1X16VA1W8Xcg8g7bYD1ktoh4/2ip/XNbnKhV1mUG1SoBgmANmrobLVC3S7bORge0P+/WNsorRhNiLSImD4qGsan15mystqX8gFtmYEG0d7TZZm6JoohwBZb3c9rJ5YNeLq5FQIGOAeuhVosF6z6Qf93YLodliSKqjCaYWi1466eTMXmw+/tWudOpNHfzQZhaLRBg/UEhCDdHI6hUbccAPD1rGDY8OBqtFhGtFhFms4h+emltxhsyx6Ch3bKNOacqUFrdCBFAkFoFlQA0tVgHvdpG0ooiEB3OjRApsDFgPeTP5r+6Gy1YvOMI3nh0EmaP6ufWtYIbVctL1xtxo8Xs8rzw4CAM7uNZ08Udw/p0OLYr/zK+uWWixR1DY3H0UseJF942lxApjW2wHvJ3B4up1YKn/pzv0EwhN6kD+eUe8N++maSzXQvYx0WBjgHroUDY7dVsEfHL//kWb+w/16Pm5avaBWyV0eT0vJ7zxNRTMWA9FAD5avfyJ2ex5Ytz/i6GbNTtWjFOldUhIapjm24P+plCPRTbYD0UaDXGwkppyxyWVjfiJ28dsXdaqQQBEVoNdi+/s8O57z91u33915uPe3NkgihaX0uN13W8Nq8UX5+/Do1KgEatQpBawIwRfTFrpOs244fSkpA2OAYCrM0FAqzNBI0tZnvHmlolIEmGccFESmLAeigQmghuJXXcqanVgkvXGx2O6Tpp4xyXGOVxeb4tq8WHx8ocjoVrNZICNmN4H2QM79j5RdTdMGA9FGD56sboAN8U3NmMsZZWCxqbW+2vty+y7XONStWho4uoO2LAeijQmgjat1t2xmmxFcgyZ3m/49AF7Dh0weW1Lz44Bo9NHSR/oYh8jAHbhUvXG+zjQG01Mts+Xdcbmv1Yso58MTXVVzhDi3oKBmwX1vz9W3xzsdrfxZDGSb5+e7kWGz4qxIyUvjhcfA2iCKcTB240m/HEO3k3Z2nZZ23d8qt+2x/PzxuFeAmztJzF/W1DYjBjRF98f8WI3cfLO7020Nq3iTzFgO0pnGSSsakV31ysRt6lajyclojdx8qd7ozbahHxWWGVpLdZOWu4pIB1Jm1QNJbPHIZ9pyu7DFg2v1JPwYBt53R5HXZ+dREAcP6avNtoK6mrplVRBN4/ehmZEwbgizNVMDZ5tgcWIL1zr8vyuLjWnem8RIGMAdtORW0T/p5/2d/FkEe7nPrf4+WYNTIOp8rqcMXgfHaUK95Mi7Xlpqv+QeYr9RScydVDON22xUlL6OeFVRgUE4ahfcM9fB/Pz7vZEdf1TXpShx31bqzBttNdv7e/Kr6On/7pG6iEm51UtZ3s4/XNxRqM7K/D+EQ9Tlyuc+t9pC5ykzG8D0KD1Wg1i/Z1awVBwP87cgnflXc96+xYSS2CNSqIomhfyzYyJAi1jc32z/tEaDGinw4pTmaREQUKBmwPcdVowgHjVcnnF1YYkRQTioxhsTh07rrk66TWYOeMjsec0fEOx/7waRGe+/CUy2t3FVzGrgLHZprJg6Idlizsq9Pi6buHMWApoLGJoBcrrb6BwgoD5rYLwq54s0yjN4Ov2i+ufdVowqHvr3lxRyLlMWDb6a5NBJ663tCCg0VVWDgpQdL5Xk1g8+LiJifjdzlclgIdA7YdJXZdDXSNLRb843g5Fk1OdLkGgFf56sW1zi9mwlJgYxtsm4NFV/FM9jG0tN4ciP/m4kmYmhyLc1X1+NEfv/Zj6RzdnRqHR28baN/s0CyKyLtQjXe+vuTxPVstIv529DKemp6MyNAgeweTKAJm0foXiwjE6bQev4fUCuzwuAj0jwqFSgDUggBBEDAsLhxhWrXDMotjE/Qel4XIFxiwbVotlg697hFaDWLCg6EPDfJTqZxLiArtsBeXRYRXAWtjarVg+cxhXt/HmfY7FXTm3+8aioWTEhUpA5EvMWDbOGsa+Kr4GqobmlFR1+SHEnXOWUeTXNNLlVwlTCOxkFyqkHoKBqyNk+/pPx447/tyeEiutmOzggErNTgZsNRTsJOrTXf6lnYWgXJlkpI985JrsL1tKAf1WKzBtpFribzdy++0B4ltub/H3joi6/qxh89dw3++fwIalYA1c1PQJ0Ir2/AyJddijdeHYHyi3mGPL7Ftj6+bnyPg2ryJPOVVwGZlZeHZZ5/FihUrsHnzZvvxwsJCrF27FgcOHIDFYsHo0aPx/vvvY+DAgS7vmZ2djUceeQSZmZnYvXu3N8Xzi7EJ+g6/4mqkbjcg0cXrjbjYtq/Wz+8a2haw8ryHNxMJXMmckIDMCdLG2xL1BB43EeTl5WHbtm0YN26cw/Hi4mJkZGQgNTUVubm5OHHiBJ5//nmEhLheQ/TSpUv4r//6L0ybNs3TYnmsu/5WagtEuRZIcbJcLBF5yKMabH19PRYvXozt27djw4YNDq8999xzuP/++7Fp0yb7seTkZJf3NJvNWLx4MdavX48vv/wStbW1nhTNLWcrDfj0uysAgPNXu8/ar7eyVTjl+vmgZA2WqLfxqAa7fPlyzJs3D7Nnz3Y4brFYsHfvXowYMQJz5sxBXFwcpkyZIulX/d/+9rfo27cvli1bJqkMJpMJBoPB4cNd567U45V9RXhlXxE+aLfFdPdhDUS5auDcroVIPm7XYLOzs1FQUIC8vLwOr1VVVaG+vh4bN27Ehg0b8NJLLyEnJwcLFy7E/v37MWPGDKf3/Oqrr/DWW2/h+PHjksuRlZWF9evXu1t8RxJDqU+EFuvnj3Y41tjcijX/822Hc//9z/kAbL3x1tlPNZ0sGyiHtw5dQN+2Nth/v2uodWaX5eYyf7YPs8XagWURRXx4rAyt7YJ02vA+0GrUGDUgUrGy+oup1Yw627+B4LiBpY1tL7L2i8oAQJWhyTqrra1DDgD6RYZwOBm55FbAlpaWYsWKFdi3b5/TNlWLxdqAl5mZiVWrVgEAJkyYgMOHD2Pr1q1OA9ZoNOInP/kJtm/fjj59+kguy7p167B69Wr75waDAUlJSe48juSOoXCtGvPG9Xc4ZmhqcRqw+9qaHHzlr9+UAgBmpcbhrZ+lS7pmz7flHQL25YfGe7zXVqA7XlKLRdv+5fK82PBg5D9/T4fjt/3u8w7H8n89G7ERnk8bpt7BrYDNz89HVVUV0tLS7MfMZjMOHjyILVu2oKGhARqNBqNGjXK4buTIkTh06JDTexYXF+PixYv4wQ9+YD9mC2qNRoOzZ89i6NChHa7TarXQar37Dy61AtL16vyBwZ1RBK3mjg8k90gHInIzYGfNmoWTJ086HFu6dClSU1Oxdu1aaLVapKen4+zZsw7nFBUVYdCgQU7vmZqa2uGev/71r2E0GvHqq6+6XSt1hzeznwItjqT/sBA71F4BIEjVc+ecSG1VbmhuRdY/CwERuCslDrcPjVV06jD1fG4FrE6nw5gxYxyOhYeHIzY21n58zZo1WLRoEaZPn46ZM2ciJycHe/bsQW5urv2aJUuWICEhAVlZWQgJCelwz6ioKADocFxu3lRCA60GK7U9sLNOrGCNdwHbYrYgSN29Q7qpxWKfHh0THozbh8Z2ei53viUpZJ/JtWDBAmzduhVZWVl45plnkJKSgl27diEjI8N+TklJCVQBUGMa1i8C6+5LdXmeLqTjzKJA+/6SGvjOaq8A8H2VEWHBagSpVfaP2PBgyStgSQlXi0Vsm8XletNDqe9LFMi8Dthba6Y2jz/+OB5//HG3rrnVzp07vSuURMl9IvDUjAiPrg20GqzUQOosYOdv+arDsVPr5yBCK9/P4Df2n8PvPy1yed4fH0vrsJ9XoAmsf30KVP6vRnZTAZavkNpH1erGVK0gmTu+pLZmKrkeglSB9u9L3RMXe/FQwNVgvWwicEbuji+p/UVyL5noye2OnK+GAAEit6UhLzBgPRRoTYSSmwicDNFyJkgtfzuo1LAKhNlkn5+pwudnqvxdDOrmGLAeEgQBS24fZG+La9+rXHy1Hl/6cFvp8tob+OfJik5ft0WW1GUTlRgRILUmyfUQqKdgwHrht5mdDyP73+NlPg3Yw8XXcbj4umz3a2w2Y9xvPgHg2Ha69+lpGBgb5tE9pQ4Fe/aDU1i/5zvrBoew/vBSCbBvdmjb+PDWn2nOMtl2zNTKJcLIPxiw1ClDU2uHY960SY7qr8P4JD1OlNZ1ed6NFjNutJg9fh9fCLAmeApQHEVAbvFm9ltTiwVLpg7CjBF9ZSwRUeBiDZbc4k3N7b6x1gVzJg6MxqPb/4VKg0mmUhEFJgasQuIjQzBndL9OX/+6+LrTX8F7g5LqRswdE4/dx8tRq+BSjkT+xoBVyJTkWExJ7nwue+YbX+FEaa3vCuRj//fTIuw7XQlBEG4ZaWH909jUipLqRjw4YQByTlWiiZ1Q1EMxYMktUsfGltfewJlKY5fn7D5ejicyhuDtwxcDYuyrO7xpi6beg51c5BapsSJ1LOt7/7qEp6a73rONqDtiwJJbpHZySZ0rYGq14L1/XcKyOwd7XCaiQMWA9ZPkPuEY0S8CyX3CMTAmDAlRoU7P66qjzB+k/mrszkLVxqZW7Pm2AosmJ3paLKKAJIg9ZMl2g8EAvV6Puro6REZ2z437Jvx2n6Re9ayFY3H/mP72Qf+i6Djbqqt/0sIKI37y1hGX7xGkFvDJyun2zRPFtj+HxUV0mEb77Icn8f0Vo33GlUoQUHSlHtfq3RuGlRAVipR+Efji7FW3rvOHe0b1Q7Dt6+DkZ47t0NN3D0dKvM5n5SLfkJo37OQKIFqJU0mD1SrowzouAi5FVFiTpPNEEUjuK22t3DMVBhSU1HpUnluV1d5Aq8WCKUNicORCtdf3U9KnEje3fPS2gQAYsL0VmwgCiNS5+q0W5Yc1ubNkoJy/Al0xmFBlNGFMD9w+nHofBmwACZa4gpUba2Z7zJ8NRxeuNQCCgMF9PFtUhihQMGADSLBGLek8X9RgAekdVUqE8amyOsRHhqBvRLD8NyfyEbbBBhCpbbAtEhfNdiYuUos1c1IAWIdc2UYF3Dr8yvZXUfTvqlH/Ol+Ne0b2w9fnr6HeFNiraxE5w4ANIFLbYL3ZsypOF4LlM4d5fL2vfVp4BfPHD0DOqQo0e/GDhcgf2EQQQCTXYH3URBAo/nGiHPPHJ0AdaPv0ELnAgA0gUgM2EHZd9bX/KbiMH05K8Hcx3Nb7/qXoVgzYACJ9mFZgfdv6qp12V0FZtwxZ6r0YsAFE6kaD3W3lKbmYLSL2nqzA/WPj/V0UyXrGPEnyFAM2gLCF0bWmFgsOn7uOu7rJtjPe7GFG3R8DNoDwW1Ga2hstKLpiRPrgaH8XxaVe+ssGtWHAdkP8tRMor2tCTWMLRvYP7Hn+UtfFpZ6JAdsN8ddOq3NV9dBq1EiKcb7UY0DgP1WvxokGAWThpESMHhAJU4sFzWYLmlstMLXtV2W2iBBhXTZwXGKUfwvazpQhsYjTaSFAsI8oyL9Ugyqj8rvGHi+tRcawPmg0mXG9oVnx97OZkBSFCK3GOhuubd8xoW2pRtvfAQEx4Zzq25txPVhSxLKdefj8TJXP3u/eUf3w1blraGj2zZTaz1ZPx7C4wG6eIOVIzRs2EVCPsO+7K46LYCvMFyuaUffHgCVF+OPXot3Hy5E5cYBPJj701rHI5B4GLCnCXy1P/5N/GQ+nKb+3F0cHkBTs5CKvna00oqG51WH4WO0N13uLeSs0SI3YtvVibe8tiiK+uVCNH6cnITuvFACgUQn282wdcdaOqJtVXdtfHZdtvPmJytaZ1Xat1GnN1LsxYMlrv/rgWxyTYU8ud909Mg5vPDrJ6Wti20aNfztaipR4HfY+M83HpSNiwFIPJQgC/s+CMbjeYEJFnbSNHonkxt9zqPty0QyqUavw+iOTkD44xjflIWqHAUteC+RFakKD1fj1vJH+Lgb1UmwiIK+5258epBYQHRYMETfbSi2i2PZ527F2Nze1zWxzfF9p76zx0dhYovYYsORzaYOikf3k7W5d86dDF/Dbj75zOCYEdN2ZiE0E5AeejNF3egnzlQIcA5a85u6Ye08mITi7hvlKgY4BSz4n1yQowVebgRF5iG2w5DW1SuiwpXZntVQR8k0zZbxSoONyhdQtNLda0NK2hJXtP6xGJSAkSO2/QlGvJTVvWIOlbiFYo+L8f+p2+D+WiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSiFcBm5WVBUEQsHLlSofjhYWFmD9/PvR6PXQ6HaZOnYqSkpJO7/PBBx9g8uTJiIqKQnh4OCZMmID33nvPm6IREfmdx6tp5eXlYdu2bRg3bpzD8eLiYmRkZGDZsmVYv3499Ho9CgsLERIS0um9YmJi8NxzzyE1NRXBwcH46KOPsHTpUsTFxWHOnDmeFpGIyK88Wg+2vr4ekyZNwptvvokNGzZgwoQJ2Lx5MwDgxz/+MYKCgryugU6aNAnz5s3Diy++KOl8rgdLRL4iNW88aiJYvnw55s2bh9mzZzsct1gs2Lt3L0aMGIE5c+YgLi4OU6ZMwe7duyXfWxRFfP755zh79iymT5/e6XkmkwkGg8Hhg4gokLgdsNnZ2SgoKEBWVlaH16qqqlBfX4+NGzdi7ty52LdvHxYsWICFCxfiwIEDXd63rq4OERERCA4Oxrx58/D666/jnnvu6fT8rKws6PV6+0dSUpK7j0JEpCi32mBLS0uxYsUK7Nu3z2mbqsVi3dIjMzMTq1atAgBMmDABhw8fxtatWzFjxoxO763T6XD8+HHU19fj888/x+rVq5GcnIy77rrL6fnr1q3D6tWr7Z8bDAaGLBEFFLcCNj8/H1VVVUhLS7MfM5vNOHjwILZs2YKGhgZoNBqMGjXK4bqRI0fi0KFDXd5bpVJh2LBhAKyhXFhYiKysrE4DVqvVQqvVulN8IiKfcitgZ82ahZMnTzocW7p0KVJTU7F27VpotVqkp6fj7NmzDucUFRVh0KBBbhVMFEWYTCa3riEiCiRuBaxOp8OYMWMcjoWHhyM2NtZ+fM2aNVi0aBGmT5+OmTNnIicnB3v27EFubq79miVLliAhIcHejpuVlYXJkydj6NChaG5uxscff4x3330X//3f/+3l4xER+Y/su8ouWLAAW7duRVZWFp555hmkpKRg165dyMjIsJ9TUlIClepm/1pDQwN+8Ytf4PLlywgNDUVqair+/Oc/Y9GiRXIXj4jIZzwaBxuIOA6WiHxF0XGwRETkGgOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFMKAJSJSCAOWiEghDFgiIoUwYImIFOJVwGZlZUEQBKxcudLheGFhIebPnw+9Xg+dToepU6eipKSk0/ts374d06ZNQ3R0NKKjozF79mx888033hSNiMjvPA7YvLw8bNu2DePGjXM4XlxcjIyMDKSmpiI3NxcnTpzA888/j5CQkE7vlZubi0ceeQT79+/H119/jYEDB+Lee+9FWVmZp8UjIvI7QRRF0d2L6uvrMWnSJLz55pvYsGEDJkyYgM2bNwMAfvzjHyMoKAjvvfeex4Uym82Ijo7Gli1bsGTJEknXGAwG6PV61NXVITIy0uP3JiJyRWreeFSDXb58OebNm4fZs2c7HLdYLNi7dy9GjBiBOXPmIC4uDlOmTMHu3bvdun9jYyNaWloQExPT6TkmkwkGg8Hhg4gokLgdsNnZ2SgoKEBWVlaH16qqqlBfX4+NGzdi7ty52LdvHxYsWICFCxfiwIEDkt/jV7/6FRISEjoE+K2ysrKg1+vtH0lJSe4+ChGRojTunFxaWooVK1Zg3759TttULRYLACAzMxOrVq0CAEyYMAGHDx/G1q1bMWPGDJfvsWnTJvz1r39Fbm5ul+2269atw+rVq+2fGwwGhiwRBRS3AjY/Px9VVVVIS0uzHzObzTh48CC2bNmChoYGaDQajBo1yuG6kSNH4tChQy7v/8orr+B3v/sdPvvssw6dZ+1ptVpotVp3ik9E5FNuBeysWbNw8uRJh2NLly5Famoq1q5dC61Wi/T0dJw9e9bhnKKiIgwaNKjLe7/88svYsGEDPvnkE0yePNmdYhERBSS3Alan02HMmDEOx8LDwxEbG2s/vmbNGixatAjTp0/HzJkzkZOTgz179iA3N9d+zZIlS5CQkGBvx920aROef/55/OUvf8HgwYNRWVkJAIiIiEBERIQ3z0dE5Deyz+RasGABtm7dik2bNmHs2LHYsWMHdu3ahYyMDPs5JSUlqKiosH/+5ptvorm5GQ899BD69+9v/3jllVfkLh4Rkc94NA42EHEcLBH5iqLjYImIyDUGLBGRQhiwREQKYcASESmEAUtEpBAGLBGRQhiwREQKYcASESmEAUtEpBAGLBGRQhiwREQKYcASESmEAUtEpBAGLBGRQhiwREQKYcASESmEAUtEpBAGLBGRQhiwREQKYcASESmEAUtEpBAGLBGRQhiwREQKYcASESmEAUtEpBAGLBGRQhiwREQKYcASESmEAUtEpBAGLBGRQhiwREQKYcASESmEAUtEpBAGLBGRQhiwREQKYcASESmEAUtEpBAGLBGRQhiwREQKYcASESmEAUtEpBCNvwtARNSZ3cfK8LuPC6FWCVAJAgQBUAkCVG1/bluShmFxOn8Xs1MMWCIKWI3NZlQZTZ2+bmq1+LA07mPAElHAEiF2+boAAQBwvd6E7yoM9s8FAda/CTfPUQlATHgwhvfzXY2XAUtEAUvsOl8hWLMTJy7X4vGdR13eb+7oeGx9LE2GkknDTi4iClgu8tUesK6C2CZI49vIYw2WiAKXi+T8z/dPQBeigeFGq6Tb5V2oxhPv5NlvLQK4d1Q//Pi2gd6W1CkGLBEFLIuLmunpcoNb96s0NKHS0ORwLCRIpVjAsomAiAKWRerv/l5Q2doZlLi3YncmIvKS2VUVVgYaFQOWiHohH1RgoVIwYNkGS0QBq1XmGuzMlL54KC3JPvpAAJAQHSrre9yKAUtEAatZ5plaw/vpMG9cf1nv2RU2ERBRwGo2m2W9n5IdWk7fz6fvRkTkBrlrsEFqBiwREQD5F3PxdQ2WbbBEFLBazPJ2cn1yuhIXrzcAAH7/8Hho1MrWMVmDJaKAJco8TutMpRH/e7wc/3u83OUsMTmwBktEAUvucbAx4cGICQ+W96Zd8KoGm5WVBUEQsHLlSofjhYWFmD9/PvR6PXQ6HaZOnYqSkpJO73P69Gn88Ic/xODBgyEIAjZv3uxNsYiInPrp7YPx2eoZ+Gz1DAT7YGUtj98hLy8P27Ztw7hx4xyOFxcXIyMjA6mpqcjNzcWJEyfw/PPPIyQkpNN7NTY2Ijk5GRs3bkR8fLynRSKiHsbVgtuBzqMmgvr6eixevBjbt2/Hhg0bHF577rnncP/992PTpk32Y8nJyV3eLz09Henp6QCAX/3qV54UiYh6oJ9MHYSZKXFtM6+sIwD+8k0JDhZddes+v394PCyiiJH9I+UvZBc8qsEuX74c8+bNw+zZsx2OWywW7N27FyNGjMCcOXMQFxeHKVOmYPfu3XKU1YHJZILBYHD4IKKeZVxiFO4b2x9zx/TH3DHxmDsmHoluTm0VBOCHaYl4eHISxiToFSqpc24HbHZ2NgoKCpCVldXhtaqqKtTX12Pjxo2YO3cu9u3bhwULFmDhwoU4cOCALAW2ycrKgl6vt38kJSXJen8iCkzOOr50IRokRIWiv75jU6QvFozpjFtNBKWlpVixYgX27dvntE3VYrEOCs7MzMSqVasAABMmTMDhw4exdetWzJgxQ4YiW61btw6rV6+2f24wGBiyRL1Cx8R8IiMZK2YPh8UiIvnZj/1QJufcCtj8/HxUVVUhLe3mpmFmsxkHDx7Eli1b0NDQAI1Gg1GjRjlcN3LkSBw6dEieErfRarXQarWy3pOIAp+zGmmgdoa5FbCzZs3CyZMnHY4tXboUqampWLt2LbRaLdLT03H27FmHc4qKijBo0CDvS0tEvd7cMfEY3CccAm5OfZ00KBqAtb113X2pEGHbc0vsPk0EOp0OY8aMcTgWHh6O2NhY+/E1a9Zg0aJFmD59OmbOnImcnBzs2bMHubm59muWLFmChIQEeztuc3MzvvvuO/vfy8rKcPz4cURERGDYsGHePB8R9TB3pcThrpQ4p68JgoCnZgz1cYk6J/tI2wULFmDr1q3YtGkTxo4dix07dmDXrl3IyMiwn1NSUoKKigr75+Xl5Zg4cSImTpyIiooKvPLKK5g4cSKeeOIJuYtHROQzgij3ZF8/MRgM0Ov1qKurQ2Skb8e6EVHvIjVvuNgLEZFCGLBERAphwBIRKYQBS0SkEAYsEZFCGLBERAphwBIRKYQBS0SkEAYsEZFCGLBERAphwBIRKYQBS0SkEI82PQxEtjVruDcXESnNljOu1srqMQFrNBoBgNvGEJHPGI1G6PWdb6TYY5YrtFgsKC8vh06ng9C2ynl3YNtLrLS0tNsus9jdn6G7lx/gM/iaKIowGo0YMGAAVKrOW1p7TA1WpVIhMTHR38XwWGRkZMD/p3Kluz9Ddy8/wGfwpa5qrjbs5CIiUggDlohIIQxYP9NqtXjhhRe69Rbk3f0Zunv5AT5DoOoxnVxERIGGNVgiIoUwYImIFMKAJSJSCAOWiEghDFiFFRQU4J577kFUVBRiY2Px5JNPor6+3v76zp07IQiC04+qqiqX9xdFEffddx8EQcDu3bu71TM89dRTGDp0KEJDQ9G3b19kZmbizJkz3eYZqqur8fTTTyMlJQVhYWEYOHAgnnnmGdTV1XWL8gPAtm3bcNdddyEyMhKCIKC2tlb2siv9DCaTCU8//TT69OmD8PBwzJ8/H5cvX1bsOdzBgFVQeXk5Zs+ejWHDhuHIkSPIycnB6dOn8bOf/cx+zqJFi1BRUeHwMWfOHMyYMQNxcXEu32Pz5s2KTg1W8hnS0tLw9ttvo7CwEJ988glEUcS9994Ls9ncLZ6hvLwc5eXleOWVV3Dy5Ens3LkTOTk5WLZsWbcoPwA0NjZi7ty5ePbZZ2Utsy+fYeXKlfjwww+RnZ2NQ4cOob6+Hg888IDs/488IpJi/vjHP4pxcXGi2Wy2Hzt27JgIQPz++++dXlNVVSUGBQWJ7777rsv7Hz9+XExMTBQrKipEAOKHH34oV9HtlH6GW504cUIEIJ47d86rMrfny2d4//33xeDgYLGlpcWrMt/KF+Xfv3+/CECsqamRo8gdKPUMtbW1YlBQkJidnW0/VlZWJqpUKjEnJ0e+B/AQa7AKMplMCA4OdlgMIjQ0FABw6NAhp9e8++67CAsLw0MPPdTlvRsbG/HII49gy5YtiI+Pl6/Q7Sj5DLdqaGjA22+/jSFDhsi+IpqvngEA6urqEBkZCY1GvmU+fFl+pSj1DPn5+WhpacG9995rPzZgwACMGTMGhw8flqn0nmPAKujuu+9GZWUlXn75ZTQ3N6Ompsb+q1hFRYXTa/70pz/h0Ucftf/n68yqVatwxx13IDMzU/Zy30rJZwCAN998ExEREYiIiEBOTg4+/fRTBAcHd6tnsLl+/TpefPFFPPXUU7KU28ZX5VeSUs9QWVmJ4OBgREdHOxzv168fKisr5XsADzFgPfCb3/ym08Z428fRo0cxevRovPPOO/j973+PsLAwxMfHIzk5Gf369YNare5w36+//hrfffedyza8f/zjH/jiiy+wefPmbvsMNosXL8axY8dw4MABDB8+HD/60Y/Q1NTUrZ4BsC61N2/ePIwaNQovvPBCtyu/pwL1GURRDIhlSzlV1gPXrl3DtWvXujxn8ODBCAkJsX9+5coVhIeHQxAEREZGIjs7Gw8//LDDNcuWLUNBQQGOHTvW5b1XrlyJ1157zeHXLbPZDJVKhWnTpiE3Nzfgn8GZ5uZmREdHY8eOHXjkkUe6zTMYjUbMmTMHYWFh+OijjxzerzuUHwByc3Mxc+ZM1NTUICoqSvJ1/n6GL774ArNmzUJ1dbVDLXb8+PF48MEHsX79esnPogg/twH3Om+99ZYYFhbWoTPBaDSKERER4uuvv+7yHhUVFeLJkycdPgCIr776qnj+/HmFSn6THM/gjMlkEkNDQ8W3337b+0K6INcz1NXViVOnThVnzJghNjQ0KFBS5+T+N1C6k8sZOZ7B1sn1t7/9zX6svLw8YDq5GLAKe/3118X8/Hzx7Nmz4pYtW8TQ0FDx1Vdf7XDejh07xJCQELG6urrDa5cvXxZTUlLEI0eOdPo+UGgUgSgq8wzFxcXi7373O/Ho0aPipUuXxMOHD4uZmZliTEyMeOXKlW7xDAaDQZwyZYo4duxY8dy5c2JFRYX9o7W1NeDLL4rWH9bHjh0Tt2/fLgIQDx48KB47dky8fv26rOVX8hl+/vOfi4mJieJnn30mFhQUiHfffbc4fvx42f8NPMGAVdhjjz0mxsTEiMHBweK4ceM6HXJy++23i48++qjT1y5cuCACEPfv39/p+ygZsEo8Q1lZmXjfffeJcXFxYlBQkJiYmCg++uij4pkzZ7rNM9hqfc4+Lly4EPDlF0VRfOGFF5yWX4nfIpR6hhs3boj/8R//IcbExIihoaHiAw88IJaUlMhefk+wDZaISCEcRUBEpBAGLBGRQhiwREQKYcASESmEAUtEpBAGLBGRQhiwREQKYcASESmEAUtEpBAGLBGRQhiwREQKYcASESnk/wPyEcS1DlFPaAAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Isolate Sheyenne National Grasslands from the grasslands_gdf\n",
    "sheyenne_grasslands_gdf = grasslands_gdf[grasslands_gdf['GRASSLANDN']=='Sheyenne National Grassland']\n",
    "display(sheyenne_grasslands_gdf)\n",
    "sheyenne_grasslands_gdf.plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "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>OBJECTID</th>\n",
       "      <th>NATIONALGR</th>\n",
       "      <th>GRASSLANDN</th>\n",
       "      <th>GIS_ACRES</th>\n",
       "      <th>SHAPE.AREA</th>\n",
       "      <th>SHAPE.LEN</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>187254</td>\n",
       "      <td>295515010328</td>\n",
       "      <td>Caddo National Grassland</td>\n",
       "      <td>68479.549</td>\n",
       "      <td>0.02694</td>\n",
       "      <td>1.159342</td>\n",
       "      <td>MULTIPOLYGON (((-95.85493 33.79814, -95.85495 ...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   OBJECTID    NATIONALGR                GRASSLANDN  GIS_ACRES  SHAPE.AREA  \\\n",
       "9    187254  295515010328  Caddo National Grassland  68479.549     0.02694   \n",
       "\n",
       "   SHAPE.LEN                                           geometry  \n",
       "9   1.159342  MULTIPOLYGON (((-95.85493 33.79814, -95.85495 ...  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "<Axes: >"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAASUAAAGdCAYAAABU/mVvAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAArdElEQVR4nO3de3xU1aEv8N/MZPKaZCZvCGRU8gIJJmKLEh/1QUDxAbbnQrkq6Adaj15pscdzrkK1gufwuM2nQC0FAbGktTWtUlpEifIQBSyhQFIDEiCChJAhkJDMTCaZ97p/ACMxr5lkHiuZ3/fzmQ/M7D1r1t4f9o+11157bYUQQoCISBLKUFeAiOhaDCUikgpDiYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpRIS6Av7idrtRX1+P+Ph4KBSKUFeHiK4hhIDZbMawYcOgVPbcFho0oVRfXw+9Xh/qahBRD86ePYuMjIwe1xk0oRQfHw/g8kZrtdoQ14aIrmUymaDX6z3HaU8GTShdPWXTarUMJSJJedO1wo5uIpIKQ4mIpMJQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqDCUikgpDiYikwlAiIqn4FEpr1qxBfn6+5/6ywsJCbNu2zbN84cKFGDVqFDQaDRITE1FUVITy8vJey125ciVGjhyJmJgY6PV6/OxnP4PVavV9a4howPPphtyMjAwsW7YM2dnZAICSkhJMnToVFRUVyMvLQ25uLlatWoXMzEy0t7djxYoVmDRpEmpqapCamtplmX/84x/x0ksv4a233sLtt9+OEydO4KmnngIArFixon9bR0QDjqK/T8hNSkpCcXEx5syZ02mZyWSCTqfDjh07MGHChC6/P3fuXBw7dgw7d+70fPbCCy/gwIED2LNnj9f1uPpbRqOxX7MEHK5thj4xFqnxUX0ug4g68uX47HOfksvlQmlpKSwWCwoLCzstt9vtWLduHXQ6HQoKCrot584778ShQ4dw4MABAMCpU6fw4Ycf4qGHHurx9202G0wmU4dXf2z/sgE/WL0PP1j9OS5Z7P0qi4j6zuf5lKqqqlBYWAir1Yq4uDhs3rwZo0eP9izfunUrZsyYgba2NqSnp2P79u1ISUnptrwZM2bg4sWLuPPOOyGEgNPpxLPPPouXXnqpx3osXboUixYt8rX63Sr5/Gscrm0BAESoOJ0uDXwnG8zYWX0BbiGgViqhUioQoVJAoVBArVQgMkKJ2MgIaGMiMEQbDW20GpooFaIjVFAqQ3cM+Hz6ZrfbUVtbi5aWFmzatAlvvvkmPv30U08wWSwWGAwGNDY2Yv369di1axfKy8uRlpbWZXm7d+/GjBkz8D//8z+47bbbUFNTg3nz5uHHP/4xXnnllW7rYbPZYLPZPO+vzmzX19O39Z+dwuIPjwEAfvlv+Zg+jlPr0sBjsTmx7ch5/Pmftfjn1819LidGrUJspAoxkVf/jEDstz6LjYzA3SNTce/Iro/ta/ly+tbvPqWioiJkZWVh7dq1XS7PycnB7NmzMX/+/C6X33XXXRg/fjyKi4s9n7399tt4+umn0dra2usk41f1t0/p0Jlm/NuazwEAEUoF1j/5Xa92NlGoWR0uHK03YusXBmyprEdTELsf/s89Wfi/D4zqdT1fjs9+T4crhOjQYvF1eVtbW6fgUalUEEKgn3npE9U1zVWnW+DZtw/hjz8aj+9cnxi0OhB153LwmKBQAHFREYhUKaGLUWP/qSa88vejaGzt/hgLpEA8OMinUFqwYAEmT54MvV4Ps9mM0tJS7N69G2VlZbBYLFi8eDGmTJmC9PR0NDU1YfXq1airq8O0adM8ZcyaNQvDhw/H0qVLAQCPPPIIli9fjrFjx3pO31555RVMmTIFKpXKv1vbg437Tnd4b3W4MXvjP/GXfy/EyKG9T3ZOoed2CzRZ7LhgtsLQYoXB2I5zLVacbmzFuZZ2nDfa4HC5MSwhBiOHxCF3aDxuGq7DTcN1SIiN7LLMs5fasGHvaXx1sRVNrXbkDInDo2OH446sFERG9G/s8bmWdrzytyNQKhRQKS//xxihVF7p67l8muRyCaiUCrx3qC6oLSBvKeD/VPIplBoaGjBz5kwYDAbodDrk5+ejrKwMEydOhNVqRXV1NUpKStDY2Ijk5GSMGzcOe/bsQV5enqeM2traDi2jl19+GQqFAi+//DLOnTuH1NRUPPLII1i8eLH/trIXja02fHS0odPnxnYHZr1VjhXTb8atI5IQoeIAeBmdbrTgt5/U4MMqA9rsrl7XN7Y7cMzQ8WptSlwUNFEqxKhViFarEKlSwi0EDtc2w31Ng/1Lgwl/r6yHJlKFu0emYvp39bg7N7VPzxq0OlzYVX3B5+/JJBD94f3uU5JFf/qUVu44gZU7Tva4zp7/ey/0SbH9qSL5WVWdEW/uPYX3/1XfITiCrSBDh+fuzcbt2SmIi/L+//kTDWZMWvFZAGsWeM8X5eD5otxe1wtqn9JgUJCR0OPyIdooZCTGBKcy1COXW+CDKgPe2P0VvjT0b2yav/yrzoin/3AIABCtViIuKgIqpQIutwCgQMSVS/GREUqolUpEqy+fotmd7tBWXFIMJfTeWTdyqJaPAg8yu9ONXdUN2PqFAV/Wm3C1IWS2OtDYKl/fylVWhxtWh7z18zdlAI6LsA8lp8uNX5Yd73GdWHXwOtzDndst8LfKc/hl2XGcN/GmbNmpAtCpFPah5Bbo9SqKup9XWch7v9hyBG/vrw11NchLgTiBCPujLTJCiXRddM/r8Kpb0Pz18LlQV4F8EIghAWF/tO2racS2I+d7XKe/41HIe1ZH75f0SR6B+P867I+28lNNva6j5g26QSGECOmlffIdW0p+5nS58acDvfdfqHn6FhTeDHwkuYT8NpPB5rzJ6tXlZZ6+9U+73YU/7P8afz18DsZ2ByJUCqhVSkSqlIiKuDylhlqlhLHdEeqqko84JMDP9tU0erUeO7q7Z2x3YO2nX+FwbTPsTjciI5SIUCoRoVJAqVDA7nTjaL0RzW0MnMEoELeZhG0oCSGwYe/p3lcEsLO6ARfMNkzKG8LpTK7x0dHz+MXfj6DBFJo71Cn0AjEZXNiG0gWzDScaWr1a98g5E46cMyE7LQ73jgxwxSR1psmCg183w+Z0w+FyY8/Ji9hxbGDfTEr9F4g7HcI2lCquTH3rC210+OwuIQRqLrRib00jdlVfwJ6T3p3qUngJxHXp8DnKrnHBbMXGz707dbuWLkYNADBZHdh17AKaLHbMuXMEAGDJh8cQFaHEY7ddh3TdwL559/OvGjGvtBIXzTwto56xo9tPPjpyHvtPXfL5e3trGvGnA7XYV9MIh0vgofx0zMHlUDrZYMYnxy8iIzEGPxx3nb+rHFQXzTYGEnmFQwL8pM3uQkpclM9TiP7+H2c6vBdC4ILJiuhIFayOy9NQBGIwWbBxVDV5S8WWkn/8+91ZmP5dPV7dchRb/lXf53I+rDqPD6t6vkVlIGrnIEbyEltKfpSoicTr/3ssHhgzFK/vPAmFQoFotRLVBjPafWgpDE+IwbwJObA6XWi3u5A3vO9P55WFjZOPkZd49S0AHrwpHQ/elO55/9yfDuODLwxefz81PmrQPSPO4WIokXd4Q24Q3DYiyaf1B+MtKGwpkbd4Q24Q3OpjKEUxlCiMcZK3IMhNi/cMkoyKUGJ0es99RIPxvjhefSNvsU8pCJRKBR7KH4bU+Cj89L5sKBUK3LpkZ7fDBwbj6RuvvpG3OKI7SJb+4KYO79fP+g62fmHo8gbewRhKVp6+kZcCMSff4DuiAmDsdYnd9h0NxgngbDx9oxAafEdUgJisnecD+sHY4Xjq9huCX5kAY0c3eSsQD9jm6VsP3G6BkxdacehMM8q7uFduWEIMxgzXhaBmgcU+JQolhlI3/vCPr/HLj47DbHV2u87h2uYg1ih4fBnRTuRvPH3rhiYqosdAykiMQXZaXBBrFDwMJfJWAM7e2FLqToE+ocN7bXQEJo9Jx22ZSbgtMxnDEwb2nEk94TglCiWGUjdGJGsQHxUBs+1ya0kXq8b/+1/5Ia5VcNjZ0U0hxNO3biiVCuTrv+nEPnupHfUt7SGsUfDwgZDkLRGAkUoMpR4UZCR0eL/fi6fpDgaBeGwOkbcYSj24PSulw/uj9aYQ1SS4XGwqkZcCMcsN+5R6UJiVjLLn74JapcTwhBhEq1WhrlLAud0CzW29PzWYCABsTv9fFGEo9UClVGDU0IE/k6QvTFYH+5TIa4H4t8LTN+qgp7FZRN+mDkAHJEOJOmjjLSbkg0DckM5Qog66uvGYqDsRKraUKMACcdsADV6BmHmSoUQdON0czU3eC8SYNoYSdeBwsalE3lOypUSBxrmUyBdsKVHAGds5cJK8xz4lCrj6Fmuoq0ADCE/fKODa7Bw8Sd4LxL3bDCXqgHMpkS+UAUgQhhJ1wCeZkC/Yp0QBx1AiX7BPiQIuEFNR0OClYihRoPGGXPIFxylRwHHwJPmCfUoUcHy8EvkiAJnEUKKOePpGvgjErBIMJeqAoUShxlCiDiwc0U0hxlCiDtjRTb7gwygpoFxuwcGTFHIMJfLgqRvJgKFEHjx1I5/x6hsFEkOJfBWIyZMZSuTB/iSSAUOJPDiam3zFwZMUUBw4STJgKJEHW0rkK45TooBqZyiRBCJCXQGSByd4o6HaaBRmJaPV5oTV4UK73YU2uwvtnr870e5weR5aGhGASboZSuRhc/DqW7i7WZ+AFT+8udf1bM7LwRSjVvm9Dgwl8mCfEsVFexcJUREqRAUoPXxqe61Zswb5+fnQarXQarUoLCzEtm3bPMsXLlyIUaNGQaPRIDExEUVFRSgvL++xzHvuuQcKhaLT66GHHurbFlGftbOlFPZiI/3f8vGVT1mXkZGBZcuWITs7GwBQUlKCqVOnoqKiAnl5ecjNzcWqVauQmZmJ9vZ2rFixApMmTUJNTQ1SU1O7LPOvf/0r7PZvHhXd1NSEgoICTJs2rR+bRX3RznvfBrxhumgAQHy0GgLiSl+QG3anq0NfUHdiJAglhRD9G/6UlJSE4uJizJkzp9Myk8kEnU6HHTt2YMKECV6Vt3LlSvziF7+AwWCARqPxuh5Xf8toNEKr1Xr9PfrG4g++xPo9p0NdDeqjSaOHYPkPb0ZcD+dVDpfb04Ftc7phc17pyL7SmZ2RGIPstHi/182X47PPZ4UulwvvvvsuLBYLCgsLOy232+1Yt24ddDodCgoKvC53w4YNmDFjhk+BRP7B20wGppv1CVjw4I347vWJUPbyeBG1Sgm1Son4aHWQauc7n0OpqqoKhYWFsFqtiIuLw+bNmzF69GjP8q1bt2LGjBloa2tDeno6tm/fjpSUFK/KPnDgAI4cOYINGzb0uq7NZoPNZvO8N5lMvm4KfQuvvg0Mo4bG4/rkWLS0OXDMYMLamd/BEG10qKvlNz6H0siRI1FZWYmWlhZs2rQJTz75JD799FNPMN17772orKxEY2Mj1q9fj+nTp6O8vBxpaWm9lr1hwwaMGTMGt956a6/rLl26FIsWLfK1+tQDjlOS360jkrB+1nehi5G3pdNf/e5TKioqQlZWFtauXdvl8pycHMyePRvz58/vsZyrLavXXnsN8+bN6/V3u2op6fV69in1w7NvH8K2I+dDXQ3qxgN5Q7Fyxs2IDsDYoEALSp/SVUKIDuHg6/Kr/vKXv8Bms+GJJ57w6nejoqIQFRXldT2pd/YrfUpxURF4fPx1KMxMhsDl0zqb0wWbww2r03Wlo9SNdsflv9ucbsRGqmB3uuFwuWFzumG2OtHSZkdzmx0tbQ60tDvgcgdi9p3BTaEAYtQq/OCW4Vg0ZQxUgXgkrWR8CqUFCxZg8uTJ0Ov1MJvNKC0txe7du1FWVgaLxYLFixdjypQpSE9PR1NTE1avXo26uroOl/dnzZqF4cOHY+nSpR3K3rBhAx599FEkJyf7Z8vIJy1tdqiUCvx0Qg5mFV6PlDj/Br7bLdDS7kB9SzsaTFbUG61osznRZnd5Bm1qY9SwO91obrMjJlKFVqsTLW0ONLfZ0dzmQEubHZcs9rDqkB8/IhnvPD0+1NUIKp9CqaGhATNnzoTBYIBOp0N+fj7KysowceJEWK1WVFdXo6SkBI2NjUhOTsa4ceOwZ88e5OXlecqora2F8lv3y5w4cQJ79+7Fxx9/7J+tIp8lxEZi3azvBqx8pVKBJE0kkjSRGDNc16+y2u0uXGqzw3LN/VlWp/vyn1dab+1XXlbH5UvgbXan5++e5XZXp8/arlwql0WEavC3jL6t331KsuA4JfIXIS4/1cXqcMHudF8Zz+OG3emG3XXlzytjfOxON6rOGbF691cBqcv9eUOwdmbg/rMIlqD2KRENNgqFAtFqldcdygX6BAxLiPEEV+2lNmw7YoDV4YLbDbw0eRRsTjfcQuDdg2fxdVNbj+UlxKrx1lPjAABJsZH93p6Bhi0logASQkCh+OYU7LzRimfePoSj9UY43QIjkjVI1ERi/uRRyEiMRWyUCi6XQKJmcIURW0pEkrg2kABgqC4af3vuDrjcAnanW4p7zWTDmSeJQkClVDCQusFQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqDCUikgpDiYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqDCUikgpDiYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqDCUikgpDiYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqDCUikgpDiYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqDCUikgpDiYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqPoXSmjVrkJ+fD61WC61Wi8LCQmzbts2zfOHChRg1ahQ0Gg0SExNRVFSE8vLyXsttaWnBc889h/T0dERHR+PGG2/Ehx9+6PvWENGAF+HLyhkZGVi2bBmys7MBACUlJZg6dSoqKiqQl5eH3NxcrFq1CpmZmWhvb8eKFSswadIk1NTUIDU1tcsy7XY7Jk6ciLS0NLz33nvIyMjA2bNnER8f3/+tI6IBRyGEEP0pICkpCcXFxZgzZ06nZSaTCTqdDjt27MCECRO6/P4bb7yB4uJiVFdXQ61W97keV3/LaDRCq9X2uRwi8j9fjs8+9ym5XC6UlpbCYrGgsLCw03K73Y5169ZBp9OhoKCg23K2bNmCwsJCPPfccxgyZAjGjBmDJUuWwOVy9fj7NpsNJpOpw4uIBgHhoy+++EJoNBqhUqmETqcTH3zwQYfl77//vtBoNEKhUIhhw4aJAwcO9FjeyJEjRVRUlJg9e7Y4ePCgeOedd0RSUpJYtGhRj9979dVXBYBOL6PR6OsmEVGAGY1Gr49Pn0/f7HY7amtr0dLSgk2bNuHNN9/Ep59+itGjRwMALBYLDAYDGhsbsX79euzatQvl5eVIS0vrsrzc3FxYrVacPn0aKpUKALB8+XIUFxfDYDB0Ww+bzQabzeZ5bzKZoNfrefpGJCFfTt/63adUVFSErKwsrF27tsvlOTk5mD17NubPn9/l8rvvvhtqtRo7duzwfLZt2zY8+OCDsNlsiIyM9Koe7FMikldQ+pSuEkJ0aLH4uvyOO+5ATU0N3G6357MTJ04gPT3d60AiosHDp1BasGAB9uzZg6+//hpVVVX4+c9/jt27d+Pxxx+HxWLBggULsH//fpw5cwaHDx/Gj370I9TV1WHatGmeMmbNmtWh1fTss8+iqakJ8+bNw4kTJ/DBBx9gyZIleO655/y3lUQ0YPg0TqmhoQEzZ86EwWCATqdDfn4+ysrKMHHiRFitVlRXV6OkpASNjY1ITk7GuHHjsGfPHuTl5XnKqK2thVL5TRbq9Xp8/PHH+NnPfob8/HwMHz4c8+bNw4svvui/rSSiAaPffUqyYJ8SkbyC2qdERORPDCUikgpDiYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqDCUikgpDiYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqDCUikgpDiYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqDCUikgpDiYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqDCUikgpDiYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqDCUikgpDiYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpMJSISCoMJSKSik+htGbNGuTn50Or1UKr1aKwsBDbtm3zLF+4cCFGjRoFjUaDxMREFBUVoby8vMcyN27cCIVC0elltVr7tkVENKD5FEoZGRlYtmwZDh48iIMHD+K+++7D1KlTcfToUQBAbm4uVq1ahaqqKuzduxc33HADJk2ahIsXL/ZYrlarhcFg6PCKjo7u+1YR0YClEEKI/hSQlJSE4uJizJkzp9Myk8kEnU6HHTt2YMKECV1+f+PGjXj++efR0tLSn2p4fstoNEKr1farLCLyL1+Ozz73KblcLpSWlsJisaCwsLDTcrvdjnXr1kGn06GgoKDHslpbW3H99dcjIyMDDz/8MCoqKnr9fZvNBpPJ1OFFRAOfz6FUVVWFuLg4REVF4ZlnnsHmzZsxevRoz/KtW7ciLi4O0dHRWLFiBbZv346UlJRuyxs1ahQ2btyILVu24J133kF0dDTuuOMOnDx5ssd6LF26FDqdzvPS6/W+bgoRScjn0ze73Y7a2lq0tLRg06ZNePPNN/Hpp596gsliscBgMKCxsRHr16/Hrl27UF5ejrS0NK/Kd7vduOWWW/C9730Pr7/+erfr2Ww22Gw2z3uTyQS9Xs/TNyIJ+XL61u8+paKiImRlZWHt2rVdLs/JycHs2bMxf/58r8v88Y9/jLq6ug5X9nrDPiUieQWlT+kqIUSHFouvy7tav7KyEunp6f2tGhENQBG+rLxgwQJMnjwZer0eZrMZpaWl2L17N8rKymCxWLB48WJMmTIF6enpaGpqwurVq1FXV4dp06Z5ypg1axaGDx+OpUuXAgAWLVqE8ePHIycnByaTCa+//joqKyvx29/+1r9bSkQDgk+h1NDQgJkzZ8JgMECn0yE/Px9lZWWYOHEirFYrqqurUVJSgsbGRiQnJ2PcuHHYs2cP8vLyPGXU1tZCqfymgdbS0oKnn34a58+fh06nw9ixY/HZZ5/h1ltv9d9WEtGA0e8+JVmwT4lIXkHtUyIi8ieGEhFJhaFERFJhKBGRVBhKRCQVhhIRSYWhRERSYSgRkVQYSkQkFYYSEUmFoUREUmEoEZFUGEpEJBWGEhFJhaFERFJhKBGRVBhKRCQVhhIRSYWhRERSYSgRkVQYSkQkFYYSEUmFoUREUmEoEZFUGEpEJBWGEhFJhaFERFJhKBGRVBhKJK12uwtf1ptgbHd4PjNbHbhkscPlFiGsGQVSRKgrQHQtl1vgk+oL+NOBWuw5eREO1+XwSYuPQmJsJE5cMEMIQBOpwuhhWpy6aIEuRo2RQ+Nhc7pR19yG2MgIREUoERmhhFqlRIRSAbcAUuOjEK1WXv5+lApCAC4hEB2hQkKsGrdcl4j8DB0UCkWI90J4YyhRyAkhcPBMM/aebMTminOovdTWaZ0LZhsumG2e9xa7C//8uhkA0GSx41SjxS91GZGiwQ/H6TGlYBiGJcT4pUzyjUIIMSjawSaTCTqdDkajEVqtNtTVIS+caDDj9//4GntONuJMU+cgCrWbhutwsz4B94xMxYQbh4S6OgOaL8cnW0oUVMcMJpQeqEX56UuoPm8OdXV6VHXOiKpzRmiiIhhKQcRQooBrt7vw8ZfnsenwOXx24mKoq+Mzm9MV6iqEFYYSBYTV4cLfKs7hgyoDDn7djHbHwD2w46N4mAQT9zb5jdXhwqEzzdh57AL+WlGHljZH718aAOKj1aGuQlhhKFG/1TW34U/ltfhjeW2HMUWDRUIsQymYGErkM7dboOZiKz4+eh77apqw/3QTBsc13K5pYxhKwcRQol653ALHDCYcOtOMj46eR+XZFrTZB24fka/YpxRc3Nvk4XS5YbG7oItRo7apDeWnm7Dz2AXs+6oRZqsz1NULmehIVairEFYYSgS3W+CDKgN++0kNTjdacEOyBjUXW3l/2RWRKt4iGkwMJULNxVb85J0Kz/vjDXIPagy2WLaUgor/BRDsTneoqyC1GIZSUDGUCDaGUo+iIxhKwcRQIjhdDKWesGctuNinFAZqLrSi5oIZsZERUKuUiFYrEa1WIVqtQlSEEk0We6irKLXzRiuSNJGhrkbYYCiFgV3VDVjyYXWoqzFg1be0Y/QwTocTLDx9CwNOXtrvl1ONraGuQlhhKIUBp4uh1B+fnWgMdRXCCkMpDLCl1D/7TzXBOEhmPBgIGEphwOXm1bX+cLoFjtYbQ12NsMFQCgNsKfUf92Dw8OpbGPheTipcLoG//6seF695Igh5bzBPzSIbhlIYuCM7BXdkp+CLOiNDqY8E20pBw9O3MMIJ8PuOLaXgYSiFEQeHBvQZ91zwMJTCiIP3uPXZIHlm64DAUAoj4TSFrb8xkoKHoRRG/vij27Dzhbuhjeb1DV9xVHzwMJTCyA0pGkQoFTCF8XzbfWW2ckR3sDCUwkz5qUuhrsKANBifZycrhlKYqTjbHOoqDEgWG1uXwcJQCjMVtS2hrsKAZGYoBQ17PMPMmie+A4vNCZvTjSPnjHh1y9FQV2lACOfn3gUbQynMjEjReP5udXCIgLdaGUpBw9O3MMZQ8l4rT9+ChqEUxvhoJe9xNHzwMJTCGA807/ER5sHDUApjvEHXewyl4GEohTE+hNJ7bt6QGzQ+hdKaNWuQn58PrVYLrVaLwsJCbNu2zbN84cKFGDVqFDQaDRITE1FUVITy8nKvyy8tLYVCocCjjz7qS7Wojx4YMxSbnr0dJbNvxW8fuwVLf3ATXn7oRuTxGWedsKUUPD4NCcjIyMCyZcuQnZ0NACgpKcHUqVNRUVGBvLw85ObmYtWqVcjMzER7eztWrFiBSZMmoaamBqmpqT2WfebMGfznf/4n7rrrrr5vDfkkITYS37m+85NfTzSYcbTeFIIayYtnusHjU0vpkUcewYMPPojc3Fzk5uZi8eLFiIuLw/79+wEAjz32GIqKipCZmYm8vDwsX74cJpMJX3zxRY/lulwuPP7441i0aBEyMzP7vjXkF0qFItRVkI6bLaWg6XOfksvlQmlpKSwWCwoLCzstt9vtWLduHXQ6HQoKCnos67XXXkNqairmzJnj9e/bbDaYTKYOL/IPZlJnfCJM8Pg8oruqqgqFhYWwWq2Ii4vD5s2bMXr0aM/yrVu3YsaMGWhra0N6ejq2b9+OlJSUbsvbt28fNmzYgMrKSp/qsXTpUixatMjX6pNXmErfxosCweNzS2nkyJGorKzE/v378eyzz+LJJ5/El19+6Vl+7733orKyEp9//jkeeOABTJ8+HRcuXOiyLLPZjCeeeALr16/vMbi6Mn/+fBiNRs/r7Nmzvm4KdUPJTOqEHd3BoxD9nHy4qKgIWVlZWLt2bZfLc3JyMHv2bMyfP7/TssrKSowdOxYqlcrzmfvK01yVSiWOHz+OrKwsr+phMpmg0+lgNBqh1fLqUX8s2FyFP5XXhroaUoiPikBmqga3XJ+IVx/JC3V1Bixfjs9+35ArhIDN1v2zxHpaPmrUKFRVVXX47OWXX4bZbMavf/1r6PX6/laP+kAVZp1KKqUC1yXFYkSKBiNSNMhKjUNmqgaZqRqkxkVBEWb7I9R8CqUFCxZg8uTJ0Ov1MJvNKC0txe7du1FWVgaLxYLFixdjypQpSE9PR1NTE1avXo26ujpMmzbNU8asWbMwfPhwLF26FNHR0RgzZkyH30hISACATp9T8AzWYzAhVo3sK4FzOXgu/12fGIvICI4jloVPodTQ0ICZM2fCYDBAp9MhPz8fZWVlmDhxIqxWK6qrq1FSUoLGxkYkJydj3Lhx2LNnD/Lyvmn21tbWQqnkPwCZDfRMykiMQXZaHLJSL78u/12D5LioUFeNvNDvPiVZsE/Jfypqm3HMYIZLCNidbrjcbjhcAk6XgNPtht3lhsMp4BYCDpcbdqcbNRdbgzqrZYRSgeuTY5GdFofcIfHITrscPpkpcYiJVPVeAAVVUPuUaPAZe10ixl6X6NN3/vLPswEJpagIJbJS45AzJA7ZV/9Mi8cNybGIULHFPRgxlMgvbP0cx5MYq0ZOWjyyrpxqXT3tGpYQAxXHKIQVhhL5hTeDC5UKQH/lKtc3rZ7Lp1yJms734FF4YiiRX1w7uFCtUiAzJQ7ZQy53NOdc6XTOTNUgWs3+HuoZQ4n8YtwNSVg78zvITovD9Uns76G+YyiRXxToE9DzbddE3uF/Z0QkFYYSEUmFoUREUmEoEZFUGEpEJBWGEhFJhaFERFJhKBGRVBhKRCQVhhIRSYWhRERSYSgRkVQYSkQkFYYSEUll0ExdcvX5ByaTKcQ1IaJvu3pcevOckkETSmazGQD4AEsiiZnNZuh0uh7XGTSPWHK73aivr0d8fPygfaKpyWSCXq/H2bNnw/4xUtwX3xgI+0IIAbPZjGHDhvX63MdB01JSKpXIyMgIdTWCQqvVSvuPL9i4L74h+77orYV0FTu6iUgqDCUikgpDaQCJiorCq6++iqioqFBXJeS4L74x2PbFoOnoJqLBgS0lIpIKQ4mIpMJQIiKpMJSISCoMJYkcPnwYEydOREJCApKTk/H000+jtbW103obN25Efn4+oqOjMXToUMydO7fHctetW4d77rkHWq0WCoUCLS0tAdoC/wnUvrDZbPjJT36ClJQUaDQaTJkyBXV1dYHaDL/wZl8oFIpOrzfeeKPHcr/66it8//vfR2pqKrRaLaZPn46GhoZAbopXGEqSqK+vR1FREbKzs1FeXo6ysjIcPXoUTz31VIf1li9fjp///Od46aWXcPToUezcuRP3339/j2W3tbXhgQcewIIFCwK4Bf4TyH3x/PPPY/PmzSgtLcXevXvR2tqKhx9+GC6XK4Bb1Hfe7gsA+N3vfgeDweB5Pfnkk92Wa7FYMGnSJCgUCuzatQv79u2D3W7HI488ArfbHcAt8oIgKaxdu1akpaUJl8vl+ayiokIAECdPnhRCCHHp0iURExMjduzY0aff+OSTTwQA0dzc7I8qB0yg9kVLS4tQq9WitLTU89m5c+eEUqkUZWVl/tsAP/JmXwghBACxefNmr8v96KOPhFKpFEaj0fPZpUuXBACxfft2v9S9r9hSkoTNZkNkZGSHmxVjYmIAAHv37gUAbN++HW63G+fOncONN96IjIwMTJ8+HWfPng1JnQMlUPvi0KFDcDgcmDRpkuezYcOGYcyYMfj8888DtDX9482+uGru3LlISUnBuHHj8MYbb/TY4rHZbFAoFB0GXEZHR0OpVHYqN9gYSpK47777cP78eRQXF8Nut6O5udlzumUwGAAAp06dgtvtxpIlS7By5Uq89957uHTpEiZOnAi73R7K6vtVoPbF+fPnERkZicTExA6fDxkyBOfPnw/sRvWRN/sCAP77v/8b7777Lnbs2IEZM2bghRdewJIlS7otd/z48dBoNHjxxRfR1tYGi8WC//qv/4Lb7e5QbigwlAJs4cKFXXZCXvs6ePAg8vLyUFJSgl/96leIjY3F0KFDkZmZiSFDhkClUgG4PD2Lw+HA66+/jvvvvx/jx4/HO++8g5MnT+KTTz4J8Zb2TtZ9IYQI+nQ3/twXAPDyyy+jsLAQN998M1544QW89tprKC4u7vb3U1NT8e677+L9999HXFwcdDodjEYjbrnllg7lhsKgmbpEVnPnzsWMGTN6XOeGG24AADz22GN47LHH0NDQAI1GA4VCgeXLl2PEiBEAgPT0dADA6NGjPd9NTU1FSkoKamtrA7MBfhTqfTF06FBPa+Pa1tKFCxdw++2392fTfObPfdGV8ePHw2QyoaGhAUOGDOlynUmTJuGrr75CY2MjIiIikJCQgKFDh/ZYblCEtEeLerRhwwYRGxvr6Zg+fvy4ANChc7epqUkolUrx0Ucf9VreQOno7oo/9sXVju4///nPns/q6+ul7ujuyrf3RVd+85vfiOjoaGG1Wr0ud+fOnUKhUIjq6mo/1LLvGEoS+c1vfiMOHTokjh8/LlatWiViYmLEr3/96w7rTJ06VeTl5Yl9+/aJqqoq8fDDD4vRo0cLu90uhBCirq5OjBw5UpSXl3u+YzAYREVFhVi/fr0AID777DNRUVEhmpqagrp9vgjUvnjmmWdERkaG2LFjhzh8+LC47777REFBgXA6nUHdPl/0ti+2bNki1q1bJ6qqqkRNTY1Yv3690Gq14qc//alnna72xVtvvSX+8Y9/iJqaGvGHP/xBJCUlif/4j/8I6rZ1haEkkZkzZ4qkpCQRGRkp8vPzxe9///tO6xiNRjF79myRkJAgkpKSxPe//31RW1vrWX769GkBQHzyySeez1599VUBoNPrd7/7XRC2qm8CtS/a29vF3LlzRVJSkoiJiREPP/xwh+/IqLd9sW3bNnHzzTeLuLg4ERsbK8aMGSNWrlwpHA6HZ52u9sWLL74ohgwZItRqtcjJyRG/+tWvhNvtDtZmdYtTlxCRVHj1jYikwlAiIqkwlIhIKgwlIpIKQ4mIpMJQIiKpMJSISCoMJSKSCkOJiKTCUCIiqTCUiEgqDCUiksr/B73blNHSrS6cAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Isolate Caddo National Grasslands from the grasslands_gdf\n",
    "caddo_grasslands_gdf = grasslands_gdf[grasslands_gdf['GRASSLANDN']=='Caddo National Grassland']\n",
    "display(caddo_grasslands_gdf)\n",
    "caddo_grasslands_gdf.plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Stored 'sheyenne_grasslands_gdf' (GeoDataFrame)\n",
      "Stored 'caddo_grasslands_gdf' (GeoDataFrame)\n",
      "Stored 'data_dir' (str)\n"
     ]
    }
   ],
   "source": [
    "%store sheyenne_grasslands_gdf caddo_grasslands_gdf data_dir"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "earth-analytics-python",
   "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.11.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}