You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
1067 lines
48 KiB
1067 lines
48 KiB
{ |
|
"cells": [ |
|
{ |
|
"cell_type": "code", |
|
"execution_count": 1, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"# import matplotlib\n", |
|
"%matplotlib notebook\n", |
|
"import matplotlib.pyplot as plt\n", |
|
"from matplotlib import cm\n", |
|
"import numpy as np\n", |
|
"import h5py as h5py\n", |
|
"import cartopy.crs as ccrs\n", |
|
"import cartopy.feature as cfeature\n", |
|
"from cartopy.io.shapereader import Reader\n", |
|
"from cartopy.feature import ShapelyFeature\n", |
|
"# matplotlib.rcParams['backend'] = \"TkAgg\"\n", |
|
"# plt.get_backend()\n", |
|
"# plt.get_backend()" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": 4, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"hdf5 = '/notebooks/3B-MO.MS.MRG.3IMERG.20190101-S000000-E235959.01.V06A.HDF5'\n", |
|
"dataset = h5py.File(hdf5,'r')" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": 5, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"precip = dataset['Grid/precipitation'][:]\n", |
|
"\n", |
|
"precip = np.transpose(precip[0])\n", |
|
"\n", |
|
"theLats = dataset['Grid/lat'][:]\n", |
|
"theLons = dataset['Grid/lon'][:]" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": 2, |
|
"metadata": {}, |
|
"outputs": [ |
|
{ |
|
"data": { |
|
"application/javascript": [ |
|
"/* Put everything inside the global mpl namespace */\n", |
|
"window.mpl = {};\n", |
|
"\n", |
|
"\n", |
|
"mpl.get_websocket_type = function() {\n", |
|
" if (typeof(WebSocket) !== 'undefined') {\n", |
|
" return WebSocket;\n", |
|
" } else if (typeof(MozWebSocket) !== 'undefined') {\n", |
|
" return MozWebSocket;\n", |
|
" } else {\n", |
|
" alert('Your browser does not have WebSocket support.' +\n", |
|
" 'Please try Chrome, Safari or Firefox ≥ 6. ' +\n", |
|
" 'Firefox 4 and 5 are also supported but you ' +\n", |
|
" 'have to enable WebSockets in about:config.');\n", |
|
" };\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n", |
|
" this.id = figure_id;\n", |
|
"\n", |
|
" this.ws = websocket;\n", |
|
"\n", |
|
" this.supports_binary = (this.ws.binaryType != undefined);\n", |
|
"\n", |
|
" if (!this.supports_binary) {\n", |
|
" var warnings = document.getElementById(\"mpl-warnings\");\n", |
|
" if (warnings) {\n", |
|
" warnings.style.display = 'block';\n", |
|
" warnings.textContent = (\n", |
|
" \"This browser does not support binary websocket messages. \" +\n", |
|
" \"Performance may be slow.\");\n", |
|
" }\n", |
|
" }\n", |
|
"\n", |
|
" this.imageObj = new Image();\n", |
|
"\n", |
|
" this.context = undefined;\n", |
|
" this.message = undefined;\n", |
|
" this.canvas = undefined;\n", |
|
" this.rubberband_canvas = undefined;\n", |
|
" this.rubberband_context = undefined;\n", |
|
" this.format_dropdown = undefined;\n", |
|
"\n", |
|
" this.image_mode = 'full';\n", |
|
"\n", |
|
" this.root = $('<div/>');\n", |
|
" this._root_extra_style(this.root)\n", |
|
" this.root.attr('style', 'display: inline-block');\n", |
|
"\n", |
|
" $(parent_element).append(this.root);\n", |
|
"\n", |
|
" this._init_header(this);\n", |
|
" this._init_canvas(this);\n", |
|
" this._init_toolbar(this);\n", |
|
"\n", |
|
" var fig = this;\n", |
|
"\n", |
|
" this.waiting = false;\n", |
|
"\n", |
|
" this.ws.onopen = function () {\n", |
|
" fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n", |
|
" fig.send_message(\"send_image_mode\", {});\n", |
|
" if (mpl.ratio != 1) {\n", |
|
" fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n", |
|
" }\n", |
|
" fig.send_message(\"refresh\", {});\n", |
|
" }\n", |
|
"\n", |
|
" this.imageObj.onload = function() {\n", |
|
" if (fig.image_mode == 'full') {\n", |
|
" // Full images could contain transparency (where diff images\n", |
|
" // almost always do), so we need to clear the canvas so that\n", |
|
" // there is no ghosting.\n", |
|
" fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n", |
|
" }\n", |
|
" fig.context.drawImage(fig.imageObj, 0, 0);\n", |
|
" };\n", |
|
"\n", |
|
" this.imageObj.onunload = function() {\n", |
|
" fig.ws.close();\n", |
|
" }\n", |
|
"\n", |
|
" this.ws.onmessage = this._make_on_message_function(this);\n", |
|
"\n", |
|
" this.ondownload = ondownload;\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype._init_header = function() {\n", |
|
" var titlebar = $(\n", |
|
" '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n", |
|
" 'ui-helper-clearfix\"/>');\n", |
|
" var titletext = $(\n", |
|
" '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n", |
|
" 'text-align: center; padding: 3px;\"/>');\n", |
|
" titlebar.append(titletext)\n", |
|
" this.root.append(titlebar);\n", |
|
" this.header = titletext[0];\n", |
|
"}\n", |
|
"\n", |
|
"\n", |
|
"\n", |
|
"mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n", |
|
"\n", |
|
"}\n", |
|
"\n", |
|
"\n", |
|
"mpl.figure.prototype._root_extra_style = function(canvas_div) {\n", |
|
"\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype._init_canvas = function() {\n", |
|
" var fig = this;\n", |
|
"\n", |
|
" var canvas_div = $('<div/>');\n", |
|
"\n", |
|
" canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n", |
|
"\n", |
|
" function canvas_keyboard_event(event) {\n", |
|
" return fig.key_event(event, event['data']);\n", |
|
" }\n", |
|
"\n", |
|
" canvas_div.keydown('key_press', canvas_keyboard_event);\n", |
|
" canvas_div.keyup('key_release', canvas_keyboard_event);\n", |
|
" this.canvas_div = canvas_div\n", |
|
" this._canvas_extra_style(canvas_div)\n", |
|
" this.root.append(canvas_div);\n", |
|
"\n", |
|
" var canvas = $('<canvas/>');\n", |
|
" canvas.addClass('mpl-canvas');\n", |
|
" canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n", |
|
"\n", |
|
" this.canvas = canvas[0];\n", |
|
" this.context = canvas[0].getContext(\"2d\");\n", |
|
"\n", |
|
" var backingStore = this.context.backingStorePixelRatio ||\n", |
|
"\tthis.context.webkitBackingStorePixelRatio ||\n", |
|
"\tthis.context.mozBackingStorePixelRatio ||\n", |
|
"\tthis.context.msBackingStorePixelRatio ||\n", |
|
"\tthis.context.oBackingStorePixelRatio ||\n", |
|
"\tthis.context.backingStorePixelRatio || 1;\n", |
|
"\n", |
|
" mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n", |
|
"\n", |
|
" var rubberband = $('<canvas/>');\n", |
|
" rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n", |
|
"\n", |
|
" var pass_mouse_events = true;\n", |
|
"\n", |
|
" canvas_div.resizable({\n", |
|
" start: function(event, ui) {\n", |
|
" pass_mouse_events = false;\n", |
|
" },\n", |
|
" resize: function(event, ui) {\n", |
|
" fig.request_resize(ui.size.width, ui.size.height);\n", |
|
" },\n", |
|
" stop: function(event, ui) {\n", |
|
" pass_mouse_events = true;\n", |
|
" fig.request_resize(ui.size.width, ui.size.height);\n", |
|
" },\n", |
|
" });\n", |
|
"\n", |
|
" function mouse_event_fn(event) {\n", |
|
" if (pass_mouse_events)\n", |
|
" return fig.mouse_event(event, event['data']);\n", |
|
" }\n", |
|
"\n", |
|
" rubberband.mousedown('button_press', mouse_event_fn);\n", |
|
" rubberband.mouseup('button_release', mouse_event_fn);\n", |
|
" // Throttle sequential mouse events to 1 every 20ms.\n", |
|
" rubberband.mousemove('motion_notify', mouse_event_fn);\n", |
|
"\n", |
|
" rubberband.mouseenter('figure_enter', mouse_event_fn);\n", |
|
" rubberband.mouseleave('figure_leave', mouse_event_fn);\n", |
|
"\n", |
|
" canvas_div.on(\"wheel\", function (event) {\n", |
|
" event = event.originalEvent;\n", |
|
" event['data'] = 'scroll'\n", |
|
" if (event.deltaY < 0) {\n", |
|
" event.step = 1;\n", |
|
" } else {\n", |
|
" event.step = -1;\n", |
|
" }\n", |
|
" mouse_event_fn(event);\n", |
|
" });\n", |
|
"\n", |
|
" canvas_div.append(canvas);\n", |
|
" canvas_div.append(rubberband);\n", |
|
"\n", |
|
" this.rubberband = rubberband;\n", |
|
" this.rubberband_canvas = rubberband[0];\n", |
|
" this.rubberband_context = rubberband[0].getContext(\"2d\");\n", |
|
" this.rubberband_context.strokeStyle = \"#000000\";\n", |
|
"\n", |
|
" this._resize_canvas = function(width, height) {\n", |
|
" // Keep the size of the canvas, canvas container, and rubber band\n", |
|
" // canvas in synch.\n", |
|
" canvas_div.css('width', width)\n", |
|
" canvas_div.css('height', height)\n", |
|
"\n", |
|
" canvas.attr('width', width * mpl.ratio);\n", |
|
" canvas.attr('height', height * mpl.ratio);\n", |
|
" canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n", |
|
"\n", |
|
" rubberband.attr('width', width);\n", |
|
" rubberband.attr('height', height);\n", |
|
" }\n", |
|
"\n", |
|
" // Set the figure to an initial 600x600px, this will subsequently be updated\n", |
|
" // upon first draw.\n", |
|
" this._resize_canvas(600, 600);\n", |
|
"\n", |
|
" // Disable right mouse context menu.\n", |
|
" $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n", |
|
" return false;\n", |
|
" });\n", |
|
"\n", |
|
" function set_focus () {\n", |
|
" canvas.focus();\n", |
|
" canvas_div.focus();\n", |
|
" }\n", |
|
"\n", |
|
" window.setTimeout(set_focus, 100);\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype._init_toolbar = function() {\n", |
|
" var fig = this;\n", |
|
"\n", |
|
" var nav_element = $('<div/>')\n", |
|
" nav_element.attr('style', 'width: 100%');\n", |
|
" this.root.append(nav_element);\n", |
|
"\n", |
|
" // Define a callback function for later on.\n", |
|
" function toolbar_event(event) {\n", |
|
" return fig.toolbar_button_onclick(event['data']);\n", |
|
" }\n", |
|
" function toolbar_mouse_event(event) {\n", |
|
" return fig.toolbar_button_onmouseover(event['data']);\n", |
|
" }\n", |
|
"\n", |
|
" for(var toolbar_ind in mpl.toolbar_items) {\n", |
|
" var name = mpl.toolbar_items[toolbar_ind][0];\n", |
|
" var tooltip = mpl.toolbar_items[toolbar_ind][1];\n", |
|
" var image = mpl.toolbar_items[toolbar_ind][2];\n", |
|
" var method_name = mpl.toolbar_items[toolbar_ind][3];\n", |
|
"\n", |
|
" if (!name) {\n", |
|
" // put a spacer in here.\n", |
|
" continue;\n", |
|
" }\n", |
|
" var button = $('<button/>');\n", |
|
" button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n", |
|
" 'ui-button-icon-only');\n", |
|
" button.attr('role', 'button');\n", |
|
" button.attr('aria-disabled', 'false');\n", |
|
" button.click(method_name, toolbar_event);\n", |
|
" button.mouseover(tooltip, toolbar_mouse_event);\n", |
|
"\n", |
|
" var icon_img = $('<span/>');\n", |
|
" icon_img.addClass('ui-button-icon-primary ui-icon');\n", |
|
" icon_img.addClass(image);\n", |
|
" icon_img.addClass('ui-corner-all');\n", |
|
"\n", |
|
" var tooltip_span = $('<span/>');\n", |
|
" tooltip_span.addClass('ui-button-text');\n", |
|
" tooltip_span.html(tooltip);\n", |
|
"\n", |
|
" button.append(icon_img);\n", |
|
" button.append(tooltip_span);\n", |
|
"\n", |
|
" nav_element.append(button);\n", |
|
" }\n", |
|
"\n", |
|
" var fmt_picker_span = $('<span/>');\n", |
|
"\n", |
|
" var fmt_picker = $('<select/>');\n", |
|
" fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n", |
|
" fmt_picker_span.append(fmt_picker);\n", |
|
" nav_element.append(fmt_picker_span);\n", |
|
" this.format_dropdown = fmt_picker[0];\n", |
|
"\n", |
|
" for (var ind in mpl.extensions) {\n", |
|
" var fmt = mpl.extensions[ind];\n", |
|
" var option = $(\n", |
|
" '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n", |
|
" fmt_picker.append(option)\n", |
|
" }\n", |
|
"\n", |
|
" // Add hover states to the ui-buttons\n", |
|
" $( \".ui-button\" ).hover(\n", |
|
" function() { $(this).addClass(\"ui-state-hover\");},\n", |
|
" function() { $(this).removeClass(\"ui-state-hover\");}\n", |
|
" );\n", |
|
"\n", |
|
" var status_bar = $('<span class=\"mpl-message\"/>');\n", |
|
" nav_element.append(status_bar);\n", |
|
" this.message = status_bar[0];\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n", |
|
" // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n", |
|
" // which will in turn request a refresh of the image.\n", |
|
" this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.send_message = function(type, properties) {\n", |
|
" properties['type'] = type;\n", |
|
" properties['figure_id'] = this.id;\n", |
|
" this.ws.send(JSON.stringify(properties));\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.send_draw_message = function() {\n", |
|
" if (!this.waiting) {\n", |
|
" this.waiting = true;\n", |
|
" this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n", |
|
" }\n", |
|
"}\n", |
|
"\n", |
|
"\n", |
|
"mpl.figure.prototype.handle_save = function(fig, msg) {\n", |
|
" var format_dropdown = fig.format_dropdown;\n", |
|
" var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n", |
|
" fig.ondownload(fig, format);\n", |
|
"}\n", |
|
"\n", |
|
"\n", |
|
"mpl.figure.prototype.handle_resize = function(fig, msg) {\n", |
|
" var size = msg['size'];\n", |
|
" if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n", |
|
" fig._resize_canvas(size[0], size[1]);\n", |
|
" fig.send_message(\"refresh\", {});\n", |
|
" };\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n", |
|
" var x0 = msg['x0'] / mpl.ratio;\n", |
|
" var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n", |
|
" var x1 = msg['x1'] / mpl.ratio;\n", |
|
" var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n", |
|
" x0 = Math.floor(x0) + 0.5;\n", |
|
" y0 = Math.floor(y0) + 0.5;\n", |
|
" x1 = Math.floor(x1) + 0.5;\n", |
|
" y1 = Math.floor(y1) + 0.5;\n", |
|
" var min_x = Math.min(x0, x1);\n", |
|
" var min_y = Math.min(y0, y1);\n", |
|
" var width = Math.abs(x1 - x0);\n", |
|
" var height = Math.abs(y1 - y0);\n", |
|
"\n", |
|
" fig.rubberband_context.clearRect(\n", |
|
" 0, 0, fig.canvas.width, fig.canvas.height);\n", |
|
"\n", |
|
" fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n", |
|
" // Updates the figure title.\n", |
|
" fig.header.textContent = msg['label'];\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.handle_cursor = function(fig, msg) {\n", |
|
" var cursor = msg['cursor'];\n", |
|
" switch(cursor)\n", |
|
" {\n", |
|
" case 0:\n", |
|
" cursor = 'pointer';\n", |
|
" break;\n", |
|
" case 1:\n", |
|
" cursor = 'default';\n", |
|
" break;\n", |
|
" case 2:\n", |
|
" cursor = 'crosshair';\n", |
|
" break;\n", |
|
" case 3:\n", |
|
" cursor = 'move';\n", |
|
" break;\n", |
|
" }\n", |
|
" fig.rubberband_canvas.style.cursor = cursor;\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.handle_message = function(fig, msg) {\n", |
|
" fig.message.textContent = msg['message'];\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.handle_draw = function(fig, msg) {\n", |
|
" // Request the server to send over a new figure.\n", |
|
" fig.send_draw_message();\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n", |
|
" fig.image_mode = msg['mode'];\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.updated_canvas_event = function() {\n", |
|
" // Called whenever the canvas gets updated.\n", |
|
" this.send_message(\"ack\", {});\n", |
|
"}\n", |
|
"\n", |
|
"// A function to construct a web socket function for onmessage handling.\n", |
|
"// Called in the figure constructor.\n", |
|
"mpl.figure.prototype._make_on_message_function = function(fig) {\n", |
|
" return function socket_on_message(evt) {\n", |
|
" if (evt.data instanceof Blob) {\n", |
|
" /* FIXME: We get \"Resource interpreted as Image but\n", |
|
" * transferred with MIME type text/plain:\" errors on\n", |
|
" * Chrome. But how to set the MIME type? It doesn't seem\n", |
|
" * to be part of the websocket stream */\n", |
|
" evt.data.type = \"image/png\";\n", |
|
"\n", |
|
" /* Free the memory for the previous frames */\n", |
|
" if (fig.imageObj.src) {\n", |
|
" (window.URL || window.webkitURL).revokeObjectURL(\n", |
|
" fig.imageObj.src);\n", |
|
" }\n", |
|
"\n", |
|
" fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n", |
|
" evt.data);\n", |
|
" fig.updated_canvas_event();\n", |
|
" fig.waiting = false;\n", |
|
" return;\n", |
|
" }\n", |
|
" else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n", |
|
" fig.imageObj.src = evt.data;\n", |
|
" fig.updated_canvas_event();\n", |
|
" fig.waiting = false;\n", |
|
" return;\n", |
|
" }\n", |
|
"\n", |
|
" var msg = JSON.parse(evt.data);\n", |
|
" var msg_type = msg['type'];\n", |
|
"\n", |
|
" // Call the \"handle_{type}\" callback, which takes\n", |
|
" // the figure and JSON message as its only arguments.\n", |
|
" try {\n", |
|
" var callback = fig[\"handle_\" + msg_type];\n", |
|
" } catch (e) {\n", |
|
" console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n", |
|
" return;\n", |
|
" }\n", |
|
"\n", |
|
" if (callback) {\n", |
|
" try {\n", |
|
" // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n", |
|
" callback(fig, msg);\n", |
|
" } catch (e) {\n", |
|
" console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n", |
|
" }\n", |
|
" }\n", |
|
" };\n", |
|
"}\n", |
|
"\n", |
|
"// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n", |
|
"mpl.findpos = function(e) {\n", |
|
" //this section is from http://www.quirksmode.org/js/events_properties.html\n", |
|
" var targ;\n", |
|
" if (!e)\n", |
|
" e = window.event;\n", |
|
" if (e.target)\n", |
|
" targ = e.target;\n", |
|
" else if (e.srcElement)\n", |
|
" targ = e.srcElement;\n", |
|
" if (targ.nodeType == 3) // defeat Safari bug\n", |
|
" targ = targ.parentNode;\n", |
|
"\n", |
|
" // jQuery normalizes the pageX and pageY\n", |
|
" // pageX,Y are the mouse positions relative to the document\n", |
|
" // offset() returns the position of the element relative to the document\n", |
|
" var x = e.pageX - $(targ).offset().left;\n", |
|
" var y = e.pageY - $(targ).offset().top;\n", |
|
"\n", |
|
" return {\"x\": x, \"y\": y};\n", |
|
"};\n", |
|
"\n", |
|
"/*\n", |
|
" * return a copy of an object with only non-object keys\n", |
|
" * we need this to avoid circular references\n", |
|
" * http://stackoverflow.com/a/24161582/3208463\n", |
|
" */\n", |
|
"function simpleKeys (original) {\n", |
|
" return Object.keys(original).reduce(function (obj, key) {\n", |
|
" if (typeof original[key] !== 'object')\n", |
|
" obj[key] = original[key]\n", |
|
" return obj;\n", |
|
" }, {});\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.mouse_event = function(event, name) {\n", |
|
" var canvas_pos = mpl.findpos(event)\n", |
|
"\n", |
|
" if (name === 'button_press')\n", |
|
" {\n", |
|
" this.canvas.focus();\n", |
|
" this.canvas_div.focus();\n", |
|
" }\n", |
|
"\n", |
|
" var x = canvas_pos.x * mpl.ratio;\n", |
|
" var y = canvas_pos.y * mpl.ratio;\n", |
|
"\n", |
|
" this.send_message(name, {x: x, y: y, button: event.button,\n", |
|
" step: event.step,\n", |
|
" guiEvent: simpleKeys(event)});\n", |
|
"\n", |
|
" /* This prevents the web browser from automatically changing to\n", |
|
" * the text insertion cursor when the button is pressed. We want\n", |
|
" * to control all of the cursor setting manually through the\n", |
|
" * 'cursor' event from matplotlib */\n", |
|
" event.preventDefault();\n", |
|
" return false;\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype._key_event_extra = function(event, name) {\n", |
|
" // Handle any extra behaviour associated with a key event\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.key_event = function(event, name) {\n", |
|
"\n", |
|
" // Prevent repeat events\n", |
|
" if (name == 'key_press')\n", |
|
" {\n", |
|
" if (event.which === this._key)\n", |
|
" return;\n", |
|
" else\n", |
|
" this._key = event.which;\n", |
|
" }\n", |
|
" if (name == 'key_release')\n", |
|
" this._key = null;\n", |
|
"\n", |
|
" var value = '';\n", |
|
" if (event.ctrlKey && event.which != 17)\n", |
|
" value += \"ctrl+\";\n", |
|
" if (event.altKey && event.which != 18)\n", |
|
" value += \"alt+\";\n", |
|
" if (event.shiftKey && event.which != 16)\n", |
|
" value += \"shift+\";\n", |
|
"\n", |
|
" value += 'k';\n", |
|
" value += event.which.toString();\n", |
|
"\n", |
|
" this._key_event_extra(event, name);\n", |
|
"\n", |
|
" this.send_message(name, {key: value,\n", |
|
" guiEvent: simpleKeys(event)});\n", |
|
" return false;\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.toolbar_button_onclick = function(name) {\n", |
|
" if (name == 'download') {\n", |
|
" this.handle_save(this, null);\n", |
|
" } else {\n", |
|
" this.send_message(\"toolbar_button\", {name: name});\n", |
|
" }\n", |
|
"};\n", |
|
"\n", |
|
"mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n", |
|
" this.message.textContent = tooltip;\n", |
|
"};\n", |
|
"mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n", |
|
"\n", |
|
"mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n", |
|
"\n", |
|
"mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n", |
|
" // Create a \"websocket\"-like object which calls the given IPython comm\n", |
|
" // object with the appropriate methods. Currently this is a non binary\n", |
|
" // socket, so there is still some room for performance tuning.\n", |
|
" var ws = {};\n", |
|
"\n", |
|
" ws.close = function() {\n", |
|
" comm.close()\n", |
|
" };\n", |
|
" ws.send = function(m) {\n", |
|
" //console.log('sending', m);\n", |
|
" comm.send(m);\n", |
|
" };\n", |
|
" // Register the callback with on_msg.\n", |
|
" comm.on_msg(function(msg) {\n", |
|
" //console.log('receiving', msg['content']['data'], msg);\n", |
|
" // Pass the mpl event to the overridden (by mpl) onmessage function.\n", |
|
" ws.onmessage(msg['content']['data'])\n", |
|
" });\n", |
|
" return ws;\n", |
|
"}\n", |
|
"\n", |
|
"mpl.mpl_figure_comm = function(comm, msg) {\n", |
|
" // This is the function which gets called when the mpl process\n", |
|
" // starts-up an IPython Comm through the \"matplotlib\" channel.\n", |
|
"\n", |
|
" var id = msg.content.data.id;\n", |
|
" // Get hold of the div created by the display call when the Comm\n", |
|
" // socket was opened in Python.\n", |
|
" var element = $(\"#\" + id);\n", |
|
" var ws_proxy = comm_websocket_adapter(comm)\n", |
|
"\n", |
|
" function ondownload(figure, format) {\n", |
|
" window.open(figure.imageObj.src);\n", |
|
" }\n", |
|
"\n", |
|
" var fig = new mpl.figure(id, ws_proxy,\n", |
|
" ondownload,\n", |
|
" element.get(0));\n", |
|
"\n", |
|
" // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n", |
|
" // web socket which is closed, not our websocket->open comm proxy.\n", |
|
" ws_proxy.onopen();\n", |
|
"\n", |
|
" fig.parent_element = element.get(0);\n", |
|
" fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n", |
|
" if (!fig.cell_info) {\n", |
|
" console.error(\"Failed to find cell for figure\", id, fig);\n", |
|
" return;\n", |
|
" }\n", |
|
"\n", |
|
" var output_index = fig.cell_info[2]\n", |
|
" var cell = fig.cell_info[0];\n", |
|
"\n", |
|
"};\n", |
|
"\n", |
|
"mpl.figure.prototype.handle_close = function(fig, msg) {\n", |
|
" var width = fig.canvas.width/mpl.ratio\n", |
|
" fig.root.unbind('remove')\n", |
|
"\n", |
|
" // Update the output cell to use the data from the current canvas.\n", |
|
" fig.push_to_output();\n", |
|
" var dataURL = fig.canvas.toDataURL();\n", |
|
" // Re-enable the keyboard manager in IPython - without this line, in FF,\n", |
|
" // the notebook keyboard shortcuts fail.\n", |
|
" IPython.keyboard_manager.enable()\n", |
|
" $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n", |
|
" fig.close_ws(fig, msg);\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.close_ws = function(fig, msg){\n", |
|
" fig.send_message('closing', msg);\n", |
|
" // fig.ws.close()\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.push_to_output = function(remove_interactive) {\n", |
|
" // Turn the data on the canvas into data in the output cell.\n", |
|
" var width = this.canvas.width/mpl.ratio\n", |
|
" var dataURL = this.canvas.toDataURL();\n", |
|
" this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.updated_canvas_event = function() {\n", |
|
" // Tell IPython that the notebook contents must change.\n", |
|
" IPython.notebook.set_dirty(true);\n", |
|
" this.send_message(\"ack\", {});\n", |
|
" var fig = this;\n", |
|
" // Wait a second, then push the new image to the DOM so\n", |
|
" // that it is saved nicely (might be nice to debounce this).\n", |
|
" setTimeout(function () { fig.push_to_output() }, 1000);\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype._init_toolbar = function() {\n", |
|
" var fig = this;\n", |
|
"\n", |
|
" var nav_element = $('<div/>')\n", |
|
" nav_element.attr('style', 'width: 100%');\n", |
|
" this.root.append(nav_element);\n", |
|
"\n", |
|
" // Define a callback function for later on.\n", |
|
" function toolbar_event(event) {\n", |
|
" return fig.toolbar_button_onclick(event['data']);\n", |
|
" }\n", |
|
" function toolbar_mouse_event(event) {\n", |
|
" return fig.toolbar_button_onmouseover(event['data']);\n", |
|
" }\n", |
|
"\n", |
|
" for(var toolbar_ind in mpl.toolbar_items){\n", |
|
" var name = mpl.toolbar_items[toolbar_ind][0];\n", |
|
" var tooltip = mpl.toolbar_items[toolbar_ind][1];\n", |
|
" var image = mpl.toolbar_items[toolbar_ind][2];\n", |
|
" var method_name = mpl.toolbar_items[toolbar_ind][3];\n", |
|
"\n", |
|
" if (!name) { continue; };\n", |
|
"\n", |
|
" var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n", |
|
" button.click(method_name, toolbar_event);\n", |
|
" button.mouseover(tooltip, toolbar_mouse_event);\n", |
|
" nav_element.append(button);\n", |
|
" }\n", |
|
"\n", |
|
" // Add the status bar.\n", |
|
" var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n", |
|
" nav_element.append(status_bar);\n", |
|
" this.message = status_bar[0];\n", |
|
"\n", |
|
" // Add the close button to the window.\n", |
|
" var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n", |
|
" var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n", |
|
" button.click(function (evt) { fig.handle_close(fig, {}); } );\n", |
|
" button.mouseover('Stop Interaction', toolbar_mouse_event);\n", |
|
" buttongrp.append(button);\n", |
|
" var titlebar = this.root.find($('.ui-dialog-titlebar'));\n", |
|
" titlebar.prepend(buttongrp);\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype._root_extra_style = function(el){\n", |
|
" var fig = this\n", |
|
" el.on(\"remove\", function(){\n", |
|
"\tfig.close_ws(fig, {});\n", |
|
" });\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype._canvas_extra_style = function(el){\n", |
|
" // this is important to make the div 'focusable\n", |
|
" el.attr('tabindex', 0)\n", |
|
" // reach out to IPython and tell the keyboard manager to turn it's self\n", |
|
" // off when our div gets focus\n", |
|
"\n", |
|
" // location in version 3\n", |
|
" if (IPython.notebook.keyboard_manager) {\n", |
|
" IPython.notebook.keyboard_manager.register_events(el);\n", |
|
" }\n", |
|
" else {\n", |
|
" // location in version 2\n", |
|
" IPython.keyboard_manager.register_events(el);\n", |
|
" }\n", |
|
"\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype._key_event_extra = function(event, name) {\n", |
|
" var manager = IPython.notebook.keyboard_manager;\n", |
|
" if (!manager)\n", |
|
" manager = IPython.keyboard_manager;\n", |
|
"\n", |
|
" // Check for shift+enter\n", |
|
" if (event.shiftKey && event.which == 13) {\n", |
|
" this.canvas_div.blur();\n", |
|
" event.shiftKey = false;\n", |
|
" // Send a \"J\" for go to next cell\n", |
|
" event.which = 74;\n", |
|
" event.keyCode = 74;\n", |
|
" manager.command_mode();\n", |
|
" manager.handle_keydown(event);\n", |
|
" }\n", |
|
"}\n", |
|
"\n", |
|
"mpl.figure.prototype.handle_save = function(fig, msg) {\n", |
|
" fig.ondownload(fig, null);\n", |
|
"}\n", |
|
"\n", |
|
"\n", |
|
"mpl.find_output_cell = function(html_output) {\n", |
|
" // Return the cell and output element which can be found *uniquely* in the notebook.\n", |
|
" // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n", |
|
" // IPython event is triggered only after the cells have been serialised, which for\n", |
|
" // our purposes (turning an active figure into a static one), is too late.\n", |
|
" var cells = IPython.notebook.get_cells();\n", |
|
" var ncells = cells.length;\n", |
|
" for (var i=0; i<ncells; i++) {\n", |
|
" var cell = cells[i];\n", |
|
" if (cell.cell_type === 'code'){\n", |
|
" for (var j=0; j<cell.output_area.outputs.length; j++) {\n", |
|
" var data = cell.output_area.outputs[j];\n", |
|
" if (data.data) {\n", |
|
" // IPython >= 3 moved mimebundle to data attribute of output\n", |
|
" data = data.data;\n", |
|
" }\n", |
|
" if (data['text/html'] == html_output) {\n", |
|
" return [cell, data, j];\n", |
|
" }\n", |
|
" }\n", |
|
" }\n", |
|
" }\n", |
|
"}\n", |
|
"\n", |
|
"// Register the function which deals with the matplotlib target/channel.\n", |
|
"// The kernel may be null if the page has been refreshed.\n", |
|
"if (IPython.notebook.kernel != null) {\n", |
|
" IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n", |
|
"}\n" |
|
], |
|
"text/plain": [ |
|
"<IPython.core.display.Javascript object>" |
|
] |
|
}, |
|
"metadata": {}, |
|
"output_type": "display_data" |
|
}, |
|
{ |
|
"data": { |
|
"text/html": [ |
|
"<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfQAAAEsCAYAAAA1u0HIAAAXiklEQVR4nO3d63KUZfq28aM7Yb0WWaooooggUVmDLCSKIAgKMozKi4gif1BoURl1BpNyB9kIvrgb5/sBEoOickvSd/fTx6/q/DJWyaVWcczTpLuJJEnqe9Q+QJIkPT6DLklSAxh0SZIawKBLktQABl2SpAYw6JIkNYBBlySpAQy6JEkNYNAlSWoAgy5JUgMYdEmSGsCgS5LUAAZdkqQGMOiSJDWAQZckqQEMuiRJDWDQJUlqAIMuSVIDGHRJkhrAoEuS1AAGXZKkBjDokiQ1gEGXJKkBDLokSQ1g0CVJagCDLklSAxh0SZIawKBLktQABl2SpAYw6JIkNYBBlySpAQy6JEkNYNAlSWoAgy5JUgMYdBW7e/du7ty545zrsd29e7f2bw+qyKCryN27dzPMrADOuR7b/PnzjfoAM+gqcufOnQA5xpl8whXnpm17OTwZplWszRkuVLljIYsD5C1OVv93UrJjnAmQO3fu1P5tQpUYdBWZCPonXMlNfnFuWvcZN7KSNQHSopVNjOQ6t7v264+wK0C28Hr1fxel+4QrBn3AGXQVMeiuGzvLpSxmaYC0GcpODqTD+Iz+mqOcCJDFLJ3xX2smZtBl0FXEoLtu7iinM4e5AfIET6bD2Iz8Ooc5HiCzmZNr/FT9n/ufzKDLoKuIQXfdXofxvMxIgCxl+V8+Pb/P+WxgU9axIbs4mIvc+Nu//1a2T8b8Et9U/+f9pzPoMugqYtBdrW1hW4Ds4I0H/vdr/JQRdk0+yf9+bdpZyhPZxEje5lTO8Xk+4Hy2sS/zmB8gS1jWt0/mEzPoMugqYtBdzc1mToaZlQ7jeZ/zkz9AB2QOczPCrlzjp3QYz1kuZYSdWc7KDDP80NgPMZTX2N2Xf2b++xl0GXQVmQj6TvbnLU7mCKfyDh/kGGfyLh/mOP/Ke5zLST7K+3ySDzifM1zIWS7mLJdyjs/zEZfzMVdynqu5wLV8ytf5jBu5xDf5gm/zJd/nCv/JVX7MV/w317mdDmMz9ptuh/EpG5vcdX6estsP7Cv+N2X/ndw1fprcVX6csh8e2BX+M2W3Jvcl30/uMt/lC77N59zMJb7JZ9zIp3ydC3yd81zLea7mY67kY77MR1zOOT7POT7PWT7LWS7mDBdymgv5gPN5n09yio9zko/yHudynH/lOGdzjA9zjDN5h9M5wgc5wqm8xcmM8l5GOZHDHM+bvJtDHMsB3skBjmQ/b+cN3so+RrOXw9nDm9nNoeziYHZxIDvZn+28ke3syzb2ZRt78jq78yq7M8KujLAjW9meV9ieLWzLFl7PZl7Ly7yalxnJJrbmJbZmI6/kRTbnRTbnBV7OBjZlAYsmfwJ+IsorWZP3Of+3/52vcCvv8mH2cjgHeCcfN+z/lBp0GXQVmQi6c7U3ws5c5cfqIe2VGXQZdBWZCPoB3pl80nt38ilv6hPevae7QxzLQY5OPtn99lR374luJweyY/Jpbu/9J7ldGWHn/ae4qU9wE09vr2QjW7KRVyb30gPbOrlND2xkci8/sFcnt5nXHrotvD5l2yb3ypRtZftDN8KOKds5Zbsm9+qUvcbuyb3O7mxjb7axL9t5Izt4Izs5kF0czG4OZg+Hsoc3s5fR7GM0+3k7+zmSA7yTgxzNIY7lTd7NYY5nlBMZ5b28PeWVlaOcnvLqytkHXmE5xceTr7Kc5kLO8OkDr7Sc44v7r7Z8mU/uv+Ly/7iWC3ydi1zPxfuvulzi5pRXXm7l/+6/+nLt/qsa17mdG4/4CkyH8QB5npeqB7TXZtBl0FXEP0N3NXeJToBsY0/1W3ptBl0GXUUMuqu5E5wLkLd4r/otvTaDLoOuIgbd1dwe3gyQf3O5+i29NoMug64iBt3V3Kb7HzDzFf+rfkuvzaDLoKuIQXc19xTPpkWr+h29OIMug64iBt3V3CKWZJhZ1e/oxRl0GXQVMeiu1v7N5QBZx4bqt/TiDLoMuooYdFdry1kZIJf5rvotvTiDLoOuIgbd1djE0/mzPp3/6Qy6DLqKGHRXYz6d//0Mugy6ihh01+35dP5oM+gy6Cpi0F2359P5o82gy6CriEF33ZxP548+gy6DriIG3XVzPp0/+gy6DLqKGHTXrfl0XjaDLoOuIgbddWs+nZfNoMugq4hBd92YT+flM+gy6Cpi0F03tpI1Pp0XzqDLoKuIQXczvQ5jadHKclZWv6WfZtBl0FXEoLuZ3lucDJBDHKt+Sz/NoMugq4hBdzO91TwdINf5ufot/TSDLoOuIgbdzfTmMi/zWVj9jn6bQZdBVxGD7mZ6LVpZw9PV7+i3GXQZdBUx6G6m16KdVaytfke/zaDLoKuIQXczvfkszFzmVb+j32bQZdBVxKC7md5a1qVFq/od/TaDLoOuIgbdzfReZ3eAXORG9Vv6aQZdBl1FDLqb6Z3i4wA5yNHqt/TTDLoMuooYdDfTu87tANnApuq39NMMugy6ihh0140NM5xlLK9+Rz/NoMugq4hBd93YEpZlFrOq39FPM+gy6Cpi0F03tp4XA6TDWPVb+mUGXQZdRQy668b2MhogZ/i0+i39MoMug64iBt11Y+f5vwDZwRvVb+mXGXQZdBUx6K4b6zAeIM+wvvot/TKDLoOuIgbddWtzmJuFLKp+R7/MoMugq4hBd93ak6xKm3b1O/plBl0GXUUMuuvWXmYkQK7yQ/Vb+mEGXQZdRQy669be5lSAHONM9Vv6YQZdBl1FDLrr1r7k+wDZwuvVb+mHGXQZdBUx6K6ba9POStZUv6MfZtBl0FXEoLtubj4LM5d51e/ohxl0GXQVMeium5vHgsxnQfU7+mEGXQZdRQy66+aGmeW3rj3iDLoMuooYdNetfcV/A2Q9L1a/pR9m0GXQVcSgu25tlBMBcoRT1W/phxl0GXQVMeiuW1vB6gC5zu3qt/TDDLoMuooYdNeNXeFWgKzmqeq39MsMugy6ihh01429zKsBcpbPqt/SLzPoMugqYtDdTK/DeIYZzjzmV7+ln2bQZdBVxKC7md6R+5/hvodD1W/ppxl0GXQVMehupreEZWnRznV+rn5LP82gy6CriEF3M7mL3AiQ53ih+i39NoMug64iBt3N5NbzYoBc4Ovqt/TbDLoMuooYdDdT6zCWNu0sZmn1W/pxBl0GXUUMupup7eVwgIzyXvVb+nEGXQZdRQy6m6nNZ0GGGU6H8eq39OMMugy6ihh0NxM7y6UA2cRI9Vv6dQZdBl1FDLqbiU18bvsVblW/pV9n0GXQVcSgu+neaS4EyDOsr35LP8+gy6CriEF3070FLEyLVq7yQ/Vb+nkGXQZdRQy6m84d4J0AGWFn9Vv6fQZdBl1FDLqbrl3n5wwznFnM8ifbp2EGXQZdRQy6m65t5jXfdz6NM+gy6Cpi0N107Eu+D7SyiCXVb2nKDLoMuooYdDcdW87KADnH59VvacoMugy6ihh097h7lg0B8hJbq9/SpBl0GXQVMejucbaRLQHyFOuq39K0GXQZdBUx6O6f7iQfBciTrKp+SxNn0GXQVcSgu3+6ecxPm3au8VP1W5o4gy6DriIG3f2THeX0/Q+Q2VH9lqbOoMugq4hBd/9kC1iYNu1c53b1W5o6gy6DriIG3ZXuOP8KkM28Vv2WJs+gy6CriEF3pVvI4rRo5yv+V/2WJs+gy6CriEF3JXuXD33PeZdm0GXQVcSgu0fdZ9xIm3aGGPYn27swgy6DriIG3T3KLnIjs5kdIGe5VP2eQZhBl0FXEYPu/m4fcD5t2gHyFier3zMoM+gy6Cpi0N1f7SBHA6TNUM5ysfo9gzSDLoOuIgbd/dkmvt98DnNziW+q3zNoM+gy6Cpi0N3v12Esa3g6QJax3LenVZpBl0FXEYPupu4K/8lCFgXIejamw3j1mwZ1Bl0GXUUMupvYea5mmFkBsoM3qt8z6DPoMugqYtDdTX7JFW6lRTstWjnK6er3OIMug65CBt3d5Jc8z0sBcpoL1W9x92bQZdBVxKC7m/ySIYayiCXV73C/zaDLoKuIQXdf8K3fntaDM+gy6Cpi0N1+jgTIKT6ufov7bQZdBl1FDLp7hvUB0mGs+i3utxl0GXQVMehuHvMzl3nV73APzqDLoKuIQR/sdRgPkKd4tvot7sEZdBl0FTHog73TXAiQvRyufot7cAZdBl1FDPpgbyvbA/jlKz04gy6DriIGfbC3nBVpM1T9DvfHGXQZdBUx6IO9YYazhCeq3+H+OIMug64iBn1wd4VbAbKJrdVvcX+cQZdBVxGDPrg7xLEAOc6/qt/i/jiDLoOuIgZ9cPcsLwTIdX6ufov74wy6DLqKGPTB3QIWZg5zq9/hHj6DLoOuIgZ9MHfvA2VaWcMz1W9xD59Bl0FXEYM+mDvLpQDZxcHqt7iHz6DLoKuIQR/MjbArQC7wdfVb3MNn0GXQVcSgD+ZWsDpt2tXvcH8+gy6DriIGfTC3iCWZzZzqd7g/n0GXQVcRgz6YW87KALnCreq3uIfPoMugq4hBH8yd4/MA2ciW6re4h8+gy6CriEEf3C1kcdq002Gs+i3ujzPoMugqYtAHd29zKkB2sr/6Le6PM+gy6Cpi0Ad7s5ntD8f16Ay6DLqKGPTB3nb2BcgoJ6rf4h6cQZdBVxGDPti7zs8ZZlaGGMo1fqp+j/ttBl0GXUUMujvBuQB5ktXVb3G/zaDLoKuIQXc3+SUvsjlAXmRz9VvcvRl0GXQVMehuYitY7Re29NAMugy6ihh0N7EOY1nAogB5hw+q3zPoM+gy6Cpi0N3UXeWHzGJ2oJWzXKp+zyDPoMugq4hBd7/fBb5Om3baDOUS31S/Z1Bn0GXQVcSgu4ftA84HyGzm+Ha2SjPoMugqYtDdn+0wxwNkMUvTYbz6PYM2gy6DriIG3f3VXmN3gKzlmeq3DNoMugy6ihh093d7jhfvf/DMqlzndvV7BmUGXQZdRQy6e5RtZMvkn6kf4pgvwXdhBl0GXUUMunvUHeRohhgKkGFmZRv7cp2fq9/V1Bl0GXQVMeiuZB3Gs4c3M4e5AdKmnad5Lgc5mqv8WP2+Js2gy6CriEF3/3RHOJVFLAkwudnMybNsyL+5XP2+fp9Bl0FXEYPuHncdxnKcs9nApixg4WTc5zIv7/NJ9fv6dQZdBl1FDLqb7l3mu2zh9bRpB8gIO6vf1I8z6DLoKmLQ3UztS77PYpYGyG4OVb+n32bQZdBVxKC7mVyHscxjfqCVL/i2+j39NIMug64iBt3N9M7xeYCsZE31W/ppBl0GXUUMuuvGnmVDgDzPxnQYq35PP8ygy6CriEF33ViH8axkzeR719eyLu/yoZ849xcz6DLoKmLQXTc3yntZwrLJt7a1aWcDm/zEuYfMoMugq4hBdzX2Ff/LPkYn4z6bOTnPtep39dIMugy6ihh0V3v3vne9lWGGc4X/VL+nV2bQZdBVxKC7XthxzgbIOp6vfkuvzKDLoKuIQXe9shWsDpALfF39ll6YQZdBVxGD7npln/J1gDzJquq39MIMugy6ihh010t7jhcC5Fk25DLfVb+n5gy6DLqKGHTXS+swNvnSO5AX2Jzr3K5+V40ZdBl0FTHorhf3by5nGcsDZAELc5Ufq9/U7Rl0GXQVMeiul7eHQwEyj/nVb+n2DLoMuooYdNfr28beABnlRPVbujmDLoOuIgbd9fo6jKXNUOazsPot3ZxBl0FXEYPu+mEv82qAvM/56rd0awZdBl1FDLrrh13jp0Ary1he/ZZuzaDLoKuIQXf9snX3v1N9UL7ExaDLoKuIQXf9skt8EyBreab6Ld2YQZdBVxGD7vppEx868yXfV79lpmfQZdBVxKC7fto5Pg+QFayufstMz6DLoKuIQXf9tvVsHIj3pRt0GXQVMeiu33adnzPMrMxidjqMV79npmbQZdBVxKC7ftw+RgNkL4er3zJTM+gy6Cpi0F0/rsN4ZjM7w8xq7FO6QZdBVxGD7vp1BzgSILs4WP2WmZhBl0FXEYPu+nX3ntLnZJjhRj6lG3QZdBUx6K6fd4hjAbKDN6rfMt0z6DLoKmLQXb9vDnMzxHA6jFW/ZTpn0GXQVcSgu37fm7wbIAc5Wv2W6ZxBl0FXEYPu+n0dxgLkKdZVv2U6Z9Bl0FXEoLsmbAGLMps51e+Yzhl0GXQVMeiuCXuRzQFyhVvVb5muGXQZdBUx6K4JO8lHAbKP0eq3TNcMugy6ihh014R1GA+0spqnq98yXTPoMugqYtBdU7aIJZnFrOp3TNcMugy6ihh015RtYmuAXOa76rdMxwy6DLqKGHTXlJ3mQoDs4VD1W6ZjBl0GXUUMumvKOoynRSurWFv9lumYQZdBVxGD7pq0xSzNcEP+HN2gy6CriEF3TdpmXguQS3xT/ZbHnUGXQVcRg+6atLN8FiA72V/9lsedQZdBVxGD7pq2Fq08yerqdzzuDLoMuooYdNe0LeGJDDFc/Y7HnUGXQVcRg+6atlfYFiDnuVb9lseZQZdBVxGD7pq2i9y4/3Wqz1a/5XFm0GXQVcSguyZuBav7/lPjDLoMuooYdNfEfcTlAH39ITMGXQZdRQy6a+rWszFAdvBG9Vv+yQy6DLqKGHTX1HUYzwIWBcgZLlS/p3QGXQZdRQy6a/K+4Nu0aWeI4Vzlh+r3lMygy6CriEF3Td8JzgXI87xU/ZaSGXQZdBUx6G4QtoQn0qKV69yufsujzqDLoKuIQXeDsGN8GCC7++i70g26DLqKGHQ3COswniGGsogl1W951Bl0GXQVMehuULaKtWnRSofx6rc8ygy6DLqKGHQ3KNvJgQA5y6XqtzzKDLoMuooYdDcom/iM961sr37Lo8ygy6CriEF3g7QhhrKMJ6vf8Sgz6DLoKmLQ3SBtGU9miKHqdzzKDLoMuooYdDdI28r2ALnIjeq3/N0Mugy6ihh0N0g7y6UA2cXB6rf83Qy6DLqKGHQ3SOswnhatvvhaVYMug64iBt0N2payPC3a6TBW/Za/mkGXQVcRg+4GbYc4FiCHOFb9lr+aQZdBVxGD7gZtHcbSZiizmN3TX9Zi0GXQVcSgu0HcKO8FyDOsr37Ln82gy6CriEF3g7pVrA2Qz3r0LWwGXQZdRQy6G9Rd4FqArGND9VseNoMug64iBt0N8paxPNDKNX6qfsvvZ9Bl0FXEoLtB3vt8EiBbeL36Lb+fQZdBVxGD7gZ981iQIYZ67nvSDboMuooYdDfoO8A7ATLKieq3TJ1Bl0FXEYPuBn3X+TlAnua56rdMnUGXQVcRg+7cL1nAwsxmTvU7ps6gy6CriEF37pe8wOYAucKt6rdMzKDLoKuIQXful5zi4wDZy2j1WyZm0GXQVcSgO9ebX6tq0GXQVcSgO3dvS1iaYYar3zExgy6DriIG3bl728K2nvpsd4Mug64iBt25e3uPcz31fnSDLoOuIgbduXu7wq0A2cr26rfcxKDLoKuQQXfut0GrZz5gxqDLoKuIQXfutw0znCdYUf2Omxh0GXQVMujO/bYhhrKC1dXvuIlBl0FXIYPu3L11GAu0soZnqt9yE4Mug65CBt25e3uZkQA5zPHqt9zEoMugq5BBd03f+5zPMpanRSstWpnH/Oxk/+T3n3cYz1a2B8hSnqh+78QMugy6ihh01+Rt5JUAadHKStbkKdZlFrMDpE07C1iUIYYCZAnL8hX/rX7zxAy6DLqKGHTX1D3HCwHyJKtylR8f+GuHOJYnWJE5zM1ilmY/b08+sffKDLoMuopMBP1tTuYTrjjXiK1ibYCsYHU+4svq9/yTHeOMQR9wBl1Ffv311wDOuR7c/Pnzc/fu3dq/TagSg65iv/76a+7cueOc67EZ88Fm0CVJagCDLklSAxh0SZIawKBLktQABl2SpAYw6JIkNYBBlySpAQy6JEkNYNAlSWoAgy5JUgMYdEmSGsCgS5LUAAZdkqQGMOiSJDWAQZckqQEMuiRJDWDQJUlqAIMuSVIDGHRJkhrAoEuS1AAGXZKkBjDokiQ1gEGXJKkBDLokSQ1g0CVJagCDLklSAxh0SZIawKBLktQABl2SpAYw6JIkNYBBlySpAf4/yPa3qGYu4iwAAAAASUVORK5CYII=\" width=\"500\">" |
|
], |
|
"text/plain": [ |
|
"<IPython.core.display.HTML object>" |
|
] |
|
}, |
|
"metadata": {}, |
|
"output_type": "display_data" |
|
} |
|
], |
|
"source": [ |
|
"fig = plt.figure(figsize=(5,3))" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": 3, |
|
"metadata": { |
|
"scrolled": false |
|
}, |
|
"outputs": [], |
|
"source": [ |
|
"# Plot the figure, define the geographic bounds\n", |
|
"\n", |
|
"# fig = plt.figure(figsize=(9,6))\n", |
|
"\n", |
|
"# fig = plt.figure()\n", |
|
"\n", |
|
"# latcorners = ([-47.72086687985285,-47.18137037025522])\n", |
|
"\n", |
|
"# loncorners = ([-16.627150415017045,-16.329425174039407])\n", |
|
"\n", |
|
"# ax = plt.axes(projection=ccrs.PlateCarree())\n", |
|
"# ax.coastlines()\n", |
|
"# ax.add_feature(cfeature.BORDERS)\n", |
|
"# ax.add_feature(cfeature.STATES.with_scale('10m'))\n", |
|
"\n", |
|
"proj = ccrs.PlateCarree()\n", |
|
"ax = plt.axes(projection = proj)\n", |
|
"ax.coastlines()\n", |
|
"ax.add_feature(cfeature.BORDERS)\n", |
|
"ax.add_feature(cfeature.STATES.with_scale('10m'))\n", |
|
"\n", |
|
"extent = [-48, -47, -16, -17]\n", |
|
"\n", |
|
"ax.set_extent(extent, crs=ccrs.PlateCarree())" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": 6, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"# Draw filled contours.\n", |
|
"\n", |
|
"clevs = np.arange(0,1.26,0.125)" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": 7, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"# Define the latitude and longitude data\n", |
|
"\n", |
|
"x, y = np.float32(np.meshgrid(theLons, theLats))" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": 8, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"# Mask the values less than 0 because there is no data to plot.\n", |
|
"\n", |
|
"masked_array = np.ma.masked_where(precip < 0,precip)" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": 9, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"cmap = 'nipy_spectral'" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": 10, |
|
"metadata": { |
|
"scrolled": false |
|
}, |
|
"outputs": [], |
|
"source": [ |
|
"# Plot the data\n", |
|
"\n", |
|
"cs = ax.contourf(x,y,precip,clevs,transform=ccrs.PlateCarree(),cmap=cmap)" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": null, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"# Set the title and fonts\n", |
|
"\n", |
|
"plt.title('August 2015 Monthly Average Rain Rate')\n", |
|
"\n", |
|
"font = {'weight' : 'bold', 'size' : 3}\n", |
|
"\n", |
|
"plt.rc('font', **font)" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": null, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"sm = plt.cm.ScalarMappable(cmap=cmap,norm=plt.Normalize(0,1))\n", |
|
"sm._A = []\n", |
|
"plt.colorbar(sm, ax=ax, label='mm/h', shrink=0.5)" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": null, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"# ax.set_extent([-180, 180, -90, 90])" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": null, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"fname = r'./layers/POLYGON.shp'\n", |
|
"\n", |
|
"# ax = plt.axes(projection=ccrs.Robinson())\n", |
|
"shape_feature = ShapelyFeature(Reader(fname).geometries(),\n", |
|
" ccrs.PlateCarree(), edgecolor='black')\n", |
|
"ax.add_feature(shape_feature)" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": null, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"precip[:1]" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": null, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"precip[732:736,3368:3373]" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": null, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"plt.show()" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": null, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [ |
|
"matplotlib.pyplot.show()" |
|
] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": null, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": null, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [] |
|
}, |
|
{ |
|
"cell_type": "code", |
|
"execution_count": null, |
|
"metadata": {}, |
|
"outputs": [], |
|
"source": [] |
|
} |
|
], |
|
"metadata": { |
|
"kernelspec": { |
|
"display_name": "Python 3", |
|
"language": "python", |
|
"name": "python3" |
|
}, |
|
"language_info": { |
|
"codemirror_mode": { |
|
"name": "ipython", |
|
"version": 3 |
|
}, |
|
"file_extension": ".py", |
|
"mimetype": "text/x-python", |
|
"name": "python", |
|
"nbconvert_exporter": "python", |
|
"pygments_lexer": "ipython3", |
|
"version": "3.5.2" |
|
} |
|
}, |
|
"nbformat": 4, |
|
"nbformat_minor": 2 |
|
}
|
|
|