From 92a1f3c822b80acdb700de7bbc412f3de35d3d5e Mon Sep 17 00:00:00 2001 From: G Paul Bailey Date: Sun, 8 Sep 2024 10:33:28 -0600 Subject: [PATCH 1/2] Add files via upload upload the code I saved from codespace --- first-cu-map.ipynb | 570 +++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 570 insertions(+) create mode 100644 first-cu-map.ipynb diff --git a/first-cu-map.ipynb b/first-cu-map.ipynb new file mode 100644 index 0000000..ec61463 --- /dev/null +++ b/first-cu-map.ipynb @@ -0,0 +1,570 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "id": "4681aaaa-53e1-4010-8fb0-39d1b85c34a8", + "metadata": {}, + "source": [ + "# Add a map to your website\n", + "\n", + "Get started with maps\n", + "\n", + "Nate Quarderer \n", + "Elsa Culler \n", + "2024-03-28\n", + "\n", + "## Get started with map-making using open-sources tools\n", + "\n", + "### Finding locations and boundaries\n", + "\n", + "[Open Street Map (OSM)](https://www.openstreetmap.org/) is an\n", + "open-source, editable map of the world – a little like a wiki for\n", + "places. They also provide a service for looking up locations using text,\n", + "which we’ll be using in this activity.\n", + "\n", + "### Mapping libraries\n", + "\n", + "You’ll need to start by importing the following libraries:" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "7ca72713", + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "94e40637", + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "07402822", + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "042060e0", + "metadata": {}, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": 1, + "id": "cd830589", + "metadata": {}, + "outputs": [ + { + "data": { + "application/javascript": "(function(root) {\n function now() {\n return new Date();\n }\n\n var force = true;\n var py_version = '3.3.4'.replace('rc', '-rc.').replace('.dev', '-dev.');\n var reloading = false;\n var Bokeh = root.Bokeh;\n\n if (typeof (root._bokeh_timeout) === \"undefined\" || force) {\n root._bokeh_timeout = Date.now() + 5000;\n root._bokeh_failed_load = false;\n }\n\n function run_callbacks() {\n try {\n root._bokeh_onload_callbacks.forEach(function(callback) {\n if (callback != null)\n callback();\n });\n } finally {\n delete root._bokeh_onload_callbacks;\n }\n console.debug(\"Bokeh: all callbacks have finished\");\n }\n\n function load_libs(css_urls, js_urls, js_modules, js_exports, callback) {\n if (css_urls == null) css_urls = [];\n if (js_urls == null) js_urls = [];\n if (js_modules == null) js_modules = [];\n if (js_exports == null) js_exports = {};\n\n root._bokeh_onload_callbacks.push(callback);\n\n if (root._bokeh_is_loading > 0) {\n console.debug(\"Bokeh: BokehJS is being loaded, scheduling callback at\", now());\n return null;\n }\n if (js_urls.length === 0 && js_modules.length === 0 && Object.keys(js_exports).length === 0) {\n run_callbacks();\n return null;\n }\n if (!reloading) {\n console.debug(\"Bokeh: BokehJS not loaded, scheduling load and callback at\", now());\n }\n\n function on_load() {\n root._bokeh_is_loading--;\n if (root._bokeh_is_loading === 0) {\n console.debug(\"Bokeh: all BokehJS libraries/stylesheets loaded\");\n run_callbacks()\n }\n }\n window._bokeh_on_load = on_load\n\n function on_error() {\n console.error(\"failed to load \" + url);\n }\n\n var skip = [];\n if (window.requirejs) {\n window.requirejs.config({'packages': {}, 'paths': {'jspanel': 'https://cdn.jsdelivr.net/npm/jspanel4@4.12.0/dist/jspanel', 'jspanel-modal': 'https://cdn.jsdelivr.net/npm/jspanel4@4.12.0/dist/extensions/modal/jspanel.modal', 'jspanel-tooltip': 'https://cdn.jsdelivr.net/npm/jspanel4@4.12.0/dist/extensions/tooltip/jspanel.tooltip', 'jspanel-hint': 'https://cdn.jsdelivr.net/npm/jspanel4@4.12.0/dist/extensions/hint/jspanel.hint', 'jspanel-layout': 'https://cdn.jsdelivr.net/npm/jspanel4@4.12.0/dist/extensions/layout/jspanel.layout', 'jspanel-contextmenu': 'https://cdn.jsdelivr.net/npm/jspanel4@4.12.0/dist/extensions/contextmenu/jspanel.contextmenu', 'jspanel-dock': 'https://cdn.jsdelivr.net/npm/jspanel4@4.12.0/dist/extensions/dock/jspanel.dock', 'gridstack': 'https://cdn.jsdelivr.net/npm/gridstack@7.2.3/dist/gridstack-all', 'notyf': 'https://cdn.jsdelivr.net/npm/notyf@3/notyf.min'}, 'shim': {'jspanel': {'exports': 'jsPanel'}, 'gridstack': {'exports': 'GridStack'}}});\n require([\"jspanel\"], function(jsPanel) {\n\twindow.jsPanel = jsPanel\n\ton_load()\n })\n require([\"jspanel-modal\"], function() {\n\ton_load()\n })\n require([\"jspanel-tooltip\"], function() {\n\ton_load()\n })\n require([\"jspanel-hint\"], function() {\n\ton_load()\n })\n require([\"jspanel-layout\"], function() {\n\ton_load()\n })\n require([\"jspanel-contextmenu\"], function() {\n\ton_load()\n })\n require([\"jspanel-dock\"], function() {\n\ton_load()\n })\n require([\"gridstack\"], function(GridStack) {\n\twindow.GridStack = GridStack\n\ton_load()\n })\n require([\"notyf\"], function() {\n\ton_load()\n })\n root._bokeh_is_loading = css_urls.length + 9;\n } else {\n root._bokeh_is_loading = css_urls.length + js_urls.length + js_modules.length + Object.keys(js_exports).length;\n }\n\n var existing_stylesheets = []\n var links = document.getElementsByTagName('link')\n for (var i = 0; i < links.length; i++) {\n var link = links[i]\n if (link.href != null) {\n\texisting_stylesheets.push(link.href)\n }\n }\n for (var i = 0; i < css_urls.length; i++) {\n var url = css_urls[i];\n if (existing_stylesheets.indexOf(url) !== -1) {\n\ton_load()\n\tcontinue;\n }\n const element = document.createElement(\"link\");\n element.onload = on_load;\n element.onerror = on_error;\n element.rel = \"stylesheet\";\n element.type = \"text/css\";\n element.href = url;\n console.debug(\"Bokeh: injecting link tag for BokehJS stylesheet: \", url);\n document.body.appendChild(element);\n } if (((window['jsPanel'] !== undefined) && (!(window['jsPanel'] instanceof HTMLElement))) || window.requirejs) {\n var urls = ['https://cdn.holoviz.org/panel/1.3.8/dist/bundled/floatpanel/jspanel4@4.12.0/dist/jspanel.js', 'https://cdn.holoviz.org/panel/1.3.8/dist/bundled/floatpanel/jspanel4@4.12.0/dist/extensions/modal/jspanel.modal.js', 'https://cdn.holoviz.org/panel/1.3.8/dist/bundled/floatpanel/jspanel4@4.12.0/dist/extensions/tooltip/jspanel.tooltip.js', 'https://cdn.holoviz.org/panel/1.3.8/dist/bundled/floatpanel/jspanel4@4.12.0/dist/extensions/hint/jspanel.hint.js', 'https://cdn.holoviz.org/panel/1.3.8/dist/bundled/floatpanel/jspanel4@4.12.0/dist/extensions/layout/jspanel.layout.js', 'https://cdn.holoviz.org/panel/1.3.8/dist/bundled/floatpanel/jspanel4@4.12.0/dist/extensions/contextmenu/jspanel.contextmenu.js', 'https://cdn.holoviz.org/panel/1.3.8/dist/bundled/floatpanel/jspanel4@4.12.0/dist/extensions/dock/jspanel.dock.js'];\n for (var i = 0; i < urls.length; i++) {\n skip.push(urls[i])\n }\n } if (((window['GridStack'] !== undefined) && (!(window['GridStack'] instanceof HTMLElement))) || window.requirejs) {\n var urls = ['https://cdn.holoviz.org/panel/1.3.8/dist/bundled/gridstack/gridstack@7.2.3/dist/gridstack-all.js'];\n for (var i = 0; i < urls.length; i++) {\n skip.push(urls[i])\n }\n } if (((window['Notyf'] !== undefined) && (!(window['Notyf'] instanceof HTMLElement))) || window.requirejs) {\n var urls = ['https://cdn.holoviz.org/panel/1.3.8/dist/bundled/notificationarea/notyf@3/notyf.min.js'];\n for (var i = 0; i < urls.length; i++) {\n skip.push(urls[i])\n }\n } var existing_scripts = []\n var scripts = document.getElementsByTagName('script')\n for (var i = 0; i < scripts.length; i++) {\n var script = scripts[i]\n if (script.src != null) {\n\texisting_scripts.push(script.src)\n }\n }\n for (var i = 0; i < js_urls.length; i++) {\n var url = js_urls[i];\n if (skip.indexOf(url) !== -1 || existing_scripts.indexOf(url) !== -1) {\n\tif (!window.requirejs) {\n\t on_load();\n\t}\n\tcontinue;\n }\n var element = document.createElement('script');\n element.onload = on_load;\n element.onerror = on_error;\n element.async = false;\n element.src = url;\n console.debug(\"Bokeh: injecting script tag for BokehJS library: \", url);\n document.head.appendChild(element);\n }\n for (var i = 0; i < js_modules.length; i++) {\n var url = js_modules[i];\n if (skip.indexOf(url) !== -1 || existing_scripts.indexOf(url) !== -1) {\n\tif (!window.requirejs) {\n\t on_load();\n\t}\n\tcontinue;\n }\n var element = document.createElement('script');\n element.onload = on_load;\n element.onerror = on_error;\n element.async = false;\n element.src = url;\n element.type = \"module\";\n console.debug(\"Bokeh: injecting script tag for BokehJS library: \", url);\n document.head.appendChild(element);\n }\n for (const name in js_exports) {\n var url = js_exports[name];\n if (skip.indexOf(url) >= 0 || root[name] != null) {\n\tif (!window.requirejs) {\n\t on_load();\n\t}\n\tcontinue;\n }\n var element = document.createElement('script');\n element.onerror = on_error;\n element.async = false;\n element.type = \"module\";\n console.debug(\"Bokeh: injecting script tag for BokehJS library: \", url);\n element.textContent = `\n import ${name} from \"${url}\"\n window.${name} = ${name}\n window._bokeh_on_load()\n `\n document.head.appendChild(element);\n }\n if (!js_urls.length && !js_modules.length) {\n on_load()\n }\n };\n\n function inject_raw_css(css) {\n const element = document.createElement(\"style\");\n element.appendChild(document.createTextNode(css));\n document.body.appendChild(element);\n }\n\n var js_urls = [\"https://cdn.bokeh.org/bokeh/release/bokeh-3.3.4.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-gl-3.3.4.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-widgets-3.3.4.min.js\", \"https://cdn.bokeh.org/bokeh/release/bokeh-tables-3.3.4.min.js\", \"https://cdn.holoviz.org/panel/1.3.8/dist/panel.min.js\"];\n var js_modules = [];\n var js_exports = {};\n var css_urls = [];\n var inline_js = [ function(Bokeh) {\n Bokeh.set_log_level(\"info\");\n },\nfunction(Bokeh) {} // ensure no trailing comma for IE\n ];\n\n function run_inline_js() {\n if ((root.Bokeh !== undefined) || (force === true)) {\n for (var i = 0; i < inline_js.length; i++) {\n\ttry {\n inline_js[i].call(root, root.Bokeh);\n\t} catch(e) {\n\t if (!reloading) {\n\t throw e;\n\t }\n\t}\n }\n // Cache old bokeh versions\n if (Bokeh != undefined && !reloading) {\n\tvar NewBokeh = root.Bokeh;\n\tif (Bokeh.versions === undefined) {\n\t Bokeh.versions = new Map();\n\t}\n\tif (NewBokeh.version !== Bokeh.version) {\n\t Bokeh.versions.set(NewBokeh.version, NewBokeh)\n\t}\n\troot.Bokeh = Bokeh;\n }} else if (Date.now() < root._bokeh_timeout) {\n setTimeout(run_inline_js, 100);\n } else if (!root._bokeh_failed_load) {\n console.log(\"Bokeh: BokehJS failed to load within specified timeout.\");\n root._bokeh_failed_load = true;\n }\n root._bokeh_is_initializing = false\n }\n\n function load_or_wait() {\n // Implement a backoff loop that tries to ensure we do not load multiple\n // versions of Bokeh and its dependencies at the same time.\n // In recent versions we use the root._bokeh_is_initializing flag\n // to determine whether there is an ongoing attempt to initialize\n // bokeh, however for backward compatibility we also try to ensure\n // that we do not start loading a newer (Panel>=1.0 and Bokeh>3) version\n // before older versions are fully initialized.\n if (root._bokeh_is_initializing && Date.now() > root._bokeh_timeout) {\n root._bokeh_is_initializing = false;\n root._bokeh_onload_callbacks = undefined;\n console.log(\"Bokeh: BokehJS was loaded multiple times but one version failed to initialize.\");\n load_or_wait();\n } else if (root._bokeh_is_initializing || (typeof root._bokeh_is_initializing === \"undefined\" && root._bokeh_onload_callbacks !== undefined)) {\n setTimeout(load_or_wait, 100);\n } else {\n root._bokeh_is_initializing = true\n root._bokeh_onload_callbacks = []\n var bokeh_loaded = Bokeh != null && (Bokeh.version === py_version || (Bokeh.versions !== undefined && Bokeh.versions.has(py_version)));\n if (!reloading && !bokeh_loaded) {\n\troot.Bokeh = undefined;\n }\n load_libs(css_urls, js_urls, js_modules, js_exports, function() {\n\tconsole.debug(\"Bokeh: BokehJS plotting callback run at\", now());\n\trun_inline_js();\n });\n }\n }\n // Give older versions of the autoload script a head-start to ensure\n // they initialize before we start loading newer version.\n setTimeout(load_or_wait, 100)\n}(window));", + "application/vnd.holoviews_load.v0+json": "" + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "application/javascript": "\nif ((window.PyViz === undefined) || (window.PyViz instanceof HTMLElement)) {\n window.PyViz = {comms: {}, comm_status:{}, kernels:{}, receivers: {}, plot_index: []}\n}\n\n\n function JupyterCommManager() {\n }\n\n JupyterCommManager.prototype.register_target = function(plot_id, comm_id, msg_handler) {\n if (window.comm_manager || ((window.Jupyter !== undefined) && (Jupyter.notebook.kernel != null))) {\n var comm_manager = window.comm_manager || Jupyter.notebook.kernel.comm_manager;\n comm_manager.register_target(comm_id, function(comm) {\n comm.on_msg(msg_handler);\n });\n } else if ((plot_id in window.PyViz.kernels) && (window.PyViz.kernels[plot_id])) {\n window.PyViz.kernels[plot_id].registerCommTarget(comm_id, function(comm) {\n comm.onMsg = msg_handler;\n });\n } else if (typeof google != 'undefined' && google.colab.kernel != null) {\n google.colab.kernel.comms.registerTarget(comm_id, (comm) => {\n var messages = comm.messages[Symbol.asyncIterator]();\n function processIteratorResult(result) {\n var message = result.value;\n console.log(message)\n var content = {data: message.data, comm_id};\n var buffers = []\n for (var buffer of message.buffers || []) {\n buffers.push(new DataView(buffer))\n }\n var metadata = message.metadata || {};\n var msg = {content, buffers, metadata}\n msg_handler(msg);\n return messages.next().then(processIteratorResult);\n }\n return messages.next().then(processIteratorResult);\n })\n }\n }\n\n JupyterCommManager.prototype.get_client_comm = function(plot_id, comm_id, msg_handler) {\n if (comm_id in window.PyViz.comms) {\n return window.PyViz.comms[comm_id];\n } else if (window.comm_manager || ((window.Jupyter !== undefined) && (Jupyter.notebook.kernel != null))) {\n var comm_manager = window.comm_manager || Jupyter.notebook.kernel.comm_manager;\n var comm = comm_manager.new_comm(comm_id, {}, {}, {}, comm_id);\n if (msg_handler) {\n comm.on_msg(msg_handler);\n }\n } else if ((plot_id in window.PyViz.kernels) && (window.PyViz.kernels[plot_id])) {\n var comm = window.PyViz.kernels[plot_id].connectToComm(comm_id);\n comm.open();\n if (msg_handler) {\n comm.onMsg = msg_handler;\n }\n } else if (typeof google != 'undefined' && google.colab.kernel != null) {\n var comm_promise = google.colab.kernel.comms.open(comm_id)\n comm_promise.then((comm) => {\n window.PyViz.comms[comm_id] = comm;\n if (msg_handler) {\n var messages = comm.messages[Symbol.asyncIterator]();\n function processIteratorResult(result) {\n var message = result.value;\n var content = {data: message.data};\n var metadata = message.metadata || {comm_id};\n var msg = {content, metadata}\n msg_handler(msg);\n return messages.next().then(processIteratorResult);\n }\n return messages.next().then(processIteratorResult);\n }\n }) \n var sendClosure = (data, metadata, buffers, disposeOnDone) => {\n return comm_promise.then((comm) => {\n comm.send(data, metadata, buffers, disposeOnDone);\n });\n };\n var comm = {\n send: sendClosure\n };\n }\n window.PyViz.comms[comm_id] = comm;\n return comm;\n }\n window.PyViz.comm_manager = new JupyterCommManager();\n \n\n\nvar JS_MIME_TYPE = 'application/javascript';\nvar HTML_MIME_TYPE = 'text/html';\nvar EXEC_MIME_TYPE = 'application/vnd.holoviews_exec.v0+json';\nvar CLASS_NAME = 'output';\n\n/**\n * Render data to the DOM node\n */\nfunction render(props, node) {\n var div = document.createElement(\"div\");\n var script = document.createElement(\"script\");\n node.appendChild(div);\n node.appendChild(script);\n}\n\n/**\n * Handle when a new output is added\n */\nfunction handle_add_output(event, handle) {\n var output_area = handle.output_area;\n var output = handle.output;\n if ((output.data == undefined) || (!output.data.hasOwnProperty(EXEC_MIME_TYPE))) {\n return\n }\n var id = output.metadata[EXEC_MIME_TYPE][\"id\"];\n var toinsert = output_area.element.find(\".\" + CLASS_NAME.split(' ')[0]);\n if (id !== undefined) {\n var nchildren = toinsert.length;\n var html_node = toinsert[nchildren-1].children[0];\n html_node.innerHTML = output.data[HTML_MIME_TYPE];\n var scripts = [];\n var nodelist = html_node.querySelectorAll(\"script\");\n for (var i in nodelist) {\n if (nodelist.hasOwnProperty(i)) {\n scripts.push(nodelist[i])\n }\n }\n\n scripts.forEach( function (oldScript) {\n var newScript = document.createElement(\"script\");\n var attrs = [];\n var nodemap = oldScript.attributes;\n for (var j in nodemap) {\n if (nodemap.hasOwnProperty(j)) {\n attrs.push(nodemap[j])\n }\n }\n attrs.forEach(function(attr) { newScript.setAttribute(attr.name, attr.value) });\n newScript.appendChild(document.createTextNode(oldScript.innerHTML));\n oldScript.parentNode.replaceChild(newScript, oldScript);\n });\n if (JS_MIME_TYPE in output.data) {\n toinsert[nchildren-1].children[1].textContent = output.data[JS_MIME_TYPE];\n }\n output_area._hv_plot_id = id;\n if ((window.Bokeh !== undefined) && (id in Bokeh.index)) {\n window.PyViz.plot_index[id] = Bokeh.index[id];\n } else {\n window.PyViz.plot_index[id] = null;\n }\n } else if (output.metadata[EXEC_MIME_TYPE][\"server_id\"] !== undefined) {\n var bk_div = document.createElement(\"div\");\n bk_div.innerHTML = output.data[HTML_MIME_TYPE];\n var script_attrs = bk_div.children[0].attributes;\n for (var i = 0; i < script_attrs.length; i++) {\n toinsert[toinsert.length - 1].childNodes[1].setAttribute(script_attrs[i].name, script_attrs[i].value);\n }\n // store reference to server id on output_area\n output_area._bokeh_server_id = output.metadata[EXEC_MIME_TYPE][\"server_id\"];\n }\n}\n\n/**\n * Handle when an output is cleared or removed\n */\nfunction handle_clear_output(event, handle) {\n var id = handle.cell.output_area._hv_plot_id;\n var server_id = handle.cell.output_area._bokeh_server_id;\n if (((id === undefined) || !(id in PyViz.plot_index)) && (server_id !== undefined)) { return; }\n var comm = window.PyViz.comm_manager.get_client_comm(\"hv-extension-comm\", \"hv-extension-comm\", function () {});\n if (server_id !== null) {\n comm.send({event_type: 'server_delete', 'id': server_id});\n return;\n } else if (comm !== null) {\n comm.send({event_type: 'delete', 'id': id});\n }\n delete PyViz.plot_index[id];\n if ((window.Bokeh !== undefined) & (id in window.Bokeh.index)) {\n var doc = window.Bokeh.index[id].model.document\n doc.clear();\n const i = window.Bokeh.documents.indexOf(doc);\n if (i > -1) {\n window.Bokeh.documents.splice(i, 1);\n }\n }\n}\n\n/**\n * Handle kernel restart event\n */\nfunction handle_kernel_cleanup(event, handle) {\n delete PyViz.comms[\"hv-extension-comm\"];\n window.PyViz.plot_index = {}\n}\n\n/**\n * Handle update_display_data messages\n */\nfunction handle_update_output(event, handle) {\n handle_clear_output(event, {cell: {output_area: handle.output_area}})\n handle_add_output(event, handle)\n}\n\nfunction register_renderer(events, OutputArea) {\n function append_mime(data, metadata, element) {\n // create a DOM node to render to\n var toinsert = this.create_output_subarea(\n metadata,\n CLASS_NAME,\n EXEC_MIME_TYPE\n );\n this.keyboard_manager.register_events(toinsert);\n // Render to node\n var props = {data: data, metadata: metadata[EXEC_MIME_TYPE]};\n render(props, toinsert[0]);\n element.append(toinsert);\n return toinsert\n }\n\n events.on('output_added.OutputArea', handle_add_output);\n events.on('output_updated.OutputArea', handle_update_output);\n events.on('clear_output.CodeCell', handle_clear_output);\n events.on('delete.Cell', handle_clear_output);\n events.on('kernel_ready.Kernel', handle_kernel_cleanup);\n\n OutputArea.prototype.register_mime_type(EXEC_MIME_TYPE, append_mime, {\n safe: true,\n index: 0\n });\n}\n\nif (window.Jupyter !== undefined) {\n try {\n var events = require('base/js/events');\n var OutputArea = require('notebook/js/outputarea').OutputArea;\n if (OutputArea.prototype.mime_types().indexOf(EXEC_MIME_TYPE) == -1) {\n register_renderer(events, OutputArea);\n }\n } catch(err) {\n }\n}\n", + "application/vnd.holoviews_load.v0+json": "" + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "text/html": [ + "" + ] + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "application/vnd.holoviews_exec.v0+json": "", + "text/html": [ + "
\n", + "
\n", + "
\n", + "" + ] + }, + "metadata": { + "application/vnd.holoviews_exec.v0+json": { + "id": "p1002" + } + }, + "output_type": "display_data" + } + ], + "source": [ + "# Work with vector data\n", + "import geopandas as gpd\n", + "\n", + "# Save maps and plots to files\n", + "import holoviews as hv\n", + "# Create interactive maps and plots\n", + "import hvplot.pandas\n", + "\n", + "# Search for locations by name - this might take a moment\n", + "from osmnx import features as osm" + ] + }, + { + "cell_type": "markdown", + "id": "f9505692-f5e7-4163-9ba7-7a08773dc238", + "metadata": {}, + "source": [ + "### Search for a point of interest\n", + "\n", + "You can use the `pyrosm` package to download and search for spatial\n", + "vector data in your area, or anywhere around the world.\n", + "\n", + "In this case, we’re looking for the location of the United Tribes\n", + "Technical College campus in North Dakota. The address in here,\n", + "`'United Tribes Technical College, Bismarck, ND, United States'`, does\n", + "not have to be complete or exact, but it should be specific enough to\n", + "narrow it down. We are also specifying that we want it to be tagges as a\n", + "`'amentity'` type, specifically a `'college'` You might have to try a\n", + "couple different searches with different addresses and/or tags to get\n", + "the address you want, just like if you are using a map website or app." + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "id": "c4a50a59", + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
nodesamenityboundaryinternet_accessnameoperatorwebsitegeometryelegnis:feature_idwikidatawikipedia
element_typeosmid
way46226108[589514700, 2188984094, 2188984095, 326889012,...universityadministrativewlanUniversity of Colorado Boulder (Main Campus)University of Colorado Boulderhttps://www.colorado.edu/POLYGON ((-105.27589 40.01047, -105.27597 40.0...NaNNaNNaNNaN
391232001[3944423259, 3944423260, 3944423261, 394442326...universityNaNNaNNaropa UniversityNaropa UniversityNaNPOLYGON ((-105.26632 40.01452, -105.26631 40.0...1620178660Q2975783en:Naropa University
\n", + "
" + ], + "text/plain": [ + " nodes \\\n", + "element_type osmid \n", + "way 46226108 [589514700, 2188984094, 2188984095, 326889012,... \n", + " 391232001 [3944423259, 3944423260, 3944423261, 394442326... \n", + "\n", + " amenity boundary internet_access \\\n", + "element_type osmid \n", + "way 46226108 university administrative wlan \n", + " 391232001 university NaN NaN \n", + "\n", + " name \\\n", + "element_type osmid \n", + "way 46226108 University of Colorado Boulder (Main Campus) \n", + " 391232001 Naropa University \n", + "\n", + " operator \\\n", + "element_type osmid \n", + "way 46226108 University of Colorado Boulder \n", + " 391232001 Naropa University \n", + "\n", + " website \\\n", + "element_type osmid \n", + "way 46226108 https://www.colorado.edu/ \n", + " 391232001 NaN \n", + "\n", + " geometry \\\n", + "element_type osmid \n", + "way 46226108 POLYGON ((-105.27589 40.01047, -105.27597 40.0... \n", + " 391232001 POLYGON ((-105.26632 40.01452, -105.26631 40.0... \n", + "\n", + " ele gnis:feature_id wikidata wikipedia \n", + "element_type osmid \n", + "way 46226108 NaN NaN NaN NaN \n", + " 391232001 1620 178660 Q2975783 en:Naropa University " + ] + }, + "execution_count": 2, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "# Search for United Tribes Technical College\n", + "ucb_gdf = osm.features_from_address(\n", + " 'University of Colorado, Boulder, CO, United States',\n", + " {'amenity': ['university']})\n", + "ucb_gdf" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "id": "8aa08622", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 3, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbIAAAGvCAYAAAA+O1BZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAA9hAAAPYQGoP6dpAABPA0lEQVR4nO3de1zUVf4/8NcwwwwXh+EWMw43EfGCiAgoQSoYirpuae0mpovdv5aakG0am/4yMxG/5Zpf11z329buWuL2NdJa3cC2SFcygjCzLbuQYIqkGRdFbnN+fygjI9fhMp+5vJ6PxzyUM+ecz3s+fpw35/M5n/ORCSEEiIiIbJST1AEQERH1BRMZERHZNCYyIiKyaUxkRERk05jIiIjIpjGRERGRTWMiIyIim8ZERkRENk0hdQDWxGAw4MyZM1Cr1ZDJZFKHQ0Tk0IQQqK2thV6vh5NT5+MuJrI2zpw5g8DAQKnDICKiNioqKhAQENDp+0xkbajVagBXd5qHh4fE0RARObaamhoEBgYav5s7w0TWRuvpRA8PDyYyIiIr0d2lHk72ICIim8ZERkRENo2JjIiIbBoTGRER2TQmMiIismlMZEREZNOYyIiIyKYxkRERkU1jIiMiIpvGREZERDaNiYyIiGwaExkREdk0JjIiIrJpXP2eyM6VX7iMlwq+gRC4+oK49uf1n2H8WbQpv/4z2ra7oY8t88bBVSmX8BOSo2MiI7JzFy83YtfHFQPWf2OLAa5gIiPp8NQikZ2TO3X9LKe+EkIMaP9E3WEiI7JzTt08lLCvDMxjJDEmMiI7N9AjshZmMpIYExmRnZMP8P9ynlokqTGREdk5nloke8dERmTnFE4D+9+8hSMykhgTGZGdG+A8BgOHZCSxPh3iWVlZkMlkyMjIMJYJIbBmzRro9Xq4uroiKSkJJ06c6LavPXv2IDw8HCqVCuHh4cjNzTV5/8MPP8Rtt90GvV4PmUyGt956q8v+Fi1aBJlMhs2bN/fikxHZj4Ge7GHgiIwk1utEVlRUhB07diAyMtKkfOPGjdi0aRO2bt2KoqIi6HQ6TJs2DbW1tZ32VVhYiNTUVKSlpeHYsWNIS0vD3LlzcfToUWOdS5cuYezYsdi6dWu3sb311ls4evQo9Hp9bz8ekd2Q8xoZ2bleJbK6ujosWLAAf/rTn+Dl5WUsF0Jg8+bNeOqpp3DnnXciIiICf/nLX3D58mW8/vrrnfa3efNmTJs2DZmZmRg5ciQyMzORnJxsMpqaOXMm1q1bhzvvvLPL2H744QcsXboUr732GpydnXvz8YjsihNHZGTnepXIlixZglmzZmHq1Kkm5WVlZaisrERKSoqxTKVSITExEUeOHOm0v8LCQpM2ADB9+vQu23TEYDAgLS0NTzzxBEaPHt1t/YaGBtTU1Ji8iOzNgI/IOCQjiZm91mJOTg5KSkpQVFTU7r3KykoAgFarNSnXarU4depUp31WVlZ22Ka1v57Kzs6GQqHAsmXLelQ/KysLzzzzjFnbILI1Az8iG9Duibpl1oisoqIC6enp2LlzJ1xcXDqtJ7vhN0AhRLuy/mjTVnFxMV588UW8+uqrPW6XmZmJ6upq46uiYuAWViWSioKnFsnOmZXIiouLUVVVhZiYGCgUCigUChQUFGDLli1QKBTGUdWNI6mqqqp2I662dDqd2W1udOjQIVRVVSEoKMgY26lTp/D4449jyJAhHbZRqVTw8PAweRHZGy5RRfbOrESWnJyM48ePo7S01PiKjY3FggULUFpaiqFDh0Kn0yE/P9/YprGxEQUFBUhISOi03/j4eJM2AJCXl9dlmxulpaXhs88+M4lNr9fjiSeewLvvvmvOxySyKwO9sgcHZCQ1s66RqdVqREREmJS5u7vDx8fHWJ6RkYH169cjLCwMYWFhWL9+Pdzc3DB//nxjm4ULF8Lf3x9ZWVkAgPT0dEyePBnZ2dmYPXs29u7di4MHD+Lw4cPGNnV1dfjmm2+MP5eVlaG0tBTe3t4ICgqCj48PfHx8TGJzdnaGTqfDiBEjzPmYRHZlwEdkzGQksX5/sOaKFStQX1+PxYsX4+LFi4iLi0NeXh7UarWxTnl5OZzaLDeQkJCAnJwcrFq1CqtXr0ZoaCh2796NuLg4Y51PPvkEU6ZMMf68fPlyAMA999yDV199tb8/BpHdGOA8xmtkJDmZ4NLVRjU1NdBoNKiurub1MrIrQzP/MWCzC/c8Eo+YYO+B6ZwcWk+/k7nWIpEDGKiFgxVOMrNmFxMNhH4/tUhE5ln79hc4caYaSoUTFE4yOMud4KxwglJ+7edrf3eWX31PIXeC8trfneVOCPUbhMlhvl0mFCcnAC09j0nuJMNNg1TQeqjg5+ECP7UKWg8X489a9dW/e7kpB/w+NaLuMJERSezrqlocLfupT32E+LpjQVwQ7ooNhMa1/dJsPVndI3ywBzb+OhJ+Hir4uKsGfJIIUX9hIiOSWH9Mjy87fwnr/vEfbMo/idlR/rgnIRgjddevKYwc7IFLDc1d9jF+iBci/DV9joXI0pjIiCTWn5eYLje2YNfH5dj1cTlG6tRwVykgw9WZix4uzpDJrm5PBpnJ35+aNQqjBnOCE9kmJjIiiQ3UDctfVnb+6KQbVdZcYSIjm8VZi0QSs4ZLUfWNZswEIbIyTGREErOG6euXmcjIhjGREUlM+jQG1Dd2PRGEyJoxkRFJbKAX9e0JjsjIljGREUlsgBbdMAsTGdkyK/gvROTYrOEaWX0TExnZLiYyIolZx6lFXiMj28VERiQx6dMYTy2SbWMiI5IY7yMj6hsmMiKJWcepRSYysl1coopIYlYx2WOAE1mLQaCpxYBmg0BziwFNLQLNBgOaWwSaDQIthqtlLQYBlcIJYVp1950SXcNERiQxK8hj+O78JTzz9olrieVaomkxoOla4mluETf8/dqfnSSn6+VX65rzHPqoQE+8teSWgfuwZHeYyIgkZg3XyM7XNeCVf38vdRgAAIM5WY8IvEZGJDlruEZmTVoMTGRkHiYyIolZwzUya8JERuZiIiOSGPOYKZ5aJHMxkRFJzBqukVkTjsjIXExkRBLjNTJTzGNkLiYyIokxkZniiIzMxen3ZJOEGddRuqva1dsyAE4892dRTGRkLiYyktzfPjqF1W99LnUYnZo8/CYsu3UYYod497hN6w3Cjc0GNLZcfTW1/r3ZgKZrf15uasFXlbUDGL3t4WQPMhcTGUnO2cpHPB+e/BEfnvwRo/UecFPK0dhsQENrMmoxoKlZGBNTY8vVcg4qeo8jMjIXExlJzlluG5dqT5ypkToEh8ARGZnLNr5ByK45K3gY0nUckZG5+A1CklPKrfvUIlkWExmZi4mMJGcrpxbJMpjHyFz8BiHJMZFRWxyRkbn4DUKSU/DUIrXRwskeZCYmMpKckiMyasPAERmZid8gJDmeWqS2OCIjc/XpGyQrKwsymQwZGRnGMiEE1qxZA71eD1dXVyQlJeHEiRPd9rVnzx6Eh4dDpVIhPDwcubm5Ju9/+OGHuO2226DX6yGTyfDWW2+ZvN/U1ISVK1dizJgxcHd3h16vx8KFC3HmzJm+fESyACYyaksI85YgI+r1N0hRURF27NiByMhIk/KNGzdi06ZN2Lp1K4qKiqDT6TBt2jTU1na+DE9hYSFSU1ORlpaGY8eOIS0tDXPnzsXRo0eNdS5duoSxY8di69atHfZx+fJllJSUYPXq1SgpKcGbb76JkydP4vbbb+/tRyQLUSp4jYxMccIHmUMmevGrT11dHaKjo7Ft2zasW7cOUVFR2Lx5M4QQ0Ov1yMjIwMqVKwEADQ0N0Gq1yM7OxqJFizrsLzU1FTU1NThw4ICxbMaMGfDy8sKuXbvaBy2TITc3F3PmzOkyzqKiIkyYMAGnTp1CUFBQt5+rpqYGGo0G1dXV8PDw6LY+9Y9TFy4h8b8/kDoMsiJfrZsBlUIudRgksZ5+J/dqRLZkyRLMmjULU6dONSkvKytDZWUlUlJSjGUqlQqJiYk4cuRIp/0VFhaatAGA6dOnd9mmJ6qrqyGTyeDp6dnh+w0NDaipqTF5keUpeGqRbmAwSB0B2RKz11rMyclBSUkJioqK2r1XWVkJANBqtSblWq0Wp06d6rTPysrKDtu09tcbV65cwZNPPon58+d3msmzsrLwzDPP9Hob1D+cOf2ebsAJH2QOs34VrqioQHp6Onbu3AkXF5dO68lueFCgEKJdWX+06UxTUxPmzZsHg8GAbdu2dVovMzMT1dXVxldFRUWvtkd9w+n3dKOWFiYy6jmzRmTFxcWoqqpCTEyMsaylpQUffvghtm7diq+++grA1RHW4MGDjXWqqqrajbja0ul07UZf3bXpTFNTE+bOnYuysjL861//6vK8qkqlgkqlMnsb1L84a5FuZA0jsitNLXBx5nU6W2BWIktOTsbx48dNyu677z6MHDkSK1euxNChQ6HT6ZCfn49x48YBABobG1FQUIDs7OxO+42Pj0d+fj4ee+wxY1leXh4SEhLMCc+YxL7++mu8//778PHxMas9SYOJjG5kyVmLlxub8W3VJXzzY+3VP6vq8M2PdZg3PhAPThpqsTio98xKZGq1GhERESZl7u7u8PHxMZZnZGRg/fr1CAsLQ1hYGNavXw83NzfMnz/f2GbhwoXw9/dHVlYWACA9PR2TJ09GdnY2Zs+ejb179+LgwYM4fPiwsU1dXR2++eYb489lZWUoLS2Ft7c3goKC0NzcjF//+tcoKSnBO++8g5aWFuMoz9vbG0ql0sxdQ5ai5GNc6AYD8UyyH2sb8O2PdVcTVVUdvv2xDt9W1eFszRV0tLkAL9d+j4EGRr8/WHPFihWor6/H4sWLcfHiRcTFxSEvLw9qtdpYp7y8HE5O17+8EhISkJOTg1WrVmH16tUIDQ3F7t27ERcXZ6zzySefYMqUKcafly9fDgC455578Oqrr+L06dPYt28fACAqKsokpvfffx9JSUn9/VGpHznLZWjidRG6prcjMoNB4PTFenzzY22bhHUJ3/5Yh58vN5nVV4jvoF7FQJbXq/vI7BXvI5NO7LqDOF/XIHUYZCUOrZiCQG+3Tt+/0tSC764lqNZTgd9W1aHs/CU0NPd97r6TDPhi7QxeI5NYT7+T+31ERtQbn6yaih9+rsfJylqcPFeLr87V4utzV7+k6ptapA6PLKz11GL15SaT0VXrCOv0xcsD+tyywRpXJjEbwkRGVsPf0xX+nq6YMtLPWGYwCFRcvIyvKmvxdVUdvrqW6L778RIaW3jXrL1a+vqnOFtdj/N1jZJs35/Xx2wKExlZNScnGYJ93BHs446U0dfLm1sM+P7CZZw8V9vmVYfvz19CM9fps3nHf6iWdPuBXp2f1iTrw0RGNkkhd8Iwv0EY5jcIvxhz/Z7FxmYDvjt/deT29bm6a6coa1H+08CeiiL7EujNEZktYSIju6JUOGGkzgMjdaYXhq80teCb1lOTVbU4WVmL4lMXUXOlWaJIyZrFD+U9qLaEiYwcgouzHBH+GkT4a4xl979ahH99WSVhVGSNfAcpMX6It9RhkBl4Jyo5rKraK1KHQFZoWrgOTk5cyNqWMJGRwzpXw/vWqL0ZETqpQyAzMZGRQ2oxCFzgDdh0Aw8XBRJCeX3M1jCRkUM6X9fAWYzUTvIoLRextkH8FyOHdK6G18eoPZ5WtE2ctUgOSePqjKVThqGhuQWNzQY0thjQ0GRAw7U/r/7cgp8vN+Grc7VSh0sWMEilQOLwm6QOg3qBiYwcUrCPO347fUS39Y58ex7z/3TUAhGR1G4d6cf1FW0UTy0SdeFyAxcsdhSzIgd3X4msEhMZURcuNXLlD0fg4cLTiraMiYyoC5c4InMIS6YM42lFG8ZERtSFyxyR2b0AL1fce8sQqcOgPuBkD5LM7qJy/OXIKbir5HBVKuCulMNVKYe7UgE3pRxuSsW1966WGd9TyeF2Q5mrcmB+m65rYCKzd09MHwGVgqMxW8ZERpIo+v4nrHrrczS19M9dyU4ywNVZDjfV9SR49c82iVF1vdyYBFVyuDpfTZgdvXeJicyujQ30xO1j9VKHQX3EREYWd+bnejyys6TfkhgAGARwqbEFlxp5TYt6btWsUZDJuECwreM1MrKo97+qwu1b/43zXOeQJDZ9tJaPa7ETHJGRRVxpakHW/v/gL4WnpA6FCM5yGZ6cOUrqMKifMJHRgPviTA0ydn+Kk+fqpA6FCACwIC4YIb7uUodB/YSJjAaMEAIvHy7Dxne/QmOzQepwiAAAOg8XPJ4yXOowqB8xkdGAOFdzBY///RgOf3Ne6lCITDx3RwTULs5Sh0H9iImM+t2B42eRmXscP19ukjoUIqM7o/2h17gieZRW6lConzGRUb+pa2jGM/tO4I3i01KHQmRC66HCprlRUodBA4SJjPrEYBA4/kM19pScxluf/oCaK7yBmKwPp9nbNyYyK1LfePUhjwICQgACVydMXP0TELha2PZnk3rCtBzooA8h0NBsQFOLAY3NBjS1CDS2tD5cUlz9s837jW3+bLrh5/N1jfj01EXUcvULsnJxIUxk9oyJbABcufZk4YuXG/Hz5SZU1zfi4uUm/Hy5CT9fK7t4uRE/1zeh+nITfrrciJ8vN/brShdEdN14JjK7xkTWz85W1yPxvz/gdHMiK+Hp5owRWrXUYdAA4hJV/WywxhV3jw+UOgwiuiY22JvrKdo5JrIBsOTWYXBx5q4lsga8Pmb/+G07APzULrgnfojUYRARgFuG+UodAg0wJrIB8nBiKAapeAmSSEpaDxXC9R5Sh0EDjIlsgHi5K3E/H59OJKnE4TdJHQJZQJ8SWVZWFmQyGTIyMoxlQgisWbMGer0erq6uSEpKwokTJ7rta8+ePQgPD4dKpUJ4eDhyc3NN3v/www9x2223Qa/XQyaT4a233mrXR2+3PVAenDwUGleu6UYklcThflKHQBbQ60RWVFSEHTt2IDIy0qR848aN2LRpE7Zu3YqioiLodDpMmzYNtbW1nfZVWFiI1NRUpKWl4dixY0hLS8PcuXNx9OhRY51Lly5h7Nix2Lp1a6f99GbbA8nDxRn/NXmoJNsmcnRyJxkmhvH6mCOQCSHMvgu3rq4O0dHR2LZtG9atW4eoqChs3rwZQgjo9XpkZGRg5cqVAICGhgZotVpkZ2dj0aJFHfaXmpqKmpoaHDhwwFg2Y8YMeHl5YdeuXe2DlsmQm5uLOXPmGMt6u+22ampqoNFoUF1dDQ+P/jmvfrmxGZOy38eFS4390h8R9UxkgAb7lk6UOgzqg55+J/dqRLZkyRLMmjULU6dONSkvKytDZWUlUlJSjGUqlQqJiYk4cuRIp/0VFhaatAGA6dOnd9nmRr3ZdkNDA2pqakxe/c1NqcAjSaH93i8RdS3Q203qEMhCzE5kOTk5KCkpQVZWVrv3KisrAQBareljErRarfG9jlRWVprdpj+2nZWVBY1GY3wFBg7Mjcy/uTkYOg+XAembiDoW4OUqdQhkIWYlsoqKCqSnp2Pnzp1wcen8i/nGu+iFEN3eWd+bNn3tJzMzE9XV1cZXRUWF2dvrCRdnOZbcOmxA+iaijgV4cUTmKMxKZMXFxaiqqkJMTAwUCgUUCgUKCgqwZcsWKBQK42joxhFQVVVVu5FSWzqdzuw2HfVh7rZVKhU8PDxMXgNl3vhA/oZIZEF6Dc+COAqzEllycjKOHz+O0tJS4ys2NhYLFixAaWkphg4dCp1Oh/z8fGObxsZGFBQUICEhodN+4+PjTdoAQF5eXpdtbhQSEtKrbVuKs9wJy5LDpA6DyGG4KuVSh0AWYtbSE2q1GhERESZl7u7u8PHxMZZnZGRg/fr1CAsLQ1hYGNavXw83NzfMnz/f2GbhwoXw9/c3XmdLT0/H5MmTkZ2djdmzZ2Pv3r04ePAgDh8+bGxTV1eHb775xvhzWVkZSktL4e3tjaCgIOP9bN1tW0q/ig7A9g++xXfnL0kdCpHdc3VmInMU/b6G0ooVK1BfX4/Fixfj4sWLiIuLQ15eHtTq649RKC8vh5PT9cFgQkICcnJysGrVKqxevRqhoaHYvXs34uLijHU++eQTTJkyxfjz8uXLAQD33HMPXn311R5vW0pyJxnSp4YhPadU6lCI7B5HZI6jV/eR2auBuI/sRkIIzHzxEL6slOYmbSJHUfBEEoJ93KUOg/pgQO8jo96TyWR4bNpwqcMgsns8teg4mMgkMH20DpEBGqnDILJrLjy16DCYyCTy2FSOyogGEkdkjoOJTCJTRvphXJCn1GEQ2SWFkwzOcn69OQr+S0vo8WkjpA6ByC5xNOZYmMgkNDHMFxNCvKUOg8ju8PqYY2Eik9hyzmAk6ncKJ/PXaSXbxUQmsZuH+iAh1EfqMIiIbBYTmRV4PIWjMiKi3mIiswIxwd5IHH6T1GEQ2Q2uV+RYmMisBK+VERH1DhOZlRgb6Impo3r+/DUi6pwAh2SOhInMijyeMhy9eCg2EZFDYyKzIqMGe+CXkXqpwyCyebxG5liYyKzMY1PDIOc9MEREPcZEZmWG3jQId47zlzoMIpvGAZljYSKzQulTw6DkgqdERD3Cb0srFODlhtTxgVKHQWSz5Jw15VCYyKzUo7cOg4sz/3mIeoPXmR0LvymtlJ+HC9JuDpY6DCKbxAGZY2Eis2KPJA3DIJVC6jCIbA5HZI6FicyKebsr8eCkEKnDILI5vEbmWJjIrNyDk4bC210pdRhENoV5zLEwkVm5QSoFFieFSh0GkU3hqUXHwkRmA9Lig6HXuEgdBpHNcOKQzKEwkdkAlUKO9KlhUodBZDOYyBwLE5mN+HVMIIbe5C51GEQ2gacWHQsTmY2QO8nw+LQRUodBZBOcmMgcChOZDfnFGB0i/D2kDoPI6smZxxwKE5kNkclkeGL6SKnDILJ6vEbmWJjIbEzi8Jtw81BvqcMgsmoGPlnToTCR2aAVMzgqI+pKfZNB6hDIgpjIbFB0kBemjtJKHQaR1brS1CJ1CGRBTGQ26onpI8CJWUQdq29kInMkTGQ2aoROjdlR/lKHQWSV6jkicyhMZDbssanD4cx5xkTtMJE5lj4lsqysLMhkMmRkZBjLhBBYs2YN9Ho9XF1dkZSUhBMnTnTb1549exAeHg6VSoXw8HDk5ua2q7Nt2zaEhITAxcUFMTExOHTokMn7dXV1WLp0KQICAuDq6opRo0bhpZde6stHtGpBPm64e0KQ1GEQWZ3GZgN+utSI6stNqLnShLqGZlxqaEZ9YwuuNLWgsdmAphYDWgwCgjMcbV6vn9pYVFSEHTt2IDIy0qR848aN2LRpE1599VUMHz4c69atw7Rp0/DVV19BrVZ32FdhYSFSU1Px7LPP4o477kBubi7mzp2Lw4cPIy4uDgCwe/duZGRkYNu2bbjlllvwxz/+ETNnzsQXX3yBoKCrX+aPPfYY3n//fezcuRNDhgxBXl4eFi9eDL1ej9mzZ/f2o1q1R28Nw57i07jEawJEJqKfzTervkx29f4zJxkgg8z0Z9n1n9uVo+N6rT87yQABQIirtwW0/VMIAYMABK79ea1MXOtTKZdBqXC6/pK3/l1ufM9Z7gS5kwzOTk5QyGXGnxXy62UKJxkUciconK6/7yyXQWF8v7Xt1b8rFU4YrHGBv6crFHLrP3EnE734daSurg7R0dHYtm0b1q1bh6ioKGzevBlCCOj1emRkZGDlypUAgIaGBmi1WmRnZ2PRokUd9peamoqamhocOHDAWDZjxgx4eXlh165dAIC4uDhER0ebjLBGjRqFOXPmICsrCwAQERGB1NRUrF692lgnJiYGv/jFL/Dss892+7lqamqg0WhQXV0NDw/bWUHj9/kn8eJ7X0sdBhHZGYWTDIM9XRDk7YYgb/drf7oh2McNQT5u8HBxHtDt9/Q7uVcjsiVLlmDWrFmYOnUq1q1bZywvKytDZWUlUlJSjGUqlQqJiYk4cuRIp4mssLAQjz32mEnZ9OnTsXnzZgBAY2MjiouL8eSTT5rUSUlJwZEjR4w/T5w4Efv27cP9998PvV6PDz74ACdPnsSLL77Y4XYbGhrQ0NBg/LmmpqZnO8DK/NfkoXjt6Cmcr2uUOhQisiPNBoGKn+pR8VM9/o0L7d7XuDoj2McNgd5uGKlVI2mEHyL8PSCz8MoqZo8Zc3JyUFJSYhwFtVVZWQkA0GpN73HSarXG9zpSWVnZZZvz58+jpaWl2363bNmC8PBwBAQEQKlUYsaMGdi2bRsmTpzY4XazsrKg0WiMr8DAwC4+ufVyVynw6K18zAsRWVZ1fRM+O12Nf3x2Fi/kn8RtWw9j4Z8/tngcZiWyiooKpKenY+fOnXBx6fxBjzdmYyFEtxm6J226q7NlyxZ89NFH2LdvH4qLi/HCCy9g8eLFOHjwYIfbzMzMRHV1tfFVUVHRZYzWbH5cEIJ93KQOg4gc3KGvz2N7wbcW3aZZpxaLi4tRVVWFmJgYY1lLSws+/PBDbN26FV999RWAqyOswYMHG+tUVVW1G021pdPp2o3Y2rbx9fWFXC7vsk59fT1+97vfITc3F7NmzQIAREZGorS0FM8//zymTp3abrsqlQoqlcqcXWC1nOVO+G3KCDy661OpQyEiB5f9zy+RPNIPYdqOJ/j1N7NGZMnJyTh+/DhKS0uNr9jYWCxYsAClpaUYOnQodDod8vOvzxZqbGxEQUEBEhISOu03Pj7epA0A5OXlGdsolUrExMS0q5Ofn2+s09TUhKamJjg5mX4kuVwOg8Ex1l37ZeRgRAZopA6DiBycEMD/FZ+22PbMGpGp1WpERESYlLm7u8PHx8dYnpGRgfXr1yMsLAxhYWFYv3493NzcMH/+fGObhQsXwt/f33idLT09HZMnT0Z2djZmz56NvXv34uDBgzh8+LCxzfLly5GWlobY2FjEx8djx44dKC8vx8MPPwwA8PDwQGJiIp544gm4uroiODgYBQUF+Otf/4pNmzb1bu/YGJlMhhXTR+I3Lx+VOhQicnClFT9bbFu9vo+sMytWrEB9fT0WL16MixcvIi4uDnl5eSb3kJWXl5uMnBISEpCTk4NVq1Zh9erVCA0Nxe7du433kAFXp+hfuHABa9euxdmzZxEREYH9+/cjODjYWCcnJweZmZlYsGABfvrpJwQHB+O5554zJjtHMDHMF3Eh3jha9pPUoRCRAztXc8Vi2+rVfWT2ylbvI7vRx2U/Ye4fC6UOg4gcWKC3Kw6tuLVPffT0O9n6b9kms00I8cakMF+pwyAiB9bSYrkxEhOZncqcOQpyPueFiCRisOC5PiYyOxWu98DdE2zzBm8isn2WXNyDicyO/TZlBDzdBnYtNCKijljyfBATmR3zdFMic+ZIqcMgIgdkyfUWmcjsXOr4IKTdHNx9RSIiG8VE5gDW3D4ak4ffJHUYREQDgonMAcidZNg6fxzC/AZJHQoROQhO9qB+5+HijD/fOx4+7kqpQyEiB8BERgMi0NsNf31gAvzU9rHiPxFZL5kF5y0ykTmY0XoN3lycgNCb3KUOhYjsGEdkNKACvNyw55EExAZ7SR0KEVGfMZE5KE83JV57KA6/uTlI6lCIyA7xhmiyCJVCjnVzxuAP86OhVvX7E32IyIHxhmiyqFmRg/HOsomI8LfdR9cQkXXhiIwsLtjHHXseSeCpRiLqH5zsQVJoPdX44rwouCvlUodDRDaMIzKS1Owof+xdOhHDtVwJhIisHxMZdWiY3yDsXTIRd0b7Sx0KEdkgTvYgq+CqlGPT3ChsuHMMVAoeKkTUczy1SFZl3oQg5C6+BUN83KQOhYhsBFf2IKsTrvfA249OxC/G6KQOhYhsgEFYbltMZNRjahdnbFsQg6dvC4ez3JInDojI1pyruWKxbTGRkdnuuyUEf18UD39PV6lDISIrVXulGZcbmy2yLSYy6pVxQV74x7KJmDKCT54moo5dvNxkke0wkVGvebop8cp9E5B15xiu1UhE7Vy81GiR7TCRUZ/dPSEIecsn49aRflKHQkRW5GeOyMiWDNa44s/3jsfvU8fC081Z6nCIyApcvMwRGdmgO8YFIP+xRMyM4DR9Ikf3MxMZ2aqb1Cq89JsYvLQgGr6DVFKHQ0QS4WQPsnkzxwzGweWTcec4rtdI5Ih4apHsgqebEptSo/DKveMxWOMidThEZEGXG1ossh0mMrKIKSP9kPfYZNw9Iciia7ARkXQs9X+diYwsRu3ijKw7x+C1B+MQ5M0FiInsHRMZ2a2EUF+8mzEZ998SAieOzojslqWeScZERpJwVcrx/24LxxsPJyD0JnepwyGiAWCpX1T7lMiysrIgk8mQkZFhLBNCYM2aNdDr9XB1dUVSUhJOnDjRbV979uxBeHg4VCoVwsPDkZub267Otm3bEBISAhcXF8TExODQoUPt6vznP//B7bffDo1GA7VajZtvvhnl5eV9+Zg0gGKCvbA/fRIWJ4VCweEZkV2RWejxmr1OZEVFRdixYwciIyNNyjdu3IhNmzZh69atKCoqgk6nw7Rp01BbW9tpX4WFhUhNTUVaWhqOHTuGtLQ0zJ07F0ePHjXW2b17NzIyMvDUU0/h008/xaRJkzBz5kyTJPXtt99i4sSJGDlyJD744AMcO3YMq1evhosLZ8tZM5VCjhUzRuKtJbdgKEdnRHbDUtfIZEIIsx9/VldXh+joaGzbtg3r1q1DVFQUNm/eDCEE9Ho9MjIysHLlSgBAQ0MDtFotsrOzsWjRog77S01NRU1NDQ4cOGAsmzFjBry8vLBr1y4AQFxcHKKjo/HSSy8Z64waNQpz5sxBVlYWAGDevHlwdnbG3/72N3M/EgCgpqYGGo0G1dXV8PDw6FUf1Dd1Dc347d+P4Z8nKqUOhYj6aGF8MNbOjuh1+55+J/dqRLZkyRLMmjULU6dONSkvKytDZWUlUlJSjGUqlQqJiYk4cuRIp/0VFhaatAGA6dOnG9s0NjaiuLi4XZ2UlBRjHYPBgH/84x8YPnw4pk+fDj8/P8TFxeGtt97qdLsNDQ2oqakxeZG0BqkU2J4WgydnjoScpxqJqAfMTmQ5OTkoKSkxjoLaqqy8+lu0Vqs1Kddqtcb3OlJZWdllm/Pnz6OlpaXLOlVVVairq8OGDRswY8YM5OXl4Y477sCdd96JgoKCDreblZUFjUZjfAUGBnbz6clSHk4Mxd/unwAfd6XUoRCRlTMrkVVUVCA9PR07d+7s8rrTjVMuhRDdTsPsSZuu6hgMBgDA7Nmz8dhjjyEqKgpPPvkkfvnLX2L79u0dbjMzMxPV1dXGV0VFRZcxkmUlDPPF249OxNhAT6lDISIrZlYiKy4uRlVVFWJiYqBQKKBQKFBQUIAtW7ZAoVAYR0w3jr6qqqrajaba0ul0Xbbx9fWFXC7vto5CoUB4eLhJnVGjRnU6a1GlUsHDw8PkRdZF7+mKNxbF4+4JHC0TUcfMSmTJyck4fvw4SktLja/Y2FgsWLAApaWlGDp0KHQ6HfLz841tGhsbUVBQgISEhE77jY+PN2kDAHl5ecY2SqUSMTEx7erk5+eb1Bk/fjy++uorkzonT55EcHCwOR+TrIxS4YSsOyOx/o4xUMp56yMRmTLr+fRqtRoREaYzUNzd3eHj42Msz8jIwPr16xEWFoawsDCsX78ebm5umD9/vrHNwoUL4e/vb7zOlp6ejsmTJyM7OxuzZ8/G3r17cfDgQRw+fNjYZvny5UhLS0NsbCzi4+OxY8cOlJeX4+GHHzbWeeKJJ5CamorJkydjypQp+Oc//4m3334bH3zwgdk7hqzP/LggjNANwiM7S1BV2yB1OERkJcxKZD2xYsUK1NfXY/Hixbh48SLi4uKQl5cHtVptrFNeXg4np+u/WSckJCAnJwerVq3C6tWrERoait27dyMuLs5YJzU1FRcuXMDatWtx9uxZREREYP/+/SajrTvuuAPbt29HVlYWli1bhhEjRmDPnj2YOHFif39MkkhMsDfeeXQiHt5ZjJLyn6UOh4isQK/uI7NXvI/MdjQ2G/D0vhPY9TFXbSGyVlZ9HxmR1K5eNxvD62ZExERGtm1+XBB2L7qZD+0kcmBMZGTzxgV54Z1HJyIh1EfqUIhIAkxkZBd8Bqnwtwfi8EhSKJ9ATeRgmMjIbsidZFg5YyT++JsYqF36fUIuEVkpJjKyOymjddi3dCJG6tTdVyYim8dERnYpxNcduYtvwZwovdShENEAYyIju+WqlGPzvHFYO3s0nOW8cEZkr5jIyO4tjB+C3YviofPgFH0ie8RERg4hOsgL7yybiElhvlKHQkT9jImMHIbvIBX+ev8EPDF9BJ8+TWRHmMjIochkMiyZMgy7HrqZpxqJ7AQTGTmkCSHe2J8+CVNG3CR1KETUR0xk5LC83ZX4873jkTlzJBQ81Uhks5jIyKHJZDIsSgzF3x+Oh7+nq9ThEFEvMJER4eqsxv3LJmFauFbqUIjITExkRNdo3Jzxp4WxWP3LcD7jjMiG8H8r0Q0emBiC/3skHoHePNVIZAuYyIg6EBngiX8sm4SZETqpQyGibjCREXXCw8UZL/0mBs/OHg2lgv9ViKwV/3cSdSMtfghyFydg6E3uUodCRB1gIiPqgdF6Dd55dCJSYwOlDoWIbsBERtRDbkoFsn8diT/Mj4YHn0BNZDWYyIjMNCtyMA5kTMb4IV5Sh0JEYCIj6hV/T1fk/Fc8Hps6nCvpE0mMiYyol+ROMqRPDcPfF93M5a2IJMRERtRHMcHeOJAxCb+MHCx1KEQOiYmMqB94uDhj6/xo/PevI+GulEsdDpFVsNRJdyYyon50V2wg3lk2CZEBGqlDIZKcTGaZVMZERtTPQnzdseeRBCyaPBQW+n9MZJWcmMiIbJez3AmZvxiFv90fBz+1SupwiCRhqQm9TGREA2himC/+mTEZySP9pA6FyOIsdUaCiYxogHm7K/HyvePxzO2joeLiw+RAeGqRyM7ckzAEe5feguHaQVKHQmQRnOxBZIdG6jywb+lEpN0cLHUoRAOO18iI7JSLsxzPzonAy/fEwneQUupwiAaMTZxazMrKgkwmQ0ZGhrFMCIE1a9ZAr9fD1dUVSUlJOHHiRLd97dmzB+Hh4VCpVAgPD0dubm67Otu2bUNISAhcXFwQExODQ4cOddrfokWLIJPJsHnz5t58NKIBlzxKiwPpk5E04iapQyEaEFY/2aOoqAg7duxAZGSkSfnGjRuxadMmbN26FUVFRdDpdJg2bRpqa2s77auwsBCpqalIS0vDsWPHkJaWhrlz5+Lo0aPGOrt370ZGRgaeeuopfPrpp5g0aRJmzpyJ8vLydv299dZbOHr0KPR6fW8/HpFF3KRW4dX7JmDt7NFwceYJErIvVn2NrK6uDgsWLMCf/vQneHldf5SFEAKbN2/GU089hTvvvBMRERH4y1/+gsuXL+P111/vtL/Nmzdj2rRpyMzMxMiRI5GZmYnk5GST0dSmTZvwwAMP4MEHH8SoUaOwefNmBAYG4qWXXjLp64cffsDSpUvx2muvwdnZuTcfj8jiFsYPwdtLJ2LUYA+pQyHqN1Z9jWzJkiWYNWsWpk6dalJeVlaGyspKpKSkGMtUKhUSExNx5MiRTvsrLCw0aQMA06dPN7ZpbGxEcXFxuzopKSkm/RoMBqSlpeGJJ57A6NGje/PRiCQTplVj75Jb8NCkEK4IQnbBUtfIzH7MbU5ODkpKSlBUVNTuvcrKSgCAVqs1KddqtTh16lSnfVZWVnbYprW/8+fPo6Wlpcs6AJCdnQ2FQoFly5b16LM0NDSgoaHB+HNNTU2P2hENFKXCCU/NCkfSCD88/vdjqKy5InVIRL1mlYsGV1RUID09HTt37oSLi0un9W48LyqE6PZcaU/adFWnuLgYL774Il599dUen5fNysqCRqMxvgIDA3vUjmig3TLMF//MmIRZY/hoGLJdThY6t2hWIisuLkZVVRViYmKgUCigUChQUFCALVu2QKFQGEdMbUdJAFBVVdVuNNWWTqfrso2vry/kcnmXdQ4dOoSqqioEBQUZYzt16hQef/xxDBkypMPtZmZmorq62viqqKgwZ3cQDShPNyX+sCAa238TA51H5784Elkrq5y1mJycjOPHj6O0tNT4io2NxYIFC1BaWoqhQ4dCp9MhPz/f2KaxsREFBQVISEjotN/4+HiTNgCQl5dnbKNUKhETE9OuTn5+vrFOWloaPvvsM5PY9Ho9nnjiCbz77rsdblelUsHDw8PkRWRtZkTokL98MtJuDrbYxXOi/mCV18jUajUiIiJMytzd3eHj42Msz8jIwPr16xEWFoawsDCsX78ebm5umD9/vrHNwoUL4e/vj6ysLABAeno6Jk+ejOzsbMyePRt79+7FwYMHcfjwYWOb5cuXIy0tDbGxsYiPj8eOHTtQXl6Ohx9+GADg4+MDHx8fk9icnZ2h0+kwYsQIcz4mkdVRuzjj2TkRuCPaH7978zi+rOz8dhYia2GpX7zMnuzRnRUrVqC+vh6LFy/GxYsXERcXh7y8PKjVamOd8vJyODldHwwmJCQgJycHq1atwurVqxEaGordu3cjLi7OWCc1NRUXLlzA2rVrcfbsWURERGD//v0IDuZSP+Q4ooO88PajE7Hjw++w5b2v0dBskDokok7JLDTdQyaEEBbZkg2oqamBRqNBdXU1TzOS1fv+/CU89dZx/PubC1KHQtShVbNG4cFJQ3vdvqffyVxKgMhGDfF1x2sP3owX7hoLb3eu2UjWxybWWiQi6f0qJgAHlyfirpgA3khNVsWqV/YgIuvi7a7Ef981FnuX3ILoIE+pwyECYKX3kRGRdYsM8MSeRxKwOTWK956R5KxyZQ8isn4ymQxzxvnjX79NxJIpoVAp+N+cpCF3ssyxxyOcyE65KRV4YvpIHFyeiBmjdVKHQw7IWc5Ti0TUDwK93bA9LQa7HroZYwM9pQ6HHIjSQmcDmMiIHER8qA/2LrkF238TgzC/QVKHQw5AKbdMiun3lT2IyLrNiNBh+mgtahuacfFSI/71ZRVeP1qOr6vqpA6N7IwzExkRDRSZTAYPF2d4uDjjvltCcN8tISj6/ifsOlqOfxw/y6WvqF84W+jUIhMZEQEAxg/xxvgh3vh/t4VjT8kP2PVxOb7hKI36gKcWiUgSnm5KPDAxBA9MDMHHZT/h9aOncODzSo7SyGxKhRU+xoWIHMuEEG9MCPHGmsuN+L/i09j1cTm+/fGS1GGRjbDUNTLOWiSibnm6KfHgpKE4uDwRG38VCd9BXKSYusfp90RkdWQyGeaOD8S/fpuExUmhcHWWSx0SWTFLXSNjIiMis3m4OGPFjJEoeCIJC+KCoLDUMudkUwa5WObqFRMZEfWan4cLnrtjDPKXJ2JW5GA+RoZMaFydLbIdJjIi6rMQX3f8YX409i2ZiElhvlKHQ1bAxdkJKoVlTj0zkRFRvxkToMHfHojDaw/GITJAI3U4JCFPV8tNCGIiI6J+d8swX+xdcgv+MD8a/p6uUodDErDUaUWAiYyIBohMJsOsyMHIe2wy7k0YYrHH3pN1YCIjIrvhrlJgze2j8cbDCRjGVfcdhgcTGRHZm5hgL/xj2UQsnTKMsxsdgKcbExkR2SGVQo7fTh+BzalRFrtZlqTBU4tEZNdmR/nj1fvGQ63icq/2iomMiOxewjBf7F4Uz3Ub7RQTGRE5hHC9B169bwLUFlrKiCyH18iIyGFE+Gvw8j3j4eLMryN74unGG6KJyIFMCPHGtgXRXHzYjvipVRbbFhMZEVmFW0dq8cLcsZyabye0Hi4W2xYTGRFZjdlR/ngkMVTqMKiPlHIneLvz1CIROajl04YjKtBT6jCoD/w8LHdaEWAiIyIro5A7Ycu8cRjEe8xs1nCt2qLbYyIjIqsT5OOGZ+eMljoM6qUIvYdFt8dERkRW6Y5xAbh9rF7qMKgXRvtb9ll0TGREZLXW3D7aopMGqH+MYSIjIrrK212J1b8cJXUYZIYAL1foLfww1T4lsqysLMhkMmRkZBjLhBBYs2YN9Ho9XF1dkZSUhBMnTnTb1549exAeHg6VSoXw8HDk5ua2q7Nt2zaEhITAxcUFMTExOHTokPG9pqYmrFy5EmPGjIG7uzv0ej0WLlyIM2fO9OUjEpHE7hgXgMThN0kdBvXQzAidxbfZ60RWVFSEHTt2IDIy0qR848aN2LRpE7Zu3YqioiLodDpMmzYNtbW1nfZVWFiI1NRUpKWl4dixY0hLS8PcuXNx9OhRY53du3cjIyMDTz31FD799FNMmjQJM2fORHl5OQDg8uXLKCkpwerVq1FSUoI333wTJ0+exO23397bj0hEVuK5OyLgrpRLHQb1wMwxgy2+TZkQQpjbqK6uDtHR0di2bRvWrVuHqKgobN68GUII6PV6ZGRkYOXKlQCAhoYGaLVaZGdnY9GiRR32l5qaipqaGhw4cMBYNmPGDHh5eWHXrl0AgLi4OERHR+Oll14y1hk1ahTmzJmDrKysDvstKirChAkTcOrUKQQFBXX7uWpqaqDRaFBdXQ0PD8vOuiGirm374Bts/OdXUodBXQjwcsWhFVMg66flWXr6ndyrEdmSJUswa9YsTJ061aS8rKwMlZWVSElJMZapVCokJibiyJEjnfZXWFho0gYApk+fbmzT2NiI4uLidnVSUlK67Le6uhoymQyenp4dvt/Q0ICamhqTFxFZp/tvCcFgjeWWPSLz3XdLSL8lMXOYnchycnJQUlLS4SiosrISAKDVak3KtVqt8b2OVFZWdtnm/PnzaGlpMavfK1eu4Mknn8T8+fM7zeRZWVnQaDTGV2BgYKcxEpG0XJzleGzqcKnDoE54uCgwb7w036FmJbKKigqkp6dj586dcHHp/DejGzOyEKLbLN2TNj3tt6mpCfPmzYPBYMC2bds63WZmZiaqq6uNr4qKii5jJCJp/SomAMO1g6QOgzqwMH4I3CVajcWsRFZcXIyqqirExMRAoVBAoVCgoKAAW7ZsgUKhMI6YbhwlVVVVtRtNtaXT6bps4+vrC7lc3qN+m5qaMHfuXJSVlSE/P7/L86oqlQoeHh4mLyKyXnInGVbOGCl1GHSDIT5uWHrrMMm2b1YiS05OxvHjx1FaWmp8xcbGYsGCBSgtLcXQoUOh0+mQn59vbNPY2IiCggIkJCR02m98fLxJGwDIy8sztlEqlYiJiWlXJz8/36Tf1iT29ddf4+DBg/Dx8THn4xGRDUgepUVciLfUYdA1Mhmw4VeRcHGWblapWeNAtVqNiIgIkzJ3d3f4+PgYyzMyMrB+/XqEhYUhLCwM69evh5ubG+bPn29ss3DhQvj7+xuvs6Wnp2Py5MnIzs7G7NmzsXfvXhw8eBCHDx82tlm+fDnS0tIQGxuL+Ph47NixA+Xl5Xj44YcBAM3Nzfj1r3+NkpISvPPOO2hpaTGO4Ly9vaFUcnUAInuR+YtRmPOHf0sdBgF49NYw3DxU2kFDv5/QXLFiBerr67F48WJcvHgRcXFxyMvLg1p9fTXk8vJyODldHwwmJCQgJycHq1atwurVqxEaGordu3cjLi7OWCc1NRUXLlzA2rVrcfbsWURERGD//v0IDg4GAJw+fRr79u0DAERFRZnE9P777yMpKam/PyoRSSQq0BOzxgzGP46flToUhzZjtA6PTQ2TOoze3Udmr3gfGZHt+P78JUz7fQGaWvgVJoWEUB/8+d7xA3pKcUDvIyMiktoQX3fcPaH7hQ6o/00I8cbL9wxsEjMHExkR2az05DCo+QBOi7p9rB5/uW8CXK1oyTAmMiKyWT6DVHg4KVTqMByCwkmG1b8Mx5a7x1lVEgOYyIjIxj0wMQR6Ll01oHwHKfG3B+LwwMQQqUPpEBMZEdk0F2c5Hk8ZIXUYdis22AtvPzoR8aHWe18uExkR2bw7o/0xWs+Zxv3JSQY8eusw7F4Uj8Eayz4o01xMZERk82QyGZ6axSdJ9xethwo7H4zD4ykjIHey/Gr25mIiIyK7kBDqi1tH+kkdhs1LCPXBgfTJSAj1lTqUHmMiIyK78btfjIKz3PpHENZqxmgdXrlvPLzdbWtJPyYyIrIbw/wG4aFJQ6UOwyZNHaXFHxZEQ6Wwrqn1PcFERkR2ZVlyGAK8rHtygrXx93TFC3eNtYnrYR1hIiMiu+LiLMczt4+WOgyb4SQDNs+LgsbNWepQeo2JjIjsTvIoLVLCO3+YL1334KShGD/Etp/vxkRGRHZpze2j4WZlSylZm2F+g7B82nCpw+gzJjIiskt6T1ekJ0v/rCxr9tycCKtZwb4vmMiIyG49MDEEI7Tq7is6oGnhWsRJ/GTn/sJERkR2SyF3wtO3h0sdhtVROMmQOXOk1GH0GyYyIrJrCaG+mBmhkzoMqzIxzBdDbxokdRj9homMiOze734xCioFv+5aJQ2/SeoQ+hX/ZYnI7gV6u+G/JnPFj1ZeNrYEVXeYyIjIISxOGobBfAAnANjFTMW2mMiIyCG4KuV40o4mOPSFvZ1mta9PQ0TUhdlR/ogN9pI6DMm5ckRGRGS71tw+Gja6Nm6/CfF1lzqEfsVERkQOJcJfgzvGBUgdhmSigzzh52Ff1wqZyIjI4SxPGQ6lnV0n6gm1SoF1c8ZIHUa/c7x/SSJyeP6erki7OVjqMCxquHYQdv3XzQjXe0gdSr9TSB0AEZEUlk4Zhr8XVaC2oVnqUAaMs1yGCH8NbovU4zc3B9vtKJSJjIgckpe7Eo8mD8P6/V9KHUqvuSnl8Pd0hb+XK/w9XaH3dEXAtb/7e7lCq3aBkwPMbGEiIyKH9V+TQ/Hl2Vq8+ekPUofSIY2rs0li8vd0RYCXm7HM3lbo6C0mMiJyaBt+FQlnuRP+XlwBISy7bd9BSvhfS0wB15LV1STlBn8vVwxS8Su6J7iXiMihKRVOyP51JO6KDcCqtz7Hl5W1/da3t7sSAV6uCGxNVl7XR1QBXm5w5ROs+wUTGRERgNgh3njn0Yl498Q5/OXI9/j4+5/a1ZHJAD+1CsE+7hji44ab1CqoXZwxSKWA2uXqa5DKGV5uzvD3coWbkl+xlsC9TER0jULuhFmRgzErcjAqfrqMH36uR92VZigVTtB6uCDIm6Moa8RERkTUgUBvNwR6u0kdBvWAfd5UQEREDqNPiSwrKwsymQwZGRnGMiEE1qxZA71eD1dXVyQlJeHEiRPd9rVnzx6Eh4dDpVIhPDwcubm57eps27YNISEhcHFxQUxMDA4dOmTyfm+3TUREtqvXiayoqAg7duxAZGSkSfnGjRuxadMmbN26FUVFRdDpdJg2bRpqazufCVRYWIjU1FSkpaXh2LFjSEtLw9y5c3H06FFjnd27dyMjIwNPPfUUPv30U0yaNAkzZ85EeXl5n7ZNREQ2TvRCbW2tCAsLE/n5+SIxMVGkp6cLIYQwGAxCp9OJDRs2GOteuXJFaDQasX379k77mzt3rpgxY4ZJ2fTp08W8efOMP0+YMEE8/PDDJnVGjhwpnnzyyT5tu63q6moBQFRXV/eoPhERDZyefif3akS2ZMkSzJo1C1OnTjUpLysrQ2VlJVJSUoxlKpUKiYmJOHLkSKf9FRYWmrQBgOnTpxvbNDY2ori4uF2dlJQUY53ebLuhoQE1NTUmLyIisi1mz1rMyclBSUkJioqK2r1XWVkJANBqtSblWq0Wp06d6rTPysrKDtu09nf+/Hm0tLR0Wac3287KysIzzzzTaVxERGT9zBqRVVRUID09HTt37oSLS+cPZpPJTBepFEK0K+tNm/6q0yozMxPV1dXGV0VFRZcxEhGR9TErkRUXF6OqqgoxMTFQKBRQKBQoKCjAli1boFAojKOh1tFRq6qqqnYjpbZ0Ol2XbXx9fSGXy7uso9PpzN62SqWCh4eHyYuIiGyLWYksOTkZx48fR2lpqfEVGxuLBQsWoLS0FEOHDoVOp0N+fr6xTWNjIwoKCpCQkNBpv/Hx8SZtACAvL8/YRqlUIiYmpl2d/Px8Y52QkJBebZuIiGxcX2eVtJ21KIQQGzZsEBqNRrz55pvi+PHj4u677xaDBw8WNTU1xjppaWnG2YZCCPHvf/9byOVysWHDBvGf//xHbNiwQSgUCvHRRx8Z6+Tk5AhnZ2fx8ssviy+++EJkZGQId3d38f3335u17a5w1iIRkfXo6Xdyvy9RtWLFCtTX12Px4sW4ePEi4uLikJeXB7VabaxTXl4OJ6frg8GEhATk5ORg1apVWL16NUJDQ7F7927ExcUZ66SmpuLChQtYu3Ytzp49i4iICOzfvx/BwcFmbZuIiOyLTAhLP4HHetXU1ECj0aC6uprXy4iIJNbT72SutUhERDaNiYyIiGwaH+PSRutZVq7wQUQkvdbv4u6ugDGRtdG6uHBgYKDEkRARUava2lpoNJpO3+dkjzYMBgPOnDkDtVqN2tpaBAYGoqKiwmYmftTU1DBmC7HFuBmzZTDm/iOEQG1tLfR6vclM9xtxRNaGk5MTAgICAFxf6soWV/xgzJZji3EzZstgzP2jq5FYK072ICIim8ZERkRENo2JrBMqlQpPP/00VCqV1KH0GGO2HFuMmzFbBmO2PE72ICIim8YRGRER2TQmMiIismlMZEREZNOYyIiIyKbZZSK7ePEi0tLSoNFooNFokJaWhp9//rnLNkIIrFmzBnq9Hq6urkhKSsKJEyeM73///feQyWQdvt544w0AwAcffNBpnaKiImNfHb2/fft2yeIGgCFDhrR7/8knnzTZVnl5OW677Ta4u7vD19cXy5Ytw7lz5ySJ+fvvv8cDDzyAkJAQuLq6IjQ0FE8//TQaGxtNttVRHy+88IJk+7kn/8aW2s+tCgsLceutt8Ld3R2enp5ISkpCfX09gL4d0wOxn3sSM2Bdx3NPYu7L8TxQ3x093de9PaZv/Gx9MjDP9ZTWjBkzREREhDhy5Ig4cuSIiIiIEL/85S+7bLNhwwahVqvFnj17xPHjx0VqaqrJ06Wbm5vF2bNnTV7PPPOMcHd3F7W1tUIIIRoaGtrVefDBB8WQIUOEwWAwbguAeOWVV0zqXb58WbK4hRAiODhYrF271qRe2/ebm5tFRESEmDJliigpKRH5+flCr9eLoKAgSWI+cOCAuPfee8W7774rvv32W7F3717h5+cnHn/8cZNtdbSvp02bJtl+7u7f2JL7WQghjhw5Ijw8PERWVpb4/PPPxcmTJ8Ubb7whrly5IoTo2zE9EPu5JzELYV3Hc09i7svxPFDfHT3d1709ppcuXdplfOawu0T2xRdfCADio48+MpYVFhYKAOLLL7/ssI3BYBA6nU5s2LDBWHblyhWh0WjE9u3bO91WVFSUuP/++zt9v7GxUfj5+Ym1a9ealAMQubm5VhV3cHCw+P3vf99pm/379wsnJyfxww8/GMuef/55AUC89957ksR8o40bN4qQkBCTshv3tZT7uSfbtvR+jouLE6tWrer0M9yop8f0QO7nnsRsbcezuftZiJ4dz0JIu697e0zv2rVLqFQqUV1d3cUe6Dm7O7VYWFgIjUaDuLg4Y9nNN98MjUaDI0eOdNimrKwMlZWVSElJMZapVCokJiZ22qa4uBilpaV44IEHOo1l3759OH/+PO6999527y1duhS+vr4YP348tm/fjn//+9+Sx52dnQ0fHx9ERUXhueeeMxn6FxYWIiIiAnq93limVCoBXF+XUoqY26quroa3t3e78rb7OisrS7L93JNj05L7uaqqCkePHoWfnx8SEhKg1WqRmJiIw4cPd9gn0PNjeqD2szkxW8vx3Jv9DPTseB7I746exN3bY3r69OloaGhAcXFxl/ugp+xu0eDKykr4+fm1K/fz80NlZWWnbQBAq9WalGu1Wpw6darDNi+//DJGjRqFhISETmN5+eWXMX369HaPhXn22WeRnJwMV1dXvPfee3j88ceRkJAgadzp6emIjo6Gl5cXPv74Y2RmZqKsrAz/+7//a9zWjdtpfezNjfFJsa+//fZb/M///A9eeOEFk/Ib93VmZmaHi6JaIuaeHJuW3M/fffcdAGDNmjV4/vnnERUVhb/+9a9ITk7G559/jrCwsA4/V0+O6YHazz2N2ZqO597s554ezwP53dGTuHt7THt5eUGpVHYan7lsJpGtWbMGzzzzTJd1Wi8+t/2NqpUQosPytm58v7M29fX1eP3117F69epO+zp9+jTeffdd/OpXv+qwj7Zt09PTsX37dgQHB0sW92OPPWb8+5tvvonz58/j5Zdfxssvv9zhtm680C9FzK3OnDmD8ePHo6amBg899BAeeughk/db2xYVFWHq1Kl499132/VhqZjb1m97TP/ud7/D7373u3b1BnI/GwwGAMCiRYtw3333AQDGjRuH9957D3/+85+RlZVl0vb06dM4cOBAh/0CltnPPY3Zmo5nc/fzmTNnMGPGDISGhnZ5PAMD+93R07h78n3b2+/knrKZRLZ06VLMmzevyzpDhgzBZ599hnPnzrV778cff2z3W0ErnU4H4OpvDoMHDzaWV1VVddjm//7v/3D58mUsXLiw01heeeUV+Pj44MUXX8TatWu7jPvMmTN48cUXO/ztxNJxA1f39ZQpU5CUlIScnByMHTsWW7ZswXvvvYe9e/cCuLqvCwsLAbT/rc6SMZ85cwZTpkzB1KlTsWbNmi6fWTRkyBCMHTsWBw4cwLlz50y2Z4mYdTqdybHZekyPHz8emZmZuPPOOy26n1vLw8PDTdqOGjUK5eXl7fpsPabff/99ODs7d7jd1pgHaj+bGzMg/fFsTsytx3N8fDyef/55/PTTTx1ut239gfru6EncNx7THW1bp9Ph6NGjJu9fvHgRTU1NncZntn650mZFWi8+Hj161Fj20Ucf9ejCZ3Z2trGsoaGh04v5iYmJ4le/+lWnMRgMBhESEtJuxlFn/ud//kcolUrJ427r7bffFgDEqVOnhBDXL9ieOXPGWOeFF14QAMS//vUvSWI+ffq0CAsLE/PmzRPNzc09+lxPPfWUACAOHTpk8Zh7cmxacj8bDAah1+vbXcyPiooSmZmZ7foz55geqP1sTsxtSXk89zTm3hzPA/nd0ZO4e3tM5+Tk9OtkD7tLZEJcnQ4aGRkpCgsLRWFhoRgzZky7qagjRowQb775pvHnDRs2CI1GI958801x/Phxcffdd7ebiiqEEF9//bWQyWTiwIEDnW7/4MGDAoD44osv2r23b98+sWPHDnH8+HHxzTffiD/96U/Cw8NDLFu2TLK4jxw5IjZt2iQ+/fRT8d1334ndu3cLvV4vbr/9dmOd1im0ycnJoqSkRBw8eFAEBASIoKAgSWL+4YcfxLBhw8Stt94qTp8+bTIdubt9HRwcLNnx0d2/saX38+9//3vh4eEh3njjDfH111+LVatWCRcXF/HNN9+Y9N2bY3qg9nN3MVvj8dxdzH05ngfyu6Mnx0dvj2lOv+/GhQsXxIIFC4RarRZqtVosWLBAXLx40aQOrt2P0cpgMIinn35a6HQ6oVKpxOTJk8Xx48fb9Z2ZmSkCAgJES0tLp9u/++67RUJCQofvHThwQERFRYlBgwYJNzc3ERERITZv3iyampoki7u4uFjExcUJjUYjXFxcxIgRI8TTTz8tLl26ZFLv1KlTYtasWcLV1VV4e3uLpUuXijNnzkgS8yuvvCIAdPjqbl+fO3dOsuOjJ//Glt7PWVlZIiAgQLi5uYn4+HiTUVSr3hzTA7mfu4rZWo/nrmLuy/E80N8d3R0fvT2m296L1ld8jAsREdk0u7uPjIiIHAsTGRER2TQmMiIismlMZEREZNOYyIiIyKYxkRERkU1jIiMiIpvGREZEZKOee+45JCQkwM3NDZ6enj1qI3rwZOikpKR2T6Juu9ZtT59o3Z2srCyMHz8earUafn5+mDNnDr766iuz+gCYyIiIbFZjYyPuuusuPPLIIz1us3HjRmzatAlbt25FUVERdDodpk2bZnyMTauHHnoIZ8+eNb7++Mc/Gt/78ssvYTAY8Mc//hEnTpzA73//e2zfvt3kCQ49UVBQgCVLluCjjz5Cfn4+mpubkZKSgkuXLpnVj10uUUVE5EheeeUVodFouq3X0ydDJyYmivT0dLNi6OiJ1idOnBAzZ84U7u7uws/PT/zmN78RP/74Y6d9VFVVCQCioKDArG1zREZE5CDMedr5a6+9Bl9fX4wePRq//e1v243YbnTjE63Pnj2LxMREREVF4ZNPPsE///lPnDt3DnPnzu2yDwAdPhm7KzbzPDIiIuqbnj7tfMGCBQgJCYFOp8Pnn3+OzMxMHDt2DPn5+R3229ETrV966SVER0dj/fr1xrI///nPCAwMxMmTJzF8+HCTPoQQWL58OSZOnIiIiAizPhcTGRGRFWn75PDOFBUVITY2ttfb6O4p2G2fSh0REYGwsDDExsaipKQE0dHRJm1bn2h911134cEHHzSWFxcX4/3338egQYPabf/bb79tl8iWLl2Kzz77DIcPHzb78zCRERFZkdYnh3dlyJAhverb3Kedt4qOjoazszO+/vprk0TW9onWO3bsMGljMBhw2223ITs7u11/bbcNAI8++ij27duHDz/8EAEBAWZ/LiYyIiIr4uvrC19f3wHpu/V0YX5+PsaNGwfg6szHgoKCDhNOqxMnTqCpqckkAf3www+YMmUKYmJi8Morr8DJyXTKRXR0NPbs2YMhQ4ZAoeg41Qgh8OijjyI3NxcffPABQkJCevW5ONmDiMhGlZeXo7S0FOXl5WhpaUFpaSlKS0tRV1dnrDNy5Ejk5uYCuHpKMSMjA+vXr0dubi4+//xz3HvvvXBzc8P8+fMBXD3tt3btWnzyySf4/vvvsX//ftx1110YN24cbrnlFgBXR2JJSUkIDAzE888/jx9//BGVlZXGa3AAsGTJEvz000+4++678fHHH+O7775DXl4e7r//frS0tBjr7Ny5E6+//jrUarWxj/r6evN2hFlzHImIyGrcc889HT5V+v333zfWgZlPhi4vLxeTJ08W3t7eQqlUitDQULFs2TJx4cIFY52ePNFaCCFOnjwp7rjjDuHp6SlcXV3FyJEjRUZGhjAYDMbYOnq1jbcn+IRoIiKyaTy1SERENo2JjIiIbBoTGRER2TQmMiIismlMZEREZNOYyIiIyKYxkRERkU1jIiMiIpvGREZERDaNiYyIiGwaExkREdk0JjIiIrJp/x8Z6dz/HV6NTwAAAABJRU5ErkJggg==", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "ucb_gdf.plot()" + ] + }, + { + "cell_type": "markdown", + "id": "eaafd823-c1f1-46f2-96e2-486a405e9ef5", + "metadata": {}, + "source": [ + "We have a map of the UTTC Campus!\n", + "\n", + "> **Warning**\n", + ">\n", + "> The Open Street Maps (OSM) database is not always complete. For\n", + "> example, try searching for UTTC with the `{'building': True}`, and\n", + "> compare it to the [map of the UTTC campus on their\n", + "> website](https://uttc.edu/wp-content/uploads/2022/12/CAMPUS-MAP-2022-23.pdf).\n", + "> What do you notice?\n", + "\n", + "### Create an interactive map\n", + "\n", + "There are lots of different ways to create maps and plots in Python.\n", + "Here, we’re going to use a tool called `'hvplot'` and `'geoviews'` to\n", + "create an interactive map, including the online `'EsriImagery'` tile\n", + "source basemap." + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "id": "8176d4d3", + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/opt/conda/lib/python3.11/site-packages/dask/dataframe/__init__.py:31: FutureWarning: \n", + "Dask dataframe query planning is disabled because dask-expr is not installed.\n", + "\n", + "You can install it with `pip install dask[dataframe]` or `conda install dask`.\n", + "This will raise in a future version.\n", + "\n", + " warnings.warn(msg, FutureWarning)\n", + "WARNING:bokeh.core.validation.check:W-1005 (FIXED_SIZING_MODE): 'fixed' sizing mode requires width and height to be set: figure(id='p1043', ...)\n" + ] + }, + { + "data": {}, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "application/vnd.holoviews_exec.v0+json": "", + "text/html": [ + "
\n", + "
\n", + "
\n", + "" + ], + "text/plain": [ + ":Overlay\n", + " .WMTS.I :WMTS [Longitude,Latitude]\n", + " .Polygons.I :Polygons [Longitude,Latitude]" + ] + }, + "execution_count": 4, + "metadata": { + "application/vnd.holoviews_exec.v0+json": { + "id": "p1101" + } + }, + "output_type": "execute_result" + } + ], + "source": [ + "# Plot UTTC boundary\n", + "ucb_map = ucb_gdf.reset_index().hvplot(\n", + " # Givethe map a descriptive title\n", + " title=\"University of Colorado, Boulder - Main Campus\",\n", + " # Add a basemap\n", + " geo=True, tiles='EsriImagery',\n", + " # Change the colors\n", + " fill_color='white', fill_alpha=0.2,\n", + " line_color='skyblue', line_width=5,\n", + " # Change the image size\n", + " frame_width=400, frame_height=400)\n", + "\n", + "# Save the map as a file to put on the web\n", + "hv.save(ucb_map, 'ucb.html')\n", + "\n", + "# Display the map\n", + "ucb_map" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "base", + "language": "python", + "name": "base" + }, + "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.11.6" + } + }, + "nbformat": 4, + "nbformat_minor": 5 +} From 9b53453f72af7a5645609b44fb6a2a020a3490d3 Mon Sep 17 00:00:00 2001 From: G Paul Bailey Date: Thu, 19 Sep 2024 09:48:59 -0600 Subject: [PATCH 2/2] Create chatgpt-first-map given the following code, use hvplot to output an html file of the university area but make it a gold outline of 4 pixels with no fill and old main point to map using esri imagery as a base --- chatgpt-first-map | 59 +++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 59 insertions(+) create mode 100644 chatgpt-first-map diff --git a/chatgpt-first-map b/chatgpt-first-map new file mode 100644 index 0000000..2b17a9e --- /dev/null +++ b/chatgpt-first-map @@ -0,0 +1,59 @@ +# modify first map using chatgpt +# given the following code, use hvplot to output an html file of the university area but make it a gold outline of 4 pixels with no fill and old main point to map using esri imagery as a base +# Work with vector data +import geopandas as gpd +from shapely.geometry import Point + +# Save maps and plots to files +import holoviews as hv +import hvplot.pandas +hv.extension('bokeh') # Ensure Bokeh extension is loaded for interactive plots + +# Search for locations by name - this might take a moment +from osmnx import features as osm + +# Search for universities in Boulder, CO +universities_gdf = osm.features_from_address( + 'Boulder, CO, United States', + {'amenity': ['university']} +) + +# Filter the GeoDataFrame to only include the University of Colorado +ucb_gdf = universities_gdf[universities_gdf['name'].str.contains('University of Colorado', case=False, na=False)] + +# Define the coordinates for Old Main, Boulder, CO +old_main_gdf = gpd.GeoDataFrame( + {'name': ['Old Main']}, + geometry=[Point(-105.2705, 40.0095)], # Correct longitude and latitude + crs="EPSG:4326" # Coordinate reference system (WGS84) +) + +# Plotting the University of Colorado area with hvplot +# Styling the university boundary as a gold outline with no fill +ucb_plot = ucb_gdf.hvplot( + geo=True, + tiles='ESRI', + line_width=4, + color='gold', + fill_alpha=0, # No fill + hover_cols=['name'], # Show name on hover + legend=False +) + +# Plotting the Old Main point +old_main_plot = old_main_gdf.hvplot.points( + geo=True, + color='gold', + size=100, + marker='*', + hover_cols=['name'] +) + +# Combine the plots +combined_plot = ucb_plot * old_main_plot + +# Save the plot as an HTML file +hvplot.save(combined_plot, 'university_of_colorado_map.html') + +# Display the plot +combined_plot