create_empty_calib_diagram.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib notebook\n",
    "import numpy as np\n",
    "import matplotlib\n",
    "from matplotlib import pyplot as plt\n",
    "import seaborn as sns\n",
    "sns.set()\n",
    "\n",
    "matplotlib.rcParams['text.usetex'] = True\n",
    "matplotlib.rcParams['text.latex.preamble'] = [\n",
    "    r'\\usepackage{bm}']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "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, fig.canvas.height);\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,iVBORw0KGgoAAAANSUhEUgAAAcIAAAHCCAYAAAB8GMlFAAAAAXNSR0IArs4c6QAAAERlWElmTU0AKgAAAAgAAYdpAAQAAAABAAAAGgAAAAAAA6ABAAMAAAABAAEAAKACAAQAAAABAAABwqADAAQAAAABAAABwgAAAADHyZ8pAABAAElEQVR4Ae2de9AVx3mnm4sulsRF2tjJlkFblXilEqBcHGnXoOy6HCF9QCzHwpFQ5Y81SEBFcgyRBGV7LZBAsMkaSTY4ThaBLs7aayHZVJJKEEhkU7FLYK2ijRMuMinHdgDHqU0icZFksLjs+Q15P80538w5M+fMpWfm6eLjnJnpefvt550zv+me7plRZ1vJkSAAAQhAAAINJTC6ofWm2hCAAAQgAIGAAELIgQABCEAAAo0mgBA2OvxUHgIQgAAEEEKOAQhAAAIQaDQBhLDR4afyEIAABCCAEHIMQAACEIBAowkghI0OP5WHAAQgAAGEkGMAAhCAAAQaTQAhbHT4qTwEIAABCCCEHAMQgAAEINBoAghho8NP5SEAAQhAACHkGIAABCAAgUYTQAgbHX4qDwEIQAACCCHHAAQgAAEINJoAQtjo8FN5CEAAAhBACDkGIAABCECg0QQQwkaHn8pDAAIQgABCyDEAAQhAAAKNJoAQNjr8VB4CEIAABBBCjgEIQAACEGg0AYSw0eGn8hCAAAQggBByDEAAAhCAQKMJIISNDj+VhwAEIAABhJBjAAIQgAAEGk0AIWx0+Kk8BCAAAQgghBwDEIAABCDQaAIIYaPDT+UhAAEIQAAh5BiAAAQgAIFGE0AIGx1+Kg8BCEAAAgghxwAEIAABCDSaAELY6PBTeQhAAAIQQAg5BiAAAQhAoNEEEMJGh5/KQwACEIAAQsgxAAEIQAACjSaAEDY6/FQeAhCAAAQQQo4BCEAAAhBoNIGxja49lc+UwLFjx9yzzz7rtm/f7g4dOhT8jR8/3k2YMMHpc/r06e66665zM2bMyLRcGduyZYtbuXJlYPfAgQMj7Gub8ig98cQTbT702neEsYqs6FbnilShzc0rr7yybTlqQcfZ5MmT3bx589zs2bOD4y4qX3hd3TiF68b3ZAQQwmScyNWDwKZNm9xDDz00IpfEUX9K+/btc5s3b3ZDQ0Nuw4YNI/KyAgJJCUjsOtPRo0eDY03HmcRNf6tXrw5EsTMvyxAIE0AIwzT4npqARG7+/PmByGlnuxpXq0/fdXJS61AnJ7W89D3qJJa64JQ7qEwrVy3UJqS61ln12rlzZ2QIrVdCF2Y61iSGL7zwQtcLr7pyigTEykgCo862UuQWVkIgAYGZM2cGJxxlVXeUrsC7JYnhtGnT3NSpU7tlS71tkO7NQfZN7Sg79E3Auka7CWHYeLjLk16IMBm+dxJgsEwnEZYTE1iyZMmwCC5cuLCnCMqwxDJrEUzscA0y6uQuQRB7UncC4W7RHTt2DN8j7r4XW5tIACFsYtQzqPOuXbucTi5KErbly5dnYBUTEMiWgMRQLUglXUTY/epsS8Fa1QkghFWPYEn+hwfGPPjggyV5QbEQ6E1g2bJlw5nUDU6CQCcBhLCTCMs9CdjgF2VUa5Cuzp7IyFAigVmzZg1Po0AISwyEx0UjhB4Hx1fX1C1qSXO1skrqttJovwULFjgNwtG9MH2qS0viO0jSCVD2bMBFN1vmx9y5c921114b7COfNPK1M5ld5bOkfLqHJ9/XrVtnq4NPs52mjlaGfLcTubqlrT72GY6LCgvv1+ZExILyhutrvnfrSlSsVLbVXXkVK4ud1ofvI0cUW9gqzWFV0nHUWacknPqJW2flVLYdF+Im3ppzK9vyIbxNbC2Zf8ZZ67M+xqysqJjqGO4WU/mvPKqPHQ9xvxcrx7dPpk/4FpEK+KPh6Jayag3qhLB06VIzG3xqcrROHjoR6G/9+vVOV/d5Jp0Iwt2+VpZERj/0Xj7oZGb3TrXv8ePHzURw0uunjvZAAhkKn8TFJ5z6mRaiE6p86rzQ0LLmfOqvV53lg/jITtg/fRcL/Wm6g92rC/tc1Hc9yMHiIl/THEdZHJs6fnWREE7GPrzOvu/Zs8e+jvjM4xgbUUhrhfzT1KhuMdW2qDz2e9m6dWsleoxoEUYdAazrSiB8cs/q5Kb5hjqxa/Spfjx6OsxLL70U/NmTaKJO2F0d7WOjRFAjW3Xilg/6lE+Wuvmgq2A72dpUkkWLFtmuwZzKfuqok7ax0DQAJX3aOvtMe1EisZO461NxlODJluqtGIS5Swzikk6GqrumxRg37W++ar9wCyfOTp7rwxcJnaLfq9xBj00JiomgjgtjLN52MaP14m5/cQ+cyOsY62RgAqeY2u+xM6b6Lej40bFjeRR/1cXSxo0b7avfn5pHSIJAGgI333zz2SuuuCL4S7Nft7ytk023zWevueaaoLwVK1ZE5nvqqae6+tRte3jbo48+Gmm/9ei4YfthH8L7isn1119/du/evZE2Bq2jjH784x8P/NBnrxT2LSqvfDWfo7Zrnepqse70X6xsm2xFJTtWFL9BU6+yutlv9WIM+/qZz3ymLWsvTp31btu5tdDr2LSYtVpOnbueDfsVV07YP4tXXsdYmpjKF9U9ym/V1eI1otIerqBF6Pd1ipfeha+uw90mgzhrV8ZxNuweT+sEEJclk/Xhq9mwQbXKrIWkbq64pOeYxrXMfKmjfFcLz1pGapnEJU0/ML+juoxtvzhuFresjhMrb5DPiRMnptrd6h+3k9Ux7tjcvXt3sGtUd6wdU8qg7sQkqahjLC6m4XEBixcvHj4+wr6H6+VT7MM+hr8jhGEafE9EYNy4ccP57GQ6vCKnL9YFq26mslL4xBDlh7pQzc9+fLR9o2z3Y6/bPibmKjNOuG1/O/HpgepxKXziC+cJi06ZJ8TwcdpL2ML+J/leZNyKPMbiYhrmlyRPq8WYBGOpeRgsUyr+ahZ+9dVXD98L01VsrxNpP7WUGOgK++DBg8HudlXdj62s9rETnuzpxNpZbw3ISJPKrKO1XqZMmdLTZaunhEz1DnOwnaPW2TYfPsMXF7rvNUhKGze1GHXvOOyDlR8W6DhRsbz6LPIYSxLTsCiG/azad4SwahHzwF918Vg3mVoJ4QEhg7inE4XsJu0iGqSsfvYNnxjCJzCzleQE60sdrXUWbrFZPTo/8+gK7ywj72U7pnTiNmFPU+YgcdNTlySEaoWrbOtZ0DGkASdKGliURFR8O8bCx0Yanr7lRQh9i0gF/JEg6Aetk4P+4loJaaoSHl6uE8Ktt97q1PJUOfqxafSZhvKXmcJdPFEnrah1YX99qqN8tRZe2Meo72HR71XHqP3LXmfHqPywbt40PmUZN40e1cWejmnjqmM8bpRop5+9+Gfpa2fZdV5GCOsc3RzrpsdWaSi3kq5qNXy636QTlQ0v1z0QX59bKj8tJbkyt7z69K2OupgJC0TY187v4Tlt4VZxZz5fl1sjX4ddS9t7kUXcbOqIBrm0RoA6dfNLBCWAEua0Pg1XpuNLFr52mGzMIoNlGhPqbCuq+xk2T0w/wM4nqKQpLdxtFSeC4bmLaWynzRtu9XXu2xpaHqzSVblOYmlS1nUclEe4e858i6qPWo02N9LiHZXP13USIR2fSv0MNDE2inm/x6bux2p//WbU8rN5hLp4zEoEVb8sfJWdJiaEsIlRz6jO+lFbC0HdlknEUFfC1pLsdMPuW3Wu17L9yKO2ZbkubuJ4+IQafohz2rKzqqN1q6Ut3/JLCK2bTa3xOL/sHpb2ixMCs+nbp45Hu5etC5dB/I/jozonOTa1fzcbWbLrVk4SX7P0pSq2EMKqRMpTP3VVa60jiaGeSSjRsBO1fpT6LoGxZynqxxj+sYbnV+lJFXYFryorrz39pAgEOnHKf/NPvshvO6FK+K01lcafrOpoFx5iqvtB8lM+6uIiyYVI2Ocnn3wyWJSt1oT4wJ6+y6bipVjaiTP8OqOwDR+/m+92T1ktsX677rOIm3Wj61Fkipn8sz/xDR/vg/DMwtdByq/yvtwjrHL0PPBdrQqdZHQS1olHJ1KJhglHlIvqorLWiLbr5K5uNxtiLuHrTMoj23kniVyc//JB93n6SVnVcc6cOcODhtSSs3ur8ik8vzOJj7qAUex0gpb4hW2F9w+/4Da8vqzvOg4k0p1J3dp2AWPb1HofpPsxi7jpmJIASvDiGMtfxUP+JplGYfULf2bha9hek77TImxStHOsq7qddO9DP2T9oMNCp+9aJwFUnqguKnWzhu/h6EetE4LsdT6/MMdqBOVJHCTM8kFJn/JNfti6fnzIoo7iaK0z46oTrQRa9tMm2QvHzfa3OmtbPy1gs5PXp8Sw808iaEx03Mj3QUTQfB8kbvIpPFjHbEZ9WstePRL9pkF87bfMOuw3So99q0NFqAMEIAAB3wjYmyIk0HqUXWdrT0IpQdeAGvVEWItWIq59SMUQQAiL4UwpEIBAAwno/XxKEsHwPbwoFOEuX7XwO0Uzah/WZUOArtFsOGIFAhCAQCwBiVyvZK1B5avLE1t61dmX7QihL5HADwhAoHYEdG9ZSd2eGigTJYgSQN0XtEFiagnq3i2pOAJ0jRbHmpIgAIEGEgg/9syqb4OuOke6amCSBkORiiWAEBbLm9IgAIGGErA5hGoVmgBKEPWnVqBEkAEy5RwcCGE53CkVAhCAAAQ8IcA9Qk8CgRsQgAAEIFAOAYSwHO6UCgEIQAACnhBACD0JBG5AAAIQgEA5BBDCcrhTKgQgAAEIeEIAIfQkELgBAQhAAALlEEAIy+FOqRCAAAQg4AkBhNCTQOAGBCAAAQiUQwAhLIc7pUIAAhCAgCcEEEJPAoEbEIAABCBQDgGEsBzulAoBCEAAAp4QQAg9CQRuQAACEIBAOQTGllMspTaZwLFjJ9zp02cGRjB+/IVuzJjRgS3ZJJVPgJiUH4MoD7KMi35zslenhBDWKZoVqYtE8NSp05l6m7W9TJ1rqDFi4mfgicvIuNA1OpIJayAAAQhAoEEEEMIGBZuqQgACEIDASAII4UgmrIEABCAAgQYRQAgbFGyqCgEIQAACIwkghCOZsAYCEIAABBpEgFGjJQZ7y5YtbtOmTW7nzp2ZeHHs2DG3ceNGt3v3bjdhwgR39OhRN23aNLdo0SI3efLkTMrACAQgAIG6EUAIC4yohEritGvXLicR3LdvX2aly9b8+fPdrbfe6rZu3TpsV+XMnDnTrV+/3s2aNWt4PV8gAAEIQOAcAbpGCzoSJH7XXnutW7BggTt48KCbPn16piVLBNX6W758eZvdefPmuYULF7qlS5e6Q4cOtW1jAQIQgAAEnEMICzoKZsyY4Q4cOBB0g0qsrr766sxKXrdunVNrU6IXlW677bZg9cqVK6M2sw4CEIBAowkghDUI/+bNm4NaxHV96v6g/tQqpVVYg4BTBQhAIFMCCGGmOIs3JnFTGj9+fNfCbbDM9u3bu+ZjIwQgAIGmEUAIKx5xG3BjQhdXnXHjxgWb9uzZE5eF9RCAAAQaSQAhrHjYTdg0XaJbmjhxYrD58OHD3bKxDQIQgEDjCCCEFQ/58ePHE9XAWoQaVEOCAAQgAIG3CSCEb7NoxDfNYyRBAALNIfC9733XnT17tjkV7qOmTKjvA5pPu6QVNh9ahFm91HP06FFBKPR52WUX+xSWxvpCTPwL/R//8cvuT/7ka27VqlWO30p0fBDCaC6szZGA3nCdZRo1alTrTfXnRDFLu9jqnwAx6Z9d1nveccftTlOs7r///kAMs/79Ze1vGfYQwjKoZ1hmr0EyVpTdS+w1zcLy5/mpN9RnkXR1qxOuun3OnKHrJwumg9ogJoMSzGf/O+64wz322GPuq1/9qps79yMDF1I3MUUIBz4kyjVgg2CSepFUOJPa6yffsWMn3KlTp/vZtW0fdYeqJSgRfPXVN9q2sVAOAWJSDvdepSouetSifneD/lbGjh3jLr30ol5FVmp7tn1Ulap6PZy1aRG9nhhz5MiRoMI+tAjrQZ5aQMA/AidOnOjqlHpQSCMJIIQjmVRqzdSpUwN/ew2asa5RPZibBAEI1I/AN7+5y33wgze4PXv+pn6Vy7lGCGHOgPM2r4d5K/UaDWotRsuft1/YhwAEiiMgEVy69E73z//8T+7OO+9ADFOiRwhTAisru4TMnisa9kGPVrNWoT1uLbxd3yWSJoRxD+bu3IdlCECgGgRMBE+ePBk4/PrrxwMx/P73v1uNCnjgJULoQRB6uSAh08t1F7TeZahXLnWmxYsXB6u2bdvWuSlYNgFdtmxZ5HZWQgAC1STQKYJWi+uu+09u0qTLbZHPHgQQwh6A8tps4iT7cS05K3vv3r321e3fv3/4u31RK0/vItRcIWv52TaJ6IoVK4JW46JFi2w1nxCAQMUJxIng0NAct3btZ9zYsUwKSBpiSCUlNWA+CdL1118fWNH3cJo7d+7wa5T05voNGzaENzvd11P3p0QuTsxWr17tNJVCttRC1D4S0E2bNrkom20FsAABCFSKACKYbbhGtSYjMxM5W6ax1kwA46Yw9Noeazi0QTbU2pRoqhwJYq9XNIV2L+Tra6+9meE8wtFOE/QHnRtVSMUbUMi5eYTEJM9Q9yOCWcZlbA3nEdIizPOI7bAdJ4CWrdd2y9ftUzYYENONENsgUF0C/YhgdWtbnOfcIyyONSVBAAIQ6JsAItg3up47IoQ9EZEBAhCAQLkEEMF8+SOE+fLFOgQgAIGBCBw48O1gsrzNEzRjjA41EoN/IoSDM8QCBCAAgdwI/MzPvMd94AMz2+wjgm04Bl5ACAdGiAEIQAAC+RHQfMA1a/57axDcrwSFIILZs2bUaPZMsQgBCEAgUwImhr/wC7/oPvKRW5ksnyld5xDCjIFiDgIQgEAeBCSG8+b9eh6mG2+TrtHGHwIAgAAEfCHw8ssvtR4O8S++uNMYPxDCxoSaikIAAj4T0BSJu+5a2HqT/EcRw4IDhRAWDJziIAABCHQSCM8T/O53v4MYdgLKeRkhzBkw5iEAAQh0IxAWQcsnMfzSl75oi3zmTAAhzBkw5iEAAQjEEYgSQeWdNWtOq5t0SdxurM+YAEKYMVDMQQACEEhCoJsIrlnD+wSTMMwqD0KYFUnsQAACEEhIABFMCKqgbAhhQaApBgIQgIAIIIL+HQcIoX8xwSMIQKCmBBBBPwOLEPoZF7yCAARqRgAR9DegCKG/scEzCECgJgQQQb8DiRD6HR+8gwAEKk7g7Nmz7vd+b73rfJ+gpkgwOtSP4CKEfsQBLyAAgZoSGDVqlFu//n+4K664criGiOAwCi++IIRehAEnIACBOhO49NJL3caNTwZiiAj6F2lew+RfTPAIAhCoIQGJ4ebNf+Auuuhi3ifoWXwRQs8CgjsQgEB9CYwfP6G+latwzegarXDwcB0CEPCLgEaH3nffJ9xbb73ll2N405UALcKueNgIAQhAIBmB8BSJEydOuN/+7Yfceeedl2xncpVKgBZhqfgpHAIQqAOBsAiqPjt37nCf+tQyWoYVCS5CWJFA4SYEIOAngU4RNC/HjBnjNHWC5D8BhND/GOEhBCDgKYE4ERwamuPWruVVSp6GbYRbCOEIJKyAAAQg0JsAItibUVVyIIRViRR+QgAC3hBABL0JRSaOIISZYMQIBCDQFAKIYP0ijRDWL6bUCAIQyIkAIpgT2JLNIoQlB4DiIQCBahBABKsRp368RAj7ocY+EIBAowi8+OJut3TpnSNepcTo0HocBghhPeJILSAAgRwJjB8/3l1wwYVtJSCCbTgqvYAQVjp8OA8BCBRB4KqrprZeo/S4s4dmI4JFUC+uDJ41WhxrSoIABCpMwMTwa1972n3ykyt4lVKFY9np+qizrdS5kmUI5EngtdfedKdOnR64iMsuu9iNGTPanT59xr366hsD28PA4ASIyeAM87CQZVzGjh3jLr30ojzcLM0mXaOloadgCEDARwI//OE/+OgWPuVIACHMES6mIQCBahHQFIkPf3i2e/LJx6rlON4ORAAhHAgfO0MAAnUhEJ4n+LnPrUMM6xLYBPVACBNAIgsEIFBvAmERtJpKDP/iL/63LfJZYwIIYY2DS9UgAIHeBKJEUHvNmjXHXXfdf+5tgByVJ4AQVj6EVAACEOiXQDcRXLOG9wn2y7Vq+yGEVYsY/kIAApkQQAQzwVgLIwhhLcJIJSAAgTQEEME0tOqfFyGsf4ypIQQgECKACIZg8DUggBByIEAAAo0hgAg2JtSpKooQpsJFZghAoKoEEMGqRi5/vxHC/BlTAgQgUDKBV1/9F3f33b854n2CmiLB6NCSg+NB8QihB0HABQhAIF8Cl132b9ynP32/GzVq1HBBiOAwisZ/4TVMjT8EAACBZhD44Ad/NajoihWfdENDs2kJNiPsiWqJECbCRCYIQKAOBCSGP/VT/9b9/M+/l/cJ1iGgGdUBIcwIJGYgAIFqELjmmv9QDUfxsjAC3CMsDDUFQQACRRDQ6NAXX9xdRFGUURMCCGFNAkk1IAAB52yKxJIlv4EYckAkJoAQJkZFRghAwGcCJoInT54Mpkkghj5Hyy/fEEK/4oE3EIBAHwTCImi7SxD/9E//2Bb5hEAsAYQwFg0bIACBKhCIEkH5PTQ0x61c+WAVqoCPJRNACEsOAMVDAAL9E+gmgmvX8j7B/sk2a0+EsFnxprYQqA0BRLA2oSy9Ighh6SHAAQhAIC0BRDAtMfJ3I1B5IXz44YfdI4880q2Oqbbt37/fXXXVVW7KlCmp9iMzBCBQDAFEsBjOTSql8kK4adMmp79e6f7773cPPPBAr2zB9rNnzzr9kSAAAb8IIIJ+xaMu3lReCJMGYsuWLU5/JAhAoJoEEMFqxq0KXjdGCKsQDHyEAATiCbzwwtdHvE9QUyQYHRrPjC3JCCCEyTiRCwIQKJnAPfd8wt16668Pe4EIDqPgy4AEePvEgADZHQIQKIaAXqr7qU+tCAo7evQILcFisDeiFISwEWGmkhCoBwETw9OnT/M+wXqE1Ita0DXqRRhwAgIQCBM4depUeLHtu8Rw7Fiu4dugsDAQAYRwIHzsDAEIZE1Ao0NvvfVX3Q9+cDhr09iDQCQBhDASCyshAIEyCNgUie9+9+/cokUfRQzLCEIDy0QIGxh0qgwBHwmYCOr1SUr/8A8/CMTwn/7p//noLj7ViABCWKNgUhUIVJVApwhaPX72Z3/OXXrpZbbIJwRyIVCbO8433nhjIkC98h09ejSRHTJBAALZEIgTwVmz5rg1a3iVUjaUsdKNQG2E8NChQ93q6TTSTKlXvq5G2AgBCGRKABHMFCfG+iRQCyHkAdl9Rp/dIFAiAUSwRPgU3Uag8kL47W9/u61CLEAAAv4TQAT9j1GTPKy8EFYxWMeOHXMbN250u3fvdhMmTHC6Lzlt2rTWCLlFbvLkyX1XaebMmcF7FOfMmeOmTp0a2FJZ6g7eu3dv8PaNZcuWuRkzZvRdBjtCYFACiOCgBNk/awIIYdZEe9jbt2+fmz9/fmvC8K1u69atw7n1iigJ2fr1692sWbOG16f5IkHdsWNH8Be137x58xDBKDCsK4wAIlgYagpKQQAhTAEri6wSQbX+li9f3mZOInXw4EG3dOlSt3PnzoFahm2GWwvjx493Dz74YN8C22mPZQj0Q+Cv//qvWsf3nSNepcTo0H5osk+WBBothK+//nrQbThu3Dg3adKkLLlG2lq3bp1TV6VELyrddtttbvPmzW7lypXuiSeeiMrSc532U1eoRPXyyy8PBJWu0J7YyFAAgfe859+7q66a6r71rf87XBoiOIyCLyUSaKQQPvfcc+6hhx5qm0qhVpMEavHixe6SSy7JJSQSOaW4rk/dH9Tfrl27At/6uV+I6OUSOoxmQODiiy9xX/jCo+5jH1sciCEimAFUTGRCoPJPllHr56qrrnJDQ0PumWee6Qnlt37rt4LuR+2naRf2p/trmzZtctdff33r+YY/6GknbQaJm5IEt1sy8du+fXu3bGyDQCUJvC2GS5ksX8kI1tPpyrcIJV4SsyNHjrjZs2d3jZJEMCwwEs/rrrsu2OeFF14IBplIEOfOnev+7M/+LNOWoQbJKJnQBQsR/6mbVmnPnj0RW1kFgeoTkBguWnRn9StCDWpDoPJC+PTTTwdPjdFgkG5dmpqqYCKop8x87WtfC6YaWCQ1ilOttttvvz24j6eu0wceeMA2D/xpwqbpEt3SxIkTg82HDx/ulo1tEPCawLe+9VfuXe+63F144YVe+4lzEBCBSneNStws9XqG6FNPPRVklQiuWrWqTQTNhu6v3XvvvUELU9MZNJgmq3T8+PFEpqxFqEE1/Sb5vmTJkuBPrdsFCxY4a5H2a5P9IJCUwNe//nV3yy23tG5B3OVOnDiRdDfyQaA0ApUWQt3nU5oyZUpPgJpfJxHUPTq1/uKSJrVbsvt6tlzkp7po+0kSPXW/btiwIfjTXEUJvARRAkmCQJ4EvvGNbwQXXnqV0osv7kIM84SN7cwIVFoIrdXU675buOU4ffr0nvBMWE1oe+6QIENaYbO6JTA9nEXzE1evXj1i0rzEXWKoaRllivuwo3ypJQFNlr/99gVtrUCJ4Ve+8j9rWV8qVR8Clb9HqFD06nbUQBhLNjjGlqM+7T5d1Daf13Wbe6gpGxJBiaEm7JeZxo/P5r7R6NHn3iiiz8suu7jMKjW+bLUEoybL33TTTa31v+nGjq3Fqaaycea30j10lT46rSWo52h2S+oWtdRrZKnyWUvQ7Nu+g3z2GiRjtk3Ue02zsPxJP21+oeomQbTlpPtnmW/MmGw7ItTlPWbMOVHM0k9sJSOge4J33HH7iCfGfOhDH3Kf//znEcFkGAvJxW8lGnOlhdBO5upGfP75590NN9wwopbqFtXJXweAujy7jSy1nS1/lkJog2CsjF6fSYWzlx3bHrangTPGzrYX+Xn69JlMitNVruKq6TNnzpzNxCZG0hFQS/D2Vneo7gmGk0Twc59b34rPaJdVvMP2+Z6OQNa/lawvZtPVJvvclRZCiYtGi+pJMffdd18gdO9+97uHKWnUp7oCLekRZr1SuPWoifpZJetutdZmnF3Nh1RK0yKUzQWtQTJK6h7tJeC9fAgM5fjfsWMn3KlTpwcuQd2haglKBF999Y2B7WEgHYG4B2ibCB47JnFsF8h0JZA7KwJZ/lbGjh3jLr30oqxc88JOpYVQBPXwagmhWoV6e4Mek6ZXEOlZm5pjqPXWGtSQ7l5Jr0dS/qxbTPJJqdegGesa1cCXpEldnSZumisZHvma1Ab5IJCGQJwI6p6gukPVEkQE0xAlb5kEKi+Eav08/vjjwUR4gQxPEVCXmZLyfPGLXwy+d/tP3aj79+8PhDDuwdjd9u+2zYS112hQEzTL382mbbNuT7Ui455jGhZgE2Xbn08IpCEQJ4JDQ3NarxHbENwTpDs0DVHylk0g21ELJdVGoqF7hOomtWeH6lPCsHDhwmBbr3uDaonpFUhKupfYa4J+2qpKjE2A4ia3SyRNCKMETduipj/Irhi89NJLsd2i9lQdMUkyYCht/cjfDALdRHDt2s8wMKYZh0Htaln5FqFFREKjSeRKEgzdk0szQEV5n3zyyeBN7mm6Ja38JJ96s4XEdtu2bcOiGN7PRE5vke9MEkl1/SpJ3MPvM1Td5b8E1sS2c/9nn302WKVH0aW5/9hph+XmEnjrrbfc2rUPjBgYo5YgItjc46IONa9Fi7AzECYMnet7LaslqKfO2IT6XvnTblcrT12ueh2TtfzMhoRuxYoVgZBF3eMLTxFR921nWrNmTSCyUa1NDRjSepUd1dLstMUyBKIInHfeee53f3eje+c73zm8GREcRsGXChMY1epCZNx5wQHUC3o1kEctRHVpSuT0Fg0JsLVqo1zSY9IkoOvXr48czCMxnT9/vtM9Q7vHqJag9vHpDfWvvfZmhqNGzw3PZ9Ro1BGTz7q///vvtXol/ot773uvHdESPDc6kZjkQ75/q1nGZezY+o0aRQj7P7YG2lOiZaM91VUp4VJLNosku9YytPuHWdjNygZCmBXJ8uz84z/+0P3ET7xzxD3BLE+45dWufiVnGReEsH7HBzUqgQBCWAL0gorM8oRbkMuNKCbLuNRRCCs/WCbLSe/hX4TmEkbdiwvn4TsE6kpAo0Nffvkld9ddS4LpRHWtJ/WCgAhUXgi5xcmBDIFsCYSnSJw6dar1Xst7EMNsEWPNMwKVF0LxVOtNaWhoyE2aNCn4zn8QgEB6AmER1N5PPLEpMIIYpmfJHtUhUHkh1Pw5eyyZnhOqkZd6pqgmjfeaRF+dMOEpBPIn0CmCVuIPf/iD1vNcz7Se6zrGVvEJgVoRqMWoUT0a7Stf+UrwzFFFx1qIGokpUYx6K0WtolixyjBYxr+AxYngrFlz3Jo1yZ8Yk+WgDP8oVdejLONSx8EytRBCOzzVMtRTW/S8URvoYqKoieRqJWb96DQrm8/kBBDC5KyKyJmVCMrXLE+4RdS9KWVkGReEsEJHzaHWJPKnnnrKPfPMM8EbKOS6RNGetamnrOQ14rRCmEpxFSEsBXtkoVmKoArI8oQb6TAr+yKQZVwQwr5CUP5OcV2nmsCurlM9Vo37icXFCSEsjnW3krIWQZWV5Qm3m+9sS0cgy7gghOnYe5lb3aZRXafcTywuXAhhcazjSspDBFVWlifcON9Zn55AlnFBCNPz93YPdZ3q1UR65qe+K4W7TvUGCFqJ+YQPIcyHa1KreYmgys/yhJu0PuTrTSDLuCCEvXlXMoeEUFMvNNDGXsyrl/1Onz69kvXx3WmEsLwI5SmCqlWWJ9zyKNWv5CzjUkchrOVrmNIexq+88krwBgh9kiBQZwInT55wnW+PTztFos58qFszCVR+Qn2/YXvuueeCFqBagkr2qDYNoNGI0quvvrpf0+wHAW8JvP/9v+wefni9u/fepa1XYb3Vej9lunmC3lYMxyAwAIFGCaFafJpSoXuDeg2SkgRQUyo0z1AjSPN6Ke8AMWJXCGRKwMTw+ed3uAceWDviVUqZFoYxCFSAQO2F8PXXXw9GiWqkqA2KsdafRoqq9adnlJIg0CQCEkP9kSAAgRq8fSIuiOr6lPjpJbVKJn7W9SkB1HNKSRCoK4GjR4+4CRMm1rV61AsCmRGoVYvQuj41JUIpLH5q9c2ZM4euz8wOHQz5TECjQ++99+Nu1ar/5mbOpMfD51jhW/kEKi+Ehw8fDqY+WNeniZ/u++mJMYhf+QcZHhRLIDxF4pOfvNf9zu84xLDYEFBaxQhUXghnzpwZTISXANpzRDXoheeIVuxIxN1MCIRFUAb1Yl2J4WOPvcv93M/9QiZlYAQCdSNQeSG0gOipMJryoBbiQw89ZKsH+nzssccG2p+dIVAkgU4RtLKvv/5GN3Uq04GMB58Q6CRQGyFUxWxgTGcl0y6rdWmvb0q7L/khUAaBOBEcGprj1q5N/j7BMnynTAiUTaDyQjhp0iREq+yjiPJLJYAIloqfwmtAoPJCuHPnzhqEgSpAoD8CiGB/3NgLAmECPGs0TIPvEKgQAUSwQsHCVa8JIIRehwfnIBBNABGM5sJaCPRDACHshxr7QKBEAohgifApupYEEMJahpVK1ZXA4cOH3NKld7qTJ0+2VZHRoW04WIBAKgIIYSpcZIZAuQTe/e5J7qMfvaPNCUSwDQcLEEhNACFMjYwdIFAeAc1vvfPOj7vFi+8KnEAEy4sFJdeHQOWnT9QnFNQEAskImBi+5z1XuF/+5Zm8TzAZNnJBIJYAQhiLhg0Q8JeAxPDGG2f56yCeQaBCBOgarVCwcLVZBF58cbf7u7/7TrMqTW0hUAIBhLAE6BQJgV4ENEViyZLfcIsWfRQx7AWL7RAYkABCOCBAdodA1gTC8wRfffVfEMOsAWMPAh0EEMIOICxCoEwCYRE0PySGTz/9v2yRTwhAIGMCCGHGQDEHgX4JRImgbM2aNcctX/5f+zXLfhCAQA8CCGEPQGyGQBEEuongmjW8T7CIGFBGcwkghM2NPTX3hAAi6EkgcKOxBBDCxoaeivtAABH0IQr40HQCCGHTjwDqXxoBRLA09BQMgTYCCGEbDhYgUAwBRLAYzpQCgSQEEMIklMgDgQwJIIIZwsQUBDIggBBmABETEEhD4KtffWrE+wQ1RYLRoWkokhcC2RFACLNjiSUIJCKwdu0690u/9P7hvIjgMAq+QKAUArx9ohTsFNpkAhdccIF7+OEN7t57l7hLLrmYlmCTDwbq7gUBhNCLMOBE0whIDB955PNu9OjRvE+wacGnvt4RQAi9CwkO1YnA2bNnnd4dGJXOP//8qNWsgwAECibAPcKCgVNccwhodKheo3T8+PHmVJqaQqCCBBDCCgYNl/0nYFMk/vIv/4+78847EEP/Q4aHDSaAEDY4+FQ9HwImgidPngwK2Lv3bwIxfOON1/MpEKsQgMBABBDCgfCxMwTaCXSKoG1997snuQsuuNAW+YQABDwiwGAZj4KBK9UmECeCQ0Nz3Nq12b5KSeNvRo8eFTsQxweSY8dyne1DHDp9UFxaY7haf2fdmTOtLySHEHIQQCADAkWIoEafnnfeGHf++WPcmDF+iozEWUmfF198QQZkMZEFgbi4SAx//OPT7q23TrvTp89kUVQlbSCElQwbTvtEoAgRlKhYC0uCeO7vHIWY2RmlIzJ/S3cEB9oIhC+iJIQXXDCq9Tc2aCG++eaP3alTzRNEhLDtEGEBAukIFCGCl1xyQdAC1OT71r8gWdeWFvTdlzRmzNtzJul28yUqrnX8vB0XiZ+SLqbUUtSfVileuuB6442TjRNDhNCfYxVPKkagSBHUVbxafjpZ+SwwY8a8HUSf/Xzby2Z8i45LuyBKLE+fdo0UQz9vNDTj2KSWFSZQhAhedNH5QUvQRFBdVohLhQ8aT11XC1H3B9UqlBiqpaiWoa9d7nlgRAjzoIrNWhN46aUX3dKld454lVLWo0M1MEbdVjohNfG+Ta0PIg8rFxZDuafjrykJIWxKpKlnZgQuv/zfuXe96yfb7OUhgipAV+e0AttQs5AjARs5quPuvPOac+cMIczxoMJ0PQn85E/+lNu8+Q/c5MmXBxXMWgRlVFfjOhmpNWiDG+pJk1r5RkBdpDruNOpXn01IzZH8JkSTOhZGwMTwy1/+oluy5N7MX6V0TgjPVUcnJhIEiiKgCy91yasnYuzYMcEcw6LKLqucUa1K8zMri35Dy33ttTdb97xaw9MGTJdddnEwmETdOa+++saA1vzZXVfh48e/I6ibvLLuKn88jPfEWrI6rWiSNskPAmniouNPA7R0X/rEibda98JPtVVC4njppRe1rav6Ai3CqkcQ/3Ml8J3v/K376Z9+T+sKubi7COoStZT0OvXqq6cGu0yaNNk9++x2273r5zPPPO1Wr14V5Nm2bXurq3dy1/y28dixY27Hju3uueeec4cPH279HXLjxo13EyZMaAn4ODdjxozg79pr/6Pt0vZpvratTLCwYMHt7p577k2QkyxZENBhGD4Ws7Dpqw2E0NfI4FfpBGyKxK/8yofcffetKlQMS698jAOPP/6Y++xnHxmx9fjxY61XTR0L1u/fv791D3Wzu+GGG90jj3x2RN4mrHjuuR3umWeecfv27Qu4XHXVFDdr1ix3++13ZFb9vMpoYh8hQpjZYYmhOhEwEdSrlLZufSaoWpPFUK3AhQvvcK+8sj9goZbnLbfc4t73vulu0qRJ7ujRo0Hr8G//9ttuy5Yt7tChQ8H6bsdEmtZrNzu+bbvnnrvd888/1+aWuOlP4rhly9OtlvP4tu1pF4ooI61PVc5fXH9PlSnhe6MIhEXQKi4x/KM/2mqLjfucN+/WYRH8tV+7Jeh+VetmypQpwUld3arTp09vieVCt3PnTrdq1apWC2h24zitWvXAsAiuXHm/27NnX/D36KObWxcGk4NuZLEcJBVRxiD+VXFfhLCKUcPn3AhEiaAK0xSJm276cG7l+mxYrQ/dB1TSfbr7738g+N7tv3nz5gUi2S1P3bbt3r3bffWr53oPJHy33PK24OkiQS1B3UsVS92f7ScVUUY/flV9H4Sw6hHE/8wIdBPBrN8nmJnTORvSide6+XSfi8Eq8cBN3NRdLOHrTOoOVYtZ6bOf7e/eaRFldPrdhGWEsAlRpo49CSCC0YjCA2MeeODcCNPonKy1C4Ybb7wxFoYGEClpYJHuo6ZNRZSR1qc65EcI6xBF6jAQAUQwGp9O1DY4Rq1B3Q8kRRPQSFlLU6dOs68jPsNTVL75zd0jtndbUUQZ3cqv8zaEsM7RpW49CSCC8YjCJ2oN/SfFE9i3b+/wxl4XDBo0o6Ru5zSpiDLS+FOnvAhhnaJJXVIRQAS74wqfqNUiJMUTOH78ePzGji2abqKUZp+0+fstI3Csgf8xj7CBQafKziGCvY+C8InaTqy990qeQ6Mnkz5lRlM2koxWTV56tjmPHDkSGNSo0F5p3LhxQRbNvUyTiigjjT91yosQ1ima1CURgTfffMN96lP35v4+wUTOeJwpfKIO39vy1WVN+r/uupGjNQf1V1MhokaBhu2GLxrC67t9T7tP2vwqu599uvlc120IYV0jS71iCVx00cXuoYc2uI99bHHrocI/CvLl8SqlWAcqskHPDrUkkRn0aShmyz6zfrKM/NMzU7M++fe656f6hC8arH5ZfxZRRtY+V8UeQliVSOFnpgR+8RevdV/4wqOBGL7//R9wTZ0n2A2qdeEpjx6unUQQutkrYltZLdfwRUPSeob5JtmniDKS+FHHPAyWqWNUqVMiAhLDL31pCyIYQ2vatLenAYRHkMZkb/RqEzV78Hg3GNZiTStsRZTRze86b0MI6xxd6taTwHvec0XmL9XtWaiHGaJOyjb5W+5u357s1U4eVq0QlyZOnDhcjrqRuyXr4jRh65Y3vK2IMsLlNek7QtikaDewrhod+od/+LXa19xGK9ozQZNUOPxkk6j7f+pmtGkTmlgfzp/EfpPyhEfVqhu5W7KHFERdfHTbr4gyupVf523cIywhurpi3LhxYzChVj8GXSGqG2rRokWJX44a53aetuPK9HV9eIqEXnB7882/5qurA/v1vve9b/iZoBKsJPfKXnnllaBcE7soJ+6++x63ePHCYNO9997jnn763EOlo/KWvU7HflmjRvV8UUua+B53PzV8MRGXx+x0fhZRRmeZTVlGCAuOtF7UOX/+fHfrrbe23nP39mt99A63mTNnuvXr1wcv8OzHrTxt9+NPmfuERVB+rFp1X+BOXcVQT36x51Dqwcy9Ho6tE7Ld9+v21BhNG1AXqWyrJfPIIw/3tF1W3MscNaqydUEhRnoQQfjNE2EexlzrwsIWzhP3vYgy4squ+3qEsOAISwTV+lu+fHlbyXptzcGDB93SpUuD97kluaJvM9BayNN2Z1k+L3eKoPn64ou73Yc//BE3atQoW1WbzxtvHBo+ET/xxOPBibjbMaTWnZK6VHu9NV1vmZ89e1bw+iDZVkoitCtWrHSbNm0O8hf1X7c65+2DLigkhLpoUOs0qrv58cfP8ZNo9uNrEWXkzclH+9wjLDAq69atC34gEr2odNtttwWrV65cGbW567o8bXct2LONcSKoeYJr1vz3WoqgheDhhx+xr63uzEUu/JBm26ATtN4vaPepwvtYnqhPvUvPulAlhhLGxx9/bPi+oeyqlalBNUuWLGm1Im8IWpxaH5d0S0A+pvmLs+XDej39xu7VLlx4R/BbD/ulF+raPVx1OUelcy3uu4OWd9T2LMqIstv0daNa907ONh1CUfW/8sorg6IOHDgQW6S6R3VC0Vu+01wx5mk71tk+N7z22pvu1KnTfe799m6XXXaxGzNmtDt9+ox79dU3avPYtNGjR7VOqBcGddPP88yZ5D9RiYpOwjaMX92a6oHQCEVt27Fjx/A2iaBakmmSTtTWKuy1n17iG9VyTPpYtSj7euO7z0ndonZPVaI4NDQUsH/++eeHRTCOi+plLe9uDxvop4zzzhsTXATqeHrrrd6/vbFjz/2uTpw41XroxFttyMeOHeMuvfSitnVVX6BrtKAI7tq1Kygpqrsk7ILEz66sNXgmScrTdpLyfcjTrSXYpMnyGoCxffuO1otfHwnelq5uOrt3aHHSvamVK+9PdaFl+0rYFi5cFNhW608jJE10deKfPHmSmzFjRksMFrt3vONi260xn+feRP+MU9ezWn/2xnoBEJ+777479v5hUkhFlJHUl7rkQwgLiqQGsij1auXZ3KI9e/Yk9ixP24mdKDHjN77xjda91Tt5dui/xkAXW3pAtf6ee25HIFbapOH3EsFeF2O9Qqn9dV8x6t5ikpaH7626XvXvtV0XI88+uz0YNCMx1AR6dSv3el6p7Gq/JGmQMpLYb1oehLCgiJuw9Zo7ZJNme81FCrudp+1wOT5+//rXv946IS9ABGOCk7brM8YMq/sgcE743p5W0YeJnrsUUUZPJ2qQgcEyBQXRHqvUqzhrEXYbZNBpI0/bnWX5tKyW4IIFiKBPMcEXCFSRAELoadTsMUx5uJen7Tz8jbKpe4JqCZ44caJtM2+RaMPBAgQgkIAAQpgAUhZZ0opPmhZhnrazqHseNr7//e/RHZoHWGxCoIEEuEfYwKCXXeXx4y8c2IW77tKoxLGtASH3B7Zuuumm1lN5NtTmAdqaQnEujWpNoxgYVykGNHCG5B+BpHHRNImLLx7tLrrofP8qkbFHCGHGQOPM9RokY/vZ/b40I/vytG1+ZfmpuX9ZpIULFwZzo15++WW3YUN9RDALNmXbqOPTe8pmmkX5aeOSNn8WPpZhAyEsiLoNgklaXFJxk708bSf1N00+TYDPIqnVdMcddwQDZvRYiKzsZuHboDbebhEOaqn4/XXy5DkdxXPvVWI/cVEc9dvqTFldzHbaLWsZISyIvE2L0GT5bunIkSPB5jQtwjxtd/O1323Hjp3I8MkyOum64Mky/frj2346Yan7WCebtE+WKbsu4W63JE8wKdvfppSfJi6twy849k6dOhM8VebkyVNtmOr4ZJls+qjaMLEQRWDq1KnB6l4DW6xrNPx28Ch74XV52g6Xw/diCIRbU03pmiqGLKUkI2D3p11kazCZjWrlQggLipceO6XUazSotRgtfxL3LG8etpOUT57sCZzrkjrbugeavW0sQqAbAR1z1h165kw2tzG6lefDNoSwoCjo0WrWcrNHonUWLSEzIYx6R5y22XNFw/tmYTtsj+/lE1C3oj1sm1Zh+fFokge6P233BtU92oSEEBYYZT2IWGnbtm2RpZrILVu2bMR2iaTeTLGg9SQVvXKpMw1iu9MWy+UT+PGPz70hQFfmVR44Uz5JPEhDwHogJIRZvCEmTdll5kUIC6SvVp7eRbh58+bhlp8VL6FbsWJF0GqMeuvE3r17LWvke+YGsT1smC/eENAI2HD3KK1Cb0JTa0dGj9YArXNdoz/+cfsgmTpXfMwDrVTnCvpWtw984APBY8H08l2d3C688EL353/+5+4Tn/iEe+973+uefPLJSJfV/al8J0+edJ/+9Kcj32LRr+3IAnNcqfebWbffIMW84x3nB60lCcaPftT+zrRB7Pq0rw1THzPm3OhYn3yL8kX+mmhnEeOoMliXnkCSuNixposwxa7zPYRWqsTyHe84zxZr8cmLeUsKo1qA6grVfT9NldCAF4ldFilP21n4l9eLebPwzTcbF154nrvggrGB4J+7d6Or9XQv6y2yThqmLyGUj0yfKJJ897Li4qJY6bhSl6hagiaCr79+IliOsjq29cSZur2YFyGMijTrciWAEKbDa2IYPmmZBZ28fErnTqrnhJAWoT+RiYqLxE9Jx5BdXClm3URQ+esohEyoV2RJEPCYgHUln3/+mGCis1yVKNqJzC/X357voZMryRcCI+NiAmifasHrWGti2BBCX45T/IBAFwIauKA/Xdmrm8u6unwTxPADwmkRdglowZs646KLFAmeRoZKAJsyTSIOO0IYR4b1EPCQgMRFj7zqfOyVL65edtnFrVbrqGCwxdGjP/LFrcb7QVy6HwJMn+jOh60QgAAEIFBzAghhzQNM9SAAAQhAoDsBhLA7H7ZCAAIQgEDNCSCENQ8w1YMABCAAge4EGCzTnQ9bcyBgT7DI0rTmNpH8IkBM/IqHeTNoXPL4/ZpvZX0yob4s8pQLAQhAAAJeEKBr1Isw4AQEIAABCJRFACEsizzlQgACEICAFwQQQi/CgBMQgAAEIFAWAYSwLPKUCwEIQAACXhBACL0IA05AAAIQgEBZBBDCsshTLgQgAAEIeEEAIfQiDDgBAQhAAAJlEUAIyyJPuRCAAAQg4AUBhNCLMOAEBCAAAQiURQAhLIs85UIAAhCAgBcEEEIvwoATEIAABCBQFgGEsCzylAsBCEAAAl4QQAi9CANOQAACEIBAWQQQwrLIUy4EIAABCHhBACH0Igw4AQEIQAACZRFACMsiT7kQgAAEIOAFAYTQizDgBAQgAAEIlEUAISyLPOVCAAIQgIAXBBBCL8KAExCAAAQgUBYBhLAs8pQLAQhAAAJeEEAIvQgDTkAAAhCAQFkEEMKyyFMuBCAAAQh4QQAh9CIMOAEBCEAAAmURQAjLIk+5EIAABCDgBQGE0Isw4AQEIAABCJRFYGxZBVMuBI4dO+Y2btzodu/e7SZMmOCOHj3qpk2b5hYtWuQmT548EKA8bQ/kWAV2zovdzJkz3ZQpU9ycOXPc1KlTgxirrEOHDrm9e/e6LVu2uGXLlrkZM2ZUgFJ5LorTpk2b3M6dOzNxIq94Z+JcUUbOkiBQAoHWie/sNddcc/Yzn/lMW+lPPfXU2SuuuOLss88+27Y+zUKettP4UcW8ebJTvBXbuL8VK1ZUEVnuPrcuEM8ePHjwrH4bN9988zC/LArOM95Z+FeUjVEqqCjRpRwIGIFrr702aP098cQTtmr4c926dW7z5s3BFW8/LcM8bQ87WdMvebKTbbU+OtP48ePdgw8+6GbNmtW5qfHLu3btcgsWLAhaz0NDQwEP/TaUDhw4EHwO8l+e8R7Er6L35R5h0cQpz0nodEKcN29eJI3bbrstWL9y5crI7d1W5mm7W7l12FYEO134rF692i1cuDD41PJLL72ECMYcQOomluCpG3T58uXu6quvjsmZfnUR8U7vVTl70CIsh3ujS73yyiuD+ne7otX9JN070gkgTaswT9t1D1re7NT6kOiR+iewfft2t3Tp0sBAt99PkhLyjncSH3zJQ4vQl0g0xA919SipO6xbMvHTDz9pytN2Uh+qmg92VY1cf34T73ZuCGE7D5ZyJrBv376gBBO6uOLGjRsXbNqzZ09clhHr87Q9orCarYBdzQLaozrEux0QQtjOg6WcCZiwabpEtzRx4sRg8+HDh7tla9uWp+22gmq4ALsaBrVLlYh3OxzmEbbzYClnAsePH09UgrUIo0YZxhnI03ZcmXVZXyQ7zYN74YUXAnS60NFFkeYPam4hqRgCRca7mBoNVgotwsH4sXfOBDTJPq+Up+28fPbFbr/sbCrAhg0bnP62bt0aTKCfO3duMKHel/rhRzuBfuPdbsXfJYTQ39jU0rO0P6g0LcI8bdcyGKFKFcFOTw3S1InOJ8foSUJap+kyNogj5BpfcyBQRLxzcDs3kwhhbmgxDAEIhAlozmDcICmbTN/P3NFwGXyHQD8EEMJ+qLFP3wR6DZIxw3YPo9c0C8uvzzxth8up4/ey2VkrUXNHaRXmf4SVHe/8a5iuBIQwHS9yD0jABsEkNZP0Byt7edpO6m9V85XNLhxnG9pfVZZV8LvsePvGCCH0LSI198emRejKv1s6cuRIsDlNizBP2918rcO2PNkp1npSkD0tqBevXsdGr/3Z3ptAnvHuXbp/ORBC/2JSa49siHyvm/XWNaoBFklTnraT+lDVfHmyU1enxE1/aZ4UVFWWVfA7z3hXof6dPiKEnURYzpWA3QvqNRrUWgWWP4lTljcP20nKr3KePNlZt6da9zYoppNV+MLITtKdSvhzGgAADBpJREFUeVjOjkCe8c7Oy+IsIYTFsaakFgGNGrQTXdy9IAmZCWHUiVPbogZUZGG7qUHKgl1cXBRvnXj1wO24UaPWUpRYzp49u6lhyLzecTHJIt6ZO1uiQYSwRPhNLXrx4sVB1bdt2xaJwEROTxvpTBJJ3WvSxGy9RqYzDWK701bTlgdh1y0uOulqcEbchY84t17EHODWewnT3BduWozS1LdbTGRnkHin8aMKeRHCKkSpZj6qlad3EeoFo9bysyrqx9t6U3nQatRE687UeqP28Kr9+/cPf7cvg9g2G039HIRdr7isWbMmeH1QlBhq7qDW65iI6gFoajyi6m0XidoWxTK8T6+YDBLvcDl1+M77COsQxYrWQS26p59+OrgyVdeZfribNm1yU6ZMCR6/FVctPY5LArp+/foRTymxffq1bfs3+bNfdr3ioouc+fPnB/M97R6VWoKKJW+ojz7ixOz6668PNup7Z7LW8/Tp0yN/M71iInv9xrvTlyovI4RVjl4NfNeP20YV6ketE2TcfaS01c3TdlpfqpY/T3aKt7Vm7P5h1fgU6a8JoIleZ9m9tnfmj1rOM95R5fm2DiH0LSL4AwEIQAAChRLgHmGhuCkMAhCAAAR8I4AQ+hYR/IEABCAAgUIJIISF4qYwCEAAAhDwjQBC6FtE8AcCEIAABAolgBAWipvCIAABCEDANwIIoW8RwR8IQAACECiUAEJYKG4KgwAEIAAB3wgghL5FBH8gAAEIQKBQAghhobgpDAIQgAAEfCOAEPoWEfyBAAQgAIFCCSCEheKmMAhAAAIQ8I0AQuhbRPAHAhCAAAQKJYAQFoqbwiAAAQhAwDcCCKFvEcEfCEAAAhAolABCWChuCoMABCAAAd8IIIS+RQR/IAABCECgUAIIYaG4KQwCEIAABHwjgBD6FhH8gQAEEhFYuXKlu/LKK4O/Xbt2JdqHTBCIIoAQRlFhHQQgAAEINIbA2MbUlIpCAAK1IjB58mSnP6UJEybUqm5UplgCo862UrFFUhoEIAABCEDAHwJ0jfoTCzyBQOUJ2H27JUuWVL4uqkDd6lOLoORQCYQwB6iYhAAEIACB6hBACKsTKzyFAAQgAIEcCCCEOUDFJAQgAAEIVIcAg2WqEys8LZHApk2b3EMPPRR48NJLL7nx48eP8KZXHtuufWXj2LFjgU3NgTt06FBgc/r06W758uXDoyFHFBJasWXLFrd9+3a3d+/ewJbsahTl7Nmz3aJFi0I53/6qMjdu3Oh27NgRlKn8U6ZMcbfddpubMWPG2xn/9ZvK0H0y81mr9+3bF9jYv3+/GxoacpdffnmQZ8TOHSueeOKJEWXIH5VhDMRBPskX1cFGhXaYChbNNy0cOHBgRJZ+eYftjjAaWnHnnXe63//93w/WRNUtlDWo34IFC4JVO3fu7Fqv8H58L4YA0yeK4UwpEGgjoBP/0qVLAwGzDRIFCZT+up0sJRY6qeoznLS/REp/Eu1OGyrTTsa2n2zoT2UuXLgwEGHbFvWpQTDKa+n48ePB1AW7MJAPlmydLXdOcZCIi0E4aR/5IzHS3/r1692sWbPCWfr6noa3/DTfu9XnhhtucF/+8peDGEp0oy4kzFltV5o6dSoiaFA8+qRr1KNg4EozCOjkKkGaNm1aIFZqzWzdujVoXRkBO3Hasn1q37lz5w6L4OrVq4PWpdmQcNgJOSw8EkcTQbXiVJ72UctUNpQ2b94ciI+V1fmp/U0E582bF+ynVpuESnb0J9tK+rR19ikRCKejR48GgiMBDvuj/FYHCWWn4IdtJPmelnea+ixevDhwwVq0Uf6ofG1XEjeSfwQQQv9igkcNIKAuP3WnWdefRGLDhg1Bi0HVf/bZZyMpqKWnE6taLBIPnVit9SIbOonLrrbZehmylpfyh8tRHq1TfiXr/g0WQv/ZyVz+Kq/EU/uZ/6Gsib+qC1eip67gsEjKJ9XB/I+7KEhcUCtjv7x7lSEGlp566in72vaplq2lcH5bx2f5BBDC8mOABw0kEHdC1D1CJQlPZ1Krzk6qaomExaMzb3ibuiDVqpKwWOsvKr9aYSpX5cQlCVTYdly+JOtN6OLyGgvdAx009cM7SZmqg9l++umnI3exmFm+yEysLJUA9whLxU/hTSVgXX+d9Z84ceLwKmv52QrrXtNy3GAYyxv+3LZt2/CiulXjknVBSniixE5dmIO0AOPKjVtvZXUT5rh9O9f3w7vTRtyyYiGxU7z0GRY8+W5cw+vjbLG+HAIIYTncKbXhBOwknwbDnj17guxp99XoTqVerb0g07/ms+/hz+uuuy68mPl3iYZE+ODBg4Ht3bt3Z1ZGWmZpCpZtCa0uVDqFUCN0lXRhEXVxkaYc8uZHACHMjy2WIZApARO0tCd1DUpRWrZsWaqWZKfzGtyTdZL46b5kuLWbdRlF2FNrT3WwFqDFKDy4qAg/KKM/Atwj7I8be0GgcAJ2T8262pI6YKNHrUWZdL/OfFZ+5/p+l9V6UletBES21fWqUa+a9qFBNFquStIgJeNjg2bC4k63qN+RRAj9jg/eQWCYgLXI0gqhJswrWYty2GCJX9Ry0kR9JQmejR6VoKg1ZaJSooupi7apFNYKNEFEBFOjLHwHhLBw5BRYVwKDtrh6cZFIWEozpUBPjVGSgGoEaRFJE+27JWstSfA0fSIq9bIRtU9e65L4YoInzrofa/c4bX1evmF3cAII4eAMsdAAAr1aKOETX144NCDDBlw8+uijw6MRo8pTt6OJTXi/FStW9NwvbYszqvykNsQtLpn/cduLXJ+kPjpGTPQ0SEZ1Y5BMkVHqvyyEsH927NkgAjb4QVWWyISTTpLz588PTnzh9Xl81z00nXB1kp05c6ZTy1DdjFqWcMg3rVe3Y/jkrf2UbD/ls+361LLu13Xul7YOxslsqjz5p6fSrFu3bthcuHWrcpXHkuqhdeafrS/jM2l9zDeb1mJzCk0YbTuffhJg1KifccErzwioVaWTok7OGuWoFpmWdaLXOn3XSa9TJLOuhsrRk10kLOZLVBnKZ/cUtV3LGoRi+9n9uc591YIJ79e5vdfynDlzgke1KZ/KCJczbty44d3ljx7DpvtpEkEJX2dSnrLFMGl9zHf5bFMptA4hNDJ+f9Ii9Ds+eOcRAT1VRSc5a5FJBHXi09NaJDLWbZm3yypT5alclWndtlovcbGRl53+dO5nfiqf9ot6NJvlSfopW/JLZckvLUsMZFuPdgsnLYcn6Wsf8dU0D9XPBxFJUx+rm+qg5IP/5hOf3QnwGqbufNgKAQhAIBUB6+pVy11CSvKfAC1C/2OEhxCAQEUIqJtXf2rdIoIVCVrLTYSwOrHCUwhAwHMC9kg1GzTjubu4968EEEIOBQhAAAJ9EgiPdtUcTQ3+UWuQ+4N9Ai1pN+4RlgSeYiEAgeoT0FSVzpGtNqiq+rVrTg1oETYn1tQUAhDImICNIpZZjZLViF0bNZpxUZjLkQAtwhzhYhoCEIAABPwnQIvQ/xjhIQQgAAEI5EgAIcwRLqYhAAEIQMB/Agih/zHCQwhAAAIQyJEAQpgjXExDAAIQgID/BBBC/2OEhxCAAAQgkCMBhDBHuJiGAAQgAAH/CSCE/scIDyEAAQhAIEcCCGGOcDENAQhAAAL+E0AI/Y8RHkIAAhCAQI4EEMIc4WIaAhCAAAT8J4AQ+h8jPIQABCAAgRwJIIQ5wsU0BCAAAQj4TwAh9D9GeAgBCEAAAjkSQAhzhItpCEAAAhDwnwBC6H+M8BACEIAABHIkgBDmCBfTEIAABCDgPwGE0P8Y4SEEIAABCORIACHMES6mIQABCEDAfwIIof8xwkMIQAACEMiRAEKYI1xMQwACEICA/wQQQv9jhIcQgAAEIJAjAYQwR7iYhgAEIAAB/wkghP7HCA8hAAEIQCBHAghhjnAxDQEIQAAC/hNACP2PER5CAAIQgECOBBDCHOFiGgIQgAAE/CeAEPofIzyEAAQgAIEcCSCEOcLFNAQgAAEI+E8AIfQ/RngIAQhAAAI5EkAIc4SLaQhAAAIQ8J8AQuh/jPAQAhCAAARyJIAQ5ggX0xCAAAQg4D8BhND/GOEhBCAAAQjkSAAhzBEupiEAAQhAwH8CCKH/McJDCEAAAhDIkcD/ByYI7q3DR1xhAAAAAElFTkSuQmCC\" width=\"225\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(1, 1, figsize=(2.25, 2.25))\n",
    "\n",
    "props = dict(boxstyle='round', facecolor='white', alpha=0.75)\n",
    "ax.plot([0, 1.0], [0, 1.0], 'k--')\n",
    "ax.set_ylabel(r'MSE')\n",
    "ax.set_xlabel(r'uncertainty')\n",
    "ax.set_xlim([-0.1, 1.1])\n",
    "ax.set_ylim([-0.1, 1.1])\n",
    "ax.set_xticks([0, 0.5, 1.0])\n",
    "ax.set_yticks([0, 0.5, 1.0])\n",
    "ax.set_aspect(1)\n",
    "ax.set_title(r'Calibration Diagram')\n",
    "textstr0 = 'UCE\\,=\\,0.0'\n",
    "ax.text(0.925, 0.075, textstr0, transform=ax.transAxes, fontsize=12,\n",
    "                verticalalignment='bottom',\n",
    "                horizontalalignment='right',\n",
    "                bbox=props\n",
    "                )\n",
    "\n",
    "#ax.annotate(r'BreastPathQ/EfficientNet-B4', xy=(0, 0.5), xytext=(-ax.yaxis.labelpad - 5, 0),\n",
    "#               xycoords=ax.yaxis.label, textcoords='offset points',\n",
    "#               size='large', ha='right', va='center', fontsize=10, rotation=90)\n",
    "fig.tight_layout()\n",
    "fig.show()\n",
    "\n",
    "fig.savefig(f\"empty_calib_diagram.pdf\", bbox_inches='tight', pad_inches=0.01)"
   ]
  },
  {
   "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
}