vs_raster_rioxry_merge_plot.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "a76c58d0-b80c-40e3-9618-35b5e73fd5e1",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAD8CAYAAABetbkgAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABaoUlEQVR4nO2dd3gc1dX/P3dmd9W7ZBVLttzkSrExGHAILdTQCYQkJCQh4Zc39c2bBmmkh7xvGkkIqSSk0AKhGELvxYC7cbdsy7JcJVuyZJXV7sz9/bE7q9nZme2SJWu+z7OPtDN37j0zO+eec8859xwhpcSFCxfHNpSjTYALFy6GHy6ju3AxDuAyugsX4wAuo7twMQ7gMroLF+MALqO7cDEOMOKMLoS4UAixWQjRLIS4eaTHd+FiPEKMpB9dCKECW4DzgDZgGfABKeWGESPChYtxiJGW6KcAzVLK7VLKQeA+4PIRpsGFi3EHzwiPNxHYZfreBiyyNhJC3ATcBKCinpRP8chQF0awqoCa6kMI07FM9B5h+R6UCh17S1EP9ca29XgYrMyhoLSfnt48cjs0BotUvL069PYnHktVkdMUqnxHADgUKEBv1pGanh7tPi9MgSpvD7t6yslt8yM1DQRo5QVIFTwHehGKwmBVHmUVPeQpAYTNE5MI2npLyW0dROb5Qg81fE/BqgKkAr6OAYJTPFTndsfSktYd2NyTA212xzOBuU9r3xKBtNxRNkbfsa6vQ0pZZT0+0oxu91vF3J+U8g/AHwCKRblcJM4dbrqGxj7tBPJv28cHa9/MWp+q5RZ/svUCZl7VihT+2LYzm9j2wQpOPmcjy5+dQ+MPliPbB0Mnk3jT1aoJ1P69j4vL16Ih+OnW8ym/qhWpx46VFAKw76LT+cZn/kmPlsddO0+n5H88aBu2wCGG6JKg9Bex/UPzuOiiZZxauI1cJYCCjo6CJkPK487BSv7w0IVMu3Mbg9NrQQiU19eizTsB9eVVqBOq2HRrPXlFA5w/aROLirZF318W2EEVyU16CvHb6Si2bXQbRVlBj4yrSYWAjGY9zfLj6jI9ZftDTct22h0fadW9DWgwfa8H9owwDXHRflIB19W8bXtOl0rUJxnYvZjvrV+PdsocUNRQv2fMp/U7p6MUFKBt2ELjt96k4wsNNP5oBTIwmBL92uRq2v2FbByo47t/+RD6I5XIwdT6MCA8oZdx4jPtrOidQoHi58zqZvob7DUsvaeHKd9fyepvzucby65Ak0rkpTdecq/QOO2Cd9j6+am0z89HCWio0yajvrQSpESbUoMQoWe2s688LbrjwcxsiWBlWOO7jhLzv/njBE0qkWdinWyyMYHFw0gz+jJghhBiihDCB1wHPDbCNMSFroakTjxmVoSOkqRUsMOUnAPs/eIgnoY6ALz7u8nbJ5H+sNSVErnsnaHvqeDtdwi+T+ON646n/idvUfHHpZDA4Cq8PtSmaZHvakU5nsZJ+N8zH5GTQ9fxFfRrXn608ULWXD2VnJfecexL+v3kPrcGj1ejS8uPvNwGqr1dVPqOsPDMTfTXSI405DEwuQwAT0M9IqBR928fvFGa+r0nCYMeM10GneYPEMPAyTCy3SRiNxEkq1lkAyOqukspg0KIzwJPAypwl5Ry/UjSkAi1d7zNkt3nUP65nVxfuzTqXKrMbZ2lNQRdWgF3/P4KJt29kWBnZ+j4lm1UbdmWtTlda2+H9vbkaCwrg4pSghWFeIqL0Xp6QJcMTK2iZ6KHzptOou6qFgpVP74nSglu35iwz8AZx3HCxBb+719X0viuVk6r2EG19zCPHzie6UXtnFy4g5MKWrjsfav5eum1zPpOKxoge3roW1BHyRdaObO0jVpfV+y9ITKSflYmjCfZzecSMaVdv6rQHfs3+ovXJpsY6TU6Usr/AP8Z6XGThQwGKVyymvYb6zPuy3gpzeuvn619D1N/txJtYCDj/jOBUlDA4KmzOPD5XhbX7WBm/kruaV1IR+dU8vL9XNb4NpXeHm5/9XwG2ivY8+9Gqv/4RlJ957zTypqXmyhqhcrcXqbkHMAnNE4t38GOvkpWHplMdzCP6pxupCrZ9K2pTHqqkZz/LKPg6bXobxTxdtNJtFySx6lnreeSijWRvuMxufGczc88WyqxHTMmWgZENAcbWoyJwAmK0NNep9thxBl9TECXDAx6s9KV8SM/3zWH13ZNpejZAvQRYHK1aRr6zrYwERoyGIw6r/f24n3lHbQTFvLquwUN0w7xxWnPRV4uRejct+8UZt/agnagHUTyL53W3k7jozX4q3J5841ZnP3eTeQqAZpy91Lm6eVAoJjHnltESTPk1goqNmjs+mCQ6uJTKXliPZ3nTePgPIHulXjTUG/NE6vVyKUibaWouZ3BkIm0h3gSP54kNy8drPSZkU1mdxndBmp9LTMrD7C2bxK1vi4OBIp5sPlEKgr7uGHSUorUxG4uM36w4WLqvxagfusm0LXUCRLCdp3taZyE7B9AHulFTKxBtu5GHwwAIHO9bP/uAqpO3E/H4UIm3y7wbGiBynK0bS0hO0BgkNqfv4H4lYenz3o3rR/XOH/GJk4p2k6+8DOl4CAbqppg/wGQKdL99jvkl5bQtLmcDWfXcXLhDhSh8/KhJla93sTkp/x43lhPwXknoOUKah/20XG8Ql57E6WPrCX30Fz2vMvDC2tnc9ziNoqUgbjPPR7DGHBi8nh92U0UThIa4qviGgJNqnHpS+Y+0sGIRsalg5F2rwEoRUWISXWI3n4GJ1fg6exHX7cVpSCftv86jptvvB+VaIPcDv8EOgKFAOz3F1HlO8LWIxPQpaD/+3V4nl8Rf1AhQCiohQWI8lL2nzcRqUDXLIlSPYAWVJC9HiYvkSDhwEIvZ122kn7Ny9r2Wi6btI7HWudxuCefnJwAFzZujLimdKlw6+pLCQZU9EGVmTetRq2egCwpDLnJDBI8HtTaGrZ+poHzz13JxsPV5N4oCe7c5UR1QqgV5ey9awLdPXl86sRX6NNyePrH76b0PxsQ+Xl0n9bIoVkqxTt1Dlzsx7Mzl8YlvXj2dbH7snoKLt7Hnh2VlNZ1c1pdC/U5nUzNORD17Df2T+Svby+mYVIH729YTpUn1gdv4OlDx1Hu6+WkghZ2Dlby2O7j+MLU56Pa6FLJ2B4TOR5m/EwYOBmpbtD7gRkrVkgpF1rPu4yeIjwT6/Deo/HB2jfRpcJf2hazZUsdM/7ux7tjPwBywI8oLED290MwiNZ12LYvpaiIgcWzaL1ARVQPkJc3yMzKA0zOP8Siom2Rl8e8FtwXLMErNErV2GAbcJYKAenhrl2L8X2pAH3NRoTHg/D50Pv6bAhT8UyoRJYUoW1uTucxRdNUWYHweNj8lSmIgMB3WFD/4zdCNMyajtA09IIcDh5XRMXqbrRCH4OlXjy9Gl0zfKh+6DhrkIaHVPLaetn2FS/fP+lRFKEzoHv5/spLmPaRdQifj+BJM2n+kJcfnP0QXhGMMK0uFW557WrmfGcfMtfHgXdXgwD9skPcPOspwJmhzExvbmM+bjx3J4ZPh9F9QkNBJyA9Ca83xncZPYvQz5zPjpugYcIh8j7rTY0ZhEDJz2fgXbORX+rguvpllHtCUWypGI5SfXHag8U88v/eg/LqqtiTihpZUqilJQwsnE7O0k2I3By0g4di26cJJTeXwdPnInSJ+vKqqOWI8PpQ66qRA370rsMopSUIRUEO+Nn7odnUXt2C/O8S9DUb8Uyso/WDjfjLJBNW6hTuOIJcMeS8UcvKaP7KLOYs3k7zwUp62/Mp3uhl4j+3hjwSYehnzOeiO1+m3ncw7bWwVfLbeVriwWniMPryiiHbihF4ZNen0c+HmpbZMrq7Rk8D3vWtaANTmFhwmAMVk5NiOeHx0HPlSex9l6DpuF3MLXqHd5dsimlnXeNlY82mS4UK9QiHvtqLp/E0Sv++NEKT8HgQBfkRhtaP9DJQ6SWvrBSZZaOhPjCA54UVqDOmEjzteBR/EGX7HrTOTmRgMLJEUPLzIRiEokL2Xj2VusfbkE/loDeHXHvB3Xuo+7+hOCvr9Kh1djLlG2/iz89nYn9HaBITAqY1IhprkMtCcQA9k3Ko9HRnZPAyGy8htd/LOm7MdyAgVRShR+wLofdDdbzGCS6jp4HO82Zw6+J/kysC/O47pXh+vADful0woRzpVQkW+uiemkfZ/StB01BKS9h77Uxu+Mx/qPEctp25jRfE6pKB2B8zlfWj+dqvND3L9y5/L2X3+UIRd0JBNNSh72iNtJHBIEX3vUnQrrMsQdu2k+ZvncBlc9fzxOZ5yL0z0fJ1ijd7KNqloQQlRSv34G+spHhnkGBLa+JOLRCqipw9hV0XFaMMQu+0AAXbvEz+2/bIveW3BzmkFVLjsV9apQIzw1uZPxt9G0yfLlxGTxHC6+Pw1UfIV0JRa59ufJE9vyljS18NU/K24xUhFdive3ny+rl0D+TQWHqIL054gFwlYNtnhMmT9P0mepGcZnkNhUBzEcLnRQYGkYFBtK3bE9zxMEDXmH3LPtZPmUddmY/Wy4McN6ONnbVl7O4ogKDCnNUC9aWVpPNqK/n5bPrFXC49aTUfLWqmW8/jN5vPpP52neDefZF2OQf6suqrhuhnn22GzwQuo6eInivm81+zn4o6VuftpK4kFOVm/qE/NfnlyP+p/NjJqn+pvKQD0suP11xIyVYQddVwNBjchODuPXj9fjz+QWY+M0jA46HW30GtqoKmEQymr1Po/f3UvqDyRME8lvTMp/ZlhYlL3kHrjTZgij4/HYFC6rydmd5OfHqyyPAaCmqCzTZ2cBk9BXjqJ1L9ue1Ue7tSus7uB7Zj0mTbpYpdgXL+8acLmPrnd9B7ekjDkz8s0DoORv6PxPVnwOBDnUmK7n+TogdESIUPBmNZQ1HZ9IVKrst72a6HYUE2fst0mBzcnHFJQXg86O86ka0/reDa6uUp7V4DYtrHuzadHXKJ0KflULnWT99Zsxm49JSs9DkmIGVMRKABcdIcrjn9LXxitEx7wwtXoseBWlVF1znTOHCZn0+f8HJExUtXBUvEuNleLxqYmbsX/Rdv0qf5eG7XTPKfz7f3n48TCI+HrR8o4KP5tlu3M0a66vVw4phgdM+UyciubrTO7K21/O89mdKvtnJTdci6biCexXw0Y3ZuyB01afohfnHnudQ+7qXk6Y1o3aEoMuHxOEq/TGHsax+u/lOFmDuDj5z7ytEmY0QxdhldCDo/cir+UkH9I7vQj9hHiqXT78AlJzPx5q1cVrk6YXM7Js/2zqN0odmszMo9R/j+qY8ycIqXv9x4Ou3L56LlSQqmHaZvSyllG6Dy9f2ZWeOFwFNTTcd7ptA9VZC/sAMpBV27Sql/VlL4anNWA3FSpW37NaVcn7t/WLo3nrkh1c2/wdGU8mM6Mk4tLUFOqkX6PMjl67IynqdxEjMe2sOphdtszyejto9WJo8H4yUclCqvHp7Jc6vmUrrWQ+29G1PTlISg/VOnMv+Gdzi9pJkCZSh5hoZCr57DL+69gsk/evuoSHi1uBjPknw+WPPWsI5ztBj9mIyM07oOg0McebronV3NjDybMNE4GCnGTsS81hcrUVs7+ITGuaUbOPfsDQyepXJHzzWRSLpk0HPtIt7/6eeo9x2KGsOgq0Dxc93VL/HowuPo3F3ChKUqSlBSsrknKox1uBCcO4WzK16LoSvbTGj0N1rW6mOa0bMNcfJxzPjOBirUI7bnnTY3DDeSZd5k2hmTgZPByHxcFZIDZwco/1ducnvohaDjREG971Bceqbn7udLM/fDTOCc0LGH9i9g4FNNBKoKEEGJd/s+9M6uxOMKgVJYiJw5GTEQhJ27QdfRe2OXckpBAVs+rXBD2D1qpm80GtCyCZfRTdh5cREfLY2uJWGst0eayVNVvdPp12kM8/FbTvsPP/71e5n8CHh6NXztvcjWPeg9PbbXTlih035ZEVUe+/NOuGLCKtb/s54pOe1oKGzuq+Gdzjpa19cy+fEAuS2H0Ft2Ran7Yv5cdlxVzIwzWjin8g16tFxWdTWws6uMkt8Xkf/GFuTkOpSD3cj8XLZ+YgJfX/jwiDL0cGkMqSLtNboQogH4G1BDKP7+D1LK24UQ5cD9QCPQAlwrpewMX3MLcCOgAZ+XUj6daJyR3L22/bbT+P6V90W+Hw0JPlwMni4MDWBA96Kh0BEo4pWO6Wx/axL1LwbIeW1DlKtOOX4Wx9+9iePz09/DbkWPlstb3VNZ215LZ0sZUx4OcmhWDmd8fBmLiuyNhu3BIp5rn82ZlVvY2V9Jsaefufm7o+4LGNE1dDytIVsahdMaPRNGrwVqpZQrhRBFwArgCuCjwCEp5W3h2mplUsqvCSHmAPcSqtZSBzwHNEkZP3XJSDL6rm+czlc//GCUO22kMNoYPBkcChbyy2cuYtbPdhFs241ncgObf1zBd09aMqzjdmn5FCn9qCJ9Q7IdoxvHR0KNN4+RTRqcGD3tt0tKuVdKuTL8fw+wkVAllsuBu8PN7ibE/ISP3yel9EspdwDNhJg+KxAeD2plRUZ9NP5rPzv8MUUubJGtuGXjM9JQ0TN+mas83Xzv4n9RdH8/W/5wMhPu7+K7Jy0ZdiYpVfsyYnLA8blblzbZ/G0SaQ9tg+U81zmHnYOVvN49I2vjQpbW6EKIRmA+8BZQLaXcC6HJQAgxIdxsImAuf9IWPpYVSE1D77E3oiUD4fGw47pqrstJLvY5XVV+tEjubEgto49LK9dw6Zlroo5nC3aW+2SwJ1BKidof5d4z+kuVvmxKeGtfxrLoh89fTuMjQfI27mXPxEb6JuYz5fvtNHizE2+QMaMLIQqBh4D/llJ2C+EYJZZUOaZwn5Haa7nkJ0eIlI4FD9SqKrSOjriFDPovWsAnr31qWGOfRwuTG8gGPSN1T6mO89z+2RR4/VxbvTzKIGb48q0TQKbjJQO7SWtp93Se3jqbWbe3o20N75XfvYcCj4d76i7gE59ZQoV6JOPJJqO7EUJ4CTH5P6WU/w4f3h9evxvr+APh40mXY5JS/kFKuVBKudBLTiYk4mmcxNYvTUctLXW6CeRpJzDx5q3UeTtHHTOOZ2TyYn9o4pusXjWNZzvnRo5pKOwPlLC1v3rE6EiE+YU7ETvzopJ/hAZVqVrdx5+bT0+ayeO9u2m/1SIkuv8MbJRS/tx06jHghvD/NwCPmo5fJ4TIEUJMAWYA9kXOEmDPl09HOWF25LuSmxtKP2SBkpvL9hvqaXhuEK2rK+a8p6GerbefwvzfruG9lWujwhezvX52J5DUkclvoAqJr7aX1549Lur4xt5aKr3pL/GyBcNGUqD4Ofvc1Yg506POS78f8fpqSn5bzPq+5Fa48SaDTN6+xcCHgXOEEKvDn4uB24DzhBBbgfPC3wmXXnoA2AA8BXwmkcXdCTVv9qEc6EQtK6P/ilPYdMc8/KfPDuU/8/oi7XouOQHf/E5y3tgYnYjQ40E5fhbbflbG9y94MKEr6GgazVykj68c9yyDdQEebT8RAE0Knt88kzJPavsihvt3f3fJFjbfWBopummG75mV3Pv2ooxpSHuNLqV8DedCvrb+MCnlD4EfpjumAc+abUgh2PPReeQe0vHtUcl5Yw1bb1uIOiBo/NabqJWVDBYqTPzqIFpfH2pxMT3vmU1nk8rgiUf4yJy3uT4nufpkmcKdII4OcpUAn1v0Ar9+5T0oJ0jOKtvMlLoOfrr6PD46900m53QcbRKBkCS+ZPEKNi2ah1i6JvqkrjHjbwHWnzKRprx9eIUW1xfvhDG3qUUtK+PweTMp/NdbMca1gUtPQRnU8T29PLoTIdjxo1P5+pUPRfnI02XAVNZsLpMffewJlPKXey+gYK9kwtM7QdPYe+VUPvX5RxMa5UYyou3bqy5l+vf96Ft2gNSHogCFQJk7k95pxbS9R/D/znohEmZsxUea3squH/1oQevspPCBN20t6LlL3o5lcsAzqZ4zz1mbFSbP9FoXI49ytRd/lU7VY1sI7t5DcN9+qv+xjv9decHRJi0Kt85/gto/7+bgh09CaqZVrZTo6zaR9+jbzPjCcv685D3sCZSm1Pe4eGN3XtfAWaWbsrrWthrsnD4ujj5ylQCLF20A1ZQPvaeH6b8I8OTB4+JcObJQ0TmnbCPVH2lBO3N+KBe9FbrG1FtX8PTnzuQ7Ky5lQPcmlQZ6zKnuqUCtrGDPB2dywcfeYH7+TpfxxjHag0U8/plzUF9aGXW876pFXPJdZ1UYRlZ9N97RVw430fb+CQR32Ke7Eh4PyrRGDi6qwl8iOHxcgKIJR1h/+fePvf3o8eCZ3MDBO3P476kP4hWay+TjHFWeHrZ/DJpeDxevCKNgySp+994zuXXxY5Gc/AaO5o6z4wvbeOFTJ1L/YiiwNK+5A615R+S8DAbRNjdTurkZFJXaeTOguRWnHf3H7Nvfem0Dn576UsyP52L84ksnP0v31QuijsnAIDP+HGBp97So40d7W+mBQDF5+wX7Fnl590/eoP2XXtTi4pALzqrS6xpiTztycNC+M45RRldLS2i4uCXjjQ8uji2Uqn1M+K8deGproo5793ax9L757A+UHCXKhmBMME25++hd0E/ugkM05e4jxxN0ZGTh8aAd6oqbmuuYZHTZOJH3VMUWMHTh4qrqFez8yNSoY8EdO6m7YwW/eek84OhLc4OGryx4hjPrQ5V6F1W1oFRVIhSB8PlQq6oiATZKSTFKXm7omAOOPUYXgkPHlVCijt+85S6coQrJoivXok6fEnVc+v0UtKpoMv3U3YrQM96+bA7DLlX7OLkwtC5/euds5IAf4fOhTJqI8HnxTKgMtT3USeCUmfCA15m2jKgabVBUDn7iVM770muOBQ1duDizdBNbP1mDyIneMDXpb9v47faz6NbzMuo/XYa3GoyN7wGp4vME6T5zKjQ10nnyBGRhPl1nNKIcP4uOT57K/J+t4kN1zpltjxn3msjJYffnTuJjH32Kam92M8O6OPbQpeVz598upf6n0Wmn1dISui6YzYRP7+CamlDwVaqqfLopyBJ5hn6z/Sw61k7gvHNWoSDJUwcpVP3U+w5FBJtTZNwx414Lnj6Xmz72BOWeo78zycXohy4Vqpf5ETk5SF2CHvLOaF2HKbr/TQa2zGXtXQ0cn78r5b3gw5Vf8LNTX4KpCZvZ4thQ3RWV1vNzXCZ3kTSK1H6OTPSF0kLrsS5YuXoDj/zntMj3sR6HMbapD8MzsZbFZ2enUouL8QGv0AhecxBP/UT7UFMpmfqvw6w8MnnkiRsGjH1GV1TarprE2aUbjzYlLsYYvjDjBSr+1cPhDy6yZXZ99QYefXFkykwPt8Yw5hm945OncM2NL7jBMS5ShiokF5Sv57KvvUjbLafFWOERAjJwtyWLkVgWjGlGV0tLqHp/K1NGKIGEi2MTk3M6eNcVq1BKYyPjtILhTxZql3o7VUt/ovZjmtGpLGdRRcvRpsLFMYDjC9voWRwdRIOUTHgzsyCaeLBj7nSi8pIp6DimGT1QU4I2CkoUuxj7KFX7KPhsW0wRkLKNR9gyUDuitAxHCG7GXCKEUIUQq4QQj4e/lwshnhVCbA3/LTO1vUUI0SyE2CyEyDi9h6/tEM/9dPExYxl1cXThUzUYjI6olMvXseTuM5JK7jDccJL4ySQ6yYY4/AKhckwGbgael1LOAJ4Pfydce+06YC5wIfBbIURGTy/Y0krJP99i1S0LUk6t48KFFe19Bej9ljLNUjLxz+v4wRNXsbavwf7CYYBdPTa7/+2+2yHTAg71wHuBP5kOj2ztNSnJOdDHPv/R32KYLrJRB81F5ij2+RHe2GBRrbubaV9ZxqoPz+GbKy5HQ8lK7b1kYfduGO9MMutzyFyi/xL4KkSNElV7DTDXXjMnUHesvSaEuEkIsVwIsTxA/CydAFs/XMyCQvuUOy6iYX5BrJ90rz9WUF/Qhch32NCia+jrNtH07cP8te30SJhrqgyfSi7BVJ5voj4zqdRyCXBASrki2Utsjtk6v1MtyZS3P1RPaywiFUZJl0mzxZTHElPbYVbhXvpOnoo6fQrqbPtqplrzDgbuqGNPIGJ6yspONSsSPWvzhJHMxJFppZbLhBAtwH2EKrb8gyzUXksV9Xe+w71tJ6d17bHw8o5EwIX5OVlf7GNFstd5u9h9tge9rBAA7ewFMaW+lIICCrcc5t6WmA1iCWFXSTWdyq7W0s4wjKq7lPIWKWW9lLKRkJHtBSnl9YxA7TUr9J4e9r2eegXmdEroHk3Y/ciZSoZ0kOnuLOva0vr/0Zg0jHG/cumj7P9mgIGJxXTOyEEpK41qJ1SVQGU+vW9VclArTHmcRDXSna6J91sPuzHOAcNee80ORS2phcCOFgmUivplvmakJyi78axFKQ0ky6yZWpKHAwWKny/PepYpP9iM6ofg3v1R5+XgIIpfY/Iv3+HX685Ke9JL1oiWrMRO1CYrb4uU8iUp5SXh/w9KKc+VUs4I/z1kavdDKeU0KeVMKeWT2RjbwGBx6tFLo0maj4X1s53K6DQBZCJ90mmbCawTjorOeWXrKbh+D2pFeVRbfWAAZflG9N4+pnwvwP9uPp+ATC6tg53Gkoy0jmd1N/pIhNHzpmcCITg82zkDphUjKS3ijZVpyKP1u7W/4awRlyrtY3Ed/7GG19l9faxRTgYDESt8zSc7+fW2s4H0jHIGjGVkIg3P7jdO5rkeExlm1FnTueqU2Jprtm2TVBHTqslt6SsVtctoZxhsEknERBNIKshEs0lFqow2OD0ng2FzhY52xmHErz3RqZSlBCEYvGAhnU1egq9A75SchAUbk6UhFW2oS8vn/t0Lad1fzoSKbmCZ7bXHBKPvPr+KDxW+YHsuXYZINAFkumZKdM1ISb90GTQdo9JoQrI0XzNtFW/OOxG5egOeiXUE9+wNM7pCME8h98ID+P8zgZ+svoDvLXgsbXriRcJZ0R4s5l9tC2hfWkv12wHyXt3E9COtCI9zFtixz+hC0DM/9Zk0XYzFl9oJ6TB5PKOcE9J9Znb9pqpBpDq2Vf2u9x2i/aQSiqoXsneej/o/HkHr7gZdI//htxC75rH7HPCtLkBbkDi3XDJ0x+vjkY757PvBNApeXk9eXygVtNHaXGoq5r4SjjoGIDyx7pmxuCYcKYyE1T7T3yAek1v/Hw6o4ViuAsXPkQuO0Nnk413XrWTvh+fhaaiPtJPL19HwyH48/UPXDmd47OrnZ5Lz5DL0vtTqFox9Rj95Hu+buyqtH96af9v4PpJxzCOJ4WbwbE2w6Uhruwk+FTqsv7sWDuRUhM5VM9bgL4NX26bSPwH2XhKK+1Jyc/HUVNMzt5Kmazejkvm7E2/JeM++RZRtTC+T0phW3ZWCArZ9UeFj+fHj3BM9fLvz6ebmTve6kUAqam+qwUQjqT1ZM7Mkkv4GrL9zonh1FYmG4MSCnUy8rpOn2+dycH0xJc9uRiPkaqMLit45wMrWBi6ufIdc4hcOsbrWnGgF6NVzWNPbwDPbZ5HzRhH1922jeN+bcft3wpgt4KDOnM6mz1bw7fMeJlc4P9xkZ1hDVdNsQvKTZdh4Y402pndiDrsc5iMdfZeM1TnVXOuJ3gPVfttFeKzQO6FLhZ2Dlfx5yXuYfvs2tP0HEF4fQlUQDXWc+EAzx+fvQhF63N87Hu3deh7/3jufzS211DzrofylFrT2jrgFFM14Tj5oW8BhdL19qWDvAWpeF+z0V9qetqpiKtLxxzQft2tnVumdpH8yWkO2lgXZ6MPpRUvEPCNhB4mnvhqfzAOMZOSvikQVeszH2lYROpN9Hdx69QN0312IOmMqMjCIPjCA3LOfZYeSS4DiRPuqvsn8+ceXo1zRQ9ONKyi+902Ce/clzeTxMGYZXevupuj+t/j7s+8GnJnRyriRH9b0ARx/XDskYvxESOc663jZsC0ky7DpuvwyeU52dGVrcjEzeRS96JEPEPU+mOEVQT475UU23lweiZzT505lx6qh/RbJ3K+GwuaBWp7vmsPvdp7JK/97KqV/Wxqy6mdZ0x6zjA6AlEx7oJedg9FS3crEQPwZ2+n/OFpAIiRzXTLGwHhMYnduuAyK2ZDc6dKUKXMb45qld6Tv8LugoOMTWtR347wBQx03/t66+DGO3FOCPO0EeqYUMOVRf8y7aAcNhWZ/Nd9663LevOFEdl+ST94VByi+N731dzIY08Y4gM45hfRouYA9cznNyvHOqUKPSjoZb/0ec61Fe7CDtZ9ExkCjr2TGj9fHSNoJ4k1OZgw3TY6GNvNkipm5FTQEanidrUnF8bnnigCfbnyJn3/zPVR96iCyKN+2nYbC33efys5XJiMCkHMYal88SNOWdehxfN/ZxJhndKHD8fmtUTO1MRvrFoVFSUIqGNeYXwSD6c2Ma/3xrUxtd71dWw1hOyHYHU9Wu3B6MY2XPl3mshqZUpHQTpNlIsNVMnTEa2cd3zq5G0yuIqPeFy3M5Mngs9Ne5C93Labnvgoeaj2R/5n+HACDUiUgPfzg1UuZc1s7k7a/MdR/Uj0nASGSUvPHPKP7SwU+Yf/YkmFss/RWhR6V88Y6UURd52TYs3n57Y7ZTR7W/q3XxUwYDpNJshpIKq5AO9tAukhVO0lG40k0HhClkkfOhZncDE2GVurGu5GIXp/Q+MSk13jpplm8+tQJ3Lr6OgZrA1S87iWnRzLriXcI9vYmTW+y6LzhNNrfHYCAQu1LCkUPLHOcQcY8o9e+cIgfnn8x35j5H4CEhhQDBmNoUolZsxnnFPTID54snCYX66RhJ1mMNuY+dJS4dDidiyeNssU4yUJD2GpAGiJtDSERbfEYM57R1e665JZsOmeVbuLF/OOZ9pWlUeeyZS1RcnORc6fBms2oVZWUbe6lt64Q32EoXb4XzaYqrIExz+j6uk1UfGce9/xyEbOL9jElp50qTzeQ3I/oJIHNzO44tg1jOknZRNpFSLLosRMPuiPTGn2alypRjG+5Lt4yIRHSNUqarzWPY448S+Z6o32qdNjer+kZm5+dmT7zxJuMBtKt5/HQngXUvpE9I6g87QSa/5+HgnU51P92NXpfH8rGHejBIHp3D2JVJ/VvhvZ5JFoKjHlGB0CTLNswlZWHm/jsxU9S4zkMxEpG6w+WjqEuNJwSYcxIe5uXIe463aF/6zLEMAw5vd9GP5pUbNVS89h2dgRDNU0WyTwXx3NxpK2T/SPZJVKi9XSMpR0ZWZ8DIEHDucyAigQH24CGQmewgCMBH/5KhXxFta257gghUIuKQtlrKis4dEYDXU0KJ1+4jo+Xv8OPSi+E34aaGjHueopLgWOC0aVXYXJjOztbK7n91fP5yTn3U6pEB/1HJF2YIVNRx60wro305eCXjb5GS/gC20v98IuVgN64WoiDMFLCsdlOtogo20UCOs3LBztjqNWTETlu5ylJcCyeHcRuUrO6y1QkXhGMbheeTAMWZo9MhmENye5XUAkF0nxuyou88f+m0/xyI9qWbTYtwxAC7az57D4jl7wDku5psOC0LRwcKGVq0UHeX/JW1IR/ZHcxen+/c39JICNGF0KUEireMI+QzPk4sBm4H2gEWoBrpZSd4fa3ADcS0jQ+L6V8OpPxI1i2joKP1zEnuAvpH+Tr3/wAX7zoCWb49gFmxjRetsTWVLOqGE8C2QVdOFv7Y49HvUg2UNGTT1VkY1gK3Ussk0WNZ7nE6nlIxqhp18bK8OloEHZwel5RBkmb52qW5BAd8YZU0DCMoBqDUrW1fySy+C8q2s7T35/N9FsaCW5vCY0zYyo731dDsFBSuhk658BHLn6R6Tn7HfsxQ6qZB89kFOsuhLgbeFVK+SchhA/IB74OHJJS3iaEuBkok1J+LVyS6V5C1VnqgOeApkQJIp1i3eNCUdnz5UV86xP/pEI9AkSvE80vhN0623zccQgbhrY7bl0+mI87GYWsE4vB6Ib2kMi1ZzdmvGui1s5hi7O1HyuSVeGNZ211XcXrMxU1PBkY9+ETGrkiEHkHDKlpfDdLc/NzMO7H8Kkn49r7x95TaXlqCv5yyQVnr+S04mbMewmSdXUOSC/fe/sSmm73I1esTziuU6x72owuhCgG1gBTpakTIcRm4Cwp5d5wXveXpJQzw9IcKeWPw+2eBr4jpVxq030EaTE6IHJy2PtfJ3HNx1/gXQVbos4NOhTMi+sGc1A90/HVm8eKkipE//AGA1rHdZpgzH0bE4LxwsbTSswTSDK+YyvtZlrtYDexmmlONKk6wWkJYdfGeCbGBijjno3nHgj7vK3P3DpBORVbTMVjoEsl6ViAtX0NrLlqCsGW1qT85U6MnonqPhVoB/4ihDgBWEGo4GJUSSYhhLkkkznGL25JJuAmgFzso40SQfr91PzyDZ7f+C46f5TP5WUrgdBDthrJrJI1yjrt4Hoz4OQWM1/jBPOLFunP9AI42hNsfu94LiMj4gvsLd/J0muFk5bgZJQ0DIbW55/oGaYyETj9Bsa63CuCkeduMG3I0BZubCyTHCZe8+9jF5Bjbmscz2RHZHNvFfJgZ8ax75nEH3qABcCdUsr5QC/hyqkOGLaSTPGQt2IHz//1VLq0fFR0vCJIrgjgFRpeoZGrBMIvQei7Ylnn2q0n7aRG0pI8JtpNx4cWE8sd1y5gYUarymuWSlYGitrIY93YY7Phx24DUKL7crrWJzR8QsNL6GMEBalCH2JCocfcn3WziQG7CcB8fWRfQ5g24zf2Cg0fOrkmpveF6THGiWhaxC4nIn2EP1bNxmky1aWSchTg8u2TU84mY4dMGL0NaJNSvhX+/iAhxh/xkkzxoHUcpO4fG3mi84QohvYx9MkVgajvXhGMO1sbEsj8AaL+N8MulNUszUMvR/R1XptoPyf12Byuafw1lidm5jBTHJrsgpGPlYkSTV6JJgAl5gmFPl4RHGJ4EYwwfRTzmyaAZDQMu7Blu+WFHjYEmie2mIldxE4mxtrcmKyMCSnZyT3yTERqm41e624Krc2P5jZVKeU+YJcQYmb40LmEqrA8xgiXZEoEva+Plp6K0P+WSCxF6FHSfGink/OPYsfMMZLTIk2AqA0SiskPH8vkQz+so3RNwAjm88bHF9ZgcpUAuSIQI22tAThOsJv4rB/zmMauMF9EogYjz928Y8xgHkOiWu/VSpuTQdRpAgrFnqsMSA8D0sMgIQkb/ZuExrPaK4zj5mdkO8Ek4VFIhtkHpcoLj52UlAEuGWTqR/8c8M+wxX078DFCk8cDQogbgVbgGgiVZBJCGCWZgmS5JFM8yECQ7SsaaKmvpNHbATivryJqL8mrWXZuJDspbvWFD7UZaj+0trNx00W5xIbU9BhV3mGtbKioxjnrC+cj5IqDIc3AeKHjuw0dfNsmrUUJ39PQufD/Ykja6FIZekZRj88SSyCdtSc7mGkfkN7IOt0wihm/jSoFAdQoQ5zZGOsVwZhoR6eJxhozkepmngf3L2TK33aRuSwPYcymkkoValUVR/5RyE+bHog554tIVkFAqhFVWjcxUzxrtHVjTLxZ3awxeMPrO7MlN8oCHqbBfM1Qu/DL6BQea3mRzEYoA9HjGgbAkAsoINXIGt86kTlF4Flf6ND9OcUH2E9Gxj1bXX7G+HbuUaf9DdZ1NoRcbPmKHx96zFo6JO29Eeu7U+yBkxfB7h2J542wi/PXpcKeQBlPfurdKK+tdrzWCcdeKqkUIVSFjldqeatvOrpU8KFHPooYYkyD+cxrdav6aMDJ8BNvL/ZQlpQh9dK8Vjf+GowXYtAhY5FdnwaddllojPNmI5TxMYxRxv0aTKlbpJZZbY5Zmpgt6A6Rf1Y12rhn6/o4+jk5r5ujAlgcmNxMj1X6G5O50WZoglciBjnzPTrByuQp+/cd1ux/bV6EZ/mmlPpKhGMiBDYZBPftp+FHB3jy/gXcfdbFTP/YZm6d+DgAAUNS2Fq6h0JngZCKZ6OaGYi3ndNgcAi9YN7wX7M7zTgXfV20ChyQKtY4e2OPVHSG1Ghjk7V/69+hewyC9KCaAkqi6LEYuhIF8BhaSjxNx2ynsD6PEE16xJhmXa7E07bM9xaRyCgMSE/MpBz5rcNjOUljq5vVzOTJ0mW9V+PvxoE6qn+RG8owm0WMG0YHQEq05h1UNO+g5/Fqrv/jR/nLcX/Di45ueunNaq+GEs20Ni/hkKROHEBiVReNFypWAmq2ue9CL3nsyk0J9+tFi7nGPJ45o6nZHhH6a+pPREv3ZKSV3TMxjx1RpU19m20SBh1O++QV66RrIEVJaqbL6N8rtMhSxvDtK4SltqX/SBCSA7OnAusS64F/n8mkV+LGkKWF8cXoJgT37af6Q/18+GP/w88//3vq1B4CKKgy/PILCVLFKint1stmyWkwo2Ji6OgXekgKmo8f0Ip4u3cayw9NYvfhEo60FyD8oTFlrs45x23kCzXP4TWtLY2xdCmi+g3RY3YriQiTWyWVWRsw36eBoTjwWOOSnVbj9Hys/mSDVrsxrTCvl51CimOusWE8zTD2GSHF4fPmdOEaIYNcgQileBqUoVHN5+ONaYdE+eMN7AqU0/hIJ/ow2M3GLaMDaN3d1N7xNret+DD7vjTIA/P/BAJ0E6MqYeZQ0WKYFcwvthxiMGlRa0VYOsS84KG/H37tE0y9C3zrd0FXO3XBPTGRUPtqqrniJ59myZl34LUauMLvnsHwXtOa32ByY11q1RyMv04vsNXYZ7cMiGgakXN67DOKZ8y03I+dodGqLidjcbfbfWemW5NqyBhKSJIbewq8YfuMdTOROaTYjqaY8eNMCnZ4oOUkKjfvSOmaZDGuGR1ABoMor62mYWM5F//4Czx+/q/IDa9NB1AJSAVNhBkhLNnsXHNewxCHRBemH9i0W2tQhF7eTf46frruPAqeKkQdhNnP7SC4dx/xfI3BffuZ/Q0Pdzx8Nl+a8FzUpKNIEVp6CCI0gNXQFF8iG5stYqzkkdDQaH+zXUy+9ZmYJ0GnvdzWa0LQo5dMcihho5WxkmGmiCRHCUvooYkiYOnD0AQGpDdKizDcbnG1CMszjmh3Saj02/0TKPhtSdbX5gbGPaMb0A4eYtbn+7jqW//Dkut/SpGQhF5/GVHpBwURxtUQdGn5DEgvewJlbBuYQK+Ww/EFu5iTu5uJ4V1zEHoBenQfj3SdxENrFzDjt4NMXrkxEvGUrK802Lab5b85jcHvvoDP/PKIkDS3IiStrJF39kE65muIuWZIKpuNjWajmZ3Kbr5/Oya3agLR56JLLVm3oBq+9JixbDYARcJYERjJJczSXrG475xpjR0nGVtMIgxKlb8+cQ5Tnhq++LFx40dPFsLjoeOGkzn1Uyv5/IQX8CIJIBiQIem+K1jO7S3ncvihOso3+fEe6kM52I3e3QOahlJchFZfxaF5RXQs1Pni2U/xxP55HPltPcWv7kBrP5ha9hEL1OoJnPncDq4tXgVAAIEWluhmNRqG1uYDkW2u5kiv2JpliWrCm336Buz8/FZrdrK+/tBYset2Jwt4vJ15xnmITrFlB2PrqoqMhA7bWc/tfOapqufG/Zifz483XkTtDXvRug6n3JcVWd+mOlIYaUY3oFZWsPv6mUy+cjtfbniKlf2N/OapC5n2rz6U1VuSVrE89ROR3T2h6htZwv7Pn84fvng7VYqfHumhT/fSEqjkvn2nsLOrjJK8Aa6pX8HivOYYI1y8wn7WicB8zI7JzW2d/OHGtcnCyUBntW+APQPG9GNZ11u3HRthrebJINEOukyY3IrmgRqe+taZ5D2aHWnuMnqaUPLzkXOmoh46QnDHzqyXykkHwuOh/6IF9NSr5B/QKdjdj+dAN9ruvcjBQYTHy/6bFvKTL/6RCjWUW8wc8Rfrgw/BWQ2NL/mthjkDicI87eMW7JcOYHEJ2vTrxPCGpLa6xABbo5uxacVpInEq6uAUV+CELq2AO399BRPuXJq198pl9HEEz8Q6Dp41CUWDmk9v49sNj0c2cIBzBJsdrBLWrJ5niriuKgeJblxndZ9Zw2SNCDdr1hhzcomQhyIYdcw4bt7AYs40E49uu2fi1LZLK+BnD13OlO+uQGaxWsu4D4EdNxCCzV+cxP6zNIoeWMbeO6fREqyIrJvjMblZNbdzJSaC3fXG8VS2Z9ppENE7C01RblKJkbwGzXaJHs17Eaxuziga4ljY4z2TZLbWDuhefvZg9pk8HlxGP9YgJdMe7GPOD/aBrlHyr+X8dNv5Mc3M8ebWtbShNg+imhg+NtLNDLO12spAdgkXEqm31sqpxl+7GHlzLHu8ffR2fvUB6cXYyRZ9/yKhO80ceZjshKhLhVtfuZKpt60dMSYH1712bOLNtRGXncjLY275XsBY4w6t0Z12T4G9Vd4Ku2CZ0LUyrrQcuj7W+GeHoTHi0xQaSzExe+xaWiE6q24otVhow495slKEHjcgxi7wKOq8GMp4a9DQp+fwvWevZPa3NqENQ4mmeHAZ/RiHUlrC/MJVke9D1vToMFojSszOCOaEeNIsmtkTR8WZJyBztlQzrfZjxEKPmkCiN8aEENq0Y6d5wFB2H6vRLZHKbnfMYPbNA7X8/cFzmfnT1WhZSA2VKlxGP8ahTSih0dcRYgpLKGoMs5vU9PjMae9Kc7KwO/VnTURhdeVZPQTWHAFR95mEpduLFjbWGUE3SiRizyu0yP/WuHzzNuKYcF8Lg1v99Qe1Qm5bdiHT/qAz6fWlwxLHngxcRj9G4L/oZHKfWxOz7tt+TTF14RJVqpQRtjFLM6fINDt3kZNUs/rZY2PkhyS3mZGdYNfG6l6zBss40WtNP6UKEcpgE450xGTkI849WMex2xdv/K+jsOTgiWz6+TyaHl45outxO7iMfowg56nlWF2l6szp3PjeUK1u8w622CCYaOlptyEl0cRghjVdtVMbQwuwS3lt7seKZIxkMceithqDKiShtFmeaNu8SWOI5ymIl/xCkwq3t5xL7hfzKFz3ZhYckZkjI6u7EOKLQoj1Qoh1Qoh7hRC5QohyIcSzQoit4b9lpva3CCGahRCbhRAXZE6+iwhsVELR2899209Cl4JBlEiaJKsVXHNa0+KQBMIGybrPrFVKDAOY3cexDxFrWTdLbmvW2ahUX8QmeYQhN501qUayKZqN6wLSw5eXX0P+JyX6uuxmickEaTO6EGIi8HlgoZRyHqAC1xHK7f68lHIG8Hz4O+GSTNcBc4ELgd8KIZzLV7pICOHxoJw4BxT7xzg4bQIfn740JljG/DFb2a3n7aBhn40nngvNSUU3awnmnXFxJalZ22AoBwAMGeHsItnMab5UpH0CCxLvKU8ULPPFl69j+n+1hCqrjCJk6kf3AHlCCA+humt7gMuBu8Pn7wauCP9/OXCflNIvpdwBNBOqw+YiTag11Wz6XD5KXq7t+b4JPso9R6K2qUbWsaYX13C7xZPs1i2tZj+7+TMUaiui+jT+T9aiD9FMH+PzN/nOrTRAdOCKaurDKOJh5LI3S3fjmlxl0NaKHqUFySE3nip0BqXKV9dczZxv7srK5pRsI+01upRytxDip4RSOvcDz0gpnxFCjIqSTOMBwT37mP3lPkef7EC5QoP3YBRzKUJHM6VMisnCgoIPLUa9Ns5ZkSzjxtummqid1chnF2tultqhLK8ag6hR/YXSSQ+NPWRlH4pvHzLo6aiOoTdETSKaVPjm+stp/HwnwX3JVUgdaWSiupcRktJTCFVHLRBCXB/vEptjtgu/bJZkOqaha2idnTGHlYIC+q84hfJr29ClYqoxZki0ocIJVpiDaZzW8cmsW81tnJYN6ez+GpRqVO51iGbyXBGIlNyyK8KhokfamK3nBYo/xs2YCMbze6LzBOq/rhHcPeyFh9JGJqr7e4AdUsp2KWUA+DdwOqOsJNN4hBwMcOCD/Xyj8QmASIink+pstyZ2KhcVIyFtYNfGur00JlOMVGImAOuEopus7TH15kwRbQbs1uBeoaHLkGFSQ0Qb4yyuO7trjZJMxnUbBiay8UfHoW3YEtN+NCETRm8FThVC5AshBKGSTBsZhSWZjjUoRUWOBjgAGRik9q5cdgfLooxc8WDUnTNg3b9uHHNCPGu5mYGsktj4RGqc2UwMMXnmiS3XZDCgkf/ebDOIXCv0SJGGAemlV8+JaAeRwg02AUXGtUa1FmO8Xj2He357AXmPLov7bEcD0mb0cHHFB4GVwDvhvv4A3AacJ4TYCpwX/o6Ucj1glGR6ihEsyXTMIRAAGZ9xc1/fxK3LLwOGDG7mj+EPN1RcgzGi96cPqd6Gvzta0oso5raqvcmqwVZreVTZI8uY5iAZM9Mb/xvtIvXOYzwBpoy8FsNevEQWBrxokayxNy+9mpq7Vo6KHAWJkFHAjJTyVuBWy2E/Ielu1/6HwA8zGdMFSWW30Xt6qHnYR8uCKhp97TGqc3SFFSMxpIwywtmFrlotz1hUaKctmskUM4DoHG92iN6sMjRh6YQSP2oIvGiRklIwVFlWNWsWpv/NMG94MUfKGXvXDePbl9a8j9nfPEBwmJI5ZhvuNtVjGAUPvc3//vtKBqQXIFZNDks+Y/1uhV1BhdDx2PW2WRo7RsNFhbIMfYw+7a6LWaeH25sLVUTW9OH7GtBDqnkgZI6LOmdW0wPSEzWxRIXOWo14YQbf5K/jU0/cyOSb9hLc1WZ7n6MRLqMfy5CSaT9Zz+f+cwPtweIoxtJRItLLUMENd5T5fzPsUjlZ+zTOGfu5bQ1vcazvEfXcNFa8icaKAENWeU0qQ1b6cL9WNT3Sp1Tw2dSkhyFp/nL3LB79xDk0fXkl2sFDjjSMRriMfoxD6+6m6Surue03H2BfsCSWyUx+aYPBzQxmTv5gbh81htlSbrPOtrZ1pDUJ9d6Qvgbzmuk3vhuS2vgY0ntQqgyGVXrDsGbQnKsMUqV2U6T0x9BqSPfNPdWoa5qP+gaVdOAy+jiA9Pup+d1yvv/kVfTqobgEg0EM9Tba2qxYGNxZHTfDzOxWKW/daWaW9naTh7lPY3kRc85iFzC72qz9GhNARIuRCgO6L6q/Hj2PAGpMbLxXBOnR8ui8YxL6CCeMyBZcRh8nkIFBZt66gW8svzLmBTfHnFuZyrwH3IpEZZHsYs+tTB2VJMJqXTczrCGtTXvNksnPFhnTZplgMLTB9L16Towxzuj/963vpviJdxKONVrhMvo4gtbdzYzP7OQrb74vJsWxVapHjFVht5ptsgmbba1OLjaDUePSZ2FsY70dUb9R7duaDIJG2WMneqJodaDH2PVmfAAOP1SHfhQyw2QLLqOPM2idncz6n538ZKvzLuF4/m8nKRov1XE81dxI0mBr0LNZfxvnjLZWA5tT0I313pwkvNGvERQD8LW1V1PzULMt/WMFLqOPQ2gdBxH3V7IvWBptjLOxRpsRlRBRRO/8sp6P56pyUvmdjtvFtdvSFxPXHifHm+VcZIebIc0JRdodDBZS86scNn17Kp6aasf+RjtcRh+nKLtnGT9+/eJoaejA7FHli03hpwbipVe2k+bW9lbPuhOs2kSi9ga9dsesywujP1XoeBna8HPPnkX49hzGUzmAHGXBMUqu/fZk27bDSIeL0QyhIPzR8fLWYBen/+2+K0QXh0g2Es4JTmo5RBvJjLENxG67tVffY46Zgn6M7wPSy+4nJ9M7s4Lpn24ddfvMldISPBPrkmrr5owbx8jZH/IrmwNFrIkpsBi2zLCLDc+Uwa39mY85SW9Hld8pQs+qYZjUdiOcFqBHyyOnU6J9pgNeS5r8EUNw/4Gk4+xdiT5OIQODTP3jdh7rmB85ZreWdtxSamLyVPeWm41vducSwWyYSxSgkyyMOPZIUA2CfMXPwQUaezdOgMqyxJ2MNFLYTOMy+jhGcO8+tv5uNl1adBafiLpsChs1doJZfeBOm0PMsMt9nqg0sRPMQTJR9dds3HeJqqlYzxkRdQZ8QuPO8+9Gz9WRe0ZP5hiRk3oyFpfRxznKH17HfXtOjjLI2W9ZFRGftlUiWy3hkWKGIn5ZI7NLzY7Jo7armiYGOy1jqOJM7HGz8c1JbTdgSHZjTK8Icu2pbyMm1tjew9GA8HhApJadx2X0cQ69p4e9SybH+KidXGV2PmhzhJm5j5ixbBg6bn52U7/xJoIIHTbuv8h3h+oyhmHPSFxh9Gse84G3QjlM1YpyR1pHEnpvb8p74I9ZRhceD+KkuUebjDGBif85wIreRntmsjHAJRWjbhPcYiDeGt2AnXrv6D93SNtsR5uRyMJIC1Wg+ClS+kOS3CbzrSYVmmbsYctNVQifL6a/sYJjltFlMIhcNXoS6I9maFu2cf+Lp8ccH5SxSRhtr0/ieCJ/t9HeGvVmDsyxLhGsy41IPzbRb8buPAitvXOVQYqUfvIVfwyt5iIOqtDZsrWOptu2Ety7L+E9jFYc2+413c1UlSwKd4aYJqqgoIheFzvBLJ0Nhra6w5KJenM6pgo9pr+YTDAJYM5Iq4gABWIwkkMu3s45L5BX2Yf0j72tqWYklOhCiLuEEAeEEOtMx1IuuySEOEkI8U743K/CCSVdjAJoZ87nmhtfsFWN46nndrDGrSejpiekz+zKs1mv263HVSQ+k3peqvaRqwSGSiKbJLwqdAKoMRoMhCaUL819juCJ0zO6h4wQJxFo0l0k0eavhEoomZFO2aU7CRVlmBH+WPt0cZTg291Fn+aLSMxUGdsOdmt566SRcDdbnJRVVm0jql/TJphcEaBC6SVf8Ud2o5kLNcQry2QsCxq8BzmwMC/Ro0gItbiYI9csSil0Va2qQinIvIhJwl9USvkKYM2bk1LZpXB+92Ip5VIZKvn5N9M1Lo4y9JY2lrTMs3Wbxb0uBXXcet7KmHb/251zMgLa7UaDkJ2hV/rwoUXlmRs0lakyYgScNBAdhf5TejOWrFpPD6Ur9iOKigBQKysSusm0g4eQ/f1x2ySDdHWqqLJLgLns0i5TO6Ps0sTw/9bjthBC3CSEWC6EWB7AnyaJLpKFDAyiv10ayZZqBMg4JXOMSuroIKXjbQcFbF1f5gQYqUTamWGt/qqj0KvncEArioxrhLkGpCdSyMEJRj+fPv5l9n1hEZ7JDY5tE0JKgttb0NrbQ7T1HMFTV5vghjSUKZNS9ptbkW2ru1PZpaTLMYFbkulooPEfrdzV9i7H88mq81EWa8suN2sOtmQzxETR4WBhN8OcZ97Y724exxzmGsqq47HVTszXzMjZx3c//Te2/Lgircg0O0i/P6kyTnpLW8a549Nl9FTLLrWF/7cedzFKENzVxuE/NdCr5zhLabuINHMYqoPhzo7pU4FdQIwTDD+5YWwbkN4IXWZpbhj2zDnkzB8jw44xWUAoau5j85YiZk1N6z7SRTaSUabL6I+RQtmlsHrfI4Q4NWxt/4jpGhejBKXrumgPFiXMDJNMwAw47x6zwqwFxOsr2d1oBgxXoVG1xZo+yy5+30y34bvXCdVS16TCU3vnoPTYr5mFx4NaPcH23NFGMu61e4GlwEwhRJsQ4kbSK7v0X8CfCBnotgFPZvleXGSInhklVHl60pa8BsxMac00k0hVt0p+J1qsjB+vwqsSrnluhLiaCzbYppQyre/VsHYwKENx/rt2VKG12hduUKY1svmWqUmtp4XXh1o2cjviEgbMSCk/4HAqpbJLUsrlwLyUqHMxouhpUCPMkE0YG2Tsos+ckGwJY7siEwYUk6ZgWNwNhtWkYptpxjqxmLUATSpUTzoEqgpBm+e0v53SjVUIVUV4PCjVVQR37optR0gd1zqjVXK1soIji6eR99iyrNdzO2ZDYF2kBrWsjOAZh9OS5jGbXOIkiow5blo/x9CUYF1vMLWxLo9XMdacIitqq6vFn26+D2u8vSp0zqxtRq2qtL+XrsNU37+B7vctRNTX0nZlQ9LWcs/kBrZ+tYl531yLcvyspK5JBS6ju4BTjqPjH5V8be7TgE2kWRLS10CyrrGo1FB2EXkObjkzDOa2U9nNTK8hIrXYzNtnjbGtkt2AXUxBvjIIPq/zfXUdpuSR1WjNO6i5fSnKcTM5cs0ix/YGBhsqOHnxJi4sXcvWr+aiFBQkvCYVuIw+zqGWlrDtCx6+POMZCpSQKplqBJttv0mute3GSxWJar+bx9Bk/GqtEJ1Oy+p2a+0vRx7ujnt9pNqtlCgHu8k9FH85ZBjwjgRzUIXOj07+N0fOz+4q99je1OIiIZq/OofvLnwgrWvt1rRO5xzzxGdo+DNgJ9WN4hOR8aUnutySpfSSFaFrovvt17wQSN6OEdy9B08cX7lnaiObP12DVhHg1uqVqEgq1CO0XSBpekRkba3uMvo4hvD6yJl9OKqIIiTHfHZWdbs2BrM5ZXaxmyxSjYpzTPZokxFWi9RKl/auw3j3IxUqc45wOL8cuuNL9WSx6bM13PieF1mQ3wIQ2Rd/2nFb6czPz1qtN1d1H8cQs6dy/fRltufs1sfJrJvjubicVOx4obKG28z6ide3Ux76RBpGvGPGmn5m/j7659XHtEsHSkEBE2a3My9vV2QsY/KZX9KKUp4995vL6OMY/Q1FVHtjc5XbbTJJRsqa93wbcGJupwkh1TZmmBk+itmt2W1NIbLGX7M2Yzc5GAxY5emm7dzsZJrxL57NVQ2rgegkH4rQqfL0IPOT3+WWCC6jj2O0H++N7M+Oh1SY3Pw9njXczJDxJHW88eKNYSwXDCmZiNkhcWRepE0Wls1y8YnU3LqNKk/PUEHLcICOD42pvgN0nFaV+UBhuIw+TqGWlTHjwm2R9Xmq0jOtMcMckiwDJ0NTIo3BzOwGDUbOOLMP3UgOmaiI5L5gCY1L0q+qKnJyOHjjaVT+704uq1xNvuKPWbKUKH5q1F4mfLQlpb3r8eAa48Ypes5q4rqq/6ChxDB7Mu6qZBFj+Sa9ScXJX570EsDiLzeMgJHJJ4mkGxqCn719PjNXrU8iA140PBPrGJwygW0fyOHmcx6lWOmPol8JewAC0sNP953PmaWb2Lq/isbBA/G6HYKIb6F3GT1NiPlzUXbtQ+s4eLRJSRnC6+PAdf2UqvaSKZ7Knankt67fU+kv1fbWcXWiUzx7RTCKsROlvPp583nMvrmNYJLFFtXKCmR9NVuvL2b6gl1cU/sMxUq/7X0Yk5GGoN1fyC9+cy2TVvclnfdQyckZ8t/bwGX0NCFXrWespp6UmkZgILWfPluqvfGSGwxHksybygSRzoQQr2S0KnQ29E8k585ygvuSq5OuFBWx/y8VXFC/jhvydzl6BczwiVAWnBtq3+CPb1Uhl72TPP0JJh+X0cch1MICpkzsGLHxsr0kSMTE5vGcPAChCLnopJN2CEiV7254LxNuzyP3JXtXpB2UynIW1bRyYsHO5NqHU1oXKH4OaYV0TyugKPnhEsJl9HEI7UgvB3pCvmA15dVm+jAkuXXNni6sE4fZHmAY88wahAFjfe6UVUaTCpsHarl320kMri9h2u3NkfRPTvBMrItkixEeD5s+V8uHi1+JSZPldB/GciIgPXzr0euYvC+7KdRcRh+HUHxeKgtDEVdmY9xII12jXDLnnP53gsHgbx+Zyr+XLGbKA51M3NyM1DS0JNbJ5pRQcuEcPnP+M7ZGyHi0qEhaA+VM/1sn+rrsFh9xGX0cQh8M0NZeBpOONiWpIRHD2u1CM0vUeH5yBZ2vrria6d/vZ/LGN9EziDFXjgwSCIfamulJRL+CToXnCFpxToa6jl3fLsYfdI38t/IZ0J23Ww7b0Cn4yFOB2SduLuDglKbKK4JRu9h+3XouM77ehbZhS+aJGA8foSNQmJI7MFcZpEDxs66/Ae/27Jd+chl9nGLiPzfzp5bFQEh9H+uwWwMbjGy3B93c5sH2hXg/l0dwR3KGs4S0HGjn0S3HpXydiuT3q84guD9J33kKSCZnnF1Jpv8TQmwSQqwVQjwshCg1nXNLMo0BaB0HKfu8ZEnHCaHvxwCzQ/QGGbM13cz0RgRce7CYL694H7t+3BSS5FmC9PvJfbvQ9pxVfY8Y4pAc1AtoeMCT9TRSkH5JpmeBeVLK44EtwC3glmQaa9C2bmfbn2aa1pOjm9kTqfwxu+0swTBG0gkVSXuwmDt+dSVTr99A7pK3s05r7RtH2DxQm9RuO4P2BztOpuD1rVmnBdIsySSlfEZKaey+f5OhnO1uSaYxhqpHNnHXzsVHm4ykkYov3pDe5u8Am/y1/OJ376Pq928j7ZI8ZgFq825a+isi361bZyHWr//aS/PQOjuHhZ5sTOEfZyh1c1ZKMrkYOWidnQT+Wk2PFto8oaFEPqMNSbnJTEY3u5RRXVo+D956ATW3Lx3WstoiLw9PEjvmjFrtANXLh8/NmdGvKYT4BhAE/mkcsmmWckkmt/bayKLk36v49cazYo6PNoZPxlIflevdVJHFwP5gCUVbu4dlHRyBorL9xkmcUbLZsYl5Q02uCHDPgVMpfi67vvMoktK9UAhxA3AJ8KGwOg5ZKsnk1l4bWUi/n/zHi9Hk6LePpuKeMydzMBJDTvZ1sO1mX9azrJohTzuOa698GV+cvf4aIrK0KBCDrHxlJlpXbBKQbCEtRhdCXAh8DbhMSmneAuWWZBqjqFqyhV83nxNz/GhFzQ0XfELjc8e/BNMyqIoaB2ppCTsuzaMp194XHhV7HzYM9kofNW8N7xapdEsy/QYoAp4VQqwWQvwO3JJMYxlax0GKf17ETn90cYLRpLpbkWrgjWGMK/ccYc/Z5ahV2cvgYmDbl+fwlcsfjTtBRko7C51cEeD57rkUvdWadVrMEHI41ypZQLEol4uEbfUnF1mG8HjY8rOFfOeCB1FF9HsxmiW71Uhn3UhijYpThc4zXfNo+WAd2tbtWaPDU1tD0YODXFq5Ju7zMnznBYqffMXPDc9+kqZPLc+K3eA5+eAKKeXCmDEz7tnFMQMZDDLre1u4fWusCj/aDHNmWKV6MtViTixsZWByallWExVG1CZWsrh0W1KTolcE8Yogh7RCpt6vD69xEJfRXVigHTyE757ySBCNGaNRqiequQaxFV0BChQ/SjA15pKahhx0rlWuth9m20D85YB5S6omFf60+wx8b25MiY504DK6ixj01iUuWzQaEMlSE0YyUXOGFT5XBOifkGLaZl2LW1AhuGsPj7+4MGnNxyc0Nu2uQe+3r7eeTbiM7iIGaj8xa/RsIJmqp6nCmnDSyuzmXW3mKqntwSKKt2Sn2koEukbTH9pZ25fYoq+Fs9MWvZ437Go7uIzuwga5Xfqw+NTNVvJUmD1hSeQEktxaGhngsJaP6HdWw9OFGAwQ1J2r1Rhquyp0dgUqqF46fL7zqLFHZBQXYwsOAsYcHpuuYc4u80uyTG/WCFLJIGO3Rq/09BCsGIagmaCGR4ndoWaXPuuO5rMQG7dlnwYbuIzuIgal67vY4Z+QsF02rPBmZrcyfjxJnkqEXESiG/nkUPjRk1egLMu+EUzvOMhjzaG96Hb0G7QM6D7Ue8sTZm/NFlxGdxEDuaGZv77xrpjjKnrMJ1U4MaaZaRPlVrOD0/rfrK7vGqzgJ80XcOsdH2HmdzYgA9lX3fWBAep+56M9WBy33fOH51DxSlvcNtmEmzPORQxkMEjpeg/EutOHDebsrckwuV0bu/X/gPRy/76T2fb0VBqe6aZ0XTPFA9uGNSe/97V1/O6f7+Xcq5axoHAn+UpoY9YO/wTu2baQ7v0h37mnbeUwUhENl9Fd2EIJHJ2IyURZUs2wBsVYr13RO4Vnf38aNfdvor7zDSRZqY+YENLvp/62t9j2txpevexkii7fS9vecqb/IUjtso3UDIMmkQguo7uwRcm2QXq0XIpU+zWkeX0+3IE0Ttlb7bK+mvHAyoXMWbKT4DAlc4gLXSPYtpsJv92NclcuM/wtIOWITDR2cNfoLmyRt/UAG/rqst6vncSOl4YZ7OuYm691uv6MOVvwN9WkSWn2oA8MjIivPB5cRndhi+CuPTzVPDspw1s2rO92vu7oMYTt/+Z0zuZrVSSLSnbg7RwZq/Zoh8voLuyhaxS8UsigTcy7Famo7vHcYU7M7qi6m9I32+Vvb/VXoBzMcvTbGIXL6C4cUfNGF+v76x3Pp+tis8Jpt5lZSpsLMyRS9Y3zD286geDuvRnTdyzAZXQXjtDXbOTJf5zOgIyt6HK0d7IZEtycKsoMDUHJs/nDmgByLMFldBfOkJL6uzbygxXvTUqFTxdOUjqZaqu26rvQ6dIKKNs8/LvCxgpcRncRF1pnJ02f38mtL15Frz48iToTqe2214STPVolucHsf9q+GHXF8GVVHWtIqyST6dyXhRBSCFFpOuaWZDrGoB08xKwvruUHr12asYU9Hfda0n2j06vn8O31l1L+bd+IxZGPBaRbkgkhRANwHtBqOuaWZDpGoQ8MMPv/DvF694xh6d/qMrP72ME4rqDTp+fwg3vez8QP7kSuWD8sdI5VpFWSKYxfAF8lOqrQLcl0DEPbso3nn5kPJJfCyQ52O86stctjDGvx9psbG2FQ+Obj76fx/1aj9/U5th+vSDev+2XAbinlGssptyTTMY5JTw1E7cxKN1tMsmmarVb1mMg4EarH8mTnccz8zV6XyR2QMqMLIfKBbwDftjttc8wtyXQMwbNiM3/ZdmrW+7X60g1Dm5PKDkNqe5eWz9pfnpC1+ubHItKR6NOAKcAaIUQLofJKK4UQNbglmY556H195P2jjC5teEoaJUrVbGe4+87bl1Hy0KphoedYQcqMLqV8R0o5QUrZKKVsJMTEC6SU+3BLMo0LFD+2mtuWhWyp6eSAS6a9NWurE+5oPZvZtx5E+l3NLx7SLclkC7ck0/iAPjDAzJ/28fjBE2xTQGULiQJmftJ8AeqXi12VPQm4JZlcpI2+qxbx4R8uoVSNznWejKEtlUnBqq4HpMotz1/DnB+2EdztuAIcl3BLMrnIOope38HqI5OAaGZMhonj5Y6zwizZA1LlllevZvbNm1wmTwEuo7tID0Ig8nIJ6KqtgSxZZrf61eNd16Pn8q1Hr2PWFzaidbvbT1OBm0rKRVpQjp+F91eHOLd0g+35VMoZJ4NXe5p441cnM+2+leiu4S1luBLdRXrY3saabQ0RCaxa6qBZkajogtPE0KUVcPPSq9l8YxNldy91retpwpXoLtKC3tPDnB908Nc/L+am+ldizjup4IkqrOwPlLI/UEy+Msg/ty6k9L5Cmh5b7UrxDOEyuou0oe/Zx4Ejk9g1WMFkX3so/jFcPNCKQ8FCXuqcyaCmMq2wgxd2zyDPG+Tc2s1MyWnnTy3vomcgB+8TpVQ/04bM8dGwYysyMHjUMqceS3AZ3UXakMEg+pOV/KHgvRx3xUauqFoZUd93+Cfwckdop1trZxkFDxdT8eh6pOZnQ3E1Ew63ITWNZQ1NvNp4KiWb91O4txWpS4JuVpisw/Wju8gcQuCZMpkN36rg5BktLF8xg5l/7EK07glNBgN+N6XTCMHJj+5KdBeZQ0qC21tounEXRwoLaOpbjh4MHm2qXJjgMrqL7EHXXP/2KIXrXnPhYhzAZXQXLsYBXEZ34WIcYNwxulpcjMhxk1m4GF8YV4zumdyANmsySmnJ0SbFhYsRxbiyugdb2xB79qMdhUL0LlwcTYwriY6USJfJXYxDjC9Gd+FinMJldBcuxgFcRnfhYhzAZXQXLsYBXEZ34WIcYNRvUxVC9ACbjzYdQCXQcbSJwKXDCpeOaEyWUlZZD44FP/pmu/21Iw0hxHKXDpeO0U6HE1zV3YWLcQCX0V24GAcYC4z+h6NNQBguHdFw6YjGaKHDFqPeGOfChYvMMRYkugsXLjKEy+guXIwDjFpGF0JcKITYLIRoFkLcPMxjNQghXhRCbBRCrBdCfCF8/DtCiN1CiNXhz8Wma24J07ZZCHFBFmlpEUK8Ex5vefhYuRDiWSHE1vDfsuGkQwgx03TPq4UQ3UKI/x6J5yGEuEsIcUAIsc50LOX7F0KcFH6OzUKIXwkh4hdbT46O/xNCbBJCrBVCPCyEKA0fbxRC9Juey++yRUfWIKUcdR9ABbYBUwEfsAaYM4zj1QILwv8XAVuAOcB3gC/btJ8TpikHmBKmVc0SLS1ApeXY/wI3h/+/GfjJcNNh+S32AZNH4nkA7wYWAOsyuX/gbeA0QvVjngQuygId5wOe8P8/MdHRaG5n6ScjOrL1Ga0S/RSgWUq5XUo5CNwHXD5cg0kp90opV4b/7wE2AhPjXHI5cJ+U0i+l3AE0h2keLlwO3B3+/27gihGk41xgm5RyZwL6skKHlPIV4JBN/0nfvxCiFiiWUi6VIW77m+matOmQUj4jpTQS1r8J1MfrIxt0ZAujldEnArtM39uIz3hZgxCiEZgPvBU+9NmwqnaXSWUcTvok8IwQYoUQ4qbwsWop5V4ITUrAhBGgw8B1wL2m7yP9PCD1+58Y/n+46AH4OCEJbWCKEGKVEOJlIcQZJvqGm46kMFoZ3W4dM+x+QCFEIfAQ8N9Sym7gTmAacCKwF/jZCNC3WEq5ALgI+IwQ4t3xSB5GOhBC+IDLgH+FDx2N5xEPTuMO93P5BhAE/hk+tBeYJKWcD/wPcI8Qoni46UgFo5XR24AG0/d6YM9wDiiE8BJi8n9KKf8NIKXcL6XUpJQ68EeG1NFho09KuSf89wDwcHjM/WE10FAHDww3HWFcBKyUUu4P0zTizyOMVO+/jWi1Omv0CCFuAC4BPhRWxwkvHQ6G/19ByFbQNJx0pIrRyujLgBlCiClhqXId8NhwDRa2hP4Z2Cil/LnpeK2p2ZWAYYF9DLhOCJEjhJgCzCBkdMmUjgIhRJHxPyHjz7rweDeEm90APDqcdJjwAUxq+0g/DxNSuv+wet8jhDg1/Nt+xHRN2hBCXAh8DbhMStlnOl4lhFDD/08N07F9uOhIC0fDApik1fNiQtbvbcA3hnmsdxFSqdYCq8Ofi4G/A++Ejz8G1Jqu+UaYts1kyZJKyMuwJvxZb9w3UAE8D2wN/y0fTjrC/eYDB4ES07Fhfx6EJpa9QICQRLwxnfsHFhKaiLYBvyEcBZohHc2EbALGO/K7cNurw7/XGmAlcGm26MjWxw2BdeFiHGC0qu4uXLjIIlxGd+FiHMBldBcuxgFcRnfhYhzAZXQXLsYBXEZ34WIcwGV0Fy7GAf4/AV0oiIWiSVUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:>"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import rioxarray\n",
    "from rioxarray import merge\n",
    "from rasterio.plot import show\n",
    "from pyproj import Transformer\n",
    "import numpy as np\n",
    "from matplotlib import rcParams\n",
    "import matplotlib.pyplot as plt\n",
    "import cartopy.crs as ccrs\n",
    "\n",
    "items = [\n",
    "r'data.2022.040506.tiff',\n",
    "r'data.2022.030405.tiff',\n",
    "r'data.2022.020304.tiff',\n",
    "r'data.2022.010203.tiff']\n",
    "\n",
    "elements = []\n",
    "\n",
    "#for val in items:\n",
    "#    elements.append(rioxarray.open_rasterio(val))\n",
    "\n",
    "\n",
    "for val in items:\n",
    "    rds = rioxarray.open_rasterio(val)\n",
    "    if rds.rio.crs != \"EPSG:4326\":\n",
    "        rds = rds.rio.reproject(\"EPSG:4326\")\n",
    "    elements.append(rds)\n",
    "\n",
    "\n",
    "merged = merge.merge_arrays(elements, nodata=0.0)\n",
    "\n",
    "image = merged.values\n",
    "show(image)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "218665c6-36b8-4442-acbc-0f4e00e2ae08",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block !important;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.DataArray (band: 1, y: 1600, x: 1500)&gt;\n",
       "array([[[-9999.      , -9999.      , -9999.      , ...,    35.593346,\n",
       "            33.86217 ,    32.04727 ],\n",
       "        [-9999.      , -9999.      , -9999.      , ...,    36.565857,\n",
       "            36.04242 ,    33.63808 ],\n",
       "        [-9999.      , -9999.      , -9999.      , ...,    36.836056,\n",
       "            35.167736,    33.3798  ],\n",
       "        ...,\n",
       "        [-9999.      , -9999.      , -9999.      , ..., -9999.      ,\n",
       "         -9999.      , -9999.      ],\n",
       "        [-9999.      , -9999.      , -9999.      , ..., -9999.      ,\n",
       "         -9999.      , -9999.      ],\n",
       "        [-9999.      , -9999.      , -9999.      , ..., -9999.      ,\n",
       "         -9999.      , -9999.      ]]], dtype=float32)\n",
       "Coordinates:\n",
       "  * x            (x) float64 -19.97 -19.92 -19.87 -19.82 ... 54.88 54.93 54.98\n",
       "  * y            (y) float64 39.97 39.92 39.87 39.82 ... -39.88 -39.93 -39.98\n",
       "  * band         (band) int64 1\n",
       "    spatial_ref  int64 0\n",
       "Attributes:\n",
       "    scale_factor:  1.0\n",
       "    add_offset:    0.0\n",
       "    _FillValue:    0.0</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'></div><ul class='xr-dim-list'><li><span class='xr-has-index'>band</span>: 1</li><li><span class='xr-has-index'>y</span>: 1600</li><li><span class='xr-has-index'>x</span>: 1500</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-976081a2-af93-4132-b13a-080673bcc520' class='xr-array-in' type='checkbox' checked><label for='section-976081a2-af93-4132-b13a-080673bcc520' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>-9.999e+03 -9.999e+03 -9.999e+03 ... -9.999e+03 -9.999e+03 -9.999e+03</span></div><div class='xr-array-data'><pre>array([[[-9999.      , -9999.      , -9999.      , ...,    35.593346,\n",
       "            33.86217 ,    32.04727 ],\n",
       "        [-9999.      , -9999.      , -9999.      , ...,    36.565857,\n",
       "            36.04242 ,    33.63808 ],\n",
       "        [-9999.      , -9999.      , -9999.      , ...,    36.836056,\n",
       "            35.167736,    33.3798  ],\n",
       "        ...,\n",
       "        [-9999.      , -9999.      , -9999.      , ..., -9999.      ,\n",
       "         -9999.      , -9999.      ],\n",
       "        [-9999.      , -9999.      , -9999.      , ..., -9999.      ,\n",
       "         -9999.      , -9999.      ],\n",
       "        [-9999.      , -9999.      , -9999.      , ..., -9999.      ,\n",
       "         -9999.      , -9999.      ]]], dtype=float32)</pre></div></div></li><li class='xr-section-item'><input id='section-6889843d-3731-449e-bb3a-a5a8173a88c5' class='xr-section-summary-in' type='checkbox'  checked><label for='section-6889843d-3731-449e-bb3a-a5a8173a88c5' class='xr-section-summary' >Coordinates: <span>(4)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>x</span></div><div class='xr-var-dims'>(x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-19.97 -19.92 ... 54.93 54.98</div><input id='attrs-62094838-ef3b-438e-8f1a-7d9ca1ff2b1c' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-62094838-ef3b-438e-8f1a-7d9ca1ff2b1c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-37e50643-6512-4c1e-8462-72374251cbf1' class='xr-var-data-in' type='checkbox'><label for='data-37e50643-6512-4c1e-8462-72374251cbf1' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([-19.975   , -19.925   , -19.875   , ...,  54.875001,  54.925001,\n",
       "        54.975001])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>y</span></div><div class='xr-var-dims'>(y)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>39.97 39.92 39.87 ... -39.93 -39.98</div><input id='attrs-d65ad0cf-b30e-4a05-b7b3-c4e97ae105d8' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-d65ad0cf-b30e-4a05-b7b3-c4e97ae105d8' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e2438886-bd48-435a-a9af-8165be80eaa1' class='xr-var-data-in' type='checkbox'><label for='data-e2438886-bd48-435a-a9af-8165be80eaa1' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([ 39.975   ,  39.925   ,  39.875   , ..., -39.875001, -39.925001,\n",
       "       -39.975001])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>band</span></div><div class='xr-var-dims'>(band)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>1</div><input id='attrs-ecf762a6-2877-4d11-a195-1e6180e01698' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-ecf762a6-2877-4d11-a195-1e6180e01698' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-53178175-6338-43a4-b387-0de56936c546' class='xr-var-data-in' type='checkbox'><label for='data-53178175-6338-43a4-b387-0de56936c546' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([1])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>spatial_ref</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>0</div><input id='attrs-7bdf2061-1c27-43a4-ad8d-c4b7434fdc7d' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-7bdf2061-1c27-43a4-ad8d-c4b7434fdc7d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-69efdcc6-e2e0-40cf-8e02-00cbcd79ddd9' class='xr-var-data-in' type='checkbox'><label for='data-69efdcc6-e2e0-40cf-8e02-00cbcd79ddd9' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>crs_wkt :</span></dt><dd>GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AXIS[&quot;Latitude&quot;,NORTH],AXIS[&quot;Longitude&quot;,EAST],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]]</dd><dt><span>semi_major_axis :</span></dt><dd>6378137.0</dd><dt><span>semi_minor_axis :</span></dt><dd>6356752.314245179</dd><dt><span>inverse_flattening :</span></dt><dd>298.257223563</dd><dt><span>reference_ellipsoid_name :</span></dt><dd>WGS 84</dd><dt><span>longitude_of_prime_meridian :</span></dt><dd>0.0</dd><dt><span>prime_meridian_name :</span></dt><dd>Greenwich</dd><dt><span>geographic_crs_name :</span></dt><dd>WGS 84</dd><dt><span>grid_mapping_name :</span></dt><dd>latitude_longitude</dd><dt><span>spatial_ref :</span></dt><dd>GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AXIS[&quot;Latitude&quot;,NORTH],AXIS[&quot;Longitude&quot;,EAST],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]]</dd><dt><span>GeoTransform :</span></dt><dd>-20.0 0.05000000074505806 0.0 40.0 0.0 -0.05000000074505806</dd></dl></div><div class='xr-var-data'><pre>array(0)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-e2328d0e-47d9-4e48-85c1-d496e3cea522' class='xr-section-summary-in' type='checkbox'  checked><label for='section-e2328d0e-47d9-4e48-85c1-d496e3cea522' class='xr-section-summary' >Attributes: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>scale_factor :</span></dt><dd>1.0</dd><dt><span>add_offset :</span></dt><dd>0.0</dd><dt><span>_FillValue :</span></dt><dd>0.0</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray (band: 1, y: 1600, x: 1500)>\n",
       "array([[[-9999.      , -9999.      , -9999.      , ...,    35.593346,\n",
       "            33.86217 ,    32.04727 ],\n",
       "        [-9999.      , -9999.      , -9999.      , ...,    36.565857,\n",
       "            36.04242 ,    33.63808 ],\n",
       "        [-9999.      , -9999.      , -9999.      , ...,    36.836056,\n",
       "            35.167736,    33.3798  ],\n",
       "        ...,\n",
       "        [-9999.      , -9999.      , -9999.      , ..., -9999.      ,\n",
       "         -9999.      , -9999.      ],\n",
       "        [-9999.      , -9999.      , -9999.      , ..., -9999.      ,\n",
       "         -9999.      , -9999.      ],\n",
       "        [-9999.      , -9999.      , -9999.      , ..., -9999.      ,\n",
       "         -9999.      , -9999.      ]]], dtype=float32)\n",
       "Coordinates:\n",
       "  * x            (x) float64 -19.97 -19.92 -19.87 -19.82 ... 54.88 54.93 54.98\n",
       "  * y            (y) float64 39.97 39.92 39.87 39.82 ... -39.88 -39.93 -39.98\n",
       "  * band         (band) int64 1\n",
       "    spatial_ref  int64 0\n",
       "Attributes:\n",
       "    scale_factor:  1.0\n",
       "    add_offset:    0.0\n",
       "    _FillValue:    0.0"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "merged"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "88bacc33-693d-47f8-a66b-09b653c0a92c",
   "metadata": {},
   "outputs": [],
   "source": [
    "merged.rio.to_raster('merged_africa_n4.tiff')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "879746b1-14af-4400-9ac6-ca83f26556df",
   "metadata": {},
   "outputs": [],
   "source": [
    "#raster_array = merged.astype('float64')\n",
    "import rioxarray as rxr\n",
    "import cartopy.feature as cfeature\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "9a1ebe72-5da0-40d7-ad08-366864446810",
   "metadata": {},
   "outputs": [],
   "source": [
    "raster = rxr.open_rasterio('merged_africa_n4.tiff', masked=True)\n",
    "raster_array = raster.astype('float64')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "914d26b0-557a-4872-9db5-72ab0b1dbd3d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block !important;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.DataArray (band: 1, y: 1600, x: 1500)&gt;\n",
       "array([[[-9999.        , -9999.        , -9999.        , ...,\n",
       "            35.59334564,    33.86217117,    32.04727173],\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "            36.56585693,    36.04241943,    33.6380806 ],\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "            36.83605576,    35.16773605,    33.37979889],\n",
       "        ...,\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "         -9999.        , -9999.        , -9999.        ],\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "         -9999.        , -9999.        , -9999.        ],\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "         -9999.        , -9999.        , -9999.        ]]])\n",
       "Coordinates:\n",
       "  * band         (band) int64 1\n",
       "  * x            (x) float64 -19.97 -19.92 -19.87 -19.82 ... 54.88 54.93 54.98\n",
       "  * y            (y) float64 39.97 39.92 39.87 39.82 ... -39.88 -39.93 -39.98\n",
       "    spatial_ref  int64 0\n",
       "Attributes:\n",
       "    scale_factor:  1.0\n",
       "    add_offset:    0.0</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'></div><ul class='xr-dim-list'><li><span class='xr-has-index'>band</span>: 1</li><li><span class='xr-has-index'>y</span>: 1600</li><li><span class='xr-has-index'>x</span>: 1500</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-aeffaed5-649f-47f1-bf76-7a7cc61617ac' class='xr-array-in' type='checkbox' checked><label for='section-aeffaed5-649f-47f1-bf76-7a7cc61617ac' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>-9.999e+03 -9.999e+03 -9.999e+03 ... -9.999e+03 -9.999e+03 -9.999e+03</span></div><div class='xr-array-data'><pre>array([[[-9999.        , -9999.        , -9999.        , ...,\n",
       "            35.59334564,    33.86217117,    32.04727173],\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "            36.56585693,    36.04241943,    33.6380806 ],\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "            36.83605576,    35.16773605,    33.37979889],\n",
       "        ...,\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "         -9999.        , -9999.        , -9999.        ],\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "         -9999.        , -9999.        , -9999.        ],\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "         -9999.        , -9999.        , -9999.        ]]])</pre></div></div></li><li class='xr-section-item'><input id='section-041ad232-5117-4222-9513-e1542704b9f8' class='xr-section-summary-in' type='checkbox'  checked><label for='section-041ad232-5117-4222-9513-e1542704b9f8' class='xr-section-summary' >Coordinates: <span>(4)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>band</span></div><div class='xr-var-dims'>(band)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>1</div><input id='attrs-753152f7-7a9b-4de8-aca5-fddb6c559ac9' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-753152f7-7a9b-4de8-aca5-fddb6c559ac9' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8452ecc1-f973-4f66-80a9-d5106430c5a3' class='xr-var-data-in' type='checkbox'><label for='data-8452ecc1-f973-4f66-80a9-d5106430c5a3' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([1])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>x</span></div><div class='xr-var-dims'>(x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-19.97 -19.92 ... 54.93 54.98</div><input id='attrs-a8634d1c-3a01-46ab-befd-5fc62be114ce' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a8634d1c-3a01-46ab-befd-5fc62be114ce' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6a73ed12-a787-496d-bbc8-7ddde8e3b1ae' class='xr-var-data-in' type='checkbox'><label for='data-6a73ed12-a787-496d-bbc8-7ddde8e3b1ae' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([-19.975   , -19.925   , -19.875   , ...,  54.875001,  54.925001,\n",
       "        54.975001])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>y</span></div><div class='xr-var-dims'>(y)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>39.97 39.92 39.87 ... -39.93 -39.98</div><input id='attrs-1608cec1-5424-48c0-8abc-40ce5521fcbf' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-1608cec1-5424-48c0-8abc-40ce5521fcbf' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-82e11985-f8d4-400c-a0ba-73c423a89b6d' class='xr-var-data-in' type='checkbox'><label for='data-82e11985-f8d4-400c-a0ba-73c423a89b6d' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([ 39.975   ,  39.925   ,  39.875   , ..., -39.875001, -39.925001,\n",
       "       -39.975001])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>spatial_ref</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>0</div><input id='attrs-817f4f8f-24ba-4bb5-8a30-ac1788a6589d' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-817f4f8f-24ba-4bb5-8a30-ac1788a6589d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c163357d-eb88-4e47-bae6-a2639a8da52b' class='xr-var-data-in' type='checkbox'><label for='data-c163357d-eb88-4e47-bae6-a2639a8da52b' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>crs_wkt :</span></dt><dd>GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AXIS[&quot;Latitude&quot;,NORTH],AXIS[&quot;Longitude&quot;,EAST],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]]</dd><dt><span>semi_major_axis :</span></dt><dd>6378137.0</dd><dt><span>semi_minor_axis :</span></dt><dd>6356752.314245179</dd><dt><span>inverse_flattening :</span></dt><dd>298.257223563</dd><dt><span>reference_ellipsoid_name :</span></dt><dd>WGS 84</dd><dt><span>longitude_of_prime_meridian :</span></dt><dd>0.0</dd><dt><span>prime_meridian_name :</span></dt><dd>Greenwich</dd><dt><span>geographic_crs_name :</span></dt><dd>WGS 84</dd><dt><span>grid_mapping_name :</span></dt><dd>latitude_longitude</dd><dt><span>spatial_ref :</span></dt><dd>GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AXIS[&quot;Latitude&quot;,NORTH],AXIS[&quot;Longitude&quot;,EAST],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]]</dd><dt><span>GeoTransform :</span></dt><dd>-20.0 0.05000000074505806 0.0 40.0 0.0 -0.05000000074505806</dd></dl></div><div class='xr-var-data'><pre>array(0)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-37a9df6a-88fe-440c-a7e1-8665d00a18c6' class='xr-section-summary-in' type='checkbox'  checked><label for='section-37a9df6a-88fe-440c-a7e1-8665d00a18c6' class='xr-section-summary' >Attributes: <span>(2)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>scale_factor :</span></dt><dd>1.0</dd><dt><span>add_offset :</span></dt><dd>0.0</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray (band: 1, y: 1600, x: 1500)>\n",
       "array([[[-9999.        , -9999.        , -9999.        , ...,\n",
       "            35.59334564,    33.86217117,    32.04727173],\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "            36.56585693,    36.04241943,    33.6380806 ],\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "            36.83605576,    35.16773605,    33.37979889],\n",
       "        ...,\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "         -9999.        , -9999.        , -9999.        ],\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "         -9999.        , -9999.        , -9999.        ],\n",
       "        [-9999.        , -9999.        , -9999.        , ...,\n",
       "         -9999.        , -9999.        , -9999.        ]]])\n",
       "Coordinates:\n",
       "  * band         (band) int64 1\n",
       "  * x            (x) float64 -19.97 -19.92 -19.87 -19.82 ... 54.88 54.93 54.98\n",
       "  * y            (y) float64 39.97 39.92 39.87 39.82 ... -39.88 -39.93 -39.98\n",
       "    spatial_ref  int64 0\n",
       "Attributes:\n",
       "    scale_factor:  1.0\n",
       "    add_offset:    0.0"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "raster_array"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "306fb33a-3fd0-416d-b63d-53c57eb58f80",
   "metadata": {},
   "outputs": [],
   "source": [
    "x, y = np.meshgrid(raster[\"x\"], raster[\"y\"])\n",
    "crs = ccrs.UTM(\"18\")\n",
    "\n",
    "transformer = Transformer.from_crs(raster.rio.crs, \"EPSG:4326\", always_xy=True)\n",
    "lon, lat = transformer.transform(x, y)\n",
    "raster.coords[\"lon\"] = ((\"y\", \"x\"), lon)\n",
    "raster.coords[\"lat\"] = ((\"y\", \"x\"), lat)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "id": "0cf0ed31-4296-4ef6-b669-6ec593485614",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/ozlem/anaconda3/envs/geocat/lib/python3.10/site-packages/cartopy/mpl/gridliner.py:451: UserWarning: The .xlabels_top attribute is deprecated. Please use .top_labels to toggle visibility instead.\n",
      "  warnings.warn('The .xlabels_top attribute is deprecated. Please '\n",
      "/home/ozlem/anaconda3/envs/geocat/lib/python3.10/site-packages/cartopy/mpl/gridliner.py:487: UserWarning: The .ylabels_right attribute is deprecated. Please use .right_labels to toggle visibility instead.\n",
      "  warnings.warn('The .ylabels_right attribute is deprecated. Please '\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgkAAAGoCAYAAADMynYuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAB7i0lEQVR4nO29ebw8RXX3/z73cr+sAorLl4gG1yRiEozEGA0GjAuiBlCikiCYmBi35DGRJ/FiHn/kMXqJwS1xSYx5FEQRRQkKGJcYzOIKCSq4sqlfWRTZl+925/z+6K47NTVVvU3PTM/Med/XvO5MdXd1dU9P16fPOXVKVBXDMAzDMIyQpWk3wDAMwzCMbmIiwTAMwzCMKCYSDMMwDMOIYiLBMAzDMIwoJhIMwzAMw4hiIsEwDMMwjCgmEkZARN4rIioi17Sx3jRp2kYReUG+nYrIgWNo12Fe/Ye1UJ+IyKki8kMR6eX1HlyyzTX5eu8ddf9Ghogc6H2vLxh1PQNE5BR3rqbdFmN+MJEwGa4EvgT8z7QbMk+IyKPym+Kza2x2FPDnwE8BV5F9L3eWbPM/+XpXNmpoA0TkXiJynoh8X0TuFpHbROSbIvJXIrLbpNoxKtPuuERkWURuFJG/m8b+DWPW2WXaDVgEVPW1wGvL1hORTaq6fQJNmheOBrYC/1Jjm4O897+oqkmB4L4PVT2mYftGYW/gCOB7wOXA/YGfBV4N3Bt48RTaNIs8AdgP+OdpNcB+18YsY5aEEnLz9ItF5BIRuUtE7sjfPy5Y7zdF5BsicqeI/LuI/Iy3bMiU75mw3ycibxSRG8meVvHMq6eKyDtE5GYRuUVE3iYim7w6jheR/86fMu8Uke+IyAdFZJ+SY9pHRD6QH8v1IvJ/AEkc+8tE5Kv50+ytIvIxEXlESf1vyNu/VUSekZf9jYhcnh/HDhG5VkROF5H9g21fnD893yUiHyfrHFMcDXxGVe8UkZ8WkU+IyA/ytt4tIpeJyCtERPK6LwL+ytv+DveEKyIX5W2+SET+XESuBW7Ilw25G0TkHvlxXiEi20TkJhH5jIjcK19+kohcmpfvEJEfichHReThRefO4wfAXqr6cFU9BHgAcHW+7PFlG3ttPkNE/jq/hq7Lz++9ReSc/Jr5tog8M9j2kXlbbxSR7SJytYicJiJ7eev45+tl+f5uF5HzRWSzd77/P2+blNtgHxF5T34d/1BE/qLguH7Wq+cIr/yJXvnPepscDdwMfE4GXRcn5cd4V36d/Jpklqkv5+flP8X7Def7eIqIfDZv590i8iX/3AX1/5mI/LOI3AW8IV/+ayLyP/nv4n/yz279U7x6NovIu/NzsV1EvifZvWBXb51Nkt0PbsmvsbcCmzCMtlFVexW8gL8DNH/dBFxG9vT6AuC9efm2vOybQC8v+y+vDrfeNV7ZNd6224CvA/+RL3P72wr8mKxzcGVvyNf5BW9fVwBfA27NPx9Qckwf8ur7DnALcEekjf6xfwO4Pn9/C/DgfJ0XeOscCJzitf1pXl2X5dt9PThPX/bWOdKr6ydk5v07vLLDvHUPzMtemH8+JP/8A+C/yTp4t93L8nXeAWzxyr8IfDFfdpH3fezMj/eq4Lt6b/55E3CJV8/3gW8D68CB+Trn523/Rn7MO711d6tx/Z0BfBm41tvfOyts59rsrqHr8s+9vE1byDpPBW4D7pVv93PA7Xm5a/96/vm/gKXgfG0H7ia7jlz73l90voGne9+fAnflx/djr+zJwfeswAvyss/knz/sHe873D6C83A18L5IXVvJ3E3u+roh3/+382MKf8PH0r9mfwB81zufx0bq30b2e/w6cBpwP++83p2f19u89U/J69jP++7uAL6a16XAx732vMHb9mrgR96x6LTvm/aan9fUG9DlV/6jdzeG88hv7sA9gQfR7/wVeGa+7E1e2e55mVvvGq9udyPYBvxCXrac/3fbfwvYk+wp/6PeDWYP4Nn0O3l3414CfgXYo+CYHuzV/+a87L5knfJGG4Njf1FetitZZ6/AP+ZlL/Dqc6LiLuApwX5/wbUz//z73nYPycs+59oA7JOXvc9b7zBv+z8h67zum3/el7yD9s6Fq+8/vPJTYjdS+p2eAkcE34f7rt6bfz7BW3fVq+NhwD3y9wcBK96yJ3nb/EaNa/BSbzslEw0rFbZzbf5Rfm4e6tVxWf5d/kbkmE/PP98J/HRe9mJvvWcG52sdODgvc9fo9RXO94FenV8gE173pt9BnxpZ7wV52TH0fzv3zr9rJ4Je7O3jUXnZsyN1fZLsd+Vfh+6afq1X5n7DV+Wf3w9IXvaPedl3I/V/E9jXXUfA/6UvKg7Jy//QW/+UvOw19EXy/nnZ4731Hk/2+787//zR/Dj2JLtfmEiwV6svczcU88v0zfBvUtWtAKp6s6pe7a13q6p+PH//Da/8vhX28W+q+rW83vVg2QWqeqeqKvDhvGw34CFkT3U3k3VMN4nIV8huWvdS1bsAROSLweuXgEd69X8o3++PyG76qWP/h9wsv5W+T/+xkWN5ef7/+ar6qWDZLwJfkczFoXlbHT+V///5/P8nVfVWv40RjgY+n7cdYAfwZ7lpdgdZ5/WEoP4qfFtV/wWi34fjV7x9/o0rVNXvqurt+ccHAv+Wm6Z7wKe97Su3R1UPBnYHfp2sI3w+8H/c8sR37POfqnoLmWhwfEpVt5F1fI775f9/Of//X6r6vfz9B7z1Dgnq/7qqXpq/d9d+leve52zNYj9uJBM1fntifIzsaX4TcDxZx7mZTDSc7a13NOmYlQvy39U1Xpn7Dfvn5b4ich+yhwKA3wZ6+TX8+3nZQ0Vkv6D+0/Pz7q4j97u7QlUvzt+fFWmXu7buBVyb7+c/veWPJRN8Lnj1w5pxJ3BBpD7DGAkLXGyHW7z3O733Q37+CNdX3MdAXap6vYgcRNZpHEJ2E/o94PdE5Fmqei79G45j76r1B5+/Snaz9bk2UscdwF7AK0XkE55Y+TWyJ1Qhe0L6Rr7ez+XbLddoF/kN+fFkoxQcb6F/0/4umWvoIWRPmkX1h1T9PpKIyIPJAuU2kZmYLyH7rR2cr1KnPeTi9N9F5INkFpSTReTU/PyWfce35XXsFJGBMrKnzo1mh7ut2LxbvPfu2q9y3TeuQ1XXReRdZE/8LyR7ggc4T1Vv9lY9mjxmJVKNOwc7I2XhefHbcjV9IeOzEnxOXUdl59Xt6w6ygNWQWypsaxitYZaEYr5C/0f9Chc4JFng34ET2P+RIrKHZHf3Z+VlW4ErReSngPuo6htU9Tmq+gj6Q/SeCKCqErwuYvDG81v58dyH7EnVxz/2s1T1se4FvIzMzxryfLIb7a8C50o/yPJX6N/Afl5VH0NmNg+5LP//FBFxnV1seONvknW0/+yVOcvGp1T14cBhwA8j27bBl/L/K8CfukIReUge3Pco+kFkT1XVXwb+us4ORORIEfl57/M96H9Hy+RPkonveBS+kv9/vIj8dP7+t73lF1OPu9wbEdlzlIYF/COZa+KR9MXh6d6+DiRzcf3zqDvKrVXX5B8vAw71fgvPAdZUtUxcfj3//1AR+cX8/XGR9b7sdgsc7+3ncLLf3EfIYpCcaH+2ZOwBPK3moRlGKSYSClDVa4C35x+PJjP/fY3M7HvYBJrwALKb01X0O8u35U+QjwC+KlnU/KUichXZkzNkQYxRVPVKshsNwJ+IyLfJnrz3DNa7Gvj7/OOpuRn/qyJyE1lH8ZRI9ZeSiZnt+fL3i8hS0J6vi8g3gf8d2f4N+f8DgatE5ErguZH1jgYuy4/F4fbxlPyYfkB2/sbBB8mCIwH+Oj833yILers3mRBzrop/EZGvk8Vr1OEpwNckG31yKZnlxrkSzlPVm0Y5gAJOJXuK3RO4XEQup/8b+Dz1Tdrf8t5fnrtEHjxqI1X1BvrX8Z5kT+6f9FY5hsz//3Ha4VX5/2cC1+WjE64l+33+SYXt3052XpeAL+TnNSa030Z27d4D+IaIfE1EvkvmWvwwWZzDXfS/k2eT3R+uoe8SMYzWMJFQzh8DLyVLqLMHWeDfN8k61nHzt2R+y33IntDfCbjhYVfly24hi0u4L5kZ/1XAu0vq/X2yju4uMt/nO4n7/l9Odvxfzet/EJlAeif9G/QAqvqvwIlkT0LHAv+gqp8mcw1cS+Zf/xbwksi25+f7/CHZjf/KcL38ienJDD8h/ilZcOkdZDfYv6G9DiJs53Yykfg3ZN/D/mTn53PAbar6LTLXz9VkFoUbiT81FvE54D/IfqMHkVli/ofs+3/eyAeRQFW/SW4JIvPxP5ys03ojmVWkV7PK88me+n8C/DR5YG1LzX279/7MIIbkaAZjVkZCVc8me1L/LNl3+nNkT/MfJt7Zh9v/KN/+q2SWoJ0Mfo935+vdSGYVezeZW+PnyFxIXwFOJh+WS5Yv451k94V7kv0e3tr8CA0jjovSNTqE9LPT/aWqnjLNtnQNETmGLKL7EFW9ZNrtMaZH7ia7nlxIqeo38vL9yDrTP1fVN06xiQOIyMNV9Tve5+fTd7s9NRLsaxhTxwIXjVnjTuAvZl0giMjv0/elh7xbVcusQQuLiNyb7Kn50WQC4QInEHL2I0ua9cEpNK+ID0mWUvvbZG10CdkuYnD0i2F0BhMJxkyRP23NwxPXAQyPTHDUSTO9iOxFFky5lexaeKG/MH9aP2XyzSrlE2TBwi6e5xtkbr6/UTPpGh3F3A2GYRiGYUSxwEXDMAzDMKLMvLtBZnTu9L333pvbbrttoOzAzXsjLeVD0UjOlrbqll1246otN5SvGHDg5iz1wTXX3zbw2S8blQdt3idarijXXH9b9Ly3TaoNV19/a7TcnQdBouu460J22Q3duXVj3dh3nNqna5PbJjzfqe/CPxZ/f3WvJb/ts8Q42+1/h6nzWeU7Tq0vu+yO7ry7eQNLuOb62z6pqkeUrzmbHHHEEXrjjTe2Utcll1wys+dq5t0NsyoSTjvtNE466aSNz6evHslyDcPOOoMj0apu67ars6+Q6/c7bKDtZZy5+oyh/Z+4dmHj/Vfdl886PU5cu3DovI+Ts1aPyvedjcw7fu38wvVPXz1y6Lz4x7PHQUez7fLyUZ3rrA98Ltuvv3/ofzfu8zJLA9eb+1znGtr9oKO4+/LzKq8fY5nloWNz5Snc+kXrpLZbZpldD3pmpXNep07/GFy7lljaeO8v7+XnPXbcg3UPj0zd86BjuPPyc0dud4oT1y68RLMZSueSQw45RC++uG7usDgiMrPnauYtCYtI7Ibg37RTN/DYdnU5fu18TjvtsMrrn7565NB+xyUQihhFFLVF2AmHxMpdB3/m6jNQdKiDWWJpoyOB8s6kqG2pdoXf3/pGx1VPKIxK6thina5PXYHQdJtR6G2I9+H9hsc96fO+yPRm/CG6DexK6wAxK8I6vYGXX5aibDkwICSqrD8qJ65duPFynxeF9fzPMcqN/fi186OuiF6F7y+0rjix4jPL30uVDj0mMMKypfzPvR8Hwx3+oNVgBzsGPsPw8ZlAmAzZlJ3t/M0ydrV1iKKOu63OvKmboiopcz+MvyOK7fv4tfPHLoRSpMz8sU66SZ2hCEmxTo/TV4/ceFUhXK+qlWFSuI4zNN+vB3+uzF8e1uEY9818OfHn4yxDs96xGPODuRumjLMitHmzrVJf032eufqM5PR2Z64+o7L/exz4pnmfaT4lH792fqFwaqNO3+8+3Gn2v+PwPBS5P1LrdoWYAGiyTUoouP91AwfDfflUdWG42IVw26auJKMpyqzH7LWBiYQpc+Laha3dgH1XQmqZYxS/Zip63z2xTtt0PQ6hUuW4wnVCcdBWu9LBmbHOpS8Gq1oHfLomDupQpWOtZokZPq9N6vPr8YMY/bqdu6HOPo3x0WsoEucJEwkdoC2hkIpAryIaqlIUuDhtcTBOmhzbOAPM0sGpw0FuTahzPXY1kK7uk3fT0RN1OvNQKFRtp28tMouCMUlMJHSEcNhZU0JrQvjfX6cu03QlzAInrl0Yfcp3wy/bwK8//D7csuPXzp+YBaBtV9k0qTJ6wieMKajaeTdZzx8qGatjnr6HrqBAr/akp/OHiYQ5wu/8i1wPdZ/82uzkZo2ikQDOvVCW46LNcxcTAGWfU7QVk9BFK8IkGOWJvq5FoIrFoug3bzTD3A0mEjpPLPdBbISCGxXh37BTN4sqN/UiV8WikOoY/XL33j9fkxRUfsIjR+o68d1afhubxpF0QRw0eZKfFEVCoElbq4xmWeTfqzEeTCR0iFin5P/oUzeAKuvE1k/d5GMuikWkzaDSNvE79LptjLm1isRDuK1f1lYsQqye1DU6LISHsxhOE78tVdoWO/ZY5kVjCqhaMiUsT0KnmLZJP7RGhImQFpE652BabpkT1y7c+O78toSEwsD9T1kSqhx7G0IytH7VFb1VmYTodcIglgMhng1yMLmZq6Nq/otYfV2w8MwDWTIlbeU1y5gloUNM4qk1vCnFTNPzFKAYntM9E8M3q1AluHTaQ0BTT/dFbaoiGMJ1Y26OUSmylLn9LHtzHGTLqnWoVfbd1rGE80U4a0IYiBhm46wrYkwMGJPAREKHiN2cUzftNkZBTCItc1vEovrDMpfMKTahVIyyuRRSFJnku0DVuTJ8UVAUf+GsCak8C3WTRaU65KJ4mkFBW08UVLVMVBEKdUYJFbkbwmWhi29UAbDMUktzvqbrXwQsmZKJhM5TdoOvS+wmmfJJt0msI/GTycRmLkx1PrFyV1ZFIJSdtzrWgPB7mLZrJvwuXSrmstEKMcFQ5mrwgzZ9Up19bAbQmEWizKJQFpDrRoA0eToP19cK7Ul9zsqaWTmKxEJMqKS+gyrntSqLIgwcmbvBMJEw59Tp7EalSfrh1LCuUVMZVx2dUWX4XxnTFgZhG6oIyKLjDl0KKXdESliGHX+V4aFVBW9MaLjt/ZiQSXdoTTrgUSwiYVmxeIrNLDmaa2NWLJCjYoGLJhJmmlRnMI0heEUUmVzHQdXRGV3o3MdFKBqq4ouA2DDJKmIkNgpCvWWp9lYVCkUCp+15MsoocovUIWURqNoZl63nz3C5gx2lFopY2+q0x5gfFst+NIe4G7c/ImESODN2lafFcQ7lsptWMW1ZRerGxbgRF+79NdffWtqWMvdGlVEmg7Nkdv/aGNesrE6YZXFH6/gTV6VmoUyNjHCxS7NwPtumlw+DHPU1y5glYU5wfthxRte7ug/cvA93/qRfXnRjqzqGvc3ocn/eClf3uBh3LEfXCK+vKiMh6tLGuZxmhzZq8GHb1204kmmZ5Y0pqWNpni1zY4bOwfDFNjCRMAecufqMsf2gUyl/q01HXTQxTr100amo8lg7qppH2xBUiyIOYnTx2Ls02iQmVIcTJxUJ7HYznvpuOH+kSD9502KLAiOOuRvmgPDH3daNctR6isawF4mB2Cu2Xnq/1TNGVhnpMC84S5N7xZZ3nVTbfbrU2cWSQ5Vd25Ns1zo9trOjU+esS/S0ndcsY5aEGSeMJHc/9vApOTSLl5nJx30jbrq9HwPR1o2tyKLQxafltigSCuEkVo4uBMWWuTiaDn+cFikL2qy0f15RtdENYCJh7vBvjlWeFqfx9DjKza/q0MY6zLMQ8KkjCKcVa9HW9ThrHax/XS9aPgKj25hImGFSqXEnGbg3DebteCZNLP9BuLyLwqmLbWqTWMyAMV3MjmAiYS6IBe61kda1SmdcZZ2ujrPu6iyP46bKUMRpUdUdNu906XeyyMx6PEEbmEiYUYo6tzZuMOO4SXXpxpd6ija6QdPvZRFFn2GMExMJM8is3Qi7JA7AhIFhGOUoFrgIJhJmjlkTCIZhGLOKuRssT8LMUGV8eIo2YhTmCXceRzmnRvew79Iw2sd6jo5TpyNzYsDPJTCvIxxGJTaboTG72HdojANt6TXLmLuho9S96RWlezWBEMc6FsMwUmQxCdNuxfRZeJHQtTHhdTquKi4EEwjV8M97l64HoxwTe4YxPhZeJMxyh5Ca9Ci23KhObIx+18ftm8gxjJaZg3kX2mDhREJXb6ajBCWGmDhohypprad5DaWuma5e4+PArAjGOLERkAsoErp40xzlRhdL4TpLk9vMOpO2MNS9VrpuATEMo9ssnEiYN8yS0A0m8fR+4OZ9uPMnY6naMIwAC1zMMJEw44SWBBMI88moZvUuuUkMY1Ywd4OJhKnThk/VhIFRl9h1N2vCweIRjLFigYuAJVOaKk1ucpY5sfvMeuflEniN8zjatowYhjEezJIwYzj3glkPuk3b+TdOXLuQc848prX6Yox7NEdYV9O6TSAYk8LcDSYSOoufTjkUBSYQZoNUZ9ako+1Kx9ikYz999Uj2rBB0mcpPEctXYRjjRjGRACYSOo+5F+aPJiMhJmFJaEKRaGjaoZdZHAzDmBwmEjqKWQsWg3lLfGQdujFPWOCiBS4aRmcoyqBona9hTB7Vdl5liMjPiMil3us2EXmFiJwiIj/0yo/0tlkVkStE5Nsi8lSv/NEi8vV82d+KiIxyDkwkTJF5eHI02iUUBCYODGP+UdVvq+rBqnow8GjgLuDcfPGb3TJVvRBARB4BPA84CDgCeIeILOfrvxN4EfCw/HXEKG0zkWAYHcSsB0YT3IOHPYC0gIKqtPKqyW8AV6rq9wrWOQr4oKpuU9WrgSuAx4jI/sDeqvoFVVXgDODoBke/gcUkTBHrBAzDaANfFLj311x/K/eZVoPmgJZHN9xbRC72Pr9LVd+VWPd5wFne55eLyAnAxcArVfVm4P7AF711tuRlO/L3YXljzJIwJQ7cvM+0m2AsOCe8/gJOeP0FU29D+Nmegqtz4tqFdr5mgxtV9RDvFRUIIrIJ+E3gw3nRO4GHAAcD1wFvdKtGNteC8saYJcEwjIkQCoIzTn76Rrl77/4b5VQRB24ds1o2Qyc/yOxpwH+r6g0A7j+AiPwjcH7+cQvwAG+7A4Br8/IDIuWNMZEwBewHa3SBtjvkIqvERRddVHl9l0Cpyu+k6nqLStuZPxeNKSRTOg7P1SAi+6vqdfnHY4DL8vcfAz4gIm8CfoosQPHLqrouIreLyGOBLwEnAH83SoPM3dCAUW5KdkMz5o1R3RYpsVI1m6ORxgIZZwcR2QN4MvBRr/gN+XDGrwGHA38CoKqXAx8CvgH8C/AyVV3Pt3kJ8G6yYMYrgU+M0q6FsCS0qaab1mU3M2MeGEcMg+9ucMR+Z27f5pIwJkI+umFiu1O9C9gvKHt+wfqvA14XKb8YeGRb7TJLQg2advRqCcCNGWfcQY5V6j7j5KebQPCoez8ya0J9tNfOa5ZZCJGQ+nHUGYs+Sh56u7EZs8o0R0A4a0Jq/9bpmYXSGD8L4W7wqRMUFcP8pMYiMA1hkHI7iAiqGl1m1MMCPWtiRuDFsCTEiAX0OMtC7EdUNh65aFvDMKoREyeqiohsLPPXMWuCiaVxobSTbXGScQ3jYO4tCS4ewM1xEctMFlLnxmM/UGPWmXZCpZDQonDGyU/fEArZ8uHfqD0hVyc8d3beEujsxxO0wdyLhDoTYNV9Khn3jytmfjWMtuiaOPDxRzK4985VGP6mzZrQvxc1OReWcMkoYmHdDYaxqHQhHXMVwiyMZR2giYUM6+xbRFt6zTAmEhoyiR+iWRGMtnDCYBbEgSOMQVBV6wBLGOX8mMgaxoZALoC7wTAWmVkSBTHqtt9iEzJOXz2SHzfYLgzkNgyzJDRgXD8e/2mvqzf3E15/Afvd/6ED7etyexeVef1OqhyTPRFnjDrT7MKfRwVU2nnNMCYSajJpdd2lG31sWt+YWOhSmxeReTz/dV1vC9/BtURqNNiinF9zN5i7oTKxZC5t09UYhCadjo3MmA7zKBBgfo9r3LQxb01suKmxOJhIqMA0UytPu6OtenOedjsXgf3u/1Ce/TvD34c794vSkfo5E4qw+IT+OQiHOTbt6BdOIMz4yIQ2MJEwo3StQzDLwfTo2rUwLnwx1ObMrvNO7DzZ+auAAjPuKmgDi0koYdGfRKriJ70xxsNPfnjFxkyI/msRCI+z6rW26B2hu3+VpYy3+5yRwiwJHSd8Qm+jI267Qw+T3hjjZ5EFmbvOqrodFp0qAmDRxVQSsySYSCiiK+o6Zspv8mTVZidugmCynPD6C7jooouiMQmLyhknP72S2dzSDhtNUYtJMJFQRJcCn5xQqCoOijrxJk+hzvpwxslPZ8+Djqm9vWGMgsW8GMZ0MJFQQBcEQtmNse6Ns4moaLovY3QW2a3gY9fe+KkazDjqCImZwQIXAQtcNHLCPPmG0RXKBEJVMT/3nVoNYlbSKimZu/DgNFF6Lb1mGBMJEcoigWeZKm4IEwrTxzJX1sv/MK+/13EQ5k5w+OcwNWzSWDzM3RBgPwRj2iy6OIDBEThVXQ0WxFgNf1gkxC0IXYrHmioWuGgiwcd+FOb7nTaLLhBC60Hd67GqX906wT6x82DnBotJyDGRkLMoP4pRAher1GsiYzQWLcVyGXVHNdSJOzCrglGGmEgwkbDoWO4Eo0uYODKMbmEiYYEY59O+CYT2Cc9plp/isKm0ZZ4x14ORxCwJJhJgccyN1pHPDrFU3BdddNGUWjMZYunC7Zo1poZigYvYEEjD6CyLNgyyDYHQVPBbDgXDiLPwloRFsSIYxixTx1W2MBkBjfHTM1OCWRIMo6Msqqm9aH6SOpaVug8AJiqMEOm185plFtaSYBYEYxbwO8wHPuV/Fy4PmVVXRVnmz7pWhaqdv90TDGMYsyQYxozwkx9eMe0mjJ061pOqMRvW+RuNcMmUFnzuhoW0JNhNw5g1Tnj9BRujG/yES+FT9SxaD5q6VRbVHWNMCkUsJsEsCYYxC5xx8tP5yQ+vGOgYw2GSsygQmjIOgWAxCYYxzEJaEgzD6B5FImdUUVAnNsEwNjBDglkSDGOWOePkp8+82b2KFaQNK4m5GY1aKEhPW3nNMqUiQURWReQrInKbiPxYRD4uIo8M1hEROUVErhWRu0XkIhE5KFjnYSLyuXydt4jIkrfsIhFRETk+2OYFInLHqAcZskhmWcPoMqHAiQmeSQohszYYxiBV3A2HAe8AvgII8H+Bz4jII1T1pnydPwNeCbwA+DbwGuDTIvIzqnp7vs7bgQ8D/wGcCjwP+IC3n63AX4nIh1V12ygHVYaIAPZkYRjTJhTsdWd9jNUzqqCwuRyMDWbcCtAGpSJBVZ/qfxaR5wO3Ao8HPi5Zj/sK4FRV/Ui+zonAj4DfBv4h3/SewCXAZcA1wL7Brs4GjgBeBrypycHUxW4Gxqwzr1axIrEwr8dsdAvJ3Q0T25/INcDtwDqwU1UPEZF7kfWNB5L1m89R1Zvz9VeBF+br/7GqfjIvfzTwXmB34ELgf6lq4wNpEpNwj3y7m/PPDwI2A59yK6jq3cC/A4/ztns1WYO3AgcBpwf13kFmpXi1iOzboF1Gh3F+Z7vBt8e8n8vU9VLkomjLLWFuB2NKHK6qB6vqIfnnVwH/qqoPA/41/4yIPILMGn8Q2cP1O0RkOd/mncCLgIflryNGaVATkfBW4FLgC/nnzfn/G4L1bvCWoaqfAu4HHKCqT1DVOyN1vwv4CfmJmAR2M5g88965TYpZD1hMEXNBuP+x3BDufdn5qHvdnbh2Iddcf2utbYw5Y/rJlI6i/0B9OnC0V/5BVd2mqlcDVwCPEZH9gb1V9Qu59eAMb5tGSB0rhIi8iUy9/JqqXpWXPQ74L+CBqvoDb933APuraqmKEZGLgMtU9eUi8ltkJ+PhwJOAt6nqXqlt99lnH33Na15T+RhCDty8T+NtR2F5931Zv/uWqex7v/s/dKTsfaO0fdR9j8o0z/uohG3f7/4PnV5janLHHXew117Jn3Ftml5DdTv9Aw44gC1btgDTu1c0ZdzX+rHH/8El3hPv3PFzDztY3/Omz7ZS16/+5n7fA270it6lqu/y1xGRq8ks9Ar8g6q+S0RuUdV9vXVuVtV7isjbgC+q6pl5+T8BnyBzSZyqqk/Kyw8F/lxVn9G07ZXzJIjIm8kEwuFOIORcn//fDPzAK78vw9aFUlT1wyJyEvCXZEGOhdx2222cdNJJdXcDZIGLz35ls0CpUdnzoGO48/JzJ75fgGf/zgXAYY2Pe5S2P/t3LuCMk/+m0bZtMM3zPip+22fNGnPRRRdx2GGHjVRHG7/Tk2paDk877bSB+8ssxTDN8rXeCRSkuSs/5MYKgurxqnqtiNyXLPD/WwXrSqRMC8obU0kkiMhbyQTCYaoaNvxqMqHwZLIREIjIbsChwPCMNNX4MzL/y01lK47CiWsXcsLrx7mHbnLGyU+fSiczax1bl4l1mHUmPpq176IrrhULdjbGhapem///kYicCzwGuEFE9lfV63JXwo/y1bcAD/A2PwC4Ni8/IFLemCp5Et4O/C5wHHCziGzOX3vlB6TAW4BXiciz8hwK7yULRPxAvNZiVPVzwL8AL2+yfVNm7cZZh/DYpnHTnYfEP12kSUBo17+LSbSvSWdvAmHB6PXaeZUgInuKyD3ce+ApZCMBPwacmK92InBe/v5jwPNEZFcReRBZgOKXVfU64HYReWw+8vAEb5tGVLEkvDT//69B+V8Cp+Tv30A23OLtZEMdvwQ8xcuR0IRXAV8F7h6hjlJEBFUtnCCnyzfTOjQdg24MkjqPMSFW5+m+jTbUpa5VadITSbV93vz0zE2Clt02JhYWgwkOgbwfcG6ew2cX4AOq+i8i8hXgQyLyQuD7wG8BqOrlIvIh4BvATuBlqrqe1/US+kMgP5G/GlMlT0LMxxGuo2SC4ZQmjVDVwyJllwHLw2tPjnnoUE0YjIdUJ+l3uimBsMwS6w1DnsfROVcVCrGhhrNsffMFg83tYEyTPM7vFyPlPwF+I7HN64DXRcovBh45vEUzbIIn+tYEn3nqWGf5Rt51UmmEfUJR0FQghPto83t19aUsJKnfwzjaMW6cFcAXBXUEglkRFgSd/XkX2sBEQsAsiYMi10iVG7dZGUaj6rkbVRT431NRZz0qTetuw6ow7euwjiUhdDmMYhkyOo6JBBMJPtO+URVR9QZcZ7rdaRzvuHz0k2aS7S/al3uSn/XzOg3B6lsERnE1rNMzoWDMLTZVdE7Xb66jtm/ax2cujzTh6IQwpiFFmXVh0rgRCU1GJnTlGOrgCwsTCPOHAKK9Vl6zjFkSZoiiG2lRIF0X6FqHNgpVj6HKzIQxcVB3P3XXnQSzJgqbBi5a3oQ5ZsITPHUVEwnMRyBS1zqJRSUUBlWtAak6Zhk/ViFLo3zYQPm8YELBmGfM3WAYLRGbmKip6T2c1KhsX10jNTNjlwXCKB29DZ+cR3RiyZS6jIkEw+gg40zENC3m5ThCzIowpyhIr9fKa5ZZeHeD/cCNcTHOYNN57XAnzSijEsJ7xzy7HcxSsrgsvEgwjDbouvnfGGY5N6T6QqFuAKMTBvOWsXGejmUUWpwFcmYxkWAYY8Ce9LvNckVP64Gb90kui80BMevWBBMHPjrzroI2WGiRMMs/ZqMbdGF2TaM6dd0L11x/ayUrwazfS0wcGCkWWiTMuuo3pk9szoMq+RGM6VBHIJy+eiR7bt6HO39Svq7fydo9ZU5QZn5kQhss9OgG+zEbbWBCYPYIxYLvfohNAOXiDozFQVDLuMiCWxIMYxyYaFhcTEjMGWZJWFyRYD9mw1gsYlN2pwIY25r8qevM87EZ7bCwImEcuBtOFb+nMGjitAliDKMedX4//nBHfxtfKKSCGm1444KiivTWp92KqbOQImFUK0LsZlJ1SFVqPZtq1jCaU/T7KfrNVa0nFqcQwy3vmqXShEEzZj2eoA0WSiS0+cP1bybhzSa80dQVECYWDKOYsg4+vnyZdcqfDIt+r7NiVThx7UJOO+2JnDQDbTW6zUKNbhj1x73M0sANJPwcrjvKfgzDSFNkNUgJBPffvQ+XVaXuw8Ykf88nrl04EyJmZrAJnhavN2r6I2ryQy8SEdPGtS32Mox5xxcLRdYFGdhmeJhkFYqsG21i4qBlVEF77bxmGOsRKtD2j3uZZWTg9lO+z9jnpu0q2s6/oZlg6D6xgLx5p+616cTAUsLCkLIkrNNDg88+dYTCuGOOTCAY42KhYhJ8wmyLk7rJ1jFtlgkFv7zKDajKMZbFV0yCsA0WoxGnLCp/HmkiDqAvEJZYokdv47+/bmhRqLKvWIyCK3P3l7a/JxMEk8NGNyywSPBp5kpId/brrEdNmXV9n9XasbQxjCsWLNmGqTMcNlZsA2led4yyUSThMS9KZwmD4/wXNeg13rnHf2fhuk449BLnrOq5DEc++P/bDJY2cTBpdOZdBW2w0CKhbsdZtZP313M3pnEIhP4+0haH9l0l1TullICp07aydRbJzB4jTAg0r0KpeIjjsFBwv7mljet1eeB/SJVRD2WEVoVQIDT9XkwcGNNkYUXCmavPKFzu30xG7eSLnmxc/W0wjQ6zaic+7c7cdTJtWkG6QGooYIpZFhBNXA3OtRD7DRcNiZzGeTIx0D1Ezd2wsI9hx6+dHy1vY4hUEUXBU+O0NkybSd10y0ZsaMG2s0iT8zprQanFomd9QwCEv6Gl/C9cP3zflkiHwYmgXKfvMjsWYQKhg6hCb72d1wyz0JaEMmtBGNzUFL+esvrCNo2TqqbXqkloivc1nU4p5ebo+hN17HyFqbyb1ls0Z0EbtHl+08G6aUEdE+IwfF0X/RZHPYaqSZdMHBhdZ2FFAgwKg9RNJxYNXVU8+Dcrv55weaquNjrnVL1ly91+/fNTtS1F68b2XecYw7aNW0iNk2kIJ9fxjVMkNbVuVNm2ibUt5taLjWxo+3ryLQqzOGFULOhyVtreGha4uLgi4fi180vjElyn7g+fcoRPK/5NJ7QcxLb3KRId4U0xNXIixSiBkykfriBRi0dTN03sGP39heuFgaEx+nUMpuoVhkd8FFkbqqxTl2ma+tu2ooQjX1Lny7eCFKdOLm6ffy2HwjsU4bFt0scxeK21dY7CEQ+zzDwcQz0ULCZhcUUCZELhrNWjkstTnXrVdeu4KqpaJ8o6yFSAVl3Kgi2blNXZt/M1j1pXtn15Yqqy7Zqa6WfJ9x8SCqyy4686qqZKsGuVc52yzGV1hJ3+dG72szDXQ6qNXW+3MRkWWiSMiyZxDG3EPkC1DtW/mdZ1oTStO6RObEZXKBtJ0GRSr64xSuc+qXb4pKwHWR3N3QfzlvioaBbLrrSxc7jAxQVnNu9kLXLc2nkct3beQFkYFb1c8NcWvkuiyl/VOsu2Dcua1B2rJ1VWZVnVfXaNWZ7/oqzdo3SYVdKQ16krlhApXCfMj+DwhWuRuFjOaxmFrnS+fmxB16aw7jJCNgSyjdcsM3t3s5bx3Q2hWIDRnmjrdmh1121DSMT2X9TZx+otutlW3Wddyo6xqpCr41Zpi5g7yG9ryn3jYkFGFa1FQ0SL1q1L2KZRhXVZkLFbFl4XqfWdYPBdW+5zG3ShQ+5CG2aXPCahjdcMY+4GMqFw3Np5hfEJKWImzVGfdItMqHWpE+ToE7oN6gqP1A091o4yd0TRfqu4SZqM5gjft9FxVAm4DM9b3f12YcRH0fkOY2bqjJgJiVkLlliKWg/88xJeL7Fz5j67gMumlpRxxySkRlAA7Ll5H+78yeD6XbFuGLPDwosEJw6cUAA4e/UYoN7TZJs357YFQmpZ2DGm2l9mKUjNT+FHnceWF52vUSwhMZHl2hCOzIi1K0ZZx+ev12Ya7i7GZhThXxMpgVg3ILVIIMQIR8c4IenqiYnKomu5Czk1yiwCVS0GsxBI2RUU0Bm3ArTBwosEGHYzPHftXM5ePabyzd5/+pv2U5xP2cQ3YYdd5ak/tm2ROTeMd4jNwDcqMfdI2Tp16i5qa+wcpc5HWFfb58Hf/7Suw9Bs79qiea5Lv21N2xjGC8XKY8t9oeDjX5NxATmZxEqpbY1pYUMgwURCkueunQtkVoXwxpF6SoXBJ8mihELuibbpjbKuSyImAIqe+vxtwja6ILSYpSDVQcdGOoy7MwuPL+XqGHUUxihDNasIhSZioguCNXSd1BXTZZYwn9h3FyZPill8wkyo4TptBVyWCYVQDMSmsu+CRcNYPBY+cLEKZdHUob/evxmGAWlFQWt1Ohnne22yvf8kVzXRTNk+3DkoWg7tDfWs0s7QZ+3+F43MiNVRFhBa9dynjr1OAGqV9kxbHMQIO9vUtVvlOvO/F0cv/wstB+4vxBcIRWiLs32EQsDN8xCzFoQCoU3MOlERnVzgoog8QET+TUS+KSKXi8j/ystPEZEfisil+etIb5tVEblCRL4tIk/1yh8tIl/Pl/2tiIykdM2SUIJzPUC5ybroia+s84S0OdRRJ3NckW98XH7uUeqtYoGps6+mpvwq81Y0rXsU90JbQaMxF1PqOqkas+LWGfeU6H67fIFXFGMAsIMdG+/LLHCxY23zCT7WQVdJSd12BkijIpNzN+wEXqmq/y0i9wAuEZFP58verKqn+SuLyCOA5wEHAT8FfEZEHq5ZEMU7gRcBXwQuBI4APtG0YWZJqMBz186t9OQx6lNy2TDDKjdg30pQ5ebeBrEhZ2WiKHxyd+/9spjVxf+fqr+twE8YNj837QSrxE00oSjrYEjd8xkeb9F6bQiEsu/Vx+/4/WXhfCzLLLPCykZZ2XlvI8tn20NLzc0w/6jqdar63/n724FvAvcv2OQo4IOquk1VrwauAB4jIvsDe6vqF1RVgTOAo0dpm1kSKuJiFD66+hygHZ9vUYBdmW8/RfikWLc9RRkYi/z3MQtGTKT45y1VX2o0RFXCbfz6RvExl7mdwn22GUCZouh6GQdFQiFF2Ma6YjomLP3//nXkX1tLLNMLfgOp79BtP7rQCX8XdSxA8eydxrTQNkc33FtELvY+v0tV3xVbUUQOBB4FfAl4PPByETkBuJjM2nAzmYD4orfZlrxsR/4+LG+MWRJq8qy1D228D5+0ynzKRX5kn7ZMtnVjFWICIea/T7XX/x8mqHHETPd1Oo/wfId1VWEUi0BVJtVpp/Zd5W9S+MGBRXEZVcpXWGGJpY3/Tgzswgq7eBaDHj12siPpigi//9Aa4f6qRCS0lWlznd7GaxyYq6EB7cUk3Kiqh3ivlEDYC/gI8ApVvY3MdfAQ4GDgOuCNbtVYawvKG2MioQHPWvtQ1AwdM51DfRN42JGGowfGcZP3663zpFfVVB0jZqkI21LUXv+/X58vUFxU+zSD+doWJE1ETnh9Fo1ISZ37Nq+9sI7YtRDbT4+eJxD835jrpgcFrW99yLaol5tBNj4PZ6lsKgx8ETCJqbuN2UBEVsgEwvtV9aMAqnqDqq6rag/4R+Ax+epbgAd4mx8AXJuXHxApb4yJhIb4FoXwCS3mWw//whS7RfhP2uGQLZ8qnWGVzrLO070LItvEbtGbetn+ip4ui/adOg/+Z3/YW+wcp8pivviYCCxbp2xfdYmNEAj/p/5SQiAmzHwzfuxaaGOESsxSVVWADMYjLLMLmzY+78ImTyz0hQMwdB5S35fvElNiLoT4Xx18i4EJhK4y0dENAvwT8E1VfZNXvr+32jHAZfn7jwHPE5FdReRBwMOAL6vqdcDtIvLYvM4TgOH5BmpgMQkt4nygZTe6mFnf3WRS/vvY9o5UR1wnmr2uD93fvkePu7lzqDzWQZf5hOs8pcY6q9ASUnbzLjI7u7pi1pVUlkm/bNSkQbH6Yx2a+x+2s6k7pix1dux7i9VRVfxWLR9cZ3ngf789oTtrmZ25qEglesq2G86TsE7PsyQ0F3n+tNfTEgTmamiAKqo7J7W3xwPPB74uIpfmZScDx4nIwWR69RrgD7Om6eUi8iHgG2QjI16m/QCKlwDvBXYnG9XQeGQDmEholfAG5ROaRncGCYlCYlaB1M3Zv9mFHYpfHnOFxOr264913v4TbdENvazzTO3fFyyxjq+KhSGsO1Weaq8fAFdUR+wYRxEFqe8ypGqHX6Vzi81l4GhyLJMYCrmTHezCCj3W2cl2T6D3xYP/e1xiaaBdSywNWCRiuI7dt0AUiaKyeqqUj0NEmEDoPqr6n8TjCZJfnqq+DnhdpPxi4JFttc1EQssM35yGI6z9ZUBueIx3QGHWuDDAKtbppgIG/Rt3WSdTJzah7YyBsXiD2OcqdTh3ht+GutYGv05/sqBxdYSjiIPQnF60j7IMlH4dqSG1sc6zrfMS1uML5+XEb2cnOwLB2bfshcfsPhddm04Q18m8WYVZnE588VB0irFMXcFEQks4MeCeZmLCwBcF7vNOthfWGxMK4bKw0xpnJ7bM8kaobMo9EXbKfpvLbqRh2/3EOWVUXc+1r2onGpryy5Je1bEutPkdpb73MlFU5bxVaWddMVhG7FjcZ2cFiMUxFCVMSl2Lse/HtzpUic2pIhQmLQ7MijAiNneDiYRRcDeqFVYGynsDnXr6InPLYv5Q/3PRDSoUBu59lZt67MnQrysViCgowmAMRugicIRuhKIbacwdUvb0WtUiEros/A4gJDZm3g+CTHVKVTvIqiIudpz++zoCI/yOyoSOo8jC5X+XZW6ZNggtYcP77o85iIn18HsKfzOjJCBLufnqMu3YBcMIMZvXCBy7dk4tE7h/w1J0yFTaJErap8r2/pNX1VwG4Q1QkIEbs6u37XbXMa3H6vOFgH/cKYEQBsH5+Amilsmy+Lm/uh1j3aC+Ouc2Vm/dTruKWyoV39KmQKhikVlnnV3YNDDCAeIBjX4gqftLWVia/A7bClT1GcXycM31t7bWjsUkS6bUxmuWMZEwIi5ngv8UA4NPNa7M3ciWNjrawU7MMaoJug2fvmuH3xkMm6lTEebDwxp9a8ly8BdrX5kAcO/9/0Xruz8ncAaHyC1vDJ3LOhw/UU/251uLfAtPLPityFJSJOSqRfUXjwiImd/DdYo6crdOmMo49jcKVbav0tnGRVE/LVH/c7WYidSxjfLbSZFKnhTLo1C8L7uNjwNFUXa28ppl7OpqgV3ZLSoKfHzTp29RyDqkYbHgdyh1btChQIgJhqYiZNjMPiyOYh1Q6E4JXRn+MYXuCb/NoXk4dSzVckEMfl/hdxP/Dp20WB5wP4TtjlkqioRBlc67+FgGXQmx5f33ywMiKKRqLoSizImp7zbWnjLWg78YO9m+cS1mbfPFQV8g+N9L+B2V/a6quLZi7S4iTMjUNNtibBuLRTDawmISWuDpa2dsvP/k6u8PdTCxGAXnbkh1Rv31q6WxLXs6LArWK/Ohhjf99Y33foR53F/vCKPKY0MrBzvr+LwRK6wM+JZjpv6iwLyeZ8Hplw2Lt5BYzIW/byd86nYmsTam9luV4sC9+DUXo+yaiPn2/eWp7UahyJ0SjmxwORJiQb2pNpX93kJBVzz0sYq7xGIPusysuwrawERCyzx17d0A/NvqH22UrbM993/GgqniF2GVm1AdwuDKWPBcXQtDaM4t6xiKhmH6wWb9kSLDAXupHAphvYPr+QGIgyNR/H2Hx+X267+PPRWnEhhVGRqaokwchEGVdUgJsbrEBNo4SV2f7qk9FkAcWmpiwad1fmP+91LkVqpSp83R0HVsCCSYu2FsHL72d9zGjSzJcv78PmiWLxp/XcV/H5LqUKoE8IXvU+0Zfrov3yZsS8xlEott6JvF/bb30+0OtmPYFdNfP3QDydC+QvN0Eb3gb4f3pFon8LTqE3WV9aoKE98lVKVzHzWQb3j/47nduADS4n2Xx61k65W79eqeO9/NNnhOl4Jthj/XjTcwgdAuFrholoSxctTaWQB8dvWlG53UDi8vQhi/4JvsU0PUXFnRGHKf4Sffajf7KjfCzBfcS95QQ4uA+7wLm4Y646J2hXEDwEYdsWXh/oePLW3JSVkFUpQJt7r1pdZdZhlF8rkEBv3poYgrcjXE2hUjdA/F64vn7AjraUNg+G2KdcBFxIZCunb557CqBaBqO8OyKoGM4XTRTihUsTqYQDDGgVkSJsAT197BsmzamNIpG+UgQ0+5LsIeijuUqkPohgPW/PrjIxNg2NRd5P+Pbeu3P9vn0sD+UkGCRYGf/nH472MZLsMRJYP19WdNDTuO2IiMso69LECt7afnlBAsC7xrsp/QMlJ2zbUhBnwrTVViwY396zAuKGIxLzA8GqZNwmsl1vH7ZZYzYbr0h6mP/pplzJIwIX799W/hs6svHSgLo7H9cj9JUYj/JOl3ymEHHXao/r7CPPcOP42tqzPWWYYWBP/zM9feX+2keHx29aWlP6YwQ2UqpmCwncWWhPS+hkVSirqdSd0n1rIn+bCtZTSJI2ieHKjdjrZsX64jXg6u6ZgoWGGFXVjZCG706ccu9IkFyDY9L1U6fhMHU2ayEzx1FhMJE+SJa+8A4DOrf4ggrETM7o6iKPSUCTvmrw/xzfT+vvp1xAVDar9+B91EHDjcuXFiIWa2TrXfb7f779/A/QmAsqeD2HE0G5lQpROsWmdq/36nlw4OTQulWHvCc1D0uSkx90DdtqXOSayDdsmRlvCzMfaHqy6zPBBDktUftwT59U8yOBOG3Q6GMU3sSpwCT1r7BwTJ3Q+bBsyboeshTNDkWI907GUdf8qsX+bXL3KBuG1c/MWoPHHtHdzJbQP7DQMDw3b7bfVjNsJ2lg1xHAX/3BTlcIi1LVZHtX2Wu45ipERSGPjpEzPBV21v3dEDVcp2sGPANeEHkfqft7F1o8wXLSussM76Rj2jtLkJZSKgrhXB4hHGg7LeymuWMUvClPj29T/il1hnWTaBbtq4kGKddfjkHHMnhOv52R0dsaBH98TkAgpDYh3rFdffxEkjWA3KuJ3b2Vfux7puH2jjStC+mFBy68fWiRFzm9ShaFRJESkxUxd/6Gj1bdqxFJTVMaqrIWVRCIMqQ/dCaMEoEgGpkTup9o8j5fKobgUTCONC6c34yIQ2MEvCFHni2jvo6Tq7yCZ2lb1Ykd0R+kMmHcPD+NyNPj5CwK23wiY2sXv+5LRpoHxpYz+bBp7OY/U5nrn2/tYsBkUcv3Y+u7Bp41xAGx1OuRWhasdZJ1CwbARK0XFVDaCrIxD62zT76dftJMNgvVSAaN1gxbA9Za4B35rgWxv8NlbLkjjadRimYq4jEMwFYUwDsyRMmcPX/g6AL//Fa9nZu4uerNPTuHpzQXuhiT0MQlxieSOpUwyX6GkXyWMidBNLrA8Mz3S3slHiDEahxzpLssy6xuMnnPVjYH3PKpIKCnWUPW0WjSBJBXJWJSYUqgw5dPSHlpYfZ3E99Tsdv60xi0RokQmPLWXFqJLroUgA+O3zvx/fzeDEQRH+8iyLZnVSQxjd+9T6/v8q+OuaFWF8qCVTAkwkdIbH/NX/4XMnv2LDvOWGSy7lnUHsCcd1EH5a6Co4YeLzHye/EjTrbJ609g8NjqBd7tJb2cTu0WWxeIPY6I2QcFx/WGed4X5VREHVp+4qAiELxnP7Hj7+SZFKaRzL1eALhdg2bv3hXB71XD9F7qXYvv0RO7ERImViJBQ+RWKgasdfdbSDLyqOXzu/Ut1Gc2Z9+GIbmEjoEL/++rcA/Sf9JVmmp8DGk1v/hve0tfe0uu9DX//GVusblaeuvZt/W/0jVmR3dup2NBBLwzf/dOKkwYDPQWL1pZ6Q/Sj5OvjxB1WTYKXaEV+v+siGJuv7+DEAZSLIb7/LMBpbpy6p0QipdatYJspGtgyKgliypOZxBXWsCE3cFIYxCubk6iCHr/3dwNO+uymtsImnrb2ndYHQVdw52EU2sSybBuIsdpW92MTu7CZ7sZLHVbjz40gnber/xTqXVKfir596ovaJiYG6AiMW+zDYuVXr7H0ffDh3Q91YgKKETrH9hdvFzkOsrI6ACK0xTqDsiORASNEsJiIW8BhGEfVfqTqqxBv4dZibYRJYMiUwS0KncZaFRebQ17+R/zj5lfR0fSOXQhlOKFywesKQJcEnZQ73iZnCXQrfsK7Yk3/dHAF11ht2saTTMqeSYcXcBFWJCaYyV0GRSyfuUusLsx3sKGzPYP6DHuHTf2wkhN9mv90x91N5UGO1cxdbLyU2YuuZQJgMCja6ARMJxgwwiiskzCHgm73LzPnpxEVxX3bzNsbbUbaP4dEt1ZM2+f/rtstnlCGBKbdFKD58l0FsGKTPdnawKTHZUzgSYg/22BAeqWyKWehaNfdCFYtAbHkY7JhaD8yCYEweEwnG3OICOi9YPWHADz/4Pp5tMPXk7a83KkX1Ncn+WIUw8j+2PIX/dF0kDorcNbHg0dRTetHTe9GyMiHhuIu7Ci096TiHonM0vmyJFqg4aWx0A1hMgjGHfGL1d/nk6u9vfI6N/ihyQZSZ64snoIrnUKj6xF1lgq9RGKXeqqM/loI/R5UhjO5zmLNghRXPzdLeuSn6XtICZXy3zVTdJhCmgcUkgIkEY844b/U4drKd7dw9UN7PGyEDnXyTYLXhIYjhSIh09r9seToRU5UbSirxUNkMmqn2NsHFClRL9hSfj6IseVGY8ChFKuI/dctOtSWk2MoQD0Yc/lwtWDU2+6MrN4FgTBMTCcZc4TJCxpJApTJU9j9Xy6IYY3Augb4P3UXa+/uIZRis2nEXrVflqaXpiAafVIBmqkMMExTFYhCKRkwMjszo5zjob9/PHQCwlW2Vjm9YVFTLuhhSNnohW6e6ULBhjt1AGc7G2fRvljGRYMwdqSyRzu0Qn8gq5mIofzKPjZBwnVjZ8MXU/otSFIfzTNTp7Ee1IFSNyagzsVXowqgyvBKyYynreNMBjIMWhVQGzFFv7nXmEBncrn+sFqg4TczdABa4aCwQF6yegJsdPuUyGOysi03h/pOuj9+J+370skDBqp1wWE/4ZN2GO8EnfJoPy0ehKLlUakhklX1uz0ctbGKF7ezYsDYUtSEc/TD6PA2DE6n5dYeELgYTB9NHsYyLYJYEY4HIrAJauE7dJ/M66xcNqQw/lw3NDIP4Yi6MphRZP8rESKpj9wnny4gPMSy/OaeSYYWUD0usN6X3KJS11wSC0TXMkmAsDJlIiKcH9il6Mq8SsxCLNUh1/FWFQ2qdVPKfUSiOeygWIU2HLVZN9Rw7xv5T+nDbticSMA0+uYejK4pzFYwTEwhdQs2SgIkEY4EIp9/OsjHG8yTUyVMQdtBLLLGDHUNxCm5ZjCYBk7GcA6k2pWjDZeDvs856VY7ZN9OX5VVIsYmVPMIgnjApq3/0czCqSLNRDN0iczdY8KiJBGNh6HdOgxMchQmUYuZ/SLsXYimaV1gZ6QZTNglTlSf6kJhwqJq7IFVnEYMipp85MeV+KHI7pOIWYoGjsVELw/XGh0PWpYkwCEWOiQOjy5hIMBaCT6z+7kanO5gnYfCptsi1UOYySLkWmlkJ1guFQlULQNjJtumWiO0nFaewlMeChIIsJhCGhzimZ8n0z0N/boZ+kOKyt5/+qIXqMy4WRbA0OY8mDmaJeEruRcNEgmFQP0WxIzUiIbQ+NBUKZftOiZaw44xlNXTL/bKUSb/KFNHuOKskQHLWlnC9OucqNQzSCYWwzMctH1cuAv88xASoCYTuY+6GDBvdYMw9n1z9/SHrQVH+g7BDGTW4sGleg2z7etkMq6Yujj3B+2W+iIgND0xRdnxu+cpGnEDxuQ5jPdx/P01zOMwwa29fJLlyJ0rqzNaYCnOtGu8RDh+1DIrGrDEVkSAiLxWRq0Vkq4hcIiKHesvuJyLni8h1IvI+Edl9Gm005od+SmZwEcvbuZud7Nj4C7Oj+e/rjEoIGWVYYpUUy2FbYhkVY0GNfuccZoeLJTVqw0Xh55bw61wu+Au399vs8IdC9use/n7CYFK3Xko0pKwMoXWmqqhcZ91GL8wUOvTbaPo3y0zc3SAizwXeCrwU+M/8/ydE5BGq+n3gtcBXgVXglcCfAK+fdDuN+cWlTW4y3K4uTdM8O4HQZAhWGwmOYu+rELo6/Dpi56LqvA9FGSjddxkKhbL2lyU4KiI12VXq+z5u7bza+zCmi6LmbmA6loQ/Bd6rqv+oqt9U1T8CrgNeki+/J3CZqn4d+A6w7xTaaMwRWabF7fRYRxOhaLEn/rInxKY3kCrbhelcY1aFUeZhqDI6oq5IKnJ1SGC4r/qUVccSUyUoM2ZBKSJ2taTiO2L06JlAMEoRkSNE5NsicoWIvGra7fGZqEgQkU3Ao4FPBYs+BTwuf/864DQR2QH8NvDmybXQ6DrnrB5be5swh3rV+RMg7XbwO+gq2Q7DYL5wu3Jf/mBw3ygCIWxTSKqjTbkhUvMtuHO3G7slxRmkO+v4fBqDs0+GLoui40rvv3j92HDJWCxH6EYBsyDMOpNwN4jIMvB24GnAI4DjROQREzi8Skza3XBvsiHTNwTlNwBPAlDVS0XkgcB9gBtUtTiPrrEwnLd6HDsSGfTKKIuabyOpUMzSEAvEczeNcQxF9Cl6si46F0XbxW54sTLf5bCDHZWeRsL9hqMlXNkygzNrrrCycV2EgaIpXB6LbJ/peR0cdbIvuraYQJhtFC29LlriMcAVqnoVgIh8EDgK+MYkdl7GtIZAhh2/+GWqug5cP9EWGZ3HTQPdlLDT8DvsMoFQ1uHGYgjKhIers45PvKw+CDMJNg+aim0bC4IM9+k638FcEYOkciHEhkXGzmMqiNEXJ8P5E4an7XbZMcvwp6KOEYouEwhGDe4P/MD7vAX4lSm1ZYhJi4QbgXVgc1B+X4atC5XYe++9ec1rXjNquybOAQccwGmnnTbtZjRi1tq+vvleG+91973hoCdvfJZcm1Z5pnfd0rB/PXsJUmhWdyyRPaVIvnfJ3/fy/0l23xs96MkD+3DrZ/X121k+Q0UxsTo0UuvSxrLBc7kUnAfZfR9WDnp6dF/uGGSjhsFzhFd3v22Dx+1/f8sbdfS3UXSjI++f+2z5LgXf2vLu+7LnQcfkbYi3O+Sq62/htNMOjS6bJOP+nZ500kljq7srtBi4eG8Rudj7/C5VfVf+PnYhdcaCPlGRoKrbReQS4MnAh71FTwY+0qTO2267bSYv1tNOO20m2w2z1faPrj4H+Un/KVcOeipy+Wc2lrtfZ8qiEHvaTUXShxTHKKxv3AVcXWUdux70ZOTyT0cDGpuOgggzIKaoGhCYGjWyctDT2XH5BYF1YWnj6X6ZZXbJrQgAO9nBTnZs5FMIcecuy5eQtW0bW4faEH5nMavDem73iLHnQcdw5+Xn5nUNXwuha2OddU7qSB6EWfqddhGl1WRbN6rqIYllW4AHeJ8PAK5ta8ejMg13w5uA94nIl4H/Al4M/BTw91NoizHHnLN6bGW/f9mcBrHhck0FQhj0ViUqv4gqIwRSs1q2OXIg5maAdACiI3QDuOVF59e5JXZhEy4oNYz3cHkRXD2x0Sp9l0RxXIKf3jl2DM5tYYmSjAZ8BXiYiDwI+CHwPLKg/U4wcZGgqmeLyH7AXwD7A5cBR6rq9ybdFmN+OWf12I0I+xC/w4tFqRcR+t5j8QhVBYL/uWowY9Xhd0XCAIZHcoTzTsQ66HBfYZChnw8hVkcoTooDSbOoghQ9emxiN3qsszNPhuW3KQuW7Hfso5iNywSCwwTCfDGpwEVV3SkiLwc+SeYt+3+qevnYd1yRqQQuquo7gHdMY9+GEeuIizrnsuyDZWmeHaG48OtrYk1IBT1WGR6ZWhZum+rIwzkfQpwLISViwpwK4WRWPXpsZeuGxSBGlvti0HrhrAe+JWErW6OiqUqiK6FcHAA8d+3c0rqM2WNdJpNMSVUvBDqZjtPmbjDmkmPXzoma1cNgs9h8BTGWIn8+LotjrPMN666SedDfxr0Pw+tGGbkwKqEFporlAQbdCK4zH3Y19MvC7I3+e1/QrLPODnZsvPfPXSgkYiKq7Ikx/M79PA0mEOYTZ0lo4zXLmEgwOsuZq88YaftQKAAbHW3YyYWdtN8phH72sqdst437K0tAVJR4JebjT1EnO2Gsrak6fVLWjth6/p8gG9aFUBzswiZgMAjT5UMI95kagplqT9HQyaLtY4QJnNZZZ2seLGkY84pNFW10llF9vH52xlg+/6IOompnWIcyE71bp2ruBr9NTeeI8InFRcRGd1TpqF1MQGzdsGwn270AxMFMl3VcMWEMgj+vg6vTJVGqgotF8Edh+McHFocw78y6FaANzJJgzC3Hrp0DlPvmY0/gRb70uh2yX1f41B7uv64LoQ2BkEqt3EY9vqskNoLEzcK5na3szDteZz3x13euhB1ekKLbp/sfJlhyHXkswNL/HpxdY3CIY/Y+lQvB3AzzTzYEUlt5zTImEoy5xgmFKtQ1RYfBirGn7tTn1P5ToiUldNzT8qhioaq7IXSNlAUmxvYRbuPHFDgxEFvuf3brbmXrxrbuXPjBi4P7rzLz5/DU0aEVoUev1nVlGLOMiQRj7nGxCb6fOuW/L+oUhzuv9bxLKZ5YqO5Tel1rgrNu1HGFVOnoXd1F7Uq5E/wn/Njx+/v1255qV9FICT8w0d82zI8QS6sc1uubl/3hlNvYxl3cxTa2mUBYECxwMcNEgrEQHLt2TmFnmHIBhGIiJhRCX7ojFpxYBb/jqrJNlYDFMOK/jLpxEbHjKxNHfifsCIVclTb766znbooVVqICwN2wU23zcyKEAmkdm/Z50TB3gwUuGgvEc9fO5aNnPj3Z4aQ6eH95GBTnB8eFJv8i83qKMH7B32/ZtkV1FZWl1vGPFwYDA4uCCf2AQ8cSw/ETMSESm/XRP88hYdmmPI1zzPLjsio6IeG3oQwLUDQWFRMJxkJx1fW3DCRJT+E6utDnHwoFR2wIpavHZ5TREXWo67KoW3dRQqqiWRX98+YSHsUEQEwshNsPWhp67MauAGxnx4ZY8AXN3uzD3dy1EfzoH8P2jSDHwREZgrBrXq+xWLQ8d8PMYu4GY+F47tq5pW6EWHnoUgjrcJ1LKoCwikugybI2qJNjITaUM8yM6IIHi+qNWRJSVoCytu7GrrgREcu5mMsmf+oPX9yFTREXgj/KIZ5G+llrH0q2xZhvzN1gIsFYUJxv2Telh8GNRXEEqfKyjjFVR5siIIyv8CmKt4it7yjqXP3tfNEQO6Ywk2KqvlgMRbiev76zDrjkRs5K4e/zdm4ZsiLERqG4QNCilNCGsSiYSDAWFl8oFFkVqgqFWEcaq9/RVsKmKqMTioSLL5RCM394DEUjQ/zRCpB2wYT7LcrTEBMEVYYyuuX+sEpXthu7DdWbGvZ45fU3F+7HmF+0JSvCrFsSLCbBMBL4/u9Y8KAfoxB2gH7gXt3hiSmym1bxbJGpznOUNqRiEMLzEy4LAx5Hpcw/HBsm6crcn29BSI3GmFTciNFtXDKlRccsCcZCkxrSlsqsWPQ0Ho7Tb0qVTjVmei966q/S8aVcL2H94Ta+i6bsWOrEPqTwhyn6Vogqk0ytM5yJMbRkOFeD5UMwDLMkGMaGUDhr9ajK26SSDIXR++HyNiiqd9QOeJTgyjqiKqwrPuTUuSPioskf2bDMUuGoCrePJZbYyrahmITYDI+GsS5mSTCRYBgRwiGOMDgsMkVsUiNX7uqNUZRRMCS13jjM5H6dRSIgTFvsb1NkaXDEjmk4TiD8Pnq4NMruvV++To9NefBhOIrBj0/wxYX7fo9aOyvZVmMxMHdDhrkbDCPnuLXzhqwCIf7wvJQpP2aaL+tsYXh65S4SHnNM/KQEA4xmVamb4taJgTCOwz+/66wPtLet+AnDmBfMkmAYHmE2wZhFwZUXbR/i6pjnTsg/V7FAydhIhSoWlJgwCMvCGRz7Vob+9xlaglLfq4tFOHv1GH5Q2DJj3jFLglkSDGMAF58QWhRiwwIdqY6/KIBxKfIX1l1mTWgjCLAOqZwKdQRQaCWpYlkIZ2YssiaEcQy+4HPJlVw7VlhhN3bbmOchdCXtxT1K22bMLzYEMsNEgmEEuDz9VU3jbXTUfs6AKvsO8zpMY9heXTO9Lw7quFXCjj/8P7iPpYFtiujRGwp29LMr7iwJhDTmmywmoZ3XLGMiwTAShB1YlSGBbY1kqBPIGFLWYaeGd4ZWjTLKEkS15VoJXQlhmcMJAz9wMdVed5x+jEmIBS8ahokEwyilyPVQhdBKMCpF+6/SMaeOoemx+dtVdcX41D0vfgCje+8LglAgFAWS+pYby4vQ5+zVY6bdhE5g7gYLXDSMIc5cfcbG+zouh9gQySrD/9p0FbjgwVEyLJYNVSzq1H0LSCros0o96f33ou/DMn9aaDcjZJi/ojcgLLKyj6/+DgDPXHt/7bbNE89dO3faTZg6CqzPdv/eCmZJMIyA49fO34hLiPnNU3MT+J2P+wvN92GGwtTwydSTb9XMiVUFQlUXQ5jWuGi9Kvhioo5bJZVYabgdvlgYTsecmizqDu7gLu6yJ2nDyDGRYBgJRk2v7BOOAEjlYCjbflZIiZSYCErPQ5EaZlrvXPijGtxkTylSabqNxcTcDeZuMIwkx6+dP+B6gOLhjg6/04tlJCzKIeATzqVQRGqGxTKKOvO67oAwT0KsTp+y+uuKgdi2WRBjb+Ac+u3wJ+k6Z/VYAItNMID+6IZFZ7YeTwyjQ6RcAzFTdtXJlmJuhaaBkk1JJT2qu++2XBPZunWnz+7HI4QCoWwGzfNWj6u1L8OYZ0wkGEYBx6+dP9CZ1e2AY7M0Vs2BEEu2VDQ7ZRv5EsJONCWEytrvAjldmzVhcg1TJqcYxargn5OV/C8W3OmCTneww2ISDMuTkGPuBsOoQBXze7UUw+VCoUoa6LozLtYlPNaw3VUsC77YEWSgbn8Gx7YJBUWRuybmHlliyVwOBjD7HXwbmCXBMEoIMzCOK7uh79Ov6p7wt61DSqSkOv+UaKjLumdhyOodzJBYJRNjFWGxks/+6AhTX4fn1g9o9LMuGsaiYyLBMBoQdsptJEoqm4La7Tc1h0IdYu2t2/E7N0TZEMpBk35qBkx/bobq1oqwDr+eFe+vfNv+vszVYIC5GxzmbjCMFqg6GqDI7F2UfCgmDNoeFtnEnbLM8kZbioZ3VqVu7gR/xsdsux57sgcwKLr8zIqp4adund3YDYBzVo/lmmaHYcwJlkzJLAmGUQnncnCkhjE2DR4sSrA0boqSC8VIuQRSCaSAaOBiVZdCWb4ENyX0MssboxnczI6xlNFlE1O59SwuwTDMkmAYjQij4mMZEmPpmUclFCejPrWXuR18039Y7kYDhLkg/GNPiY1Ux++sAsPzL8RTMDuBsJT/h8zN4I+uSOFbb/x02haTYIDlSXCYJcEwGhLrCGNpjmMJlZoSCoS6ddaxGMTEQvg+JVLKOumiIY3hhE1l7MVeGyLGDWEsy6wYfj9+Wy0mwQCLSXCYSDCMioQuhxhFcyy0NQKhrM4200k76iZTChMxCTIwf0J6FMXSwMsvS7GDHSyxxG7slrSMpM5HLFNlW7N1GsY8YO4Gw6hAmJ4Zyp/kR3UxpPIu1I3+Dyc48uuoSlEnu4MdQ+XhBE4uJsFPWFREOCHTHuyx8cTvP/kvs8wOdgwEUPqWmyKLQpHr5uzVY2wmRGPC0UHdxESCYYxAKBTK5lmoO4VzakRBLBV0nTraICVgfPESm+2xSlvWWd8IPtyVXVlhZUNYhE/7fkyCczn4y0LxAGmBEIopY3HJpoqW0vXmHXM3GEaLtD0ywXW6yyyj6FA8QZEpfRxuhyptdTQRJ+sMx1yssMImdt8oD0VZLBuk79IYd3bKtlFVVG3sXRfoQkyCiPyNiHxLRL4mIueKyL55+YEicreIXJq//t7b5tEi8nURuUJE/lZEJC/fVUTOzsu/JCIHlu3fRIJhVCAVj1A3cLBK55R68vZTG9eljSfj+AiDYetA3VwHRftxQxh3sn0jMVJovUnFIWxl64Z7ospohzba3gYiQn5PNwyATwOPVNVfAL4DrHrLrlTVg/PXi73ydwIvAh6Wv47Iy18I3KyqDwXeDPx12c5NJBjGCLSV1yBlFfBjCFKTJFWte5KEgYuDy3rJV7Z+/7a0Nf/bxlZ69NglmNXR7ct9DzvYMWB9qdI+6Is9Z/3o0eOs1aOanwBj5lGgp9LKa6R2qH5KVXfmH78IHFC0vojsD+ytql/QzCR1BnB0vvgo4PT8/TnAb0iJIjWRYBg1SHW2o4iFMt/+JP3jdYcf+ttl/9OBlo6i1MzLXmft/9/BDu7iLraxdcCNEIoBF4PgB0YuBXWG7x1V0jcbi0WL7oZ7i8jF3utFDZv0e8AnvM8PEpH/EZHPicihedn9gS3eOlvyMrfsBwC58LgV2K9ohyYSDKODhNM1j1pXlYmT+uv78yiUz7vQrE3DwxrXA+uAP6GWsxBsYxvQ7/hdufvv/vxcCVXFW0xYPHjzvnNlUYiN0jEmwo2qeoj3epe/UEQ+IyKXRV5Heeu8GtgJvD8vug54oKo+CvhT4AMisjdE/ZLODFm0LIqJBMOoiD8bZFX/dpXgwljnHZrrqw93TMcN+JTVFzP9F1FHQITuhXSdyxtP9+4c+IGKfk6KshiFbL/luRLCxFiOeelcq+T6MDK0JVdDFXeDqj5JVR8ZeZ0HICInAs8Afid3IaCq21T1J/n7S4ArgYeTWQ58l8QBwLX5+y3AA/I6dwH2AW4qapuJBMOoQXiTDYVAkSiILS/zm5eJjJBUZ+3XUzTqwJ8wqWpnXpeybIvu/3Z2sJWtA8t3Y7eNoZE+4bTPg3UOjgjx//vbFomxFVbmRigY1enI6IYjgD8HflNV7/LK7yMiy/n7B5MFKF6lqtcBt4vIY/N4gxOA8/LNPgacmL8/FvisEx0pTCQYRk18i0KMaib94XUmFXtQFtjnCwV/dsUysVC0XL116tSzno9KcOf6Lu5KbTaUJjtljYkNIXUuCoA98lkkw3TaFrNgTIm3AfcAPh0MdXwC8DUR+SpZEOKLVdVZBV4CvBu4gszC4OIY/gnYT0SuIHNRvKps55ZMyTBapupT/6RGHBSJj5RVYXnAfD8oFIoCD4vcFOXuhaXoeoNxEcsDnXcqY2LR0Mzh/fbPj4ttcPX79R63dt7wxsbc4kY3TJt8uGKs/CPARxLLLgYeGSnfCvxWnf2bSDCMKTBugeA6f///6HWWuwnqLiuqNy0a1gesAkWpnt3xVznf4SgIfxsTCItJF0TCtDF3g2HMGXVHRjQVEH78QFZP8URMdeoMLRluSugwHsG5AHpeW8JREuX7HIxPcKKiR28jy2WTeIQzV59hcQzGzGOWBMPoKOkJnnoDT9nDQwkHcyz4AYvDdVV1jQx23qFAKNqmKqn13TH6MQHhKIThYMXi0Rmxc5uawbPp3BebuT/buLv2dkY3yKaKNkuCWRIMoyaTejoMOyYdWBYG9w2b5P3/o7obqlgJ6g6b9LcrikNwLHl/LndCmDipqjWjKOlT6J5o6q550to/sC0YnWHMFj1t5zXLmEgwjBnBPdPEggrLCDvFOsMqq9Ikp0L1bfqddCxAMbaPOnW77I3hvkblWWsfGvj80dXn8NHV57RWv89Zq0eZe8NoHXM3GEYNTl89cqhsXFMxV8XPaVD+tD/8VBwL1Ku6z/R+yl0NRW0N63dDFF3ipDABUhioWdaeMP1zLCvjOusIMrJo+Ojqc+jR49i1c3jW2oc4Z/XYkeoLOWv1qKlef/NKV0Y3TBuzJBhGRWICoYhx5j0IO9hYh1sl8+C4BE7VkRB14hbCjtxlXwxHcMRzUAxaFqrONeECF0fhWWsf4ti1czY+++/b4Li18wasIfOURnqqTDDjYpcxkWAYFXACoYl5fFI0GVkQyydQbbvYk/fw03qTNvWf8sNkTvH00r6roKitsbiHVFKrWJrscZjyz149pvU6DaNNzN1gGCX4FoR1emxiZSiwLaSNp/NR0yHH4hBGJebSKHM9xK0c6W3KAiCdy8GfsyEM0PQtJL7oGN5XempuYMDdMI5cCc9dO7e1uvzzcfbqMdlUf8ZIzLoVoA1MJBhGAb5AWA46pC7hz7eQynbYRg6DorqrxkW4bZq1YX1DqEGWIbFocqyytqf2MUs4C0dsSmyjORaTkGEiwTAS1I1BGDdlI6mKche0leQonNPBlRXFGTTdt7+vKrNblk1cVbQft45vkaiS0nmanL16zECchi8OlljiQZv3mUaz5go1kWAiwTBSnLh24ZCrYdTOtg6xtMTNn8D9jrD+cfgdtW85qLPvOvuI1e+3ezs7oiMgYkLB7zy35zkVYi6T2LZtBC62jQtMdO6FUBzE3htGU0wkGEaClCWhbbGQchW0PUWzv68m2xV9boMiC0SYYTIVt1CcJCmWnbJeCucusZRnoYwN38wmwjJGxdwNNrrBMGoxjs4xfEqvMp3yKLSVLnkaDAqH5Vpt8+MmwvpSGSt9uuZ+gsG5Jox2URsCCZhIMIxatDGJUareSW0fe5ouGiY4bVJ5D4ZHb9Rva8qyMhwMOf3zAMMpwd1EVA7/vdi8A0YLmEgwjBao04mMq8OpHr3fi1osQqEwzo6xrtiKT041aH1JxTEU7Tc1PHI4riHbrgvWBNe2Xv7nhEKYidIYHVVp5TXL2BVlGDUJn7zLOtWwQ46ZvV35qG1qy8oxKSGT+hwfoRGmk07HF7jPRfXNGoNBtPEkWG46bScWLPviaJi7wQIXDSNKUdBi+L4sL0HsCXecoyTGHfw4CqmRGqlz6FNl/H9KFBRNrV2VZZY5fu38RtuOSux6LItDEJbYGxsGaYyGWRIMYwSqdG5uvcEUw+MJgJwVYk/9o5yX8PyG+RWcwHBJmJokHNruTUk9ScJhuEXnKBQOu7DJrAkNUczdAGZJMIwh6vidi2/YszuKoEv04w3ST87LLLOdHXnK7EELhTPBZ3XVm6vCBf9N47spug5jVqwd7GAP9mBf2cwSy9zEEnvIPvTURj40ZdZdBW1gIsEwpsS4O555Eh2xTtEv6ydXWh8aJglZB+o+150F0+3nxLULWziSalQVquF3vI1t3KY/5mZuYm8OYqdu5xZuGkcTjQXBRIJheHQher0pdTMhziKpOJB+/Ed/gqb4iIhmk15NQiC0de3tZAfHr53Px878LZT1sUxMtRDMgaugDUwkzAmnrx450Scdo3sUBUrOG75YcGLApSj2J31ys0WOwjh/VzFhMOr350Y2fOf6H3PS2jsa12OA9kwkmEiYE0wgGGF657KJnuZFSDghEIqBFW9K75RlIYZzQSyzjJZOq1WfcVqr3Lk4a/Uofji2vRiLhIkEw5gjqgiFWScUOst5QqFsvoL+8bo4hLrBitk+MkFx9fW3ttHkWtT9zpw1xVlSNrEbu7BiImFE3OiGRWd2xkwZxpiZ5XgEn1Ryp9BvP4v4iahiox7cKAZH0YyQKdw2s+TLX2GFZZbZhRVW2MTT186YdpPmAhsCaSLBMDaYZ5fNLOVQKCIc3hjO4rgjyGUQn/dhuXDZ8Wvnj00gOCHq8ji0NRfIOuussMKu7M4ym/jM6h+OXKdhgLkbDGPhWM+z/XedVKZK31ISG+7oExvN4I+AcNutsz61bIpN8S1Dzt2wzjpPWvu7PIHSw6bbwFlHLXARTCQYxlzjm+enPY9BUaxEanKmYWtBfB0/2HB4v4O5EHxxMEmXwumrR7Zm0fHPpe9i2cl2IHOVnHbaoa3sa3GZfVdBG5hIMIw5JnzanqYFIYyViKWprhpsGXM1pPc7HJcwrXiDcXwfThw554VhtImJBMOYY4ZHAkx3xEPsSTqVRTFcHi4LXQ2htSAWe9CFYMRhN0n8+MpSfvetCb18ZMe6iYSWMXeDiQTDmGu6IAhCd0HK5J5qa9+VEBMYfXEQczn4pvjnrp1b/yDGTCgWyqYdD9mUTwu9xBK7sGksbVxYFDB3g4kEw1gkxikaYtMxV5kiOyYC/I4+5UqIiQJfNLj/Syx1TiCUTTNeRniuVlhhN/ZCa+SCMIpRQLsf3zt25mNclGG0xDwPgxw3qbiCMh98WTxBKALC0QkhrtN063VNILSJS8G8C5tYkmWeaGmYjZYxkWAYASYUmpOKG5DIOsX19F0I7n9RzgPod5hL3m2tCzEIjjavK3cOXYbJnWxnq97RWv1GhiVTMpFgGFFMKDTHz4joPleZAcEXD34sgW8xWNoQH33R4Hzy7r3rOJdm4PbW1P3jxNeu7LpRdhM/bqVNRo4CPWnnNcNU/hWJyKqIfEVEbhORH4vIx0XkkcE6IiKniMi1InK3iFwkIgcF6zxMRD6Xr/MWEVnylj1IRM4UkS0isi1f5wIRedToh2oYxiSpkxMg7Cx3sCNP+jTobqgyo+M66zx37dyZczMUBXX6+Ou4DJM9ep2ymhjzQx2pfRjwDuBxwBOBncBnRORe3jp/BrwS+CPgl4EfAZ8WkXt467wd+DDwNOBngOcBiMgK8GngPsBzgIcDxwJfBvx9GMZEMGvCaFR9Si5aL3QxOOHgCK0Fz107d+YyJ4bJrmAwbXNs/eEyC1hsn3ZcDQvjblDVp6rqe1T1MlX9OvB8sg798ZBZEYBXAKeq6kdU9TLgROAewG97Vd0TuAS4DLgG2DcvPwh4CPAyVf28qn4v//+XqvqvIxyjYRhTwnV+VW+T/Q4zPd4/XDZLAYpFwjMmFoAhseC/79FjiSWetfah9htrQK+l1wjk1vkfisil+etIb9mqiFwhIt8Wkad65Y8Wka/ny/42758RkV1F5Oy8/EsicmDZ/kdx2t0j3/7m/PODgM3Ap9wKqno38O9k1gfHq4ELga1kwuD0vPzHZKfz2SJiQzONTmDWhHaIxSSkhgGmZqsMXQ++FWEWBEJV6sYphJNaGXPJm1X14Px1IYCIPILMEn8QcATwDhFxP5B3Ai8im8DjYflygBcCN6vqQ4E3A39dtuNRRMJbgUuBL+SfN+f/bwjWu8Fbhqp+CrgfcICqPkFV78zLfwj8MfAa4JY8buG1YUyDYUwaEwqjETePF+dPWPL++usMWxd69Dh27ZyWWtod/MBPHzdPRfZ+cKin0TLdD1w8Cvigqm5T1auBK4DHiMj+wN6q+gVVVeAM4GhvG/dgfg7wG87KkEKyOuohIm8iUzC/pqpX5WWPA/4LeKCq/sBb9z3A/qp6RLSy4brvARwO/ArwTODngN9T1ffF1t9nn330Na95Te1jmDYHHHAAW7ZsmXYzGrGIbT9w8z5jaE09lnffl/W7b5l2MxqRaru7O2lQJoGDQtGhMoCrrh+us03avtbrXEepO7fkf4qyxBKKcuX1Nw+tN+7f6UknnXSJqh4yth1MmZV77qv3fGI7k2T9+KPnfw+40St6l6q+q8q2InIK8ALgNuBi4JWqerOIvA34oqqema/3T8AnyNz4p6rqk/LyQ4E/V9VniMhlwBGquiVfdiXwK6p6Iwlqm/VF5M1kAuFwJxByrs//bwZ+4JXfl2HrQhJVvR34GPAxEfkL4JPAa4GoSLjttts46aSTqh9ARzjttNNmst2wuG0/ffXI8pXGyJ4HHcOdl8+mWd1ve1EwXphJcImlAXN6mEjppDEHKbZ9rVe5hvwZHuPLs/OyB3twF3clXS2z/DudQ24sElQi8hk8i7vHq8lcB68l09KvBd4I/B5xHakF5ZQsi1JLJIjIW8kEwmGq+q1g8dVkQuHJwFfy9XcDDgX+d539OFRVReRbwC812d4wjG6S6gTDqZ/DqP0wydIsUUdkxlJcZ+X9497E7pWGhBojMKEcB+6pvwwR+UfAKeMtwAO8xQcA1+blB0TK/W225LF/+wA3Fe2zTp6EtwO/CxwH3Cwim/PXXpB16MBbgFeJyLPyHArvBe4APlCh/oNF5DwROVZEHiEiDxWRF5Ipptl8fDLmCotNaEYoBqpMl7y9JBjPCYizVo9qp5EdIjXJk588ag/2YDt3c9TaWRNv38KgdGV0w/7ex2PIRgZCZnF/Xj5i4UFkAYpfVtXrgNtF5LF5vMEJwHneNifm748FPqslMQd1LAkvzf+HwxH/Ejglf/8GYHeyXAj3BL4EPCV3IZSxBbiKLHDxQDIB833gNODUGu00DKNDlAUuhnM91Inud0Kh64mE6rqqUudgKbey7MImtrO1jaYZ3ecNInIwmWy5BvhDAFW9XEQ+BHyDLG/Ry1TVmd5eQvaQvjtZnMIn8vJ/At4nIleQWRCeV7bzyiJBK2SEyBXJKfRFQ2XywIk/qbudYUySE9cunHpswjyRGv5Yvt3ykCti0ZiFtNMzTwcSIanq8wuWvQ54XaT8YuCRkfKtwG/V2b/lIzAMYyqEFgQoT00czggJ3e8sRxWV7pwsscRu7MY662xnK3dxV0stNJJYyEfHf12G0UEsNqEd/BENYRriVJ4A33rg3rvgvbNXjxl3k2vTptWpR48d7NiY/KrrLhZjPjBLgmEYE6d8mF/x80voajhr9ai5yrro485F1y0mc4cLXFxw7KozjAaYNaE5ZSMbyvCHQa57f12jDSuCGw7qMis6oWCuhgnRgdEN08ZEgmE0xIRCM9KuhHp3U18sLLM8d8MhY4mlXByGuRomhLb0mmFMJBjGCJhQaJdYIGMRs5hUqSk9emxlq03oZEwUEwmGYUyEMDgxtTyVTGh4/UGB0DWf/SgC0p8JMzwumxZ6QnR/gqeJYIGLhmGMlWWWBhLGh4mUXFkYq+AEQ2pKad/d0DWBMCqxoM5wVkxjAvRm3FfQAnbFGYYxdgSJxiL4ZbHlMcvDoMgYzpvQJdpwR/Xosc76RgCjYUwSsyQYhjF2tGb0VpNMjF0c4dCUvoVlfWjmS2NCKMiMj0xoA7MkGIYxdto02sbM8EBn8yTUtSakZn0EG/o4cWx0g4kEwzDGT93QrSrDIV0H2vXpkptMEZ3Chj4ak8ZEgmEYjak6IVMTioSCP8GTyx9wzuqxY2vLpIgdswUrThFLpmQxCYZhNCcceeCPSvDXGbfFdZ31zlsURmGZZcuPMGksLTNglgTDMEYgZklIWQCaWB3KrAnOirDM8lyY4lPnaJ31zsZcGPONiQTDMGrjD02s2vmn8h0U1R9b158i2lkPjl07p1IbpkHTYZB+zIVZEaZET9t5zTAmEgxjRCw1c0bRrI7heinqWhvWWe+0QGiLLuaAWASk185rlrGYBMNYYML4gaL1UvjioE5dVSd5CjvIddY5fu380v10iRPXLqw9K6SfF2HWjteYH0wkGMYC02Sq5pDQ7RDLmhi+j4uBYsEy6zEHTYSCiYMpooDOtqugDUwkGMYC4zrmWOfddHhjOAcDkJy7gQrlsy4OfHzXVEwwmCjoGDPuKmgDEwmGscCEQxXHtQ9l0OLgBzGmLA/z3mE6y4LFtBhdxkSCYSwosWmZfYtC1UDEsn34uDrrxDDMMyYQuo3M+MiENrDRDYbRArN2s/ef6uPDDOPlTanibvA5fu18zlo9qrX9G0ZtVNt7zTBmSTCMBSMlClK0EadQtn1oWThz9Rkj7ccwWsEsCWZJMIy26LI1oSg5UZ062iZs0zjngjAMoz5mSTCMlqg7vG0WCAMbm3TilgjImEUEi0kAsyQYRmt02ZLgAhHr4lsgRnnKlzwqwc2zUL7f+ZiLwZhhXJ4Ei0kwDMOYDH4WQcMwuo+JBMNYECbt719meUAUuM9heWw7w+gE5m4wd4NhtEnXXQ5Vytok7PCrCIB11s3VYHSAlmaAnHGhYSLBMFpm1oTCOPDFwDrrKMq691d1W8Mwpou5GwxjwfBHKbTtgoh18HU6fRMIRmdQkBkPOmwDsyQYxhjosjUBxjlPQ9pa4AuAmBvClS3ZbcnoCuZusF+jYYyLRRUK/fqHBYMvBvz3bp0llnju2rljbZdhGNUxd4NhjBEnFLqWaGmcIx1SoxfcyIbY+pAJBHM3GJ1CF3sCMjCRYBgLSRszPKbrrp4LIXQtHLt2TtvNMYxm6Oy7CtrARIJhTIAT1y7sjDWhjjgoerKvkvOgal3mYjCMbmIxCYYxIboeo1CHMM5gcNmwiVYGJosexASC0Vm0185rhjGRYBhTYJqzHdaZx6Eor0HKMhA7NpcnwceCFI2uIz1t5TVSG0TOFpFL89c1InJpXn6giNztLft7b5tHi8jXReQKEflbEZG8fNe8vitE5EsicmDZ/k0kGMYEcdYEv5OehemR14O/sKwOS/mfCQTDKEdVn6uqB6vqwcBHgI96i690y1T1xV75O4EXAQ/LX0fk5S8EblbVhwJvBv66bP/dvzsZxpwzqSyITWhrtIEgNnLBmDG0U+6G3BrwHOCskvX2B/ZW1S+oqgJnAEfni48CTs/fnwP8hrMypDCRYBhTYlxZD6tSRZwUxR6MglkRjM6jtJlM6d4icrH3elGDFh0K3KCq3/XKHiQi/yMinxORQ/Oy+wNbvHW25GVu2Q8AVHUncCuwX9FObXSDYUwYf6TDNF0N4953auSDCQRjZmgv6PBGVT0ktVBEPgNsjix6taq62c6OY9CKcB3wQFX9iYg8GvhnETkIolHCLjCiaFkUEwmGMSVcroJpuRtS+02JB9fhV53yOZU8yTCMQVT1SUXLRWQX4FnAo71ttgHb8veXiMiVwMPJLAcHeJsfAFybv98CPADYkte5D3BT0b7N3WAYU6DrwyGLhIsfuFiHJZZsCmhjhuhUTMKTgG+p6oYbQUTuIyLL+fsHkwUoXqWq1wG3i8hj83iDEwD3w/sYcGL+/ljgs3ncQhITCYYxJU5cu3CjM+7SCIdxWTbMzWDMFAra01ZeLfA8hgMWnwB8TUS+ShaE+GJVdVaBlwDvBq4ArgQ+kZf/E7CfiFwB/CnwqrIdm7vBMDpA2DFP0g3hBEq4v3GmbjYMozqq+oJI2UfIhkTG1r8YeGSkfCvwW3X2bXcAw5givtsh7JAn2UHHREqb+7/q+ltaq8swJkZ33A1Tw0SCYUwZP8GS3zGPy5IQdv6x/fhlLkNjnfYss8wKKxy3dp7FIRgzSqdiEqaGiQTD6BC+UPA78zaf6qu6EeoKgxCLQTCM2cdEgmF0gBPXLkya+Mfhdmja+Zdtt5z/mfXAmAt6vXZeM4wFLhpGx3Edc2Hu1BqEQZFVgySLxIqfD8EEgjEXKOiMuwrawESCYXQYvwNvZSAV8VEM4b7C9+V1rnP82vkttdAwjK5gIsEwOsCZq8+YdhMGYhXCwEUbCmksHjrzQYdtYCLBMKZMFwSCI+V2CAWELxrMgmDMLSYS7PHAMLrIOOY8GMUakNrWBIJhzDdmSTCMKXP82vmF1oQqgYWprIk+dYMTU26GZZZMHBhzj6IWuIiJBMPoBE4o+BaE1CyLMdEQBhkWCYIi0ZHK0+C30zAWBhMJJhIMoyuELgZfIBRlYmySW6FK1kUfEweGsZiYSDCMjrDEEj16LA114JlYEJrFFaQsEoPrWCCiYQxgeRIAEwmG0Rl69Dhu7TzOXj2G3ggZEf1Jo4rw4yBMGBhGiMUkgIkEw+gM29kBMCQQqlgCHFUFApgwMIxSTCSYSDCMruA6eJfW+KzVowCGBIJ17oZhTAoTCYbRUfw5EM5aPYqrr7+Vk0wgGMaEMHcDWDIlw5gJbNIkw5g8qr1WXrOMiQTDMAzDMKKYu8EwDMMwAtSGQAImEgzDMAwjgsUkgLkbDMMwDMNIYJYEwzAMw4igDZOazRMmEgzDMAxjCHM3gLkbDMMwDMNIYJYEwzAMw4hglgQTCYZhGIYxhGIiAczd0CofX/2daTfBMAzDMFrDLAkt8sy190+7CYZhGEYrKD102o2YOiYSDMMwDCPEMi4C5m4wDMMwDCOBWRIMwzAMI4JZEkwkGIZhGMYQitKzjIsmEgzDMAwjhlkSLCbBMAzDMIwEZkkwDMMwjAg2BLKhJUFEThYRFZG3eWUiIqeIyLUicreIXCQiBwXbPUxEPpev8xYRWfKWPUhEzhSRLSKyLV/nAhF5VPPDMwzDMIz6KEpPe628ZpnaIkFEHgv8AfC1YNGfAa8E/gj4ZeBHwKdF5B7eOm8HPgw8DfgZ4Hl5nSvAp4H7AM8BHg4cC3wZuFfdNhqGYRiGMTq13A0isg/wfuCFwGu8cgFeAZyqqh/Jy04kEwq/DfxDvuo9gUuAy4BrgH3z8oOAhwBHqOoVedn3gM/XPB7DMAzDaAU1d0NtS8K7gHNU9bNB+YOAzcCnXIGq3g38O/A4b71XAxcCW8mEwel5+Y+BHvBsEbE4CcMwDGPqmLuhhkgQkT8AHgr8n8jizfn/G4LyG7xlqOqngPsBB6jqE1T1zrz8h8Afk1knbsnjFl4bxjQYhmEYhjE5RLXcnCIiPwP8J3Coqn4rL7sIuExVXy4ijwP+C3igqv7A2+49wP6qekSlxmTxC4cDvwI8E/g54PdU9X0F21xcpW7DMAyjVW6sem+fRUTkX4B7t1TdzJ6rqiLhBcB7gHWveJlsyu0emevgW8BjVPUr3nYXkJ2cE2s3LItz+CTwcFU9sO72hmEYhmGMRlV3wz8DPw8c7L0uBj6Yv/8OcD3wZLeBiOwGHErD4EPN1Mu3gL2abG8YhmEYxmhUChJU1VuAW/wyEbkTuElVL8s/vwV4tYh8i0w0/AVwB/CBsvpF5GDgL4H3Ad8AtgO/DvwecFaVNhqGYRiG0S5tjiR4A7A7WS6EewJfAp6iqrdX2HYLcBVZ4OKBZBaO7wOnAae22EbDMAzDMCpSKSbBMAzDMIzFwyZ4MgzDMAwjiomEBojIE0TkYyLyw3wOixcEyxvPYyEip4rId4N1D8j38/Gg/Ml5+YPHdKiIyEtF5GoR2Soil4jIod6y+4nI+SJynYi8T0R2H1c7gjaN7fznyy7K6w1fHxyx3asi8hURuU1EfiwiHxeRR85I218mIl/L236biHxBRJ7e9XZHjmMc886M65yfEqnz+q6325gvTCQ0Yy+y1NL/C7g7srzxPBbAZ4GHisgB3rqHk8VoPEFElr3yw4Dvq+pVox5QDBF5LvBW4PXAo8hGqnxCRB6Yr/Ja4KvAU8iGx/7JONoRYZzn3/EeYP/g9Ycjtvsw4B1kWUifCOwEPiMi/vwkXW37FuDPgV8CDiG7Tv9ZRH6h4+3eQMYw78wE2v7toM6fn5F2G/OCqtprhBfZCI4XeJ8FuA54tVe2O3A78Ide2VeAXyXLN/FO4KV5+R5kozue7637/8g65C3AL3vlnwfeM8Zj+xLwj0HZd4G1/P2HgePy9ycDb5j1858vuwh42wTavheZuHrmrLU939dNZB1K59sN7ANcSSbONvbV5bYDp5AlrIst62y77TVfL7MktM9I81io6l1ks18e7q17ONkP+nOuXET2Int6+LdxHISIbAIe7R9HzqfoH8frgNNEZAfZRF5vHkdbajLqPCKT5B5k1ryb888z0XYRWRaR55GJnM8zG+0e17wz4+bBkrnVrhaRD0rftdj1dhtzgomE9hlpHoucf6MvBg4EforsZrwhEoBfIxvCGt702uLeZE8fyeNQ1UuBBwIPAH5eVa8bU1vq0Mb5B3iRiNwRvF7aclvfClwKfGEW2i4iPy8idwDbgL8HjlHVr89Au8c278yY2/4l4AVkroI/yNvzeRHZr+PtNuYIm3FxfIRjSyUsU9XtDP/IIRMJf5ELhMOBL6vq3ZLNl/FGyWbKPBz4rqpuabvhAYXHoarrZNk2u8Yo5x/gbLIEXz4/bqdpICJvIhN6v5afw4GmhauHZVNq+7fJMqzuCzwbOF1EDvObFaw/9XZLNu/M68nmndlesGrn2q6qnxhokMgXyfLJnAh80a0WbDb1dhvzhYmE9nEd5mbgB175fUn/UEM+T/a0dlj+ughAVb8jIreTBY4dxphcDTk3kvnLNwfldY5jGrRx/gFuVdUrWmuVh4i8mSyA7HAdDDrtdNvzzsbVe7GI/DJZsOrr8rIutvtXyaxil4mIK1smCwJ+MZkJHrrZ9gFU9Q4RuRx4GFmqfJiBdhuzjbkb2udqRpzHQlW3kpmgD6cfj+D4HHAUWbzAuFwNrkO4BO84cp5Mw/k4JsTI53+ciMhbyeI3nqj5jKoenW57hCVgV7rd7n9mwvPOjIu8XT9LFrDY5XNuzBFmSWhAHjT40PzjEvBAyeafuElVvy8jzGPh8W/ASWQ3Yf9H/zmyFNjLDIqHcfAm4H0i8mWyqcBfTBYf8fdj3m8hEzr/e4hIaEXZrqo3jdDutwPPB44Gbvbqv0NV71BV7XDbTwUuIHtqvQeZ0DkMeHqX261jnnfGYxzn/DTg42TDn+9LFlOxJ3B6l8+5MWdMe3jFLL7Ibo4aeb03Xy5kw5euI4sq/hzwyJr7ODSv8z+C8p/Ny6NDo8ZwrC8FriFzf1wCPGHezz+Z+IrV/58jtjtWpwKneOt0te3vBb6XXwc/Aj4DPLXr7S7Y19u63nYya8e1ZEOifwh8BHhE19ttr/l62dwNhmEYhmFEsZgEwzAMwzCimEgwDMMwDCOKiQTDMAzDMKKYSDAMwzAMI4qJBMMwDMMwophIMAzDMAwjiokEwzAMwzCimEgwDMMwDCPK/w/22hZSAUvpCgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1152x432 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig = plt.figure(figsize=(16, 6))\n",
    "\n",
    "\n",
    "#greyscale = raster.mean(dim=\"band\")\n",
    "\n",
    "# Plot on a map\n",
    "ax = plt.subplot(projection=ccrs.PlateCarree())\n",
    "raster.plot(\n",
    "    ax=ax,\n",
    "    x=\"lon\",\n",
    "    y=\"lat\",\n",
    "    transform=ccrs.PlateCarree(),\n",
    "    cmap=\"cubehelix\",\n",
    "    shading=\"auto\",\n",
    "    add_colorbar=True,\n",
    ")\n",
    "\n",
    "\n",
    "\n",
    "g1 = ax.gridlines(draw_labels = True)\n",
    "g1.xlabel_style = {'size': 14, 'color': 'k'}\n",
    "g1.ylabel_style = {'size': 14, 'color': 'k'}\n",
    "\n",
    "\n",
    "g1.xlabels_top = False\n",
    "g1.ylabels_right = False\n",
    "\n",
    "\n",
    "\n",
    "plt.title(\"chirps-dekad/africa_3-monthly/merged\",fontweight='bold', size=14)\n",
    "\n",
    "plt.tight_layout()\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "id": "98287fdf-1a6e-44da-8bd8-68ca6865c460",
   "metadata": {},
   "outputs": [],
   "source": [
    "#geom_dem, info = mask (raster_array, [geom],  crop=True, invert=False, all_touched=False) \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d8903db9-bdf4-442b-8b32-bbfe8d9d9b89",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}