examples/GulfStreamImageOverlay.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.3.0+207.g853b5c0.dirty\n"
     ]
    }
   ],
   "source": [
    "import os\n",
    "import folium\n",
    "\n",
    "print(folium.__version__)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import cartopy.crs as ccrs\n",
    "import matplotlib.pyplot as plt\n",
    "import cartopy.feature as cfeature\n",
    "\n",
    "LAND = cfeature.NaturalEarthFeature(\n",
    "    'physical', 'land', '50m',\n",
    "    edgecolor='face',\n",
    "    facecolor=cfeature.COLORS['land']\n",
    ")\n",
    "\n",
    "\n",
    "def make_map(bbox, projection=ccrs.Mercator()):\n",
    "    fig, ax = plt.subplots(\n",
    "        figsize=(8, 6),\n",
    "        subplot_kw={'projection': projection}\n",
    "    )\n",
    "    ax.set_extent(bbox)\n",
    "    ax.add_feature(LAND, facecolor='0.25')\n",
    "    ax.coastlines(resolution='50m')\n",
    "    return fig, ax"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import numpy.ma as ma\n",
    "from scipy.io import loadmat\n",
    "\n",
    "fname = os.path.join('data', 'mercator_temperature.mat')\n",
    "merc = loadmat(fname, squeeze_me=True)\n",
    "\n",
    "x, y = merc['x'], merc['y']\n",
    "bbox = [x.min(), x.max(), y.min(), y.max()]\n",
    "\n",
    "level = 10  # 10 meters temperature.\n",
    "data = ma.masked_invalid(merc['temp'][level, ...])\n",
    "\n",
    "fig, ax = make_map(bbox=bbox)\n",
    "cs = ax.pcolormesh(x, y, data, transform=ccrs.PlateCarree())\n",
    "fig.savefig(os.path.join('results', 'GS.png'), transparent=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from PIL import Image, ImageChops\n",
    "\n",
    "\n",
    "def trim(img):\n",
    "    border = Image.new(img.mode, img.size, img.getpixel((0, 0)))\n",
    "    diff = ImageChops.difference(img, border)\n",
    "    diff = ImageChops.add(diff, diff, 2.0, -100)\n",
    "    bbox = diff.getbbox()\n",
    "    if bbox:\n",
    "        img = img.crop(bbox)\n",
    "    return np.array(img)\n",
    "\n",
    "\n",
    "with Image.open(os.path.join('results', 'GS.png')) as img:\n",
    "    img = trim(img)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "image = np.zeros((61, 61))\n",
    "image[0, :] = 1.0\n",
    "image[60, :] = 1.0\n",
    "image[:, 0] = 1.0\n",
    "image[:, 60] = 1.0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><iframe src=\"data:text/html;charset=utf-8;base64,PCFET0NUWVBFIGh0bWw+CjxoZWFkPiAgICAKICAgIDxtZXRhIGh0dHAtZXF1aXY9ImNvbnRlbnQtdHlwZSIgY29udGVudD0idGV4dC9odG1sOyBjaGFyc2V0PVVURi04IiAvPgogICAgPHNjcmlwdD5MX1BSRUZFUl9DQU5WQVMgPSBmYWxzZTsgTF9OT19UT1VDSCA9IGZhbHNlOyBMX0RJU0FCTEVfM0QgPSBmYWxzZTs8L3NjcmlwdD4KICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2Nkbi5qc2RlbGl2ci5uZXQvbnBtL2xlYWZsZXRAMS4yLjAvZGlzdC9sZWFmbGV0LmpzIj48L3NjcmlwdD4KICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2FqYXguZ29vZ2xlYXBpcy5jb20vYWpheC9saWJzL2pxdWVyeS8xLjExLjEvanF1ZXJ5Lm1pbi5qcyI+PC9zY3JpcHQ+CiAgICA8c2NyaXB0IHNyYz0iaHR0cHM6Ly9tYXhjZG4uYm9vdHN0cmFwY2RuLmNvbS9ib290c3RyYXAvMy4yLjAvanMvYm9vdHN0cmFwLm1pbi5qcyI+PC9zY3JpcHQ+CiAgICA8c2NyaXB0IHNyYz0iaHR0cHM6Ly9jZG5qcy5jbG91ZGZsYXJlLmNvbS9hamF4L2xpYnMvTGVhZmxldC5hd2Vzb21lLW1hcmtlcnMvMi4wLjIvbGVhZmxldC5hd2Vzb21lLW1hcmtlcnMuanMiPjwvc2NyaXB0PgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL2Nkbi5qc2RlbGl2ci5uZXQvbnBtL2xlYWZsZXRAMS4yLjAvZGlzdC9sZWFmbGV0LmNzcyIgLz4KICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9tYXhjZG4uYm9vdHN0cmFwY2RuLmNvbS9ib290c3RyYXAvMy4yLjAvY3NzL2Jvb3RzdHJhcC5taW4uY3NzIiAvPgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL21heGNkbi5ib290c3RyYXBjZG4uY29tL2Jvb3RzdHJhcC8zLjIuMC9jc3MvYm9vdHN0cmFwLXRoZW1lLm1pbi5jc3MiIC8+CiAgICA8bGluayByZWw9InN0eWxlc2hlZXQiIGhyZWY9Imh0dHBzOi8vbWF4Y2RuLmJvb3RzdHJhcGNkbi5jb20vZm9udC1hd2Vzb21lLzQuNi4zL2Nzcy9mb250LWF3ZXNvbWUubWluLmNzcyIgLz4KICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9jZG5qcy5jbG91ZGZsYXJlLmNvbS9hamF4L2xpYnMvTGVhZmxldC5hd2Vzb21lLW1hcmtlcnMvMi4wLjIvbGVhZmxldC5hd2Vzb21lLW1hcmtlcnMuY3NzIiAvPgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL3Jhd2dpdC5jb20vcHl0aG9uLXZpc3VhbGl6YXRpb24vZm9saXVtL21hc3Rlci9mb2xpdW0vdGVtcGxhdGVzL2xlYWZsZXQuYXdlc29tZS5yb3RhdGUuY3NzIiAvPgogICAgPHN0eWxlPmh0bWwsIGJvZHkge3dpZHRoOiAxMDAlO2hlaWdodDogMTAwJTttYXJnaW46IDA7cGFkZGluZzogMDt9PC9zdHlsZT4KICAgIDxzdHlsZT4jbWFwIHtwb3NpdGlvbjphYnNvbHV0ZTt0b3A6MDtib3R0b206MDtyaWdodDowO2xlZnQ6MDt9PC9zdHlsZT4KICAgIAogICAgICAgICAgICA8c3R5bGU+ICNtYXBfMTlkZGYzN2Y3M2UwNDI2YmE2OGZmOWQ4YzI4ODc4NmMgewogICAgICAgICAgICAgICAgcG9zaXRpb24gOiByZWxhdGl2ZTsKICAgICAgICAgICAgICAgIHdpZHRoIDogMTAwLjAlOwogICAgICAgICAgICAgICAgaGVpZ2h0OiAxMDAuMCU7CiAgICAgICAgICAgICAgICBsZWZ0OiAwLjAlOwogICAgICAgICAgICAgICAgdG9wOiAwLjAlOwogICAgICAgICAgICAgICAgfQogICAgICAgICAgICA8L3N0eWxlPgogICAgICAgIAogICAgPHN0eWxlPgogICAgICAgIC5sZWFmbGV0LWltYWdlLWxheWVyIHsKICAgICAgICBpbWFnZS1yZW5kZXJpbmc6IC13ZWJraXQtb3B0aW1pemUtY29udHJhc3Q7IC8qIG9sZCBhbmRyb2lkL3NhZmFyaSovCiAgICAgICAgaW1hZ2UtcmVuZGVyaW5nOiBjcmlzcC1lZGdlczsgLyogc2FmYXJpICovCiAgICAgICAgaW1hZ2UtcmVuZGVyaW5nOiBwaXhlbGF0ZWQ7IC8qIGNocm9tZSAqLwogICAgICAgIGltYWdlLXJlbmRlcmluZzogLW1vei1jcmlzcC1lZGdlczsgLyogZmlyZWZveCAqLwogICAgICAgIGltYWdlLXJlbmRlcmluZzogLW8tY3Jpc3AtZWRnZXM7IC8qIG9wZXJhICovCiAgICAgICAgLW1zLWludGVycG9sYXRpb24tbW9kZTogbmVhcmVzdC1uZWlnaGJvcjsgLyogaWUgKi8KICAgICAgICB9CiAgICAgICAgPC9zdHlsZT4KPC9oZWFkPgo8Ym9keT4gICAgCiAgICAKICAgICAgICAgICAgPGRpdiBjbGFzcz0iZm9saXVtLW1hcCIgaWQ9Im1hcF8xOWRkZjM3ZjczZTA0MjZiYTY4ZmY5ZDhjMjg4Nzg2YyIgPjwvZGl2PgogICAgICAgIAo8L2JvZHk+CjxzY3JpcHQ+ICAgIAogICAgCgogICAgICAgICAgICAKICAgICAgICAgICAgICAgIHZhciBib3VuZHMgPSBudWxsOwogICAgICAgICAgICAKCiAgICAgICAgICAgIHZhciBtYXBfMTlkZGYzN2Y3M2UwNDI2YmE2OGZmOWQ4YzI4ODc4NmMgPSBMLm1hcCgKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICdtYXBfMTlkZGYzN2Y3M2UwNDI2YmE2OGZmOWQ4YzI4ODc4NmMnLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAge2NlbnRlcjogWzM2LjYzMTY1MjE2NTYsLTY5LjUwMDAwMDY2NjZdLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgem9vbTogNiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG1heEJvdW5kczogYm91bmRzLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbGF5ZXJzOiBbXSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHdvcmxkQ29weUp1bXA6IGZhbHNlLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY3JzOiBMLkNSUy5FUFNHMzg1NwogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB9KTsKICAgICAgICAgICAgCiAgICAgICAgCiAgICAKICAgICAgICAgICAgdmFyIHRpbGVfbGF5ZXJfNTg0ZWRmMjNlYzVmNGNjOTg1NzgxODMwMGM1ZTA3MjAgPSBMLnRpbGVMYXllcigKICAgICAgICAgICAgICAgICdodHRwczovL3tzfS50aWxlLm9wZW5zdHJlZXRtYXAub3JnL3t6fS97eH0ve3l9LnBuZycsCiAgICAgICAgICAgICAgICB7CiAgICAgICAgICAgICAgICAgICAgbWF4Wm9vbTogMTgsCiAgICAgICAgICAgICAgICAgICAgbWluWm9vbTogMSwKICAgICAgICAgICAgICAgICAgICBub1dyYXA6IGZhbHNlLAogICAgICAgICAgICAgICAgICAgIGF0dHJpYnV0aW9uOiAnRGF0YSBieSA8YSBocmVmPSJodHRwOi8vb3BlbnN0cmVldG1hcC5vcmciPk9wZW5TdHJlZXRNYXA8L2E+LCB1bmRlciA8YSBocmVmPSJodHRwOi8vd3d3Lm9wZW5zdHJlZXRtYXAub3JnL2NvcHlyaWdodCI+T0RiTDwvYT4uJywKICAgICAgICAgICAgICAgICAgICBkZXRlY3RSZXRpbmE6IGZhbHNlLAogICAgICAgICAgICAgICAgICAgIHN1YmRvbWFpbnM6ICdhYmMnCiAgICAgICAgICAgICAgICAgICAgfQogICAgICAgICAgICAgICAgKS5hZGRUbyhtYXBfMTlkZGYzN2Y3M2UwNDI2YmE2OGZmOWQ4YzI4ODc4NmMpOwoKICAgICAgICAKICAgIAogICAgICAgICAgICAgICAgdmFyIGltYWdlX292ZXJsYXlfMmJmNDRhYTg2ZTRiNGQ3NDkyYzRkNzQxYTViYTAzNTMgPSBMLmltYWdlT3ZlcmxheSgKICAgICAgICAgICAgICAgICAgICAnZGF0YTppbWFnZS9wbmc7YmFzZTY0LGlWQk9SdzBLR2dvQUFBQU5TVWhFVWdBQUFEMEFBQUE5Q0FZQUFBQWVZbUhwQUFBQVcwbEVRVlI0MnUzUG9RMEFBQXpEc1A3LzlFYUc5MEFOSW9VNmsweGJ1VWxSME5EUTBORFEwTkRRME5EUTBORFEwTkRRME5EUTBORFEwTkRRME5EUTBORFEwTkRRME5EUTBORFEwTkRRME5EUTBORFEwTkRRMEMrNnFnV3pZV3BwQTlnLzdBQUFBQUJKUlU1RXJrSmdnZz09JywKICAgICAgICAgICAgICAgICAgICBbWzM0LjUsIC03MS43XSwgWzM2LjYsIC02OV1dLAogICAgICAgICAgICAgICAgICAgIHsKICAiYWx0IjogIiIsCiAgImNsYXNzTmFtZSI6ICIiLAogICJjcm9zc09yaWdpbiI6IGZhbHNlLAogICJlcnJvck92ZXJsYXlVcmwiOiAiIiwKICAiaW50ZXJhY3RpdmUiOiBmYWxzZSwKICAib3BhY2l0eSI6IDEuMCwKICAiekluZGV4IjogMQp9CiAgICAgICAgICAgICAgICAgICAgKS5hZGRUbyhtYXBfMTlkZGYzN2Y3M2UwNDI2YmE2OGZmOWQ4YzI4ODc4NmMpOwogICAgICAgICAgICAKICAgIAogICAgICAgICAgICAgICAgdmFyIGltYWdlX292ZXJsYXlfZjk1MzJiN2Q2NWM2NGIxYzg3OGU0MGEzMzQ4YmE4OTUgPSBMLmltYWdlT3ZlcmxheSgKICAgICAgICAgICAgICAgICAgICAnZGF0YTppbWFnZS9wbmc7YmFzZTY0LGlWQk9SdzBLR2dvQUFBQU5TVWhFVWdBQUFaa0FBQUZIQ0FZQUFBQnprVEdSQUFCY2NrbEVRVlI0MnUyOUM3UmU4NTMvMzdRNmJkRUxSdXNTY3IvZmhDQVM5MHNTZzFKMUswSXFya1hkYVVKZEVpR1I1SlNwQ1NhVFFhV1ZRVnN5U3VxSEtxMmZ2N1orREtVbVRWSjFxY21NWmJFc0s1WXgrKys5VjU4enovZTFuL044enM1enpzbTV2UGRhN3hYT3Z1OW5QNS9YOC8xK2JwLzR6R2MrazMzaUU1K3dMTXV5ckxiV256NmgvOWg3NzcwdHE4dHFyNzMyeXJiZmZ2dHM1TWlSV1d1WC9mYmJMOXRtbTIyeTdiYmJMcHM2ZFdyejMvZnZkVVNpeVZ1Y2ttaks0RXZxcTgrNXFRWmNtR2p5bU84bW1qVHV5a1FIN0h4Vm92MTNtNTFvMzcydlNiVFAvdGNtMm12S3ZJSTZjamxnL0t4RSswK1luV2l2eVhNVDdYYlVna1RqcGkxTU5QYTBwa1RETDBrMTlMdlFaZVUwNk9wVS9lY3ZURFRxdktaRXZMNGRUMjVLdE90eEN4UHgvaWh1djlQMC85WFkwNXNTN1hCR0txNGY5ODJGaVhpdUNVZWtHditOVkx5V25VOU1OZWFzcGtTanY1MXE1SGxGL1JVMGhvelY5YlhUVGp0bFgvemlGMXR0RE45NDQ0MU1vL2l0dDk0Nkd6Tm1qQ0ZqeUJneWhveGx0YXc5OTl3eisrUW5QNW05L2ZiYnJUYUl4eHh6VFBiNXozOCsyMnl6elF3WlE4YVFNV1FzcTc0RWpMdnZ2cnRWeHZEZGQ5L05kdHR0dCt6TFgvNXlEcWQxNjliVjNHN1NKaWNrbXRMdmdsUkR2NU5xdTNOUzlmNTJxbUV6RWtYUW9aSGViODg1QlhYbWhSQ2MrTFg1aVhZK1lXR2lIYjdWbEdqNGQxSU51elFWMXhNaWc2OUtOZVNLVkFPdlRjWDFPNTZTcWhvQ0VxOS90eU1YMUZjQW1iR25OalVyTXVLanowNUZLQkFxZlBhRXpDN0hMMHpFZTZ1K05vblBocCtkWk1oWUcwekRoZzNMRGZ6dXUrL2Vac2VVZjBWNjhjVVhXelI2SDMzMFVYYmxsVmRtRzIrOGNYNytYWGZkTmZ2c1p6K2JQZjMwMDRhTUlXUElHREpXZDFHZlBuM3lGL0FyWC9sSzg5LzIyR09QYkpkZGRza2QrZXM3WmRhM2I5L3NiLzdtYjdJVFRqZ2hlL3p4eDV1TjNXT1BQWllkZHRoaDJSZSs4SVY4eENNL1RQVUlhUG55NVlhTUlXUElHREpXZDlFT08reVErMEkrL2VsUForUEhqOCtucmpiWlpKUDhwWnd3WVVLcklzbzBDaEZZdUc3SEhYZk10dDEyMjl5eHIyTis2VXRmeXNHalVjN1lzV01MRU50ODg4Mnp4WXNYR3pLR2pDRmp5RmpkUllLS2ZDRUhISEJBRGh0TldXMjU1WmE1d1cvTmlFWGh4NS82MUtmeWZlcUJTRkZuQWhyQm92L1hOWXdiTnk2UFNydnBwcHRxR3NuSm01MmNxT0RJSDNWcElxNmYwdmU4UkpPSHowdzBhY2ZMRXhFcWRLUkxYV21aOFBYNWlXaklhRWdKaVdFelVvMjRLQldoTTNoV3FvSFgxQmNkLzd5ZW5VNWFXRmU3VElWd2YxeS82N0VMNmlvQkxoejlJODlQeFd1bDBXY2dRT0ZhSVVJbHVoZENyZGJ6TVdTc0RTWkJRa2IrL2ZmZnp3WU5HcFI5NHh2ZnlIcjM3cDM3YWlJNGFXU2lmZTY0NDQ1czAwMDNiUldVQkpPaFE0Zm1veG1CVElEU0tFcFRaUUxjTysrOFk4Z1lNb2FNSVdOMUI0MGFOU3FmeGxxN2RtMXVpSTQ2NnFqOFpaVGhsMSttcGYwMERTWWdISHJvb2RtSEgzNll2ZkRDQy9rMFdLMFJqSUFpYUFsSUdqRUpSc09IRDgrT08rNjRiTkdpUmRtYU5XdGFaU1FOR1VQR2tERmtyQzRtWmVjZmZ2amh1UkZhdFdwVnR0RkdHMlc5ZXZYS0lkUFNQb0tHUmg1ejVzeEp3cEMxWC9WVW1QdzVXMnl4UlI1UThNMXZmak83N2JiYnNqZmZmSE85amFRaFUzK0pjaXNNR1VQR2tMRTZYQnBoL04zZi9WMCtqWFhhYWFkbFcyMjFWVE5vT0pMUk5vS1NSajRQUHZoZ1l1QWVlZVNSN0hPZisxd09HWTF5TkJVbUVNblAwMUxlUzZOTHdYRS8rckpFaGVUTVFSY25vcU9meVpWTVZwUTY4eEpCaG81eU9xc0ppV0V6VTlHdzB0RS81SExveWxTRDVrQ3pVNDI0TUJVTk01TXZxUUpFOER3S0NZOEVDN2F2WGtjanoyZkpaRXdtcm5JOWt5a0pzUUpBY1g1ZU81OFZqeWNaTXRZR2thQ2ljaTU2OXdRRytXRTBwYVdwTUVXTVNmcTduUHFDaUtCVWEzcExQaGxObHdsQWtxYlJubnp5eVhZMXFvYU1JV1BJR0RKV0o1ZW1zd1FIdlh1Q2cwWWgrbTlCUmFISEdvMGNlT0NCMmZ6NTg3TVZLMWJrL3BkYXl3Y2ZmSkFkZWVTUjJkS2xTL01reTQ1WURCbER4cEF4Wkt3dVVQNmx1aHk0RWpPVmdUOWp4b3c4V2ZLZWUrN3B0RWJWa0RGa0RCbER4dXJFbWpoeFl1NS9PZXVzczdKcDA2YmxPU3BEaGd6Sm5uLysrZElHN29ZYmJzaHpiSlIwZWVxcHA3WTQ0bW5QaFk3OHlTTWhRR2kvUGE1T05QSHcrYW0rVmxSblhxS0VQQmFZSkVSWTRMSldrY1ZxTVRDQTBCbCtjU3BDak5DaTRhYlJqeUJTdWdBbURIa0JERlVxSkg0QzRLRWpubEJDTXVXWU0xTXhrSUEvQ0FySm9SZWtZdkZTeVpDeE9yenZ5OS8rN2QvbXp2eEdsMXR2dlRXZmNwTS9Sd21YQ2xYV2RGc2xMTG8xaTZiYkRCbER4cEF4Wkt4dUFCZUJRQW1Za2hJdzEzZlJhRVZoeVFMTTZOR2preEdTL0R1WFgzNTVxNDR6ZCs3Yy9BdlFxQy9Ia0RGa0RCbER4bXBIZUNndlJmWEFWTHBGeFMxbDdQVjMvYmV5OHBXdklrZStraUZWbzJ6bHlwWHJiZENVNzZJSU5LbC8vLzc1S0VhUlp6cTJjbXdVdGZhTFgveWlWY2VxSklBcWFNQ1FNV1FNR1VQRzJnQytFMlhtVnhlclZFa1haYzBydkZnT2VobDIrVmRVcGtWbDh6V1NxT1M3Q0FTcWRIelJSUmRseno3N2JKc1lOSTA2RGpua2tEemNXU01aUWV6RUUwL01lOGkwZG5Ta05nRG5uMzkrdnIrbTdrYU1HTkcyZ1FCQndVc1dqS3hWUkxFckxUUkVoRXJCVVk4Q2xEUktkRTdUTWMva1RCbytRbzdIbzVHbk01dVFHSDlNS3E2UGtrc0pPU1ovc2tCbmRaQkM5S3hHblp1S253V2ZSU0VRZ0pBaFZNNUpSZUR6c3lIZ0pVT21oMmpublhmT25ldXNiaXpmaVBxNTZPL0tUZEYyR3BVSUlvSkZKVWxTSXdWRmYybEVvaEdFU3JOb0JQREtLNi9VTkR3eStPMFpVcXpSVU5uam4zVFNTWG1RZ01DbzBaWGdVZ2xDVUN0bVE4YVFNV1FNbVc1YmpWakd2VDFMNnV0WCs0QUJBM0tEcXYrV29kVm9vMUl2VEg5VDBxTmdJa004YWRLazdMbm5ubXQyanYva0p6L0pzKzAzUlBSV1d5MGFVUW1ldWwvZG8vck82UDgxS2xMbEFFUEdrREZrREpsdUtZMHc5QmtJQm0zZDgxNCtDd0hrKzkvL2ZtNFFWR2xZMDBXdnZmWmE5dDU3NzNWWUFtTm5XdDUrKyszc3hodHZ6UGJiYjcrOGlabW0ydjduZi83bmY2R3gxYmRTc1IvTXgrckprT0c4dmlGanlCZ3luVnlxMWFYUFFDT0wxalRyaWh6eDhvUElHYTZwTFdYV0wxdTJMUE5TWW1SU0VqSlJJQUFocys5ZWN3cnFTZ3NoVTBpV2hDRXNKRXRDaEJJTkk0L1BBcG1GaEVSY0g1M3BMQUpKaUJBYU81N2NsSWlHbGdVNSs5MndJTkgyaTY5TDFPZWY1Nlc2ZFc2cXFuV0ZmZjloZnFMK0N4WW1ZcUlwbnpXREpnZ3g3azhnaHBwVGxDSFRTYVE2WHBxMjBmVE4rc0pGVllybE4xRytpQ0tuZnZXclg3V0pVZEYwbVVEVkZqa2xob3doWThnWU10WUcwTWlSSS9NUFE5Rk9aZUdpS1RGRmRjbVpmY3N0dDdTWk1kR1UycFFwVS9Jb01ZMktWSkxma0RGa0RCbER4cERwb3RGZitoeFVZcVcxY0ZFb3NSejFBcE44REczcFg3bjMzbnVidzVFRm1QdnZ2OS9UWlQwRU1uc2VPQzlSdlI3MGhvd2hZOGgwRVNsSlVaK0RwcnNpdUNqMFZnVW01VzlwYW1wcVU3am9XR2VjY1VhelAwYzF3ZGpEUlkzQ3FoZEZuRDN3d0FQWjZhZWYzaTBDQ1NaL2FYcWlLWDk3V2tIZGVZa2dRMmN3b1VKSE53MFoxeGVLTnRKeEh6aW5DMUNCV0xDUlVPTDkwQkQzbjc4d1VaK2I1NmU2Zlc2cTIxSnQvMC96VWdFY2hmVkwvbGNFVXZXNlhJVFFUZk1UOWYzK2drUUVZTisvVDlYblJnaFFLMmhSTEVPbUV5VSs2bk5RR0hGTGNOR1VtaElnRlg0N2I5NjhOamZveWhWUjhxVkdSd29jRUdBZWV1aWhIQ0lhS1NsQVFldDBuVWNmZlhSZU8wekFxNjZtM0VpNUdFUEdrREZrREJtcm5hUndZMzBPeXBxdnRVNDlWalI2bVQxN2RwdkNSUUJSUG94R0w0S0t6aVBBYUtyc3NjY2V5MnVFYVVTamMvZnIxeSt2UDZaUmwyQllBWXZBcCsyNnkyTElHREtHakNIVDdhU3MrMHBmRlk1d05ISlJJdVYvL01kL3ROcFFLQk5mSXcyVlQ1aytmWHAyekRISDVPVllsQnN5YnR5NFBJcE5vNUpLZVJaVkx4WkFCQmxsK1QvOTlOTjV0SnFjL2dvczBFaXFBcnhLRkpzaTRwWXNXZEx0akd4M2g4eWVCMStYaUQ0SlE4YVFNV1M2b1NxZElUWDlWQTBlQlFMbzMrcndZWTFrbm5qaWlkemdhM3JyekRQUHpNT1ZwZHR2dnoySGhSejJDbVZXeFdNQlpmdnR0ODhCcHBJd3F2ZWwvUlJzb0Ntd3l2a0dEaHlZNy9Qd3d3L253UVFDVG1WOWRZaTBJTlNXVVd4ZU9oZGs2T2d1UUlVRktKRWN5WUtYaFlSQ0pHY1NLbEdSUm00ZjdjOUFnY0ZYcGVyWHREQlJuMXV1U3hWQmhJYWVFQ3FybTBxbzdMRUlBYTdIdlJPQWhhQUZQcHZicnkzSWtPa2trdUhYNTZES3haVVJnNHBPQ2hpYTBqcjU1Sk9iYTRkVnlzQUlHb01IRDg1SEZKcnEwaFNXb0tSeUtXV1RPZ1VUaldvMGNoR2dCQ09CcFZJMW9CSWl2V2pSb2g1WkpjQ1FNV1FNR1VPbVM1ZkxyL2czOU44eStKcStVcmtaVFVkTm5UbzFuOW9TY0FTUFNpbjk5cGl5MDFTYWpxOFJsYUFpNE9qZnRnNlI5bUxJR0RLR2pOWEJpWmdhbFFnaW1ySVNZRFNGcHRHRjFtbDZxeU5nVngzQkpuK082cHg1NmRyTGhLL1BUMlRJR0RLR1RBOGJ4VlFpdFRSNnFUVGUwclNZcHNFMHF0QzBtUDV0ei9CcCtXTVVRU2JBekpvMXk2T1diZ3daT3NhWkRGa29lbmhKZlJFeWhBcU5Pby9QWk1sQ2djM2dmTkg1Qjh4TFJlZDF3ZkVlT09yNy9DTkVLTkhRNDN4MHRsTjB6aWVPZWpqeUM0NzRBQm9GU1B6ZzJrUUQvbVYyb3FFL3VUTFJxT1hmVFRUMlo1Y21HdmZBaklJTW1RMHNWVjZ1VEpWVklLTVE0b3B2cHIzQ3BUVXlVakNBL0Q2YUV0TUk1cWFiYnVyU3BmeTlHREtHakNGandkaHJXcW9DR1RuejI3T3ZqSTR0c0ZSQ2xvODg4c2pzaGh0dXlHdVVlVEZrREJsRHhwRHBadE5rQ2kvV1ZOaXh4eDZiVDQvSjZkNGVVMkdLWEZQVW1TTFFEampnZ0R4RTJVdm5YUTRZUHl2UmZudGNuV2pmdmE5SnRQdWgxeVZpMzNkRHhwQXhaSHFZMU9aWWVTaWFGdnZOYjM2VGh3MnJCMHhiQWt6VFlaWFd5Y3B2VVkyemRldlcyWUwzUU1nd0daRkd2RkNBRWhDZzQ1eU9meVl2RWlvVUF3bVlIRWxvY1Arb0tHUC9oUXNURVJKOWwxNlQ2bzVVQldkM0JCVUt6dlVDTEw2M01CVUNENUxrU0VDamNHMXd2UGY3NFp4RUErK2FuYWd0b0ZGV2hzd0c4TUZvMUtJRXlZTVBQamdmV1NqcHNSSkNyTkl0bXRKU1MrYldoQndMSmhvTnliY2lxT2k0R3JHb0ZNd3BwNXlTclZtenhsYmJrREZrREJsRHBxZFVXaFpnbEZRcG40aEdNQUtDOGw4MCtsQStpcElnQlprS2VMaS9FaU1GRWtXQ0tVaEFFV2dISFhSUWR2MzExemVYajJIVlpDK0dqQ0ZqeUJneVBVQWFhUWdxR3EzSTZTOEh2UHd5QW9XZ1V5a2hVeWs2V2ZtNzhtWXFqY05Vdyt5NDQ0N0w3cnp6enV5OTk5N3Jsa1oyNzBsekUrMDFHWm95TDlIZUI4eE5SS1BNZmk1U1J5NzdUNWlkaU5mTHZCVkR4cEF4Wkt6MXlxYlh5RU93RUdCVWJWa08rU3V1dUNJYlBYcDBkdkhGRitlakZJRkVJNW5ERGpzc216RmpSdTVMV2JwMGFWNlhyS2ZrcnZSMHlOQlJ6eVphTzUyVWlvNzhRa0ZKR0hVMjNTcEFCTW1TTEdnWk5iNHFGTVFFTkFnMXJpZEVCbDZiaW81OEdtMGE0Z2dxQldmNFAxNVhUb1FNaWt3V29ISjlxb0x6dnFxQkdhKzEvN0tyRXczNThWV0podDk3UlNKQ1pQUy9YcFpvelAycHVEMlBOK3luOVRYeXZzc0xNbVE2U0pYaWtocWhhS3BMeVk4Y2lhZ1hpNU1nRFJsRHhwQXhaS3oxZ295bXc4YU9IWnRkZSsyMWhva2hZOGdZTW9hTTFiWmh5L0xCcU5Oa2RTaXhZUFBTU3kvMUtJanNkdVNDUkRTYUhRR1o2dU54ZS9wRXBIckx2dnRjazJpUHIxNlh5SkF4WkF3WlE2WkRraThWRmFZcE14VzlWS2Evd3BrcnBXVFdybDFyeUZRMHZaeDJQakhWK0c4c1NEVHg4UG1KQ0szMmhnd2hNZTZicVhnL0xQaFljTlJma0twUXNCSlFHRHdyMWNCclV0Rm9ENWliaXZzUHVSd0traStIWEFsaGZ4Ni9BQmxjTHgzbFJXYzRpalN5b09XUytvM0JxbzE4TGpZT1kzSm1CQms0L3V0QkpWZFZjbVFFbFFnaVhFOG9ETHBuVmlJbVl6S1FnTUR1ZitlY1JEeWVaTWhzb0NBQVZWaFcxV1U1K3dXYTg4NDd6NUF4WkF3WlE4YVFzZHFuWVpucWloa3lob3doWThnWU1sYWJTeVZsbEVQVFhaWmRqMTJRQ2ozaXV5SmtKaHl4b0Zuc3YyTElHREtHakNIVHFhV01meVZjOWhUSU1QbU9qbWc2cWlPTk9qY1ZtMWFOT1NzVnowK2pUNGhVQTZZV1pMZy9DejdTTWMvMXZGNUNoTW1SaEFpTk1pRVJRWVJHbkpDSUhQMU1ob3lPeC8wTFRjNEFJUjZ2VUlDU0VDRjAySWlMa01EK0JVZ0JBaUZrQ0JGQXA3QS9FaXFySGYwMCtyV01lRjFJL0doT0lrSWl2TmVvV0NpZUhjOHZHVEtkUU1PSEQ4K256QXdaUThhUU1XUU1HYXZOcFFLWGt5Wk5NbVFNR1VQR2tERmtyTGFYa2pRWExWclVaYUZDSTJ2SUdES0dqQ0ZqeUhRU3FhUy9lcjZvcEV4M2dRd0xNdExJMDVFOSt1eFVoQUNOZEJRSXNPUEpUWW00ZitINmtMd1lyZWY5OFBxNWZlRjZBQmxDbG81OGF0RHNWSFNjRTBLRnhsNHpvV0E5b1VDakgwSUtFSW5PVDZnUm9tenlGVFVGQ3h0L3NRaGxFQmdRTmpFakJLUGoxNEZjNUhndkhlUVFCQ21FejQ3UFBnQ29aTWhzWUNsUFJuWE11dkppeUJneWhvd2hZOGgwVXFuaTh1R0hIMjdJR0RLR2pDRmp5Rmh0TDdWaG5qOS9mcGVDQ3Z1VkdES0dqQ0ZqeUJneW5WVHFML1BrazA5MmFjZ1VqQ2dnUWFQTlpFT3VweEVteEpoOEdScnRTeUFVZ0l5TWVtSDdZSDhXZ0N4QUNZRUxOTElzV0VralRjaHcrMElCVEFRTzhQcWpwbU4wN1BQOGtlT2UyL2VmdnpBUmt5c0hYSmNxYWxMRzdXa0lhV2hwcUduSTZZaVBJQk9wRUhoQXcxOHYrWk1BS2h1RUFNQ3hvVnBwUUFlUUtSei9ZeGt5Rzdob1pxOWV2YkozMzMzWGtERmtEQmxEeHBEcGFRRG8zNzkvL205N1JaV3A5TCtlZlZkYkRCbER4cEF4WkF5WkJxVlMvSG91N1FVWjFTcjc2OFB2Y2dVdERSbER4cEF4WkF5WkJqUmt5SkQ4d1NqeXF6Mk92OGNlZStROVpIUU90V1B1Y3BCQndVY2FmVUtIVUtHam5KQ0pIUGxNVnFUanUxRFFFUVVhYVVRTFlzRkdHcldTeVl5RUFPK2ZVQ0lrYUtSNVB3WElzWUFtSUZ0d3ZFZUJBQkFoUWlOUFJ6aUxQQTc5eVpXSm9xS05CV2MzQzFvQ0NvUHVucFdJVGJTaW9wSnMxRFh3cnRtSkNvWWY5eHNsTUJhU1JmKzVaUkZJQlVWQTQvcW9HQ2lER0c1Y1VGK0V5aTNGYXpCay9xb0pFeWJrT1N2YmJydHQ4d2hETUdpUGM2bFoyZWFiYjU2ZlF4MHpEUmxEeHBBeFpBeVpidXAzR1RWcVZENmFxSUNsSXBYZmI2L3pEaGd3SUQrdnpuUE9PZWNZTW9hTUlXUElHRExkVFJxbEtFZWxBcFZOTnRra2I0ODhjdVRJdkhObGU1NWIvcGlERHo0NFArL2RkOSs5d1NHeTM1NXpFckhudkNGanlCZ3lob3doVTBKNzdybG5YcFJTOTYzV3grUEhqKy93M0ppbXBxYjgvRys5OVZhWGc4ek9KeXhNeEdUSGdxT2JCU3VSckJnbFh6SjVNMnB5RlJuUnlGRmRLQUNKNC9kZnVEQlJ2Kytsb3RIbDlSSUNmRDZFRXFGRzZFVUZMUXZKbVJDZkI2KzMvNEtGaVdpNFdCQ1JScG9RWWVPdHdSQWh4T05GeDQ4Z1J1Z1FNbVB1VDhYMXZQN3FKbU5TQVNvRXd6L1ZWMTFIL3ovTks2VUNWSUpBZ3JEWUp3TUpDTFVhMTlBaklhTmFZYnJucmJmZXV0Mml4K3I1ZmxRUWM4MmFOZGtGRjF6UUthYkREQmxEeHBBeFpBeVpOaHpGZlBLVG44eHZYTGtxSFgxK1RjZHRzODAyK2ZtWExGbGl5Qmd5aG93aFk4aDBKOG52b3Z2dDNidjNCZ2swVVA1TkpVZG0zYnAxaG93aFk4Z1lNb1pNZDVGYUhILys4NS9QcDZ0MjJtbW5Eai8vaUJFanNrMDMzVFIvNktvbTBGbVd2U2ZOVFVTb0ZBcGd3cEhQcG1FRkFSSjBiSE05OTZkUkxoaEZOc1VDTkFwTnR3aVZrb0VCVENaa3dVVW1yTkZJMDVFZkdYMUNLMHFRaTFSdzNyS2dJd3dUa3hYcFdLZGpQSFRVMDdsY3Nxa1lqMTlvN0lYckl4UWlpQkVxWTM5MmFWMHhrSUNCRHFIenZnNkVpc1V6b1Z2TEtTeVllZFA4dWdxVFQzc3laSlNQOHVsUGZ6cTc4TUlMYzBQZjBiNFlUZE50dlBIRzJXYys4NW44b1QvMzNIT0dqQ0ZqeUJneWhreFh5WGVaT0hGaUNCa1plZmxoTklybzZHdlVPU3VsYW80NDRvaE9sUWRqeUJneWhvd2hZOGpVa0thOFZQcEZVMkRLYzJsTktaYytmZnJrbWYwZEJULzlPM2JzMk9aaW1OS0diclc4ejc3WEpqSmtEQmxEeHBBeFpHcE1QMzN1YzUvTHBreVprbjM5NjEvUHk3UzBwdWlsZ0NUZlRFZGNvNmJsTkhvU1lJNC8vdmo4WWQ5MDAwMGJmT1JDeU94KzZIV0p4azFibUtpUVBBaWp6d0tNTkpwTVZxU1JwaEV0clhsTnFhTDFjS1JIMjBmSm5qVDZoVVpRZ0JDaHcrTVJpb1ZreDhnWkd5VUhCazI1UXNNVUpPejF1MzVoSXQ1ditQa0VueWNMYWhZZ2lnVEJBclJ3ZjMzdnVDWVJvVVVvRVVMakhwaVJhTWNIWmliaS9vVkFnWHFmRHgzNUJNZHRnU0xBOHdjR0ZieEwwUStHTHAyTXFRUksrVGUwekp3NU0vZTNLTUZTeVpVdDdhUHlNVi82MHBmeUFwZ2RjWTI2bHFPT09pcTc3YmJic3ZmZWV5OWJ2SGh4cDVnZU0yUU1HVVBHa0RGa0Fza0hvK2t2TGErKyttbytzaEZvNm8xU05LcFFDUE9Ydi96bERuSDhmK1VyWDhsT1BmWFVUbGVMekpBeFpBd1pROGFRcWVIZjZOdTNiM05kTVNVMXF2YVlGbVhPcTFTTE12bnJ3VVA3QzBMSzlHK3ZNdjRWNlRvRndkZGVlODJRTVdRTUdVUEdrT25NaFN3RkJrV0ZWWklvQnc4ZW5JOWFGaXhZa0k5aWxNR3YwY25PTys4Y0hrdFRiUFBuejgrbnp0clRONk1xenJxbXpyanNjY2gxaVZqd2tzbVNUQTRzTkwyaW94cEdna2E0dEdDMG1QeFlFQjNqTVBLRm9uNzhZc0ZvTS9teWJGT3VLRENBRUNha0NjV3cwUlNlUndIcWM1dEthZUMxMERVbGhmMGppQkJLaEVyNCtUZjZmdkI5Z09HbVk1L0puVHMvT0NQUmhJY3VUa1FJRlFwdzFtdVFWaklaTTJ4cVJoRXFVVU8yQ0ZJM2Q0R21aVU9IRHMwdlVsRmhsZWdzalVoKzlLTWZaZC8rOXJlelk0NDVKZzhQcnF5VFg2YmU4WVlQSDU0Ny94OTk5TkY4OUtQL2I0OVJsMFpMQnh4d2dDRmp5Qmd5aG93aDA1a2gwNjlmdjl4Wlh5bG91WExseXV5ZGQ5N0p3NUFGRkUxN3ljOGlDS2xQUzJ2NndLaXNpMFkxOTkxM1h4NzVKZEMwaFk5R2ZpS050SFM5dXJhMWE5Y2FNb2FNSVdQSUdES2RGVElheFdoYVRORlpDbEdlUFh0MmJpUTFjcEZUZlgwN1Z3b0dtalpUT1BIU3BVdnpxYk10dHRnaW5HNkxrajExVEZWWnZ1R0dHN0tYWG5xcDAwQ0Z0Y2dNR1VQR2tERmtlaVJrbExSWWJlaGwvQmQrL0tXdVhoWXRXcFJIaVNtYXJKRnpxVE9sanZQUlJ4L2x5WkhUcGszTFJ6WHJHOTZzRVl3YWtYWEdoWkRaOWRnRmllallaeE1zTnEyaVk1b0ZKQXRHQWFMUkt4aEJHQ0VhS1JhSUxEUVZJK1FhZFJUalM4UHpGUXB1SXJtVDF4ODV4bms4R20wNjJtbWt1WDNoZVFSR3R3QjFmbjVCZ2REQzg2VWhnOUh1OTZNNWlkaTBqTW1TTk9wY3orUkxRb0dPOXpEUWdjNzIyNjlORkVGbjR2KzVLTkdVWDM0NzBkNlBuSitvZXQrbytHYkI4UjRwZ0VycC9SbDAwdGtoVTZuclZjbGxVUmEvZkMrVlJVQlFSRmxiWmV3cllLQjZTdXZKSjUvTUV6dzFnaXA3TEYyWDhuVU1HVVBHa0RGa0RKbE9DaGxWS1phUmw3R1dmMFhKakI5ODhFR3pzWnczYjE2YkZyZVU0MTlncVY0VUVDRFliYnZ0dHJtdlIwRUdtZ3FMcHZVVVRQREdHMjhZTW9hTUlXUElHREtkMlNlamhsNHk2clVXK1dMV1o1VFJrdVNjdi9QT093dm5FV2lPUHZyb2JQcjA2Ym5mUjM0aEJSM1VLc0twdjJuOWd3OCsyR2tkKzRhTUlXUElHREtHVEZVU28wWXorKzY3Yng1RlZyMG9BYk10UzhKb3BLVFJVYlM4OWRaYnVVTmZJeCtPYWpUYVVRNVBaNDRlRzMvTWdrUXNjRG44RWdpT2ZUYjlvdU83a0d3SkkxUkl4Z3VTNzBLb3dOSE84OVBvRlJ6OS9HSXNxaTg2VTJtVW80S2RoSERoZm5IOUVjUUptUWpTaFdSUjNCOGQzd1hJUUh4K0xGakpwbCtUSGpzbjBiRlBuWnpvck44ZG0raWlaNDlJZE5tL0haYm9pdWUvV2xjem4vdGFvZ3VlUFNyUktiODVJZEVSdno0OUVZMCtreTFaRUpOUTRmM3Y5ZkFGaVhpKzQ1K2FudWl3WDMycldUdzNnUndtUnphcUNDcjh2blNWWkV3WmRQazRWR1c1ZXBrOGVYSStmYVd3WlkxNEdqMlBwc1NVYTlQYVJSRnVPK3l3US9QK3UreXlTejZLVVZpMUlXUElHREtHakNIVHlTR2ppREdOQ3ZUZk11Yk1tRC90dE5PYVMrWXI4cXpSOHdsWStyYzFpM3hES2hOVEtXc2o2ZnFVRE5yWjgyQU1HVVBHa0RGa0RKbVBwYW13U25hL1JqSWFLVlF2eTVZdHk2ZW5OQUpSbmJLMjZwYTVldlhxMEhETFQ2TnB2T3F5TVFxQjVwU2VJV1BJR0RLR2pDSFRTU0VqSjdxdXBSS21YR251OWNBREQyU25uMzU2ZHY3NTUyZjMzbnR2dHRGR0d5VWppa2FuekNaTm1oUWE3alBQUEROMy9sZkt4aWlhN1BMTEwrK1UwV1NFeXRoVG14SkZ5WlZNQnFRak96S2FJV1NpWkw4SU1ndzBvQkZrOHR5ZGMrcXE0T2hsQWNvZ1dYUElGVTJKQ28zT2dnS2JoRFlETGZqOG95WmtoZVJDRnJBTUhQMjhQeGFNcEJHbTBhYlJiM3J4Z0VUWHY3UmZvdmt2VGtwMDlRc0hKU0prQ0JIcXF1Y1BTVFR2OTFNUzhmelI5ZkQ0Wi96dXVFU0U1dUcvUGlNUjEwLy96WWwxMVJKd0pDWnlzcUZjYWVpVVRNWU1remtSRkNGMU9wK01Sakdjd3FxTWJGUU9Scmt0R3MyMFpWK2ExdmhWVkpSVGtXM3lFeWtxVGRla3ZCMUR4cEF4WkF3WlE2WUxRVVlaLzhwVCtmRERENXVOcGtZdWdzSG8wYVB6b0lCR3MvMHBRWXRUYzlXTFlLSnJVRFNhZkVKbm4zMTJjbjJHakNGanlCZ3loa3dYZzB3bENWUC82dnJrLzFBRkFFbkFhZXV3YVlVbnEyUk5TOHNKSjV5UWgxVzN4bit6b2FGaXlCZ3lob3doMDJNZ295bXVNbERRTkZSMWlYeU5HRFJ5VWV0aWpTamtySzlFaGJXbFZNdE1qdjAzMzN3ejYyb0xvYkxqS1UySjJHU01qbVVhdlVKQncwWUxHR0o5SVhrd0tNQklxQlNTQWZGRlk1T29ZVCs5b3E3b3VDNDROd2tkZkxFSzBBQVVDQjJ1SjZRaWFQRjU4Zk1vZkQ2QUNBMFJuOGZCajUrVjZOei9kM1FpUXVDYTN4K1lpT3N2ZWU3cmlVNzc3ZkdKYUlUcEdLZlJqaFE1MWdtSkNJbzMvMkhQUk45L2FaOUV2SDlDaVlFTnZENG1aKzc3NkhuTll0QUVBYy9QamdZK2hBU2h3NlpuVEZ6RjhkbndyZGIzcTkwaG8yUExuMUVHTW9jY2NraUxCdlhrazAvT050dHNzM2E1VmdVQnRGUnR3SkF4WkF3WlE4YVE2V1NRVVJTV2psMG1TMS85WTZaTW1aSVlVWlhOVitTWjhsSTAybWlyQXBtMW90czBiWGJ1dWVkMldxZStJV1BJR0RLR2pDRlRsYm12WTZ2WVpSa24vTlNwVXhNamV2amhoK2NSWFNvcjA5Wk9mMHJSWTByMFZEVGJQZmZjMHltaE11R0lCWWtNR1VQR2tERmtlaVJrTkUxV2RycE1aZjdwZ0ZkRE1SMUhFVjZWRU9mMmxFWmdsY1pvWFFFeU81MjBNTkhJODVzU1JVYXFiQUhLZ3BGRE1tQUVIVUtDZ1FOTXRtU2dBQjMxTEpCSXlEQlpidlMvWHBhSWptMDYrZ3VPZGpyR2NYMTA1QmVjL1lRT29GSW91SW4xVWNGTUZtemM0K0VMRTAxN2VscWl5SEVlSlVOeWV4cHhRb0JHbFZDakVUN20vNTZTNktnblQwdEU1emlQUjBQTjgvTjRoQUx2ajVCaDRBQ2ZEKytIVGN6Ry8veVNSTHV1K0U2ejJBQ054VGJIL3V6U1JBeHFxVzZBbG90UXdYcis0SXFDWnZqZFkvS28xRzZRVVpUWVh3OWVzM3B4U3hKSXFydEtLbjlGeDlCVTJRc3Z2SkNYK3BlVHZyMUJNMjdjdUR6aDBwQXhaQXdaUThhUTZXU1FrYzlFU1pNNnJwenBaZmJWVkZYMU5OWFRUeitkanlnZWYvengvUCsvOTczdjVTSE9na0I3UWtabGEzUWZob3doWThnWU1vWk1KNExNb0VHRDhrS1N2WHIxeXN1d2xQV2hDQ2hYWEhGRlhTT3JjaTVLaXF5dWlOeldFdXh1dlBIR3psR0w3S3ZYSlRKa0RCbER4cERwa1pDUjcwWFJYOHFOMFFpbU5SMHNCU0VGQ0ZUK1gyWDhXNnBzdkdyVnF1ekNDeS9NYnIzMTF1eTY2NjdMeThIb1BDcFcyZFlqTVNWL2RwYXNma0ptM0xTRmlkaDByTkEwQzBhcGtQeUhBcGdNRE9EMkJjRUlGZ3BZQm81OUdtRkNyWkJNaUNLQkxGaklGNTlRS1h3eDJkU0tCU0NESm1kY1g0QUduaWVUWFFtVnd2MERlblNzL3QwdnowNUVJMC9IZXlSQ2hZNXRKamN5R1pQWHMvOHZ6azMwMVNmT1RCUWxLNTd3LzUyVWlORGk4YWlvU1Jpdmo4Y25kQmpZd09STlFweWdvT281OXdrZ0FvckptWU4vZkZXaUNCSUQ3NXFkYU5EZHMxS2hnUnkvYTd3WFhxL1VacERScUVMaHZ5cTVvbi9WVWJKTWNFQUZTQm9KeVovRFpkYXNXYm0vUmlNZGpUSXEwM0dTcHMvYXNnVzBqcW1SMHFtbm50b3BRcGtOR1VQR2tERmtlanhrTkFJNTZLQ0Rjc2U4UW8xYnU5L1lzV1B6aTZqNFdPVFUxK2lFaTR4L2RiNk5Tc0ZvQkNTMUZtaXRrV3FZNlZ3YUhXbFV0bVRKRWtQR2tERmtEQmxEWmtORFJxTUxqV2FVWDFLMkNyTE9yeVRNaXBIWHlLUlM1cit5YVBwS0k1eU9yS09tL0I0VjVlem9aYmVqRmlReVpBd1pROGFRNmRHUWtWOUYwMWNDVEpsUlRIVlZBRWxkTWZYL3lvczU4Y1FUazR4L1RaVzF4c2ZUMXNta3VxK09uaklqWkFpVk1XYzJKUnArY1NvYUtSb3hOaWtMa3pNQmhVS2dBS0FTT2JZTHlZWTRYNkhKRmhMQWFHU2pwbEpEUHY2eVZZdUJBSVJXd2ZIUGdwTkJjaVlESWFLQ21ReVVZRUljNzQ5R09XcnlSWWpRMFI4WlRUcTZhWFJwSkt1VENXc2xGQkpDTk9xRUN1KzNYdktpUkdod1BadUlFVG9NSk9EOUVycHNpa1lROEVkUHRXTmY0djNVQXlJZC80UU0zMjBHd1VTSnl2eEJGZ1VlTUtpRXoxWnFFOGhVY2txVVNLbU1lYlpUcm5iczE1S2kwQVNSU3JpelJpeHk2bGZ5WmQ1OTk5M2NSNkpSVGtlM2cxYVUzSVlleVJneWhvd2hZOGowYU1pb21aZjhKYjE3OTg1OUtocHg2RzhDaCtBZ1lOUWJoY2ovVVhIa3kvR3ZiWlhaWHoxVnBjck1PbjVIQVVaQTB6M3B1Z3daUThhUU1XUU1tUTBJR1RuOUR6end3THlDc1lBaS80eHFqYWx0c2hZNTBDditGRG5WYStYT2FKcE4xeUdqcmlyTE12S0MxSXN2dnBnZjQrV1hYODZQM1JHakdWVW9FQndya1d2dDd0Zy81THBFaG93aFk4Z1lNb1lNd3BkbGpKOTQ0b25zK3V1dnowY2NBb1VjNTBxY1ZOU1p5c0pVU3YrM05DSVJnR1RjTlFxcWpJRE9QLy84Wm1OOHhCRkg1UDZhbHNyVUNISTZ0cWJvR2gyWmFVU2xhMVhrV2tkRFp1enBUWWxHblF1ZGw0b1FvUkZqa3pLS2p2OG9jS0JRaEJGR3RkQ1VEUHN6K1pLT2N4YmQ0NHZQTHg1ZmRINHg2THdzSkt5eDBSS0xDTjQwdjY0S2phS3dmeUZ3QUpDS0NsYlNDRWRRNFhwQ2hjbURpMS9lUGRHU2x5Y200dkZwQ0NPalRVVVFZdklpb1VNSVJPY24xTGcvejBmSFBaTXJHZmpBNU04b09Yak0vWmNsNHZ0YUx4bVRQemdZeEVMb1JJNzVDQnFSYWtHbFhTQlR5Vy9SeUVNNU1xcWtyS2t1SlRVS1BzcWMxOGltWXJnbEFhSFdjWVlQSDU0RDVweHp6c21hbXBxeWQ5NTVwOWtZNjcrMW45WXJ2MGJUYWdwbEZsamtDOUs1ZGV4R08yY0tqZ0xaVTA4OTFTRk9mMFBHa0RGa0RCbERwcFdSWXZTOWFCcE00Y2YzMzM5L0Rod0ZCNmgwLzVaYmJ0bmljUVFhYmJ0aXhZcWFSbm41OHVYNXlFaWpHZ1VMVktiWk5LSVM2QnE5RC9tREJNUWJicmloWS9KZ0RCbER4cEF4WkF5WnhqcE9LdW55d1FjZnpPNjQ0NDQ4V2t6UmFCcng2TysxZkRTVlVaQ3FNTmRhMXExYmx6YzMwekZ1di8zMjdPS0xMODVISFlKTUZNMFc1ZTNvMmlvaklrUEdrREZrREJsRHBwTkRSbE5hY3Y3UG5qMjcyYkErOHNnaitick5OOTg4TDZpcEhCdU5IcXF6OXpVTnhpNlo5UllCU1pEUitkYm5PdnYyN1p2N2dUVEswak41NXBsbk5rZzBHWnVNamJpd3ZncE5zSkFzV2RpSGtFRnlZS0ZKR0tCQ28xa29nTW1tV1ZoUEkwM0hPNytFZkxGcEZGamdrRWFNamxaQ3A5Q1lDVVVEV1pDejBBanFId1BkRWdpUTQvVkdUYmFpNUVvNjd1LzY0NDUxeFNaZGZONDBJbVVkNzVFSU1SYklaTEltcjQvbkoxU1lQRXJIUGFIS0ptYmZmdWFZUkRUMGRMNFRETFVTRmxzeS9KR1JKeVFZZE1CM0tWSUVvZlZSaDBDbVVqNUd2aFE1OGxYb3NucDU1WlZYc3JsejUrYVEwVGJ5aVFnMmxXaTExMTU3clZYR1dsTjEyMnl6elhvblhncFF6ei8vZkxQZmFFT0ZMQnN5aG93aFk4Z1lNcTJJMEpJdlJnRUExZjRaRmJpc0ZMbzg5dGhqbS92RVZQdGJCQnRObFNsY1dZbWFMVlZscmw0MERTZS96UHBHbG1rVW8rQ0JtMjY2S1g4b2dvNGhZOGdZTW9hTUlkTkpJU1BuZldWRUlMOUo5VFNZL2x0Sm01b08wOGhGUGhBQjUrMjMzMjQydkdlY2NVWU9LWUZwazAwMkNRMzEyclZyOHhCb1RiMlZhZmRjbWM3VGlPblJSeC9OcjFYWHJIRHM5bHIybWp3M2tTRmp5Qmd5aG93aHM1NVRaQ3JMb3VNTEZMVVNLVFhGcGUwMHNsRUk4c3laTTV1TjhhR0hIcHFQZWhUKzNOcmxWNy82VlQ0S0VuQmFHeEduYyt2YW5uMzIyZnhhMVJXelBSZENoc21WZE1RUC8wNHFRb0dPZmtLbXNQK01WR0Z5SUpNUGthekk1TWtDVkZoUUVvMlI2QmlOaXZEUmlOQlJlL3hUMHhQUjhjc3ZBUjN0dkw0K3Q4OU5WSEQwNDNrVUlCS3M1L2xvZE9uSUoyUzRuazIybHEwY2wyakY2bUdKZnZ6SEhSTHhlQkZFMmxvTUJHREJ5OU4rZTN5aWMvL2YwWW40K1VkTnlaaU1Tc2MrUldpVmhVcmt6SzgyK2xHUUJOY1RRb1JPMldUS0NHcFJBVStwM1gweXl2UVhQT1Q0VnlXQWVwMHlkOXh4eDN4a282bXZ5cUpnZ2VPT082NlVFZGYyQWtkcnJrK2xZNVE0cWh3Y1hhdWV4Y01QUDJ6SUdES0dqQ0ZqeUhRRnlGVENnaldTMFhuazFLKzNyYWJIQktRMWE5YXNsd0ZYTjB1ZHF6WGRNZ1U4alhydXZmZmU3SmUvL0dWK2ZTcVIwOTZMSVdQSUdES0dqQ0hURG1Yekt6NmFxQytNREwybXlPYlBuMTg2NDE3QWtDK250VzJXRlkybWMyaVVwV3Q3N3JubjJqMFB4cEF4WkF3WlE4YVFhUWRWV2h2TDBWK3ZtNlg4Sk5wV1hUYVZzMUlwdE5tYVJYMW9LbjRjblVlSm9MVkdOVHFId0NLUVhYWFZWZmwxVFo0OHVVT1NMVm53a2dVdVdiQ3lrR3daRkZ5TW1wSVZvSEw5d2tRMGdnVkhOd3RJTXRrUVRjYTRQNDBxR3lleDhSSkZDREg1anNsN2RKelRrY3o5Q1oxK1A1eVRpSUVCQlVjL2t6VVpPQUJJOGY1cFZBbVZLUG1TUnZGbnEwY2tJbVJ1Ky9mZEV0R0lSbzc2MWlUazFRc1VvT0drbU96SUh4VU1kT0R6SUpRSUdUNWZRcHJyZVQ5TUZxWXpuWUVCZkg4WmlGTDlyak81c3RIaW54RkV5Z1lhUklFSEhRNFpTZVZnZEQ2Vmw0bVNKZ1VDUllySlQzUDMzWGUzeXFBcjUwWStuYnZ1dWl2ZlI4VTVGZHBNWDVCSzBBaENmL25MWDVwSFdJcFFNMlFNR1VQR2tERmt1akJrVkhoU1BoZUZHaXZidjk2SXB0cFBJOUFJSEdVWGhVSnZ0ZFZXTldGMyt1bW5aL3ZzczAvK0VNNCsrK3dPS3h0anlCZ3lob3doWThpMG85UkpVM2twbWc2VHNhOFhjVlpSdjM3OThnb0FaUmRGbWJFbHRNS1ZkZjVLeUxJS2JYN3d3UWR0QnBXSmg4OVBaTWdZTW9hTUlXUElkQ0JrS3Y0WlJYWnBOS09zL25xZE02WCsvZnVYTGxpcFV2MkNDU0dtaEV2VlJGTWttclpwNjRXUUtUUWRPNnNwRVIzeGhFclU5S3NBSllqSmxRVW9sRTBtUkNCQTJOU0xqbThrSHhhYWh0R0lGL1pQSGZGMHZQS0xSVWN3aytrSUpYNnhDTG13b09hdGN4TU51bnRXSWtLTGdROXMwa1hIZEFRWk5oMGpaSmF2R3BXSVRiZ2lxRVRPNVNqaGo5dno4eUprbUp6SjVFc21VeElTMGZPS2tsa3AzbjlVOEpJL1d2aWpndTlEOVRvVzIrU3pJblNpNHJHRVF2UURvQzIwd1NCVEdaM0lMeUxZeVBDM0JCcEJRcmsyWmFiTDNuLy8vYnhrdjhyRmNMcXVVcU9zdlJaRHhwQXhaQXdaUTZZVFFFWlFVWjJ5U2xVQVFhZldkcW9wSmhpMWRsRklzaHo3cXV6TVVZeEdSS3FuMXA2TElXUElHREtHakNIVENTQlQ3YU5SQkppdVE5Q3BOYlVtT0xSMlVjdG4rWHNZVkZBWkVTbVh4cEF4WkF3WlE4YVE2U0dRa1RTNjBIWEloNklLenRYcjVMeFhOODNXTEdxSTFsTEdmeVg1c3IyWGNkOWNtR2o0SlUySldMQ3lBQmxBZ283N1FrSEwyYW1pZ3BRRmlOQ0ljM3NtVXhJcWNQUVgxdk40VEViRTlrd0daYkluQXduNkw3czZFWjJsTkZxRUNvMFVqVHlQVjBqV0RKSXJDY0Vva0lHT2FrS0FUY2hvTkFtbGUxZU5TY1RrU3hyaHlIRlBTTkNvMHVqUzhWMHYrYkJXQWlJaFEvRkhCTytIZ1JCOG5tVUxqdkw2ZUwrOFAvNklZTk04cW5wYnZqdUVCaUVUZlhZUnBLSWdqRWkxUHA5T0F4bEpBRkR5cGNLVmxVeXB5TE5LVW1WMWhlWjZ5OVNwVTJ1MmR0YlVuSTQxYjk0OFE4YVFNV1FNR1VPbUowSkd5Wm55b3dnU3ltRlJJN1BiYnJzdFQ3QnNiU0ttQUZXcmJJMThOSnBDVTBTWklXUElHREtHakNIVEF5RlRHWEdvaUtZaXpqUk5ObVBHakh3S3JEWExrVWNlMldKblRGVmJibStIdnlGanlCZ3lob3doMDhraFUrMmtsOE5maVpJYW5jeVpNeWMwN3NxNWFhbjZjcVZKV2tjc0xFako1RWV1andwWU1ybXkvL3lGaVFxTzlrV0JDSlYvcUM5Q2dJRURoU1ptUWNGTU5qVXJORG5EL2dWSXdiRk95RVRKbW5RY00xbVQwT0VYaVk1OEpsdjJ2ZU9hUkV5Mm85R2hVV0pnUWxuSFA5Y3oyVEp5L05QdzBMQ3hLQ0t2dnl4a0tFS3NOWWFzV3Z3UlFTZ1RPZ3lVNFBQbS9peTR5dnZqanhJK0R3WUM4RWRIOVRvbWV0SnhIeVZqUnNWTnkwS0U5OTRhZFZySVZLUlJqY3JRUkQxZW5ubm1tWnFKbDVxQ1V5Q0JnZ0VVYldiSUdES0dqQ0ZqeVBRQXlPeTIyMjUxRXpDcjgxcVVQS21SaU1yQXRMU2NmUExKaGFreVRZOXBGS1Q4Ry9sM3lyWU5NR1FNR1VQR2tERmt1aWhrMUhKWjUyVkdma3RGTmVXZjJXKy8vVm8wN0NycFg5MXlXVDRZK1hXV0xGblM3bERaNmFTRmlRd1pROGFRTVdRTW1RME1HWTFnTkwzVlVnSm1yVXJNTFhXdFZGU1pxanBYa2k4cmJad2ZmL3p4RGhtNUVES0Q1alNsUWpMa2lJdFNFVElGUno2TkxpQlFNT3JCOWdYbzBIRWZRUVBxdjNCaG9nSHptbExOaGJDZWtDa0VEaUE1azhscy9KTHppM2I0cjg5SVJNYzZqVG9kdmR5ZVJwSEpwRXl1SXdTSGZHeElxa1hIT1p0d1JRVXllYjAwaW9RSW01ajljT1V1aVhpL2hEUWhROGM5bncrTmJpUWVqK2Rud2N3SVFnd01pSkkzQ1IwK1B4WllqU0FUUFEvdVh3MFZBaXFDVE9UWWp3cHNjdi9vZUFRS2k3dEtHOVFubzNJeGxWNHUxYU9RYXVkL1pUcE5JY2pLbTZtMVhIenh4Zm02NnVSTmhVQjMxR0xJR0RLR2pDRmp5SFJDeUZTeThDdWdFVWpZODBXakhmV0RrZFRBck5ZaTU3NENCQ29qSlBsd1huamhCVVBHa0RGa0RCbERwcnRBUnNaZFVWekRody9QSGU2YXNsTGZGam40bzM0eHFwWmNLU2xUblVpcDdwbHFqWHpoaFJmbTAyVXROUmFUWDZmU00wYm5WVnZsamx3TUdVUEdrREZrREprT2dFeGxSQ0xKNmE2d1lmbEdOTEpRSkZtdFRQekt2c3J5MTM3eXJRaFUrcHNnTTNQbXpORElxK3lNaW14V290RTBzdW5JWmRSNVRZbG9SRm53a283OWdkZW1pcElkYWRTcFFnRk1RaWFDQ2h6MVRJNWs0TUhBYTVvU0VhcUVKdGNYQWhrQVFTYXIwVEhNTHoyTk5JM0k5Uy90bDJqeHk3c25vbU9ZUnA2QkExRVRxcWdwRzQxcTVNZ25SQ0xJOFBycDJGN3k4c1JFMzM5cG4wU0VEZzBOcjUvUWpJeHNsSnpKSHcyRURBMWRCSjJqbmp3dEVRTS8rTHdKR2I1ZkVZUUpCbTdQOWZVQXpDQ01zc1VxeXpZeDQ3VkY1MjkzeDcraXdEUXFxU1JRVmxvc0R4dzRNRGY4QXM5bW0yMldUNUZ4ZEtQL0Y0aDBMWUtTUmpYYTlvRUhIZ2pMK3F1Q3M0NHZNS2wwektKRml3d1pROGFRTVdRTW1lNEdtVnJRVWFpeTRMSHh4aHRuTjk1NFkzYkZGVmZrK1N5Q2tLYTVOS1dtanBjVEprekk5NUhmUllEUmlPUzU1NTRMRGJ4cWtRbG1PcGVtNlZRRXM2UHlZUXdaUThhUU1XUU1tUTBJR2ZhRUVXaVV2YTlRWTBWL0NTeUNUV1ZxcmJLdHdDTXdhVVJ6N3JubmhrWmU3WnMxeGFiaW1xM1p2dEZsdDZNV0pESmtEQmxEeHBBeFpNN2I4SGt5bXM2U1gwWVorSExPYTVwTC94NTAwRUU1YUFTWDZ1MzEvL3I3bFZkZVdkZm9MMTY4T0IvOWFOdldGdE5zUzhpd29DV2hFVFlSaTZBQ0kweUk5YnQrWVNMdVgwaXVKRlFpaU9CK2VQMkRyMG9WTlZuai9vWHIvOEcxaWZpbGpRcEkvc01mOWs3RVpFUTZ1dS82NDQ2SjZCaW4wYVVScDFFamRKaDhTY2p3aXg4NTlnbWhDRHBSTWllM0o0U2JYandnRVkweWpYaFVRSk5RNGZxeWtLR2ptZGNUUWVld1gzMHIwYlNucHlXS25qZTM1L0Y0dlJRTmRmVzd6bWZEZDZXdElSTkJKMUtuYTFxbUVZMml3VFRkcGFtdVdtVm1sRXNqZ0R6NjZLTjFEZitrU1pPeWFkT21kY2owbUNGanlCZ3lob3doMHdVZ1U1RkdNMWRmZlhYdTZLOEZHdVhBS0VTNXN5eUdqQ0ZqeUJneWhrd1hnb3pBY3RSUlIyV3Z2dnBxbnEydndJQnEwQ2d3UUw2YmpuYm9HektHakNGanlCZ3kzUUF5R3FVbzlQaTN2LzF0OXRaYmIrVTVNd3FGcm9CRzZ6V2wxdG9Xek8yOTdIeml3a1JEcm14S1JFYzNqU3dkOVFWSE41SVRhZVFKaGJMSmxqVHFBNjZEa0R3WlFnWDN6NEtlM0w3ZzZFZXlKWTBRb1VLamQvTWY5a3hFUno0ZDI0VE9ELzU5MTBTRURxRkVpREdaazE5VUp0U3g2Q0VERitpb3A1R0xqRjZVak1uMVovM3UyRVFNbENnYkdFREhMdzBSSVJTSlJxdXRJUk1sRkRKUVlQcHZUa3hFNko3eW14TVNzV0FuUDI5Q3FUcW9Jbm9XalVJbFNyNHNDN1ZhQ1orZEFqS0NpWEpwNUxoZnVuUnB0bkxseWp6WHBqS3EwZCtYTFZ2V2FhYkxEQmxEeHBBeFpBeVpMZ1NaNnU2Vm1oWlRscjlHTGZ2c3MwOTI5TkZIWjJ2V3JNbU4rNXR2dnRrcHBzd01HVVBHa0RGa0RKa3VDSmxLN1RHTlhLb1hkY1ZVSGt5dlhyM3lQQnVWa0Rubm5IT3kyMisvUFZ1OWVyVWhZOGdZTW9hTUlXUElsT3N6YzhvcHArVGRMSlhFV2FrQW9OSXh5clBSMUpxbTBWVGJUTDRhUmFmSmw5TlJ5OWhUbXhKRmtHSFRzb0tqRzVBb0ZKU0VDZ1V3MllRTWtLR2p2MURRRXNjdkpGY1Ntb1FLSUJwQlpmdkYxeVdpRWVhWG1GQ2hrWXNndy9XUmFGUVpDRUFvOFhwb1pDSWp4K1MrcUNsWkJKV29TUmVQUjBnd1VJRDNUd2dUWWpTcU5PSTBYSVFRRFdIVU9JdlFpU0FUN1YrcnlHTzFDQVVXWUQzMnFaTVRzVUFuSWNQazN1cHJaZEJMVk9BeWdrelpKbVpSZzdqV3FOTkJSdEtvUlJCUkJuOEZMdldnSk4rTm90TU1HVVBHa0RGa0RCbERwbDFhT1N2SHhwQXhaQXdaUThhUU1XVGFwWEtBU3RTODg4NDdIUUtaMFdjM0pTSkVxRUovR1BwVTZKTUJCT2lUb2RFdTVOa1FNbGpQOHpQdkpvSUtJVkxJZTBIdHN6NjN6MDNFUElqSTUwSWpIdVhCRUJMY1Bzb0Q0ZmtJRlo2UDR2SFl2eVNDQ3ZOV3FHaDdyaWRrNklNaEZQbDhWcXdlVmxlRU5PK1BQZzBhb2Nqd1JVYU0wSWkyYnhRcVpmTmdDQ0VlbjlkZi9Tem9rNG1lWFZtZlNWbC9GZTh0ZWpaU3Q0Qk1wV3FBd3A4TkdVUEdrREZrREJsRHBzMmw1TTE5OTkzWGtERmtEQmxEeHBBeFpOcGVsWXJONzcvL3ZpRmp5Qmd5aG93aFk4aTBmWDZOcHN6bXpadlg3cEFaZVY1VG9xaDJGNkZUY09TemZ3dFVjTlF6cjJWZVUxMkZ0Y2lDUUFVR0FoUUNEMjZkbTJqWVQ2OUl4QmVUZVJpUkl6N3EvMElqeiszcG1LWmptOWRERVRvUjlIajloRTdaV21KY0grWFI4UHA1UGpxaUdXakIrL25aNmhHSmxxOGFsWWozeSt1aEl6eHkxRWVPKzhpUlh6WVBoa1l4Z2c2aHdmM0xRcENHdjE1UVJEMGcxWUpJRkFnUTFTS0w4bXI0V1hSYnlDaDM1cTgza28wWk04YVFNV1FNR1VQR2tERmsyazVxMlR4czJMRDhacFNjYWNnWU1vYU1JV1BJR0RKdFd2Tk1sWnZWZDBiMzh2VFRUeHN5aG93aFk4Z1lNb1pNMjJqbzBLRjVxWm5wMDZmbk4zVCsrZWUzTDJUT2Iwb1UxZTRxOUZTaEk1MkJBbkNrRTBxc2JVWkhQcE05bVl6SndBQW1SL2E1L2RwRTdGblBCREU2ZG1sa2FJU1l6RWRvRUJJMDZwR2puVWFWa0tFalBJSlF0UC84RnljbDR2VndmUVFaSmxNU0toRmtlSHdHSHRBSU1IbVMxOHRBQ241ZWZENE1MR0R5WWRua1NScktzcEFwQzRrb01DQkt4dVQrOVNCU3k5QlhPL3BwNU10Q2h1THhHRmpBN3phdmpkZk8vV3VCcWx0QVJpVm9Uai85OU95blAvMXBma003N0xDRElXUElHREtHakNGanlMU05WSmw1eFlvVjJicDE2L0liMm1pampRd1pROGFRTVdRTUdVT21iYVJDbXBkZWVta09BSTFxMnRzdlk4Z1lNb2FNSVdQSTlDRElqQjQ5T205MmR2Lzk5MmNISDN4d3UvdGw2UGduUklaK04xVUJRaXlvT1NkVm9kOE1DbWh1djJSZW9qNC91RFpSdngvTlNUVG83bG1KUnQ1M2VTSW1SL0pMU3Njd2pSaU4wSS8vdUVPaWgxWVBUVVRJRUNJMDBsSEJTeG81UWlwS3ZxUlJKVVRLUWlrS0pDQlVJc2MrajA4SVJaQ2hvNStmYi9RakllcEh3MzRxZFBSSGp2SElNSlpOenFRaUtCQWlVY0loQzFwUzlmckRTRFRDTk96VlJwdEE1YlUxQ2huYWd2NTN6a25FSUI1ZU82OXY0dis1cUtCdWs0d3A1Ny82eTl4d3d3M3Q3cGN4WkF3WlE4YVFNV1I2R0dUVVYyYmpqVGR1aHN4blAvdFpROGFRTVdRTUdVUEdrRmsvb0VqOHUvSmtOSFdtRVUxNyttVU1HVVBHa0RGa0RKbHVESm5ldlh0blgvemlGL042WmRWL256aHhZaDVacHVUTTl2VEwwUEUvOUxKVWhBNGQrWVhrUmpqcWk0NzR0TUJrdngvT1NVVEgvTmlmWFpxSUw5cTBwNmNsb3RHZ1VXYlRLa0xqMTMvcWwranBWL29rZW16Tm9FUXNPQmxCSW5MMEV4SVJaS0prU3E1dkZGSmxBd3NvUW9vRktBa0ZRb3BHbGU4RERSY0RBWGo4Yno5elRDSkNqT2Vqc3pwS3JveUtPRWFPZnQ0UElSQTU2cU1DbDRReXhlM0xGclVzYzY5bEMyQVNFZ1ArWlhhaVB2ODhMeEdoTStUSFZ5VnFUVkhPRFFvWitWRlVQYm5zZnR0dnYzMjIzWGJiNVZXWHYvemxMK2ZIbVRCaFFyNXVpeTIyeUE0NDRJRDh4Z1FiUThhUU1XUU1HVU9taDBKRzU5MTg4ODFiWFFTek1uSlI3NWhERHowMGUrT05ON0x6empzdkd6Um9VRDZDcWNCSGZXWDAvenIreXBVckRSbER4cEF4WkF5Wm5ncVpUVGZkdEZYYjl1M2JOOS8rSzEvNVN1NXorZGEzdnBVWS9wZGZmam4zeC9UcTFTc2Z6Ung3N0xINTlvY2Njb2doWThnWU1vYU1JZE5USVNQVmN1Slhlc1I4L3ZPZno2ZkZOdHRzczN6YlBmYllJMCs0dlB2dXV3dkcvNk9QUHNwMjJtbW5IRFJQUGZWVXZyMysrei8vOHovYkZqSVhOQ1VhZG1tcVFvRkpGS1FzTkMyN2VYNmk3ZjlwWHFLK1M2OUpOT2IreXhMeFM4RW1WSFFzMDRnU0lpeUlTSWo4N3BYdEVqMzdTdTlFaEF3REFXaTBvNEtVaEF5Tk9JMTIxUFNMUmpnS05JaWFvRVhRaVVTbzhQcVlUTWxrUjk0ZmZ6VFF5RVpGRStrSUozVDRmdEhSei8walNKUnRPbGJXMGM5azBBaUNaUjMvaEV3VUNCQkJyUjZRb3dLWVVaQUV0NDhjLzRRS05lNkJHWWxxQlJ0c1VNaW9jckxPclJGSXJmVUN5d2tubkpCSGlWMTU1WlU1UUpZdVhWb1hBTys5OTE3dWkzbjIyV2R6WDQyT2Y5SkpKeGt5aG93aFk4Z1lNajBOTXZLaEtFcE1PUzBqUjQ1TTFnMFpNaVNmU2xNOXN2VmQ3cnZ2dnZ3R1ZRMmdrZU1ZTW9hTUlXUElHREpkRURJQ2llQnl5eTIzNU5OaUd0bklpYi8xMWx2blJTOGZmdmpoaG9GUUdjM01tREhEa0RGa0RCbER4cERwU1pCUlIwdjVYYlM4Kys2NzJaRkhIcG1OSFRzMk8rV1VVN0lYWDN5eFlOeFZhWG55NU1uWldXZWRsZDExMTEzWjIyKy9IUUxoa1VjZXlXOVMxUURrczJtTFpjU0ZUWW1ZZk1rbVlHd3kxdmY3Q3hLeFlDVUxWSTVhL3QxRWZQR1lMRW1qVFNPNVl2V3dSRS85cVcraWYvdnp0b24rOE9ldDY0cmJFMUkweXBFalBITDhjM3M2b3ZtbDU1ZWNScEZHbkpENzRjcGRFaEhLZkw2RUlEOFBRb1JRNFBVdzJaSDNTOGpRVVI5QmhxTFJqWXdxb1JJNTRpUEhlMVF3azlmSDR4RXEvTHdqeUVSTnh4cE56b3pPRndWR2xJRk9CR1FtVHRMV1JDS2srSU5GMnFDUTZkZXZYN2JWVmx0bDU1NTdibmJycmJkbTgrZlB6NlpObTVidHYvLysyWUVISHBoOTdXdGZ5NDQrK3VoczExMTN6VWM1bi83MHAvTXBOam4rdi9DRkwyU2YvT1FucytYTGx6Y2IvM2ZlZVNmZlIrSE9EejMwVU82ZitlMXZmNXVOR1RNbXY5Ry8vL3UvTjJRTUdVUEdrREZrZWdwa2xJaXAwR1JCUTJISkFvbnlYQVFmL1YzNU1JS0t1bDZxNEtVaXk2cjNWemp6MVZkZm5Sdit4WXNYNXo0Y0pXZnFHQXBqWHJac1dSNWR0dnZ1dStjM3FrQUNROGFRTVdRTUdVT21oMENtTGNyTEtGOUcwV01LSGhneFlrUnpiVE1GRk15ZVBUc0hVU1cvUnYvV0NuMDJaQXdaUThhUU1XUU1tWnFRbVQ1OWVqN1NHVDU4ZUdFcWJyLzk5bXV1eWp4czJMRDhYMDNQTmJvTXY3Z3BFUXRlRnBJdHY3Y3dFWnVJTWZseTRGMnpFN0dBSlY5YUdpa2FRVUprMWF0YkpmcnphNmxlZjIzclJGeFB5RERaa283enNvNXZPdElqeURBWmtGOGNCa3JzL09DTVJEUUtoQVFodzBDR3FBbGJGSmpBNUVvNjlpUElSTW1TVVZIRnFBaGpXVWQzbEx3WUpTdEdUY2Q0UEg3K0ZDSEQ2eTNiWkl6cm8wQ0dzc21tdFF4MVMrZU9JQk1CT2tyV1pQRk9mbmVpQXBwZHZtbVpPbUl1WExnd254cFR3RUQxT3VYS2FOcE16bjZWcnJub29vdXlLVk9tNURmODRZY2ZHaktHakNGanlCZ3loa3g5S1dMczBVY2Z6YWZLeG84ZlgvRDNhRDByQXFneTgzLy85MzhiTW9hTUlXUElHREtHVEgwcHV1eTExMTdMbmZ0NzdybG5zazVCQXZvN3c1Wi8vdk9mMS94N3ZXWFlqS1pFaG93aFk4Z1lNb1pNTjRlTW9DSllLS3RmR2YxY1AyN2N1RHphak1zOTk5eVQzL1NTSlV2V0d6SVVDMkQybjc4d0VRdGdGaHovdDg1Tk5QUW5WeWJpaDNiYWI0OVBSQ05JUnp5aDhWK3ZiNXRvN2V2YjFCV2h4R1JMUWkyQ0NpRlNGakpjVDBmdjZIKzlMQkdidURIaGpGOWtPdXA1ZjRRTXhRS2doQmFoRUVHRkVHSmdBSTlIb3hvWm5yTEpqalRLaEFCL0JFVk55NmpJa1I5QmhlTCtFZlI0ZjN4ZVVlQUUzNmRhaHJkYTlmYm50bXc2Vmc5SXRhQVVCVFZFOTFyZFVFMGlaQWdscWN0Q1JnM0tLaVZqbE5YZlVqVUJMczgvLzN5KzN6SEhIR1BJR0RLR2pDRmp5Qmd5TGJkZ2xyTy9wYVpuMjI2N2JWNWNrNHVjL3JwZkZkczBaQXdaUThhUU1XUU1tZEx3R1RCZ1FGNzdUTldiYXkyNlgxVU1NR1FNR1VQR2tERmtESmxTZ0ZIWXNoSXhWZXEvcGVWVG4vcFVmdU90aGd3S1lETDVrc21XTElCSnlKUk52dVNYNE9ZLzdKbUlVUG5MeHlDcEZxRVNRWVpRYWpUWmtrWTJLbkFaSldOU2RIVHYrTURNUkgxK2NHMHFOSUZqb0FBaFE4ZCs1UGlQbXJJUkdoRlVHTmhCOFVkSDVHaU9JRkcyeVZoa3VLSUNsSFRFRTVLRUJndDJVcEhqdjJ6Z1F0bm5RME5Qdzh4QUhocmxDQnoxVktZaFdtdCtZUERkNFhlRnlaZ01ESkE2RFdRRUNKVi9XZC85RlUybWtqSXFVYU9JczNwTHBZK05JV1BJR0RLR2pDSFRoU0VqY0tpRWYvLysvVU9BeVBBcldrd2prWlk2WmRZTEFsRHRNem42VlJReldwUS9ZOGdZTW9hTUlXUElkSU9SalByRDZQaXFvQ3hmQ1l0Y1ZrWWhtc0xTRkpmQ2pwbTlINDFnQkJqOTI5cE0va3JiNTFiN1lBd1pROGFRTVdRTW1jNEpHWTFLVkM5TWlaTTZqMFlyR3JWb2RDT1lLUHhZeFNzMTFhVkZVVjhDVTJ1UHJYMVZsNnkxZ0tsRWw5VnovQk15UXk5THhhWmsvUmN1VEZSdzlDUDVrcERwOThNNWlmaFNzdWtZa3lzSkZVS0Q2eWttVzY1NU5kWGphd1lrb3BHUGtpc2pxTkNSejBDRzVhdEdKYnAzMVpoRWhCeU5CcDJUZys2WmxZaEdoSkFqVkhoOWhBcVRNUWt0T3ZvSlNUNGZpdnZUNkVVRkwybGtvMlRGeUpGUGxYWGNsNFhHdEtlbkpXS0JXSzR2R3dqQSs0K2FzUEY1UndtTmJQUTE5bWVYSnFyZXRteVFRYlMrYkhJbnhYdmpEN29ONXBNUkREU1ZwVnBqZ295bTBGVEFVaU1RNWJPb1NkbHp6ejJYRy9pWk0yYzIvMTFGTDlWTHBsNzdablhSVkIrWjFpNVBQUEZFR01Kc3lCZ3lob3doWThoMEljaXdwcGpnb09rejlZeFJnY3ZxVWNnSEgzeVE5NUhSQ0VmVmxaVnNxV2d4VGJVSkRKb1dtekJoUXU2N0VheldyRmxUcXU3WVpaZGRsdC8wMldlZmJjZ1lNb2FNSVdQSWREZklWUHRTTkZyUnRKVWM4ZC80eGpleTFhdFhGd3ordW5YcjhvNldBb3hnbzJtM2pUYmFLQzhiVXhuOWxGblVkVlAzKytDRER4b3lob3doWThnWU10MFZNdFZUYVR2dXVHTWVlcXhSaTFvdzExdmVmZmZkN0kwMzNsanZDc3JxanFuNzFZaXBSY2pNYkVwVUtJQjViU3BDcGQ4TkN4SVZJSFB6L0VTRGYzeFZJbjdKNkdpblk1N1FlTDFCc1FBbWpUd0REd2dKUW9YUUlaU1k3TWlDbTcvK1U3OUVYTS9yNC9YUTZOREkwSkVlRmNRa1JOamtqT2RuUVV3NitybWV6NnNzVkdnRWFmVDVma1ZHdXIzRjgwZlh3L1duL09hRVJGeGZGaklSZE9nYzUvT1BFaUFwR3Y1cUE4MGZuR1doVXpaNXMrenhlWDJkT2s5R1VpQ0FRTE5vMGFLc3JSZFZYYjdsbGx0YWxlMXZ5Qmd5aG93aFk4aDBROGhvQ20yVFRUYkp6am5ubkRhRHkvdnZ2NS9ObURFais4eG5QdE1jdXF4dW1ZYU1JV1BJR0RLR1RBK0RqSXBhamhrenBsU3ZsNWFXeW5UWXBFbVRtdUdpNElGWnMyWVZtNUI5cHltUklXUElHREtHakNIVEJTQlRKbk5mam4yTk50NTg4ODMxQnNzenp6eVRSNlJWSjF1dVhiczJXN0ZpUmQwOEdrSW1jdlF6K1pJUW9kaTBiUHNsOHhMUm1VWWo5dSt2YnAySVVHRUJTenJ1STZqdytJK3NHWnlJam0wNjJxT0NtSVRPa3BjbkpscXhlbGdpQmhvdzBJRWloSGc4bm8vWHovdGpZRUZVOFBLSzU3OWFWM1RzUjVDSm1vNUZ6bGxDaUVZMU10Sk03aXdySGk5UzJmMmo4L1ArR05qQTUxRVdPbEZCVFFaQ2xHMzZWdjNaTlZKY3MxYUJUU1pUUnNlUElNWHRPeHd5U3NKVXVISnJZQ1BILzlTcFU5Y0xMaSsrK0dMV3UzZnZacmdvTlBxRkYxNW9mVHRsUThhUU1XUU1HVU9tNjBGR2pud2RXOW44dGNySlZFdFJYNHNYTDE0dnlGU214TlF0TXlxT2FjZ1lNb2FNSVdQSWRCUElWUHdzT3I1eVhKUkUyUktNbFB1eVBvQm9WWUZMK0ZqWS84V1FNV1FNR1VQR2tPbWlrQkZZbE4ydk1qSnFKRmJkeVZLakcyWDNLMng1enB3NVdYc3RFV1FHekV0RnFIQjdRb05OeUFxQkFOQ2d1MmNsWWhPcTM3MnlYYUlJS25UVVI5QjU5cFhlaVJoWVFLTWNOUm1qMFkwQ0EyaTA2YmgvYVBYUVJJUWVSU2hSUEY2VTNFbkhQNitmbnhlaEVCWEFqQXBpMGpFZkZhaWttRHhJeHplTmNGbUk4UG9iaFJLdmgvZGZOdm15TEZUS0tvSk9WSEF6MnI1ZWd6Y2FmVUlqY3VSSCs5T1JUMFVRWXFLcDFDSFJaU3BnS2FlK2dDT2d5RStqYkgrRksrdmZsU3RYSmxDUTgvL3h4eDgzWkF3WlE4YVFNV1FNbWRacHhJZ1IrWWhHTmNkVUtGTkJBWmRmZm5uMjFsdHZaVXVYTHMxT091bWtiTlNvVVhrOU1wV08wYlk3N0xCRG0weWhHVEtHakNGanlCZ3kzUnd5a3VxTktmSkxrRkgvbU05KzlyTTVVSlNCcitpd3dZTUg1NkhNY3VDcnlabjhPWnBpZS9ubGwzTlluSEhHR2ZuZm9tWFV1VTJKREJsRHhwQXhaQXlaSGdDWmFna2s4cy9vMzNyVGJBS1E2cFhwLzNXZHFzN2NLR1FJbFJBeWMxSVJLbjF1VE1XQ21YMytlVjRpdm1RMG1uVE1FeEpzSWtiSVVEeCtCSkdveVJnaHcvMmo0M005a3pHWi9FZ0lFQktFSkIzLzNENXFNc1lDbVlSQW8wYVZScHFRb2xFcjJ6U3NiUEpsMmVzdHF3Z3FrYU05S3FnWkpYTkdFQ09VSXBVTkZDaGJVTFJhM0paQkJsRmlKNkVRSldNU0d0eWYwQ0drV0RCVDZsUVovOVhhYmJmZDhpbXpuL3prSjlrUlJ4eVJYNmlnMUpyRmtERmtEQmxEeHBBeFpNS1d5b2NlZW1oMjVwbG5OaWRZdHJiY2pDRmp5Qmd5aG93aFk4aTBDSmpOTjk4OGQvb3JNRURYSnVEVUs4MXZ5Qmd5aG93aFk4Z1lNalZybXlucVRNNThSWlZwZWt6dG1WVWs4K2FiYjI0dUUxT212YktXa1JjMEpXS0JTMEtEaWdwaTlybHBmaW8wSWV0enkzV0orRUhSeUVYSmtiLy84emFKNkxqbmVpWW5Ncmt6U3A0a0ZDZ2FZWXFPLytqNExEREovUWtsTmsyanVEM0Y1MEdJRWlwUk1pQ05LZ3RhUmthWTI5UFFSRkNoNGFHUmk0eHkyV1RIUmlGRGFORHh6ZXZsODJIeTZrWFBIcEdJMi9QK0dDaFFGanFORnRpTUhQOWxGSjJiNzBaWkNERndnRkJoNHVXWSt5OHJxQ0hJYkwzMTF0bkVpUlBYdXppbVlLTG9zUk5QUERHNzY2Njc4bndaRmJKVVVjc0tZTlluWDhhUU1XUU1HVVBHa09rR2tORitLbXhKdUNpTFh5TVN0VlhXMUplbXV4UWxwdVJMaFMwcklWTWh6TnJtNUpOUHp2OVZqc3dUVHp5UlEyTDgrUEg1c1MrOTlOTDF5b3N4WkF3WlE4YVFNV1M2Q1dTVTUxSVp0ZWhmVFgwcEQwYVZrUjkrK09HOEc2V1NMZSs5OTk3c29ZY2V5cDU2NnFtOFFyS1NMT1ZuVVZNeDFUV3JqRndxeDFYT3pQb3Vob3doWThnWU1vWk1ONENNQUtQUmlvQ2cwWWV5K0RXeVdaK2VNSm9XMjNubm5aczdXSmIxdzFRdmZmN3h1bFF3K25UVTA1RS9ZRzVUS2dRRzhIaDk3N2dtRVQ4WVF1WFByMjJWS0lKS0JCa21iL0o4Tk9vMCtvUUVqWDRFbFFnU1BCK2JoRVc2N044T1M4VGo4ZnhSNEFMRjVNcEdIZmswY3BGb0JKbE1HVUVtU3I2TWtodkxGclJzdEFBbXI0ZUdrZnNUS253L0dEZ1JRU3hTbzRFQVpTSFVsb3F1amU4S2JWVlo2QkF5dysrOW9xQ0dJWFAxMVZmbmNGRVZaVUhpN2JmZmJxZ0V6Sm8xYTBwRmtoa3lob3doWThnWU10MFVNb29JdSsyMjI3TE90aGd5aG93aFk4Z1lNbDBjTXZLL0tMcE1KZnczT0ZSdW01dktrREZrREJsRHhwRHB1cEJSeVJkMXUxVEUyTEpseXpvZFpQcmZPU2RSMzZYWEpHS0J5dEgvZWxsZDdicmlPNG40VWw3OXdrR0oyRlFyS21qSjVNcXlrR0ZCU0RibEltUzRubzd3eUZITy9RbVZhSHRDSTJyNlJTUEQvWmtjeXMrRHg2T1JhclRnSTQxaVdlaHdmeWJnRVRKc1NoWTUrbWwwSThoRWp2R3k0dkVaZUVCRHlQMFppRUZIUDQ5ZnRpQmxXZWkwZFlITk1rRURaUU1CNmhYZnJQVnVSUVV6b3lacGJaS01XU25aUDNueTVJYWM4NGFNSVdQSUdES0dqQ0dUVEk4TkdEQWd6M1BwREtNWFE4YVFNV1FNR1VPbW0wQkc1ZmlWLzdMWlpwdGx6ejc3N0lhSHlxMXpVeGt5aG93aFk4Z1lNbDBUTWlxdkw3aW9Ddkw2NUw5MEJHVDZMN3U2cmtZdC8yNGlPazc1cGFkam1ZNTBRb1NPZTRwUUlTUUlHU1puRWxLRURndEs4dm9KQVJhRXBLTGtSaDR2Z2d5UEZ6VkJJMFJvVkFnTlF1V3E1dzlKUkVqeDgyN3JaRU1xZ2d5M3A2SGcrMHBERVVHbHJORnZ0SUFrRlowL09sL1U1STMzUzZqUXlKYzEzTzN0MkcvUFpNeEdDMnlPLy9rbGlRaWhDRHBTWGNpb2s2Vkt3RXlkT3JYaDNCVkR4cEF4WkF3WlE4YVFTU0xJbEgzLzNlOSt0L1Bsd1JneWhvd2hZOGdZTWwwWE1uTHliN25sbHRsQkJ4M1VLYURTNzBkekVoa3lob3doWThnWU1sMFlNa09HRE1uOU1PKzk5MTZuaE15Z2UyWWxHdkF2c3hNTitmRlZpZmdnYVlUWTFJb1FpS0JSMWxILzlDdDlFaEVxNzd5K1hTSTJMV1B5SlpNajZXaW5rUzhMbWNpUlQwancrVEw1TVhMTU16bVNSb1dPK3loUUlJSkM1S2huc2w4RWpVYVRNMmxrYVJpNFBuTFVSd1V4dVQrUFR5Tk9SWTV4cm8rYWl2SDQvTHlpSm05bElSTUZCa1QzRTBHbjdQbXE5MjAwcUlEbllqSW1nMGoyLzhXNWlRZ1ZCa0V4TUtBdVpGUnFYeU1ZRmJxVUgyYjU4dVdkWm5yTWtERmtEQmxEeHBEcDRwQ1IxUEpZcFdJT1BmVFFUdVdETVdRTUdVUEdrREZrdWdGa2xHeXBiUDROblF0RG40b2hZOGdZTW9hTUlkTU5JTFBwcHB2bUdmMGJlb2tjK1FQdm1wMklVT0ZOMGlqOWJQV0lSRkV5WkxRK1NxNWtzaVNoOC9icnFTTElyRmc5TEJFZDZZUUNJZkxEbGJzazRubzYrZ2tWUWlKS2xvd1VOWjJLSE9rME9vMG1TeEpxa1hqK3N1TDVDUWs2YnlNakhVR2xiUEpsVlBBeWdoUU5IeUVTQlFMdzgrTHo0L240dk1vNjNxT2lrN3orUmdNTDZwMC9PbmVqRGRVSUlVS0h0cFJRWVlGTUJncElDV1NrWWNPR0dUS0dqQ0ZqeUJneWhrejdRT2FpaXk0eVpBd1pROGFRTVdRTW1iYUZqTHBjNnQ5VnExYTFPMFRZQTNyd3gyQ29saUZqeUJneWhvd2gwODBnbzlMOStyY2psZ2d5VEo0Y2VkL2xpYmllemltK3REU3FOUEtOUW9XQkFYVHNFd284L250dmJKL292MTdmTmhIUHgrUkhRaVpLdm93Z1EwangrSVIyNU1nblJHZ2thRVFhYllwRkkxczJPVEtDUWxrSWxZVVNqVGFOVUFRSktrb1c1Zk9LQ21oRyt4TVNrYk02T244RUdXNGZGY2lNa2tlalFJZXlnUVNFSHAzdlZQVzJVUkJDQkEyZXUyeXgwQWc2dEwxMUlhUGFaR3ZYcmpWa0RCbER4cEF4WkF5WnRvZE1SeTZHakNGanlCZ3lob3doWThnWU1vYU1JV1BJR0RKZER6SmNDSTA5SHI0d1VkUlloOGxxZEN3L3NtWndva2FoRWtIbXNUV0RFckhBSm85UHlLeDlmWnRFREJ5SW1ud3hPVE5LeGlSMDZQaG53VTFDaGV2TEprTkdCUnNqbzFmV1NKVk5saXdMRGE1bmNta0VHUnBwR2hMZWI1VHNHRDMvUmdNSG9tVEt5RER5ZnZuNVJaOEgxMGROMGNvYTJyTFFvRDJpdUgxMC9IcmlzWmg0SGpXOEsvc3NlSHdXeUtUdGxnd1pROGFRTVdRTUdVUEdrREZrREJsRHhwQXhaQXlabWd0dmpEZHV5Qmd5aG93aFk4Z1lNdTBHR1RxSzZGemk5a3dJNG9OaVFjWmYvNmxmb3FqQVpkUmtqQVVxS1RZUm94RW5sTjUvbzI4aUZzemsvalJhZEt6ei9nbU5xRUFtOWYyWDlrbkVRQUNxTEdRaVNKUnRlaFZCcHF6alA0SkVvMURpL2pUU2tWR21FWTN1djJ5Z1JGbEhmNVNNR1NWblJwQnBWRHcrRFRmdEZSVkJKRkpVd0xPTUk3OHM0TXBDaHA4VnowZmJYTGRBcGlGanlCZ3lob3doWThnWU1vYU1JV1BJR0RLR1RNK0dESHNTR0RLR2pDRmp5Qmd5aGt5N1FZYk9JdDVJVk9TTkQ0cU9KVDRJR2owbVIwYU9mRUtGeVpSc2NrYkhQbzA2dC8vTHgrQ29GaUd6NnRXdEV0R1JIem1hNlpqbjlaUnRZc2JrVENaN0VrSk0xb3lTTGNzNm9pUEg5NGFHVEtNRk1xTUNqcnhmcnFmUmppQVRKYXMyV2hDekxHU2l3SVZHQzN3Mm1nd2JHZmJJc1I5QnZmcFpSb21sWlFNRHFETEZPaVVlUHdvRWtBd1pROGFRTVdRTUdVUEdrREZrREJsRHhwQXhaSG9JWlBnaDhNSU5HVVBHa0RGa0RCbERwdDBndzRmT0M0MmFTUEdoODhiM2V2aUNSSHp3TkxKUC9hbHZJanJlQ1JrbVA5SklselhhRER4Z0FVd21Zekk1bFBjVEdjRkdJY05BQTI0Zk5UVnJGREpSOGwvaytJNk1iS09RYWV1bVpUUTZVZE90S0xtUVJwL1BOM29la2RGdEZDcU5RaWFDU2xRd2xPOVRsRHpMNTlOb01qQ1R4WmxNWFM5SWdjOHFna3FaSUlQV3FNTWMvNGFNSVdQSUdES0dqQ0ZqeUJneWhvd2hZOGdZTXAwZk1qeXdJV1BJR0RLR2pDRmp5TFFiWktJdkFUOTB2clJNMXFPNFBSOHNiNHpYeDVlRUJTc2p5TkRSejBBQ0dsa2FaU1pqTXZDQUJURGZmcjEzSWlaM1h2UDdBeFB4UzhGa3pPaDZJK2hFa0tFWUNNQXZFYjlrWlIzN1pZMGMxMGRHTlZyUDk2bFI2TkRvODR2UDk1bnJlWCtSSVlrS2hFWUZKeU9vdHpka29zQURxdXo3eGVOSDM2L0lmbEg4UGx6MS9DR0orUDJ1L29ISUgyeThWLzZnNTJjVEpjcEdqdjZvaVJtRHZEZ0FrQXdaUThhUU1XUU1HVVBHa0RGa0RCbER4cEF4Wkxvd1pPcTk1SWFNSVdQSUdES0dqQ0hUYnBDSmlxcnhRK05EWnJMZXNwWGpFa1hKakhTTVJjNHZPdVllV2owMEVaTWRDUmttVDdMZ1pHU1VDUmtHSHZEOExJaEp5TkN4enBlVzE4ZmtTYXJSZ3BsUklBQUREL2lsall3R1A5OG9XWmZySThkdVdjaEVqdUdvSUdia3VPYjdTOFBCKytQemlaSXhvMENIeU5FZE5RR0xBaE02RzJTaXdKT3k3d01oUkRCRXFnZWR5UFpGUVZuUnUxSTJxSVFzWU9JbGl4ZExob3doWThnWU1vYU1JV1BJR0RLR2pDRmp5Qmd5aG93aFk4Z1lNb2FNSVdQSUZDRlQ3VWhuQVVvNmQzaGlmbWcwZWpUYUxFQVpGYXhrY2lTTkhEODBPdEtZM0xqbVk4TmVMWjd2M2xWakVrVk53QWdaUXZPUk5ZTVRNUkNBMTBjbzhYaDB0RWNGTExtK2JJSE1zdHV6b0NkL2RFUkd1aXdFb3VUSnlJZzJDcDJvNEdYazZPY1hOeW9veS8yNW5zN2U2SDc1K1VRL0NuaitDUHFSSVlzS1dyYTM0ejlLL2kzN1B2QjVFaFJSTWlhaFV5MCtDLzVBaVg2Z1JldjVMdkVIRHhNdENaUnhEOHdveUpBeFpBd1pROGFRTVdRTUdVUEdrREZrREJsRHhwQXhaQXdaUThhUU1XUU1tU0prcXBObmRsM3huVVJSUVVvNnZtaWttZnpJcGwxMGZQLzV0YTBTRVFwMG5QUDRGSS9IL1gvLzUyMFNNVEFoY3F3VE10eWZ6K05ucTBja2lwcWk4WGhSTWlTTlBOZFRVWEptcEtqSkdiOG9VWEliMTBkUTRucnVIMEdNa0lpU0Y2TmtTeG9wZnBHanhrLzhZdlA3UmtNUUpXTkdrT0dQTkQ0L2JoOEZZdENRUlJDS0NsaEd5YnlOT3Y3TEtub2YrUHpxUVVNcUF4bCtWNktnaXFqaEdoUGI2Y2puZ0dQOHp5OUpSRmJVa2lGanlCZ3lob3doWThnWU1vYU1JV1BJR0RLR2pDRmp5Qmd5aG93aFk4Z1lNdlVod3hQVE1jY2JwOUdrVVdkQnlDZ1prcEJZVlFNc2pZakhwK09kUm5QK2k1TVNNVGt5U3M2TWtrKzVQb0lLejgvckl4UzVQaEloeFlLYmtYaDlYQjhadFFnNlpadUNSVkNKdGkvcitLVlJpcjdvaEU1WlJRVTFvK3VQUGcvZVgxUXdrK3NKblNqNU1qTHFVUk83dGs3T2pOYnplVWFRWWJKeUJLRjZRT0s1K1lPanJPTytVVEZJVERKa0RCbER4cEF4WkF3WlE4YVFNV1FNR1VQR2tERmtEQmxEeHBBeFpBd1pRNllJbWVvUGtRK0JScEhKaElRS1JZalEwVTdIUGdNQkNLa29jQ0FLRk9EMU1EazBNdUkwbXBHalBTcFl5ZTE1L3JKR244bW8wVXZPN2N0Q0t0cWVCVko1UGlhbjBlaVZoVTdrNkc0VU9wR1JpWnJ1bFcwS0dEbjJJNlBkcU9PZm90R1BJRk0yZVRKNnZoU1BWeFk2WmFGR1JRVkh5elFwS3lzZW05Zk9kNm1zbzU4UTRYcSt1L3pzSkVQR2tERmtEQmxEeHBBeFpBd1pROGFRTVdRTUdVUEdrREZrREJsRHhwQXhaSXFRcVU0V1pJRkdHbVVhY1RyU0NSR0tScC9RZVB2MTNvbmVlMlA3Uk8rOHZsMGlicy9qUlFVeGYvMm5mb25LUWlVU2pTeEZxRFFxSHA4UW9QT1FMMjZqa09HWHFGRW9SaEJxRkZJMENoRjBhbjJSeWtDbjBZS1BaWnRzUllFTTBmT0pubGZVS0t0c3NpUHZMeXFBdXI3UDg1ajdqc3dHVGh5WUhYN0xsTHJKdDJXVE14a293ZmV6ekkrOHN0ODlIcHVmSFFFY05jempldTdQejRLMlJUSmtEQmxEeHBEcGNaQ1o5cXZqczgyMi9GSnVBUGM2ZHpkRHhwQXhaQXdaUThhUWFSdklUUDNGc2RrMlE3ZkpqZC9tVzIyZW5mYWJieGd5aG93aFk4Z1lNb1pNWTVDWit1dXBXZDl4ZlN0R0w5ZFhiejRvckZWbnlMUVJaS29oRVRuaUc0VU1Bd1ArNi9WdEU3My9SdDlFNjk3b2w0anJDUjFDSnJvK0pwTkdqdXVvNENTTmZoUUlFQldvWkhJb2t6MHBicy96ZHpSMEltaEh4NHVnRTUyUDY2TXZKaUVWRmRnc202eFhOcG0wVVVWUUpHUWlSWUVYVVNCQ2xPeElvMDFEeHMrSDYxdjZQTDc1aTZPelhyMTY1VVp2NEs3OXNxT1hIbHdUeWxGQnpySlE1UE1yOC8yS2ZzQkZ3Q3I3QTQ3WHltZGROb2hJTW1RTUdVUEdrT2tSa0dudHlNK1FNV1FNR1VQR2tERmtEQmxEeHBBeFpBd1pROGFRTVdUK0NwbDZqblpDb0N4VTZHZ250SGg4bnA5UWlhNHZjdlF6K1pKR25vNzdzbzczQ0NwbEljSmsyQi8vY1lkRVVVSE5LQm15bHJPdVd0dytna1lFa2JJRk4va2k4L3FpODBWZnZHajd5S2lXZFJSSFJpOVNCSTJ5UnBLR0pBcWNpS0RENjJHVHNnZ3FVWEpqZEwzUjUxRTJVS0pzMHpNZVB5cVkyY2dQckxLS0lGRVdXdEVQT01tUU1XUU1HVVBHa0RGa0RCbER4cEF4WkF3WlE4YVFNV1FNR1VQR2tERmtESmtpWkQ3OHk0Qm0wYWhIQlM4cE5nRmo4dVhyTllwYVZ1c3ZnWGc5RWZSNGZZK3NHWndvYXJyVjBWQWhOQWdWaXNlTGtqRWpvMTNXOFY4V0dwRUkrZWlMMFpvWHZaN3pOTm8rQ2dTZ0lpTlgxbWhGaXBweU1TbVIrNWROWm8yTVBJMHFyNGZRNFhvK1AwS2R6enRxc2haQkpYTFVSNEVLMFk4SVhnKy9iL1YrNUVTZlRWbEZueVhQeDJzdEN5bkprREZrREJsRHhwQXhaQXdaUThhUU1XUU1HVVBHa0RGa0RCbER4cEF4WkF5WkltUSsrc3VnWnJIZ1pGbEhQN2VuSTU2UUlJUjRQRUxyZDY5c2w0aU8vQWg2Tk9JMGNvUkdsSXdaUWFXc0k1L3JxNHVYU3J6K3lMRmZ0a2xTV2FoRVRkWEtGaFRsK2NvNkx5T0lST3Q1Zkg0eEk4aVVoVTVab3hVcGdoQ3ZweXhrb2tBUkhvOFFJUFFJSFJaaGpDQWFKVU9XTFJES3p6OUs5b3grTlBEOGZENzFuUDBSQk5vYVFtWFBGMEVwS1pCcHlCZ3lob3doWThnWU1vYU1JV1BJR0RLR2pDRmp5Qmd5aG93aFk4Z1lNb1pNVGNqUTBVNmpUeU1mSlZzU0tsRUJUa0tCVGRTb0NESnNTa1lvUkpDaG83OXNnY3ZJa2I5ODFhaEVoQW9MZVBMNGJkMDBqQzk5MUlTdHJDT2YrMGVPL3VqRmJ4UXlFWVJwaEdoMElzZzBtcnhaTnBremdnNnZKeXFTV0JieVVRSmgyY0FBaWdVNENhV282UnMvSDE1dkJNMnl5YTZFREo5dlBjaDBORlRhUTRhTUlXUElHREtHakNGanlCZ3lob3doWThnWU1vYU1JV1BJR0RLR2pDRmp5QlFoVTUyTVNRalFTRE9aOGVsWCtpU2k0NStPL2JMSm1ZVGE0MnNHSklxZ3crc2xKR2preWpZWmk2REU3UWtaUWlRU3owY2p5ZnVKcmpkcUdoWTU5cWtJSXR5K1VhaVVUY1lzV3hRd0NnUW9DeDBhblVhVC9jb1cwT1QxUmdsNjBmT0tBZ0VpbzgzejgzNElrU2hRZ0pEaDhTSUk4dmxFQlZLajV4MGxqOVlMUEdrVU1oME5wVm9GVkEwWlE4YVFNV1FNR1VQR2tERmtEQmxEeHBBeFpBd1pROGFRTVdRTUdVUEdrQ2xDcHJvSkdCMzlkRVN2V0Qwc0VTRkVJOC9qUlpCWisvbzJpYUpBQkVLSDBIdHN6YUJFTENnWk5TbUxraTJqNUUwYWVUN1BoMVlQVGNUblMwWEpsN3hlSm5kRzk4K1hQb0pJRkFnUUpUdnlSYVZSaXE0blN1YU1ITk5SRWNDb3lSdU5SZ1Nkc2xCb3RPa1pqU0t2TjRKSzlQbTFKaUd2RWVqdyt1bG9aeUFBSWNQNzQvMUVnUWxSSUVkVWNKUEg0LzNWKzFIVDFwQnA3NEtiaG93aFk4Z1lNb2FNSVdQSUdES0dqQ0ZqeUJneWhvd2hZOGdZTW9hTUlXUEl0QTR5Ly9YNnRzMmlJNTBGSG1uMHVEMlRLWmtjR1VIbUx5WEY0eE02aEdDVTNGaTJvR01VT0VBalQ2Z3dNSUhyS1JwWkdnVWFmVUtGaWlBWk9mTHJOVjJxNVRpUG1xU1ZMY2daTlpYajloRTBvdXNyR3dnUU9Zb2pTTFExVktLQ2tCRlV5djRJaUpKL294OGQwZk9ObmtmWlFKSG8vV2kwRVZpWndJb0k0Tkd4eXhhNExBc1J2bHUxWk1nWU1vYU1JV1BJR0RLR2pDRmp5Qmd5aG93aFk4Z1lNb2FNSVdQSUdES0dUQkV5MVFVc0NSRTZrbW1rYWNRSm1haXBHUXRvTmdvZG5wL1FvV2prbVN4Sm94dzFKZVB6WWVBQno4LzFFV1FpSXhBbFowYkpvMlVkK1kwMlJlUDV5a0lqU2g3bCswdnhlaUxJUkk1L09vb2pTSlJ0c2hXSjU0OEtYa1pHbitMejV1ZkQ5NnZzKzhicmlReGZGSGdSL1lpSmtvVjVQV1dkNTVHaHJwY01YRGFvb2xHSVJFRXJVYUpxclhNYU1vYU1JV1BJR0RLR2pDRmp5Qmd5aG93aFk4aDBZY2hZbG1WWlZqdm9ULzgvR3grUWZBUHVsZzhBQUFBQVNVVk9SSzVDWUlJPScsCiAgICAgICAgICAgICAgICAgICAgW1szMS4wMjgzOTI3OTE3NDgwNDcsIC03OC4wXSwgWzQxLjk3ODc5NzkxMjU5NzY1NiwgLTYxLjBdXSwKICAgICAgICAgICAgICAgICAgICB7CiAgImFsdCI6ICIiLAogICJjbGFzc05hbWUiOiAiIiwKICAiY3Jvc3NPcmlnaW4iOiBmYWxzZSwKICAiZXJyb3JPdmVybGF5VXJsIjogIiIsCiAgImludGVyYWN0aXZlIjogZmFsc2UsCiAgIm9wYWNpdHkiOiAwLjk1LAogICJ6SW5kZXgiOiAxCn0KICAgICAgICAgICAgICAgICAgICApLmFkZFRvKG1hcF8xOWRkZjM3ZjczZTA0MjZiYTY4ZmY5ZDhjMjg4Nzg2Yyk7CiAgICAgICAgICAgIAo8L3NjcmlwdD4=\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x7fea5242a9b0>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from folium import plugins\n",
    "\n",
    "zoom_start = 6\n",
    "\n",
    "m = folium.Map(location=[y.mean(), x.mean()], zoom_start=zoom_start)\n",
    "\n",
    "plugins.ImageOverlay(\n",
    "    image=image,\n",
    "    bounds=[[34.5, -71.7], [36.6, -69]],\n",
    "    colormap=lambda x: (1, 0, 0, x),\n",
    "    zindex=1,\n",
    ").add_to(m)\n",
    "\n",
    "plugins.ImageOverlay(\n",
    "    image=img,\n",
    "    bounds=[[y.min(), x.min()], [y.max(), x.max()]],\n",
    "    opacity=0.95,\n",
    "    pixelated=False,\n",
    ").add_to(m)\n",
    "\n",
    "\n",
    "m.save(os.path.join('results', 'GulfStreamImageOverlay_0.html'))\n",
    "\n",
    "m"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Bypass the figure creation entirely:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from matplotlib import colors\n",
    "\n",
    "\n",
    "carto_temps = [\n",
    "    '#009392',\n",
    "    '#39b185',\n",
    "    '#9ccb86',\n",
    "    '#e9e29c',\n",
    "    '#eeb479',\n",
    "    '#e88471',\n",
    "    '#cf597e',\n",
    "]\n",
    "\n",
    "cmap = colors.ListedColormap(carto_temps, name='carto_temps')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def colorize(array, cmap):\n",
    "    normed_data = (array - array.min()) / (array.max() - array.min())\n",
    "    cm = plt.cm.get_cmap(cmap)\n",
    "    return cm(normed_data)\n",
    "\n",
    "\n",
    "colored_data = colorize(data, cmap)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><iframe src=\"data:text/html;charset=utf-8;base64,PCFET0NUWVBFIGh0bWw+CjxoZWFkPiAgICAKICAgIDxtZXRhIGh0dHAtZXF1aXY9ImNvbnRlbnQtdHlwZSIgY29udGVudD0idGV4dC9odG1sOyBjaGFyc2V0PVVURi04IiAvPgogICAgPHNjcmlwdD5MX1BSRUZFUl9DQU5WQVMgPSBmYWxzZTsgTF9OT19UT1VDSCA9IGZhbHNlOyBMX0RJU0FCTEVfM0QgPSBmYWxzZTs8L3NjcmlwdD4KICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2Nkbi5qc2RlbGl2ci5uZXQvbnBtL2xlYWZsZXRAMS4yLjAvZGlzdC9sZWFmbGV0LmpzIj48L3NjcmlwdD4KICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2FqYXguZ29vZ2xlYXBpcy5jb20vYWpheC9saWJzL2pxdWVyeS8xLjExLjEvanF1ZXJ5Lm1pbi5qcyI+PC9zY3JpcHQ+CiAgICA8c2NyaXB0IHNyYz0iaHR0cHM6Ly9tYXhjZG4uYm9vdHN0cmFwY2RuLmNvbS9ib290c3RyYXAvMy4yLjAvanMvYm9vdHN0cmFwLm1pbi5qcyI+PC9zY3JpcHQ+CiAgICA8c2NyaXB0IHNyYz0iaHR0cHM6Ly9jZG5qcy5jbG91ZGZsYXJlLmNvbS9hamF4L2xpYnMvTGVhZmxldC5hd2Vzb21lLW1hcmtlcnMvMi4wLjIvbGVhZmxldC5hd2Vzb21lLW1hcmtlcnMuanMiPjwvc2NyaXB0PgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL2Nkbi5qc2RlbGl2ci5uZXQvbnBtL2xlYWZsZXRAMS4yLjAvZGlzdC9sZWFmbGV0LmNzcyIgLz4KICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9tYXhjZG4uYm9vdHN0cmFwY2RuLmNvbS9ib290c3RyYXAvMy4yLjAvY3NzL2Jvb3RzdHJhcC5taW4uY3NzIiAvPgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL21heGNkbi5ib290c3RyYXBjZG4uY29tL2Jvb3RzdHJhcC8zLjIuMC9jc3MvYm9vdHN0cmFwLXRoZW1lLm1pbi5jc3MiIC8+CiAgICA8bGluayByZWw9InN0eWxlc2hlZXQiIGhyZWY9Imh0dHBzOi8vbWF4Y2RuLmJvb3RzdHJhcGNkbi5jb20vZm9udC1hd2Vzb21lLzQuNi4zL2Nzcy9mb250LWF3ZXNvbWUubWluLmNzcyIgLz4KICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9jZG5qcy5jbG91ZGZsYXJlLmNvbS9hamF4L2xpYnMvTGVhZmxldC5hd2Vzb21lLW1hcmtlcnMvMi4wLjIvbGVhZmxldC5hd2Vzb21lLW1hcmtlcnMuY3NzIiAvPgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL3Jhd2dpdC5jb20vcHl0aG9uLXZpc3VhbGl6YXRpb24vZm9saXVtL21hc3Rlci9mb2xpdW0vdGVtcGxhdGVzL2xlYWZsZXQuYXdlc29tZS5yb3RhdGUuY3NzIiAvPgogICAgPHN0eWxlPmh0bWwsIGJvZHkge3dpZHRoOiAxMDAlO2hlaWdodDogMTAwJTttYXJnaW46IDA7cGFkZGluZzogMDt9PC9zdHlsZT4KICAgIDxzdHlsZT4jbWFwIHtwb3NpdGlvbjphYnNvbHV0ZTt0b3A6MDtib3R0b206MDtyaWdodDowO2xlZnQ6MDt9PC9zdHlsZT4KICAgIAogICAgICAgICAgICA8c3R5bGU+ICNtYXBfMDRiOWM5YTI4NGVkNGRmZGJkYzBhMDNhODUyMzBhNTMgewogICAgICAgICAgICAgICAgcG9zaXRpb24gOiByZWxhdGl2ZTsKICAgICAgICAgICAgICAgIHdpZHRoIDogMTAwLjAlOwogICAgICAgICAgICAgICAgaGVpZ2h0OiAxMDAuMCU7CiAgICAgICAgICAgICAgICBsZWZ0OiAwLjAlOwogICAgICAgICAgICAgICAgdG9wOiAwLjAlOwogICAgICAgICAgICAgICAgfQogICAgICAgICAgICA8L3N0eWxlPgogICAgICAgIAogICAgPHN0eWxlPgogICAgICAgIC5sZWFmbGV0LWltYWdlLWxheWVyIHsKICAgICAgICBpbWFnZS1yZW5kZXJpbmc6IC13ZWJraXQtb3B0aW1pemUtY29udHJhc3Q7IC8qIG9sZCBhbmRyb2lkL3NhZmFyaSovCiAgICAgICAgaW1hZ2UtcmVuZGVyaW5nOiBjcmlzcC1lZGdlczsgLyogc2FmYXJpICovCiAgICAgICAgaW1hZ2UtcmVuZGVyaW5nOiBwaXhlbGF0ZWQ7IC8qIGNocm9tZSAqLwogICAgICAgIGltYWdlLXJlbmRlcmluZzogLW1vei1jcmlzcC1lZGdlczsgLyogZmlyZWZveCAqLwogICAgICAgIGltYWdlLXJlbmRlcmluZzogLW8tY3Jpc3AtZWRnZXM7IC8qIG9wZXJhICovCiAgICAgICAgLW1zLWludGVycG9sYXRpb24tbW9kZTogbmVhcmVzdC1uZWlnaGJvcjsgLyogaWUgKi8KICAgICAgICB9CiAgICAgICAgPC9zdHlsZT4KPC9oZWFkPgo8Ym9keT4gICAgCiAgICAKICAgICAgICAgICAgPGRpdiBjbGFzcz0iZm9saXVtLW1hcCIgaWQ9Im1hcF8wNGI5YzlhMjg0ZWQ0ZGZkYmRjMGEwM2E4NTIzMGE1MyIgPjwvZGl2PgogICAgICAgIAo8L2JvZHk+CjxzY3JpcHQ+ICAgIAogICAgCgogICAgICAgICAgICAKICAgICAgICAgICAgICAgIHZhciBib3VuZHMgPSBudWxsOwogICAgICAgICAgICAKCiAgICAgICAgICAgIHZhciBtYXBfMDRiOWM5YTI4NGVkNGRmZGJkYzBhMDNhODUyMzBhNTMgPSBMLm1hcCgKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICdtYXBfMDRiOWM5YTI4NGVkNGRmZGJkYzBhMDNhODUyMzBhNTMnLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAge2NlbnRlcjogWzM2LjYzMTY1MjE2NTYsLTY5LjUwMDAwMDY2NjZdLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgem9vbTogNiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG1heEJvdW5kczogYm91bmRzLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbGF5ZXJzOiBbXSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHdvcmxkQ29weUp1bXA6IGZhbHNlLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY3JzOiBMLkNSUy5FUFNHMzg1NwogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB9KTsKICAgICAgICAgICAgCiAgICAgICAgCiAgICAKICAgICAgICAgICAgdmFyIHRpbGVfbGF5ZXJfZGUxNzlmMTRmZTI3NGI0ZmE3ZDk5MmQ2MmIxNjdiOWUgPSBMLnRpbGVMYXllcigKICAgICAgICAgICAgICAgICdodHRwczovL3tzfS50aWxlLm9wZW5zdHJlZXRtYXAub3JnL3t6fS97eH0ve3l9LnBuZycsCiAgICAgICAgICAgICAgICB7CiAgICAgICAgICAgICAgICAgICAgbWF4Wm9vbTogMTgsCiAgICAgICAgICAgICAgICAgICAgbWluWm9vbTogMSwKICAgICAgICAgICAgICAgICAgICBub1dyYXA6IGZhbHNlLAogICAgICAgICAgICAgICAgICAgIGF0dHJpYnV0aW9uOiAnRGF0YSBieSA8YSBocmVmPSJodHRwOi8vb3BlbnN0cmVldG1hcC5vcmciPk9wZW5TdHJlZXRNYXA8L2E+LCB1bmRlciA8YSBocmVmPSJodHRwOi8vd3d3Lm9wZW5zdHJlZXRtYXAub3JnL2NvcHlyaWdodCI+T0RiTDwvYT4uJywKICAgICAgICAgICAgICAgICAgICBkZXRlY3RSZXRpbmE6IGZhbHNlLAogICAgICAgICAgICAgICAgICAgIHN1YmRvbWFpbnM6ICdhYmMnCiAgICAgICAgICAgICAgICAgICAgfQogICAgICAgICAgICAgICAgKS5hZGRUbyhtYXBfMDRiOWM5YTI4NGVkNGRmZGJkYzBhMDNhODUyMzBhNTMpOwoKICAgICAgICAKICAgIAogICAgICAgICAgICAgICAgdmFyIGltYWdlX292ZXJsYXlfYWEyZGQ2YjA4OGMxNDgzYjgyNGZjYWMxNjc3NzUyMWUgPSBMLmltYWdlT3ZlcmxheSgKICAgICAgICAgICAgICAgICAgICAnZGF0YTppbWFnZS9wbmc7YmFzZTY0LGlWQk9SdzBLR2dvQUFBQU5TVWhFVWdBQUFEMEFBQUE5Q0FZQUFBQWVZbUhwQUFBQVcwbEVRVlI0MnUzUG9RMEFBQXpEc1A3LzlFYUc5MEFOSW9VNmsweGJ1VWxSME5EUTBORFEwTkRRME5EUTBORFEwTkRRME5EUTBORFEwTkRRME5EUTBORFEwTkRRME5EUTBORFEwTkRRME5EUTBORFEwTkRRMEMrNnFnV3pZV3BwQTlnLzdBQUFBQUJKUlU1RXJrSmdnZz09JywKICAgICAgICAgICAgICAgICAgICBbWzM0LjUsIC03MS43XSwgWzM2LjYsIC02OV1dLAogICAgICAgICAgICAgICAgICAgIHsKICAiYWx0IjogIiIsCiAgImNsYXNzTmFtZSI6ICIiLAogICJjcm9zc09yaWdpbiI6IGZhbHNlLAogICJlcnJvck92ZXJsYXlVcmwiOiAiIiwKICAiaW50ZXJhY3RpdmUiOiBmYWxzZSwKICAib3BhY2l0eSI6IDEuMCwKICAiekluZGV4IjogMgp9CiAgICAgICAgICAgICAgICAgICAgKS5hZGRUbyhtYXBfMDRiOWM5YTI4NGVkNGRmZGJkYzBhMDNhODUyMzBhNTMpOwogICAgICAgICAgICAKICAgIAogICAgICAgICAgICAgICAgdmFyIGltYWdlX292ZXJsYXlfYTExOTA3OGZiM2YyNGRkYjk5NWEzNzlmZjNlZjA5Y2YgPSBMLmltYWdlT3ZlcmxheSgKICAgICAgICAgICAgICAgICAgICAnZGF0YTppbWFnZS9wbmc7YmFzZTY0LGlWQk9SdzBLR2dvQUFBQU5TVWhFVWdBQUFHY0FBQUJUQ0FZQUFBQ0NxaGxDQUFBRWprbEVRVlI0MnUyYlRXN1hNQkRGY3pGdXc1RzRBdUlLWEFDeFFCd0FJWlpzdTRZTmY5UktrVkxYbm5sdlBodzdjYVJaRkNVV25WL202MDI2YmJOZkg1OGVxTDM3K3VPTmxXZlY3dm4wKytjcisvdjRWN1h5UHE5dGw3aElJRytnRkpjR3B3YnErRy9JODhnOTIyVXVFQkI2bkFabUIxS0xHQmJPdFNOSGdHTUJzOFBSb2tWTFowamtTYzl1bDdvRU9CRlFhbkNRbElVQ090NTNQVGdGSUV2VVdOOXlOSjNkcnlHb3dMR21Nd3VZRmhCdk4zZXRhS20weEpsZ3lsUWt0ZGtTcU5ZOTE0UkRYaEZwQ0lYQzJIYm5TNEtDcEtFTUlNLzIrUGJseFc0TlJuTjhEeWc3aUpyZEtrTFE5SVhXaGQySmtVQ21nNU10d1dnUkV3RUdCVElOSEZTWE9vSmlvNlBXUVIxL2JqazZJMUttZ01PKzdaYVVwWUZCMjEwRWxnWmhtb1lnV25iWGFna3F5MFFXZmZTWjI0Q1JPakF0Z2pTblo3WFUwN1MybmdrY2lacnlQaW45b09scGVqamVGaGZkVEtJcHkxcThNeUJOVi95WlFzMmtLY1RCdlFFTkxaOUVXUStwSlFQU2tIVW1Ha0ltSEFUVWRIRFlMMXM4Z0R6bk1wRVJIVVdYUzJkZU9FZG4vdm53R2JJV2pObzlVOER4UW1uOXN0YTVoUVZpTVJUUWFRcUJkWE9JekIrV2hxQVhFS1FUTEgrLzA4QkUxZ0FKaGpSZ0lvN1ZuRzVKZDJpdE9oVU1xd0tnaFZhYit0RzNIWFU4QzJZb1ZicVd5dEFKSCttRTBJaXpwQ0cyOEVlQjZRS25WZnladW5NV0dBMlFKWjBORFllVjQydk9ZQUJKanZ6MS92dUxJVTVsYTlId2NMeFJJemtKQmRTQ3M0UFo0V2lPOHM0OFE4R1JaaGx2T21zNWcwbHBUTlI0MCtKUWNEenJZbTNRUktkdnpiRk0xTEJkM3RCd01qL0VRK1FSeEtrV09KNjVaaWc0MFh2NGFEZzdJSXRUTStFY3owb0R3KzdqSTBSSjFwbHM5Q0FEYWhTZ1ZEaFpZQ3g1UDZyVFFoU0VxSlQyZkZaS0krRGRvYlQyK3ExV1dPdTZ2RUtsMU1abk5nbXBLUzBTanBTV0pFQlJTckkwWTJVQlNvT1REUWFkVzFBZzZOeGpnWDhaT05aZHkrN1lvNE1aSytHY01aaVd6NlkxQTcyZ1JGa05EcU8zUlRZRjZYQWk5dmJXS01qZVpQWlNDOUtVQVk5ZVZpdnlXWTczZG01c1MzODZIQlNNMUhtZERTZTYxa2hudGw3YWJuQVE2YjRuR004d0daVWFoNGdjTmxwcVlEUllGcURvTWcxZFNVZlhuUzZ0ZFBtZmw2QzBIS3dObVphSVk1Um1Cdkp3cWpTYTBqUXdMSnpXdVJFYTNKbUx0aTV3SXFLbUpmRWZaeUxtSEUyTzZWMWJVdUgwaUJwTm5tSEJaQXk2dzhMUnB2MElNRkt6d0theW5nUHJhWENRMlNZYXp2RmVhNTBaUVdWSWc0TitQZU9wTXdqUTBjRjRBS1ZFeTI0UlEyWUVtQkhnV0NDbFJrMzA5Tzk1ZGdRd0xDZ1RIRVJoam5qYmUwazFvd0l5Z2RIbWd0SEFJQXIwOUhBa1pkbFN2RDFhMm9KenVGcVJ3a0N4dHNCWE5qZWNLREFMRGdlSWhoSUJac0hCSU5GUlk2a3hDdzYrUjRMZ2VEc3lSaE9iSFZCa202N0NrZlliNkJiVHFqYXYxTmFBdzdTamFNcGFjQUxnc0xNQzBvV3g0dVlDVTlRY2FVaGpwUHdJT0d2bUtlQll2eU5qNFhpV2FWZlUxTXh3b29aTGkrUHZJdHNnM3k5c0VXQ1FobURWRlA3RHhjMExaVGs5RnM2clB3SHhRRmx3Y3EwTG5BWDBKRGdMektSd1Zob2NFTTZxVHdQQVdSRXpJWndselF5WTFyUmhOT29idEFVbkliVXRPRVk0WHBVQU9XTTVuTFAvVE5SKzVlSUdXTndBQUFBQVNVVk9SSzVDWUlJPScsCiAgICAgICAgICAgICAgICAgICAgW1szMS4wMjgzOTI3OTE3NDgwNDcsIC03OC4wXSwgWzQxLjk3ODc5NzkxMjU5NzY1NiwgLTYxLjBdXSwKICAgICAgICAgICAgICAgICAgICB7CiAgImFsdCI6ICIiLAogICJjbGFzc05hbWUiOiAiIiwKICAiY3Jvc3NPcmlnaW4iOiBmYWxzZSwKICAiZXJyb3JPdmVybGF5VXJsIjogIiIsCiAgImludGVyYWN0aXZlIjogZmFsc2UsCiAgIm9wYWNpdHkiOiAxLAogICJ6SW5kZXgiOiAxCn0KICAgICAgICAgICAgICAgICAgICApLmFkZFRvKG1hcF8wNGI5YzlhMjg0ZWQ0ZGZkYmRjMGEwM2E4NTIzMGE1Myk7CiAgICAgICAgICAgIAo8L3NjcmlwdD4=\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x7fea523d63c8>"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "m = folium.Map(location=[y.mean(), x.mean()], zoom_start=zoom_start)\n",
    "\n",
    "plugins.ImageOverlay(\n",
    "    image=image,\n",
    "    bounds=[[34.5, -71.7], [36.6, -69]],\n",
    "    colormap=lambda x: (1, 0, 0, x),\n",
    "    zindex=2,\n",
    ").add_to(m)\n",
    "\n",
    "plugins.ImageOverlay(\n",
    "    image=np.flipud(colored_data),\n",
    "    bounds=[[y.min(), x.min()], [y.max(), x.max()]],\n",
    "    opacity=1,\n",
    "    pixelated=True,\n",
    "    zindex=1,\n",
    ").add_to(m)\n",
    "\n",
    "\n",
    "m.save(os.path.join('results', 'GulfStreamImageOverlay_1.html'))\n",
    "\n",
    "m"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}