diff --git a/02_classification.ipynb b/02_classification.ipynb index 1c735e1..5d37750 100644 --- a/02_classification.ipynb +++ b/02_classification.ipynb @@ -12,7 +12,7 @@ }, { "cell_type": "code", - "execution_count": 2, + "execution_count": 15, "outputs": [], "source": [ "import numpy as np\n", @@ -30,22 +30,22 @@ }, { "cell_type": "code", - "execution_count": 8, + "execution_count": 16, "outputs": [], "source": [ "train_from_existed = False # 是否从现有数据训练，如果是的话，那就从dataset_file训练，否则就用data_dir里头的数据\n", "data_dir = \"data/dataset\" # 数据集，文件夹下必须包含`img`和`label`两个文件夹，放置相同文件名的图片和label\n", "dataset_file = \"data/dataset/dataset_2022-07-20_10-04.mat\"\n", "\n", - "color_dict = {(0, 0, 255): \"yangeng\", (255, 0, 0): 'beijing',(0, 255, 0): \"zibian\"} # 颜色对应的类别\n", - "# color_dict = {(0, 0, 255): \"yangeng\"}\n", + "# color_dict = {(0, 0, 255): \"yangeng\", (255, 0, 0): 'beijing',(0, 255, 0): \"zibian\"} # 颜色对应的类别\n", + "color_dict = {(0, 0, 255): \"yangeng\"}\n", "# color_dict = {(255, 0, 0): 'beijing'}\n", "# color_dict = {(0, 255, 0): \"zibian\"}\n", "label_index = {\"yangeng\": 1, \"beijing\": 0, \"zibian\":2} # 类别对应的序号\n", "show_samples = False # 是否展示样本\n", "\n", "# 定义一些训练量\n", - "threshold = 5 # 正样本周围多大范围内的还算是正样本\n", + "threshold = 2 # 正样本周围多大范围内的还算是正样本\n", "node_num = 20 # 如果使用ELM作为分类器物，有多少的节点\n", "negative_sample_num = None # None或者一个数字，对应生成的负样本数量" ], @@ -70,7 +70,7 @@ }, { "cell_type": "code", - "execution_count": 9, + "execution_count": 17, "outputs": [], "source": [ "dataset = read_labeled_img(data_dir, color_dict=color_dict, is_ps_color_space=False)\n", @@ -99,7 +99,7 @@ }, { "cell_type": "code", - "execution_count": 10, + "execution_count": 18, "outputs": [], "source": [ "if len(dataset) > 1:\n", @@ -130,7 +130,7 @@ }, { "cell_type": "code", - "execution_count": 11, + "execution_count": 19, "outputs": [], "source": [ "# 对数据进行预处理\n", @@ -147,13 +147,84 @@ }, { "cell_type": "code", - "execution_count": 12, + "execution_count": 20, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ - "1923it [00:00, 5114.42it/s] " + "\n", + " 0%| | 0/9019 [00:00 39\u001B[0m \u001B[38;5;28;01mreturn\u001B[39;00m \u001B[38;5;28;43mopen\u001B[39;49m\u001B[43m(\u001B[49m\u001B[43mfile_like\u001B[49m\u001B[43m,\u001B[49m\u001B[43m \u001B[49m\u001B[43mmode\u001B[49m\u001B[43m)\u001B[49m, \u001B[38;5;28;01mTrue\u001B[39;00m\n\u001B[0;32m 40\u001B[0m \u001B[38;5;28;01mexcept\u001B[39;00m \u001B[38;5;167;01mIOError\u001B[39;00m \u001B[38;5;28;01mas\u001B[39;00m e:\n\u001B[0;32m 41\u001B[0m \u001B[38;5;66;03m# Probably \"not found\"\u001B[39;00m\n", + "\u001B[1;31mFileNotFoundError\u001B[0m: [Errno 2] No such file or directory: 'dataset_2022-07-21_17-19.mat'", + "\nDuring handling of the above exception, another exception occurred:\n", + "\u001B[1;31mFileNotFoundError\u001B[0m Traceback (most recent call last)", + "Input \u001B[1;32mIn [22]\u001B[0m, in \u001B[0;36m\u001B[1;34m()\u001B[0m\n\u001B[1;32m----> 1\u001B[0m data \u001B[38;5;241m=\u001B[39m \u001B[43mscipy\u001B[49m\u001B[38;5;241;43m.\u001B[39;49m\u001B[43mio\u001B[49m\u001B[38;5;241;43m.\u001B[39;49m\u001B[43mloadmat\u001B[49m\u001B[43m(\u001B[49m\u001B[43mdataset_path\u001B[49m\u001B[43m)\u001B[49m\n\u001B[0;32m 2\u001B[0m ground_truth \u001B[38;5;241m=\u001B[39m data[\u001B[38;5;124m\"\u001B[39m\u001B[38;5;124mx\u001B[39m\u001B[38;5;124m\"\u001B[39m][data[\u001B[38;5;124m\"\u001B[39m\u001B[38;5;124my\u001B[39m\u001B[38;5;124m\"\u001B[39m]\u001B[38;5;241m.\u001B[39mravel()\u001B[38;5;241m==\u001B[39m\u001B[38;5;241m1\u001B[39m]\n", + "File \u001B[1;32m~\\miniconda3\\envs\\deepo\\lib\\site-packages\\scipy\\io\\matlab\\mio.py:224\u001B[0m, in \u001B[0;36mloadmat\u001B[1;34m(file_name, mdict, appendmat, **kwargs)\u001B[0m\n\u001B[0;32m 87\u001B[0m \u001B[38;5;124;03m\"\"\"\u001B[39;00m\n\u001B[0;32m 88\u001B[0m \u001B[38;5;124;03mLoad MATLAB file.\u001B[39;00m\n\u001B[0;32m 89\u001B[0m \n\u001B[1;32m (...)\u001B[0m\n\u001B[0;32m 221\u001B[0m \u001B[38;5;124;03m 3.14159265+3.14159265j])\u001B[39;00m\n\u001B[0;32m 222\u001B[0m \u001B[38;5;124;03m\"\"\"\u001B[39;00m\n\u001B[0;32m 223\u001B[0m variable_names \u001B[38;5;241m=\u001B[39m kwargs\u001B[38;5;241m.\u001B[39mpop(\u001B[38;5;124m'\u001B[39m\u001B[38;5;124mvariable_names\u001B[39m\u001B[38;5;124m'\u001B[39m, \u001B[38;5;28;01mNone\u001B[39;00m)\n\u001B[1;32m--> 224\u001B[0m \u001B[38;5;28;01mwith\u001B[39;00m _open_file_context(file_name, appendmat) \u001B[38;5;28;01mas\u001B[39;00m f:\n\u001B[0;32m 225\u001B[0m MR, _ \u001B[38;5;241m=\u001B[39m mat_reader_factory(f, \u001B[38;5;241m*\u001B[39m\u001B[38;5;241m*\u001B[39mkwargs)\n\u001B[0;32m 226\u001B[0m matfile_dict \u001B[38;5;241m=\u001B[39m MR\u001B[38;5;241m.\u001B[39mget_variables(variable_names)\n", + "File \u001B[1;32m~\\miniconda3\\envs\\deepo\\lib\\contextlib.py:135\u001B[0m, in \u001B[0;36m_GeneratorContextManager.__enter__\u001B[1;34m(self)\u001B[0m\n\u001B[0;32m 133\u001B[0m \u001B[38;5;28;01mdel\u001B[39;00m \u001B[38;5;28mself\u001B[39m\u001B[38;5;241m.\u001B[39margs, \u001B[38;5;28mself\u001B[39m\u001B[38;5;241m.\u001B[39mkwds, \u001B[38;5;28mself\u001B[39m\u001B[38;5;241m.\u001B[39mfunc\n\u001B[0;32m 134\u001B[0m \u001B[38;5;28;01mtry\u001B[39;00m:\n\u001B[1;32m--> 135\u001B[0m \u001B[38;5;28;01mreturn\u001B[39;00m \u001B[38;5;28;43mnext\u001B[39;49m\u001B[43m(\u001B[49m\u001B[38;5;28;43mself\u001B[39;49m\u001B[38;5;241;43m.\u001B[39;49m\u001B[43mgen\u001B[49m\u001B[43m)\u001B[49m\n\u001B[0;32m 136\u001B[0m \u001B[38;5;28;01mexcept\u001B[39;00m \u001B[38;5;167;01mStopIteration\u001B[39;00m:\n\u001B[0;32m 137\u001B[0m \u001B[38;5;28;01mraise\u001B[39;00m \u001B[38;5;167;01mRuntimeError\u001B[39;00m(\u001B[38;5;124m\"\u001B[39m\u001B[38;5;124mgenerator didn\u001B[39m\u001B[38;5;124m'\u001B[39m\u001B[38;5;124mt yield\u001B[39m\u001B[38;5;124m\"\u001B[39m) \u001B[38;5;28;01mfrom\u001B[39;00m \u001B[38;5;28mNone\u001B[39m\n", + "File \u001B[1;32m~\\miniconda3\\envs\\deepo\\lib\\site-packages\\scipy\\io\\matlab\\mio.py:17\u001B[0m, in \u001B[0;36m_open_file_context\u001B[1;34m(file_like, appendmat, mode)\u001B[0m\n\u001B[0;32m 15\u001B[0m \u001B[38;5;129m@contextmanager\u001B[39m\n\u001B[0;32m 16\u001B[0m \u001B[38;5;28;01mdef\u001B[39;00m \u001B[38;5;21m_open_file_context\u001B[39m(file_like, appendmat, mode\u001B[38;5;241m=\u001B[39m\u001B[38;5;124m'\u001B[39m\u001B[38;5;124mrb\u001B[39m\u001B[38;5;124m'\u001B[39m):\n\u001B[1;32m---> 17\u001B[0m f, opened \u001B[38;5;241m=\u001B[39m \u001B[43m_open_file\u001B[49m\u001B[43m(\u001B[49m\u001B[43mfile_like\u001B[49m\u001B[43m,\u001B[49m\u001B[43m \u001B[49m\u001B[43mappendmat\u001B[49m\u001B[43m,\u001B[49m\u001B[43m \u001B[49m\u001B[43mmode\u001B[49m\u001B[43m)\u001B[49m\n\u001B[0;32m 18\u001B[0m \u001B[38;5;28;01mtry\u001B[39;00m:\n\u001B[0;32m 19\u001B[0m \u001B[38;5;28;01myield\u001B[39;00m f\n", + "File \u001B[1;32m~\\miniconda3\\envs\\deepo\\lib\\site-packages\\scipy\\io\\matlab\\mio.py:45\u001B[0m, in \u001B[0;36m_open_file\u001B[1;34m(file_like, appendmat, mode)\u001B[0m\n\u001B[0;32m 43\u001B[0m \u001B[38;5;28;01mif\u001B[39;00m appendmat \u001B[38;5;129;01mand\u001B[39;00m \u001B[38;5;129;01mnot\u001B[39;00m file_like\u001B[38;5;241m.\u001B[39mendswith(\u001B[38;5;124m'\u001B[39m\u001B[38;5;124m.mat\u001B[39m\u001B[38;5;124m'\u001B[39m):\n\u001B[0;32m 44\u001B[0m file_like \u001B[38;5;241m+\u001B[39m\u001B[38;5;241m=\u001B[39m \u001B[38;5;124m'\u001B[39m\u001B[38;5;124m.mat\u001B[39m\u001B[38;5;124m'\u001B[39m\n\u001B[1;32m---> 45\u001B[0m \u001B[38;5;28;01mreturn\u001B[39;00m \u001B[38;5;28;43mopen\u001B[39;49m\u001B[43m(\u001B[49m\u001B[43mfile_like\u001B[49m\u001B[43m,\u001B[49m\u001B[43m \u001B[49m\u001B[43mmode\u001B[49m\u001B[43m)\u001B[49m, \u001B[38;5;28;01mTrue\u001B[39;00m\n\u001B[0;32m 46\u001B[0m \u001B[38;5;28;01melse\u001B[39;00m:\n\u001B[0;32m 47\u001B[0m \u001B[38;5;28;01mraise\u001B[39;00m \u001B[38;5;167;01mIOError\u001B[39;00m(\n\u001B[0;32m 48\u001B[0m \u001B[38;5;124m'\u001B[39m\u001B[38;5;124mReader needs file name or open file-like object\u001B[39m\u001B[38;5;124m'\u001B[39m\n\u001B[0;32m 49\u001B[0m ) \u001B[38;5;28;01mfrom\u001B[39;00m \u001B[38;5;21;01me\u001B[39;00m\n", + "\u001B[1;31mFileNotFoundError\u001B[0m: [Errno 2] No such file or directory: 'dataset_2022-07-21_17-19.mat'" + ] + } + ], "source": [ "data = scipy.io.loadmat(dataset_path)\n", "ground_truth = data[\"x\"][data[\"y\"].ravel()==1]" @@ -241,25 +332,8 @@ }, { "cell_type": "code", - "execution_count": 8, - "outputs": [ - { - "data": { - "text/plain": "", - "application/javascript": "/* Put everything inside the global mpl namespace */\n/* global mpl */\nwindow.mpl = {};\n\nmpl.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(\n '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\nmpl.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 = document.createElement('div');\n this.root.setAttribute('style', 'display: inline-block');\n this._root_extra_style(this.root);\n\n parent_element.appendChild(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 (fig.ratio !== 1) {\n fig.send_message('set_device_pixel_ratio', {\n device_pixel_ratio: fig.ratio,\n });\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\nmpl.figure.prototype._init_header = function () {\n var titlebar = document.createElement('div');\n titlebar.classList =\n 'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n var titletext = document.createElement('div');\n titletext.classList = 'ui-dialog-title';\n titletext.setAttribute(\n 'style',\n 'width: 100%; text-align: center; padding: 3px;'\n );\n titlebar.appendChild(titletext);\n this.root.appendChild(titlebar);\n this.header = titletext;\n};\n\nmpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n\nmpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n\nmpl.figure.prototype._init_canvas = function () {\n var fig = this;\n\n var canvas_div = (this.canvas_div = document.createElement('div'));\n canvas_div.setAttribute(\n 'style',\n 'border: 1px solid #ddd;' +\n 'box-sizing: content-box;' +\n 'clear: both;' +\n 'min-height: 1px;' +\n 'min-width: 1px;' +\n 'outline: 0;' +\n 'overflow: hidden;' +\n 'position: relative;' +\n 'resize: both;'\n );\n\n function on_keyboard_event_closure(name) {\n return function (event) {\n return fig.key_event(event, name);\n };\n }\n\n canvas_div.addEventListener(\n 'keydown',\n on_keyboard_event_closure('key_press')\n );\n canvas_div.addEventListener(\n 'keyup',\n on_keyboard_event_closure('key_release')\n );\n\n this._canvas_extra_style(canvas_div);\n this.root.appendChild(canvas_div);\n\n var canvas = (this.canvas = document.createElement('canvas'));\n canvas.classList.add('mpl-canvas');\n canvas.setAttribute('style', 'box-sizing: content-box;');\n\n this.context = canvas.getContext('2d');\n\n var backingStore =\n this.context.backingStorePixelRatio ||\n this.context.webkitBackingStorePixelRatio ||\n this.context.mozBackingStorePixelRatio ||\n this.context.msBackingStorePixelRatio ||\n this.context.oBackingStorePixelRatio ||\n this.context.backingStorePixelRatio ||\n 1;\n\n this.ratio = (window.devicePixelRatio || 1) / backingStore;\n\n var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n 'canvas'\n ));\n rubberband_canvas.setAttribute(\n 'style',\n 'box-sizing: content-box; position: absolute; left: 0; top: 0; z-index: 1;'\n );\n\n // Apply a ponyfill if ResizeObserver is not implemented by browser.\n if (this.ResizeObserver === undefined) {\n if (window.ResizeObserver !== undefined) {\n this.ResizeObserver = window.ResizeObserver;\n } else {\n var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n this.ResizeObserver = obs.ResizeObserver;\n }\n }\n\n this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n var nentries = entries.length;\n for (var i = 0; i < nentries; i++) {\n var entry = entries[i];\n var width, height;\n if (entry.contentBoxSize) {\n if (entry.contentBoxSize instanceof Array) {\n // Chrome 84 implements new version of spec.\n width = entry.contentBoxSize[0].inlineSize;\n height = entry.contentBoxSize[0].blockSize;\n } else {\n // Firefox implements old version of spec.\n width = entry.contentBoxSize.inlineSize;\n height = entry.contentBoxSize.blockSize;\n }\n } else {\n // Chrome <84 implements even older version of spec.\n width = entry.contentRect.width;\n height = entry.contentRect.height;\n }\n\n // Keep the size of the canvas and rubber band canvas in sync with\n // the canvas container.\n if (entry.devicePixelContentBoxSize) {\n // Chrome 84 implements new version of spec.\n canvas.setAttribute(\n 'width',\n entry.devicePixelContentBoxSize[0].inlineSize\n );\n canvas.setAttribute(\n 'height',\n entry.devicePixelContentBoxSize[0].blockSize\n );\n } else {\n canvas.setAttribute('width', width * fig.ratio);\n canvas.setAttribute('height', height * fig.ratio);\n }\n canvas.setAttribute(\n 'style',\n 'width: ' + width + 'px; height: ' + height + 'px;'\n );\n\n rubberband_canvas.setAttribute('width', width);\n rubberband_canvas.setAttribute('height', height);\n\n // And update the size in Python. We ignore the initial 0/0 size\n // that occurs as the element is placed into the DOM, which should\n // otherwise not happen due to the minimum size styling.\n if (fig.ws.readyState == 1 && width != 0 && height != 0) {\n fig.request_resize(width, height);\n }\n }\n });\n this.resizeObserverInstance.observe(canvas_div);\n\n function on_mouse_event_closure(name) {\n return function (event) {\n return fig.mouse_event(event, name);\n };\n }\n\n rubberband_canvas.addEventListener(\n 'mousedown',\n on_mouse_event_closure('button_press')\n );\n rubberband_canvas.addEventListener(\n 'mouseup',\n on_mouse_event_closure('button_release')\n );\n rubberband_canvas.addEventListener(\n 'dblclick',\n on_mouse_event_closure('dblclick')\n );\n // Throttle sequential mouse events to 1 every 20ms.\n rubberband_canvas.addEventListener(\n 'mousemove',\n on_mouse_event_closure('motion_notify')\n );\n\n rubberband_canvas.addEventListener(\n 'mouseenter',\n on_mouse_event_closure('figure_enter')\n );\n rubberband_canvas.addEventListener(\n 'mouseleave',\n on_mouse_event_closure('figure_leave')\n );\n\n canvas_div.addEventListener('wheel', function (event) {\n if (event.deltaY < 0) {\n event.step = 1;\n } else {\n event.step = -1;\n }\n on_mouse_event_closure('scroll')(event);\n });\n\n canvas_div.appendChild(canvas);\n canvas_div.appendChild(rubberband_canvas);\n\n this.rubberband_context = rubberband_canvas.getContext('2d');\n this.rubberband_context.strokeStyle = '#000000';\n\n this._resize_canvas = function (width, height, forward) {\n if (forward) {\n canvas_div.style.width = width + 'px';\n canvas_div.style.height = height + 'px';\n }\n };\n\n // Disable right mouse context menu.\n this.rubberband_canvas.addEventListener('contextmenu', function (_e) {\n event.preventDefault();\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\nmpl.figure.prototype._init_toolbar = function () {\n var fig = this;\n\n var toolbar = document.createElement('div');\n toolbar.classList = 'mpl-toolbar';\n this.root.appendChild(toolbar);\n\n function on_click_closure(name) {\n return function (_event) {\n return fig.toolbar_button_onclick(name);\n };\n }\n\n function on_mouseover_closure(tooltip) {\n return function (event) {\n if (!event.currentTarget.disabled) {\n return fig.toolbar_button_onmouseover(tooltip);\n }\n };\n }\n\n fig.buttons = {};\n var buttonGroup = document.createElement('div');\n buttonGroup.classList = 'mpl-button-group';\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 /* Instead of a spacer, we start a new button group. */\n if (buttonGroup.hasChildNodes()) {\n toolbar.appendChild(buttonGroup);\n }\n buttonGroup = document.createElement('div');\n buttonGroup.classList = 'mpl-button-group';\n continue;\n }\n\n var button = (fig.buttons[name] = document.createElement('button'));\n button.classList = 'mpl-widget';\n button.setAttribute('role', 'button');\n button.setAttribute('aria-disabled', 'false');\n button.addEventListener('click', on_click_closure(method_name));\n button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n\n var icon_img = document.createElement('img');\n icon_img.src = '_images/' + image + '.png';\n icon_img.srcset = '_images/' + image + '_large.png 2x';\n icon_img.alt = tooltip;\n button.appendChild(icon_img);\n\n buttonGroup.appendChild(button);\n }\n\n if (buttonGroup.hasChildNodes()) {\n toolbar.appendChild(buttonGroup);\n }\n\n var fmt_picker = document.createElement('select');\n fmt_picker.classList = 'mpl-widget';\n toolbar.appendChild(fmt_picker);\n this.format_dropdown = fmt_picker;\n\n for (var ind in mpl.extensions) {\n var fmt = mpl.extensions[ind];\n var option = document.createElement('option');\n option.selected = fmt === mpl.default_extension;\n option.innerHTML = fmt;\n fmt_picker.appendChild(option);\n }\n\n var status_bar = document.createElement('span');\n status_bar.classList = 'mpl-message';\n toolbar.appendChild(status_bar);\n this.message = status_bar;\n};\n\nmpl.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\nmpl.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\nmpl.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\nmpl.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\nmpl.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], msg['forward']);\n fig.send_message('refresh', {});\n }\n};\n\nmpl.figure.prototype.handle_rubberband = function (fig, msg) {\n var x0 = msg['x0'] / fig.ratio;\n var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n var x1 = msg['x1'] / fig.ratio;\n var y1 = (fig.canvas.height - msg['y1']) / fig.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,\n 0,\n fig.canvas.width / fig.ratio,\n fig.canvas.height / fig.ratio\n );\n\n fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n};\n\nmpl.figure.prototype.handle_figure_label = function (fig, msg) {\n // Updates the figure title.\n fig.header.textContent = msg['label'];\n};\n\nmpl.figure.prototype.handle_cursor = function (fig, msg) {\n fig.rubberband_canvas.style.cursor = msg['cursor'];\n};\n\nmpl.figure.prototype.handle_message = function (fig, msg) {\n fig.message.textContent = msg['message'];\n};\n\nmpl.figure.prototype.handle_draw = function (fig, _msg) {\n // Request the server to send over a new figure.\n fig.send_draw_message();\n};\n\nmpl.figure.prototype.handle_image_mode = function (fig, msg) {\n fig.image_mode = msg['mode'];\n};\n\nmpl.figure.prototype.handle_history_buttons = function (fig, msg) {\n for (var key in msg) {\n if (!(key in fig.buttons)) {\n continue;\n }\n fig.buttons[key].disabled = !msg[key];\n fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n }\n};\n\nmpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n if (msg['mode'] === 'PAN') {\n fig.buttons['Pan'].classList.add('active');\n fig.buttons['Zoom'].classList.remove('active');\n } else if (msg['mode'] === 'ZOOM') {\n fig.buttons['Pan'].classList.remove('active');\n fig.buttons['Zoom'].classList.add('active');\n } else {\n fig.buttons['Pan'].classList.remove('active');\n fig.buttons['Zoom'].classList.remove('active');\n }\n};\n\nmpl.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.\nmpl.figure.prototype._make_on_message_function = function (fig) {\n return function socket_on_message(evt) {\n if (evt.data instanceof Blob) {\n var img = evt.data;\n if (img.type !== 'image/png') {\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 img.type = 'image/png';\n }\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\n fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n img\n );\n fig.updated_canvas_event();\n fig.waiting = false;\n return;\n } else if (\n typeof evt.data === 'string' &&\n evt.data.slice(0, 21) === 'data:image/png;base64'\n ) {\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(\n \"No handler for the '\" + msg_type + \"' message type: \",\n msg\n );\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(\n \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n e,\n e.stack,\n msg\n );\n }\n }\n };\n};\n\n// from https://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\nmpl.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 }\n if (e.target) {\n targ = e.target;\n } else if (e.srcElement) {\n targ = e.srcElement;\n }\n if (targ.nodeType === 3) {\n // defeat Safari bug\n targ = targ.parentNode;\n }\n\n // pageX,Y are the mouse positions relative to the document\n var boundingRect = targ.getBoundingClientRect();\n var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n var y = e.pageY - (boundingRect.top + document.body.scrollTop);\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 * https://stackoverflow.com/a/24161582/3208463\n */\nfunction simpleKeys(original) {\n return Object.keys(original).reduce(function (obj, key) {\n if (typeof original[key] !== 'object') {\n obj[key] = original[key];\n }\n return obj;\n }, {});\n}\n\nmpl.figure.prototype.mouse_event = function (event, name) {\n var canvas_pos = mpl.findpos(event);\n\n if (name === 'button_press') {\n this.canvas.focus();\n this.canvas_div.focus();\n }\n\n var x = canvas_pos.x * this.ratio;\n var y = canvas_pos.y * this.ratio;\n\n this.send_message(name, {\n x: x,\n y: y,\n button: event.button,\n step: event.step,\n guiEvent: simpleKeys(event),\n });\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\nmpl.figure.prototype._key_event_extra = function (_event, _name) {\n // Handle any extra behaviour associated with a key event\n};\n\nmpl.figure.prototype.key_event = function (event, name) {\n // Prevent repeat events\n if (name === 'key_press') {\n if (event.key === this._key) {\n return;\n } else {\n this._key = event.key;\n }\n }\n if (name === 'key_release') {\n this._key = null;\n }\n\n var value = '';\n if (event.ctrlKey && event.key !== 'Control') {\n value += 'ctrl+';\n }\n else if (event.altKey && event.key !== 'Alt') {\n value += 'alt+';\n }\n else if (event.shiftKey && event.key !== 'Shift') {\n value += 'shift+';\n }\n\n value += 'k' + event.key;\n\n this._key_event_extra(event, name);\n\n this.send_message(name, { key: value, guiEvent: simpleKeys(event) });\n return false;\n};\n\nmpl.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\nmpl.figure.prototype.toolbar_button_onmouseover = function (tooltip) {\n this.message.textContent = tooltip;\n};\n\n///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n// prettier-ignore\nvar _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\nmpl.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\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n\nmpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n\nmpl.default_extension = \"png\";/* global mpl */\n\nvar 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.binaryType = comm.kernel.ws.binaryType;\n ws.readyState = comm.kernel.ws.readyState;\n function updateReadyState(_event) {\n if (comm.kernel.ws) {\n ws.readyState = comm.kernel.ws.readyState;\n } else {\n ws.readyState = 3; // Closed state.\n }\n }\n comm.kernel.ws.addEventListener('open', updateReadyState);\n comm.kernel.ws.addEventListener('close', updateReadyState);\n comm.kernel.ws.addEventListener('error', updateReadyState);\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 var data = msg['content']['data'];\n if (data['blob'] !== undefined) {\n data = {\n data: new Blob(msg['buffers'], { type: data['blob'] }),\n };\n }\n // Pass the mpl event to the overridden (by mpl) onmessage function.\n ws.onmessage(data);\n });\n return ws;\n};\n\nmpl.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 = document.getElementById(id);\n var ws_proxy = comm_websocket_adapter(comm);\n\n function ondownload(figure, _format) {\n window.open(figure.canvas.toDataURL());\n }\n\n var fig = new mpl.figure(id, ws_proxy, ondownload, element);\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;\n fig.cell_info = mpl.find_output_cell(\"
\");\n if (!fig.cell_info) {\n console.error('Failed to find cell for figure', id, fig);\n return;\n }\n fig.cell_info[0].output_area.element.on(\n 'cleared',\n { fig: fig },\n fig._remove_fig_handler\n );\n};\n\nmpl.figure.prototype.handle_close = function (fig, msg) {\n var width = fig.canvas.width / fig.ratio;\n fig.cell_info[0].output_area.element.off(\n 'cleared',\n fig._remove_fig_handler\n );\n fig.resizeObserverInstance.unobserve(fig.canvas_div);\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.innerHTML =\n '';\n fig.close_ws(fig, msg);\n};\n\nmpl.figure.prototype.close_ws = function (fig, msg) {\n fig.send_message('closing', msg);\n // fig.ws.close()\n};\n\nmpl.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 / this.ratio;\n var dataURL = this.canvas.toDataURL();\n this.cell_info[1]['text/html'] =\n '';\n};\n\nmpl.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 () {\n fig.push_to_output();\n }, 1000);\n};\n\nmpl.figure.prototype._init_toolbar = function () {\n var fig = this;\n\n var toolbar = document.createElement('div');\n toolbar.classList = 'btn-toolbar';\n this.root.appendChild(toolbar);\n\n function on_click_closure(name) {\n return function (_event) {\n return fig.toolbar_button_onclick(name);\n };\n }\n\n function on_mouseover_closure(tooltip) {\n return function (event) {\n if (!event.currentTarget.disabled) {\n return fig.toolbar_button_onmouseover(tooltip);\n }\n };\n }\n\n fig.buttons = {};\n var buttonGroup = document.createElement('div');\n buttonGroup.classList = 'btn-group';\n var button;\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 /* Instead of a spacer, we start a new button group. */\n if (buttonGroup.hasChildNodes()) {\n toolbar.appendChild(buttonGroup);\n }\n buttonGroup = document.createElement('div');\n buttonGroup.classList = 'btn-group';\n continue;\n }\n\n button = fig.buttons[name] = document.createElement('button');\n button.classList = 'btn btn-default';\n button.href = '#';\n button.title = name;\n button.innerHTML = '';\n button.addEventListener('click', on_click_closure(method_name));\n button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n buttonGroup.appendChild(button);\n }\n\n if (buttonGroup.hasChildNodes()) {\n toolbar.appendChild(buttonGroup);\n }\n\n // Add the status bar.\n var status_bar = document.createElement('span');\n status_bar.classList = 'mpl-message pull-right';\n toolbar.appendChild(status_bar);\n this.message = status_bar;\n\n // Add the close button to the window.\n var buttongrp = document.createElement('div');\n buttongrp.classList = 'btn-group inline pull-right';\n button = document.createElement('button');\n button.classList = 'btn btn-mini btn-primary';\n button.href = '#';\n button.title = 'Stop Interaction';\n button.innerHTML = '';\n button.addEventListener('click', function (_evt) {\n fig.handle_close(fig, {});\n });\n button.addEventListener(\n 'mouseover',\n on_mouseover_closure('Stop Interaction')\n );\n buttongrp.appendChild(button);\n var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n titlebar.insertBefore(buttongrp, titlebar.firstChild);\n};\n\nmpl.figure.prototype._remove_fig_handler = function (event) {\n var fig = event.data.fig;\n if (event.target !== this) {\n // Ignore bubbled events from children.\n return;\n }\n fig.close_ws(fig, {});\n};\n\nmpl.figure.prototype._root_extra_style = function (el) {\n el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n};\n\nmpl.figure.prototype._canvas_extra_style = function (el) {\n // this is important to make the div 'focusable\n el.setAttribute('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 } else {\n // location in version 2\n IPython.keyboard_manager.register_events(el);\n }\n};\n\nmpl.figure.prototype._key_event_extra = function (event, _name) {\n // Check for shift+enter\n if (event.shiftKey && event.which === 13) {\n this.canvas_div.blur();\n // select the cell after this one\n var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n IPython.notebook.select(index + 1);\n }\n};\n\nmpl.figure.prototype.handle_save = function (fig, _msg) {\n fig.ondownload(fig, null);\n};\n\nmpl.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.\nif (IPython.notebook.kernel !== null) {\n IPython.notebook.kernel.comm_manager.register_target(\n 'matplotlib',\n mpl.mpl_figure_comm\n );\n}\n" - }, - "metadata": {}, - "output_type": "display_data" - }, - { - "data": { - "text/plain": "", - "text/html": "
" - }, - "metadata": {}, - "output_type": "display_data" - } - ], + "execution_count": null, + "outputs": [], "source": [ "%matplotlib notebook\n", "model = AnonymousColorDetector(model_path)\n", diff --git a/04_multi_classification.ipynb b/04_multi_classification.ipynb index 3f30c7d..1ce3ea4 100644 --- a/04_multi_classification.ipynb +++ b/04_multi_classification.ipynb @@ -26,7 +26,7 @@ "import pickle\n", "from sklearn.tree import DecisionTreeClassifier\n", "# %matplotlib notebook\n", - "from main_test import virtual_main\n", + "from main_test import pony_run\n", "from models import AnonymousColorDetector\n", "from utils import lab_scatter" ], @@ -39,11 +39,12 @@ }, { "cell_type": "code", - "execution_count": 18, + "execution_count": 25, "outputs": [], "source": [ - "img_path = r\"C:\\Users\\FEIJINTI\\Desktop\\721\\zazhi\\Image_2022_0721_1351_38_946-002034.bmp\"\n", - "model_path = [\"models/beijing_dt_2022-07-21_16-44.model\", \"models/tobacco_dt_2022-07-21_16-30.model\",\"models/zibian_dt_2022-07-21_16-45.model\",\"dt_2022-07-21_17-19.model\"]" + "img_path = r\"C:\\Users\\FEIJINTI\\Desktop\\7261\\tobacco\\Image_2022_0726_1410_31_993-000254.bmp\"\n", + "# img_path = r\"E:\\zhouchao\\725data\\14.png\"\n", + "model_path = [\"models/background_dt_2022-07-26_15-56.model\", \"models/tobacco_dt_2022-07-26_15-57.model\",\"models/zibian_dt_2022-07-21_16-45.model\",\"models/dt_2022-07-21_17-19.model\"]" ], "metadata": { "collapsed": false, @@ -54,10 +55,23 @@ }, { "cell_type": "code", - "execution_count": 19, + "execution_count": 25, + "outputs": [], + "source": [], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 26, "outputs": [], "source": [ "img = cv2.imread(img_path)[:, :, ::-1]\n", + "# img = cv2.imread(img_path)\n", "img = cv2.resize(img, (1024, 256))\n", "b_detector = AnonymousColorDetector(file_path=model_path[0])\n", "b_result = b_detector.predict(img).astype(np.uint8)" @@ -71,7 +85,7 @@ }, { "cell_type": "code", - "execution_count": 20, + "execution_count": 27, "outputs": [], "source": [ "t_detector = AnonymousColorDetector(file_path=model_path[1])\n", @@ -87,7 +101,7 @@ }, { "cell_type": "code", - "execution_count": 21, + "execution_count": 28, "outputs": [], "source": [ "z_detector = AnonymousColorDetector(file_path=model_path[2])\n", @@ -102,7 +116,7 @@ }, { "cell_type": "code", - "execution_count": 22, + "execution_count": 29, "outputs": [], "source": [ "s_detector = AnonymousColorDetector(file_path=model_path[3])\n", @@ -117,7 +131,7 @@ }, { "cell_type": "code", - "execution_count": 25, + "execution_count": 30, "outputs": [], "source": [ "result = 1 - (b_result | t_result)" @@ -131,22 +145,16 @@ }, { "cell_type": "code", - "execution_count": 26, + "execution_count": 31, "outputs": [ { "data": { - "text/plain": "", - "application/javascript": "/* Put everything inside the global mpl namespace */\n/* global mpl */\nwindow.mpl = {};\n\nmpl.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(\n '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\nmpl.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 = document.createElement('div');\n this.root.setAttribute('style', 'display: inline-block');\n this._root_extra_style(this.root);\n\n parent_element.appendChild(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 (fig.ratio !== 1) {\n fig.send_message('set_device_pixel_ratio', {\n device_pixel_ratio: fig.ratio,\n });\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\nmpl.figure.prototype._init_header = function () {\n var titlebar = document.createElement('div');\n titlebar.classList =\n 'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n var titletext = document.createElement('div');\n titletext.classList = 'ui-dialog-title';\n titletext.setAttribute(\n 'style',\n 'width: 100%; text-align: center; padding: 3px;'\n );\n titlebar.appendChild(titletext);\n this.root.appendChild(titlebar);\n this.header = titletext;\n};\n\nmpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n\nmpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n\nmpl.figure.prototype._init_canvas = function () {\n var fig = this;\n\n var canvas_div = (this.canvas_div = document.createElement('div'));\n canvas_div.setAttribute(\n 'style',\n 'border: 1px solid #ddd;' +\n 'box-sizing: content-box;' +\n 'clear: both;' +\n 'min-height: 1px;' +\n 'min-width: 1px;' +\n 'outline: 0;' +\n 'overflow: hidden;' +\n 'position: relative;' +\n 'resize: both;'\n );\n\n function on_keyboard_event_closure(name) {\n return function (event) {\n return fig.key_event(event, name);\n };\n }\n\n canvas_div.addEventListener(\n 'keydown',\n on_keyboard_event_closure('key_press')\n );\n canvas_div.addEventListener(\n 'keyup',\n on_keyboard_event_closure('key_release')\n );\n\n this._canvas_extra_style(canvas_div);\n this.root.appendChild(canvas_div);\n\n var canvas = (this.canvas = document.createElement('canvas'));\n canvas.classList.add('mpl-canvas');\n canvas.setAttribute('style', 'box-sizing: content-box;');\n\n this.context = canvas.getContext('2d');\n\n var backingStore =\n this.context.backingStorePixelRatio ||\n this.context.webkitBackingStorePixelRatio ||\n this.context.mozBackingStorePixelRatio ||\n this.context.msBackingStorePixelRatio ||\n this.context.oBackingStorePixelRatio ||\n this.context.backingStorePixelRatio ||\n 1;\n\n this.ratio = (window.devicePixelRatio || 1) / backingStore;\n\n var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n 'canvas'\n ));\n rubberband_canvas.setAttribute(\n 'style',\n 'box-sizing: content-box; position: absolute; left: 0; top: 0; z-index: 1;'\n );\n\n // Apply a ponyfill if ResizeObserver is not implemented by browser.\n if (this.ResizeObserver === undefined) {\n if (window.ResizeObserver !== undefined) {\n this.ResizeObserver = window.ResizeObserver;\n } else {\n var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n this.ResizeObserver = obs.ResizeObserver;\n }\n }\n\n this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n var nentries = entries.length;\n for (var i = 0; i < nentries; i++) {\n var entry = entries[i];\n var width, height;\n if (entry.contentBoxSize) {\n if (entry.contentBoxSize instanceof Array) {\n // Chrome 84 implements new version of spec.\n width = entry.contentBoxSize[0].inlineSize;\n height = entry.contentBoxSize[0].blockSize;\n } else {\n // Firefox implements old version of spec.\n width = entry.contentBoxSize.inlineSize;\n height = entry.contentBoxSize.blockSize;\n }\n } else {\n // Chrome <84 implements even older version of spec.\n width = entry.contentRect.width;\n height = entry.contentRect.height;\n }\n\n // Keep the size of the canvas and rubber band canvas in sync with\n // the canvas container.\n if (entry.devicePixelContentBoxSize) {\n // Chrome 84 implements new version of spec.\n canvas.setAttribute(\n 'width',\n entry.devicePixelContentBoxSize[0].inlineSize\n );\n canvas.setAttribute(\n 'height',\n entry.devicePixelContentBoxSize[0].blockSize\n );\n } else {\n canvas.setAttribute('width', width * fig.ratio);\n canvas.setAttribute('height', height * fig.ratio);\n }\n canvas.setAttribute(\n 'style',\n 'width: ' + width + 'px; height: ' + height + 'px;'\n );\n\n rubberband_canvas.setAttribute('width', width);\n rubberband_canvas.setAttribute('height', height);\n\n // And update the size in Python. We ignore the initial 0/0 size\n // that occurs as the element is placed into the DOM, which should\n // otherwise not happen due to the minimum size styling.\n if (fig.ws.readyState == 1 && width != 0 && height != 0) {\n fig.request_resize(width, height);\n }\n }\n });\n this.resizeObserverInstance.observe(canvas_div);\n\n function on_mouse_event_closure(name) {\n return function (event) {\n return fig.mouse_event(event, name);\n };\n }\n\n rubberband_canvas.addEventListener(\n 'mousedown',\n on_mouse_event_closure('button_press')\n );\n rubberband_canvas.addEventListener(\n 'mouseup',\n on_mouse_event_closure('button_release')\n );\n rubberband_canvas.addEventListener(\n 'dblclick',\n on_mouse_event_closure('dblclick')\n );\n // Throttle sequential mouse events to 1 every 20ms.\n rubberband_canvas.addEventListener(\n 'mousemove',\n on_mouse_event_closure('motion_notify')\n );\n\n rubberband_canvas.addEventListener(\n 'mouseenter',\n on_mouse_event_closure('figure_enter')\n );\n rubberband_canvas.addEventListener(\n 'mouseleave',\n on_mouse_event_closure('figure_leave')\n );\n\n canvas_div.addEventListener('wheel', function (event) {\n if (event.deltaY < 0) {\n event.step = 1;\n } else {\n event.step = -1;\n }\n on_mouse_event_closure('scroll')(event);\n });\n\n canvas_div.appendChild(canvas);\n canvas_div.appendChild(rubberband_canvas);\n\n this.rubberband_context = rubberband_canvas.getContext('2d');\n this.rubberband_context.strokeStyle = '#000000';\n\n this._resize_canvas = function (width, height, forward) {\n if (forward) {\n canvas_div.style.width = width + 'px';\n canvas_div.style.height = height + 'px';\n }\n };\n\n // Disable right mouse context menu.\n this.rubberband_canvas.addEventListener('contextmenu', function (_e) {\n event.preventDefault();\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\nmpl.figure.prototype._init_toolbar = function () {\n var fig = this;\n\n var toolbar = document.createElement('div');\n toolbar.classList = 'mpl-toolbar';\n this.root.appendChild(toolbar);\n\n function on_click_closure(name) {\n return function (_event) {\n return fig.toolbar_button_onclick(name);\n };\n }\n\n function on_mouseover_closure(tooltip) {\n return function (event) {\n if (!event.currentTarget.disabled) {\n return fig.toolbar_button_onmouseover(tooltip);\n }\n };\n }\n\n fig.buttons = {};\n var buttonGroup = document.createElement('div');\n buttonGroup.classList = 'mpl-button-group';\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 /* Instead of a spacer, we start a new button group. */\n if (buttonGroup.hasChildNodes()) {\n toolbar.appendChild(buttonGroup);\n }\n buttonGroup = document.createElement('div');\n buttonGroup.classList = 'mpl-button-group';\n continue;\n }\n\n var button = (fig.buttons[name] = document.createElement('button'));\n button.classList = 'mpl-widget';\n button.setAttribute('role', 'button');\n button.setAttribute('aria-disabled', 'false');\n button.addEventListener('click', on_click_closure(method_name));\n button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n\n var icon_img = document.createElement('img');\n icon_img.src = '_images/' + image + '.png';\n icon_img.srcset = '_images/' + image + '_large.png 2x';\n icon_img.alt = tooltip;\n button.appendChild(icon_img);\n\n buttonGroup.appendChild(button);\n }\n\n if (buttonGroup.hasChildNodes()) {\n toolbar.appendChild(buttonGroup);\n }\n\n var fmt_picker = document.createElement('select');\n fmt_picker.classList = 'mpl-widget';\n toolbar.appendChild(fmt_picker);\n this.format_dropdown = fmt_picker;\n\n for (var ind in mpl.extensions) {\n var fmt = mpl.extensions[ind];\n var option = document.createElement('option');\n option.selected = fmt === mpl.default_extension;\n option.innerHTML = fmt;\n fmt_picker.appendChild(option);\n }\n\n var status_bar = document.createElement('span');\n status_bar.classList = 'mpl-message';\n toolbar.appendChild(status_bar);\n this.message = status_bar;\n};\n\nmpl.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\nmpl.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\nmpl.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\nmpl.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\nmpl.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], msg['forward']);\n fig.send_message('refresh', {});\n }\n};\n\nmpl.figure.prototype.handle_rubberband = function (fig, msg) {\n var x0 = msg['x0'] / fig.ratio;\n var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n var x1 = msg['x1'] / fig.ratio;\n var y1 = (fig.canvas.height - msg['y1']) / fig.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,\n 0,\n fig.canvas.width / fig.ratio,\n fig.canvas.height / fig.ratio\n );\n\n fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n};\n\nmpl.figure.prototype.handle_figure_label = function (fig, msg) {\n // Updates the figure title.\n fig.header.textContent = msg['label'];\n};\n\nmpl.figure.prototype.handle_cursor = function (fig, msg) {\n fig.rubberband_canvas.style.cursor = msg['cursor'];\n};\n\nmpl.figure.prototype.handle_message = function (fig, msg) {\n fig.message.textContent = msg['message'];\n};\n\nmpl.figure.prototype.handle_draw = function (fig, _msg) {\n // Request the server to send over a new figure.\n fig.send_draw_message();\n};\n\nmpl.figure.prototype.handle_image_mode = function (fig, msg) {\n fig.image_mode = msg['mode'];\n};\n\nmpl.figure.prototype.handle_history_buttons = function (fig, msg) {\n for (var key in msg) {\n if (!(key in fig.buttons)) {\n continue;\n }\n fig.buttons[key].disabled = !msg[key];\n fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n }\n};\n\nmpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n if (msg['mode'] === 'PAN') {\n fig.buttons['Pan'].classList.add('active');\n fig.buttons['Zoom'].classList.remove('active');\n } else if (msg['mode'] === 'ZOOM') {\n fig.buttons['Pan'].classList.remove('active');\n fig.buttons['Zoom'].classList.add('active');\n } else {\n fig.buttons['Pan'].classList.remove('active');\n fig.buttons['Zoom'].classList.remove('active');\n }\n};\n\nmpl.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.\nmpl.figure.prototype._make_on_message_function = function (fig) {\n return function socket_on_message(evt) {\n if (evt.data instanceof Blob) {\n var img = evt.data;\n if (img.type !== 'image/png') {\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 img.type = 'image/png';\n }\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\n fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n img\n );\n fig.updated_canvas_event();\n fig.waiting = false;\n return;\n } else if (\n typeof evt.data === 'string' &&\n evt.data.slice(0, 21) === 'data:image/png;base64'\n ) {\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(\n \"No handler for the '\" + msg_type + \"' message type: \",\n msg\n );\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(\n \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n e,\n e.stack,\n msg\n );\n }\n }\n };\n};\n\n// from https://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\nmpl.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 }\n if (e.target) {\n targ = e.target;\n } else if (e.srcElement) {\n targ = e.srcElement;\n }\n if (targ.nodeType === 3) {\n // defeat Safari bug\n targ = targ.parentNode;\n }\n\n // pageX,Y are the mouse positions relative to the document\n var boundingRect = targ.getBoundingClientRect();\n var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n var y = e.pageY - (boundingRect.top + document.body.scrollTop);\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 * https://stackoverflow.com/a/24161582/3208463\n */\nfunction simpleKeys(original) {\n return Object.keys(original).reduce(function (obj, key) {\n if (typeof original[key] !== 'object') {\n obj[key] = original[key];\n }\n return obj;\n }, {});\n}\n\nmpl.figure.prototype.mouse_event = function (event, name) {\n var canvas_pos = mpl.findpos(event);\n\n if (name === 'button_press') {\n this.canvas.focus();\n this.canvas_div.focus();\n }\n\n var x = canvas_pos.x * this.ratio;\n var y = canvas_pos.y * this.ratio;\n\n this.send_message(name, {\n x: x,\n y: y,\n button: event.button,\n step: event.step,\n guiEvent: simpleKeys(event),\n });\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\nmpl.figure.prototype._key_event_extra = function (_event, _name) {\n // Handle any extra behaviour associated with a key event\n};\n\nmpl.figure.prototype.key_event = function (event, name) {\n // Prevent repeat events\n if (name === 'key_press') {\n if (event.key === this._key) {\n return;\n } else {\n this._key = event.key;\n }\n }\n if (name === 'key_release') {\n this._key = null;\n }\n\n var value = '';\n if (event.ctrlKey && event.key !== 'Control') {\n value += 'ctrl+';\n }\n else if (event.altKey && event.key !== 'Alt') {\n value += 'alt+';\n }\n else if (event.shiftKey && event.key !== 'Shift') {\n value += 'shift+';\n }\n\n value += 'k' + event.key;\n\n this._key_event_extra(event, name);\n\n this.send_message(name, { key: value, guiEvent: simpleKeys(event) });\n return false;\n};\n\nmpl.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\nmpl.figure.prototype.toolbar_button_onmouseover = function (tooltip) {\n this.message.textContent = tooltip;\n};\n\n///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n// prettier-ignore\nvar _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\nmpl.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\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n\nmpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n\nmpl.default_extension = \"png\";/* global mpl */\n\nvar 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.binaryType = comm.kernel.ws.binaryType;\n ws.readyState = comm.kernel.ws.readyState;\n function updateReadyState(_event) {\n if (comm.kernel.ws) {\n ws.readyState = comm.kernel.ws.readyState;\n } else {\n ws.readyState = 3; // Closed state.\n }\n }\n comm.kernel.ws.addEventListener('open', updateReadyState);\n comm.kernel.ws.addEventListener('close', updateReadyState);\n comm.kernel.ws.addEventListener('error', updateReadyState);\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 var data = msg['content']['data'];\n if (data['blob'] !== undefined) {\n data = {\n data: new Blob(msg['buffers'], { type: data['blob'] }),\n };\n }\n // Pass the mpl event to the overridden (by mpl) onmessage function.\n ws.onmessage(data);\n });\n return ws;\n};\n\nmpl.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 = document.getElementById(id);\n var ws_proxy = comm_websocket_adapter(comm);\n\n function ondownload(figure, _format) {\n window.open(figure.canvas.toDataURL());\n }\n\n var fig = new mpl.figure(id, ws_proxy, ondownload, element);\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;\n fig.cell_info = mpl.find_output_cell(\"
\");\n if (!fig.cell_info) {\n console.error('Failed to find cell for figure', id, fig);\n return;\n }\n fig.cell_info[0].output_area.element.on(\n 'cleared',\n { fig: fig },\n fig._remove_fig_handler\n );\n};\n\nmpl.figure.prototype.handle_close = function (fig, msg) {\n var width = fig.canvas.width / fig.ratio;\n fig.cell_info[0].output_area.element.off(\n 'cleared',\n fig._remove_fig_handler\n );\n fig.resizeObserverInstance.unobserve(fig.canvas_div);\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.innerHTML =\n '';\n fig.close_ws(fig, msg);\n};\n\nmpl.figure.prototype.close_ws = function (fig, msg) {\n fig.send_message('closing', msg);\n // fig.ws.close()\n};\n\nmpl.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 / this.ratio;\n var dataURL = this.canvas.toDataURL();\n this.cell_info[1]['text/html'] =\n '';\n};\n\nmpl.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 () {\n fig.push_to_output();\n }, 1000);\n};\n\nmpl.figure.prototype._init_toolbar = function () {\n var fig = this;\n\n var toolbar = document.createElement('div');\n toolbar.classList = 'btn-toolbar';\n this.root.appendChild(toolbar);\n\n function on_click_closure(name) {\n return function (_event) {\n return fig.toolbar_button_onclick(name);\n };\n }\n\n function on_mouseover_closure(tooltip) {\n return function (event) {\n if (!event.currentTarget.disabled) {\n return fig.toolbar_button_onmouseover(tooltip);\n }\n };\n }\n\n fig.buttons = {};\n var buttonGroup = document.createElement('div');\n buttonGroup.classList = 'btn-group';\n var button;\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 /* Instead of a spacer, we start a new button group. */\n if (buttonGroup.hasChildNodes()) {\n toolbar.appendChild(buttonGroup);\n }\n buttonGroup = document.createElement('div');\n buttonGroup.classList = 'btn-group';\n continue;\n }\n\n button = fig.buttons[name] = document.createElement('button');\n button.classList = 'btn btn-default';\n button.href = '#';\n button.title = name;\n button.innerHTML = '';\n button.addEventListener('click', on_click_closure(method_name));\n button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n buttonGroup.appendChild(button);\n }\n\n if (buttonGroup.hasChildNodes()) {\n toolbar.appendChild(buttonGroup);\n }\n\n // Add the status bar.\n var status_bar = document.createElement('span');\n status_bar.classList = 'mpl-message pull-right';\n toolbar.appendChild(status_bar);\n this.message = status_bar;\n\n // Add the close button to the window.\n var buttongrp = document.createElement('div');\n buttongrp.classList = 'btn-group inline pull-right';\n button = document.createElement('button');\n button.classList = 'btn btn-mini btn-primary';\n button.href = '#';\n button.title = 'Stop Interaction';\n button.innerHTML = '';\n button.addEventListener('click', function (_evt) {\n fig.handle_close(fig, {});\n });\n button.addEventListener(\n 'mouseover',\n on_mouseover_closure('Stop Interaction')\n );\n buttongrp.appendChild(button);\n var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n titlebar.insertBefore(buttongrp, titlebar.firstChild);\n};\n\nmpl.figure.prototype._remove_fig_handler = function (event) {\n var fig = event.data.fig;\n if (event.target !== this) {\n // Ignore bubbled events from children.\n return;\n }\n fig.close_ws(fig, {});\n};\n\nmpl.figure.prototype._root_extra_style = function (el) {\n el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n};\n\nmpl.figure.prototype._canvas_extra_style = function (el) {\n // this is important to make the div 'focusable\n el.setAttribute('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 } else {\n // location in version 2\n IPython.keyboard_manager.register_events(el);\n }\n};\n\nmpl.figure.prototype._key_event_extra = function (event, _name) {\n // Check for shift+enter\n if (event.shiftKey && event.which === 13) {\n this.canvas_div.blur();\n // select the cell after this one\n var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n IPython.notebook.select(index + 1);\n }\n};\n\nmpl.figure.prototype.handle_save = function (fig, _msg) {\n fig.ondownload(fig, null);\n};\n\nmpl.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.\nif (IPython.notebook.kernel !== null) {\n IPython.notebook.kernel.comm_manager.register_target(\n 'matplotlib',\n mpl.mpl_figure_comm\n );\n}\n" + "text/plain": "
", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMcAAAD8CAYAAADDuLCoAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABLH0lEQVR4nO2dd5hdVdW433XObXOnZUqmJJMySSa9QUIqSocACviJSBFU1KCComBD8RNFPgugoiI/6VIEkSZCQu+QkE56zySZZCaZ3m87Z/3+uDdTMj0zl5kk532eeebOPvucs++es85ea+219xJVxcHBoS1GfzfAwWGg4giHg0MHOMLh4NABjnA4OHSAIxwODh3gCIeDQwfETThEZIGIbBGR7SLyk3jdx8EhXkg85jlExAS2AmcBRcBy4DJV3djnN3NwiBPxGjlmAdtVdaeqhoAngQvjdC8Hh7jgitN1hwJ7W/xdBMxuWUFEFgILAVwJrhnTCsw4NaXvsLAxe/g+sVE2VGaRnVJNlhlqdSyCzabqLFL9jbjEIqwmtTV+PIkhRvkqccXuFVILj/SufxrUZtfBbLIzK8k0I7261tHCyrXBMlUdfKTnx0s4pJ2yVvqbqt4L3AswaHyWLnk5GVOOTf/AsmCYKx+7hTuuuIdP+5rL764axj9vOZ/iTyn5z0fYc7aH7HFVNGwbxJLL/4xX3BRF6jj3zz/ilevvIM3096odLzX4+P5T3+fSC17hhvSdvfxWAx8zd/vu3pwfr6exCBjW4u88YH9HlX1m5JgVDIBZXje/u/RRrl7yVYIabipfUTOSpGuKeP6zd2E2RohkhFl04n389wt34hU3AIMMF+5Pl+Pug/453x/gnS/dzmeS1/XoPEttLLXblNfZAartxl63a6ASL4PcRdQgPwPYR9Qgv1xVN7RXP33CYC17a1Cft2OgEVYLt5hYanNf9TBGeQ5yRkKQCBbjF32bM6du5L5hH/R3M9tQZweosCMMdyW1KrfUxkZx91Llixdm7vaVqjrzSM+Py+taVSPAdcArwCbgqY4EA8DW9rSw/qHBDrUpu7H4RH5ROqndt2dPOPQQ/bpsMr9793wa1IspBgGNkLzZjddotgV+WTqR2ytGt7pnnR3odRuOhCTD10YwAEwxBqxg9AVx02VUdZGqjlXV0ap6W7zu09f4DU+bspuy3uPZndPYEekbFeKqQcvY8Jm/clFiHQCpRgL+Mw7y5gszCKsFwOzEHTz9+7MZ987VWGrzdqPBvD/cwGuNCQAENczNB6dQZwf6pE0ObYmLWtVT0sZnafnbqf3djH4lqGGeqB3KlcklTfbXqw1url12BQW5BznwxAiqTwmw6dT7adAQJ7xxLQsmbuQvQz48pu213jAg1aqekmgG+7sJ/covSidRbYf4SsrBVg/6GQlBln3qb+z4YAQVc8JsOvV+DliNnPj098nIqHMEI84MiJ6N6LGrt3aHJzfP4O6KWa1UpLurhvHtffNJM/28fOXtvHfmn7Cx+dQb15MwvJYPT3iCB2ryOvUWXbBtAQs2n98vdsqxwIAQDp/R1gg+nojs8/PIkvmcte4KGuwQbzSa3PX8Z7gs4yMA8t1J5LmSuLX0REyvxerZj+AWk9+8ez4VltXhde8a+TTbVw/jlHUXOwJyBAwI4XBLx//geNKeZ6o/sJMtMJXat7M5/cff5Xdf+hL+YqEwnNmq3k8yl7PhlPuaPERXzfmQPFdCh9fNdyex8ot/pKQihX1WQ1y/w7HIgDDIC6Yk6JbX8vq7GXGlwQ616wkDOGPjBRysTcKyDNwuiwRPmPem/euYdpN+EvTWII9X+EiPqLJ6FxYxENgTqeOX+xeQ66vmMylrWBkYSaIR5LLkfXjFzTf3nslofxm/GNw2MPn6Ea9z/VtXcM/pj7DAf8g50XPBsNQmgtU0u+7QOwbEyJE9MV2L38zo72b0irBaFEUa2RlJ5R8H57O6JI9w2GR89kEmphQzP3kr9xSdxoy0Pczw7+J7Sy7lWye+w8JB69kfUWrVzQyP2S3v06G5kKJIIz6B3NgE3Z5IHReu/jrLZj7ujDocIyOH3W6c4sCjM9XILSb57iTy3RZnjHiXsrx6nq0r4PLkneyN2DxTcyI1QR+PbziJS+at4LdznuHmZy7nqW1nE/EL1137LLO8JR3eO6wWD1QP56uphU0jwyv143i3ciz/zH8LgOGuJE7J206DhkiVjm0Rh+5xXBvkPaUjwWiPTDORhan7STJ8TPD4uTlzM+9OeY7X5v+VZMNipLuMxH2Cp9amdm4jv37nAraG66m2G/lF6aRWAYphtRi76Jvc/XDrJTGDXbWsfHs8O8J1TWV/yl1BqtEzwXA8We0zIISjLuLt7yYcEY/XZvBAdU636zfYIT4O5TDclcQsr5sVN/2V1/74Fxad/Fc8aQG+sPrrmAiPLJ/LhNe+SXGkLnafXHLfMLE8tLInAupm+MsBzn38h/yydGKP2v7d/Sdx0qpLmPjhlxj75tf4f1VDe3R+TxkonsGeMCBsjoQxQ7T+/cT+bka7BDXMj0vmkuWuZWHaKtKMhCa74IclJ/Dc5mlsP/XhLq/zaoOb3xWeyzeGvcelyZVN5X+qHMm9T5yHWHDWxcv4U+4KghrmrPWXUPlGLrd+4xH2hjIY5T1AhlHPVI/VagT78YHpPLVkFunDqihIL+PJ/De79b2WBqKj9Rh3gGTDc0wa8ceEzeF1DdyVaQYGqa5GHvh4HvfVngIu5bEz/858n8GpKZt4tmR2l9d4pCaTR/fN5d6CfzLa3RzdWhSp4/E7zyVRlVCKUJBwAIiODu9OeY57hw8hx6zmorToCNJgW0xafC2bz/tb08P8u+w13PjZD0gzfGwIRYDujcJzfIcM9oH5UhoIDAzhMAaucLjF5JeDN/DLMzZgqU2jhpoezJN9lfzPqR91eY2rUsq4KuW/QOuw748CQ3A3KJWfrycUdFEd8XN31TCuHRRdYbwwdT8tXbqLGzJJXeemdkEIr9n8ps8yow/4dO/R66GqswMkGb6uK36CDAibI+koCTw0xSDJ8DW5SVONBG7PWX3E17sosYqFtzxLfmY5P5/5EgvTVnHn6+exNtQ2DP17xTP5+aNfIpwEye04BtYEg5yx8QKWBcNtjh0NDDTBgAEyctRaA69j4oWlNnUa5PWGbF6smEZhbTp3FfwLE+XDwGBGTCjh0hVfZ+O8x5rqz13zRcoK03n7G78nWQy84m9zzTfrJ5DqaSTHDALHnv3QHwyIkSNkH73qQE/5IGjwh/KZhNXk/4YsJsNXzw93XMzla67mhue/zPdHvkZ4ZzJ7Yp6qoEao+3Aw/3vGcwx3JbW7ycKGcIhH/r6AG/NeaXfFnsORMSCEw5Sjw8/eF/MBn/bBjzNWc2lyJbmuJP416lUGJ9Rhf5BG+sQyLkhs4IIzPuLDxuj+FH7Dw2cv/pDb/3ExlR0ED071+Djv6ve5+olrOWjVN5U/X5/Eqw3OKHKkDAjhiBwlI0eEvpms/OmBeTxeGw2XMcUg3VNPY7bNm9MeBeDO3FVN7t49kTr+L2sVnjkVzH7/Wx1e8/+y1zLu5F283TikqWxJ3Rh+/Idv8FLD8aO29iUDYp5j+OQULXw9u7+bccTsCtexKjiERCNIulnHSFeItBaG++E02CH+Z+vnWJC9gS8mr+fp2kkE1MUP03e0qhdWi8kPX8drV95OuuHi6brhfCXlYI/a9sOSE3j+1Tmsv/LPx+RcRmccE/Mc4aN8JeA7jaO4a8vpGIZNoidMY9hNdZ2PnLRa7hn3TyZ5WodzXL/vNLbszKXm3mHce/HJhHckYw8NcPr8TczwNnuibGzEAp8ISYavlWA02CFs7C69PLfnrOanV3zQxoh36JoBMXIMnTRI975xxLs2fmJ0Fnh4OGG1OGA1kmyY7cY6fRCw+dZfr6NhqI23zMBbpVTNDnHGhM2t9q7aE6lrZWS/2uDm48BwysLJ/HfHZDbOe4xKq4HlwVRm+2pIEq+zrjzGMTFypJpHxyq1ngQeusUk7zDPUVGkjgOWhxleD7O8ipow+l/1BDN8BL9bgb/RR3W49UgQ1qhQrgh5+Mbyq/CsTELmVXLlmGX4X05mTPCruLclEMiO4E4L8J0p7/CdtF7tgukQY0AIx9GqVoXVotIONM1Qd0Wu6eeJmgImezZjYHDKJStZeXoegZeSqSpOw50Q5uGRiwEPQQ1zZ/lkHn7pdOacuoH1/5jEqEc+RkMh6rfP4KHL5+L1Cr+c+QKP5s5hdkYhzz94ChNPKorvlz6OGBDjrx4l6zkOZ1ckwOl3/5BLdp7RrfqmGJyXtJ79kSBuMfnr0I84LWcb1ScFmDtuBxs+/SBecfHrsvFMfOo7PLJ5FmZQ2Hj/JJJKLCo+Pw0zJ5uiz1h8MPte7rnxL1ySdJBhiVW8UTyOiB9O9kV3Rbyx+MRWoewOPWdACEdA+8+LYqnNjnAdHwRsiiI9e5jGuhOZ+JktbPvnuKbVeV2R54KPQ81h7hcNWsmDJz/MoyPf4KWGVE5YfgX/veM0MlcJKYkBAjkRBu0MUjLbpHISJP+rka3n/J0SC+rVg1tMLspYSUPIzQ1XPQvAmRs/x7a6LPJc8VkKENRwt7/v0cyAUKuE/nMKBDXCXaWn8fK2iSwo2Mifhyzv0fkPjnyJnT96Cbd0by7BViVguwmrxe/KJ7F4/0TGDiplWO5i/lv+aX4x8SV+tPNycj9QvH9IofQSpfHHVfwq/y1ueusL/HHYC1yx60JWrBqDuhVcyuyJOzghax9fSy3h12VTmD94J78c/DFmHFy3ltrHjUt4QHirRk5O0p2v5/Z3M/qENcEgW8NZfD6xsl2v0a5wHYmGcNO+c7AR3nt/EukTyll2wr8BuHzXaRz8WT6NWR5q8wyu+OprvFc+hn3Vqfxx8lN8c+UVfHbMep7fPI3ZIwv5fu6r3Lrns9ww7BU+7Ys+vPHwVn13/0m89O4MvMPrmuK+Bjpx3w5URIaJyFsisklENojI9bHyW0Rkn4isif2c1+Kcm2KJMreIyDld3cMeGNpdr7DU5ht75/O9a6/jzzddyn8bUpiz5uJWxx+pyeTn+z5DiWWydN8IDJTsj5TGkJughvnG3vkU/mUcOy4z8VRHsD3w4PNnsue/+dhvpvPNx77Jr6a/wG+zV/LCvL+x4rWJbAgO4fmCV5qS4sTLjZtkBkkZXcWNk17v9bX2ROqOirweXY4cIpIL5KrqKhFJBlYCFwGXAHWqesdh9ScCTxDNCzgEeB0Yq9qxkpo7KU33vZHZ0eEBT1DDnLHui4SezKZsho272iA0JEziVg8PXnMXthr8du95fLxhBIm7XfzrW3cywiWsCPq5+UffoHySiZrAhFp0SxLzzlrPiqen0DDUJnOVkLqzEctrsvtcD578WgINHhI2+xAbLA+8+fXfN+1AcjRwKEYt3vMxcR85VLVYVVfFPtcSzbfR2YLjC4EnVTWoqruA7UQFpUO8MnAXO3WHz22Nbnzw0C1/YP3n/sKnz1pL2jI3asDlT38XQ2zuHfUsWy+8h/TTi3m6egZJho/5vjDF84XQIJucj8LkPujDiAjr/z6Z0y5djq/UIJApNGZ7CaW6yFyjjBlcxtun/pmCc3aQu6SRcLLymdt+yC9KJx0167RNMY6KicoetVBERgInAIeWv10nImtF5EERSYuVtZcss40wichCEVkhIivqK4+OxU4d8ciYf/POlKeZ5EnAb3i4Y+jr5F22i8ZcCyMCJ3gMssxE3GLy7pTnmjZ2c4vJc//zJyacuJshP9tOySw3kQQltTDIf1dOJ5ykZK4NIrZSm2ciCvVhD0/XTKX0ryMxghZZK5TEAxaWGj2apGyJpfZRI1ifJN02yEUkCXgHuE1VnxWRbKCMaCLMW4mqXleLyN3AElV9LHbeA8AiVX2mk2vXAlt691WOCzKJ9rlDx7TsoxFxzyYrIm7gGeBxVX0WQFUPtDh+H/Bi7M8eJcuMsaU3uuHxgoiscPqpc/qyj7rjrRLgAWCTqv6hRXlL3+vngPWxzy8Al4qIV0TygQJgWV801sHhk6Q7I8d84EpgnYisiZX9FLhMRKYTVasKgWsAVHWDiDwFbAQiwLWdeaocHAYqcZsEFJEFwF1E95a5X1V/20ndhap6b1wacgzh9FPX9GUfxSsPuUk0D/lZRG2Q5cBlqtp2/30HhwFKvJzNs4DtqrpTVUPAk0TnPxwcjhriFXjY3lxHq30zRWQhsBBAvJ4ZSa5MRow6iLfFuuty20Wd5WWEu57uUhLxUlOYBA1HHp4gIoSy/bgP1HMoJlISfDDcQreG6UmcpLhM3GOU+ooE3Adbfw/xekgYFcIlFuXlKagJqan1VDX48ZVE0GAodg0XGmmeKJUEH9rYs/zj9iA/OXkV7C/JwFXWvf4UEXC70FDnG8WJYYDXA8EQatvYaYmYdUE03HZyVzweIsluxIa8vIMkSc+XK6yrzSQ/qbTNuQctD+WhRCIRk2x/DXs31JbF3ZV7BLT3jVs9UjG98F4RWeAekrl4nusKFj73WlPieoAyq54bis7lkRHvdnqz2ytGc+2gTfgNDz8+MJ1lN52E55UV0AuVMTjjJLwvrwLbAhHs2dNJu20PNadWt3pQu8SC8IiZnPP0Oyy66TR8LzY77iq+NJcPb/0rrzYmMsVTxmNVM3jj4DjOyNrC+2eNwDpwELNgFLOf3sxbN83HuygaMWz6UrEC1T37QtVQc/4crvjnuzz96KkMuePDTqvbJ0+nbJqfIc/sIFJyoPVBkVZ9Kx4v5ZefSNYbRRAMseuvg0l5LonUHQ2wdG3rc8NgZoyi+KwcXrnpdoBuLxY7xJi3v8I/5/6d6d7WIfkvNfh4/MBclu7I5/sz3uB7E99o2rFbRGYADwMJwCLgeu3CpoiXWtWtuY6YbXI3gLWjkF///spWezMFVLG16yb6JEKdhvlF6STOTVlL0VU92BJTBNewPFwjh7cqTnhzXVQwAGPaBAAi9pF1l/vVFTz3uzPwfX8/ZmZzBqvKicqVhWdx+3evpNTy8ND6uZg/T+PVmz6NdSC6mYLUNfDo+tm4GpodflZVDwUjRso/l/LRpZMY+k5tl3U9uw6S+2YpkZIDiMsVFQjAHJSKMW0CrqFD2HHnHPb/YB6B06fialSqZg9FgyFOHbGdA5+yifjdYLRd5bnv/BzMYPS57KlgALhcFslG2//xivpRNERaRQkMjz1jAPcQ1VQKYj8LurpPvIRjOVAgIvki4gEuJTr/cTiziMZeYXi9VI2F/9SPbDqYbSbwoyEvd3mz76TtJstM5McZq9kaymHYg9H+EFf7A6ORmIi4o50YOnsGG3+ZQyRnUKs6dqBZbZHdxQTT3KR5uqeqGcnJbcpSH1uKeUkDVnkFRmIiNZfNYfAqqDk3gn/Jdr6y5iv41viRJR/jfal5TUmkuITRV6zGfHtVt+7dFdambejydV3Wi+zbj7VpGwAycQxmakr0/Kpq7DUbscsrGPKeTdJ+G/+SrZR82qZqjEHF+eN4/Y0TmDJhD76NRU3ntWToi8UEMoVVwfTWbVO7V9G6NREfflerMJggMCs2J5eiqktio8UjRINnOyUuwqGqEeA64BWigYpPqeqGdqo22SZ2IEDB49UUhdKbojbdYjKqB4qf3/DwqYTt7PmqjZmczLbbZ8KsKbhGDMOVE90Xy0hOpujaaZjDopufJewo5wdzXiGQGY35tk49EXPs6FbX3XzrWKq/Uttm84N2MUx0zHDMglGYKa0fDKu8AsRgyz3jOO1HHxJKFqyaGiQliS+OXoXRIsRs/w/nUXPZHMzBgzESo29XMyWlSeDNjHSs004Ew8SYOr77nXQE2Gs3txmt7ECAhOeXkfLPpdgFw/nmyW9x1WWv4bnyAEm74azMTZSfmc++L09ocz1r+y4S9yuj3BWtyk0x8MuRxYcdwpBWmlKI6DM2lKg2c4h24/3aXKtXLekEVV2kqmNVdbSq3tZBNQHI9tey9+Z5lE9L5YGPPkWdNj8lG8M923whw1Qm5+1HQyEQqPtVPQfOysOO/XPtujqGPbCZSGHUX2Bt38UfXj2fSEJUbYgkmGy9ZnCzGpGRzrj7ajhz2BZ8Zue2hisnG9fIYejqDWiCB7LbCcO3Lcb+JczH5+aS+fcl0TaUHOSBlfMZelFh06iTutPi+lv+ReM//UhuFge+O4+xbzXQeO6J0Y5zu9lxmYk5YQyyp7hHfdTX6Ir1/P2t0/lB+hauHvEBtlt4oWQqNfkGSfvb30LVU2eTZ7ZdUdibRJ8WBmWBNmqa0g0buD36e5lsETAsywxx6aUvs70xCwMlSZoNrVneni3JzDITGZdygI/Vy7jf7AQgMialSU0KnH8S3oog8mHsrSWC7bFp6j8DJs4sxBqSS2TffqrPGEtCaYgP75jF0G9v7/TekYNliBG9jjQEEav9B0OXr6OlmGkoxIQf7OLgRePIkugAm/jMR/zN+gLu60rY8otB/HHOw9yy6bPkvL8Nyc2h+O+pyG4D3bUXu6F/tzYyEhL43PzlNGqIV8snkbk2wPbJObgmNJD8m/XtPoXlk0xsjmzv4XDIRbntZfRh5SY2HrNVMIaHqK1bRNTuPUR34v36fQnecqLGEScnbuE3uW/xvzlvNsX6h9XqcPPkzqiPeBGXi6rTRmEdLEWWNOvYJVcG2HWdIDFPh8ycTPKQWtz10X+Uf1sFp2dupmpe1J+Q9NRSXO+txVNnH67PtsGVPRijIB/75OlY23cR2dW9/aP2/O9cyEgj44ElWDU10Xa5PSS9sw3vN02+d+Kb/LHwLHKuLsMeNYTAox5q6xIY//PN/SsYs6ZgZqRjNzbyZtFYkgwf/8x/i8LPejATI8wYsQemjcM1YlibU8WCj0O9U6EOx8KgJujD5Wl69fiAZapaDNSKyJxYrOBVwH+6ul6/CkcL24RZXjepRkIr74VbTFJ6mNQkrBZvPj8DuzFA2tL9UZfjIa+T388FBeuIBFxgRcuqCxJxvTKIxFV7AJBQmP/sn0bluObhXSMRjJBNWaDz1XaRkgPonv1UTOxBNlfDJH1uCYWXZOMa1vxy23f9TLb87zjsvfv5x5/Ow/f9BLBskv50gKGJVYz9cRl7F05qsqXaY+/P57Hnf+dh+OOzFagRjCCpKRheL5lJ9U07vOdP34fbEyHHV4OV4AazraqUvMdmmufI1pDYdW72htvmrXeLxeiUMlyuphFpd4u4vm8B9xN1AO0AFnd1n/4eOVDVRZ0d7+mKsV2RAAll0YE8Urin9UHbpt7yIrUu1FYQoWKikFxkNfnyI4V7cP06ndSdrYd82x1th9qdqKqqYFlE/N2f2DI8bs7M3cL/ffkRCv+UiuHzYZ88nRu//jSZqwUzN5uwX7A3bGHbT8azuzqd8s8lYJeWMeGCLdg5bR+SQwx7rY7TPruKzX+a1PcCIoLsKaFqZg4NZ07ltlHPNh26dvhbJLyRxOXpS6gendDuCDpoYy17I0emVhmN7T8TtgqlwSSyU2sxomktag4dU9UVqjo5ZgNf19UcBwwA4ehr8l0+ak+rJ3LadHT+9FbHjLRBzEreyciXIqA2VVfOIWulTdK6klb13Ot3c3DWYa5gAZ8ZbrIpOsPqwXZRqsr+wCAGmzU8PeM+8t+D3z56L1WWn1CSQMQiY1MIVEnfAINvaBbkjS+OQ3YW4crJbnJNt2LpWgq/MoIRI0sJzmvrNeoNOncqdl09Sf/+iMR3N/Or3Rfgjb3I5vkOkHjQ5rHyeVRNgLqL266Slu17+PneC47o3nZCx0Lld4VIcPVN6rf+Nsj7HAPB47HwLi/EbmhADROdPRlZ8jFq29xXeDIJIRsxTdKfWYtdX8/hPiirvIKxP1uH3WIm3HegsVsTknYwyIhHdra5ppmWhlVV1WbWXoNBXl89DXu6kOmt463CAhavnsKEnxWSU7sKTU7G++56bCDz/WLs/VFBtgMBhv72Qyyg5vyJpK1KxtrS1mFgbdiC73wPGunb/XPlw4+bDG2rpoZNuwrYP1KxzQayzETKJ5m8t2806VNL8S9Ojn7/yuYU03ZtLat2jKWNVd1LAhE3PtMRjnYxxaD+YCJWXT06ZzLGsg3Iik3Rf2QkQmlVEmPWFWJFIp2Ggdj1reOP7AQXPleYjiKawmfPZP/JbvzFEEkQGrPzcdcI+U/sJ7JrN1tvHoenwmD475ahMXvnkKBkfmRy/YLX+eJDNzDitmWorVi2hZGcjPi8WKWl0Tak+KHUjcTafaj9Kf9ajqU2RmIiRtogIkX7WrVNw/FdH25OHIur1MMET7PqFhgTQFenM+6szVSEEhFv25HNKOv7zeGS3EFCtskwd3mvr3XMqVWW2k06qSxZi0YiTQ+HNgbwLUuKTsbFcOXmEDljRpfXDQ5yk2h2/JD5lu9g5H9qSNpvUTfKYviMfcy/4GP0gRBb/99MUrcI7pMqKbpxFhVfmcOOx6dR/aU5uIYOIfFghG/+5Hvk/3lT9IGPORDs2loie5vnruyPN2HX1mIOzcXMyYY5UzEnjo3WV8WaXkDqvxowM9I7amZckIYARhBebmjWJ78wdRWJe5Usby2uqoa28VlA1vLotkZ9hSHKIHcjAcvdJ1vMHnPCYYrBOSevwUxKbKPC2I2NJO89bFGiaWK7urYjPNURGq2OO9yqrGTvOalEfAYF31mG91suVt8/Ffmqi4JHwjRmCwmeMOdcspSUPSHSX0ngJ794jCvfXMJDf/sjdUOMVmpHEyJgmJgZ6ey4fTau/BFEdu8lUrSPgzOT2PzN6KYv4nLh2ryHpSvGIcmf3B5WOn866jLxVgr+FlP8F6SuInNNDe8UjaFmYvvCOmh1KVvCfbtI1BCbDXtyORAe1Ptr9b45A4ughln88WTadUaokvTvj1oVRYr24Xmt61zitqfrrhr2SjVmWKNv8W07ybhvCZHCPcgHaxh264cMunAP6781CXdVAH9phJ3BLNxi8c3tl3L11xYRPnMGZkY6VVfNbQ4ZGTQIZk6k+vFB/P1z96Fmczv8Byx+u+BJdN40zOws7OG5FFy/tK2XrhPE7cHw+9uEunQLw2T3uQlgGOT9t4T/Vp3QdKjA3YiV5CH9nkQC6QbmoFTMQamtTteiYu4tPaXn9+2EiJqo9s2u/ceccLxYn8H4P9dh13YdedqEfdjbS6SVauIaNRLLa1DeIjRB3B4OXjuPXb+Z2xRqsvszqez/dMf/GDMznboRfnTVRrwvr+L1eXk8eMo89OYMllXlIz8pRRISSPhSMYUPj8KVNxSrspLa/ET+NfERtgVzoKxZJbQ8wqd8+/jyQy+y8eY8jG6GkRiJiey4Yw7Fz08g8NIQap/NYcePJnXr3FbYFr5SYfvVWdROGcwVaUubDmWZiZRNScD/8V7MAOy5ZhK5L1vISVOaT29s5KW1U9q7cud4LTLMtjvih9VkU1U2WZk1eNuJ2u0pA0I4wkcYRtAeBZ6DqLt3yXDE40HzopNr5ph8ttyaxoGZLvISq5pm1s2hOVSdEOLET28BMRCXi8CQMGZDW+EQr5f6z8+m9P5kGrLMpolJq6aGA+fnc/XDL3BL3ovsPpBBpGgfiVc1EqjzMPGF/UROn0HjYIPP3/xDHv3FZ9AWISmpTyznzP/3I6osPyNGH2TznfltXLr2KSe0CpMHsBsDjHwpRNafEvBf7yblS9Xk33pkUb9DXy0lcUIl7jqLb2y4sqn8oFVPXZ5iZ6dz8NQweW/W8tbWsRT+EMwJBdFQ9tlT+NTErT2+pzshTHY7wuEWi0mDSkj2Bkkxer8X74DwVrn7UEbXBYdi1AZ6nRRZ3dFo1503u7GqXQyZU8w3st7hO/+6jKS/pmK8vY7c14Zgj48KgzF6JPnPKL6VW9Dk5FYjV8k1M6iZEWDMr/3IkuZFRnLSFH75k4c4O6EegwS866Mz65GSA4xbWMnbV84hVYNk/TV6jisnm8BzafiuTCRSXAK2Rd5vP+Kl5+bgLykj/4QMOCxXetmUBIZsP2zixbYw34oKQ2/7yd6+m+qqyVDg4Z6JTxDdTwNCqmRPP4DevIeJ/5uFdbCMsV8DMyeLwNhsvMYYdM1WPnrzBKyvvN2jyd5IyEWF3TZywlbhYDCJ/VUpGMN6/8IdECMH0Cq5fG8Iqwlm97+WmZKCWTCqVZkGg+iK9ahpMnd4IWMeieD+fTp/Kj6LrO8E8b23ETsQIO2DvWx4aRyojbVlO+5XV2DX1FF9XmsVJeeeFRR8ZTWy5ONW5RWTkjjRU8b/q4rePzK9jpLvzwMRNBwi/cElTQ8xIpSfmc93RryJJrWY7bYtrE3bsCorcb25so17Ouehj4nsP7KoXVdOdpe2iIZDDHnBjadGufLp66izo87uPFcS3x31BkZGOlZWGmZmOhoOEdlbhOuNlVgbtmAHAuQ/V0ux1cDWcH23kwdJuYeSyKA25YYoaZ4GGuu9x5ZBntFOxtUj4YPqAiit6LjC4WuWPW6stA5Wo63fRmkgCTUFFfAYEay9+5rmQCJF+8j7zYdtJ/YO0+o0HGp3ye7gl3cyb9ENFAYyGPPSNYy4S/ifr7xNcEF0wz5xuXANy0O8Xvb+bC7fvvlpyiNJSKCtS1lmTGL3L+cROqf1Zn92fX2ny4XN7CxcuTntHtNwuHlOphOSt9cgl5WS92aE3ZHme33Kt49wwRDMknLsrLR2zzX3lbElnEqe6SbX9BPUcJcLnoxgx3ZddyZqu8uAEY6+2nX7iswl6NCs6DUHD26lg4vbg86Z2uYc6SBeSsMhGm8dQvFcH4FMN25pMVTHlte2h7TzPInLhSs3BzMtjfDZM9n127m4/wV3nP4kH945i9FPWFSN9fP4xpNouK6qKXRFkxLY8rcpnHjeRh7/ynk8e/6cNpN84vZQOiMF2614y3qmazfMGMHey0e1e8wqr2gzGdoe4fQEXIaNpyrIbfua0rTgE4OqggSs0jKqJqW2e65dWcWiqmn4DQ+mGHjF3W5q6lbtSuhY2A3pO/t1QAhHQO2m1X+9xRC7ySAPj8/DaOE+1HAIY/lhCxItG+nE1+7buI8RzxzAWxkh3OKtJB4PBxa0DcXGEBoz2+lWMcDjpuyC8dg/KMNXKgS+l8kP3v4ilRfVM/H2dfz+538n6V0/XxixCnv25OgEptfNL+a/wIbHJ8LStUR2FrYZCax5k7j0u68yYnEAXdnegsuO8S/ZTt4j23p0Tquv5fXiqgvT8J9sCj+TyMq3m1clLg+m4gooGomQ9sbOds+3AwH+u7VnHiv//vYf27CabKjouwxhA0I4Ki0/+45g3UZ7vFs3HtkSjSMy3lvdFHoBsZFj5sRW9a2qKnTjjg6vFxk2GCprUFMItJgE1GCQzIfa5g8U0yTUzktSwyEiu/eS9shS/N8Whty1DF29gfHfXUeo1M9vcz7gmhVfIuvBVTz473OoHpUQjXwNhrn9Hxfz5W8vinp52sH8YB0PPn0Oebdvb7XEN3jeSdReOqfD7xY92aD6tJ4FOBnJyZQtnIvh86HBILJhB7X5EB5k4y1vVnlO8lZjfPkg4nI1bRjRHnZRzyKGpYOBwy0W4wZ1fJ+eMiCEI9cVYHgfZSZaXzsEbWxftdBwqM1WMa4huYRO6fjN5dpbCumpmEGbzWVZTSHrxrQJmHlDWtWtu2QO4o+qBGZmBsbk2Fv0sHmTTTcMZtctJ4FhYgcCJO00+emBeRibktBwhOG/WsKgx5aCKtIYxF+iPLRtLtV/tDB8bb00Gokw4rYVvPfeZErnZzWVl09yc+Ov/knD/8xua2vFsMrKSXpqabvHOkJDIbKf2sje608EEYyMdL549vtkjSonaZ/dFBKSYvj4fN4ajLT27Y1DJBf2bNIumN6xWmX3YdruASEcfYnRw8y0muClcXDHYSGR4hKsrTvwlDcyfFBV8zLYkvKmB86YPJ7I6TOwPGCNGgIKWlePvSGacsRISKDoKzFBUWXi7QfIn7eH0oWzcI0cjq9c+eiOmQz/1UdNcVKoYg5KJefJCkIXVjHkFiH478MWNok0eZM0HGL0j5aS/nDzgz70zyv5+aNf4sxfvIcrr8v9BJove9KUDoUJQKeNRQNBkoqiqrBVXMIT62ZyYE86DdkGuyNRh4EpBlelriM8MQ8zLQ3XqJFtrmVOGkfV1J6lbg4P6rhusqtnm911xjEnHOmehqh+310OlpO2srTzOqqE03wkuZtjh6wDB5sW8ZTOTWPnJSYHTrbZdqWfceduY8dD49h6/wzCZ88EEayYX8Dw+dj442yuznufl356O+Y/Qlxww1uELq+AWZNaP5RuD8MTKqgtSkFFyHhgSdNaeDMlBTMjndCMMa3a2dIe0WCQYb9ewkcXjyeyr/vuXNm4s/MN8Zatww4ESH08OrppJIJWeEGUmhODrAs26/0+Mdl7egJkDGLLta29YkZyMrsuzmDN+X/u1cYKhwiryZry9p0kR8KAmATsT+xgEKOypst6RsCiMuBHaFs3474lZNwXq+fzsevrJzL67QqKFiSSsGUvtiqe6qjxaofCjP/eWv5w8WWc/ps7yfNX8fmUVVwydSV3/vlM3vvvXIbdGp30s8orWPL1GRiXGNSNTiJxNdR9YTaBdINAutCQH2bct1djZKQTGTcMd3EVVlFx6xD1WJxXj/qkGx6qwxmxKELFt+r5asESTkkoBqLu8STDhz22Hi0+iHt4EkZiYtP1g3PHYQbB7CNVyFahurFny6o745gbOXqKBoNYZV1nEjMbw13qs0ZiIpgmGesDSCjM8KeLiOzei4YjpG8OUnvhCZR8dzZ2IEDGG4Usqh/B34YuZZjL4ItrvsYHRaOwp9Q27xJoW5j7y0neLVSNNtlx+1zOvvk9QguqCWTZ3H/Gg9S+OJyxr9Zg/aqS0P02lZd2HX4fDzzVIZJ8QUyU8GGjjhiKhsLYW5IwMpttL19RDe468IqbTaH+3UGlPY5J4bBnTuhUZ26Ja1gewXO7zpJVPi2V7ITWwYxmSkp0Y7VD9yoYgQzJpnKcF2vbTiJ7YvMRauOuDpD62hZcZ5RF9X+Xi1/99wvcXRV1B589fDMJnjCjv18eFYrsLPb9eB6ffnkbky7fyH++/Xu2X3EPvxi8kcUz7uU7Z79MjllH3eIcXlg6g7IX8zg7eyO+yvjmCTL8/qZJQ3PwYMq/MTfqzi0qp6Qslfu3zWO/5WmyIYIa5jNj1mMOyWbwaht1Nysr1satDF5VT6UdIKcHWpX42g887GuOSeFwbd6DmZnZ4XagLbHLK0hc2+UWRvjLIgSs1teLTM7nc3e/xvZHpsOsKey8OJW6iZkE0wQzPQ1zfNRFqpEIsrkQq7ISeT6DstOGE9lbxOgfLuHO187n5gMn8+8PZjP4yqjts/uX8zjtjR0MOq0EQ2zuHr6YZENYGQwR1DBf3no5r1w2h+uu+Q5D/7EBX7EJp1Ziq0HSit0Yfj/FN85rFQHbJ4hQfsk09v/PKDBMrLIysl7YgYZCRIr2kfG6j9oqP89UzWRZbBbbK24m+vejHjdJL65p89Jy7yljXSiFNLMH7tx6F+VWW+9mWTAJ0zjGJgH7knrLgzY2UnXG6G4JhzE4g5Lzh2NMn4gxbUK7rlKAuhwXHsOKRrjG1B5ZspY/vfAZ3jnlL9SN8JO0F8KJBmoAOYOhuNnQP6RnD15RRdW45uvm/yfM+h9MZfzPt7LnmgnMXLSbJVffwffStlLVkMDGuiGYRPeVrbITmPTkd/BeESAwJBl1CYjB4I8jFGSU8swdZ6KqbP7rRGonBzF2dS30PUKVtIeXMHhNA+aYkVFjfHA6ZnrUVZu+oQ5x2ZyavImZ3uYR7Espe9nxlWykIJ+dX2phlBsmjeNzWFLf/vxNh80w23cWVAT9jM3owrnSA44p4QhqmAQzTONpUyibLq3CuztC3S6qxivbfuRhyw0JGEPajzMyQ9AQ8WBnpiGHVANVCu7czlev/A7Jz68mY20DroCCgr11V3RDhcMRQSzBHJTKwWvnsfMSk6LTvISm5VM/PMLzD5/CO4EsKu0Ab550L6P9pdxePoNvvXUVd+87ndH/qmfT/w3n13/7O+N+sZ6KxzMonepi5ZaRpD2ylJ3fHoOn2M34azdhlR3ZOmpz3BjMDuYmxO2heL6f8jnR+RR7/ebmZcdrtpCwPoFnKmbSYDevpwirRSRB0c3bGfVoC6+ZbWGGbVZVtRNp0AnisxjsausYURV216QhhuKWHqSJ6IBjzlu1eN1kxr+1jqSRJyA+b5ebC1jbdzHm+7ua/u6oS5OLgrgMi8Dm7a0iX63SUox3SlHAtbOYhqmj0U56tWJKKp4a2PqziVhpIeZO2MGmsmzMW7Yx9m2LfT+ex7n+SrySyO/KC3jx96dS+7lavj33TR677xxm/eVjdr01lbv2n8XK5QWkrReCk21y3nCBKvl3rMOub8A+fAFXD7AL97ZOPCNC+dfmUHNGA+EGNyP/HcZd3bZfNRxixJNFbF82kTlnTOftK28n15VEQC2GTjoAYmBt39XqHNs0WFM4jKV5FlM9Fn6j610QNWhQGkmBwzyHeYlVlLqSqKrzM9LT+3TtXQqHiAwjumV7DmAD96rqXSJyC/AN4NA49tNDG7SJyE3A14guF/iuqr7S65Z2A1dsLYEdCDD4niWo19smK1KPMUywLcRW7E6WXxrJyWz+eT5LL7yDkCrnRH5Ew8gw465dHbU53B7qLjyByCXlpD6UxpDzi9i6MY8/DX+R7/EZqtIHsefr47h/4V94ui6Hi5NKyHZX8+2bn+a8xN34xc3FN3zM6YtvYOyzdXxsjeV7Fy3iH+vPY8z3myf+erQCsgM0GGxdIAa2CyIVPjLWGHheWdrhPEikcA+uPfsYHp7Kb845jROTdrOqbgQlq3IYRcxBMWsKtteF8d5qvGsLoWoMo9wB/Eb3cnWI127XIN9Rk0nQMnG7Iww2e98P3Rk5IsCNqrpKRJKBlSLyWuzYH1X1jlYNF5lINB/HJGAI8LqIjP0k0i2bYiANzbk5Cn92IjlLIngXt42B6i7lX5vF4EdXYwM+M0JweB7WnqJWAmdmZ1H09wyWzvgDN+07h7yESh5aeBc7Q1n8/YzPYwZtChfazBqxlaEJVazdPwXzy8Kw6Urh+R4mJhez48XB+ELbKAxn8vOlF3LKaX9mpLuMIa5akg0PXnHzeMVUTp2+iWkP7eWkhJ0kShjvD8I8v+RkrI09X1HXbWyLwf9vCd3NH2bPn0o42cU7j53EW56TMEKQFFAiJ0/GfGsV5rYiTNPAIjr3kzS8BjdCcaSO3G6EEWnAjBnkrUeOqen7eKdoDLZtEOhs+O4mXdocqlqsqqtin2uJ5tvoLBbhQuBJVQ2q6i6ie5O23fIuDuyJ1KEJ0TRlaisjn6sh4f3NvbpmzWjANHFVB6PeqsPfmIbJljuHkptcy+yXvs+OX07g3/8+hVleN2f6i7jmz88QTnIx4gGTqv9x8/SKmZRNSySyt4jEJTtYXDONzXU53Jy7mDtHPMcU736Wnv4XTn3vO6wJDGdzaDABjWCpzdnJ63h77Xj+8soCvr7yy0z3esly1XDuvz9i/w/m9ep79iXG+2tIeHMdQ96swFuh1A9Vct4tx/VhNGLYqqxssoe0tg5Vocq2uyUY3SESMdr1ZvWUHhnkIjISOAE4tIXHdSKyVkQeFJFDFlx7yTLbCJOILBSRFSKyorS8d4NKUMMENcxP9n6W5I2epjXaunpDr9UMd3UslioUYW9NGpHdrUcNbItxt9biuiKMJERIKK4n/7Eipi27jGTDQ42dEHWvBi2sAwcZ+81VZN4bVYOs8gre+94cliwbzzn//gEHrASeqJqFpcrUYUUMc1dwQWIDqUYCy4PKEDPI2NHFJBww+OMJTwFworeEy1I29mgL0jYYZnT/q75CFTsQwF67meAgYeKsXTTkp7ZV14i6uS/MX0e+u/cP88FgctT2c7WvdvWUbuchF5Ek4B3gNlV9VkSygTKiSUBuBXJV9WoRuRtYoqqPxc57AFikqs90cu1aYEvvvspxQSbRPnfomJZ9NCLu2WRFxA08Azyuqs8CqOqBFsfvA16M/dmtZJmHsUVVu56mPs4RkRVOP3VOX/ZRl2pVLNnHA8AmVf1Di/KWS64+B6yPfX4BuFREvCKSTzQ5zTIcHI4yujNyzAeuBNaJyJpY2U+By0RkOlG1qhC4BkBVN4jIU8BGop6uaz8JT5WDQ1/TbZujxxcWWQDcRXQjo/tV9bed1F2oqvfGpSHHEE4/dU1f9lFchCOWGH0rcBZRG2Q5cJmqbuzzmzk4xIl4xVbNArar6k5VDQFPEp3/cHA4aohXbFV7cx2zW1YQkYXAQoBEv8wYP6ZvM4sebShKSG3CGCR1cy2KQ+esXBssi7sr9wjoMil6TC+8V0QWjB/jWbzslZ5FZh5LhNVi3DPf5sYzF/HllG0k9TCDrkP7mLnbmxKeiMgM4GEgAVgEXN9V0sx4qVXdmuuI2SZ3x6kNRw1htRj+is3ig5MdwehbhseeMYB7iGoqBbGfBV2dHC/hWA4UiEi+iHiIBiK+0E69WURjr45rvOKi9pvVfCm3Z/tHOXRJEJgVm5NLUdUlsdHiEeCirk6Oi1qlqhERuQ54hagr90FVbW+fysNtk+OWb495l9P9RRzatcOhTwgRfcbCRLWZQ7Qb73c4cVsJqKqLVHVsLCn6bR1UcyxPoqH2RaF0Tn7/263KK60Gnq9vDshrsEN9tqfwcYTSDRu4Pfp7JeDhtslxy82Z6zllVuvw+jTTz5ZALhvca5nkSejWKjmHVniI2rpFRO3eQ3Qn3q/f15AvJ2ocHfeYYnBqQttR4Yb0zWwLH7E38njHByxT1WKgVkTmxGIFrwL+09XJ/SocqhoBruvPNgx03GLyWX/XOzI6tMvuFnF93wLuJ+oA2gEs7urk/larUNVFM6c57suWhNVqtXdsXyX2OQ5pequo6gpgck9Odnp9AFFm1TP348+zJNj7TZVb4hjyR0a/jxwOzWSaidw/4TFSDQvom/XUgGPIHyGOcPQxh6tEPWWSp28Shzr0HketigOWNmc3cjh6cYSjj3GL2ZQVtTMstVkTDBJWq0dZjRw+ORzh6CciWDxRNZtn6jJp0M63LAV4tcHdxqiuswPsCMd/K/7jFUc4+gmvuPld9houTa5sk3f7oFXPM3UpTapZWC3GeyqxD4t4CKjF3kjKJ9bm4w1HOAYgWWYiAXXz+W0XUBypwy0muWZCG0M/00xsd1bdoW9wvFUDlCuSyynwvIA7tiqwLxJKOvQMRzgGMLO8bqBzw94hfjhqlYNDBzjCMUCoswO9cuk68yp9jyMcA4TXGzOP+FxLbR6qHkm13diHLXJwhGMAUGcHmO0tOWKj2xSDLyZvY32oN3kIHA7HEY4BwPawUGX37l+RZvqZ73P+nX2J05sDAEOUCttZ0zLQcFy5A4BxbpMGOwD0IFG9Q9xxRo4BgFfcpJl+6uxAj85zPFTxxRk5BhDt7XYYVgsDwRSDBjuEjY2Bgd/wUGuHcBmms4w2Tji9OsD5S2UBHwSj/6Y6DVNhR6iwo1G8AVWC2py8c1nQGUn6EmfkGOCck7SBHeEMoIEss/VuiHmHpSaOhps49BXOyNFLqu1GKq2GuF1/rNvDS5XTnA0S+gFHOHpJknhJicPO6IdCSdxi8vshb7axK8Jq8UHApjhS5xjmccIRjl5iinFEBnFXI4FN83G/tN09ZFM4TLoRINNMoNQKfiKhI2VWfau9e491upNqeZiIvCUim0Rkg4hcHyu/RUT2icia2M95Lc65SUS2i8gWETknnl/gaCSoYXZEOn+YW65Bby+sZILbTa26cYtJniupzWrCeJBmJHCK72Dc7zNQ6I5BHgFuVNVVIpIMrBSR12LH/qiqd7SsLCITiebjmAQMAV4XkbFOuuVmvOJmrLt7xvPdVcOY4StkhKuR3BYGuFtMZnk/2QVQphikmcfPRGWXI4eqFqvqqtjnWmATnec2uBB4UlWDqrqL6N6ks/qisccCh+82ctCqb3W8zKpvOr4jXMd0325O8NrkupKoswOt7AtL7Sb1zDHY+54eKcsiMhI4AfgoVnSdiKwVkQdFJC1W1l6yzDbCJCILRWSFiKwoLT9+BhW3mK3UJK8YvNzQHE2baSY2Hd9nJeEWq0nFSjJ8uIgeezcAOyKNbA0HeLgmi5oezq47dE2385CLSBLwDnCbqj4rItlAGdEkILcCuap6tYjcDSxR1cdi5z0ALFLVZzq5di2wpXdf5bggk2ifO3RMyz4aEfdssiLiBp4BHlfVZwFU9UCL4/cBL8b+7FayzMPYoqozu9vo4xURWeH0U+f0ZR91x1slwAPAJlX9Q4vy3BbVPgesj31+AbhURLwikk80Oc2yvmisg8MnSXdGjvnAlcA6EVkTK/spcJmITCeqVhUC1wCo6gYReQrYSNTTda3jqXI4Gum2zdHjC4ssAO4imk32flX9bSd1F6rqvXFpyDGE009d05d9FBfhiCVG3wqcRdQGWQ5cpqob+/xmDg5xIl7hI7OA7aq6U1VDwJNE5z8cHI4a4hWy3t5cx+yWFURkIbAQwMSc4cfZENmhb6mlsizurtwjoMuk6DG98F4RWZDoyVg8K3IKxMn+cTg+eV2frjz0WURmAA8DCcAi4HrtwqaIl1rVrbmOmG1yN5btCIZDPBgee8YA7iGqqRTEfhZ0dXK8hGM5UCAi+SLiIRqI+EI79aK2ieV4egc6RnJyfzfhSAgCs2JzcimquiQ2WjwCXNTVyXFRq1Q1IiLXAa8QdeU+qKob2ql6uG3S74jbg4a7zrR0vGHX1vZ3E46EENFnLExUmzlEu/F+hxO3NeSquoiobtcZ7dkm/YojGMccSjds4Pbo75WAh9smDg59iYeorVtE1O49RHfi/fpdOJYTNY4cHOKBD1imqsVArYjMicUKXgX8p6uT+1U4VDUCXNefbXA4ptndIq7vW8D9RBff7QAWd3Vyv+9bpaqLUiS9v5vhcGxSc+iDqq4AJvfk5P5WqxwcBiyOcDg4dIAjHA4OHeAIh4NDBzjC4eDQAY5wODh0gCMcDg4d4AiHg0MHOMLh4NABjnA4OHSAIxwODh3gCIeDQwc4wuHg0AGOcPQEGXALFx3iiCMcPUDMTzaTkkP/4ghHD9BIpL+b4PAJ4giHg0MHOMLh4NABjnA4OHSAIxwODh3gCIeDQwc4wuHg0AGOcDg4dIAjHA4OHdCdVMvDROQtEdkkIhtE5PpY+S0isk9E1sR+zmtxzk0isl1EtojIOfH8Ag4O8aI7Ox5GgBtVdZWIJAMrReS12LE/quodLSuLyESi+TgmAUOA10VkrJNu2eFoo8uRQ1WLVXVV7HMtsInOcxtcCDypqkFV3UV0b9JZfdFYh84Rt6e/m3BM0SObQ0RGAicAH8WKrhORtSLyoIikxcraS5bZRphEZKGIrBCRFWGCPW+5Qxuc3CJ9S7fzkItIEvAOcJuqPisi2UAZ0SQgtwK5qnq1iNwNLFHVx2LnPQAsUtVnOrl2LbCld1/luCCTaJ87dEzLPhoR92yyIuIGngEeV9VnAVT1QIvj9wEvxv7sVrLMw9iiqjO72+jjFRFZ4fRT5/RlH3XHWyXAA8AmVf1Di/LcFtU+B6yPfX4BuFREvCKSTzQ5zbK+aKyDwydJd0aO+cCVwDoRWRMr+ylwmYhMJ6pWFQLXAKjqBhF5CthI1NN1reOpcjga6bbN0eMLiywA7iKaTfZ+Vf1tJ3UXquq9cWnIMYTTT13Tl30UF+GIJUbfCpxF1AZZDlymqhv7/GYODnEiXuEjs4DtqrpTVUPAk0TnPxwcjhrilROwvbmO2S0riMhCYCGAiTnDT0qcmjIAEIEuRmgxjWg9lwsN9P28j4iAYYAqatt9fv2BSC2VZXF35R4BXSZFj+mF94rIAj8pi2fLGXFqygChk119XDnZbPrZSIYvsvEv3Q5JglVeETtP0LlTkSVruxSwrjGR6eMxduzFqqnpuvpRzuv6dOWhzyIyA3gYSAAWAddrFzZFvNSqbs11xGyTu+PUhqMGq7yS007aQCjVxKqsxKqobD6oirFycx8IBmBb6JqN2MHjJiJheOwZA7iHqKZSEPtZ0NXJ8RKO5UCBiOSLiIdoIOIL7dSbRTT26rhGI2FWPTqVlG21sYLWgqB9+TCr9u31BjZBYFZsTi5FVZfERotHgIu6OjkuapWqRkTkOuAVoq7cB1V1QztVD7dNjlsSSyyM3SW0mhDqhq3i0Ckhos9YmKg2c4h24/0OJ26LnVR1kaqOVdXRqnpbB9Wc/TUBVHEFbPZdOa5VsTl6JDpvWvQPEazTTsTw+fqhgUc1Sjds4PaIl0HeXQ63TY5bvItWkLdpBC33VLS27yI85iQSJo0jnOnHu6uMSCDQb208CvEQtXWLiNq9h+hOvF+/L5NdTtQ4Oq5xDR0CYhDZWdjmmPf11QRzknCt2k6kcM8n37ijGx+wTFWLgVoRmROLFbwK+E9XJ/ercKhqBLiuP9swELDTUnANb18F1kgE15ursGtr2xwzUzqfGzIHpWKOG9MnbTxK2d0iru9bwP1EHUA7gMVdndzfIwequqi/29Dv7NxDaEQGZkZ69O/DUx10YJSHThyDuDrWjK3qGqzUBMyxo/uqpUcbTZM5qrpCVSfHbODruprjgAEgHA5gNzTgqmxE/H7KFs5FZkzCHDemy2WvnvW7O9/5XRWWrcNO9iFebx+3+tinvw1yhxi6ZReWreS8aqCNAawDB7s8xyor7961V7bnRXfoCkc4+ghxuTBiapFW12D30KskHg8aCjlG9wDCEY4+QiMRrNJyjAQfRkY6hm2j9Q3txzAZJobHjaqioRCoYtfV9XjCz/D5UMt2NlaIE45w9CW2hV1fj11fDyKIp63NIC4XalnYwSCu/BFYGcmYxRVEivb1+Hbi8SBeL1ZpaV+03uEwHIM8XnQQw6SRSDSYMCGB8rk5VI1NBFeLXIOG2a6L1pw4to0XSy2LwLThzqx5nHBGjn7Cbmgg9fGl0c9+f/MBtdH8odjJozCXb0KDQcyMdGyvG3G5W6lQ4vFghm0iJ03AtWZ7u3MhDkeOM3J0A8Pvx0hMjEuqZdfQIZRePg2M2Oihiv3xJsRSir81AyM5Ga1vQDbtAG29SMmqrMR4dw2uNdHA5s7mPBx6jtOb3UBDIYzkZFzpaVglB/vMABa3h4pThpO9eDeWaaJ2c0yuLPmYvAMjsS2rlefLHJSKVVXdonHqjBhxwhGObqCRCFZlJVRWdl25B4jPi78kTKT4ANhtdy86PNbK8Pmwaur6tA0OHeOoVf2IXVuL661V7QrG4bhycwjPnfgJtMrhEI5w9DexuQ3xejsM8TASE7HyBuN6f32HguTYG32PIxwDBD1hHMawIe0a/RoKoyvWd2jrGD4f5VeehNHS6+XQaxzhGAAYycmYNQGobn+WvCsHgB0IMPj5zVAwotnrRTSk3ZkDOXIc4RgAGKnRSb92Z7pjD7uYZttjLbAqKzHKazASmoVBIxEYO7LP2nm84QhHf3LoLe9xdxhXJdPGR/eu6iw0PUakaF80dCWG3dCAvc5Je3KkOMLRXxgmZkoS4vViJydgJfswkpPbVFOvCdLOv0mke5OSzu4lR4wjHP2FbWHX1WPmZqOGQSjdhwzLjc7Et2TZhlYeKnG5EK8Xc/wYGi88qbleHGbvj3cc/18/o4kJiCq+DzZjtRO2LqaJqo2RkAAFI1DDwCyrRiuqifgycOUNxR6UDDv3YKSmoI2NrWfQHY4YZ+ToRzQSwdq0HaO6Hru+oX1P1YzxuPKGIvnDqJySSiDXD6aB+Lw0ZhpoXR32+s3YDQ0cXJCPWsfHJtGfBM7I0UvMQanYdfXdMpg7IlK4p41gGMnJ2PUNGIEIjeOycdVHSF+8FcKRaLxVfT05D1VgxQxwIzERNXHirPoQZ+ToLbZijDvy3T2MBF+rRVHicuEaORxr8qiowV64j4ZsN+79FVgVldjBIHZjNBDRbmhovpBlkfWc45nqS5yRo5dYdfXI9sIjOteVkw2AXVvXtDBKLQuraD+uqmrsxgBGgo+M9/dhlRyMrjOP1TMmjwcDxFKkqpbIvv3g7IbYpzgjR2+xrbYr/rrhOXLl5oDLRahgCOJufkeJy40xKBVJSsLMzIjdQzEzM1rdpyE/BbEUe3shmuRvrhsnzIJR1H1hdvPevccB3Um1PExE3hKRTSKyQUSuj5XfIiL7RGRN7Oe8FufcJCLbRWSLiJwTzy8w0BCvFzO1850IjcREwiOzsUoO4F5f2ORdErcH8bghEkH9PnRQMpKYiHWwFLuyqtU1EhavonriIAyvF2vL9m5v03OkSF0DSXsbcW3rcovZY4buqFUR4EZVXSUiycBKEXktduyPqnpHy8oiMpFoPo5JwBDgdREZezykWzazs9DaulauVMPna71YKSMdq6ISWboWVQVDcI0YhpU1CNlciDEoFbu8AqmsQevriUwvgJFZuIvKoUyaZsA1EkEUgicV4HprTbfC3ntDpLgEig9LkXCM0+XIoarFqroq9rkW2ETnuQ0uBJ5U1aCq7iK6N+msvmjsQMc6WNraSDZMqi+a3jqc3LIwBw1q8k7ZdfUQDFGf5ycybTTYNkZ6Glpbi5E9GCNs4zpQjdbUUX/mJIyp45svH1Ysj0H512Y5EblxoEc2h4iMBE4APooVXScia0XkQRFJi5W1lyyzjTCJyEIRWSEiK8IcI5mGDnfHetx4aq1Wbl6rqjq6qrAFVt5g/M8vw71+F9rQGN3vasZ4AvmZGPVBKK9CIxGS3t+OHKiIrv2YOZnkFUWEk02MEM4MeRzodh5yEUkC3gFuU9VnRSQbKCOaBORWIFdVrxaRu4ElqvpY7LwHgEWq+kwn164FHD9k12QS7XOHjmnZRyPink1WRNzAM8DjqvosgKoeaHH8PuDF2J/dSpZ5GFtUdWZ3G328IiIrnH7qnL7so+54qwR4ANikqn9oUZ7botrngPWxzy8Al4qIV0TyiSanWdYXjXVw+CTpzsgxH7gSWCcia2JlPwUuE5HpRNWqQuAaAFXdICJPARuJerquPR48VQ7HHt22OeLaCJGFqnpvf7djoOP0U9f0ZR8NCOFwcBiIOOEjDg4d4AiHg0MH9LtwiMiCWAzWdhH5SX+3pz8RkUIRWReLVVsRK0sXkddEZFvsd1qL+sd8DFtsgvmgiKxvUdbjPhGRGbG+3S4if455YTtHVfvtBzCJpr0dRTSh+sfAxP5sUz/3RyGQeVjZ74GfxD7/BPhd7PPEWH95gfxYP5r9/R3i0CefBk4E1vemT4hOJ8wFhGia5XO7und/jxyzgO2qulNVQ8CTRGOzHJq5EPhH7PM/gItalB/zMWyq+i5QcVhxj/okNieXoqpLNCopj7Q4p0P6Wzi6FYd1HKHAqyKyUkQWxsqyVbUYokGgQFas/Hjuu572ydDY58PLO6W/VwK2p/cdz77l+aq6X0SygNdEZHMndZ2+a0tHfXJEfdXfI8eRxGEds6jq/tjvg8BzRNWkA4dCdWK/DyUoP577rqd9UhT7fHh5p/S3cCwHCkQkX0Q8RBdJvdDPbeoXRCQxtpgMEUkEziYar/YC8OVYtS8D/4l9Pp5j2HrUJzHVq1ZE5sS8VFe1OKdjBoA34jxgK1HPws/6uz392A+jiHpaPgY2HOoLIAN4A9gW+53e4pyfxfptC93wvhyNP8ATQDEQJjoCfO1I+gSYSfRlswP4K7HokM5+nPARB4cO6G+1ysFhwOIIh4NDBzjC4eDQAY5wODh0gCMcDg4d4AiHg0MHOMLh4NAB/x/KQQwQ66qlRwAAAABJRU5ErkJggg==\n" }, - "metadata": {}, - "output_type": "display_data" - }, - { - "data": { - "text/plain": "", - "text/html": "
" + "metadata": { + "needs_background": "light" }, - "metadata": {}, "output_type": "display_data" } ], @@ -168,7 +176,58 @@ }, { "cell_type": "code", - "execution_count": 24, + "execution_count": 41, + "outputs": [], + "source": [ + "s = cv2.cvtColor(img, cv2.COLOR_RGB2HSV)[..., 1]\n", + "high_s = s > 175" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 42, + "outputs": [ + { + "data": { + "text/plain": "", + "application/javascript": "/* Put everything inside the global mpl namespace */\n/* global mpl */\nwindow.mpl = {};\n\nmpl.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(\n '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\nmpl.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 = document.createElement('div');\n this.root.setAttribute('style', 'display: inline-block');\n this._root_extra_style(this.root);\n\n parent_element.appendChild(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 (fig.ratio !== 1) {\n fig.send_message('set_device_pixel_ratio', {\n device_pixel_ratio: fig.ratio,\n });\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\nmpl.figure.prototype._init_header = function () {\n var titlebar = document.createElement('div');\n titlebar.classList =\n 'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n var titletext = document.createElement('div');\n titletext.classList = 'ui-dialog-title';\n titletext.setAttribute(\n 'style',\n 'width: 100%; text-align: center; padding: 3px;'\n );\n titlebar.appendChild(titletext);\n this.root.appendChild(titlebar);\n this.header = titletext;\n};\n\nmpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n\nmpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n\nmpl.figure.prototype._init_canvas = function () {\n var fig = this;\n\n var canvas_div = (this.canvas_div = document.createElement('div'));\n canvas_div.setAttribute(\n 'style',\n 'border: 1px solid #ddd;' +\n 'box-sizing: content-box;' +\n 'clear: both;' +\n 'min-height: 1px;' +\n 'min-width: 1px;' +\n 'outline: 0;' +\n 'overflow: hidden;' +\n 'position: relative;' +\n 'resize: both;'\n );\n\n function on_keyboard_event_closure(name) {\n return function (event) {\n return fig.key_event(event, name);\n };\n }\n\n canvas_div.addEventListener(\n 'keydown',\n on_keyboard_event_closure('key_press')\n );\n canvas_div.addEventListener(\n 'keyup',\n on_keyboard_event_closure('key_release')\n );\n\n this._canvas_extra_style(canvas_div);\n this.root.appendChild(canvas_div);\n\n var canvas = (this.canvas = document.createElement('canvas'));\n canvas.classList.add('mpl-canvas');\n canvas.setAttribute('style', 'box-sizing: content-box;');\n\n this.context = canvas.getContext('2d');\n\n var backingStore =\n this.context.backingStorePixelRatio ||\n this.context.webkitBackingStorePixelRatio ||\n this.context.mozBackingStorePixelRatio ||\n this.context.msBackingStorePixelRatio ||\n this.context.oBackingStorePixelRatio ||\n this.context.backingStorePixelRatio ||\n 1;\n\n this.ratio = (window.devicePixelRatio || 1) / backingStore;\n\n var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n 'canvas'\n ));\n rubberband_canvas.setAttribute(\n 'style',\n 'box-sizing: content-box; position: absolute; left: 0; top: 0; z-index: 1;'\n );\n\n // Apply a ponyfill if ResizeObserver is not implemented by browser.\n if (this.ResizeObserver === undefined) {\n if (window.ResizeObserver !== undefined) {\n this.ResizeObserver = window.ResizeObserver;\n } else {\n var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n this.ResizeObserver = obs.ResizeObserver;\n }\n }\n\n this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n var nentries = entries.length;\n for (var i = 0; i < nentries; i++) {\n var entry = entries[i];\n var width, height;\n if (entry.contentBoxSize) {\n if (entry.contentBoxSize instanceof Array) {\n // Chrome 84 implements new version of spec.\n width = entry.contentBoxSize[0].inlineSize;\n height = entry.contentBoxSize[0].blockSize;\n } else {\n // Firefox implements old version of spec.\n width = entry.contentBoxSize.inlineSize;\n height = entry.contentBoxSize.blockSize;\n }\n } else {\n // Chrome <84 implements even older version of spec.\n width = entry.contentRect.width;\n height = entry.contentRect.height;\n }\n\n // Keep the size of the canvas and rubber band canvas in sync with\n // the canvas container.\n if (entry.devicePixelContentBoxSize) {\n // Chrome 84 implements new version of spec.\n canvas.setAttribute(\n 'width',\n entry.devicePixelContentBoxSize[0].inlineSize\n );\n canvas.setAttribute(\n 'height',\n entry.devicePixelContentBoxSize[0].blockSize\n );\n } else {\n canvas.setAttribute('width', width * fig.ratio);\n canvas.setAttribute('height', height * fig.ratio);\n }\n canvas.setAttribute(\n 'style',\n 'width: ' + width + 'px; height: ' + height + 'px;'\n );\n\n rubberband_canvas.setAttribute('width', width);\n rubberband_canvas.setAttribute('height', height);\n\n // And update the size in Python. We ignore the initial 0/0 size\n // that occurs as the element is placed into the DOM, which should\n // otherwise not happen due to the minimum size styling.\n if (fig.ws.readyState == 1 && width != 0 && height != 0) {\n fig.request_resize(width, height);\n }\n }\n });\n this.resizeObserverInstance.observe(canvas_div);\n\n function on_mouse_event_closure(name) {\n return function (event) {\n return fig.mouse_event(event, name);\n };\n }\n\n rubberband_canvas.addEventListener(\n 'mousedown',\n on_mouse_event_closure('button_press')\n );\n rubberband_canvas.addEventListener(\n 'mouseup',\n on_mouse_event_closure('button_release')\n );\n rubberband_canvas.addEventListener(\n 'dblclick',\n on_mouse_event_closure('dblclick')\n );\n // Throttle sequential mouse events to 1 every 20ms.\n rubberband_canvas.addEventListener(\n 'mousemove',\n on_mouse_event_closure('motion_notify')\n );\n\n rubberband_canvas.addEventListener(\n 'mouseenter',\n on_mouse_event_closure('figure_enter')\n );\n rubberband_canvas.addEventListener(\n 'mouseleave',\n on_mouse_event_closure('figure_leave')\n );\n\n canvas_div.addEventListener('wheel', function (event) {\n if (event.deltaY < 0) {\n event.step = 1;\n } else {\n event.step = -1;\n }\n on_mouse_event_closure('scroll')(event);\n });\n\n canvas_div.appendChild(canvas);\n canvas_div.appendChild(rubberband_canvas);\n\n this.rubberband_context = rubberband_canvas.getContext('2d');\n this.rubberband_context.strokeStyle = '#000000';\n\n this._resize_canvas = function (width, height, forward) {\n if (forward) {\n canvas_div.style.width = width + 'px';\n canvas_div.style.height = height + 'px';\n }\n };\n\n // Disable right mouse context menu.\n this.rubberband_canvas.addEventListener('contextmenu', function (_e) {\n event.preventDefault();\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\nmpl.figure.prototype._init_toolbar = function () {\n var fig = this;\n\n var toolbar = document.createElement('div');\n toolbar.classList = 'mpl-toolbar';\n this.root.appendChild(toolbar);\n\n function on_click_closure(name) {\n return function (_event) {\n return fig.toolbar_button_onclick(name);\n };\n }\n\n function on_mouseover_closure(tooltip) {\n return function (event) {\n if (!event.currentTarget.disabled) {\n return fig.toolbar_button_onmouseover(tooltip);\n }\n };\n }\n\n fig.buttons = {};\n var buttonGroup = document.createElement('div');\n buttonGroup.classList = 'mpl-button-group';\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 /* Instead of a spacer, we start a new button group. */\n if (buttonGroup.hasChildNodes()) {\n toolbar.appendChild(buttonGroup);\n }\n buttonGroup = document.createElement('div');\n buttonGroup.classList = 'mpl-button-group';\n continue;\n }\n\n var button = (fig.buttons[name] = document.createElement('button'));\n button.classList = 'mpl-widget';\n button.setAttribute('role', 'button');\n button.setAttribute('aria-disabled', 'false');\n button.addEventListener('click', on_click_closure(method_name));\n button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n\n var icon_img = document.createElement('img');\n icon_img.src = '_images/' + image + '.png';\n icon_img.srcset = '_images/' + image + '_large.png 2x';\n icon_img.alt = tooltip;\n button.appendChild(icon_img);\n\n buttonGroup.appendChild(button);\n }\n\n if (buttonGroup.hasChildNodes()) {\n toolbar.appendChild(buttonGroup);\n }\n\n var fmt_picker = document.createElement('select');\n fmt_picker.classList = 'mpl-widget';\n toolbar.appendChild(fmt_picker);\n this.format_dropdown = fmt_picker;\n\n for (var ind in mpl.extensions) {\n var fmt = mpl.extensions[ind];\n var option = document.createElement('option');\n option.selected = fmt === mpl.default_extension;\n option.innerHTML = fmt;\n fmt_picker.appendChild(option);\n }\n\n var status_bar = document.createElement('span');\n status_bar.classList = 'mpl-message';\n toolbar.appendChild(status_bar);\n this.message = status_bar;\n};\n\nmpl.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\nmpl.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\nmpl.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\nmpl.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\nmpl.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], msg['forward']);\n fig.send_message('refresh', {});\n }\n};\n\nmpl.figure.prototype.handle_rubberband = function (fig, msg) {\n var x0 = msg['x0'] / fig.ratio;\n var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n var x1 = msg['x1'] / fig.ratio;\n var y1 = (fig.canvas.height - msg['y1']) / fig.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,\n 0,\n fig.canvas.width / fig.ratio,\n fig.canvas.height / fig.ratio\n );\n\n fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n};\n\nmpl.figure.prototype.handle_figure_label = function (fig, msg) {\n // Updates the figure title.\n fig.header.textContent = msg['label'];\n};\n\nmpl.figure.prototype.handle_cursor = function (fig, msg) {\n fig.rubberband_canvas.style.cursor = msg['cursor'];\n};\n\nmpl.figure.prototype.handle_message = function (fig, msg) {\n fig.message.textContent = msg['message'];\n};\n\nmpl.figure.prototype.handle_draw = function (fig, _msg) {\n // Request the server to send over a new figure.\n fig.send_draw_message();\n};\n\nmpl.figure.prototype.handle_image_mode = function (fig, msg) {\n fig.image_mode = msg['mode'];\n};\n\nmpl.figure.prototype.handle_history_buttons = function (fig, msg) {\n for (var key in msg) {\n if (!(key in fig.buttons)) {\n continue;\n }\n fig.buttons[key].disabled = !msg[key];\n fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n }\n};\n\nmpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n if (msg['mode'] === 'PAN') {\n fig.buttons['Pan'].classList.add('active');\n fig.buttons['Zoom'].classList.remove('active');\n } else if (msg['mode'] === 'ZOOM') {\n fig.buttons['Pan'].classList.remove('active');\n fig.buttons['Zoom'].classList.add('active');\n } else {\n fig.buttons['Pan'].classList.remove('active');\n fig.buttons['Zoom'].classList.remove('active');\n }\n};\n\nmpl.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.\nmpl.figure.prototype._make_on_message_function = function (fig) {\n return function socket_on_message(evt) {\n if (evt.data instanceof Blob) {\n var img = evt.data;\n if (img.type !== 'image/png') {\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 img.type = 'image/png';\n }\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\n fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n img\n );\n fig.updated_canvas_event();\n fig.waiting = false;\n return;\n } else if (\n typeof evt.data === 'string' &&\n evt.data.slice(0, 21) === 'data:image/png;base64'\n ) {\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(\n \"No handler for the '\" + msg_type + \"' message type: \",\n msg\n );\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(\n \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n e,\n e.stack,\n msg\n );\n }\n }\n };\n};\n\n// from https://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\nmpl.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 }\n if (e.target) {\n targ = e.target;\n } else if (e.srcElement) {\n targ = e.srcElement;\n }\n if (targ.nodeType === 3) {\n // defeat Safari bug\n targ = targ.parentNode;\n }\n\n // pageX,Y are the mouse positions relative to the document\n var boundingRect = targ.getBoundingClientRect();\n var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n var y = e.pageY - (boundingRect.top + document.body.scrollTop);\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 * https://stackoverflow.com/a/24161582/3208463\n */\nfunction simpleKeys(original) {\n return Object.keys(original).reduce(function (obj, key) {\n if (typeof original[key] !== 'object') {\n obj[key] = original[key];\n }\n return obj;\n }, {});\n}\n\nmpl.figure.prototype.mouse_event = function (event, name) {\n var canvas_pos = mpl.findpos(event);\n\n if (name === 'button_press') {\n this.canvas.focus();\n this.canvas_div.focus();\n }\n\n var x = canvas_pos.x * this.ratio;\n var y = canvas_pos.y * this.ratio;\n\n this.send_message(name, {\n x: x,\n y: y,\n button: event.button,\n step: event.step,\n guiEvent: simpleKeys(event),\n });\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\nmpl.figure.prototype._key_event_extra = function (_event, _name) {\n // Handle any extra behaviour associated with a key event\n};\n\nmpl.figure.prototype.key_event = function (event, name) {\n // Prevent repeat events\n if (name === 'key_press') {\n if (event.key === this._key) {\n return;\n } else {\n this._key = event.key;\n }\n }\n if (name === 'key_release') {\n this._key = null;\n }\n\n var value = '';\n if (event.ctrlKey && event.key !== 'Control') {\n value += 'ctrl+';\n }\n else if (event.altKey && event.key !== 'Alt') {\n value += 'alt+';\n }\n else if (event.shiftKey && event.key !== 'Shift') {\n value += 'shift+';\n }\n\n value += 'k' + event.key;\n\n this._key_event_extra(event, name);\n\n this.send_message(name, { key: value, guiEvent: simpleKeys(event) });\n return false;\n};\n\nmpl.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\nmpl.figure.prototype.toolbar_button_onmouseover = function (tooltip) {\n this.message.textContent = tooltip;\n};\n\n///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n// prettier-ignore\nvar _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\nmpl.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\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n\nmpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n\nmpl.default_extension = \"png\";/* global mpl */\n\nvar 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.binaryType = comm.kernel.ws.binaryType;\n ws.readyState = comm.kernel.ws.readyState;\n function updateReadyState(_event) {\n if (comm.kernel.ws) {\n ws.readyState = comm.kernel.ws.readyState;\n } else {\n ws.readyState = 3; // Closed state.\n }\n }\n comm.kernel.ws.addEventListener('open', updateReadyState);\n comm.kernel.ws.addEventListener('close', updateReadyState);\n comm.kernel.ws.addEventListener('error', updateReadyState);\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 var data = msg['content']['data'];\n if (data['blob'] !== undefined) {\n data = {\n data: new Blob(msg['buffers'], { type: data['blob'] }),\n };\n }\n // Pass the mpl event to the overridden (by mpl) onmessage function.\n ws.onmessage(data);\n });\n return ws;\n};\n\nmpl.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 = document.getElementById(id);\n var ws_proxy = comm_websocket_adapter(comm);\n\n function ondownload(figure, _format) {\n window.open(figure.canvas.toDataURL());\n }\n\n var fig = new mpl.figure(id, ws_proxy, ondownload, element);\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;\n fig.cell_info = mpl.find_output_cell(\"
\");\n if (!fig.cell_info) {\n console.error('Failed to find cell for figure', id, fig);\n return;\n }\n fig.cell_info[0].output_area.element.on(\n 'cleared',\n { fig: fig },\n fig._remove_fig_handler\n );\n};\n\nmpl.figure.prototype.handle_close = function (fig, msg) {\n var width = fig.canvas.width / fig.ratio;\n fig.cell_info[0].output_area.element.off(\n 'cleared',\n fig._remove_fig_handler\n );\n fig.resizeObserverInstance.unobserve(fig.canvas_div);\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.innerHTML =\n '';\n fig.close_ws(fig, msg);\n};\n\nmpl.figure.prototype.close_ws = function (fig, msg) {\n fig.send_message('closing', msg);\n // fig.ws.close()\n};\n\nmpl.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 / this.ratio;\n var dataURL = this.canvas.toDataURL();\n this.cell_info[1]['text/html'] =\n '';\n};\n\nmpl.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 () {\n fig.push_to_output();\n }, 1000);\n};\n\nmpl.figure.prototype._init_toolbar = function () {\n var fig = this;\n\n var toolbar = document.createElement('div');\n toolbar.classList = 'btn-toolbar';\n this.root.appendChild(toolbar);\n\n function on_click_closure(name) {\n return function (_event) {\n return fig.toolbar_button_onclick(name);\n };\n }\n\n function on_mouseover_closure(tooltip) {\n return function (event) {\n if (!event.currentTarget.disabled) {\n return fig.toolbar_button_onmouseover(tooltip);\n }\n };\n }\n\n fig.buttons = {};\n var buttonGroup = document.createElement('div');\n buttonGroup.classList = 'btn-group';\n var button;\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 /* Instead of a spacer, we start a new button group. */\n if (buttonGroup.hasChildNodes()) {\n toolbar.appendChild(buttonGroup);\n }\n buttonGroup = document.createElement('div');\n buttonGroup.classList = 'btn-group';\n continue;\n }\n\n button = fig.buttons[name] = document.createElement('button');\n button.classList = 'btn btn-default';\n button.href = '#';\n button.title = name;\n button.innerHTML = '';\n button.addEventListener('click', on_click_closure(method_name));\n button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n buttonGroup.appendChild(button);\n }\n\n if (buttonGroup.hasChildNodes()) {\n toolbar.appendChild(buttonGroup);\n }\n\n // Add the status bar.\n var status_bar = document.createElement('span');\n status_bar.classList = 'mpl-message pull-right';\n toolbar.appendChild(status_bar);\n this.message = status_bar;\n\n // Add the close button to the window.\n var buttongrp = document.createElement('div');\n buttongrp.classList = 'btn-group inline pull-right';\n button = document.createElement('button');\n button.classList = 'btn btn-mini btn-primary';\n button.href = '#';\n button.title = 'Stop Interaction';\n button.innerHTML = '';\n button.addEventListener('click', function (_evt) {\n fig.handle_close(fig, {});\n });\n button.addEventListener(\n 'mouseover',\n on_mouseover_closure('Stop Interaction')\n );\n buttongrp.appendChild(button);\n var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n titlebar.insertBefore(buttongrp, titlebar.firstChild);\n};\n\nmpl.figure.prototype._remove_fig_handler = function (event) {\n var fig = event.data.fig;\n if (event.target !== this) {\n // Ignore bubbled events from children.\n return;\n }\n fig.close_ws(fig, {});\n};\n\nmpl.figure.prototype._root_extra_style = function (el) {\n el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n};\n\nmpl.figure.prototype._canvas_extra_style = function (el) {\n // this is important to make the div 'focusable\n el.setAttribute('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 } else {\n // location in version 2\n IPython.keyboard_manager.register_events(el);\n }\n};\n\nmpl.figure.prototype._key_event_extra = function (event, _name) {\n // Check for shift+enter\n if (event.shiftKey && event.which === 13) {\n this.canvas_div.blur();\n // select the cell after this one\n var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n IPython.notebook.select(index + 1);\n }\n};\n\nmpl.figure.prototype.handle_save = function (fig, _msg) {\n fig.ondownload(fig, null);\n};\n\nmpl.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.\nif (IPython.notebook.kernel !== null) {\n IPython.notebook.kernel.comm_manager.register_target(\n 'matplotlib',\n mpl.mpl_figure_comm\n );\n}\n" + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "text/plain": "", + "text/html": "
" + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "%matplotlib notebook\n", + "fig,axs = plt.subplots(3, 1)\n", + "axs[0].imshow(img)\n", + "axs[1].imshow(s)\n", + "axs[2].imshow(high_s)\n", + "plt.savefig(\"2.png\",dpi=900)" + ], + "metadata": { + "collapsed": false, + "pycharm": { + "name": "#%%\n" + } + } + }, + { + "cell_type": "code", + "execution_count": 8, "outputs": [], "source": [], "metadata": { diff --git a/config.py b/config.py index 9702068..9c0b574 100644 --- a/config.py +++ b/config.py @@ -4,22 +4,26 @@ import numpy as np class Config: # 文件相关参数 - nRows, nCols, nBands, threshold, = 256, 1024, 22, 3 - nrgbRows, nrgbCols, nrgbBands, rgb_threshold = 1024, 4096, 3, 2 + nRows, nCols, nBands, spec_size_threshold, = 256, 1024, 22, 3 + nRgbRows, nRgbCols, nRgbBands, rgb_size_threshold = 1024, 4096, 3, 4 + # 需要设置的谱段等参数 selected_bands = [127, 201, 202, 294] - bands = [127, 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, - 294] + bands = [127, 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, + 211, 212, 213, 214, 215, 216, 217, 218, 219, 294] is_yellow_min = np.array([0.10167048, 0.1644719, 0.1598884, 0.31534621]) is_yellow_max = np.array([0.212984, 0.25896924, 0.26509268, 0.51943593]) is_black_threshold = np.asarray([0.1369, 0.1472, 0.1439, 0.1814]) black_yellow_bands = [0, 2, 3, 21] green_bands = [i for i in range(1, 21)] - # 机器学习模型相关参数 + # 光谱模型参数 blk_size = 4 pixel_model_path = r"./models/dt.p" blk_model_path = r"./models/rf_4x4_c22_20_sen8_8.model" - rgb_model_path = r"./models/dt_2022-07-20_14-40.model" - rgb_tobacco_model_path = r"models/beijing_dt_2022-07-21_16-44.model" - rgb_background_model_path = r"models/tobacco_dt_2022-07-21_16-30.model" + + # rgb模型参数 + rgb_tobacco_model_path = r"models/tobacco_dt_2022-07-21_16-30.model" + rgb_background_model_path = r"models/beijing_dt_2022-07-21_16-44.model" + threshold_low, threshold_high = 5, 255 + threshold_s = 175 diff --git a/main.py b/main.py index 0a1cfe6..30d152e 100644 --- a/main.py +++ b/main.py @@ -1,20 +1,20 @@ import os import time import numpy as np +import scipy.io + from config import Config -from models import ManualTree, AnonymousColorDetector +from models import RgbDetector, SpecDetector, ManualTree, AnonymousColorDetector +import cv2 -# 主函数 def main(): - threshold = Config.threshold - rgb_threshold = Config.rgb_threshold - manual_tree = ManualTree(blk_model_path=Config.blk_model_path, pixel_model_path=Config.pixel_model_path) - tobacco_detector = AnonymousColorDetector(file_path=Config.rgb_tobacco_model_path) - background_detector = AnonymousColorDetector(file_path=Config.rgb_background_model_path) - + spec_detector = SpecDetector(blk_model_path=Config.blk_model_path, + pixel_model_path=Config.pixel_model_path) + rgb_detector = RgbDetector(tobacco_model_path=Config.rgb_tobacco_model_path, + background_model_path=Config.rgb_background_model_path) total_len = Config.nRows * Config.nCols * Config.nBands * 4 # float型变量, 4个字节 - total_rgb = Config.nrgbRows * Config.nrgbCols * Config.nrgbBands * 1 # int型变量 + total_rgb = Config.nRgbRows * Config.nRgbCols * Config.nRgbBands * 1 # int型变量 if not os.access(img_fifo_path, os.F_OK): os.mkfifo(img_fifo_path, 0o777) if not os.access(mask_fifo_path, os.F_OK): @@ -25,6 +25,67 @@ def main(): fd_img = os.open(img_fifo_path, os.O_RDONLY) fd_rgb = os.open(rgb_fifo_path, os.O_RDONLY) data = os.read(fd_img, total_len) + # 读取(开启一个管道) + if len(data) < 3: + threshold = int(float(data)) + print("[INFO] Get threshold: ", threshold) + continue + else: + data_total = data + rgb_data = os.read(fd_rgb, total_rgb) + if len(rgb_data) < 3: + rgb_threshold = int(float(rgb_data)) + print(rgb_threshold) + continue + else: + rgb_data_total = rgb_data + os.close(fd_img) + os.close(fd_rgb) + # 识别 + t1 = time.time() + img_data = np.frombuffer(data_total, dtype=np.float32).reshape((Config.nRows, Config.nBands, -1)).transpose(0, + 2, + 1) + rgb_data = np.frombuffer(rgb_data_total, dtype=np.uint8).reshape((Config.nRgbRows, Config.nRgbCols, -1)) + # 光谱识别 + mask = spec_detector.predict(img_data) + # rgb识别 + mask_rgb = rgb_detector.predict(rgb_data) + # 结果合并 + mask_result = (mask | mask_rgb).astype(np.uint8) + mask_result = mask_result.repeat(Config.blk_size, axis=0).repeat(Config.blk_size, axis=1).astype(np.uint8) + t2 = time.time() + print(f'rgb len = {len(rgb_data)}') + + # 写出 + fd_mask = os.open(mask_fifo_path, os.O_WRONLY) + os.write(fd_mask, mask_result.tobytes()) + os.close(fd_mask) + t3 = time.time() + print(f'total time is:{t3 - t1}') + + +def save_main(): + threshold = Config.spec_size_threshold + rgb_threshold = Config.rgb_size_threshold + manual_tree = ManualTree(blk_model_path=Config.blk_model_path, pixel_model_path=Config.pixel_model_path) + tobacco_detector = AnonymousColorDetector(file_path=Config.rgb_tobacco_model_path) + background_detector = AnonymousColorDetector(file_path=Config.rgb_background_model_path) + total_len = Config.nRows * Config.nCols * Config.nBands * 4 # float型变量, 4个字节 + total_rgb = Config.nRgbRows * Config.nRgbCols * Config.nRgbBands * 1 # int型变量 + if not os.access(img_fifo_path, os.F_OK): + os.mkfifo(img_fifo_path, 0o777) + if not os.access(mask_fifo_path, os.F_OK): + os.mkfifo(mask_fifo_path, 0o777) + if not os.access(rgb_fifo_path, os.F_OK): + os.mkfifo(rgb_fifo_path, 0o777) + img_list = [] + idx = 0 + while idx <= 30: + idx += 1 + fd_img = os.open(img_fifo_path, os.O_RDONLY) + fd_rgb = os.open(rgb_fifo_path, os.O_RDONLY) + data = os.read(fd_img, total_len) # 读取(开启一个管道) if len(data) < 3: @@ -40,36 +101,40 @@ def main(): continue else: rgb_data_total = rgb_data - os.close(fd_img) os.close(fd_rgb) # 识别 t1 = time.time() - img_data = np.frombuffer(data_total, dtype=np.float32).reshape((Config.nRows, Config.nBands, - -1)).transpose(0, 2, 1) - rgb_data = np.frombuffer(rgb_data_total, dtype=np.uint8) + img_data = np.frombuffer(data_total, dtype=np.float32).reshape((Config.nRows, Config.nBands, -1)). \ + transpose(0, 2, 1) + rgb_data = np.frombuffer(rgb_data_total, dtype=np.uint8).reshape((Config.nRgbRows, Config.nRgbCols, -1)) + img_list.append((rgb_data.copy(), img_data.copy())) + pixel_predict_result = manual_tree.pixel_predict_ml_dilation(data=img_data, iteration=1) blk_predict_result = manual_tree.blk_predict(data=img_data) rgb_data = tobacco_detector.pretreatment(rgb_data) - rgb_predict_result = 1 - ( - background_detector.predict(rgb_data) | tobacco_detector.swell(tobacco_detector.predict(rgb_data))) + # print(rgb_data.shape) + rgb_predict_result = 1 - (background_detector.predict(rgb_data, threshold_low=Config.threshold_low, + threshold_high=Config.threshold_high) | + tobacco_detector.swell(tobacco_detector.predict(rgb_data, + threshold_low=Config.threshold_low, + threshold_high=Config.threshold_high))) mask_rgb = rgb_predict_result.reshape(Config.nRows, Config.nCols // Config.blk_size, Config.blk_size) \ .sum(axis=2).reshape(Config.nRows // 4, Config.blk_size, Config.nCols // Config.blk_size) \ .sum(axis=1) mask_rgb[mask_rgb <= rgb_threshold] = 0 mask_rgb[mask_rgb > rgb_threshold] = 1 - mask = (pixel_predict_result & blk_predict_result).astype(np.uint8) mask = mask.reshape(Config.nRows, Config.nCols // Config.blk_size, Config.blk_size) \ .sum(axis=2).reshape(Config.nRows // 4, Config.blk_size, Config.nCols // Config.blk_size) \ .sum(axis=1) mask[mask <= threshold] = 0 mask[mask > threshold] = 1 - mask_result = (mask | mask_rgb).astype(np.uint8) + # mask_result = (mask | mask_rgb).astype(np.uint8) + mask_result = mask_rgb mask_result = mask_result.repeat(Config.blk_size, axis=0).repeat(Config.blk_size, axis=1).astype(np.uint8) t2 = time.time() - # print(f'shibie time is:{t2 - t1}') print(f'rgb len = {len(rgb_data)}') # 写出 @@ -78,6 +143,14 @@ def main(): os.close(fd_mask) t3 = time.time() print(f'total time is:{t3 - t1}') + i = 0 + print("Stop Serving") + for img in img_list: + print(f"writing img {i}...") + cv2.imwrite(f"./{i}.png", img[0][..., ::-1]) + np.save(f'./{i}.npy', img[1]) + i += 1 + print("save success") if __name__ == '__main__': diff --git a/main_test.py b/main_test.py index 880f8a6..2a329b0 100644 --- a/main_test.py +++ b/main_test.py @@ -10,20 +10,27 @@ import cv2 import matplotlib.pyplot as plt import numpy as np -from models import Detector, AnonymousColorDetector -from utils import read_labeled_img +from config import Config +from models import Detector, AnonymousColorDetector, ManualTree +from utils import read_labeled_img, size_threshold -def virtual_main(detector: AnonymousColorDetector, test_img=None, test_img_dir=None, test_model=False): +def pony_run(test_img=None, test_img_dir=None, test_spectra=False, test_rgb=False): """ 虚拟读图测试程序 - :param detector: 杂质探测器，需要继承Detector类 :param test_img: 测试图像，rgb格式的图片或者路径 :param test_img_dir: 测试图像文件夹 - :param test_model: 是否进行模型约束性测试 + :param test_spectra: 是否测试光谱 + :param test_rgb: 是否测试rgb :return: """ + if (test_img is not None) or (test_img_dir is not None): + threshold = Config.spec_size_threshold + rgb_threshold = Config.rgb_size_threshold + manual_tree = ManualTree(blk_model_path=Config.blk_model_path, pixel_model_path=Config.pixel_model_path) + tobacco_detector = AnonymousColorDetector(file_path=Config.rgb_tobacco_model_path) + background_detector = AnonymousColorDetector(file_path=Config.rgb_background_model_path) if test_img is not None: if isinstance(test_img, str): img = cv2.imread(test_img)[:, :, ::-1] @@ -31,42 +38,45 @@ def virtual_main(detector: AnonymousColorDetector, test_img=None, test_img_dir=N img = test_img else: raise TypeError("test img should be np.ndarray or str") - t1 = time.time() - img = cv2.resize(img, (1024, 256)) - t2 = time.time() - result = 1 - detector.predict(img) - t3 = time.time() - fig, axs = plt.subplots(3, 1) - axs[0].imshow(img) - axs[1].imshow(result) - mask_color = np.zeros_like(img) - mask_color[result > 0] = (0, 0, 255) - result_show = cv2.addWeighted(img, 1, mask_color, 0.5, 0) - axs[2].imshow(result_show) - axs[0].set_title( - f' resize {(t2 - t1) * 1000:.2f} ms, predict {(t3 - t2) * 1000:.2f} ms, total {(t3 - t1) * 1000:.2f} ms') - plt.show() if test_img_dir is not None: - image_names = os.listdir(test_img_dir) + image_names = [img_name for img_name in os.listdir(test_img_dir) if img_name.endswith('.png')] for image_name in image_names: - img = cv2.imread(os.path.join(test_img_dir, image_name))[..., ::-1] - if test_model: - data_dir = "data/dataset" - color_dict = {(0, 0, 255): "yangeng"} - dataset = read_labeled_img(data_dir, color_dict=color_dict, is_ps_color_space=False) - ground_truth = dataset['yangeng'] - world_boundary = np.array([0, 0, 0, 255, 255, 255]) - detector.visualize(world_boundary, sample_size=50000, class_max_num=5000, ground_truth=ground_truth) + rgb_data = cv2.imread(os.path.join(test_img_dir, image_name))[..., ::-1] + # 识别 + t1 = time.time() + if test_spectra: + # spectra part + pixel_predict_result = manual_tree.pixel_predict_ml_dilation(data=img_data, iteration=1) + blk_predict_result = manual_tree.blk_predict(data=img_data) + mask = (pixel_predict_result & blk_predict_result).astype(np.uint8) + mask_spec = size_threshold(mask, Config.blk_size, threshold) + if test_rgb: + # rgb part + rgb_data = tobacco_detector.pretreatment(rgb_data) + background = background_detector.predict(rgb_data) + tobacco = tobacco_detector.predict(rgb_data) + tobacco_d = tobacco_detector.swell(tobacco) + rgb_predict_result = 1 - (background | tobacco_d) + mask_rgb = size_threshold(rgb_predict_result, Config.blk_size, Config.rgb_size_threshold) + fig, axs = plt.subplots(5, 1, figsize=(12, 10), constrained_layout=True) + axs[0].imshow(rgb_data) + axs[0].set_title("rgb raw data") + axs[1].imshow(background) + axs[1].set_title("background") + axs[2].imshow(tobacco) + axs[2].set_title("tobacco") + axs[3].imshow(rgb_predict_result) + axs[3].set_title("1 - (background + dilate(tobacco))") + axs[4].imshow(mask_rgb) + axs[4].set_title("final mask") + plt.show() + + mask_result = (mask | mask_rgb).astype(np.uint8) + # mask_result = rgb_predict_result + mask_result = mask_result.repeat(Config.blk_size, axis=0).repeat(Config.blk_size, axis=1).astype(np.uint8) + t2 = time.time() + print(f'rgb len = {len(rgb_data)}') if __name__ == '__main__': - model = AnonymousColorDetector(file_path='dt_2022-07-20_14-40.model') - virtual_main(model, - test_img=r'C:\Users\FEIJINTI\Desktop\720\binning1\tobacco\Image_2022_0720_1354_46_472-003051.bmp', - test_model=True) - virtual_main(model, - test_img=r'C:\Users\FEIJINTI\Desktop\720\binning1\tobacco\Image_2022_0720_1354_46_472-003051.bmp', - test_model=True) - virtual_main(model, - test_img=r'C:\Users\FEIJINTI\Desktop\720\binning1\tobacco\Image_2022_0720_1354_46_472-003051.bmp', - test_model=True) + pony_run(test_img_dir=r'E:\zhouchao\725data', test_rgb=True) diff --git a/models.py b/models.py index bb4eae9..5716f15 100644 --- a/models.py +++ b/models.py @@ -9,17 +9,19 @@ import pickle import cv2 import numpy as np import scipy.io -import tqdm from scipy.ndimage import binary_dilation from sklearn.tree import DecisionTreeClassifier from sklearn.metrics import classification_report from sklearn.model_selection import train_test_split from config import Config -from utils import lab_scatter, read_labeled_img -from tqdm import tqdm +from utils import lab_scatter, read_labeled_img, size_threshold -from elm import ELM +deploy = False +if not deploy: + print("Training env") + from tqdm import tqdm + from elm import ELM class Detector(object): @@ -87,7 +89,7 @@ class AnonymousColorDetector(Detector): y_predict = self.model.predict(x_val) print(classification_report(y_true=y_val, y_pred=y_predict)) - def predict(self, x, threshold_low=10, threshold_high=170): + def predict(self, x, threshold_low=5, threshold_high=255): """ 输入rgb彩色图像 @@ -98,9 +100,11 @@ class AnonymousColorDetector(Detector): x = cv2.cvtColor(x, cv2.COLOR_RGB2LAB) x = x.reshape(w * h, -1) mask = (threshold_low < x[:, 0]) & (x[:, 0] < threshold_high) - mask_result = self.model.predict(x[mask]) - result = np.ones((w * h,)) - result[mask] = mask_result + result = np.ones((w * h,), dtype=np.uint8) + + if np.any(mask): + mask_result = self.model.predict(x[mask]) + result[mask] = mask_result return result.reshape(h, w) @staticmethod @@ -133,6 +137,12 @@ class AnonymousColorDetector(Detector): bar.close() return negative_samples + def pretreatment(self, x): + return cv2.resize(x, (1024, 256)) + + def swell(self, x): + return cv2.dilate(x, kernel=np.ones((3, 3), np.uint8)) + def save(self): path = datetime.datetime.now().strftime(f"{self.model_type}_%Y-%m-%d_%H-%M.model") with open(path, 'wb') as f: @@ -292,6 +302,137 @@ class BlkModel: self.rfc = pickle.load(f) +class RgbDetector(Detector): + def __init__(self, tobacco_model_path, background_model_path): + self.background_detector = None + self.tobacco_detector = None + self.load(tobacco_model_path, background_model_path) + + def predict(self, rgb_data): + rgb_data = self.tobacco_detector.pretreatment(rgb_data) # resize to the required size + background = self.background_detector.predict(rgb_data) + tobacco = self.tobacco_detector.predict(rgb_data) + tobacco_d = self.tobacco_detector.swell(tobacco) # dilate the tobacco to remove the tobacco edge error + high_s = cv2.cvtColor(rgb_data, cv2.COLOR_RGB2HSV)[..., 1] > Config.threshold_s + non_tobacco_or_background = 1 - (background | tobacco_d) # 既非烟梗也非背景的区域 + rgb_predict_result = high_s | non_tobacco_or_background # 高饱和度区域或者是双非区域都是杂质 + mask_rgb = size_threshold(rgb_predict_result, Config.blk_size, Config.rgb_size_threshold) # 杂质大小限制，超过大小的才打 + return mask_rgb + + def load(self, tobacco_model_path, background_model_path): + self.tobacco_detector = AnonymousColorDetector(tobacco_model_path) + self.background_detector = AnonymousColorDetector(background_model_path) + + def save(self, *args, **kwargs): + pass + + def fit(self, *args, **kwargs): + pass + + +class SpecDetector(Detector): + # 初始化机器学习像素模型、深度学习像素模型、分块模型 + def __init__(self, blk_model_path, pixel_model_path): + self.blk_model = None + self.pixel_model_ml = None + self.load(blk_model_path, pixel_model_path) + + def load(self, blk_model_path, pixel_model_path): + self.pixel_model_ml = PixelModelML(pixel_model_path) + self.blk_model = BlkModel(blk_model_path) + + def predict(self, img_data): + pixel_predict_result = self.pixel_predict_ml_dilation(data=img_data, iteration=1) + blk_predict_result = self.blk_predict(data=img_data) + mask = (pixel_predict_result & blk_predict_result).astype(np.uint8) + mask = size_threshold(mask, Config.blk_size, Config.spec_size_threshold) + return mask + + def save(self, *args, **kwargs): + pass + + def fit(self, *args, **kwargs): + pass + + # 区分烟梗和非黄色且非背景的杂质 + @staticmethod + def is_yellow(features): + features = features.reshape((Config.nRows * Config.nCols), len(Config.selected_bands)) + sum_x = features.sum(axis=1)[..., np.newaxis] + rate = features / sum_x + mask = ((rate < Config.is_yellow_max) & (rate > Config.is_yellow_min)) + mask = np.all(mask, axis=1).reshape(Config.nRows, Config.nCols) + return mask + + # 区分背景和黄色杂质 + @staticmethod + def is_black(feature, threshold): + feature = feature.reshape((Config.nRows * Config.nCols), feature.shape[2]) + mask = (feature <= threshold) + mask = np.all(mask, axis=1).reshape(Config.nRows, Config.nCols) + return mask + + # 预测出烟梗的mask + def predict_tobacco(self, x: np.ndarray) -> np.ndarray: + """ + 预测出烟梗的mask + :param x: 图像数据，形状是 nRows x nCols x nBands + :return: bool类型的mask，是否为烟梗, True为烟梗 + """ + black_res = self.is_black(x[..., Config.black_yellow_bands], Config.is_black_threshold) + yellow_res = self.is_yellow(x[..., Config.black_yellow_bands]) + yellow_things = (~black_res) & yellow_res + x_yellow = x[yellow_things, ...] + tobacco = self.pixel_model_ml.predict(x_yellow[..., Config.green_bands]) + yellow_things[yellow_things] = tobacco + return yellow_things + + # 预测出杂质的机器学习像素模型 + def pixel_predict_ml_dilation(self, data, iteration) -> np.ndarray: + """ + 预测出杂质的位置mask + :param data: 图像数据，形状是 nRows x nCols x nBands + :param iteration: 膨胀的次数 + :return: bool类型的mask，是否为杂质, True为杂质 + """ + black_res = self.is_black(data[..., Config.black_yellow_bands], Config.is_black_threshold) + yellow_res = self.is_yellow(data[..., Config.black_yellow_bands]) + # non_yellow_things为异色杂质 + non_yellow_things = (~black_res) & (~yellow_res) + # yellow_things为黄色物体(烟梗+杂质) + yellow_things = (~black_res) & yellow_res + # x_yellow为挑出的黄色物体 + x_yellow = data[yellow_things, ...] + if x_yellow.shape[0] == 0: + return non_yellow_things + else: + tobacco = self.pixel_model_ml.predict(x_yellow[..., Config.green_bands]) > 0.5 + + non_yellow_things[yellow_things] = ~tobacco + # 杂质mask中将背景赋值为0,将杂质赋值为1 + non_yellow_things = non_yellow_things + 0 + + # 烟梗mask中将背景赋值为0,将烟梗赋值为2 + yellow_things[yellow_things] = tobacco + yellow_things = yellow_things + 0 + yellow_things = binary_dilation(yellow_things, iterations=iteration) + yellow_things = yellow_things + 0 + yellow_things[yellow_things == 1] = 2 + + # 将杂质mask和烟梗mask相加,得到的mask中含有0(背景),1(杂质),2(烟梗),3(膨胀后的烟梗与杂质相加的部分) + mask = non_yellow_things + yellow_things + mask[mask == 0] = False + mask[mask == 1] = True + mask[mask == 2] = False + mask[mask == 3] = False + return mask + + # 预测出杂质的分块模型 + def blk_predict(self, data): + blk_result_array = self.blk_model.predict(data) + return blk_result_array + + if __name__ == '__main__': data_dir = "data/dataset" color_dict = {(0, 0, 255): "yangeng"} diff --git a/utils.py b/utils.py index a8249f0..25607e1 100644 --- a/utils.py +++ b/utils.py @@ -96,6 +96,14 @@ def lab_scatter(dataset: dict, class_max_num=None, is_3d=False, is_ps_color_spac plt.show() +def size_threshold(img, blk_size, threshold): + mask = img.reshape(img.shape[0], img.shape[1] // blk_size, blk_size).sum(axis=2). \ + reshape(img.shape[0] // blk_size, blk_size, img.shape[1] // blk_size).sum(axis=1) + mask[mask <= threshold] = 0 + mask[mask > threshold] = 1 + return mask + + if __name__ == '__main__': color_dict = {(0, 0, 255): "yangeng", (255, 0, 0): "bejing", (0, 255, 0): "hongdianxian", (255, 0, 255): "chengsebangbangtang", (0, 255, 255): "lvdianxian"}