Plot_Shapes.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "%reload_ext autoreload\n",
    "%autoreload 2\n",
    "# %matplotlib widget\n",
    "%matplotlib notebook"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import geopandas\n",
    "import os"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support. ' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            if (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>');\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option);\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width / mpl.ratio, fig.canvas.height / mpl.ratio);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function() {\n",
       "        comm.close()\n",
       "    };\n",
       "    ws.send = function(m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function(msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data'])\n",
       "    });\n",
       "    return ws;\n",
       "}\n",
       "\n",
       "mpl.mpl_figure_comm = function(comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element.get(0);\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
       "        return;\n",
       "    }\n",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable()\n",
       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
       "    fig.close_ws(fig, msg);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message(\"ack\", {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>');\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items){\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('tabindex', 0)\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager)\n",
       "        manager = IPython.keyboard_manager;\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.find_output_cell = function(html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i=0; i<ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code'){\n",
       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] == html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "}\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel != null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4nOzdd3BTZ9o2cO9kdr+Z3Zl3Zvfd2ZmdnUSQhUACqYQE2DdhyWaXtE0npJBK2rK7CaTJ2MammlANphoDpmObZowsy5Z7kXvv3Zbcq1zVz/X9IetgudsqR9K5fzPPTAJGOjo61rn0lPtxAyGEEEII4RU3rg+AEEIIIYTYFwVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBICCGEEMIzFAAJIYQQQniGAiAhhBBCCM9QACSEEEII4RkKgIQQQgghPEMBkBBCCCGEZygAEkIIIYTwDAVAQgghhBCeoQBoAYPBAIVCAaVSid7eXmrUqFGjRo2aEzSlUgmFQgGDwcB1lOAMBUALKBQKuLm5UaNGjRo1atScsCkUCq6jBGcoAFpAqVSyFxDX32Z6e3tRpWiDV2gGHnK/gXs3huLejaEIkBZyflxct9Oxxbh3Yyie3hqO7h7+9tZ+d0mGezeG4gmvMJfqtb6YWMpe7497hSG7spHzY5puu5NVjZ03s/H9JRnePhyDL08n4UHh3d/fqdqD7jfw4l4JNl1MRYC0EAlFdWhq75rw+ZRKJaoUbQhNrcAnAfEQbLr7WM/sEGHjhVTsu5OLy0lliC+sQ0vHxI81urV2dOMJrzDcuzEUP1yWzfqcyFs78ZjnLazYfgf3bgzFfSOO8SmfcM7fM2q2b4fEeWbX+fvHYlEpb7Pqc5g6cJRKJddRgjMUAC3Q29sLNzc39Pb2cn0oZlRaPQRCEQRCEV47msL14XCusWeIPR+NPUNcHw5nsuq62PMgKW7h+nCsokDRgwc8xRAIRfjqQjY6+9VcH9KUqtr6cTKhGi/7J7Hvx+i2YncsXjuaggVeYsz3FOP/9sTis6BMhGTJkSfvQb68B90DGjAMY/GxfBaUibnu4x/H/Zsj8Ip/MnZFlKK8pW/Sx9LqDdgWXgKBUIRXjyRbdFzdAxp43CzEYh8JeywPbYlEZFGzRY9LnEdmXRe+C8nHnzdHQCAUYcmOaMi7Bq32+I56/7YnCoAWMF1An55KQFGj0uIPY2u5mtHAfmieTKjm+nA4F1XcjEe3RmHF7lho9fyd7wEAm28WQiAU4aldUrT1qrg+HIu53yiAQCjC3w8mQK3Tc304k+oe0GDd6XSzgDXPIwKfBWVi881C+Ekr4B1WhJMJ1VBp7ftaelVaSEtasT+qHJtC8vBuQBqe3hUzJhD+80gyDkRXIKWqA/KuQfSrdbiWrcD6c5lY7H03rO2KKLXKcTEMg85+NeRdg3Y/J8QxVLX14x8HEyEQirDmhMxqj0sBkAKgRUwX0KsHoiEQivD2iVSHuKnuuGP8Fv7TtQKuD8UhPLc/HgKhCCGZDVwfCud2ikrMbtKlzb2oautziOt2Ng5JK9meKlGBY/cOud8whu+57iKsO52Oc6l16HDwHsumniGE5TXiywtZbE/MZG3JjmjsFpdRWCNW1dgzhHkexusvtbrDKo9JAZACoEVMF9B3F1Mxf3gY6pk9cbiYVo/2Pu4+2I/GVUEgFGFjcB5nx+BIHh4eRqpqm3wIiw88hnsAR7f5HmIUNTrfXJghjR7rz2UNhw+pw/TCj3ZlRK98bFkr14czK+19alzNaMA3V3PxzJ44zPcwfuY9ti0K+yTlKFQooTc45vknzs/zViF7vYkKmmGw8FqjAEgB0CIjL6CK1j4s8707ZDLXXYT3A9Os9m1lJu4UNEEgFOH1YzT/DwD7nnAZyh2FclAL/5hKfHQmA/84mIglO6LZ8+Mrts6wnb0MqHU4kVCNvx9MgEAowgOeYocMIAYDww6P+lppaNQRMAyDAbUOOp5PqyD20afS4tWjKezn1YuHknAyoWrW8wIpAFIAtMjoC6hPpUVAYjX+eSTZrHclssh+E+4NBgabgvMgEIqw4VKO3Z7XkZl6Kuo6Brg+FIe0Yncse60609DdyPl08zwicEFWx/UhjUvRPcgep7MOtRPiCFRaPQ5EV2CxjwRrA2RmYfClw0lY4CXGMt8YfBeSj4SK9kl7CSkAUgC0yGQXkLxrEF9eMA5NrfZLtMvxZNd34bUR35BiSp1zqMnaPjmbgdePpSAsr5HrQ3FIF2R17DXz3yu5Fg+t2MurI75oOXKwOpdqPL+r9sc77BA1Ic6kT6XFzRwF3g1Im3AFu2lKlndYEfZJynEgqhyZdV3sY1AAdOEA6OPjM6bg44IFC9i/V6lU2LBhA373u9/hN7/5Dd588020ts4sME11AZW39LEX4oHoCqveWBmGQbNyCDkN3TidXIs3jt0Nfou8Jbid32S153J2eyLLjD2il6lHdCKBSTXs9ROSKef6cKYltbpjRC+74y4A2RVRSiWZCLGR9j41RAXNuJBWj5KmXqRWdWBLWJFZCaGR7aXDSfghNB+bgzMoAHJ9ALbi4+ODRYsWoaWlhW0dHXfn43399de49957ERsbi+zsbCxbtgwrVqyY0XNM5xvEgahy9sL7ObJs1q/HhGEY3MhR4K/74sedyP/TtQK0OnBvCBcKFUp2XmZiRTvXh+OQdHoDex0di6/i+nCmJTy/kT3m0CzHDa2mhTdLd0ppvhwhdjKo0SE8vwnb75TA42Yh/nMl16y38N6NoRQAuT4AW/Hx8cGjjz467t8plUr88pe/xLVr19g/Kysrg5ubG9LS0qb9HNPtQjb1rizwEk/YCzik0aNXpR0zRKTS6pHb0I3QLDncbxSw9ZAEQhH+vDkCK3bHYt3pdAQm1Tj0MBjXTPMi57iL4HWrCE08Lgg9HoZh2DILG4PzcCS2EmsDZDibUutww5btfWp8cjYDi30keHx7NNackGFI47hzFy+m1bO/s19dyHa480kIXzT1DOFOQROOxFZi6/UsCoBcH4Ct+Pj44Ne//jX++Mc/Yu7cuXj//ffR0GCsAxcbGws3Nzf09PSY/Zv77rsPBw8enPAx1Wr1uFvJTHUBjZwEPnInimblENJqOiG8XsB+M1myIxprA2RYGyDDC4eS2Npbf/n57kT9+Z5iHIiuQL9aZ8EZ4he1To/vQvLNFg18eCYD34fmI6q4hc4lwF5jzx9IMOtZ/vflHFzNaED3gIbrQ4ROb8Bqv0T2PQxKrXX44t4Mw5jNs6xp7+f6kAjhPZoD6MIBUCwWIzQ0FAUFBZBIJFi+fDnuu+8+9PX14fLly/jVr3415t8sXboUP/3004SPOd68wulcQAzDYOlOKQRCEfykFQCMBXk/DcrEi4cSzW6247XHt0fj4zPp2BtZhtv5TQ5xI3ZWqVUdZqvHBEIR/rovHm8dT+X60DhX3KTEjjslEF4vwIuHkvDI1iiz8/TULin6VFpOjzG9phMCoQiLfSSobHWuuo6LvSV4dm8cTtDuPIRwjgKgCwfA0Xp6evA///M/OH369KwD4Ex6AG/kKJBe04kBtQ6BSTVsT97VDGMv5OPbo7HAS4y1J2X45GwG4svbMKDWIbehG7dyGxGe3wRpSSsU3YM0ZGQDFa19OJtSC+H1guGFM5EUrEdhGAZ+0gq8c1KG+0fsAvHa0RRkjVhNZ08HoysgEIrwQWA6J89viU0heexcVGcLr4S4GgqAPAqAAPDkk0/C3d191kPAo010Aam0enYF0gOeEWwR2A2Xc9hCtaaq5lvCiix/YRZiGAZVbX2IK2+DuLAZ+fIeDGr4MSSq1unhc7vY2KvkLcH2OyV0cx5HWF6jWdFogdC4ndm60+n4+GwGfgjNx1cXsrHhUg6+C8mHr7gUQSm1CM5swMW0eqvUFxQVNLNTJS6nO9+2fgzD4MMzGWyA1egce+iaEFdGAZBHAbC/vx+//e1vcfjwYXYRyPXr19m/Ly8vt9oikM5+NX66VoBFw8Hv8/OZuJzeYLYApEWpQkxpK6e9e10DGpxLrcUze+LGDDvPdRfhZf8k7I8qR4vStReXJFW244VDSewikcXeEvzjYCJOj1pY0+ng+7baGsMwyG3onnLKwnjt9WMpkFV3AjAuilhzQoadohJsv1MyrfJIdR0DeHBLJARCEdxvFDhNrcLR8uU9bG/qv0d8ISSE2BcFQBcOgN9//z0SEhJQV1eH1NRUPP/88/j973+P9nZjGZCvv/4a9913H+Li4pCdnY3ly5dj+fLlM3qOqS4gjc4AedegQ5Z+KG5SYrlvDNaclLELS146nIQ3jqWM6elZ7CNBdj03Q372YjAwiC1rhcfNQrNSAXOEIrzin4wr6fX4+4EE6h0EkFjRjq3hxTgvq8PNXAWuZDTgaFwVLsjqEJRSi+Px1fC5XYwNl3LG7IqzakT5or/tN/73xbT6KZ/TtN/v2gCZ04em8PwmNgS+dyrN4RexEOKKKAC6cABcu3Yt/vjHP+JXv/oV/vSnP2Ht2rWorr47+dpUCPq3v/0tfv3rX+ONN95AS8vMtmxz5gvoP1dyjZXSf47FrdzGMUO+LUoVbuYq8LJ/EhsQ1wbIcC1b4fQ34Kl09KtxXlaHV/yTx/RkzXE3lknpovmC05ZZ14W/jbOy2DT0LhCKcDC6wmyF/Eh1Hf1471QanvaVosJFAvil9LulYTxvFXJ9OITwjjPfv63FZQOgPTjzBdSv1mFLWBGUg5Ov6uxTac223BIIRXjzeCrSazqh1jlu7TVrkXcNwi+6gl3FPbIG43cheTiZUI2GztltRs4neoOxFIrXrSJIS4xTHzr71XjAU8ye02f2xI07JcLXRXfS2D9cJP7pXTFcHwohvOPM929roQBoAb5cQHqDce7X0bgqdl6jQCjCo9uiEJTqeIWCbaWqtQ8/XsvHU8Nh8KkRofD1Yyk4mVBNxbhnQKXVY6eohD2Hy31jxvQux5S2Ys7wkHxEoeNu9zYb606nQyAU4csLWVwfCiG8w5f792QoAFqAjxeQonsQG4Pz2CC4LjAd8z3FWHc6HfujynmxUIJhGOTWd+FoXBXWnU5nA4pAKML9myOw/lwmYstanXahgj1UtfWPWXx0LrXO7GdyG7qxwMvYQ+hx0/WGSU11FpMqaXtCQuyNj/fv0SgAWoDPF5DewOBIbOWYOXLLfWMQltfIq5p6rb0qXEirx5vHU83OxXP74xGe30RBcBSDgWHD35IdUnxzNRc3chRmPcnyrkE8vt24GOmjMxkuuVDC9PoceR9jQlwVn+/fJhQALUAXENAzqMHt/CYEJtVg5d67PTp/3hyBRd4SCK8XILKoGdElrWjsGUJwZgO7AOD70Hy8dTwVLxxKwrdXc11i+LSqrQ877pSwtR9NoXivpAxVba6xgMFSDMOw5+Z0cu24P/Pt1Vz2Z5RD3O4+Ym06vYHdGm6uu8jlyywR4ojo/k0B0CLOcAHp9AbElbchp6Hb5s+lHNLCO6xo3LqC02mPb4/GlxeycF5WZ5XCwVzqU2lxOKYSD/tIzF7jwz4SfH0xG9U83Q9WOaQ125M5qnj8lfemsi+mc/ZZUCY2Bufh35dz4B1WhE+DMrEpOA+3chudbtrBrhHzHr++mM314RDCS85w/7Y1CoAWcPQLKKehG6v9EtmehtqOgVk9Tq9KizPJtdOuZ2gwMChp6sXRuCp8GpSJ14+lYNVwzbcHPMXwvFWIreHFOBxTifD8JlzJaBhTe/AvP8dia3gx+lRap15kotLqISpoxmdBmex2gAKhCIu8Jbzr+VF0D+Kvw3UA57iLsD+qfML3trNfjR+v5ePRbVFTfnF471QadkWUOs3uNd+H5uO5/fH46Ew6BtTOccyEuBpHv3/bAwVAC5guoAOiPGTVdXEeVKra+nAgugK+EaV4a9R8NIFQhOBM4/ZZegMDjc4And6AtJpOiAqaUdLUO+HxHxguVyEtaZ31sTEMg0KFcsK5gSqtHrFlrTiRUI0lO6Rjjv1atmLWz+0oBtQ65Ml78OzwUPnPkWVcH5Ld9Km07JeAFbtjp72XsN7AIKehG+dldfCPqURAYjV+vJaP7XdK8HNkGVb7JWKZbwz7uM5Qksd0TYsKXGtVMyHOhAIgBUCLmC6gezeGDm9RVcjJZPX6zgGEZMrN5p2ZVqR+H5qP/w4XfZ7nEYG/7otndyGYP6IGm0AowpqTMtzKbTQrcpw3vPXXHHcRylrs84syqNHBf5zh07/8HItChZLzoG2JnkENu9PI5+f5U/5jT2QZBEIRlvnGWL3nM6a0lQ2B606nW/WxrU2t07PXc0kTf288hHCNAiAFQIuYLqBPAuLxbkAa5ghFeOFQEnLtMN/OJCilli1Dcr+7CA9uicS28BJcTKtH0/DOCspBLdackI07fPbotii8cSwF8z3Mw+DTu2Lw3qk0CIQiPH8gATvulNjtNZnoDQwii5rx5E4pu2WdQGjcozi1usPux2MN1e397OsI4cnqz/Y+NbuP70Rz/iyVUNHOnldFt2P3App21zkcU8n1oRDCWxQAKQBaZOQFJC1pMZuv9GlQJi7I6nAtW4FbuY24lduIKxkNOJ1ci6NxVdhxpwSfnM3AmhMyfHM1F/Kumd+0ihqVbHD7+8EE7BaXTvg4DMOgtLkXyZUdaOwZQotShbqOAXZeX7NyCHuGh9RGh8TFPhJ0cDzR/nJ6A172T2J3jpjrLkLmNIcRHQnDMHhoOAyVunAPEMMwKGpUYldEKf7ycywEQhFePZJss95b5ZCWrU25wEuM87I6aHSOWTrmaFwVBEIRNgXncX0ohPAWBUAKgBYZfQF19qvxXUg+O8Q3k/bx2YxpP69GZ8Dp5Fq2SO7rx1KsemNVDmqRXd+Fa9kKRBY1sz2JjqCzX81uTffjtXyuD2fGGIbBOydkuN9dhOJGJdeHMyMDah1b01CnN6BXpUWfSmu2e0dtxwBOJlSP+SKx2EeCIhu/3sy6Ljw3PM9QIBRhtV8i6ma58MmWvMOKIBCKsFNk/151QogRBUAKgBaZ6AKqbu/HXkkZ1p/LxEdnMvB+YBreD0zD+nNZ+M+VXPx0rQA+t4txNqUWH53JYOffTUdkUQu7iEAgNBbJ5VPRZQA4HFPp0HPoTL1fwZkN2BpejG+u5uKL81nYFJyHA9EV7HsXltfI9aFOimEYZNd3Y39UOV46nIT5HhFYuTcOH53JMJtvOsfdOE3gbwcS8MaxFPbP53uKseFSjl0Lg/epjKWITLtsPLVLiqo2xyq5s+NOCTtfkQ/7aRPiiCgAUgC0iKUXkMHA4O3h1bo/XS+Y9GcZhsFucRl7c12yQ4rL6Q1j9k7lg6y6LvY82LpXaabkXYN4fUQImqw58iKAnIZuvHbU/HWYel4na6/4J2Pd6XRcTKuHcpC7As5tvSr8/WACBEJj6aFNwXkO05Nd2tzLLsRaGyDDkIZCICH2RgGQAqBFLL2ADo7oDfr35ZxJfzY4s4H9Wd+IUt7XD3vpsHEi/R4HK6UysoDxB4Hp2HGnBIFJNbiYVo9j8VVwv1EI9xuFiCmdfUkdW9DpDWjrU6G+cwCnk2vZFeILvMT4z5VcXMtWIKe+CxfT6nE1owG5Dd1QafUY0ujRrBxCXHkbJMUtDrWbS0e/Gu+MWDz0133x065laWsjF63sFjvWNUwIH1AApABoEUsuIIZh2JWR8z3EqGydeJuwKxkNZr0sBPhmxFZh1uxJs2QuZXx5G7siO768zWrHZIkhjR6y6k4oh4wFta9lK3AmuRYNnYO4lq3Aj9fysWp/vFmRalP7/HwW2vuca5eN0QwGBrLqTrzsn4yX/ZOwW1yKfgf58rT+XCYEQhGW7pRyfSiE8A4FQAqAFrHkAlJp79YDWxsw8fw/g+HuvqmvHkl2uJ4jrgSl1JqFldnucjKarLpzVv9OXNjMrsjefLPQKscyXQzDQDmkRVuvil2UUdrci0PSSqzaFw/PW4XwulXElvURCI1bkI03hPuApxgv+ychKGX6O784gy/O3+2ZXewjwYHoCrvOv2tWDqGqrZ89p4MaHf5xMJF9Lwgh9kUBkAKgRaZzAVW29mHDpRzUjLP3q2m3jm+u5k7471OrO9gblzOWPbGV0XvKCoQiq8ylKm5SzriOY768hx0y/epCtt32MVZp9dh8s3Dc3rvx2lx348KDoJRanEqswVvHU+ErLkVsWSsae4Zcej6p3mDs/Vw1YpXw0p1Sm85TVA5q4R9TabYyeYGXGG8eT8W/L+VAIDTW4Rz55aVrQIOchm4oh7TQ6AzsqmtCiHVRAKQAaJGJLiCd3oCM2i4ciK7AXknZhDtomOaLnUutm/A5Vg3vnfrYtii6GYyiNzD4LCjTLOREFlleaPj6JNvOtfWpEJolh3dYEbbfKcEr/sl4/kACW8rH1iGqpr0fB6Mr8FlQJnttTNae3CnFD6H5uJhWj1YHmp/HFYOBwaX0evb8PLotCudldWa731iiqWcIQSm18LhZOGannbHvTTQ+PpOBwKQatPWpEJIpx4rdsWY/89i2KNzOb7LKsTk7V+qRJtyjAEgB0CKjL6ACRQ8Ox1RiT2QZUqo6pgwDLxwyLmT4dpIewK8uZLM3KlfuobHEyAUyAqFoxvPWlENalDbf/RA4nVxr9vc6vQFbwoqwckT5HYFQxC4wWOgViY3BeWjrs27A0uoNiCxqxj5JOc6l1iEgsXraNSbnuBtbecvEc0v57Fh8FR7fHm12zlb7Jc56VXm/WocTCdVs+RlrtoVekS4b3tU6PWTVHVAOapFV34UHt0TitaMpZoXns+u78dGZDMz3FE/65YyQmaAASAHQIpZcQN0DGvYD/kRC9YQ/V9ykZH+OaoZNbORQ+ZvHUrFTVDLh8F5bnwpRxS04l1qHkEw53j6RijeOpWBreDE8bhZinkcEApNqUNsxALVOj/KWPrZe4xx3EV46nARfcSmCMxsgLmxGs9Ly8iIGA4OwvEZsDM7DB4Hp+O+V3DG9QZO1p3ZJcSC6An7SCtQNH3dFax/nO7g4MpVWjzeHp2GMHBaeTk8TwzDQ6Ayoae9HYFINluyInnGwm267f3MEChQ9djgj9qMc0mKfpBwLhntJF3qJ8fn5LHZhnGmrx7KWXnZu7SNbo1AxyWI5QmaCAiAFQItYugjEtJPHREPEAFDV1geB0LhS2FG3tnIUuQ3deGZPHLsLxf/9HIu3T6SiT3U3CG4Z3oVhdHt6VwwCEqtxPL4aQSm1OJNcC9+IUnxzNRdFjUpUtvYhprQVyiHrzBkbudo4vaaT3R92tm2OuwgZtTRHdKaUg1r2HJrmUu6KKJ1wHmehQom3T6Sy29vZo/3rkustEjmTXDvpa37YR4L3AmRY5hsDgVCElXvjHK7mJ3FuFAApAFrEdAH5R85uSzLTHK70molXnur1BizdKZ3y54gRwzCIKW3FC4eS8PLhJDyxPRqhWXIAxpWYC70ix73hLNlhPMcvHEqyyyrevZIyhOU1mpWzsbQ9sjUKP17Lx6X0epovOk0tShV7/jxuFLL//W5AGgY15uVi8uU9Y4aN7dEe3BJpsz2U7UlW3YHNNwpxJrkWap0e/5xGYXFTMHfkounEOVEApABoEdMFdO/GULx6NGVGPTClzb2YM/wBJ52ktMupxBo8sjUKT2yPRs4MV6fyGcMw+C4kD3PdRVjkLUFYXiM+CEyf8CazNkCGf1/OwU5RCQ5ElSMopdYmK0QTKtqx4XIO/nEwkR1WtkVbc0JmlaFpV9Y9oGFLsbzinwytTo9r2QosGt7m7t2ANHZleU17v9kKYnu3ThcYyg/OlLOvJ6asFUfjqqb12o/HVSKuvA059d3UC0ishgIgBUCLmC6g+7+/DoHQOFfnvKwO1e39k35jb+tVYbHP3b1Uc+onDnbrz2Xh//bEOtyOF86gX63DutPpeMAzAmsDZH+bgF0AACAASURBVNO+2S70isTSnVIs9IqE+41CSIpbEJIpx25xGTbfLMR3IflY7ZeIdafTkVHbheImJeRdg1Mej364puMze+LGzD2zZntoSyS+vpiNb6/mmg1/k7u0egP7HizdKUVD5933L6ehG4u8JXj9WApOJ9Xgs6BMdh4aF23F7liXWAAm7xpkX9P9m2d+Hky/w8fjqjjdZpC4BgqAFAAtYrqAapvasWG4rpepvXk8FY0T7D06ci9bgVBktgJ1NNOqwq3hxbZ6GS6to0+NF4bnBNqyTWeeFsMweHpXjM2OYcmOaJQ09bJD3mRiJxKqIRAai0KP3oXHYGCQXNkOcWEzZ6FvZDuVWMPRWZodlVaPkKwGfBqUiTUnZXh8ezS8bhUiLK9xRq/bNEd6jrsIz+2Px5oTMty/OQIPeIrhG1Fq9px6A4PajgEqFUOmjQIgBUCLjLyAtHoD/GMq2Q3oTW3pTime2x+PfxxMREBiDaQlrdDo9Fi1L479cOuZpAbZh8PDhO8HptnxlTm/nkENbuYq8H97bDtZ/2EfCXzFpdPukZhq8rsl7UhspY3PqmtQDmrZcDF6BX5jzxC7Y4qpTBOX7RM71Ja0JoOBwbOjyiUJhCIcjaua9pw/U/vmai40OgNb/SCj9u4X52f2xLF1OF88lITlw4tF3jyeivD8JoRkyc2mQHQPaCAtacWl9Hqr1Xwkzo0CIAVAi0x0AZW39OHdgLQxH2imxRyvHknG/qhy9s+jSyaeA/iKv/FD86sLrrcS0NqCMxvgdasIX1/MnvbuGLNpj26Lgp+0AkmV7TPe9UNvYGwy/PtDaD6tEp8medcgu2dzUkU7GjoHEZIpx+abhXhoy/iLhLhoD/tI0DPoXGFlZDkmS1t5Sx/SajrZ61o5pMVfp1H83NTmeUTgZf+kMf9msbcEsWW0pSbfUQCkAGiRqS6gnkENChQ9SKvpRExpK7aFF0+4CrVpguFi07fm9ecyXWIl4EzUdgxAOaSd9HVrdAa09amgNzD41oorasdrD22JxJHYSvSrdRMez3R09qutuqBg/bks3l0blvrP5Vy84p+MOe4ifBCYPuWuHVy0A9EVXJ+mGVt/LtPq5yE4s4F9/D6VFoeklXjreCo+OZuB8PwmxJW3IamyHaKCZrx6NAUvHErCi+P03j436nfuTkETfWniMQqAPAmAu3fvhpubG7799lv2z1auXAk3Nzez9tVXX83ocWdzAbX3qXEioRrvnUrD0l1S9sPobErtuKU7Iota2J/hy7dWg4HBppA8CIQibArOw8dnM9iFNVVt/Uiv6URYXiO7E4dAKMLj26MhvF5gs5vxxuA8q+7G0NQzNOHw9BPbo7HhUg4+P59lVgx65MKhke2NYylWOy4+6FfrzILKE9ujJ/xixkVbuTcOh6SVeGqXFKeTa2e8sw1XmpVD096pZibtP1dyZ1zWiGEY5DR0I6a0FcmVHegeHvaNKm4xe+yX/ZPstnc3cSwUAHkQADMzMzFnzhw88sgjYwLgF198gZaWFrbN9EKw9AJSDmqx2PvuTX2xtwSvHU3Bx2cz8OGZDDyzJw4bg/PYIsHrz2XN6nmczcjJ4vM87g7lLpogAI1sD0/jZ2ba5nlE4Fh8FXwjSuFzuxhpVqrH2DWggfuNAnx8NgPXsxVo7VWhUKE0G/ZT6/RYtT8eS3dKcT1bgcy6Lrx1PBVLdkjx7dVcHIuvQmpVh1WOhy9OJdZwFu6+upCNo3FViCtrQ1uvCoruQVS390PeNYiOfjX61To09QwhsqgFuyJK8VlQJt47lTZpsXhHoNEZZrRzzUzbf6/kWq2XO62mEz63i3H/8DSR1X6JbLkfwh8UAF08APb392P+/PmQSqVYuXLlmAA48v9nwxoXUF3HADYF55kFwYnaXHcRL4b6Pgsy9s7sjypHvrwHH5/NmNacvkXeEra2oj1u5Nn1Xciu77bpJH1fcSkWe0uwU1SChIp2tifEmRYGOBKVVo/nDyTY5RoZ3V72T3LZ9+1E/PRq+lnSLqXXI6O2y6xkjyXC85vYx37taAqSKtuhpVXEvEEB0MUD4EcffYSNGzcCGBv4Vq5cid///vf43//9XyxatAju7u4YHJz8g0WtVqO3t5dtCoXCaheQVm9ARWsfxIXNCMmS40pGA47GVeGrC9lssVqBcPJt41wBwzBs6ZuR+5+qtHoUKnqw/U4x/jZqLo9p2Mm0bZS922q/RAzMcl6grLoT689lwjusCPWdA2AYBgNqHYY0euTLe/DkTinmuovw8dkM7JOUI6GinUpdWMA3otTq7/9cd+M1sMBLzG5DOLo9vj0a2ZPU+3RmOr1hwukJtmietwqt1mOXWtVhduyLvCVYsTsWC7zEWLpTik/OZiC5knrYXREFQBcOgFevXsXixYuhUhnnbY0OgAEBAZBIJCgsLMSlS5fwpz/9CW+88cakj+nj4zNm3qC9LqBn9hhLK/zjYCI8bhaibcR8tGblEG7lNiK+rA0ZtV1OvQ1YaXMvBELjsOtEc3NalEN481gKvr2ay67mFAhFWHPSdsWVp2oRhc0zep0d/Wocj682m3s2zyNiWjfSpTulCEqptcbp5pU7BU1m14u12t8OJOD5Awl4bFvUuHsE//NIstX2kHZEVzLq7fq79tLhJNR3Dljt+BXdg/C6VYTHtkVN+Jz7JOW0YMTFUAB00QAol8vxhz/8AQUFBeyfTTXkGxsbCzc3N1RXV0/4M7bsAZzKxTTzD9mVe+PgfqPQbGWbqdzBar9EXM9WQFZ9t4RCdn0XDkkrEZhUY9XFDNa2J7IMAqEILx5KnHAXi9qOAQiEIlS39+NmrgIbLuXg2T1xNpmAPt021V6lVW39CM2S41xqHcLzm/DkTqnFzxlX1maLt8DlMAwDP2kFWyvOFu3DMxNPU1i1L97sC5urseZ+1pO1tQEyZNZ12WwajMHAoKylF4kV7ajtGECevAceN+/uDz3fU4xXj6agqq3fJs9P7IsCoIsGwFu3bsHNzQ333HMP29zc3PCLX/wC99xzD/T6sT1LAwMDcHNzg0Qimfbz2PsCKm3uHffDdq67sZfh/cA0PDiqjtmj26LGrTsXkumYu0UkVbaz9RLfPD7+6lbvsCIIhCKb9ObMpr14KGnC11PS1DtuSQprtAe3ROJAdAXSajqRVNmOuHIKhKPVtPfjtaMpnF8ji30kiCpu4fp02MS5VNsVN1/oFYl/HknG9WwFZ/Ofr2cr2M8kgdD4Rdta8xAJdygAumgA7OvrQ1FRkVl78sknsW7dOhQVFY37b1JSUuDm5mbWazgVri6glKoOeIcVYX9UOSTFLWbDS/KuQfjcLsa7AWmTDmms9ku06zFPl1qnxwMjarJVtfWN+ZkNl3NsdsOZTbuV2zjua8lp6Majk7wHtmgHoit4sVBoukwLiuzRHtoSiYSKdnQNaDCo0WFQo0NtxwB+vJaPeR4ReHKnFPujypHb4FpzAY/F2W4ByGRF8u1JqzcgrrwNj2+PhkAowtO7YhCUUsuWlyHOhwKgiwbA8YwcAq6ursb27duRnZ2Nuro63L59G/fffz+effbZGT2mo19A/WodNgXn4fPzWYgpbTUboop3sN4ird6A2/lNYybRByaZ74OaWddl91A1Vbt/cwR8bhejuEmJ9j41Muu68F1IPmc9lP++nDNhYXG+2R9VbrfFQe+clE14HBqdwWWD+elk25XVcbTtDVuUKrPpG0/vikF5y9gvqcTxOfr92x54GQDlcjmeffZZ/O53v8P/+3//D/PmzcOPP/5o9zqA9vZDaD77wfXeqTRUttr/g0ul1aOkqRexZa24lq1AdHErfCNK8cTwN+uR7aEtkWZDLW29Kocq2Ovo7dm9cfj8fBb8YyqdppiwLTAMg5SqDrNhPFu0+R5i+NwuRnh+E4oalewetq4uubLdZuf0cnrD1AdgZ/KuQRyIKmcX5j3gKTarWECcg7Pdv22BNwHQFpztAtIbGByOqTT7gN0UnGfTsiLdAxqE5TViU3AentsfzxZfnazNEYrw5YUsVIwKqBWtfZyHKmdtj2yNQkiW3GV7oaaje0CDdafT7XbOH9sWhTw5P4LByMUS1mgbLudAVt3p0HUTO/vVZjUlW3tV0OkNkHcNOvRxEyNnu3/bAgVACzjrBZRR22X2YVsg70GvSmvV8jGNPUPYElY07h6rD/tI8NLhJLwfmIZV++Ox0EuMR7dG4auL2QjJkqNrgnk18q5BzoOUs7f3TqWhrsN6JTScjXJIa9eadR+fzeD6JdtFa6/KaudstV+ixftt20tZS6/Z9WT67xW7YxGWN/7cYOIYnPX+bU0UAC3grBdQecvdnrSRPXJz3Y015t47lYYtYUU4m1KLOwVNiC1rRXBmA27mKhBX1oa6joEJew0ZhhkuzZLNPu4/DiZit7gM8eVtaO1VzboXKrW6g/MA5QrtAU8xjsdXO3W9SEt43iocc+3bsrnaoo/xJFRYbxj4REI1lIPOUzcxparDbOHayHY8vprXve6OzFnv39ZEAdACznoB9am0WLk3Dh8Ezmw4bOGIEjN/3hyBVfvj8VlQJn4IzYfP7WK43yjA34aHROZ5RGDNSRlSqzus9gFounFTs0775mouL4vbJlS0Y81JGd47lWaX8/yKf7LLDwk29gxhodf4IUggFOGve+Ontd0l+1njFTllbU1HYtypSInEinb0DGrwrxFfgHeLy7g+PDIOZ71/WxMFQAs48wWk0RlQ1zGA+s4BaPUGqLR6tPWqkNvQjdAsOXzFpfj6YjbWnJThxUNJ+DQoE58GZWC1X+KE33ZN7QFPMfZEllm9h8me87f40t44lgJ5F79qmrX1qeArtv6WcJO1vZIyl+8JautTIamyHZKiZuwedX7fOSmDWqdHXHnbuHUZ57obF30JhMZRiLeOp044FcQZqLR6tqi9QCjCzVwF14dERnHm+7e1UAC0AF8vIIOBQWPPEJIrO3AxrR7H4quwP6oc+yTlEBc2o2fQNh/cIydcU7NeW+wt4dV8JYZhsFNUYvfzvDW8mOuXbjeigmaz1z5yNe/3w9UIFnlL8NGZDGTXd0OnN8BgYNA1oHGpoLw/qpw9B+k1nbydduGI+Hr/HokCoAXoArKtpp4hHIuvwl5JGXaLy+x2o+Zr+yAwHdn1XVy/7Xax/pz9CkSb2kKvSK5ftt0kV96dr/vl+Syzv+voV6NQoXSpoDeR7gGN2YjJ68dSbFp1gUwf3b8pAFrEdAF5hGTQL7WVdQ1obLaFGrXJ2ydnM1x+pfCP1/Ltfl6X+8Zw/bLtRm9g8NQuY93FH0LzuD4cThUoevB/e2LZ6yCtppPrQyKgAAhQALSI6QJ6ZX8U3g9MQ3h+Ey++1drD68e437+Vz22+pxgHospR097vkgsYuOgBfNhHAs9bhfC4Wcjp3rb2Epolh0Aowqp98VwfCudUWj3+eSQZAqEI7wem8WJluKOjAEgB0CKmC+gxzzC8dyoNH57JwL8uZdP+kBbS6Q1m29ZR47Yt9IrEe6dc66b1in8y5+f1REI116fBpvLkPRAIjV8mhjT82BVlMtKSVrPtIc+m1HJ9SLxGAZACoEVMF9D6Uwn424EEPLYtCt+H5uN4fPWYXSzI9PWrdZzfnKmN34Jc5KbVr9bhv1dyOT2X8z3EEBU0I7u+G6lV1iuX5CiGNHp2327/GMfa05criRXtWLLDODQ+x13EuxX4joQCIAVAi5guoHs3hrIf6s8fSMC60+l4N8DYIygpbqGVXzPAMAwKFD2cBx1q47f7N0dg3el0nJfVYcBJdmuYiEqrd6ie5uf2x+OpXVJ8eCbDZcLg1vBi9vU5c1kXa+rsV+MvPxvnBPqKS7k+HIdnMDBQDmoh7xpEaXMvVFrr9CZTAKQAaJHxAuB4bbVfIqKKW1zmQ91WVFo9NlzK4fxGTG167buQfK4vGYv5SSs4P4/jtYZO1+gZ6h7QsK9ptV8iOvvVXB+SQ9hxx1iG6JuruVwfikPo6FcjrrwNgUk1cL9RiC8vZEF4vQCr/RLNhs2/D81Hh5WuIQqAFAAtMt0AaGqfnM1AaTN/L7bJqLR6rDkp4/zGS236baFXJAY1zt0LCACX0us5P5ejW44LzbfcJylnb+LP7o1DZFEL1Dp+zwk0FbXfFMy/FdIDah2KGpUIy2uE8HoBntkTN+Xvwxx3ET4LykRMaStae1VWOQ4KgBQALTLTAGhqK3bHolCh5PrwHU5ZSy/ePJ7K+c2X2vRbeH4T15eNVVxIc5wQ+Ni2KIdbNGEwMBYNvZW19GLF7rulUB7wFOPDMxkQFTTzcopMcGYDBELjsD/f3ClowkKvyGn/PsxxFyG1usPqx0EBkAKgRWYbAE1t/blMmgQ8inJIi1X74jm/CVObXvvyQtbUb6qTGLlrA5ftTLJjLbRp61PhLz/H4pGtUQjJks/6cboHNNgrKcOSHdFmr/fDMxlOP590pr66YNwreLVfIteHwpnuAQ3iytrgHVbEzokc3Vbti4eooNkmz08BkAKgRSwNgAKhCI9vj0ZmHT92X5guRfegWeFUao7b5nuKrTYkwzWGYfBZkP3rA45sC7zEY8KQwcCguEmJokYlbuYqcCCqHL7iUiRVttvlvGy+Wcge3wOeYovrQjIMg8rWPuyTlGOBl3GXjI/OZKBZOWSlI3Zsap2ePZ/H4qu4PhyHwDAMWntVKG3uRVVbHxTdgzbbUtSEAiAFQItYIwAKhMa5VE09/Pjwm66sui7Ow81kbckOKbxuFeFqRgPOy+rgH1OJHXdK8ENoPr44n4V3Tsrw6tEUvOKfjBcPJWG1XyI+OpOBqOIWdPSr8Y+DiZy/Bms19xuFXF8uVtPWp2JLl3DRlu6Umh1Pa68Ka06MPzd2sbcEyiGtTc/HgFqHB7eYD9dZsxTQyC3jluyQusyXicmY5pwu3SnlXc/nbIRkybE1vBjnZXXQ6Ky34xYFQAqAFrFWABQIjbsE+EaUQkKrhQEAGbWOFQCf2B6NT4MycUhaibSaTot7Qb68kMX5a7JWW3c63UrvumNo71PD/UYhJ+dyvqcY2hHbSr53Km3Snz8Sa9v6eilVHWOe8+ldMYgotN6wXFx5G9vj//aJVJcvF/Pq8I4gpx1kqL+hawDLfGOwyFuCxT4SrNgdi9ePpmDF7lg8uCUSK3bH4turuTaZhzcdnwZlsNfehbR6qz0uBUAKgBaxZgAc2aJLWrl+aZzT6g3424EEToPNQq9I/BxZhqJG629cr9LqsS28BIt9JOzzLdkhxRvHUjB/xObxztAuWvFD2ZFIS1ox38P+70VsmfH3v7hJOeXPPr0rxqaLKDJqu7D+3NgvK1+ct+7cz5yGbran8aldUtwpcN1tNR/Zauxh5nqzAL2BQWiWHG8eG7vwbvQ8TVO7IKuz6zGWNPXiiR3ReMg7Eg96iVHUaL3FkxQAKQBaxBYBcLGPBFVttIsIANR2DHA2HPfF+Sy7LNBRafVIruyAvGuQveGZ9lB1huYbUeqyN2oAiC1rxf12Lha9NkAGwDgvauTK2YmarecCGgwM9knK8faJVHjdKkJMaatNQmd5S5/ZArAPAtNR3zlg9efhmikAljRxFzxKmnrx4qEkCIQiLPAU4/PzWYgvb0NRoxLZ9V3IqO1EbkM3ajsGkFrVwe6as8hbgl6VcdrBoEaHhs5B5Mt7kF3fhcy6LmTUdqGoUQl51yA7SiLvGoS0pBWRRc2ILG7GeVkdtoWXYGNwHr44n4UPAtOx5qQMbx1PxWvD02ZW7Ytnd0wRCI3lg8pbrHu+KABSALSILQJgQKJr7w86U009Q9gUkocH7NArtthbgq3hxajt4Pamo9Mb8OzeqWtjcd32RJa5dPgzCeEgkN/Ob8KQRo+57lP/7IuHkqActO1cQHtRafXwk1awveAPbolEqAUrjx3Rq0dTIBByV0KpuEnJjjws9pHgZEL1lCV+OvvV7GfwA55iPLs3bsprc5G3BE/vijH7s9eHX/t02xx3ETZcykF7n/ULiFMApABoEVsEQFsteXd2fSotQrLkWBtg/WLRy3xjcF5Wh34HmpBd2dqH54eHwF87mgL/mEr890ou5nlwv3XZB4HpSK/p5PoU2dU+iX1LxDzgKcat3MZp//xy3xhEFjW7TCCv7xww+11/fHs0Vu2Lx+abhah04n3Wm3qG2NeUUmX/OXUMw+DDM8Y5dW8dT53Rzizx5W14br95ia6FXpFY5huDZ/bE4a/74rFqn3E7w5HTWP68OQIvHU7CW8dT4XmzEOvPZWGnqAQBidW4lF6PW7mNEBU0I7KoGdElrYgta0VaTSdKmnptOh+UAiAFQIvYIgAeiCrn+mU5vMaeIZxOrsU7J2XYKynDtvASfH4+C6v9ErHI2/jNdo67CE/ulOKfR5In7U2b5xGB4ibHLMptMDAoa+k1WxTQ0DmIlRz0Dq72S4T7jULkutAOFTOh1unx+Pbx50XZqr0wPEQ3k/bW8VSXKTJvMDDYFVE65jXOcRdha3ix2e+Fsyhr6cVbx1MxzyMCb59Ixb8uZUNS3Iy2PtuvfjYYGBwYrnU51100q5EOhmFQ3d6P1OqOScv26Ic/u/LlPQ7bO00BkAKgRUwX0KaLqXh4xGR+S9qldNecUD9TDMOgUKGEqKAZBYqeMVtHheU1jvuhyTAMlEPaMTeH1l4VwvObsC28BD63i7FPUo7j8dXIl/fY9HXYSteABkmV7QhKqcVucRk2XM6Z1oKFx7dH4+0TqfjpWgH8pBV471TahEM5c9xFeDcgDbfzXXdC/kxsCy+xe/CeTXvAU8zp/DJr6+hXo6SpF7Flrfji/N0FKbvFZVwf2oz1qrR42d882D+7Nw4PbYmEpLjFJs/ZNaDBzVyFWTmhwKQamzyXM6EASAHQIiMvoF6VFj9Hllk8YdxRA+CQRo+kynb0DU8A7lfrEJBYDe+wIiRWtFu8YELeNYjvQvLx2tEUvHNSNuZDcp6HcRjhp2sF2HA5B49vj0Yaz4Yhp9LYM4Tr2QrsFpdh/blM/P1gAl47moItYUWQlrROOJzSNaDBtWwF/KQV8LldjP9cycWZ5Fq08aAm20z0q3UT1uSzVdsaXoyndkln/O9c+QZ/M1cBgdBYNNsZS8YYDAy7F65/TAVePpzEvp4yCxc6VLX143BMJd4PTMPK4WA58rpY6BWJqxkNVnolzo0CIAVAi4x3AeXLe/Dx2YxpTd4e3R7aEomOGczJsJe6jgGsGp778ezeOPz3Si6eGGc47K/74vHhmQy43yjEnYKmadfK0+oNWH9u5jswzHEX4YfQfKRWGVfRWrNIKCHjGdTo8G7A5LX5rNmKm5So7RjAp0GZWLIjGi8eSsI3V3Px53G+aAqvF+BMci2OxlXZZUiRKwzDsCtYz6XWcX04FtPpDfjncG3AB7dE4kpGw4x73LPquvD6sYkXWKz2S8ReSRkaOmnrURMKgBQALTLZBdSsHMKR2MoZbWkmvF7AwasYq6ylF68fS8H7gWn4/HwWu/DgmT1x6FVpEZ7fNK3X89O1qV9Pa68KH53JmPY5Gtke3RaFV/yTseFSDnwjSnE6uRYpVR1jhosJsSa1Tg+Pm7YvFL3YWzLhl6jIomYs9pbgHwcT4R1WhP1R5S5bj3E8Z1NqIRCKsGp/vEv8vvcMasy+WHwfmj+tENir0sL9RgH77+Z5ROCDwHRcSq9Hek0najsGHGpxmyOhAEgB0CLTuYBUWj3WnU6f1gf+x2cz7DLXakijxwVZHVqU4/cS6PQGbArOG3N88zwi8NO1Aqw5Ob1hMNMcsg/PZOCjMxn455FkfBCYjtZeFdJqOifsyZiqbQsv4cWWUeSuitY+bAkrQqMDbZmYVtMJ34hSeIcVsSu2rdl23Cnh+iU6LOWQlq0T985JmcMuNJgJvYHBqcQa9v2fapFGek2nWZkV4fUCm5RLcUUDah1uZ1ZRAOT6AOxh9+7dcHNzw7fffsv+mUqlwoYNG/C73/0Ov/nNb/Dmm2+itXVmO3CYAuCybeF4Zk8c6ib4hTUYGJyX1WGhV+SkH/gPbYlEqh1KAzy+PRqPbouatPbTlYwGrNgdO6Oh7EPSSrON40e2dafTEVPairjyNjx/IAEPbYnEmhMy7I8qx0uH7873e2RrFEqaeuErNq7+e8BTPGYey2q/RJee40TG0ukNOC+rw4NbIvFBYDpyHGw1MsMwCM9vwgIv69Wr5HqnCEeXUNHOfqZuvuk6+1H//aDxy0RolhwMw6CqrR+S4hakVXciraYTl9Lr8eO1fHx6NgN/3hyBlXvjeFeWyRK3chux2EeCezeGUgDk+gBsLTMzE3PmzMEjjzxiFgC//vpr3HvvvYiNjUV2djaWLVuGFStWzOixR5eBWewjgXdYEW7kKJBY0Y7iJqVZxfyK1j62Cvx47f3AtCkLclpqUHN3c/e//ByL87I6dkNyU+9janUHXj+WAp/bxbiRozBbeTebdiqxhg3Bpm2uRursV+NkQjV+jizD0bgq1HcOoLVXBeWQlj0mjc6A9j41NDoDDAbGpttfEcd1MLoCAqEI92+OwMkEboumD2n0SKhox9mUWpxIqEZYXqNV93ge1NDQ3VRSh/cq/vPmCJc5X0fjqiAQGveFHlnC6sEtYzsQPjqTjiGN8w+B28OQRo8DUeVsp8ZTPuEUALk+AFvq7+/H/PnzIZVKsXLlSjYAKpVK/PKXv8S1a9fYny0rK4ObmxvS0tKm/fjTqQM4evPqqYq7fhqUadO6dHHlbWOec5G3BHHlbVixOxYv+ychprQVP0eWsXP/Pj+fNe0aaA9tiURAYjViSlvZX7Tylj4cj692yW2d+KhvSIv3TsmQXGn/QrZ1HQP46754PLc/Hh+fzeDkiwDDMPAVl457Q7ZmK2+hHsCpMAzDrpJ2laoAQxo93hixoGO+hxgv+yfhywtZ+MvPsfj4bAZ23ClBdEkroSX5wAAAIABJREFUdE5YC9HetHoDQjLlZtsqbricg+4eJQVArg/Alj766CNs3LgRAMwCYGxsLNzc3NDTY14D7r777sPBgwcnfDy1Wm0s+TLcFArFlAFwsY/EbEVeRGHztD78PwhMt0nR3diy1nGf74fQfITlNeLRbVFY7C1BWUuv2XyUqZppVfD9myPw/IEE9ua4wEuMb67mOtTcLTIxg4FBUmU7jsVXsTUTD0krcSS2Ej63i7HmhAwPeIrxsI8ET+2S2rzH2hFFl4z/O2TttiuilOuX6hQ2heRhtV+iwyyiswad3oDUqg4kVrSze++SmdHoDLic3oC//Hw3+C33jYG40LhjDi0CceEAePXqVSxevBgqlTF8jQyAly9fxq9+9asx/2bp0qX46aefJnxMHx8fuLm5jWlT7QTyaVAmO5z56gz2QlyyI9omda6OxFaO+3zvnUqDr7gU8zwiMM8jYsqN6B/2kUCjM6BA0YOAxOpxf8bndjEVEXYi7jemt7p1V0Qp8py0iLalPG/ZfgWwQGhcdCUubMbfDiTgnZMy/PdKLufD3o6ovKUPAqFx0VlcWRvXh0M41qvSIiCxGst87y6QWbJDisCkGrPhcgqALhoA5XI5/vCHP6Cg4O43QmsEwNn0AJraY9ui8FnQ9GrdzfOIgPuNAuwUleCQtNL6JwjAIen4IXAmzbRR9ydnM7BqX/y4P3NBVmeT4yezp9bpkV3fDT9pBU4mVCOisBlNPUNgGMbsGt0tLsOeyDK43yjAT9cK4CsuxbVsBara+ng7B9NgYGZU2snS9uROqdnCkv9eyeX6FDgk4XVjKZQ/b45AQGI1DY3yiE5vgLxrEJLiFrjfKGS3AxUIRXhqlxRnU2rHHamgAOiiAfDWrVtwc3PDPffcwzY3Nzf84he/wD333IOYmJhZDQGPZou9gAVC45wPaUkrIota8PeDCTYp3qnW6c1KCNiqjZ4DSeyLYYy7DhySGncHWLpTijkTrOxeulPKrrhe5hvD9aE7JNOiA67aq0dT0O2Eu1/Ymlqnx3+v5JqdJ1coDUPGKm/pwyFpJd45KcOSHdHjfp49fyABwZkNk05RoQDoogGwr68PRUVFZu3JJ5/EunXrUFRUxC4CuX79OvtvysvLbbIIZLZtrrsI4sJmzHU3fouxRWX/fHkPnts/fs+dtZrpQ7i1V4WY0lacS63Dz5FlUHRTRXpbS6psn7A+3aPborD+XCY2BufhZf+kMVsYHo+nocbxmFZo2rttuJyDwKQa/PtyDvWqT8BgYBCc2cBWWvj6YjbXh0SsKKehGx9OsGnAfA8xXjiUBI+bhUip6pheEW0KgK4ZAMczcggYMJaBue+++xAXF4fs7GwsX74cy5cvn9Fjmi6giJwaq9b/MrWndknxQ2g+lvnGYFu4bYrCDmn0Nt3aarVfIl4/lmL2LW1LWBFvhxDtpV+tw8M+xqGQOe7Gm+GFtHoUKHrQ0a8e8wE5pDEOC4dmyRFX3kbzNifAVQ8gLaKavuTKu+/R2ydScS1bQdtEOrFBjQ4/Xstn39O57iKsP5eFi2n1KGpUor1PPe1tR0eiAMjjAGgqBP3b3/4Wv/71r/HGG2+gpaVlRo858gJ6xT/ZpjeA+Z5i+Ekr0DNoneGfjn41fgjNx+GYSrT1qrDcd+xw8OgCzDNp7wakYcedEraUxXchxl/gNSdkFC7soLVXxZbxWbk3Dil2KDDu6hiGwcdnZ7dtoaXN/UYh/n4wAZtC8lDWwt8b1nQwDDOmHuOK3bGILGqmzx4nI6vuZPd9FgiNW+RZa0oUBUAeBUBbGHkBTXd7NEvbar/ECbdwm4l+tQ7P7InDPI8I1HcOsMe/zDcGX5zPwmq/RCzZMb3af+O1kCw5BtQ6XJDVoWdQA5VWjw2XcyAQinA4xjYLW7igHHLceUZx5W1siJ/nEYGGLqrDOBsGA4OeQQ3y5D2chL/Rbb6HGPsk5VyfFofGMAxKmnpxLL4KS3dK2XP3QWA6TT9xEoFJd8uQPbE92up1HikAUgC0yMgLyE9aYbcbwMM+ElzPVuCCrA7nUutm3SOQUNGOJ7ZH46sL2fh0xOpPS4LfyDbXXYQ9kWXsqtEdd0ogEIrwydkMK78T3NDpDVixOxZ9DlynS941iC8vZGHpzmiccKASIv1qndPUEAzONNYSm+4qfmu38eY9LfASU6/uNA1qdDgQVY75nsZpOo9ti6Kt0xycVm9gv7x+dSHbKp0eo1EApABoEdMF9PWZpHGHUG3Z3jkpQ2RRM3xuFyOyaPpD1x39agRnNuBIbCWOxlUhs64LgPEb84ZLOZM+57H4Kvxtlpvej9xTeJlvDLaEFTl9qYYzybUQCEWQdzlmj4JWb0BYXiP+6Z+MhV6ROJPM7f7JBgMDr1tFeHKnFAu9IrEtvMTh54LqDYzZXtVctImC54rdsVyfHqdS3znATtV5Ynu0y2wd54qKm5Tsda7W2eaLIgVACoAWseUq4MnapfR6fHw2A9KSsfvqTodyUIvTybXYcCkHC7zEeMU/GQkV7ehX6xCQWI21AeMPZ89xn/hmNJO2an88RAXNVn437Eut02PxcL2pZ/fG4WaugutDMtMzqMFqv0T2nG+4lINbuY12PYYWpQo3chTYE1kGj5uF2HA5B0t2RMM/ptKhe01Hqmnv5zT8mYLeeH++dKeU69PjdIY0evZ8bgrOozmBDqhPpWWrF3wWlGmz56EASAHQIlwEwAVeYii6B7FidyzmuouQVNlu0WuQdw3ibEotHtsWhb2SMrZHpq5jANKSVviKS2d9rPM97q6M/vhsBvvBez3bscLSbN3MVbAh8PPzWVwfjpktYUXsud9+p4Tt6bUXjc5gtgWTqS3ZIXWqm26zcojzADhRe+ekjOvT45Rk1Z1s2aOD0RVcHw4Z5Vi8sdTSU7ukNhn6NaEASAHQIlz1AJrKewiEIrx5PBWAcajKklIHuQ3d+C4kH+43zPfTNMxyCKyoUYnuAQ1bYia+vA3Vw70py3xjXKasRfeABhfS6m36QTUbP4TeLZvwyNYovHAoCTtFJWhWDo0JYCGZcvjHVCK6pNVqxXNDs+R463gqHt0Whbnu5tes160iqzyHPSiHtJwHvYnaxuA8rk+P07ogq2PP4z5JObROPB2lZ1CDn64VYMkOKR4bru/pzHMcPwhMh0Aosvm2hxQAKQBaxN4BcJ5HBF46nITj8dV42f9uKLud34Tw/CbsFht78Gra+5FS1YH9UeXwj6lEanUH8uQ9+OlaAXIauid9TfKuQfwcWYbd4jL4x1RiUKPDRxMU35ysLd0pRb9aB++wIjy0JRIhmcZVwd9ezYX7jQLsuFOC6JJWh58D5qz6VFrsFpeNu6Dn1SPJZvMv15y4O+QfktmAAoX19vg1GBho9QbkyXvw2LYoLPONQUmT83zg6vSGCXdO4bqtP2e74TE+GLkn+kuHk6x63duLckiLZ/fGjXt9bApxviHufrWOHVU5l1pn0+eiAEgB0CL2DoC7IkrR2a+GqKAZj22LGvdn5nuK8dCWSCRXdqCzX40LsroxHxDPH0hAeH4T+tW6CT8gWpQq/HStAB+dycDt/KYZH+sPofnQ6g1Y7ZeItQEytPY6Vg8ZX+iGw1dIptxsAc9bx1PZXo8LafXYFJxntoemLXoQDAZmVgVbTVKrO3AsvgpBKbV27bF5apd0xtc/BUDncCu3kd05RCA01i+NLXOOL6YMw+A/w9vfLfeNQVJlO/LlPfC8VcgOcc9kgSDXFN2DeO1oCgRC4/Bve5/aps9HAZACoEXsHQAXeUuwKSQPaTWdOJdaN+bvEyrakdvQjdoO83pvegOD7PouHI+vxvuB5rt+LPaRjPl5k/rOAQiEItzKbURKVQcCk2omnJA+us11NxYgNv3/nzdH4MsLWVi5N87m3+zIxEYu8Bnd49HRrzZ7Dx2t9yAwqQarhrcu3BNZZrfnfet4Kudhb7zmTEPpjqytT4WNwXlm2yGuOSFD3QSfi45ApzfA53YxBEIR7t8cgex68zm+/7qUDYFQhP1RjlcvkmEYyKo7EZIlZ6ekiAub2SC+2EeC3ClGqqyBAiAFQItwNQdwdJvjPrMbYmVrH94PTENmXRf0BgZtfSqkVnXgdn4TgjMb2Bt/z6CGfY5/XcpGVHELntolxatHUyweFnt2bxy8bhUhlWqZ2U1Ilpwtx/PioaRxe9HiytrY9yi6pBVtfSqotHqHCIOZdV1Y6GWsDWarrRHH4z1iQY0jNT8pLWCwpqaeIfhGlOLB4fpzC70ixwQrRzCk0eP/s3fe0U1c6fv3nu9J9pxfNmfPZpN/tkRAIA1II2wIKYSQXkjZsGlsCqQQspuETZFcwICNCcWG0MH03jEg2XLvXe5d7pKbXOUiS1aZ5/fHWGPLlm1JM+r3c877B7YYzUjXc5+5932fd+nukc5T57Ibx73muW0jPd4Xb09ESZPSCWc6HuWglsnxM8bTW0YWFd7cneowWy0iAIkAZIWzBOCXJ3IhuFLETIZrI0psnqBPZtQzYu69/em4LyASP14sxOAQPemP3h4xRmhMFfIbuyfchrYm5gSKOf5WCOaQdamYFY6Vx3MnLPbQGyhmAhwdM/3oFdzlh7MQLCzDrjjaR3J3vBQRBU3oHuCmReFkJFYqmGuIKHCcpU2KtN3pYs9ckG4g9kHWpWLyYhdtTUBVW59LPAAZKZKPdKSZyNpp6fBWqjFcZSVwe3QleHx6R+he/xGXiAfWRmGruMJunn/mIAKQCEBWOEsAzg+ORVVbH24WNTPu9j9eLLTpGgwGyqRC0xgv70hGZ79mwgTj78/lm+1QMFncFxCJa/m0WEiqasfSPWlYfTqP42+FMJrOfg32JFQz38HCzfFT5jcVyXuw+nQe5gaKTVZ6l+6ZuN/148Gx6Oi3b84OACYfdfnhLOxJqHbIygZFUXhnb5pVY90R8culoqlPnmATykEtntg0Yu6/cHM8QmOqnGL63twziKiSFpzJasQliRxNPYPMSv5EhSud/RpczZcz559e4xo7LZ8N99Lel1iDIZ0BSVXtuJgrc8gD5FiIACQCkBXO3AKe5R9pkrMyXUDbQrQop7ZXoSgKA5oRF/zj6fWMXcvoMK4wzlknxoGkGrMi8XBq3YQ2MR+FZ+JIah1WHjdtzP7uvnQ8ujEGL4YlQdFHikO4ZHBID7VWj6iSVoRGV46rArZWMBkMFHrVWjT1DELS0IUTGfUIFpaBf7kI/MtF+PlSISMS9yXav9Xc2DzFZQ7ywpN1qcz+jTgzPjiY6ZBr91aaewbx6dFs3DPmPvvNacmUbgpcoNMb8L8LhRN+/9MFwkntp+TdKvD49Oq9tStr/RodkqraEVPWxqk4418uAo8vxFax43J4J4IIQCIAWeEqOYCjY26geMpm561KNeYFxeCni4WoVvQxP99wo8zsMWf4ijAvKIbpzcjjC/FiGF1ROnudGI8Hx2JzZMW4vMCXwpKZFcqJYpOo3N5fk1eg0xvw/bl8s5/xktAkhKfUoqHTPkntxkpEHp82J86osZ8HGUVReGzjiKidLhA6ZDI2Im3rw6HkWgQLy/DTxUJ8eCjTZKXI0X/rrrQ16akMDulxvbB5XAHdktAkrDlfgPCUWhxPr8eO2CoE3SxDeEotUqUdVlcSG4v1jB6pN0a5L7yxKxUrj+cyRVD3+IpwIqN+0uMlViqYe7Wl5DV2w/dqMe4PGLlv29pxaiwURTFbwGsuON/DkghAIgBZ4YoCkMcX4j9n86c890+Hl+Jf2ZmCs9mN+OBgpknnjrGx7ECGyUrf48Ej1hj3+IpszgdcuDkeEQVNbt8X2NmEp9SafK6PbIjGz5cKcTa7EYND9s2ryaztHGeVcji1zm7v98P5AuZ9Zq8Tu0QbvtF5Y46MiSr4CfahsrUPP10snPReaYyPwjOxI7YKnx3NxiZROc7nNJrsvBjpGhhC4PVSJvd2hq8Iv8VJcSiZ/pt+/+DIKjdFUZB1qdBlwaqcsVf51yclU762rLkX746qdn8pLBnPbEnAKztTONmlERa1mOT87YqTsj4mW4gAJAKQFa4oAGf5RSLTAg+3zn6NyYqeJfHMFvP5gMaYs258LqGl8fz2RCgH3aM/rCvy+XCP5hBRORS9alZdYWyluWeQ2eKZLhBOuRJtKxk1nZguEOKtPWlQax2XND4VTT2DVv9NsY0LuTJnX7ZX0jUwhLjyNvwWJ8WqUxJ8c1oC36vFCBaW4ZvTEtwXYF4gfnEil1m1lTR04YODmSapPObC1pZ/Adfo6vWpHCKSq9ox00/EzB8/XSxEkbyH09Xl0ZX0yw5kuMS9nghAIgBZ4YoC8JEN0RZvPYTFVFl17FWnJON+9unRbAReL8WqUxL0a3RT3syMYawuniYQIjylFpm1nW7djsnZCK7QwuuBtVFo7HR8oroRiqKwbH8G/rEpljEDtweVrX3oUTk+cXwq0qo7MD/YccbRpB2ca1LR2gv+5SJ8djQbAddKsPzwiPXJd+fysTaixMTK5fVdKQhPqUVtez8u5MqYFcbpAtvNnM/nyPDKzmSsjTD1i6QoCnmN3Qi6WYZPjmQzKRWfH8uxm2G/cTXyiU1xLpO2QAQgEYCscEUByOMLsX+4h2Jtez/e3ps2YeVar1rLFHpYGrP8I/HtmTzm3+/tT8fjwbH47lw+LkvkJq+9d5L8vx8vFmJBCJ07teZCAd7dl46Pw7Ow/HAW/rkvHa/sTMGy/RnIrXc9Hy5XpKa9n1l9cnaCdWVrL7Od5X+t2OuEfY9qCBtvlk2Z/8omZvlF4vHgWMwLinWZCZUwOUfT6kzypJ/ZEo/ntiWafWCrVvRhd7wUkgbb81svDd+Pl4QmMWOkUNaDN3ebVvM/ujEGv8VJ7bqabjStdqXKdSIAiQBkhasKwHt8RUioVOBmEZ1E/GRIHMpbRga5Tm9gbgij86m4iic2xWHN+QIsHmVGOjamCYT434VCXJbImfY/xoltfnAs3j+YAUWfGhRFoVrRz7ot0JDOgM2RFQgRlWNXnBRnshoRU9aGaoVjPL7aetVTVmh3Dwwh4FoJ3t6bxiSCW8OVPDne2ZOGr0/mOj2n0njDN8aGG2VeJwRblWr4XS2e9EHIkvjHplh8eSIXexKqEV/RhobOAci6VAgWluHDQ5ms2usRHEtGTSe+OS3BVydzcTKjAX1q+22F9qq1zEPIgaQaCK4UM/Yx9/pHYs2FApzNbkR2XZfd74FGMTo3UGxSeOhMiAAkApAVrioAeXx662C0w/rTW+Kh1upxJU+OaQIh9iRUAwB+uVRk0fH+czYfIaJys797a0+aidWLJdvAS0KT8MyWBMZiYfG2RIiKW5BQqYDgSjH+ezYfzT2D2JtYzQjGFcdyEFdOT4DRpdZtixTLlROey6s7U3AlT26XvLmOfg1+vDhi5fDtmTzkNXaPu+Fuiaow6cV7cooKP3OM9gd77bcUdDrAl28ihnQGbLxpWlUe6iJmtI6mT61FREET1t8oxXv7080abRtjbqAYH4dnYau4AuLS1kltPgiEqRj7N8jj0x6ujvDsHI1Ob2DMqR9YG4UNN8rs3ut3KogAJAKQFa4sAM1Fi3KQsayY4SuCuLTVxCTYXrH8cBZTpGBJfHsmDzN8RVPmUllTZDC6xdlEMT84FluiKlAo4yYB+mZRM+aYMdnm8emil2BhGc5kNWLNmFXYt/em2SzeLuTKmK3gYKHj2qVNxJmsRua6Pj+W43WrgObQGyhUK/ogKm6BsKgFUSWtiCtvQ217v9XWIQTCZKi1eoREluPlHclYcSwH2XXOS6np7NfgXwdGKuUf2xiD3PouGAwUVEM6JFQocCarEedzGpFYqbB7IRsRgEQAssLdBOBYU9EZviIIrhRZZGkwVTy3LXHCyjdbc6Fm+UfiyZCJPdasETg/XZzYUNVcfHNawuoGdC670aZ+yT9dLGSdi2O8yQZeL2V1HC7Ib+w2MdJ9ZWcKpG2usQVEIBAci8FAIbFSYdKreJZfpNl75cLN8XbtvEIEIBGArHA3ATg2QqMrzfb6tSXe2JU6oQC0Np7fTm8HD2h0jJWBuZgTKEZMWduEq3W9ai0aOgeg1RvG+dRZEj9fKrR6JXBIZ8CuOKlN1/35sRzWK4+yLhV+uliI+cGxqHGRXJsieQ98rxYzY+2pX+PtmvtEIBBcG2OXldFzxlO/xmPl8Rx8djSbSZN4fnsiEioVdjkHIgCJAGSFuwvAk5kN45L1nR3PbUtk8lPKW3oZf6rJ4uUdydgmrkR+Yzeq2vpwKrMBX53MZW4uCyZZRZwq/K8VT2qkrDdQUGv1qGrrw29xUhODbGvj6S3xrLdIlw2v/r2xK5XVcexBe58GCzfTeamhMVXOPh0CgeBkdHoD5N0qtPdpTB5+qxV9TBvL6QKhXXwDiQAkApAV7i4AXTHe2ZsGAKhq6zPJF3FmzA+ORYioHAkVClzNl+OH8wVYeTwHL4QmcW71cb2wmdWYnDucc/hCaBJKmpQwGChQFAVFnxoFsh7UtPej30w3AkdhbG/1ZEic086BQCC4Psb89Ht8RXbZMSACkAhAVhAByH18HJ6FSxK5XT3UXDluFrETgBdyZSY5d3MDxWarTucHx+I/Z/NxOqvBbh07zNGv0THn4IpGzgQCwXJq2vux4UYZUqUdnB/71Z1069ENN+xTzEYEIBGArCACkH08v31ir0Bvi69O5nIyLstberH6dB7uC4jESzuSwOPT2yhPbIoz265vmoDu8uKIvrLxFW3g8WkfMo3Oddq4EQgE6/n3kWzmPnK9sJlTM+n39tO9iVedmrqXsS0QAUgEICs8VQCyNa61NGavE8PvarHTr9eeMV1At+ez5LXncxo5HZ86vQEtykHUtPebiK1+jQ5ZtZ3YEVuFZftHttmf2ZIA1ZB9t4d9h7/vf+5Lt+v7EAgE26EoCpKGLsSUtSFV2gFZl8qsRdGirab94e/1j8Rbe9LwvwuFUxrfT8X57EbcFxAJwZUiu1jCEAHowQJw3759mDt3Lm6//XbcfvvtWLBgASIjI5nfL1q0CD4+Pibx9ddfW/UenigA39ydilOZDfjiRC6W7knDiYx6s6tGXMTyw1k4lFzr9Gu2Rzy9JR6746UIi6my+POLr2jj+s/AIsqae/HwsEh9e2+aXVcCjT1RSREIgeCa6A0UVh4f79s6N1CM/5zNx4mMekQUNGFtxIhDw9cnJeMsuxZP0ObOUkKjK8Hj03nh9uhUQgSgBwvAGzduQCQSQSqVoqqqCn5+frjllltQWkp7oy1atAhffvklWltbmbB2IHiiAJQ0mBqFdvRrTJJx/3UgAz9dLMTqM3kmzcxticc2xmCFFQbRrhIPTtDJYfG2RPhfK0ZkcQsSKhR4MSzJquM+vCGa9VOzrUgaupkCkln+kTiZ2WCX91k3PGl8eCiT9LAlEFyQnPou5p70xq5ULAlNmtQrds2FAgD0qqG0rQ/ncxoZy6d7/SOx/HAWvjuXj2pFv1Xn0a/RYcONMnQP2CdXmAhADxaA5vjTn/6Ew4cPA6AF4Pfff8/qeJ4oAI1FCGMn585+zTiLEoqiUNKkRHp1B85mN+Jwah2Sq9rxzt6R3r7TBfSqkrOvi8tQ9Kmh6FOjWK5EqrQDGTWdTMsu1ZAOX57ItfnYbKuA2dDYqWJW6Hh8IT44mMl5np6sS8VMJvsSazg9NoFAYIesS8W0bFtxLIf5uU5vgKShG9ujK7HyeA7eP5iBny8VIqFCYfZBTtalwgcHM03ubc9sSeA0R5AtRAB6iQDU6/U4d+4cbr31VpSV0RVFixYtwp133ok///nPmD17NgQCAVQq65arPVEAThcI8dD6aMzwFeHpLfH45Eg2/K8V40BSDVKlHRY1nl95PAebIyuQXt2BakWfy3kN2vKZ/HSxEP7XirE/aXLRUtI0cc9hS+JMFrd5gNZCURS+P5fPnM+VPDnn73EyswE8Pl18klPvvNZUBAJhhNNZDUz+99xAMWrbrVuxGwtFUchr7EZ4Si3uD6B3Td7bn47mHufscoyFCEAPF4DFxcW47bbb8H//93/44x//CJFIxPzu4MGDEIvFKC4uxunTp/HXv/4V77zzzqTH02g06O3tZUIul3ucAJwq5gSKISxqGffZUBSF7oEhlDX3IlXagRRpOyQNXUisVLAyYnZ2bLhRhvzGbqvGXWhMlc3vZ8/WR9ZgFO2PbYxhPRGYw7g68OnRbM6PTSAQrCOxcqRX+vsHMzhv15hd18WkzswPjoWiT83p8W2BCEAPF4BDQ0Oorq6GRCKBQCDAnXfeyawAjiU+Ph4+Pj6oqZl4hScwMHBc4Yi3CUAen+68seZCAVYez8WyAxl4fnsi84TnabE7XmrRWLuSJ8dTv8bj2a0JrFriNbnI07GiT43FwxY97+1P57wK73h6PTMZEAgE57J9uOCCxxfarU1jXccAZg8XxB1JrbPLe1gDEYAeLgDHsmTJEnz11VdmfzcwMAAfHx+IxeIJ/z9ZAfS+mLNOjPY+zYRjomtgCCGics7ez9zqqrNo6BxgDLlf2ZmCsmbubpQtykHGsDqmzDnVzwQCgd69eWYLbecy00/E5DPbg6CbZeDxhVi8PRHlLc4VXkQAepkAXLx4MT799FOzv0tLS4OPjw+KioosPp4n5gCSGB+PbowB/3IRzuc04lRmA7oHhvDVyVwsCU3CNAG37/XNafuYntpKUlU7U9H3yIZodI2pyFNr9cio6USKtB0DVraY+zWqAjw+nRxuSW4pgeAJ6A108Vx5Sy/CU2rR1uvc7dA+tRbv7U/HTD8REisVdn2vtl61Sb/05YezUO8AA3pzEAHowQJQIBAgOTkZ9fX1KC4uhkAgwO9+9zvExMSgpqYGGzduhEQiQX19Pa5fv44ZM2bg2Wefteo9iAD0zjA+LdsjBFeK7fQXYTsd/Ro8upFuzP7Z0Wx8FJ6J9/bxAi5pAAAgAElEQVSn4/VdKZgTOOJx+ND6aCRXtVt8XNWQjvEfFJe22vEKCATn0qIcRHhKLVYez2Hslt4/SJuw/3Sx0KnnpuhVg8cXYhpfaHcjeABo7hnEN6clmDG8AzBnnRjp1dy3kpsKIgA9WACuWLECPB4Pt956K+666y4sWbIEMTExAACZTIZnn30Wd9xxB37/+99j5syZ+Pnnn4kPIAmnxQxfkUuLoLCYKqY109j4x6ZYPLEpDjMEQry+i+7f+fCGaIuKZzbcoLeE/K66nvAlENgibevDZ0ezMX3MTsHsdWK8sSsVr+9KcXrlPwA8sYku1HNkVX5D5wDe3TdyT1l1SoLkqnbo9Nx3/TAHEYAeLAAdARGAJLiK8JRaZw/nSTEYKKTVdOCH8wXYcKMMV/PluCyRo0DWA4OBwpDOAP8xbf34l03TKZQqLfYkVJv8LLK4BTy+EAs3x9st+ZxAcDQURSGioAkz/UTM38OyAxk4mFyDInmPw0SOJRgMFF7ZST+4RRQ0OfS91Vo9frlUZJJK8+jGGGyJqkBMWZtd+4UTAUgEICuIACTBVbiCLQJbNDo9LkvkzDV9djTbRNR19msww1eEzNpO5meqIR3TK/mNXSl4ITQJay4UoFDW44xLIBA4IUXajme2JGD2OjFWHMuxa3tFW9DpDUiVdmB3vBSCK0W4b7jYq7RZ6ZTzqWjthd/VYibVxNhF6fntiaji2JLGCBGARACygghAElzE4m2JHtUWTdLQhfnBsZgx/FRvTCynKAoPro3CTD8RLuTKmNfn1nfhmS0J+Ch8pAvJPb4iFMudMxkRCGxoUQ4yubHv7U93ub9t5aAWrw6v+Bnj9d9cYyuaFqbtWBdRgnlBtBicFxRjF39UIgCJAGQFEYAkuIisUStinsLZ7AaTa0yvoZO8V5/OY342uiMARVFo6h5EYqWCmTxdwSuMQLCUFuUgruTJTapcXW0Mx1e04YFhQ+bpAiF+OF+AjTfLUNLkeg9bXQNDjFB9ZksC5yKQCEAiAFlBBCAJtvHBwUxnD2O7oDdQJv6IgddLAQA9qiF8cDAT84JioRwcn/NnMFCMWayxLzWB4Mrk1ncxFb3GeHlHst22Lm2lrmPAJCdxV5xlJvfOpFWpZlwXHt0YgwQObWqIACQCkBVEAJJgG0nDtimxZW1YF1ECg4f54e1NrGau9auTuRBcKUJceRsaO1Vmt8aMRSGz/CJRreC+BR2BwBUURWFzZAUzvqcL6Ny1HbFVVntiOgL+5SLw+EK8ujMFEQVNnHf3sRdtvWq89hu9Ejg2h5gNRAASAcgKIgBJsIm396aBoijIu1VYtJV+yr2QI5t64LkRegPFuP+PjQ8OZsL3arGJEey+xBrw+LRYJBBcFVmXCl+eyGXG8i+XikxSGlyNCzkyxnfPkVYvAJBe3YHwlFqkVXdgcMi2ql6NTo/lh+kc4ee2JULSwP4aiAAkApAVRACSsDWmCYRMZeDnx3Lw0PporDye41L2EFySU9+FTaJyfByehXlBsZjpJ8Ky/fS22Sy/SFzLb4LBQCG9poP5jFKljjeHJRCmoq1Xjfmj8vxc3cKpUNbDnOt35/Lt9j4GA4WYsjbsS6xBYyedr1et6DO57z20Phq746U2rT42dA4wXYl4fPoBMkVK+wba0kmICEAiAFlBBCAJW2P54Sx6DKm1jAdWbbv3bHlqdHqkSNvxrwMjuVOLtiagtEnJmMM+tD7arn1JCQRraVMO4o1dqeDxaYsSrrYj7cmh5Frw+EK89lsKtHZ4wFSqtONyIO/xFeGtPWlYEppk9v73zt60SXusT4S8W4VfLhWZ5DJOF9AP1J8cyUZlq+V5l0QAEgHICiIASdgaxi4ZBcNP5//YFOvk0ex4BjQ6qLV6fHUyl7mhP7EpDopeNZ7fnsj8O6KgifQKJjiVFuUg1kWU4O09aeDx6b7Yzuphaw1avQErj+eAx7dftx1j2oYxzPVHv9c/EnUdA7iaL2eq/D8/lmOzRU5TD/19zB3VipLHF+L+gChcL7SseIwIQCIAWeGqAnBOoBj3+Iqcfh4kzMeLYUnMGMpr7AaPL8TTW+KdOJKdw/zgWGYS7R4YwnPbaNG34UYZ6joGsHhYBPL4QiwIiUPAtRIiBAkORzFmy/e/Z/MgdbEKX3M09wwyhso8vhDZddzn/tV1DOCxYfNm36vF6OzXgKIoVCv6ISxqwc2iZiRVtaOjf2S1r7K1j3ngY7vroTfQ71XarGRyBHl8Ic7nTO1pSAQgEYCscDUBOMs/EvsSa6DTG9Cr1qKhcwAtykG0KtWILm3FZ0eznX6OJEwLPVKlHXh4fTS+OJHjxJHseIrk9Mrn6FWAFGk7eHwh5gaKoRzUol+jw644qclT/oNro1w+54rgOVAUxfTAXrw9EWnVHS5n7GwOiqKYVIrHNsaYGK9zdfyokhbGrHnR1gSrtnRf3pEMHl+INRcKUN7Sy8lnqjdQCLxeSs+FfpHIrJk8h5gIQCIAWeFsAfjm7lREFDShvKUXla196LfAeqBFOYjrhc348WLhuOVzEvaPhzdEQ60dqYSLLGnBC6FJ+O/ZPHsOVZcju66L+Uy+PZOHwSE9KIrCk5vjcH9AlEkXkMEhPQKulZh8jueyG9GjGgJAt5jbJq7EyuO5OJJaZ9ZfkECwhdiyNvD4QjywNsrl2rlNRlOPCm/sSsUsPxHnBsrybhVWHMth/hZf2ZlidT7flTy5yVbxS2HJuJInZy0EKYrCVydz8cauVDweFDPOSkqt1aOpZxCyLhVkbZ1EADr7BNwZZwnAlcdzkVPfxfqPRa3VY3t0pdNFkTfFuogSk+8gPIVO0LZndZ4rQlEUTmU24LltiXg8OBa9ai3UWj3zOZmbUJSDWiwdzsEyxhOb4jDLzzTdYYavCK/sTMH6G6VeVVhD4BaKovDW8HgLiSx39ulYRURBE3h8utiCKzr6NdgkKsf9AXQnkZl+ImyPrjR5oLWGhEoFVh7PwX0Bkczf7qdHs9HYyU5oK1VDzArjIxtiECwswxcnck228Xl8If7+w0UiAJ19Au6MowXg4u2JnFtj1HUMOF0UeVMYjZ+NbImijWSNnTK8DYOBgmpoZOX6H5vom3RkcYvZ1w/pDNiXWGNSXTjLLxL3+kdi/Y1Sk5wnHp9OCneHZH2Ca6Ea0iFYSPtX3usfaVPFqjMJjakCjy/EqlMSVsfR6Q1Ir+6A39ViE6G2bH8GZ3mQykEtdsdLMcs/Ei/vSMa9/iJsE1eiT237Sn7XAC0CHxxuezc6ZvlF4r6ASCIAQQQgKxwpAFccy7GLu3xVWx9m+UXa/fxJ0BPJ2KflY+n1eDIkDoIr9qnQcze2iekV6U+OZE/5WuWgFteHVzp4fCFyhw1uW5VqiIpbmJZyp7Ma7H3aBA8ivqLNJD3G3caPclCLeUH0g1REQZNNx6AoCpcl8nGrZkv3pCGhQmGXPEhpWx/+ezafeS+2htXdA0Mokvdg9ek8HEquhaShGz2qIebcyRYwEYCscJQAFFwpsmv1o7StD/7XirHiWA5e3ZnCVHWR4DY+PTpe1Gy8WYbXf0vB6Uz3mmTsRUmTEjy+EHPWiS2eZIxegvzLRSY/N64mbo+utMepEjyQ9JoOTB/OTZsfHIuDyTXOPiWrEVyhW749ty0RGp3127M6vQErj490OXlkQzR+vlSIVKn9C2AoioK4tBUbb5bZ9X0AUgQCEAHICkcIwDd3pzqlO0SLchC/XCoy6+lEwrYwt61p9Og6nl7v8O/YFRnQ6HCvP70iHV/RZtH/MW53jc2jXH06j/nsz2VPbQtB8G6aewaZh4b/ns13m165oymW92Da8JjPstGk2ti68V7/SOxPqnHLz8ESiAAkApAVjhCANyw0tbQXnf0aVLb2oVDWA0lDF+o7BlDR2ovLEjme3hLvdFHlLvHEpjizQv7z4Wq6s0SgMGwcnoA+Cs+06PWbROXg8YV4ZkuCyWes1uqxLoKuHp4moKuNuwaG7HXaBDemrLkXCzfT97Pntyea5KW6CxRF4e29dPeNDTeszymmKAp7E6qxcHPchA+srsTgkB4XcmXYm1htU8cgIgCJAGSFIwSgotd1W2FVK/qY7RISk8dExqRfDDeU3xUndfC357pIGmiLmEc3xli05SRt62Ma3SdUKEx+p9MbTNrNhYjcq5qTYF8oisKehGrGOP+5bYlo7hl09mnZRFo13Uf7Xv9IKPosnzc6+zW4mi/Hl8P3okc2RONUZr39TpQDCmU9WLQ1gfm7fmxjDKJKWpBT32WRHRpABCBABCArHCEAa1zcxsKYb0Ji4vjwUKZZIaPTGzAnUIxZfiJUuUFnAUfR3qdhCpNCIsstSoH4+VIheHwhwmKqxv1OpzcwBrEPrY922wmewD1F8h68uZvu7fv1SYlJxwp3w9iP11JHge6BIfhdLTbpqztNIER4Sq1LmV33qIaQUKlAbFkb0qo78NXJkfzEf2yKZVZujXGPrwj8y0Vo7Jzc/5AIQCIAWeEIARhbZlkelLPoHhjC4u2JuD8gCm/vTUPg9VJcy29CXmM3LubKEHi9FMv2Z5gtx3dGPLctEfOCYvHjxUIcTatj3PKNMXtdFB7g8Fwf2RA94Y3IuNI1yz/SLk3a3RljA3seX4ilu1Phe7UYx9LqJhRvxtevPG6+o4pOb8Dru1LA49PGtYNDtnmXETwL4wOsu/tw1rT3g8envflalFM/4Egauk0qfF/7LQXrIkpQ0eoYMSTrUkGjpc3flSotcuq7xt0D1Vo91t8onXCXafG2RPSqtWjrVeP7c/lYtDUBC0LiTF7zxKY4/OtABiIKmsaJWiIAiQBkhSME4H/OuseNaaonRoOBQluvGjn1XTiX3Yi1ESV4Z28aYypqz3ghNAlPbo7DvsRqhEZXYne81GRrXdrWh52xVXhgbRRmCET47lw+05eWTTweHGvS0WIsJzPqGVFKGM/1wma8vCMZD2+IZj7TWf6RCLxeisRKhcnKYGkzXT08yy8S1Qrzq6nybhXTuupoWp2jLoPgolAUhaXDq3/uXoR1PqcRPL4Qyw5kTPnazn4N5gzb3Dy/PRGZNhaL2Eptez9m+UXis6PZ+Oe+dOyKl4LHF+KpX+NRIOtBqrQD7+1PN7HiWbwtEW/uTsXCzfF4Y1cq9iXWTNjxJ7e+Cx+HZzFpIcbYl2ha0U0EIBGArHCUDUxipWLqk3FT9AYKlyRypvqOq5jpJ8J0gRBbxRV4dNjWxpyfF0VRiClrw/PbacFnFAiLtydifnAs7vUX2fT+7x/MmHKrMUhYhhfDkrBJSPLSJkKnN+CSRI5t4kpmBY/HF+Lfh7Pw3LZE7EmoRnxFG1qVg/j3EbrX9T82xU5oUntiWHS/uTvVwVdC4AIu2/ydzmpgHipsKSJwJdYOFzttsiDHdX9SDfNgzFWxS0PnAC7kyhBb1jbl6nqhrMfi++iTIXHj8notRTmoRaGsh0n/4PGF+N+FQsbUmwhAIgBZYRxAS0Nj7Gqm/PVJdm7u7oBqSIfTWQ3YFSfF2ogSm4tL5gSKcSqzwcRw2ejoP8svEtfyR4xRDQYK/MsjOYyPbIhGdVsfSpqUeDEsiXlKfiLEVJw+sDYKr+ykxch0AT2BzPKPxIthSRBcKUKRvMeia1Zr9VgenmWzXYO3oTdQiChogu+VYjy0Pnrcd//SjmRGwK84Zn4ruEhOTz4LNsU6+OwJbEmVduDRjTG4mi9nfSy9gWJyx8JTajk4O+dirJyfqmVdU88g89o9CdWs37epZxDfnskzsQt7ZEM0wlNqUdvejwGNDi3KQVQr+lHb3g9ZlwrKQS0EV4qZB/Xlh7OQIm3Hd+fycV9AJGb4irD6TB6yajs5SY2hKAo7Y6XM+c1eJ8bPlwoRcDGbCEBnn4A7M/oJQq3VI6e+C9vElXghNMliwWJJmDMQ9nSqFf0Ql7biTFYjdsVJEXCtBM+OqvqaSCibq5pWDmrxcXgWIwKNT4AhkeWMiPO9WoyGUT0oFb1q/HC+ADy+EDtjq3AhVwZhUQvk3Spo9Qbo9AZIGrpsMlodTWOniuT/2UBHvwbH0+ux+kweXgpLHudXSW/51+C/Z/Ox/HAWPjmSjRXHcpgeobP8RCQP0M0IHe5bfl8A+9ZskoZu8PhCzA0U29zL1pXYPbyN+suloglfQ1EUU+n7QmgSq+YCHf0a7IitMknheXdfOp4ck4NnSUgaRjp+aPWGKb8PW1vESRq68MauVOZ97/npMhGAzj4Bd2ayJeSa9n7sipMyq0hsQuTifkyOgqIoSBq6TAx+pwmE+OxoNrLrJm8bZDBQWDyc1xdT1ob2Pg2zyjh6VXA0xm0V0knC9eno1+BIah3mB8dO+QC2YFMcHlwbhVWnJG7p9+at6A0Us/I+ka2SpRi3Qb84kcvR2TkX3ytFU3YUOpRci38dyMBMPxEKZZbtUpgjVdoxri9wWTM9B+r0BpzJasRbe9KY19zjK8JD66MxN1Bs8v94fCGWH86y+CG6T02vHD71a7zFVi+lzUrkNXYzOeoGA4WoklasjShBeHwJEYDOPgF3xpIcgq6BIZzNbjTxIrMmVh7PdamSfFehoXMACZUKq1YCjHYPJzPqmZZj84JiJnz9+ht07sivURVcnDLBAai1ehxJrcM/96XjkyPZTG/Rx4Nj8P7BDAQLy/D9uXwmQfylsGTWq7gEx2Hs+rLmfIHNx6AoisknPeYhxUAvhdEr2xN1vMmt72LG/MVcmc3vUyxXMiJuQUgcrhc2Tzg/6Q0UBjS6cb9Xa/VQ9KrN/m4iBjQ6ZgWYx7e8QcLexGrw+EKsPp037nckB5AIQFZYM4AoikJ4Si1jOGpJ/OtABtmm4pA1F+gt3a3iCrT1qpkVRI1OD4OBQl5jN8JiqvDD+QJ8FJ7JVKGZ85YjuAcURaFFOYghnQFKlRYp0nbsTazGimM5zIQ40QowwfW4kicHjy/Ex+FZFv8frd6AuPI2CK4UYc2FAuY+cK9/pFv7/o3mw0OZ4PHH5/Xp9AaExVQxW7Pfncu3eUFB0avGE5virF65Y0u1oh9v7UkDjy/Eku1JyKiZOmda0tCNHbFV+ORINv51IAMvhiWNew0RgB4sAPft24e5c+fi9ttvx+23344FCxYgMjKS+b1arcbq1atxxx134LbbbsO7776LtjbrPPdsGUB5jd3jjCvNxdt70yxe5iZYxiUJPXm8vCMZ28T00+TDG6KRW9+FJRNsG94XEIlNQvs3Jic4ntHVgU9sikPg9dIJ84/KW3ptzj0icEdSVfvwyn2sRav/mbWdE6YEHEyumfL/uwsHhre0vzk9UjA4oNHhk+HKeB5fiEPJNRiwcU4paVLiqV9HWuX1OuBvoaRJidd+SzH5zq7kTV4ARFEUs+o3Oq4XjH/IIwLQgwXgjRs3IBKJIJVKUVVVBT8/P9xyyy0oLaVd0letWoW///3viI+Ph0QiwYIFC7Bw4UKr3sPWAdSjGsKaCwXjfIqMW5K746UekZjsahhX/UbHD+cLmFXZOevEWH0mD/uTanAlT46t4gp8cCgDj2yIhqxrcld5gvsxoNHhu1HbwcaHA0lDt8nrhEUt+N+FQvAvT5xgT3AMAxod8wC9JDQJp7Maxhmta3R6xJW3YcVwn20en24rGHi9FDtjpRBcKcLZ7EaPSq0xVtUau4BQFIWVx+nrvz8gCtfym6wu+pB1qXA8vR7LDmQw+dJPb4lHfcfA1P+ZA4xFf9MFtBn8VnEFDGOuYUCjw43CZmyJqoDgSjG+PJGL14bzRD88lIkvTuROuF1MBKAHC0Bz/OlPf8Lhw4ehVCpxyy234NKlS8zvKioq4OPjg8xMyxrQA+wHkKxLhdCYKvz3bD4EV4oRUdBE8pHsiMFAjVvpM7ZB+ua0BEqV6VOtVm9g8gbf2pPmURMGYYRetRYXcmTMg8D9AVEmgv+tPWn44XzBlK2lCI6hvmPApIsFj0+3+Fu8LREvhCaZpNlMEwjhd7V43N+2p2HsKX5quAjEuFU+yy9y3APNVEgauhkj/H+O6pT09cnx90h7USjrxtO/xuPfR7LGCc4+tRYlTUpk1HRiXtB4/9gH10XhpAW9jIkA9BIBqNfrce7cOdx6660oKytDfHw8fHx80NNjWgl19913IywszOLjkgHkfigHtXhl2ArEGB8czBz3ZGmkVanGPzbFYvnhLJQ2T9zVg+D+1LT3Y846MVMoZMTYe3S6QIj/ns1n+nMPDulxLrsRb+9NY5VYT7CeFuUgfouTYtn+DLN51Y9tjEHQzTKX76XOFWuGLat2x0tNHA+s8fpr7FThP8NFU8b4/nwBjqTWQd7tuIefM1mNuNefLjT5fIyfp7CohfmdMRaExMHvajHCYqoQLCxDenWHRe9D5m8PF4DFxcW47bbb8H//93/44x//CJFIBAA4c+YMbr311nGvnz9/Pn755ZcJj6fRaNDb28uEXC73+gHkjgzpDLiQ04j7hm8kufWTW8gYW9Z5ckcWb+dQci32JdbA72oxs+J7KrMem0TlzCqxMR5YG4VfLhWZtKp6YxfpLOIsBjQ6VLX1IbO2E0lV7WjqGfS61fpdcbQP4JoLBUiRtjMpLZbkkUcUNOH9gxkm5vuLtyVCXNrqgDM3JaqklTmHT45km+QaKlXaca1D392Xju6BIZveiwhADxeAQ0NDqK6uhkQigUAgwJ133omysjKbBWBgYCB8fHzGhTcPIHfmmS10jslUnTjEpa3EC9DDWRASh6iSVry5OxWz14lt6lEdV25dERmBwBXGjkYhonLmIcb3arHZ11IUhRMZ9WaN9T85kj1p/3J7otbqmbadgddLx+3KxFe0gccXYuHmeOgNFOuiLCIAPVwAjmXJkiX46quvbN4CJiuAnsUHB2nrhMlsQAwGCu/uS8fDG6KhIx07PJYTGfVQa/WoVvRbLfyMca9/pNMmT4L3UtqsZLbB4yvamLzlyxLzFbObROUm43aagO6RO7oTkjO4POzSsHBzPIZ04++1xlXO13elcPJ+RAB6mQBcvHgxPv30U6YI5PLly8zvKisrHV4EQnAuS4e9pW4WTWwqKmnoYm6U+Y3WJVMT3IMBjQ617f24PlxNaKsA5PGFmB8cS3JFCQ4l4BrdsWjl8RwcTK5hRJ2xeMJgoFDZ2ofkqnbsSRixSNmXWAN5N92b1xUwVm1vE4/fbclv7GbOO5QjX1Yyf3uwABQIBEhOTkZ9fT2Ki4shEAjwu9/9DjExdOeHVatW4e6770ZCQgIkEgmefPJJPPnkk1a9BxlA7kvXwBCT8zLZk++ARodXh20Fpmq0bgnxFW1kJdGFaO4ZZL5fa0zaJ4sH1kahgEWrLQLBGozdbgKvlzLFH89uScDK4znYFScd56XH49M+s66WJ2k8z403x/uuhqfUTro6aAtk/vZgAbhixQrweDzceuutuOuuu7BkyRJG/AEjRtB/+tOf8P/+3//DO++8g9ZW65JeyQByXzJrO8HjC/HMloQpX9uv0WHp7lTc4ytiNbHn1HfhuW2JNv9/Avfsjpfi65MSLNtvW6vGiWLx9kRnXxrBSzCu+n10KAsPr482Ox6nCWiPy8+OZuNIap1Ldpgyehl+ejR73O/SqzvA4wsxe52YM49cMn97sAB0BGQAuS9Gt/ile9Isev31wmbw+HS1p60dWp7fnsiqhynBfpgzCWcTi7ZO/WBBIHDBqlMSWuQNj72ZfiIE3SxDeEotvjuXj3URJW5hZH8tvwk8vhDvH8wY9zuDgcKc4ar7ytY+Tt6PzN9EALKCDCD3pKqtDw+spas8L+RY5t9W2qxkJndbejQ39QyCx5+4WTvB+Rj7jXIRs/wisTNWalG7MgLBVlRDOjweFMOMu2e2JCC5qt3Zp2UT4lLaAua138wXeTw2XCFc0sRNji2Zv4kAZAUZQO6HWqvH89vpPJkPD2Va3B6pvKXXZILfHFlh1ftGFNBPt95iTOuOnMlq5HQVkMc3n9BOILBFqzcgNKYKj2wY2fJdfjhrQkN7d6BFOcjkZY/tukNRFNOysUU5yMn7kfmbCEBWkAHkfuyIrQKPL8TjwbHo7Ld8dYaiKHx3bsQl/+EN0VYlI6+NKMG8oBiXS7wmjEBRlEn/WC6Cf7mItHckcMqARoflh7OYMbYgJA6HU+s8orjs7b30Kry5/r0LQuLA4wuRXmNZp4+pIPM3EYCsIAPIvdAbKMwb3i65PkGD8MkY0hmYHpk8vtDighCKorB4eyK+PJFr9XsSHIu8W8VZNbAxXgxL8ojJmeB8ypp7mar1B9ZG4Wx2I2dFEa7AyuN028UjqXXjfmdsdxckHF8lbAtk/iYCkBVkALkX0rY+8PhCPLg2yuYJ+Vz2yDahuZuUOdJr6Ao2c/YGBNfjm9MSzreCxza0JxCsoba9H//cl86Mp0c2RHukL+n+JLqi+dmtCWjuMd0GTqxU4NMj2VgQEmfV7s1EkPmbCEBWkAHkXiRUKsDjs7PooCgKG2+WgcenLQmK5FOvAn51kn6qPZ3VYPP7EhyHrEuFWX6RnArAK3nmuzIQCFPR1DOIh4dz/Wb4irD6dB4UvWpnn5ZdGNDo8MO5fDy0PhqPbYxBc89Ivp9Ob8DLO5LB4wsRnlLL+r3I/E0EICvIAHIvvjhBCzHBFfM9Mi3FYKCYrYqHN0QjRTpx1Z3RcPq+gEiXcdwnTM1nR7M5FYALN8ez7l1K8C70BgpH0+qYnuWv/ZbCWQGEKyPrUjF/N2OFnrGTybL9GawLXsj8TQQgK8gAci+MPlKpUvZJxEqVlskHnCaY2N6lpIm2j1kbUcL6PQmOI7mqnfNt4PU3Sp19WQQ3QNGrxp6EaqYilsenWwyOrYz1ZHyv0qbQ353LN/n56NV5wZUiVvmPZP4mApAVZAC5F8YcmnUR3EzEij41XhlOyJ7pJ0Jufde418SWtQ2/Zwl+vlSIXXFSnMtuRPfAECfnQGHI1TwAACAASURBVLAPFEXZZRXQG1ZwCLZjNKg3xj2+IoRElnudn6RRAP54sXDc7y5J5Mzn89y2REQWt9jkrkDmbyIAWUEGkHtRrejDs1vp7ZQ2jnJo1Fo9BFeKwOPTVXmj7QsoamSreGw8GRKHDg4SmQn2o0+tZTwjuRSBtcQLkmAGo1k8jy/EnHVinMps8Nq0gQ036Dzrj8IzzRbsJVYqMD84lvm8XtmZgiOpdShtttwkmszfRACyggwg92Ph5njw+JZbuFhCR7/GZJIPuFaCFuUgLuTKJhUDy/ZncNbYnGAfTmY2cL4VPGedGOtvlCI0uhL7EmuQUKkg/pBeCkVRKJD14I1dqSZjxBV79TqSYvlI56WvT0rMvqZPrUVoTBXT1ckYJzLqobXA5YHM30QAsoIMIPdjSWgSeHzuzEQBehVw7CQ/TSC0yE/u1yjrOooQHMtUIp6r+OxoNtJrOty6kwPBOhIrFXhh+H40NsiDIfC/C4XM5zEZ3QNDOJRci8WjPFrnBcXifM7kbTfJ/E0EICvIAHI/nt5CrwBy2S9TrdWbtGSyJu71j0Sr0jMtHTyFpKp2fByeZZKUb694YlPcpFXlBM9Ab6Awd7go7R5fEV7ekYwXw5LwcXgWYsvanH16LkGBrIf5uxjQ6KZ8vUanh+/VYuZz5fEnN/wn8zcRgKwgA8j9MD5xrzlfYHEfYEs4m217H9mfL41PdCa4HoeSax2yGjjLP9Jjfd4INNcLm5nvu6qtz9mn45LoDRTzGe1JqLb4b0KrNzA5hP/YFDthK0YyfxMByAoygNwPoxk0j08nGHNViKHVG0z6c1obv0ZVkDwwF6ejX8OsINs7Vp2SWLTqQXA/2vtGcobnBcV4VCs3rjF2BuHxhZguoHtrN/VMXUmv0enxxCa6d3DQzVLozGypk/mbCEBWkAHknkQUNOH+ADpx+Klf4yFt5eYJXKc3YDGLqtGokhZOzoNgP2RdKnx4KBMz/bjtF2wungyJQ3Rpq90fDBR9alyWyHEyox6Shm6ohojwtCdrI0qY77haQVb/JkNvoBhLmNGxPbpySuEcUdAEHl+Ih9dH4x+bYlHSZFohTOZvIgBZQQaQ+yJt68OLYUl4/bcUPLs1Af0crbbsjJXaPOG/sjOFrAa4CXHlbZhhQZEPF/HyjmQcSa1DF4fekRRFQdalQtDNsnHXMV0gxIpjOaRzjR04nl7PfM6kN7jltCgHIbhiKgQfD47F9cLmSbeGT2c14Ifz+eDxhVi0NcGkOpjM30QAsoIMIPemWzWEJ0PobYIfzhdwstIyeovZlvj0aDapAHQToktbmW4wjohZfpHYk1DN6pwpikLg9VKsOJaDFcdycDKjHvmN3XgpLHnc+723P90iOw2CZai1embn4ceLheSztYHmnkGsjSjBg6OsX2b4irAzVjrhw3NnvwbzgmLA4wtxJmukMpjM30QAsoIMIPcnp76LWQE5mlbH+njybhXriX53vJSDKyM4Ar2BwmWJnGlPZe+YJqAT4q0VD61KNUqalAi6WYbZ68Q4kVEPAOhVa/He/vQJ3++SRG6HT807aVWqmc+V5PuyQ9rWhx8vFjLtPXl8IRaExEHWZb5d3uHUOvD4tAdnfccAADJ/A0QAsoIMIM9gdPuloJtlrG7OBgNlclOyJWb6iTjpV0xwDB39mgk7vtgrntmSgC9O5GLZ/gy8visFP10sREKlYlxle79Gh3URJbjXPxLLD2fhZEY9Igqa4He1GEv3pDHG6BPFyuM5TvpUPY9qRR8eXBuF57YmOPtUPIpLEjnjufr6rhTkmGnJqdUbmAedxdsT0darJvM3iABkBRlAngFFUfg4fKSCN6qkldXxjP1/2cTsdWJiBeImdA8MOWwFcKp4MiQO1/Kb0NZLF3aMFnifHs1mWiFaGi+GJTn74/UYiuQ9zHdE4Jbyll7MXjfy4P3duXw0dA6YvKZVqWZSfhZtTUBtU7vXz99EALKACEDPoV+jY7qEPLMlAY2d5rcSLOWzo9msJ3PBlWKOro5gb/Iau/HYxhinC0CuY26g2NkfrcdgFIDzgmLIFrAdaOgcMHkQmxMoRrXCtO+2rEuFp36lH4reDov1+vmbCEAWEAHoWSh6R54Q7wuIxJU8Ov9J1qXCxVwZRMUtFptHlzX3sp58pwuEkBKTWLfhu3P5zHfnqAphruL+gKgJja4lDeO31AjW06vWYpY/LVCIz6N9qFb0YZOonBm7T2+JR217P1qUg4zormnvx4Nro/D3Hy56/fxNBCALiAD0PJp6BjFn1FbCws3xmDaqBdjOWMsLNEYLAltjxTGSg+UODGh0zMrChVwZPj+WM+67fDHMfN9XZ0fQzTL0qIZwNV9u9vf/OZvv7I/XYzib3cg8HKyNKHH26XgsRfKecbnYr+9KYQpAtkdXEgEIIgBZQQSgZ2KsMBt98zCKwpd3JFt8nIbOAeaJn00s258xYXUbwTV47bcU5vsKjaky6WAwLygWhbIeGAwU/nuW/UMBl/HQ+mhmJcpYKTk2Xt2Z4uRP17O4WUS3gZsmIFXW9kTWpcInR0xTcRaExGFvYjVuFjVj+b4Er5+/iQBkARGAnk2BrAc3Cpsh61Ix/n6LtydadYwjE0yq1sbJYdsOgmtyYJTg+/FiIUJjqph/Xx41yVMUhd3xtpuFcx2jzYhHd6gYG5KGbmd8rB7LulGfdXoNqfi3J0qVFhdzZXhwbRTeP5jBfO5kBZAIQFYQAegdVLb2YdFw9eRPFwut+r8GA4V/H2FXEPLhoUws3Z1KjGNdmH6NDguG80d5fCEeHS4ImbNOjMEhU4Patl41q/HAVSw7kMGs/ml0+km3qF8KS8bJjHpk1nZy2pHEW6Eoivls39uf7uzT8QoUvWoIi1sY83YiAIkAZAURgJ6PolfNbAc/GRIHebf1W7EtykGbe8fe4yvCyuN0PtnLO5LR3qexw1USuGBPQvW472/NhYJxrzuZ2eBU4TfTT4R9iTUmBU2lzUqL//80gRB+V4vx/bl8pFd3QKPz3vaFbLr2fH1SAh6fdh0gOA6KonAsrQ4zfrrs9fO3xwrAkJAQPP744/jDH/6Au+66C2+99RYqKytNXrNo0SL4+PiYxNdff23xexAB6PnEldOefktCE9HZb7v4GptTaE3MD47FimM5eHVnCr4+KeHw6ghcotbqsfpMnol4rzJTxc1mLLCNZ7cmoKRJafb8l+5Otfkh5bGNMbg/IAoPb4jG0j1p2BkrhWpopMpVb6DGrYS6O+Eptfj8WA5WnZKg24ZV0WUH6O3ILVEVdjg7wlS0d3V7/fztsQLw5ZdfxrFjx1BaWorCwkK89tpruPvuuzEwMGIOuWjRInz55ZdobW1lwprBQASg5yProlu7zfAVoba9f+r/MAF1HQM2rwIaY/XpPDy0Phqi4hYOr5DAJRRF4T/DhR4TifXj6fVOEX9PTtIqCwCW7c/g9P3mB8fiP2fz8XF4FuYEijFnnRhrI0oQGl2JqJJWRJW0oEDWA7VWD7VWj58uFtpthZuiKJQ0KSEsakG1og8lTUpk1HRO+X9alINIq+6ApKEb5S29OJFRj69PSvDqzhTMDRRjxrD43SaunPRY5jD2X46vaLP1sggsIPO3BwvAsbS3067fyckjVZyLFi3C999/b/MxyQDyDj48lAken87rYiO+gm6WcTKxzguKNVldIbgWufVd4PGF+Dg8y+zvFX1qE2shR8Qsv0g09wxOet7Gce7omOErwuPBseDxhXj/YAbO5zSiqq0PBgs8N42rihqdHqLiFpzOakBChQIVrb3Iqu1EclU7mnsGmcrb0TFdIMTS3anYJCpHQoUCWr0BAddKsHRPGt7YlYqHN0RbdQ3p1dYVcxjHwMVcmVX/j8ANZP72IgFYXV0NHx8flJSMeC8tWrQId955J/785z9j9uzZEAgEUKksz/EiA8g7UPSq8fKOZOZmf72w2abjKAe1eMSKSWWyOJpWx/FVErhCrdVjll8kFm6On/A1r+5M4WQcWBL3B0QhRdo+5XkvP5zlsHOyJOYEivHFiVycymxAXHkb3t6bBsGVIhxOrcPh1Dr4XyvGLP9IBF4vZYq02MRclj28BVeKLB4jla19zP87ldlg8f8jcAeZv71EABoMBrz++ut46qmnTH5+8OBBiMViFBcX4/Tp0/jrX/+Kd955Z8LjaDQa9Pb2MiGXy71+AHkLQzoDfrlUBB5fiH8dyLD5OKP94djEJ0eyObw6AtcYCz0mWnUb/UBhz1i8LREFsh6LzvnbUfmLJKyPrWLLcvlqhztR8Ph0z2/loNbicUXgDiIAvUQArlq1CjweD3L55Kab8fHx8PHxQU1NjdnfBwYGjisa8fYB5E1UtNLt3easE9vcy1M1pMO8IPY9Y3fFWd6RhOAcvj2Th6Sq8Stv0aWtdhcj0wVChESWQ621vPBC7IDz8uT434VC1HUMTNouUm+g8MyWkdXKxEqFTWOLwB4iAL1AAH777bf429/+hrq6qbfMBgYG4OPjA7HYfAN0sgLo3ai1eubGzeapnQtz6NBo65POCY5FrdWPswmhKMru27/TBELElFlfWKDVG5hcPBK2x30BkXh7bxqqFeOLxgplPczr8hqJubYzIQLQgwUgRVH49ttv8Ze//AVSqWWrJWlpafDx8UFRkWW5HGQAeReKvhEDXzZFGGqtHk9simM90SSbWV0iuC5KlRY/X7K/BczK47k2n+P5nEanCyhPCXNbwmez6c93+WHzBUIEx0Hmbw8WgN988w3++Mc/IikpycTmZXCQzsmpqanBxo0bIZFIUF9fj+vXr2PGjBl49tlnLX4PMoC8C2N150w/kc1bwEa4mGif3Zpg1RYfwflM1G+Xi3hzdyqW7c9An9r21WmDgWL86Uiwi7F9wzNqOpnfkRV850Pmbw8WgOZy9Xx8fHDs2DEAgEwmw7PPPos77rgDv//97zFz5kz8/PPPxAeQMCGKUS28LLGomAyKorBNXMl6ktkcSUxk3QWKovDmsNnyLP/ICb/T2evEWLw9cdzPl+5Jm7Kt4NnsRtbn2aMawisOrFL25GhR0gsOBgNl9ucE50Hmbw8WgI6ADCDvora9Hzy+EA+ujeLkeBRFQXClmNUEM00gxOu7UhASWU6qCd2E9j4NCmQ9WH3afNXtdAHt2zf6Z49ujMGa8wWMcHx3XzpCY6pwfwBdTfrm7lREFDRNWoBgDbIuFWtbFBJCpFd3QDWkw+fHcpifnc4iti+uAJm/iQBkBRlA3kVmLb2Fs3hbImfHVGv1CLpZhukcGQNvE1d6XMstT+VG4Xhz4me3JuBAUg3yG7tRIOtBbn0XEioUCE+pheBKMb4+KcGHhzIhuFKEwOulCBGVI7O2k3VKgjmy67rw6Eb2FeveHN+dy8e/xmyprziWjTQrTaMJ3EPmbyIAWUEGkHeRW9+F5YezOPfg61EN4WByzbhVH1tjSWgSqhXje9ASXAud3oDwlFoEXi/F+hulEJe2craCxxWyLhXe25/udCHlzvFSWDIeGPb94/GFmOUvwqJtCc7+ar0eMn8TAcgKMoC8i9JmJXh8IR7ZEM3pce0xwc5ZJ0ZFKxmXBPZQFIXyll7sipNi6e5Uh7ex84R4KSwJHxzMwBfHc/D89kR8eSIXulEWQZ39GnT226cPMsE8ZP4mApAVZAB5F1VtdPum2evM+0TaysLN8XaZdJ7YFIcOMqkQOKZfo0NWbSfOZjdik6gcT2+xz/j15HhofTSq2vqgHNQi6GYZoktbnf21eh1k/iYCkBVkAHkXV/Lk4PGFeHtvGqfHtWdbsI03yzg9VwJhLDq9AYmVCuTUdyFV2oEXQpOcLrDcIT45ko03d6eisdPy/vME7iDzNxGArCADyLuILG4Bjy/EC6FJnB532X77+a7d4ytCeQsZnwTHodHpyaqgFRFZ3OLsr8wrIfM3EYCsIAPIu5A0dIPHpy06uKy6XH+j1K4TzNLdqax9CwkEayhpUuKRDdFOF1fuEIuIobtTIPM3EYCsIAPIe1AN6RgT35XHczg9ttFexp5hS29YAoENxpQJElPHxptldrHy8WS0egNq2vtR2qyErEuFQlkPtkdXYsWxHLwUloxntiRg8fZELDuQgXPZjeM+XzJ/EwHICjKAvAOKovDliVzw+HQFcH3HAOfHX2LnvKk1Fwo4PWcCYSp0egNWnZI4XVy5S/znbD56VENWfcaJlQr0smj95450DwzB72oxY4JuaVyWyJljFMh6sOZUhtfP30QAsoAIQO/gsoReyZjlHwlJQ7dd3uNirsyuk8s9viIkVbXb5dwJhInQ6Q34ZIr2dSRG4sG1Ufg1qmJCSxiKopDX2I1t4kqExlQh6GaZRxu/Dw7psStOis+OZmNPQjXe3J2Ke3xFzOf1wNooPB4ci3v9IzFNIMTH4Vk4mVGP5Kp2SBq6kVnbifnBseDx6YfgAlkPqtr6MMs/En//4aLXz99EALKACEDPJL2mAyuP52D54SzsiK3CE5viwOMLsd3ODdzDYqrsOrksP5yFAY3OrtdAIIylRTmIOaStnFVxf0AU/nehEKcyG1DSpIRGp8eNwmYTy6jXfkvB0CgvQU9CpzfgRmEzFm8b3xObxxfixbAkpNd0MNu6FEVNmOd8LrsRPL6QEYLGmC246vXzNxGALCAC0POILm0125btyRDHeOpl1HTi50uFdpswR2+DEAiOIqOmE/f6c9PphgTdL7pYrnT218oper0BWbWdWBdRgsdHibV5QbFYeTwH/9yXjlWnJMio6bSqqE1voLAnoRpBN8tMPsMXfo3y+vmbCEAWEAHoWaTXdGDusPBafToPIZHleDEsCd+fy4eiV+3Qc+lVayd8+mUT/9yX7tDrIBCMlDX3Yv2NUrz2W4rTBZS7xr3+kfj6pATnshud/XVyhqxLhTXnC/DIhmj8Y9OI8Ht0Ywx2xFZBOchdjmNyVTvTyeb7E+leP38TAcgCIgA9A4OBwoGkGswYzi15e28aNDrn59VI2/owex23K4EPro0i1YYEp/NSmP3Mz909lu5ONRHJT2+Jx4pjObiYK0Obgx9EuUarNyC9ugOHkmtxOLUOxXIlk2LD4wvx+dFs/O9CIRIqFHbb3pY0dOFCrgytHV1eP38TAcgCIgDdH4OBwprzBcwN6IfzBS7lyVXb3o83dqVyNrk8ty3R2ZdE8HIoijKZ9EmYRnpNB5SDWvheLUZsWZtdPDxblWrIulToVWtR2qxEP8e5wQYDBUWvGnmN3biQK0PQzTIsP5yFh9ab94acFxSDzNpOaPWOy2kk8zcRgKwgA8j92SauBI9PV8mezmpwydWxIZ0Bv0ZVMFsXbCKItIYjOJnESoXTRZYrx/8uFNr9PtQ1MIRZwzmZv8VJkVPfZdV7anR6JFQqcDVfjou5MuyIrcJPFwvx4aFMLNqagFl+E+d7ProxBqtOSZiCln/uS0dVW58dr9Y8ZP4mApAVZAC5N31qLR5YS3tJXciROft0puRAUg3rySU8pdbZl0HwYtRaPeYFxbIex54eexKq7fo9iEtbwePTK29j6RoYwpHUOry5OxUPro3CL5eK0NA5gBRpO7ZEVeCni4V46ldavL0YNrF/6XQBXTz3UXgmAq+X4mx2I/Ibu6EbXuWjKMqpqTZk/iYCkBVkALk36yJKwOMLsXh7okuu/I2Fi44hqdIOZ18GwYs5O2zJQWLqCLpZxogle30HM/1ECIupgn54m7mxU2Xx+c0LisW6iBIsP5wF/uUi7IqT4rJEjszaTsi6VA7dzrUFMn8TAcgKMoDcl2v5TcyNzF0MkkublawnFbICSHAW5S29zLYjCcvis6PZnK+ShafUmthMlTWPzF+Shi4mP/OZLQkQFrUwBSmPbozBz5cKGaHn7p6iZP4mApAVZAC5J0XyHryzNw08vhBboiqcfToWo9bqWU+g84JiXKrIheA9kI4gtsXnx3I4F1ujOw+tOJaDb8/k4Vx2I6ra+qDoVSOjptPk9Vq9wS7FKM6EzN9EALKCDCD3Q63VY0FIHKYJhFh7rYTZ+nAX/nehkPWEIixqcfZlELyM7Loupwspd45Xd6agVcmdBYyouIU59oNroxBT1sbZsd0FMn8TAcgKMoDcix7VEP49vApxX0Ak+tywibqiT83aG/B6YbOzL4PgJVAUhaiSFjy6McbpIsrd46WwZPRycM/SGygTn8Glu1Pd8l7IFjJ/EwHICjKA3AeNTo8XQumKtWkCIU5mNjj7lGzmRmEzK0uYakW/sy+B4AWohnT44kSu04WTJ8VH4ZmsDZLPjSnEed5NiuC4hszfRACyggwg9+Fg8oiFiqShy9mnwwq9gUJoTJXNk8h/zuY7+xIIHo5aq8cHBzOdLpg8MVafzrNJBGp0elzNl+PBYeur0RGeUmu3zhuuCpm/iQBkBRlA7kGPagj3BdDFE75Xi519Oqyp6xjA/GB2XmqFsh5nXwbBgwkWljldKHlyvLk7FXUdAxZ9F809g9gqrsBjU2zDF3jZPYHM30QAsoIMINcnpqyNufEt3BxvF18tR6NUaVlPIM9sSYCiz737ihJck9z6Lk661pCYPO4PiMK57MYJt28VvWqEiMox009k0fFOZtQ7dqA4mZiCOq+fv4kAZAERgK4NRVF4egvtWP/g2iiPMUGmKIrpYMImXtmZ4vZeXgTXQm+gmFxbEo6JF8OScDKzAdK2PtS09+N8TiM+Cs+0WoT/eLHQ2cPHYSh61Zjx42Wvn7+JAGQBEYCuTVlzL3Nza+/TOPt0OOWzo9x4qoVGVzr7UghuDEVRUPSOrCRfzZeDxxdi+eEsi1eeSLhGrD6T58SR5FhKmpR44dcor5+/PVYAhoSE4PHHH8cf/vAH3HXXXXjrrbdQWWk62anVaqxevRp33HEHbrvtNrz77rtoa7PcD4kIQNfGKADv8RUhu869Cz/G4nu1mJOb/n0BkahWOL4RO8EzCLxeio/CMwHQRQbPbk3A/QFRKGlSYvG2RKeLGhKWx4GkGiePJsdC5m8PFoAvv/wyjh07htLSUhQWFuK1117D3XffjYGBkcTZVatW4e9//zvi4+MhkUiwYMECLFy40OL3IAPItaEoyqRZub0brDuSpbtTObvxPxkS53UJ4AT2pEjbwePTubXrb5Tig4OZmC6gWyt+Sexf3CqmCYRoUQ46e0g5FDJ/e7AAHEt7ezt8fHyQnJwMAFAqlbjllltw6dIl5jUVFRXw8fFBZmamRcckA8j1qW3vZ+woHg+OdfbpcMabHApA4wSw7EAGAq6VmGzpEQjm6OzX4OEN0ePG0dxAMRIqFZhOikDcKt4/mOHsIeVwyPztRQKwuroaPj4+KCkpAQDEx8fDx8cHPT2mKx933303wsLCzB5Do9Ggt7eXCblc7vUDyB3oGhhibnSeUAUMAMsOZNhtMnhofTQiCpq80hyWYBlJVe1OFy0kuIuz2Y0TfteShm4EXCtBsLAMb+9Nw7dnJBAWtUA15N4FZEQAeokANBgMeP311/HUU08xPztz5gxuvfXWca+dP38+fvnlF7PHCQwMhI+Pz7jw5gHk6lS19THt3x5cG+UxoibweqndJ4Vvz+RB6yGCmcAtW8UVThctJLiLmvaJuwMdT683ee2/hh8+7wuIhP+1YrfdOiYC0EsE4KpVq8Dj8SCXy5mf2SIAyQqge0FRFLZEVWC6QIjpAiGOp9c7+5Q441p+k0MmhvU3SmEweIZoJnCDwUDhmS0JThctJLiJe3xFk+6M5DV2Y3t0JdZGlGB3vBT7EqtNvv8H10bhVGaD23USIQLQCwTgt99+i7/97W+oq6sz+bktW8BjIQPItaEoCvOC6I4ZR9Pqpv4PbkR6dYfDJoiXwpKJXyCBobyl12Fjj4RjoqrNOicAiqKQUdOJd/amMceYFxSL/57Nx8VcGdRavZ1GH3eQ+duDBSBFUfj222/xl7/8BVKpdNzvjUUgly9fZn5WWVlJikA8DGNLqu/PcdP/VtLQja9PSvDJkWzsiP3/7d1/TJT3AcfxhyI/gloKgkXFO3/EqClq02EC2rRpaNqkNjWaqrGa2ujSYE2FzUSwtblNBzbrj8SsZWvnQud0zB9trYoVLYpp1anMCmiLYKWKikWFA52FyfHZHwtPe+NQOC139vt+Jd8/fJ6Hu4ePX7gPd/d872SPf3HeKZ9W1PXqA8SGw2cD8n0i+Bw4dTnghYVxZ8c/Dnf9HsCbafO06y+fn+700ZQz/nRAO8ovBPWrBzx+/4wL4MKFCxUdHa2SkhLV1dXZ4/r1H96vkJ6eLofDoT179qi0tFSpqalKTU3t9n0wgYLfodNX5MzarodW7Lrt29p5vK7T4rYz/rg/IOvoVV1s7tUHiDHLPw1Y2UVwudZyQyOXscjzz2lMz9t/W3Oi9YZHn1dd0lu7TmrRun/ZtztpVbF+veGY0v9WqsUFR4NqzVEev3/GBdDXxRqWZSk/P98+pmMh6JiYGEVFRWnatGmqq6vr9n0wgYLf1ZYb9i+jpu//4/ftHKm5otHLd8iZtV0LPjiitQdqNO3dLzQiu1CTVhXLfd3/2/aHx9Ou8b/pvAzHTzkWrivt1e8RwesXK3cFvLQw7tyY+Lvd2lv5ndruwDN2Ddda9VZRpca5dvq8rwm/LdKs9w7o7JV/34GZ6D8ev3/GBbA3uN1uWZal2tpar4tDGME1xmZ/qKGZG3X01Pkef+3Zi5c1+91iDc3cqKGZGzX73WI1NLrV1NSkC/VX9NDyLRqauVEz/1Dc69/X3Lw99nn11iitOhfw/09G4MejOYW9PvcYP90YsWSTGhrdutzQeMfmyHeXG7T5YJV+v/WofvnnfZq8Ypscv/K+3z/tLlfNhUtaW/KVlv79oBb/9QvNydujp98s0of/rOrW/TQ0unXgq7OqPHOxR+fXcRGn2+0OdJUIGArgbeiYQAwGg8FgMO6+8ePVQUxDyswK9AAACZpJREFUAbwNHo9HtbW1crvdP9lf2h0lk2cZyYzcgnOQG5mRW3APX7m53W7V1tbK47m7lq+5kyiAQa6pifcp9BSZ+Yfc/ENuPUdm/iE3/5CbbxTAIMfE7Tky8w+5+Yfceo7M/ENu/iE33yiAQY6J23Nk5h9y8w+59RyZ+Yfc/ENuvlEAg1xLS4tcLpdaWloCfSp3DTLzD7n5h9x6jsz8Q27+ITffKIAAAACGoQACAAAYhgIIAABgGAogAACAYSiAQSA3N1fJycnq16+f4uPjNXXqVFVWVnod0/G5xbGxserbt6+mT5+uixcvBuiMg0NeXp7GjRun/v37q3///kpJSdGOHTvs/WR2a6tWrZJlWcrIyLC3kVtnLper0ycIjB492t5PZl07d+6c5syZo9jYWEVGRiopKUlHjhyx97e3t+u1115TQkKCIiMjlZaWpqqqqgCeceA5nU6fn1rx0ksvSWK++dLW1qbly5dr2LBhioyM1IgRI7RixQq1t//w+cbMNW8UwCDw5JNPKj8/X8ePH9exY8f01FNPyeFw6Nq1a/Yx6enpGjp0qIqLi1VaWqqUlBRNmjQpgGcdeFu3blVhYaGqqqp08uRJvfLKKwoLC9Px48clkdmtHD58WMOGDdP48eO9CiC5deZyufTAAw+orq7OHpcuXbL3k5lvDQ0NcjqdeuGFF3To0CGdPn1aRUVFOnXqlH3M66+/rujoaG3ZskVlZWV65plnNHz4cH3//fcBPPPAqq+v95pru3fvlmVZ2rt3ryTmmy85OTkaMGCAtm/frpqaGm3atEn9+vXT6tWr7WOYa94ogEGovr5elmVp3759kiS3262wsDBt2rTJPubrr7+WZVk6ePBgoE4zKMXExGjNmjVkdgtXr17VqFGjtHv3bj366KN2ASQ331wulyZMmOBzH5l1LSsrSw8//HCX+9vb25WQkKA33njD3uZ2uxUREaGCgoLeOMW7QkZGhkaOHKn29nbmWxemTJmi+fPne22bPn265syZI4m55gsFMAhVV1fLsixVVFRIkoqLi2VZlhobG72OczgcevvttwNxikGnra1NBQUFCg8P14kTJ8jsFp5//nllZmZKklcBJDffXC6XoqKiNGjQIA0fPlzPPfeczpw5I4nMbmbs2LHKzMzUs88+q/j4eD344IN6//337f3ffPONLMvSl19+6fV1jzzyiBYvXtzbpxuUWltbNWDAAOXk5EhivnUlJydHTqdTJ0+elCQdO3ZMAwcO1Lp16yQx13yhAAYZj8ejKVOmaPLkyfa29evXKzw8vNOxEydO1NKlS3vz9IJOeXm5+vbtq9DQUEVHR6uwsFASmd1MQUGBkpKS7Jc9flwAyc23HTt2aOPGjSorK9POnTuVmpoqh8Oh5uZmMruJiIgIRUREaNmyZTp69Kjee+89RUZG6oMPPpAk7d+/X5Zl6cKFC15fN2PGDM2cOTMQpxx0NmzYoNDQUJ0/f14SP6Nd8Xg8ysrKUkhIiPr06aOQkBDl5uba+5lrnVEAg0x6erqcTqdqa2vtbfzAd621tVXV1dUqLS1Vdna24uLidOLECTLrwtmzZzVw4ECVlZXZ2yiAPdfY2Kh7771Xa9asIbObCAsLU2pqqte2l19+WSkpKZJ4UO6OJ554Qk8//bT9b+abbwUFBUpMTFRBQYHKy8u1du1axcbG8sfGTVAAg8iiRYuUmJio06dPe23nKf/uS0tL04svvkhmXfj4449lWZZCQ0PtYVmWQkJCFBoaqs8++4zcuik5OVnZ2dnMtZtwOBxasGCB17a8vDwNHjxYEi/L3cq3336re+65R1u2bLG3Md98S0xM1DvvvOO1beXKlfbV+sy1ziiAQaC9vV2LFi3S4MGDfV6S3vGm382bN9vbKisrjX/Try+PPfaY5s2bR2ZdaG5uVkVFhddITk7W3LlzVVFRQW7ddPXqVcXExGj16tVkdhOzZ8/udBFIZmam/axgxxvz33zzTXt/U1OT0W/M/zGXy6WEhATduHHD3sZ88y02NlZ5eXle23JzczVq1ChJzDVfKIBBYOHChYqOjlZJSYnXpf/Xr1+3j0lPT5fD4dCePXtUWlqq1NTUTi+tmCY7O1v79u1TTU2NysvLlZ2drZCQEO3atUsSmXXXj18ClsjNlyVLlqikpEQ1NTXav3+/Hn/8ccXFxam+vl4SmXXl8OHD6tOnj3JyclRdXa3169crKirKfmO+9L+lOe677z598sknKi8v19SpU41emqODx+ORw+FQVlZWp33Mt87mzZunIUOG2MvAfPTRR4qLi/N6WZy55o0CGAR8LfhpWZby8/PtYzoW/oyJiVFUVJSmTZumurq6wJ10EJg/f76cTqfCw8MVHx+vtLQ0u/xJZNZd/18Aya2zWbNmadCgQQoPD9eQIUM0a9Ysr7XsyKxr27ZtU1JSkiIiIjRmzBivq4ClHxbnvf/++xUREaG0tDT7Sk6TFRUVybIsn1kw3zprbm5WRkaGHA6HvRD0q6++qtbWVvsY5po3CiAAAIBhKIAAAACGoQACAAAYhgIIAABgGAogAACAYSiAAAAAhqEAAgAAGIYCCAAAYBgKIAAAgGEogAAAAIahAAIAABiGAggAAGAYCiAAAIBhKIAAAACGoQACAAAYhgIIAABgGAogAACAYSiAAAAAhqEAAgAAGIYCCAAAYBgKIAAAgGEogAAAAIahAAIAABiGAggAAGAYCiAAAIBhKIAAAACGoQACAAAYhgIIAABgGAogAACAYSiAAAAAhqEAAgAAGIYCCAAAYBgKIAAAgGEogAAAAIahAAIAABiGAggAAGAYCiAAAIBhKIAAAACGoQACAAAYhgIIAABgGAogAACAYSiAAAAAhqEAAgAAGIYCCAAAYBgKIAAAgGEogAAAAIahAAIAABiGAggAAGAYCiAAAIBhKIAAAACGoQACAAAYhgIIAABgGAogAACAYSiAAAAAhqEAAgAAGIYCCAAAYBgKIAAAgGEogAAAAIahAAIAABiGAggAAGAYCiAAAIBhKIAAAACGoQACAAAYhgIIAABgGAogAACAYSiAAAAAhqEAAgAAGIYCCAAAYBgKIAAAgGEogAAAAIahAAIAABiGAggAAGAYCiAAAIBhKIAAAACGoQACAAAYhgIIAABgGAogAACAYSiAAAAAhqEAAgAAGIYCCAAAYBgKIAAAgGEogAAAAIahAAIAABiGAggAAGAYCiAAAIBhKIAAAACGoQACAAAYhgIIAABgGAogAACAYSiAAAAAhqEAAgAAGOa/MZj3UGPeOlYAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "df = geopandas.read_file(os.path.expanduser(\"~/Downloads/datasets/natural_earth/ne_50m_ocean/ne_50m_ocean.shp\"))\n",
    "df_river = geopandas.read_file(os.path.expanduser(\"~/Downloads/datasets/natural_earth/ne_50m_rivers_lake_centerlines_scale_rank/ne_50m_rivers_lake_centerlines_scale_rank.shp\"))\n",
    "df_lake = geopandas.read_file(os.path.expanduser(\"~/Downloads/datasets/natural_earth/ne_50m_lakes/ne_50m_lakes.shp\"))\n",
    "df_lake_historic = geopandas.read_file(os.path.expanduser(\"~/Downloads/datasets/natural_earth/ne_50m_lakes_historic/ne_50m_lakes_historic.shp\"))\n",
    "base_map = df.plot()\n",
    "# df_river = df_river[:100]\n",
    "river_plot = df_river.plot(ax=base_map)\n",
    "lakes_h_plot = df_lake_historic.plot(ax=base_map)\n",
    "lakes_plot = df_lake.plot(ax=base_map)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>scalerank</th>\n",
       "      <th>featurecla</th>\n",
       "      <th>name</th>\n",
       "      <th>name_alt</th>\n",
       "      <th>min_zoom</th>\n",
       "      <th>wikidataid</th>\n",
       "      <th>label</th>\n",
       "      <th>name_ar</th>\n",
       "      <th>name_bn</th>\n",
       "      <th>name_de</th>\n",
       "      <th>...</th>\n",
       "      <th>name_pl</th>\n",
       "      <th>name_pt</th>\n",
       "      <th>name_ru</th>\n",
       "      <th>name_sv</th>\n",
       "      <th>name_tr</th>\n",
       "      <th>name_vi</th>\n",
       "      <th>name_zh</th>\n",
       "      <th>wdid_score</th>\n",
       "      <th>ne_id</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <td>0</td>\n",
       "      <td>3</td>\n",
       "      <td>Historic Lake</td>\n",
       "      <td>Aral Sea</td>\n",
       "      <td>None</td>\n",
       "      <td>5.0</td>\n",
       "      <td>Q35883</td>\n",
       "      <td>Historic Aral Sea</td>\n",
       "      <td>بحر آرال</td>\n",
       "      <td>None</td>\n",
       "      <td>Aralsee</td>\n",
       "      <td>...</td>\n",
       "      <td>Jezioro Aralskie</td>\n",
       "      <td>Mar de Aral</td>\n",
       "      <td>Аральское море</td>\n",
       "      <td>Aralsjön</td>\n",
       "      <td>Aral Gölü</td>\n",
       "      <td>Biển Aral</td>\n",
       "      <td>鹹海</td>\n",
       "      <td>5</td>\n",
       "      <td>1159110427</td>\n",
       "      <td>POLYGON ((60.21527 46.72382, 60.31408 46.73384...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>1</td>\n",
       "      <td>5</td>\n",
       "      <td>Historic Lake</td>\n",
       "      <td>Lake Chad</td>\n",
       "      <td>None</td>\n",
       "      <td>6.0</td>\n",
       "      <td>Q125309</td>\n",
       "      <td>Historic Lake Chad</td>\n",
       "      <td>بحيرة تشاد</td>\n",
       "      <td>চাদ হ্রদ</td>\n",
       "      <td>Tschadsee</td>\n",
       "      <td>...</td>\n",
       "      <td>Czad</td>\n",
       "      <td>Lago Chade</td>\n",
       "      <td>Чад</td>\n",
       "      <td>Tchadsjön</td>\n",
       "      <td>Çad Gölü</td>\n",
       "      <td>Hồ Tchad</td>\n",
       "      <td>乍得湖</td>\n",
       "      <td>5</td>\n",
       "      <td>1159110451</td>\n",
       "      <td>POLYGON ((15.29243 13.23020, 15.28137 13.12633...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2 rows × 31 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "   scalerank     featurecla       name name_alt  min_zoom wikidataid  \\\n",
       "0          3  Historic Lake   Aral Sea     None       5.0     Q35883   \n",
       "1          5  Historic Lake  Lake Chad     None       6.0    Q125309   \n",
       "\n",
       "                label     name_ar   name_bn    name_de  ...           name_pl  \\\n",
       "0   Historic Aral Sea    بحر آرال      None    Aralsee  ...  Jezioro Aralskie   \n",
       "1  Historic Lake Chad  بحيرة تشاد  চাদ হ্রদ  Tschadsee  ...              Czad   \n",
       "\n",
       "       name_pt         name_ru    name_sv    name_tr    name_vi name_zh  \\\n",
       "0  Mar de Aral  Аральское море   Aralsjön  Aral Gölü  Biển Aral      鹹海   \n",
       "1   Lago Chade             Чад  Tchadsjön   Çad Gölü   Hồ Tchad     乍得湖   \n",
       "\n",
       "  wdid_score       ne_id                                           geometry  \n",
       "0          5  1159110427  POLYGON ((60.21527 46.72382, 60.31408 46.73384...  \n",
       "1          5  1159110451  POLYGON ((15.29243 13.23020, 15.28137 13.12633...  \n",
       "\n",
       "[2 rows x 31 columns]"
      ]
     },
     "execution_count": 39,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_lake"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "ename": "AttributeError",
     "evalue": "'GeoDataFrame' object has no attribute 'properties'",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mAttributeError\u001b[0m                            Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-36-393fdbaba6da>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mdf_river\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mproperties\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;32m~/.pyenv/versions/3.7.3/envs/georasters/lib/python3.7/site-packages/pandas/core/generic.py\u001b[0m in \u001b[0;36m__getattr__\u001b[0;34m(self, name)\u001b[0m\n\u001b[1;32m   5177\u001b[0m             \u001b[0;32mif\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_info_axis\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_can_hold_identifiers_and_holds_name\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   5178\u001b[0m                 \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 5179\u001b[0;31m             \u001b[0;32mreturn\u001b[0m \u001b[0mobject\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__getattribute__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mname\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   5180\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   5181\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m__setattr__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvalue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mAttributeError\u001b[0m: 'GeoDataFrame' object has no attribute 'properties'"
     ]
    }
   ],
   "source": [
    "df_river.properties"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}