Azimuth.ipynb

{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "name": "Azimuth.ipynb",
      "provenance": [],
      "collapsed_sections": [],
      "include_colab_link": true
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    },
    "widgets": {
      "application/vnd.jupyter.widget-state+json": {
        "a252f72c51754c19b56345fd62e8621d": {
          "model_module": "@jupyter-widgets/controls",
          "model_name": "ButtonModel",
          "state": {
            "_view_name": "ButtonView",
            "style": "IPY_MODEL_e5e18e85d0a1447d898ad373c514a7e2",
            "_dom_classes": [],
            "description": "Calculate!",
            "_model_name": "ButtonModel",
            "button_style": "",
            "_view_module": "@jupyter-widgets/controls",
            "_model_module_version": "1.5.0",
            "tooltip": "",
            "_view_count": null,
            "disabled": false,
            "_view_module_version": "1.5.0",
            "layout": "IPY_MODEL_9d56bd267d76484da2fc96a7bc76be6e",
            "_model_module": "@jupyter-widgets/controls",
            "icon": ""
          }
        },
        "e5e18e85d0a1447d898ad373c514a7e2": {
          "model_module": "@jupyter-widgets/controls",
          "model_name": "ButtonStyleModel",
          "state": {
            "_view_name": "StyleView",
            "_model_name": "ButtonStyleModel",
            "_view_module": "@jupyter-widgets/base",
            "_model_module_version": "1.5.0",
            "_view_count": null,
            "button_color": null,
            "font_weight": "",
            "_view_module_version": "1.2.0",
            "_model_module": "@jupyter-widgets/controls"
          }
        },
        "9d56bd267d76484da2fc96a7bc76be6e": {
          "model_module": "@jupyter-widgets/base",
          "model_name": "LayoutModel",
          "state": {
            "_view_name": "LayoutView",
            "grid_template_rows": null,
            "right": null,
            "justify_content": null,
            "_view_module": "@jupyter-widgets/base",
            "overflow": null,
            "_model_module_version": "1.2.0",
            "_view_count": null,
            "flex_flow": null,
            "width": null,
            "min_width": null,
            "border": null,
            "align_items": null,
            "bottom": null,
            "_model_module": "@jupyter-widgets/base",
            "top": null,
            "grid_column": null,
            "overflow_y": null,
            "overflow_x": null,
            "grid_auto_flow": null,
            "grid_area": null,
            "grid_template_columns": null,
            "flex": null,
            "_model_name": "LayoutModel",
            "justify_items": null,
            "grid_row": null,
            "max_height": null,
            "align_content": null,
            "visibility": null,
            "align_self": null,
            "height": null,
            "min_height": null,
            "padding": null,
            "grid_auto_rows": null,
            "grid_gap": null,
            "max_width": null,
            "order": null,
            "_view_module_version": "1.2.0",
            "grid_template_areas": null,
            "object_position": null,
            "object_fit": null,
            "grid_auto_columns": null,
            "margin": null,
            "display": null,
            "left": null
          }
        },
        "59c5228d07634c0faac1818ad593b5ed": {
          "model_module": "@jupyter-widgets/output",
          "model_name": "OutputModel",
          "state": {
            "_view_name": "OutputView",
            "msg_id": "",
            "_dom_classes": [],
            "_model_name": "OutputModel",
            "outputs": [
              {
                "output_type": "display_data",
                "metadata": {
                  "tags": []
                },
                "text/html": "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe src=\"about:blank\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" data-html=PCFET0NUWVBFIGh0bWw+CjxoZWFkPiAgICAKICAgIDxtZXRhIGh0dHAtZXF1aXY9ImNvbnRlbnQtdHlwZSIgY29udGVudD0idGV4dC9odG1sOyBjaGFyc2V0PVVURi04IiAvPgogICAgPHNjcmlwdD5MX1BSRUZFUl9DQU5WQVM9ZmFsc2U7IExfTk9fVE9VQ0g9ZmFsc2U7IExfRElTQUJMRV8zRD1mYWxzZTs8L3NjcmlwdD4KICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2Nkbi5qc2RlbGl2ci5uZXQvbnBtL2xlYWZsZXRAMS40LjAvZGlzdC9sZWFmbGV0LmpzIj48L3NjcmlwdD4KICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2NvZGUuanF1ZXJ5LmNvbS9qcXVlcnktMS4xMi40Lm1pbi5qcyI+PC9zY3JpcHQ+CiAgICA8c2NyaXB0IHNyYz0iaHR0cHM6Ly9tYXhjZG4uYm9vdHN0cmFwY2RuLmNvbS9ib290c3RyYXAvMy4yLjAvanMvYm9vdHN0cmFwLm1pbi5qcyI+PC9zY3JpcHQ+CiAgICA8c2NyaXB0IHNyYz0iaHR0cHM6Ly9jZG5qcy5jbG91ZGZsYXJlLmNvbS9hamF4L2xpYnMvTGVhZmxldC5hd2Vzb21lLW1hcmtlcnMvMi4wLjIvbGVhZmxldC5hd2Vzb21lLW1hcmtlcnMuanMiPjwvc2NyaXB0PgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL2Nkbi5qc2RlbGl2ci5uZXQvbnBtL2xlYWZsZXRAMS40LjAvZGlzdC9sZWFmbGV0LmNzcyIvPgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL21heGNkbi5ib290c3RyYXBjZG4uY29tL2Jvb3RzdHJhcC8zLjIuMC9jc3MvYm9vdHN0cmFwLm1pbi5jc3MiLz4KICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9tYXhjZG4uYm9vdHN0cmFwY2RuLmNvbS9ib290c3RyYXAvMy4yLjAvY3NzL2Jvb3RzdHJhcC10aGVtZS5taW4uY3NzIi8+CiAgICA8bGluayByZWw9InN0eWxlc2hlZXQiIGhyZWY9Imh0dHBzOi8vbWF4Y2RuLmJvb3RzdHJhcGNkbi5jb20vZm9udC1hd2Vzb21lLzQuNi4zL2Nzcy9mb250LWF3ZXNvbWUubWluLmNzcyIvPgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL2NkbmpzLmNsb3VkZmxhcmUuY29tL2FqYXgvbGlicy9MZWFmbGV0LmF3ZXNvbWUtbWFya2Vycy8yLjAuMi9sZWFmbGV0LmF3ZXNvbWUtbWFya2Vycy5jc3MiLz4KICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9yYXdjZG4uZ2l0aGFjay5jb20vcHl0aG9uLXZpc3VhbGl6YXRpb24vZm9saXVtL21hc3Rlci9mb2xpdW0vdGVtcGxhdGVzL2xlYWZsZXQuYXdlc29tZS5yb3RhdGUuY3NzIi8+CiAgICA8c3R5bGU+aHRtbCwgYm9keSB7d2lkdGg6IDEwMCU7aGVpZ2h0OiAxMDAlO21hcmdpbjogMDtwYWRkaW5nOiAwO308L3N0eWxlPgogICAgPHN0eWxlPiNtYXAge3Bvc2l0aW9uOmFic29sdXRlO3RvcDowO2JvdHRvbTowO3JpZ2h0OjA7bGVmdDowO308L3N0eWxlPgogICAgCiAgICA8bWV0YSBuYW1lPSJ2aWV3cG9ydCIgY29udGVudD0id2lkdGg9ZGV2aWNlLXdpZHRoLAogICAgICAgIGluaXRpYWwtc2NhbGU9MS4wLCBtYXhpbXVtLXNjYWxlPTEuMCwgdXNlci1zY2FsYWJsZT1ubyIgLz4KICAgIDxzdHlsZT4jbWFwXzdlMjJiYzQ1NmE4ZTRmODFiNzY5MmFkMGZmNGIxNzJjIHsKICAgICAgICBwb3NpdGlvbjogcmVsYXRpdmU7CiAgICAgICAgd2lkdGg6IDEwMC4wJTsKICAgICAgICBoZWlnaHQ6IDEwMC4wJTsKICAgICAgICBsZWZ0OiAwLjAlOwogICAgICAgIHRvcDogMC4wJTsKICAgICAgICB9CiAgICA8L3N0eWxlPgo8L2hlYWQ+Cjxib2R5PiAgICAKICAgIAogICAgPGRpdiBjbGFzcz0iZm9saXVtLW1hcCIgaWQ9Im1hcF83ZTIyYmM0NTZhOGU0ZjgxYjc2OTJhZDBmZjRiMTcyYyIgPjwvZGl2Pgo8L2JvZHk+CjxzY3JpcHQ+ICAgIAogICAgCiAgICAKICAgICAgICB2YXIgYm91bmRzID0gbnVsbDsKICAgIAoKICAgIHZhciBtYXBfN2UyMmJjNDU2YThlNGY4MWI3NjkyYWQwZmY0YjE3MmMgPSBMLm1hcCgKICAgICAgICAnbWFwXzdlMjJiYzQ1NmE4ZTRmODFiNzY5MmFkMGZmNGIxNzJjJywgewogICAgICAgIGNlbnRlcjogWy0yMC4wNjY2NjQzNDcwMDcyODMsIC00NC4wMjU0NTQyMjU4NDQ5NV0sCiAgICAgICAgem9vbTogMTMsCiAgICAgICAgbWF4Qm91bmRzOiBib3VuZHMsCiAgICAgICAgbGF5ZXJzOiBbXSwKICAgICAgICB3b3JsZENvcHlKdW1wOiBmYWxzZSwKICAgICAgICBjcnM6IEwuQ1JTLkVQU0czODU3LAogICAgICAgIHpvb21Db250cm9sOiB0cnVlLAogICAgICAgIH0pOwoKCiAgICAKICAgIHZhciB0aWxlX2xheWVyXzQzYjM5YjgzNTIyMzRjMmViMTk2YjQwM2U5ZDZjYWE1ID0gTC50aWxlTGF5ZXIoCiAgICAgICAgJ2h0dHBzOi8vc3RhbWVuLXRpbGVzLXtzfS5hLnNzbC5mYXN0bHkubmV0L3RlcnJhaW4ve3p9L3t4fS97eX0uanBnJywKICAgICAgICB7CiAgICAgICAgImF0dHJpYnV0aW9uIjogbnVsbCwKICAgICAgICAiZGV0ZWN0UmV0aW5hIjogZmFsc2UsCiAgICAgICAgIm1heE5hdGl2ZVpvb20iOiAxOCwKICAgICAgICAibWF4Wm9vbSI6IDE4LAogICAgICAgICJtaW5ab29tIjogMCwKICAgICAgICAibm9XcmFwIjogZmFsc2UsCiAgICAgICAgIm9wYWNpdHkiOiAxLAogICAgICAgICJzdWJkb21haW5zIjogImFiYyIsCiAgICAgICAgInRtcyI6IGZhbHNlCn0pLmFkZFRvKG1hcF83ZTIyYmM0NTZhOGU0ZjgxYjc2OTJhZDBmZjRiMTcyYyk7CiAgICAKICAgICAgICB2YXIgbWFya2VyX2ViMDE2OGQwMzU4ODQyYWZiMDNiYmE3NDQ3MjBlOTdmID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFstMjAuMDU4NzIwOTIyMTAxMjQ4LCAtNDQuMDEwMzg5MDQyOTg3MjldLAogICAgICAgICAgICB7CiAgICAgICAgICAgICAgICBpY29uOiBuZXcgTC5JY29uLkRlZmF1bHQoKSwKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgKS5hZGRUbyhtYXBfN2UyMmJjNDU2YThlNGY4MWI3NjkyYWQwZmY0YjE3MmMpOwogICAgICAgIAogICAgCgogICAgICAgICAgICAgICAgdmFyIGljb25fNWIxN2JlZjMwNWE5NGE2NGI2NDJhMDdmMGZiNjNmNjYgPSBMLkF3ZXNvbWVNYXJrZXJzLmljb24oewogICAgICAgICAgICAgICAgICAgIGljb246ICdpbmZvLXNpZ24nLAogICAgICAgICAgICAgICAgICAgIGljb25Db2xvcjogJ3doaXRlJywKICAgICAgICAgICAgICAgICAgICBtYXJrZXJDb2xvcjogJ2JsdWUnLAogICAgICAgICAgICAgICAgICAgIHByZWZpeDogJ2dseXBoaWNvbicsCiAgICAgICAgICAgICAgICAgICAgZXh0cmFDbGFzc2VzOiAnZmEtcm90YXRlLTAnCiAgICAgICAgICAgICAgICAgICAgfSk7CiAgICAgICAgICAgICAgICBtYXJrZXJfZWIwMTY4ZDAzNTg4NDJhZmIwM2JiYTc0NDcyMGU5N2Yuc2V0SWNvbihpY29uXzViMTdiZWYzMDVhOTRhNjRiNjQyYTA3ZjBmYjYzZjY2KTsKICAgICAgICAgICAgCiAgICAKICAgICAgICAgICAgdmFyIHBvcHVwXzdiMmQyNzNmYTYzNTQ1YmM5Yzc1YmNlOGJlODNlMjhlID0gTC5wb3B1cCh7bWF4V2lkdGg6ICcxMDAlJwogICAgICAgICAgICAKICAgICAgICAgICAgfSk7CgogICAgICAgICAgICAKICAgICAgICAgICAgICAgIHZhciBodG1sXzg5ZGM0YjNhOWU0ZTQ4ZDg5MzJkYmI4OGRiODY4MDA4ID0gJChgPGRpdiBpZD0iaHRtbF84OWRjNGIzYTllNGU0OGQ4OTMyZGJiODhkYjg2ODAwOCIgc3R5bGU9IndpZHRoOiAxMDAuMCU7IGhlaWdodDogMTAwLjAlOyI+PGNlbnRlcj5Qb250byBGZXJydWdpbm9zbzwvY2VudGVyPjxicj48Y2VudGVyPi0yMC4wNi0tNDQuMDE8L2NlbnRlcj48L2Rpdj5gKVswXTsKICAgICAgICAgICAgICAgIHBvcHVwXzdiMmQyNzNmYTYzNTQ1YmM5Yzc1YmNlOGJlODNlMjhlLnNldENvbnRlbnQoaHRtbF84OWRjNGIzYTllNGU0OGQ4OTMyZGJiODhkYjg2ODAwOCk7CiAgICAgICAgICAgIAoKICAgICAgICAgICAgbWFya2VyX2ViMDE2OGQwMzU4ODQyYWZiMDNiYmE3NDQ3MjBlOTdmLmJpbmRQb3B1cChwb3B1cF83YjJkMjczZmE2MzU0NWJjOWM3NWJjZThiZTgzZTI4ZSkKICAgICAgICAgICAgOwoKICAgICAgICAgICAgCiAgICAgICAgCiAgICAKICAgICAgICB2YXIgbWFya2VyX2UyZjZlMDMzZGI5NTQ1MTFiYjJmNmYxMTU2YWZlNmYyID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFstMjAuMDYyNTIzODY3NTUwMTI3LCAtNDQuMDM4MDEwMjg0NDkwMjNdLAogICAgICAgICAgICB7CiAgICAgICAgICAgICAgICBpY29uOiBuZXcgTC5JY29uLkRlZmF1bHQoKSwKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgKS5hZGRUbyhtYXBfN2UyMmJjNDU2YThlNGY4MWI3NjkyYWQwZmY0YjE3MmMpOwogICAgICAgIAogICAgCgogICAgICAgICAgICAgICAgdmFyIGljb25fNGZkNzUwOGE2MmQwNGVlMmFhZDNjYjAzMTU0NDdiMjIgPSBMLkF3ZXNvbWVNYXJrZXJzLmljb24oewogICAgICAgICAgICAgICAgICAgIGljb246ICdpbmZvLXNpZ24nLAogICAgICAgICAgICAgICAgICAgIGljb25Db2xvcjogJ3doaXRlJywKICAgICAgICAgICAgICAgICAgICBtYXJrZXJDb2xvcjogJ2JsdWUnLAogICAgICAgICAgICAgICAgICAgIHByZWZpeDogJ2dseXBoaWNvbicsCiAgICAgICAgICAgICAgICAgICAgZXh0cmFDbGFzc2VzOiAnZmEtcm90YXRlLTAnCiAgICAgICAgICAgICAgICAgICAgfSk7CiAgICAgICAgICAgICAgICBtYXJrZXJfZTJmNmUwMzNkYjk1NDUxMWJiMmY2ZjExNTZhZmU2ZjIuc2V0SWNvbihpY29uXzRmZDc1MDhhNjJkMDRlZTJhYWQzY2IwMzE1NDQ3YjIyKTsKICAgICAgICAgICAgCiAgICAKICAgICAgICAgICAgdmFyIHBvcHVwX2EzMjE4NTc0Nzg0YTQzNjFiY2UyY2Q4NGI4M2NjZGFmID0gTC5wb3B1cCh7bWF4V2lkdGg6ICcxMDAlJwogICAgICAgICAgICAKICAgICAgICAgICAgfSk7CgogICAgICAgICAgICAKICAgICAgICAgICAgICAgIHZhciBodG1sXzU2Yjk2ZWE5MzQyYjRjNGJiZGY1NDQ3NDJhN2NiZmIxID0gJChgPGRpdiBpZD0iaHRtbF81NmI5NmVhOTM0MmI0YzRiYmRmNTQ0NzQyYTdjYmZiMSIgc3R5bGU9IndpZHRoOiAxMDAuMCU7IGhlaWdodDogMTAwLjAlOyI+PGNlbnRlcj5Qb3N0byAzOTwvY2VudGVyPjxicj48Y2VudGVyPi0yMC4wNi0tNDQuMDQ8L2NlbnRlcj48L2Rpdj5gKVswXTsKICAgICAgICAgICAgICAgIHBvcHVwX2EzMjE4NTc0Nzg0YTQzNjFiY2UyY2Q4NGI4M2NjZGFmLnNldENvbnRlbnQoaHRtbF81NmI5NmVhOTM0MmI0YzRiYmRmNTQ0NzQyYTdjYmZiMSk7CiAgICAgICAgICAgIAoKICAgICAgICAgICAgbWFya2VyX2UyZjZlMDMzZGI5NTQ1MTFiYjJmNmYxMTU2YWZlNmYyLmJpbmRQb3B1cChwb3B1cF9hMzIxODU3NDc4NGE0MzYxYmNlMmNkODRiODNjY2RhZikKICAgICAgICAgICAgOwoKICAgICAgICAgICAgCiAgICAgICAgCiAgICAKICAgICAgICB2YXIgbWFya2VyX2ExMTc1NWRhYWRmODRkMWI4OTgxNTVmMTdjNTFjM2RhID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFstMjAuMDc4NzQ4MjUxMzcwNDc3LCAtNDQuMDI3OTYzMzUwMDU3MzM0XSwKICAgICAgICAgICAgewogICAgICAgICAgICAgICAgaWNvbjogbmV3IEwuSWNvbi5EZWZhdWx0KCksCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICkuYWRkVG8obWFwXzdlMjJiYzQ1NmE4ZTRmODFiNzY5MmFkMGZmNGIxNzJjKTsKICAgICAgICAKICAgIAoKICAgICAgICAgICAgICAgIHZhciBpY29uX2RiNjBhNGVmNWE0ODQwOWRhZWQ4YjM3YzEyZjQzNTY1ID0gTC5Bd2Vzb21lTWFya2Vycy5pY29uKHsKICAgICAgICAgICAgICAgICAgICBpY29uOiAnaW5mby1zaWduJywKICAgICAgICAgICAgICAgICAgICBpY29uQ29sb3I6ICd3aGl0ZScsCiAgICAgICAgICAgICAgICAgICAgbWFya2VyQ29sb3I6ICdyZWQnLAogICAgICAgICAgICAgICAgICAgIHByZWZpeDogJ2dseXBoaWNvbicsCiAgICAgICAgICAgICAgICAgICAgZXh0cmFDbGFzc2VzOiAnZmEtcm90YXRlLTAnCiAgICAgICAgICAgICAgICAgICAgfSk7CiAgICAgICAgICAgICAgICBtYXJrZXJfYTExNzU1ZGFhZGY4NGQxYjg5ODE1NWYxN2M1MWMzZGEuc2V0SWNvbihpY29uX2RiNjBhNGVmNWE0ODQwOWRhZWQ4YjM3YzEyZjQzNTY1KTsKICAgICAgICAgICAgCiAgICAKICAgICAgICAgICAgdmFyIHBvcHVwXzJjZmVhODRkMjhlODQzNDBhNzc2ZTViYmQwNjE2Y2MxID0gTC5wb3B1cCh7bWF4V2lkdGg6ICcxMDAlJwogICAgICAgICAgICAKICAgICAgICAgICAgfSk7CgogICAgICAgICAgICAKICAgICAgICAgICAgICAgIHZhciBodG1sXzI0MTc0MThlNmFhYjQwMWFiMDg5YmQ1ZmEyYjM2ZDAyID0gJChgPGRpdiBpZD0iaHRtbF8yNDE3NDE4ZTZhYWI0MDFhYjA4OWJkNWZhMmIzNmQwMiIgc3R5bGU9IndpZHRoOiAxMDAuMCU7IGhlaWdodDogMTAwLjAlOyI+PGNlbnRlcj5JbnRlcnNlw6fDo288L2NlbnRlcj48YnI+PGNlbnRlcj4tMjAuMDgtLTQ0LjAzPC9jZW50ZXI+PC9kaXY+YClbMF07CiAgICAgICAgICAgICAgICBwb3B1cF8yY2ZlYTg0ZDI4ZTg0MzQwYTc3NmU1YmJkMDYxNmNjMS5zZXRDb250ZW50KGh0bWxfMjQxNzQxOGU2YWFiNDAxYWIwODliZDVmYTJiMzZkMDIpOwogICAgICAgICAgICAKCiAgICAgICAgICAgIG1hcmtlcl9hMTE3NTVkYWFkZjg0ZDFiODk4MTU1ZjE3YzUxYzNkYS5iaW5kUG9wdXAocG9wdXBfMmNmZWE4NGQyOGU4NDM0MGE3NzZlNWJiZDA2MTZjYzEpCiAgICAgICAgICAgIDsKCiAgICAgICAgICAgIAogICAgICAgIAogICAgCiAgICAgICAgICAgICAgICB2YXIgcG9seV9saW5lXzlmMGJhNTMwMzAxZTQ4NzY5MmZkYmY5YzA2ZmRmZmU5ID0gTC5wb2x5bGluZSgKICAgICAgICAgICAgICAgICAgICBbWy0yMC4wNTg3MjA5MjIxMDEyNDgsIC00NC4wMTAzODkwNDI5ODcyOV0sIFstMjAuMDYyNTIzODY3NTUwMTI3LCAtNDQuMDM4MDEwMjg0NDkwMjNdXSwKICAgICAgICAgICAgICAgICAgICB7CiAgImJ1YmJsaW5nTW91c2VFdmVudHMiOiB0cnVlLAogICJjb2xvciI6ICJibHVlIiwKICAiZGFzaEFycmF5IjogbnVsbCwKICAiZGFzaE9mZnNldCI6IG51bGwsCiAgImZpbGwiOiBmYWxzZSwKICAiZmlsbENvbG9yIjogImJsdWUiLAogICJmaWxsT3BhY2l0eSI6IDAuMiwKICAiZmlsbFJ1bGUiOiAiZXZlbm9kZCIsCiAgImxpbmVDYXAiOiAicm91bmQiLAogICJsaW5lSm9pbiI6ICJyb3VuZCIsCiAgIm5vQ2xpcCI6IGZhbHNlLAogICJvcGFjaXR5IjogMSwKICAic21vb3RoRmFjdG9yIjogMS4wLAogICJzdHJva2UiOiB0cnVlLAogICJ3ZWlnaHQiOiAyLjUKfQogICAgICAgICAgICAgICAgICAgICkKICAgICAgICAgICAgICAgICAgICAuYWRkVG8obWFwXzdlMjJiYzQ1NmE4ZTRmODFiNzY5MmFkMGZmNGIxNzJjKTsKICAgICAgICAgICAgCiAgICAKICAgICAgICAgICAgICAgIHZhciBwb2x5X2xpbmVfYjcwMjY5MGFhYmZmNDNkNmE0MjNhMDVlODA4ZTAzOGQgPSBMLnBvbHlsaW5lKAogICAgICAgICAgICAgICAgICAgIFtbLTIwLjA1ODcyMDkyMjEwMTI0OCwgLTQ0LjAxMDM4OTA0Mjk4NzI5XSwgWy0yMC4wNzg3NDgyNTEzNzA0NzcsIC00NC4wMjc5NjMzNTAwNTczMzRdXSwKICAgICAgICAgICAgICAgICAgICB7CiAgImJ1YmJsaW5nTW91c2VFdmVudHMiOiB0cnVlLAogICJjb2xvciI6ICJncmVlbiIsCiAgImRhc2hBcnJheSI6IG51bGwsCiAgImRhc2hPZmZzZXQiOiBudWxsLAogICJmaWxsIjogZmFsc2UsCiAgImZpbGxDb2xvciI6ICJncmVlbiIsCiAgImZpbGxPcGFjaXR5IjogMC4yLAogICJmaWxsUnVsZSI6ICJldmVub2RkIiwKICAibGluZUNhcCI6ICJyb3VuZCIsCiAgImxpbmVKb2luIjogInJvdW5kIiwKICAibm9DbGlwIjogZmFsc2UsCiAgIm9wYWNpdHkiOiAxLAogICJzbW9vdGhGYWN0b3IiOiAxLjAsCiAgInN0cm9rZSI6IHRydWUsCiAgIndlaWdodCI6IDIuNQp9CiAgICAgICAgICAgICAgICAgICAgKQogICAgICAgICAgICAgICAgICAgIC5hZGRUbyhtYXBfN2UyMmJjNDU2YThlNGY4MWI3NjkyYWQwZmY0YjE3MmMpOwogICAgICAgICAgICAKICAgIAogICAgICAgICAgICAgICAgdmFyIHBvbHlfbGluZV8wYjExNzM3ZmM4NTM0M2UxODc3M2M1YTJhZjlhN2NmOCA9IEwucG9seWxpbmUoCiAgICAgICAgICAgICAgICAgICAgW1stMjAuMDYyNTIzODY3NTUwMTI3LCAtNDQuMDM4MDEwMjg0NDkwMjNdLCBbLTIwLjA3ODc0ODI1MTM3MDQ3NywgLTQ0LjAyNzk2MzM1MDA1NzMzNF1dLAogICAgICAgICAgICAgICAgICAgIHsKICAiYnViYmxpbmdNb3VzZUV2ZW50cyI6IHRydWUsCiAgImNvbG9yIjogImdyZWVuIiwKICAiZGFzaEFycmF5IjogbnVsbCwKICAiZGFzaE9mZnNldCI6IG51bGwsCiAgImZpbGwiOiBmYWxzZSwKICAiZmlsbENvbG9yIjogImdyZWVuIiwKICAiZmlsbE9wYWNpdHkiOiAwLjIsCiAgImZpbGxSdWxlIjogImV2ZW5vZGQiLAogICJsaW5lQ2FwIjogInJvdW5kIiwKICAibGluZUpvaW4iOiAicm91bmQiLAogICJub0NsaXAiOiBmYWxzZSwKICAib3BhY2l0eSI6IDEsCiAgInNtb290aEZhY3RvciI6IDEuMCwKICAic3Ryb2tlIjogdHJ1ZSwKICAid2VpZ2h0IjogMi41Cn0KICAgICAgICAgICAgICAgICAgICApCiAgICAgICAgICAgICAgICAgICAgLmFkZFRvKG1hcF83ZTIyYmM0NTZhOGU0ZjgxYjc2OTJhZDBmZjRiMTcyYyk7CiAgICAgICAgICAgIAo8L3NjcmlwdD4= onload=\"this.contentDocument.open();this.contentDocument.write(atob(this.getAttribute('data-html')));this.contentDocument.close();\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>",
                "text/plain": "<folium.folium.Map at 0x7fc6b1e38f60>"
              }
            ],
            "_view_module": "@jupyter-widgets/output",
            "_model_module_version": "1.0.0",
            "_view_count": null,
            "_view_module_version": "1.0.0",
            "layout": "IPY_MODEL_a19b2113063e48d282e9fcd80969ca5f",
            "_model_module": "@jupyter-widgets/output"
          }
        },
        "a19b2113063e48d282e9fcd80969ca5f": {
          "model_module": "@jupyter-widgets/base",
          "model_name": "LayoutModel",
          "state": {
            "_view_name": "LayoutView",
            "grid_template_rows": null,
            "right": null,
            "justify_content": null,
            "_view_module": "@jupyter-widgets/base",
            "overflow": null,
            "_model_module_version": "1.2.0",
            "_view_count": null,
            "flex_flow": null,
            "width": null,
            "min_width": null,
            "border": null,
            "align_items": null,
            "bottom": null,
            "_model_module": "@jupyter-widgets/base",
            "top": null,
            "grid_column": null,
            "overflow_y": null,
            "overflow_x": null,
            "grid_auto_flow": null,
            "grid_area": null,
            "grid_template_columns": null,
            "flex": null,
            "_model_name": "LayoutModel",
            "justify_items": null,
            "grid_row": null,
            "max_height": null,
            "align_content": null,
            "visibility": null,
            "align_self": null,
            "height": null,
            "min_height": null,
            "padding": null,
            "grid_auto_rows": null,
            "grid_gap": null,
            "max_width": null,
            "order": null,
            "_view_module_version": "1.2.0",
            "grid_template_areas": null,
            "object_position": null,
            "object_fit": null,
            "grid_auto_columns": null,
            "margin": null,
            "display": null,
            "left": null
          }
        }
      }
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/assismauro/CursoPython/blob/master/Azimuth.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "YjKSN-T2KiNL",
        "colab_type": "text"
      },
      "source": [
        "# Azimuth\n",
        "\n",
        "This Colab Notebook calculates the crossing point of two lines of sight, plot the corresponding azimuths and show the crossing point distances.\n",
        "\n",
        "To run:\n",
        "\n",
        " * Click \"Run\" button (small arrow) left to the next cell.\n",
        " * Fill data\n",
        " * Click \"Run\"button again\n",
        " * click in \"Calculate!\" button below the form cell.\n",
        " "
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "CeIo2Gem9lb2",
        "colab_type": "code",
        "cellView": "form",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 1000,
          "referenced_widgets": [
            "a252f72c51754c19b56345fd62e8621d",
            "e5e18e85d0a1447d898ad373c514a7e2",
            "9d56bd267d76484da2fc96a7bc76be6e",
            "59c5228d07634c0faac1818ad593b5ed",
            "a19b2113063e48d282e9fcd80969ca5f"
          ]
        },
        "outputId": "2ba6b7fc-f3c3-4706-ec66-eea07abd4bea"
      },
      "source": [
        "!pip install simplekml\n",
        "import ipywidgets as widgets\n",
        "from IPython.display import display\n",
        "from shapely.geometry import Point, LineString\n",
        "import numpy as np\n",
        "import matplotlib.pyplot as plt\n",
        "from matplotlib.pyplot import figure\n",
        "import math\n",
        "import folium\n",
        "import simplekml\n",
        "\n",
        "#@title {run: 'auto'} ## Dados de Entrada\n",
        "\n",
        "#@markdown #### Torre A \n",
        "nomeA = \"Ponto Ferruginoso\" #@param {type:\"string\"}\n",
        "xA = 603484.00 #@param {type:\"number\"}\n",
        "yA = 7781714.00 #@param {type:\"number\"}\n",
        "azA = 220 #@param {type:\"number\"}\n",
        " \n",
        "#@markdown #### Torre B \n",
        "nomeB = \"Posto 39\" #@param {type:\"string\"}\n",
        "xB = 600593.00 #@param {type:\"number\"}\n",
        "yB = 7781310.00 #@param {type:\"number\"}\n",
        "azB = 150 #@param {type:\"number\"}\n",
        " \n",
        "button = widgets.Button(description=\"Calculate!\")\n",
        "output = widgets.Output()\n",
        " \n",
        "# trigonometry\n",
        "def sincos(azim):\n",
        "    az = np.radians(azim)\n",
        "    sinaz = np.sin(az)\n",
        "    cosaz = np.cos(az)\n",
        "    return sinaz,cosaz\n",
        " \n",
        "def calcIntersection(xA,yA,azA,xB,yB,azB): \n",
        "    pointA = Point(xA,yA)\n",
        "    sinA, cosA = sincos(azA)\n",
        "    pointB = Point(xB,yB)\n",
        "    sinB, cosB = sincos(azB)\n",
        "    length = 10000\n",
        "    endlineA = Point(pointA.x+(length*sinA), pointA.y+(length*cosA))\n",
        "    endlineB = Point(pointB.x+(length*sinB), pointB.y+(length*cosB))\n",
        "    line1 = LineString([pointA,endlineA])\n",
        "    line2 = LineString([pointB,endlineB])\n",
        "    intersection = LineString([pointB,endlineB]).intersection(LineString([pointA,endlineA]))\n",
        "    return  ((pointA.distance(pointB),\n",
        "             pointA.distance(intersection),\n",
        "             pointB.distance(intersection)),\n",
        "            intersection)\n",
        " \n",
        "def plotResult(p1,p2,p3,nomeA,NomeB,distances):\n",
        "     figure(num=None, figsize=(12, 9), dpi=80, facecolor='w', edgecolor='k')\n",
        "     x_values = [p1[0], p3[0]]\n",
        "     y_values = [p1[1], p3[1]]\n",
        "     plt.plot(x_values, y_values,label='Visada '+nomeA)\n",
        "     x_values = [p2[0], p3[0]]\n",
        "     y_values = [p2[1], p3[1]]\n",
        "     plt.plot(x_values, y_values,label='Visada '+nomeB)\n",
        "     title = 'X: {:10.2f}, Y: {:10.2f}\\nDistância entre torres: {:10.2f}m\\n'.format(\n",
        "             p3[0],p3[1],distances[0])+\\\n",
        "             'Distância {0} à interseção: {1:10.2f}m\\n'.format(nomeA,distances[1])+\\\n",
        "             'Distância {0} à interseção: {1:10.2f}m\\n'.format(nomeB,distances[2])\n",
        "     plt.title(title)\n",
        "     plt.legend()\n",
        "     plt.show()\n",
        "\n",
        "# converte utm para latlong\n",
        "def utmToLatLng(zone, easting, northing, northernHemisphere=True):\n",
        "    if not northernHemisphere:\n",
        "        northing = 10000000 - northing\n",
        "\n",
        "    a = 6378137\n",
        "    e = 0.081819191\n",
        "    e1sq = 0.006739497\n",
        "    k0 = 0.9996\n",
        "\n",
        "    arc = northing / k0\n",
        "    mu = arc / (a * (1 - math.pow(e, 2) / 4.0 - 3 * math.pow(e, 4) / 64.0 - 5 * math.pow(e, 6) / 256.0))\n",
        "\n",
        "    ei = (1 - math.pow((1 - e * e), (1 / 2.0))) / (1 + math.pow((1 - e * e), (1 / 2.0)))\n",
        "\n",
        "    ca = 3 * ei / 2 - 27 * math.pow(ei, 3) / 32.0\n",
        "\n",
        "    cb = 21 * math.pow(ei, 2) / 16 - 55 * math.pow(ei, 4) / 32\n",
        "    cc = 151 * math.pow(ei, 3) / 96\n",
        "    cd = 1097 * math.pow(ei, 4) / 512\n",
        "    phi1 = mu + ca * math.sin(2 * mu) + cb * math.sin(4 * mu) + cc * math.sin(6 * mu) + cd * math.sin(8 * mu)\n",
        "\n",
        "    n0 = a / math.pow((1 - math.pow((e * math.sin(phi1)), 2)), (1 / 2.0))\n",
        "\n",
        "    r0 = a * (1 - e * e) / math.pow((1 - math.pow((e * math.sin(phi1)), 2)), (3 / 2.0))\n",
        "    fact1 = n0 * math.tan(phi1) / r0\n",
        "\n",
        "    _a1 = 500000 - easting\n",
        "    dd0 = _a1 / (n0 * k0)\n",
        "    fact2 = dd0 * dd0 / 2\n",
        "\n",
        "    t0 = math.pow(math.tan(phi1), 2)\n",
        "    Q0 = e1sq * math.pow(math.cos(phi1), 2)\n",
        "    fact3 = (5 + 3 * t0 + 10 * Q0 - 4 * Q0 * Q0 - 9 * e1sq) * math.pow(dd0, 4) / 24\n",
        "\n",
        "    fact4 = (61 + 90 * t0 + 298 * Q0 + 45 * t0 * t0 - 252 * e1sq - 3 * Q0 * Q0) * math.pow(dd0, 6) / 720\n",
        "\n",
        "    lof1 = _a1 / (n0 * k0)\n",
        "    lof2 = (1 + 2 * t0 + Q0) * math.pow(dd0, 3) / 6.0\n",
        "    lof3 = (5 - 2 * Q0 + 28 * t0 - 3 * math.pow(Q0, 2) + 8 * e1sq + 24 * math.pow(t0, 2)) * math.pow(dd0, 5) / 120\n",
        "    _a2 = (lof1 - lof2 + lof3) / math.cos(phi1)\n",
        "    _a3 = _a2 * 180 / math.pi\n",
        "\n",
        "    latitude = 180 * (phi1 - fact1 * (fact2 + fact3 + fact4)) / math.pi\n",
        "\n",
        "    if not northernHemisphere:\n",
        "        latitude = -latitude\n",
        "\n",
        "    longitude = ((zone > 0) and (6 * zone - 183.0) or 3.0) - _a3\n",
        "\n",
        "    return (latitude, longitude)\n",
        "\n",
        "def plotMapFolium(z, p1, p2, p3, nomeA, nomeB, distances, northernHemisphere=False):\n",
        "    lat, lng = utmToLatLng(z, p1[0], p1[1], northernHemisphere=False)\n",
        "    g1 = (lat, lng)\n",
        "    lat, lng = utmToLatLng(z, p2[0], p2[1], northernHemisphere=False)\n",
        "    g2 = (lat, lng)\n",
        "    lat, lng = utmToLatLng(z, p3[0], p3[1], northernHemisphere=False)\n",
        "    g3 = (lat, lng)\n",
        "\n",
        "    points = [g1, g2, g3]\n",
        "    #print(points)\n",
        "\n",
        "    # Load map centred on average coordinates\n",
        "    ave_lat = sum(p[0] for p in points)/len(points)\n",
        "    ave_lng = sum(p[1] for p in points)/len(points)\n",
        "    #print(ave_lat, ave_lng)\n",
        "\n",
        "    title = 'Intersection coordinates (X: {:10.2f}, Y: {:10.2f})\\nDistance between towers: {:10.2f}m\\n'.format(\n",
        "             p3[0],p3[1],distances[0])+\\\n",
        "             'Distance {0} to intersection: {1:10.2f}m\\n'.format(nomeA,distances[1])+\\\n",
        "             'Distance {0} to intersection: {1:10.2f}m\\n'.format(nomeB,distances[2])\n",
        "    \n",
        "    mapa = folium.Map(\n",
        "        location=[ave_lat, ave_lng],\n",
        "        tiles='Stamen Terrain',\n",
        "        zoom_start=13\n",
        "    )\n",
        "    #print(title)\n",
        "    folium.Marker(\n",
        "        location=g1,\n",
        "        popup=\"<center>{}</center><br><center>{:.2f}-{:.2f}</center>\".format(nomeA, g1[0], g1[1]),\n",
        "        icon=folium.Icon(color='blue')\n",
        "    ).add_to(mapa)\n",
        "    folium.Marker(\n",
        "        location=g2,\n",
        "        popup=\"<center>{}</center><br><center>{:.2f}-{:.2f}</center>\".format(nomeB, g2[0], g2[1]),\n",
        "        icon=folium.Icon(color='blue')\n",
        "    ).add_to(mapa)\n",
        "    folium.Marker(\n",
        "        location=g3,\n",
        "        popup=\"<center>Interseção</center><br><center>{:.2f}-{:.2f}</center>\".format(g3[0], g3[1]),\n",
        "        icon=folium.Icon(color='red')\n",
        "    ).add_to(mapa)\n",
        "\n",
        "    folium.PolyLine([g1,g2], color=\"blue\", weight=2.5, opacity=1).add_to(mapa)\n",
        "    folium.PolyLine([g1,g3], color=\"green\", weight=2.5, opacity=1).add_to(mapa)\n",
        "    folium.PolyLine([g2,g3], color=\"green\", weight=2.5, opacity=1).add_to(mapa)\n",
        "        \n",
        "    mapa.save(\"/content/mapa2.html\")\n",
        "    return(mapa)\n",
        "\n",
        "def create_kml(z, p1, p2, p3, nomeA, nomeB, distances, azA, azB, northernHemisphere=False):\n",
        "    lat, lng = utmToLatLng(z, p1[0], p1[1], northernHemisphere=False)\n",
        "    g1 = (lng, lat)\n",
        "    lat, lng = utmToLatLng(z, p2[0], p2[1], northernHemisphere=False)\n",
        "    g2 = (lng, lat)\n",
        "    lat, lng = utmToLatLng(z, p3[0], p3[1], northernHemisphere=False)\n",
        "    g3 = (lng, lat)\n",
        "    kml = simplekml.Kml(open=1)\n",
        "    pntA = kml.newpoint(name=nomeA, description=\"Ponto de Observação A\",\n",
        "                   coords=[g1])  # lon, lat optional height\n",
        "    pntB = kml.newpoint(name=nomeB, description=\"Ponto de Observação B\",\n",
        "                   coords=[g2])  # lon, lat optional height\n",
        "    pntC = kml.newpoint(name=\"Cruzamento\", description=\"Incêndio\",\n",
        "                   coords=[g3])  # lon, lat optional height \n",
        "    # Create a linestring that will hover 100m above the ground that is extended to the ground\n",
        "    name1 = f'Linha de visada de A\\nDistância: {distances[1]:.2f}m\\nAzimute: {azA}{chr(176)}'\n",
        "    linestring1 = kml.newlinestring(name=name1)\n",
        "    linestring1.coords = [g1, g3]\n",
        "    linestring1.altitudemode = simplekml.AltitudeMode.clamptoground\n",
        "    linestring1.extrude = 1\n",
        "    # Create a linestring that will hover 100m above the ground that is extended to the ground\n",
        "    name2 = f'Linha de visada de B\\nDistância: {distances[2]:.2f}m\\nAzimute: {azB}{chr(176)}'\n",
        "    linestring2 = kml.newlinestring(name=name2)\n",
        "    linestring2.coords = [g2, g3]\n",
        "    linestring2.altitudemode = simplekml.AltitudeMode.clamptoground\n",
        "    linestring2.extrude = 1     \n",
        "    kml.save('linhas_pontos.kml')     \n",
        "\n",
        "def on_button_clicked(b):\n",
        "  with output:\n",
        "    distances, intersection = calcIntersection(xA,yA,azA,xB,yB,azB) \n",
        "    if not intersection.is_empty:\n",
        "      #plotResult([xA,yA],[xB,yB],[intersection.x,intersection.y],nomeA,nomeB,distances)\n",
        "      create_kml(23, [xA,yA],[xB,yB],[intersection.x,intersection.y],nomeA,nomeB,distances, azA, azB, northernHemisphere=False)\n",
        "      display(plotMapFolium(23, [xA,yA],[xB,yB],[intersection.x,intersection.y],nomeA,nomeB,distances, northernHemisphere=False))\n",
        "    else:\n",
        "      print(\"There's no crossing point!\")\n",
        "  \n",
        "button.on_click(on_button_clicked)\n",
        "display(button, output)"
      ],
      "execution_count": 5,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "Requirement already satisfied: simplekml in /usr/local/lib/python3.6/dist-packages (1.3.5)\n"
          ],
          "name": "stdout"
        },
        {
          "output_type": "display_data",
          "data": {
            "application/vnd.jupyter.widget-view+json": {
              "model_id": "a252f72c51754c19b56345fd62e8621d",
              "version_minor": 0,
              "version_major": 2
            },
            "text/plain": [
              "Button(description='Calculate!', style=ButtonStyle())"
            ]
          },
          "metadata": {
            "tags": []
          }
        },
        {
          "output_type": "display_data",
          "data": {
            "application/vnd.jupyter.widget-view+json": {
              "model_id": "59c5228d07634c0faac1818ad593b5ed",
              "version_minor": 0,
              "version_major": 2
            },
            "text/plain": [
              "Output()"
            ]
          },
          "metadata": {
            "tags": []
          }
        }
      ]
    }
  ]
}