Planning_Response_Phase_Facility_Flooding_and_Set_Cover_Optimization.ipynb

{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "name": "Planning_Response_Phase_Facility_Flooding_and_Set_Cover_Optimization",
      "version": "0.3.2",
      "provenance": [],
      "collapsed_sections": [],
      "toc_visible": true,
      "include_colab_link": true
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/bwsi-hadr/2019-student-final-exercise/blob/master/Planning_Response_Phase_Facility_Flooding_and_Set_Cover_Optimization.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "9bmz5lfibhG0",
        "colab_type": "code",
        "outputId": "f732e118-8de0-4bc9-8758-38cd81633014",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 1000
        }
      },
      "source": [
        "# need to specify location of some certificates for rasterio\n",
        "!export CURL_CA_BUNDLE=/etc/ssl/certs/ca-certificates.crt\n",
        "!sudo mkdir -p /etc/pki/tls/certs\n",
        "!sudo cp /etc/ssl/certs/ca-certificates.crt /etc/pki/tls/certs/ca-bundle.crt\n",
        "try:\n",
        "  import rasterio\n",
        "  import rasterio.plot\n",
        "  import rasterio.merge \n",
        "  import rasterio.mask\n",
        "except:\n",
        "  !pip install rasterio\n",
        "  import rasterio\n",
        "  import rasterio.plot\n",
        "  import rasterio.merge\n",
        "  import rasterio.mask\n",
        "  \n",
        "try:\n",
        "  import rasterstats as rs\n",
        "except:\n",
        "  !pip install rasterstats  \n",
        "  import rasterstats as rs\n",
        "  \n",
        "try:\n",
        "  import pyproj\n",
        "except:\n",
        "  !pip install pyproj\n",
        "  import pyproj\n",
        "  \n",
        "import networkx as nx\n",
        "try:\n",
        "  import osmnx as ox\n",
        "except:\n",
        "  # osmnx depends on the system package libspatialindex\n",
        "  !apt install libspatialindex-dev\n",
        "  !pip install osmnx\n",
        "  import osmnx as ox\n",
        "\n",
        "try: \n",
        "  import geopandas as gpd\n",
        "except: \n",
        "  !pip install geopandas \n",
        "  import geopandas as gpd\n",
        "  \n",
        "try:\n",
        "  import contextily as ctx \n",
        "except:\n",
        "  # install dependencies for contextily\n",
        "  !apt install libproj-dev proj-data proj-bin\n",
        "  !apt install libgeos-dev\n",
        "  !pip install cython\n",
        "  !pip install cartopy\n",
        "  # install contextily\n",
        "  !pip install contextily==1.0rc1 --no-use-pep517 --no-cache-dir\n",
        "  import contextily as ctx\n",
        "  \n",
        "import fiona\n",
        "from shapely.geometry import Point, LineString, Polygon\n",
        "  \n",
        "import gdal\n",
        "import pandas as pd\n",
        "import numpy as np\n",
        "import matplotlib.pyplot as plt\n",
        "import pathlib\n",
        "from google.colab import drive\n",
        "drive.mount('/content/drive')\n",
        "%cd /content/drive/'My Drive'/BWSI-Remote-Sensing/'Final_exercise'"
      ],
      "execution_count": 1,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "Collecting rasterio\n",
            "\u001b[?25l  Downloading https://files.pythonhosted.org/packages/a2/d6/f5dc5b31b5bbb8c9cf0ae2f38a31e4406994db48f9f14d9434588c935c7c/rasterio-1.0.24-cp36-cp36m-manylinux1_x86_64.whl (19.7MB)\n",
            "\u001b[K     |████████████████████████████████| 19.7MB 2.6MB/s \n",
            "\u001b[?25hRequirement already satisfied: click<8,>=4.0 in /usr/local/lib/python3.6/dist-packages (from rasterio) (7.0)\n",
            "Requirement already satisfied: attrs in /usr/local/lib/python3.6/dist-packages (from rasterio) (19.1.0)\n",
            "Collecting cligj>=0.5 (from rasterio)\n",
            "  Downloading https://files.pythonhosted.org/packages/e4/be/30a58b4b0733850280d01f8bd132591b4668ed5c7046761098d665ac2174/cligj-0.5.0-py3-none-any.whl\n",
            "Requirement already satisfied: numpy in /usr/local/lib/python3.6/dist-packages (from rasterio) (1.16.4)\n",
            "Collecting snuggs>=1.4.1 (from rasterio)\n",
            "  Downloading https://files.pythonhosted.org/packages/58/14/8e90b7586ab6929861161e73e9fd55637a060e4d14dd1be14a4b8a08751f/snuggs-1.4.6-py3-none-any.whl\n",
            "Collecting click-plugins (from rasterio)\n",
            "  Downloading https://files.pythonhosted.org/packages/e9/da/824b92d9942f4e472702488857914bdd50f73021efea15b4cad9aca8ecef/click_plugins-1.1.1-py2.py3-none-any.whl\n",
            "Collecting affine (from rasterio)\n",
            "  Downloading https://files.pythonhosted.org/packages/56/5d/6877929932d17850fa4903d0db8233ec8ed35aab7ceae96fa44ea6d479bd/affine-2.2.2-py2.py3-none-any.whl\n",
            "Requirement already satisfied: pyparsing>=2.1.6 in /usr/local/lib/python3.6/dist-packages (from snuggs>=1.4.1->rasterio) (2.4.2)\n",
            "Installing collected packages: cligj, snuggs, click-plugins, affine, rasterio\n",
            "Successfully installed affine-2.2.2 click-plugins-1.1.1 cligj-0.5.0 rasterio-1.0.24 snuggs-1.4.6\n",
            "Collecting rasterstats\n",
            "  Downloading https://files.pythonhosted.org/packages/43/30/7f0b3ed0debd10a6687110dd9d46b308d33ff4696b88c2461758d86ff428/rasterstats-0.13.1-py3-none-any.whl\n",
            "Requirement already satisfied: rasterio>=1.0 in /usr/local/lib/python3.6/dist-packages (from rasterstats) (1.0.24)\n",
            "Collecting simplejson (from rasterstats)\n",
            "\u001b[?25l  Downloading https://files.pythonhosted.org/packages/e3/24/c35fb1c1c315fc0fffe61ea00d3f88e85469004713dab488dee4f35b0aff/simplejson-3.16.0.tar.gz (81kB)\n",
            "\u001b[K     |████████████████████████████████| 81kB 5.9MB/s \n",
            "\u001b[?25hRequirement already satisfied: shapely in /usr/local/lib/python3.6/dist-packages (from rasterstats) (1.6.4.post2)\n",
            "Collecting fiona (from rasterstats)\n",
            "\u001b[?25l  Downloading https://files.pythonhosted.org/packages/89/4a/193cd6a75e51062c85f4e1cd6f312b3bbda6e26ba7510f152ef5016f0b16/Fiona-1.8.6-cp36-cp36m-manylinux1_x86_64.whl (17.9MB)\n",
            "\u001b[K     |████████████████████████████████| 17.9MB 8.2MB/s \n",
            "\u001b[?25hRequirement already satisfied: numpy>=1.9 in /usr/local/lib/python3.6/dist-packages (from rasterstats) (1.16.4)\n",
            "Requirement already satisfied: cligj>=0.4 in /usr/local/lib/python3.6/dist-packages (from rasterstats) (0.5.0)\n",
            "Requirement already satisfied: affine in /usr/local/lib/python3.6/dist-packages (from rasterio>=1.0->rasterstats) (2.2.2)\n",
            "Requirement already satisfied: attrs in /usr/local/lib/python3.6/dist-packages (from rasterio>=1.0->rasterstats) (19.1.0)\n",
            "Requirement already satisfied: snuggs>=1.4.1 in /usr/local/lib/python3.6/dist-packages (from rasterio>=1.0->rasterstats) (1.4.6)\n",
            "Requirement already satisfied: click-plugins in /usr/local/lib/python3.6/dist-packages (from rasterio>=1.0->rasterstats) (1.1.1)\n",
            "Requirement already satisfied: click<8,>=4.0 in /usr/local/lib/python3.6/dist-packages (from rasterio>=1.0->rasterstats) (7.0)\n",
            "Collecting munch (from fiona->rasterstats)\n",
            "  Downloading https://files.pythonhosted.org/packages/68/f4/260ec98ea840757a0da09e0ed8135333d59b8dfebe9752a365b04857660a/munch-2.3.2.tar.gz\n",
            "Requirement already satisfied: six>=1.7 in /usr/local/lib/python3.6/dist-packages (from fiona->rasterstats) (1.12.0)\n",
            "Requirement already satisfied: pyparsing>=2.1.6 in /usr/local/lib/python3.6/dist-packages (from snuggs>=1.4.1->rasterio>=1.0->rasterstats) (2.4.2)\n",
            "Building wheels for collected packages: simplejson, munch\n",
            "  Building wheel for simplejson (setup.py) ... \u001b[?25l\u001b[?25hdone\n",
            "  Created wheel for simplejson: filename=simplejson-3.16.0-cp36-cp36m-linux_x86_64.whl size=114021 sha256=6b754396a8df6d18ba573033e02b02156c001059af0842b706c4d61e5d8d3ab6\n",
            "  Stored in directory: /root/.cache/pip/wheels/5d/1a/1e/0350bb3df3e74215cd91325344cc86c2c691f5306eb4d22c77\n",
            "  Building wheel for munch (setup.py) ... \u001b[?25l\u001b[?25hdone\n",
            "  Created wheel for munch: filename=munch-2.3.2-py2.py3-none-any.whl size=6613 sha256=8e6db14330c81901ee82c1ee8c82b98b8b4d40f1e23f375611656bebf4ca5004\n",
            "  Stored in directory: /root/.cache/pip/wheels/db/bf/bc/06a3e1bfe0ab27d2e720ceb3cff3159398d92644c0cec2c125\n",
            "Successfully built simplejson munch\n",
            "Installing collected packages: simplejson, munch, fiona, rasterstats\n",
            "Successfully installed fiona-1.8.6 munch-2.3.2 rasterstats-0.13.1 simplejson-3.16.0\n",
            "Collecting pyproj\n",
            "\u001b[?25l  Downloading https://files.pythonhosted.org/packages/16/59/43869adef45ce4f1cf7d5c3aef1ea5d65d449050abdda5de7a2465c5729d/pyproj-2.2.1-cp36-cp36m-manylinux1_x86_64.whl (11.2MB)\n",
            "\u001b[K     |████████████████████████████████| 11.2MB 2.7MB/s \n",
            "\u001b[?25hInstalling collected packages: pyproj\n",
            "Successfully installed pyproj-2.2.1\n",
            "Reading package lists... Done\n",
            "Building dependency tree       \n",
            "Reading state information... Done\n",
            "The following package was automatically installed and is no longer required:\n",
            "  libnvidia-common-410\n",
            "Use 'apt autoremove' to remove it.\n",
            "The following additional packages will be installed:\n",
            "  libspatialindex-c4v5 libspatialindex4v5\n",
            "The following NEW packages will be installed:\n",
            "  libspatialindex-c4v5 libspatialindex-dev libspatialindex4v5\n",
            "0 upgraded, 3 newly installed, 0 to remove and 4 not upgraded.\n",
            "Need to get 555 kB of archives.\n",
            "After this operation, 3,308 kB of additional disk space will be used.\n",
            "Get:1 http://archive.ubuntu.com/ubuntu bionic/universe amd64 libspatialindex4v5 amd64 1.8.5-5 [219 kB]\n",
            "Get:2 http://archive.ubuntu.com/ubuntu bionic/universe amd64 libspatialindex-c4v5 amd64 1.8.5-5 [51.7 kB]\n",
            "Get:3 http://archive.ubuntu.com/ubuntu bionic/universe amd64 libspatialindex-dev amd64 1.8.5-5 [285 kB]\n",
            "Fetched 555 kB in 1s (660 kB/s)\n",
            "Selecting previously unselected package libspatialindex4v5:amd64.\n",
            "(Reading database ... 131289 files and directories currently installed.)\n",
            "Preparing to unpack .../libspatialindex4v5_1.8.5-5_amd64.deb ...\n",
            "Unpacking libspatialindex4v5:amd64 (1.8.5-5) ...\n",
            "Selecting previously unselected package libspatialindex-c4v5:amd64.\n",
            "Preparing to unpack .../libspatialindex-c4v5_1.8.5-5_amd64.deb ...\n",
            "Unpacking libspatialindex-c4v5:amd64 (1.8.5-5) ...\n",
            "Selecting previously unselected package libspatialindex-dev:amd64.\n",
            "Preparing to unpack .../libspatialindex-dev_1.8.5-5_amd64.deb ...\n",
            "Unpacking libspatialindex-dev:amd64 (1.8.5-5) ...\n",
            "Setting up libspatialindex4v5:amd64 (1.8.5-5) ...\n",
            "Setting up libspatialindex-c4v5:amd64 (1.8.5-5) ...\n",
            "Processing triggers for libc-bin (2.27-3ubuntu1) ...\n",
            "Setting up libspatialindex-dev:amd64 (1.8.5-5) ...\n",
            "Collecting osmnx\n",
            "\u001b[?25l  Downloading https://files.pythonhosted.org/packages/69/77/98cbee7f94abe2c60013fd8751f2b2d6054d48078c61f9774cff0b341f43/osmnx-0.10-py2.py3-none-any.whl (70kB)\n",
            "\u001b[K     |████████████████████████████████| 71kB 2.7MB/s \n",
            "\u001b[?25hRequirement already satisfied: networkx>=2.3 in /usr/local/lib/python3.6/dist-packages (from osmnx) (2.3)\n",
            "Requirement already satisfied: descartes>=1.1 in /usr/local/lib/python3.6/dist-packages (from osmnx) (1.1.0)\n",
            "Collecting geopandas>=0.5 (from osmnx)\n",
            "\u001b[?25l  Downloading https://files.pythonhosted.org/packages/21/80/da2a33c9201cd4ce693f4aa6189efc9ef1a48bec1c3b02c3ce9908b07fec/geopandas-0.5.1-py2.py3-none-any.whl (893kB)\n",
            "\u001b[K     |████████████████████████████████| 901kB 8.0MB/s \n",
            "\u001b[?25hRequirement already satisfied: numpy>=1.16 in /usr/local/lib/python3.6/dist-packages (from osmnx) (1.16.4)\n",
            "Requirement already satisfied: matplotlib>=2.2 in /usr/local/lib/python3.6/dist-packages (from osmnx) (3.0.3)\n",
            "Requirement already satisfied: Shapely>=1.6 in /usr/local/lib/python3.6/dist-packages (from osmnx) (1.6.4.post2)\n",
            "Collecting Rtree>=0.8.3 (from osmnx)\n",
            "  Downloading https://files.pythonhosted.org/packages/33/0d/0bcb0fdd759570e85697189f675fa462131e5636482f2832288bbe4bb6f7/Rtree-0.8.3-py3-none-any.whl\n",
            "Requirement already satisfied: pandas>=0.24 in /usr/local/lib/python3.6/dist-packages (from osmnx) (0.24.2)\n",
            "Requirement already satisfied: requests>=2.21 in /usr/local/lib/python3.6/dist-packages (from osmnx) (2.21.0)\n",
            "Requirement already satisfied: decorator>=4.3.0 in /usr/local/lib/python3.6/dist-packages (from networkx>=2.3->osmnx) (4.4.0)\n",
            "Requirement already satisfied: fiona in /usr/local/lib/python3.6/dist-packages (from geopandas>=0.5->osmnx) (1.8.6)\n",
            "Requirement already satisfied: pyproj in /usr/local/lib/python3.6/dist-packages (from geopandas>=0.5->osmnx) (2.2.1)\n",
            "Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.6/dist-packages (from matplotlib>=2.2->osmnx) (0.10.0)\n",
            "Requirement already satisfied: kiwisolver>=1.0.1 in /usr/local/lib/python3.6/dist-packages (from matplotlib>=2.2->osmnx) (1.1.0)\n",
            "Requirement already satisfied: python-dateutil>=2.1 in /usr/local/lib/python3.6/dist-packages (from matplotlib>=2.2->osmnx) (2.5.3)\n",
            "Requirement already satisfied: pyparsing!=2.0.4,!=2.1.2,!=2.1.6,>=2.0.1 in /usr/local/lib/python3.6/dist-packages (from matplotlib>=2.2->osmnx) (2.4.2)\n",
            "Requirement already satisfied: setuptools in /usr/local/lib/python3.6/dist-packages (from Rtree>=0.8.3->osmnx) (41.0.1)\n",
            "Requirement already satisfied: pytz>=2011k in /usr/local/lib/python3.6/dist-packages (from pandas>=0.24->osmnx) (2018.9)\n",
            "Requirement already satisfied: idna<2.9,>=2.5 in /usr/local/lib/python3.6/dist-packages (from requests>=2.21->osmnx) (2.8)\n",
            "Requirement already satisfied: certifi>=2017.4.17 in /usr/local/lib/python3.6/dist-packages (from requests>=2.21->osmnx) (2019.6.16)\n",
            "Requirement already satisfied: urllib3<1.25,>=1.21.1 in /usr/local/lib/python3.6/dist-packages (from requests>=2.21->osmnx) (1.24.3)\n",
            "Requirement already satisfied: chardet<3.1.0,>=3.0.2 in /usr/local/lib/python3.6/dist-packages (from requests>=2.21->osmnx) (3.0.4)\n",
            "Requirement already satisfied: attrs>=17 in /usr/local/lib/python3.6/dist-packages (from fiona->geopandas>=0.5->osmnx) (19.1.0)\n",
            "Requirement already satisfied: click<8,>=4.0 in /usr/local/lib/python3.6/dist-packages (from fiona->geopandas>=0.5->osmnx) (7.0)\n",
            "Requirement already satisfied: six>=1.7 in /usr/local/lib/python3.6/dist-packages (from fiona->geopandas>=0.5->osmnx) (1.12.0)\n",
            "Requirement already satisfied: cligj>=0.5 in /usr/local/lib/python3.6/dist-packages (from fiona->geopandas>=0.5->osmnx) (0.5.0)\n",
            "Requirement already satisfied: click-plugins>=1.0 in /usr/local/lib/python3.6/dist-packages (from fiona->geopandas>=0.5->osmnx) (1.1.1)\n",
            "Requirement already satisfied: munch in /usr/local/lib/python3.6/dist-packages (from fiona->geopandas>=0.5->osmnx) (2.3.2)\n",
            "Installing collected packages: geopandas, Rtree, osmnx\n",
            "Successfully installed Rtree-0.8.3 geopandas-0.5.1 osmnx-0.10\n",
            "Reading package lists... Done\n",
            "Building dependency tree       \n",
            "Reading state information... Done\n",
            "proj-data is already the newest version (4.9.3-2).\n",
            "proj-data set to manually installed.\n",
            "The following package was automatically installed and is no longer required:\n",
            "  libnvidia-common-410\n",
            "Use 'apt autoremove' to remove it.\n",
            "The following NEW packages will be installed:\n",
            "  libproj-dev proj-bin\n",
            "0 upgraded, 2 newly installed, 0 to remove and 4 not upgraded.\n",
            "Need to get 232 kB of archives.\n",
            "After this operation, 1,220 kB of additional disk space will be used.\n",
            "Get:1 http://archive.ubuntu.com/ubuntu bionic/universe amd64 libproj-dev amd64 4.9.3-2 [199 kB]\n",
            "Get:2 http://archive.ubuntu.com/ubuntu bionic/universe amd64 proj-bin amd64 4.9.3-2 [32.3 kB]\n",
            "Fetched 232 kB in 1s (370 kB/s)\n",
            "Selecting previously unselected package libproj-dev:amd64.\n",
            "(Reading database ... 131342 files and directories currently installed.)\n",
            "Preparing to unpack .../libproj-dev_4.9.3-2_amd64.deb ...\n",
            "Unpacking libproj-dev:amd64 (4.9.3-2) ...\n",
            "Selecting previously unselected package proj-bin.\n",
            "Preparing to unpack .../proj-bin_4.9.3-2_amd64.deb ...\n",
            "Unpacking proj-bin (4.9.3-2) ...\n",
            "Setting up libproj-dev:amd64 (4.9.3-2) ...\n",
            "Setting up proj-bin (4.9.3-2) ...\n",
            "Processing triggers for man-db (2.8.3-2ubuntu0.1) ...\n",
            "Reading package lists... Done\n",
            "Building dependency tree       \n",
            "Reading state information... Done\n",
            "The following package was automatically installed and is no longer required:\n",
            "  libnvidia-common-410\n",
            "Use 'apt autoremove' to remove it.\n",
            "Suggested packages:\n",
            "  libgdal-doc\n",
            "The following NEW packages will be installed:\n",
            "  libgeos-dev\n",
            "0 upgraded, 1 newly installed, 0 to remove and 4 not upgraded.\n",
            "Need to get 73.1 kB of archives.\n",
            "After this operation, 486 kB of additional disk space will be used.\n",
            "Get:1 http://archive.ubuntu.com/ubuntu bionic/universe amd64 libgeos-dev amd64 3.6.2-1build2 [73.1 kB]\n",
            "Fetched 73.1 kB in 0s (178 kB/s)\n",
            "Selecting previously unselected package libgeos-dev.\n",
            "(Reading database ... 131375 files and directories currently installed.)\n",
            "Preparing to unpack .../libgeos-dev_3.6.2-1build2_amd64.deb ...\n",
            "Unpacking libgeos-dev (3.6.2-1build2) ...\n",
            "Setting up libgeos-dev (3.6.2-1build2) ...\n",
            "Processing triggers for man-db (2.8.3-2ubuntu0.1) ...\n",
            "Requirement already satisfied: cython in /usr/local/lib/python3.6/dist-packages (0.29.13)\n",
            "Collecting cartopy\n",
            "\u001b[?25l  Downloading https://files.pythonhosted.org/packages/e5/92/fe8838fa8158931906dfc4f16c5c1436b3dd2daf83592645b179581403ad/Cartopy-0.17.0.tar.gz (8.9MB)\n",
            "\u001b[K     |████████████████████████████████| 8.9MB 2.8MB/s \n",
            "\u001b[?25h  Installing build dependencies ... \u001b[?25l\u001b[?25hdone\n",
            "  Getting requirements to build wheel ... \u001b[?25l\u001b[?25hdone\n",
            "    Preparing wheel metadata ... \u001b[?25l\u001b[?25hdone\n",
            "Requirement already satisfied: numpy>=1.10 in /usr/local/lib/python3.6/dist-packages (from cartopy) (1.16.4)\n",
            "Requirement already satisfied: setuptools>=0.7.2 in /usr/local/lib/python3.6/dist-packages (from cartopy) (41.0.1)\n",
            "Requirement already satisfied: six>=1.3.0 in /usr/local/lib/python3.6/dist-packages (from cartopy) (1.12.0)\n",
            "Requirement already satisfied: shapely>=1.5.6 in /usr/local/lib/python3.6/dist-packages (from cartopy) (1.6.4.post2)\n",
            "Collecting pyshp>=1.1.4 (from cartopy)\n",
            "\u001b[?25l  Downloading https://files.pythonhosted.org/packages/27/16/3bf15aa864fb77845fab8007eda22c2bd67bd6c1fd13496df452c8c43621/pyshp-2.1.0.tar.gz (215kB)\n",
            "\u001b[K     |████████████████████████████████| 225kB 40.1MB/s \n",
            "\u001b[?25hBuilding wheels for collected packages: cartopy\n",
            "  Building wheel for cartopy (PEP 517) ... \u001b[?25l\u001b[?25hdone\n",
            "  Created wheel for cartopy: filename=Cartopy-0.17.0-cp36-cp36m-linux_x86_64.whl size=9714102 sha256=ed285d3bc65815f651683e3e693e269f249dac14dee6a9af753c4740b18edd27\n",
            "  Stored in directory: /root/.cache/pip/wheels/cd/cf/40/539f798f94e921e94fd376a5f9d213a6febe77754c0b187c73\n",
            "Successfully built cartopy\n",
            "Building wheels for collected packages: pyshp\n",
            "  Building wheel for pyshp (setup.py) ... \u001b[?25l\u001b[?25hdone\n",
            "  Created wheel for pyshp: filename=pyshp-2.1.0-cp36-none-any.whl size=32608 sha256=30ee11b98d2641938baf116e0f5f8f418a363e73b5e042f830cc0132bebf0ee6\n",
            "  Stored in directory: /root/.cache/pip/wheels/a6/0c/de/321b5192ad416b328975a2f0385f72c64db4656501eba7cc1a\n",
            "Successfully built pyshp\n",
            "Installing collected packages: pyshp, cartopy\n",
            "Successfully installed cartopy-0.17.0 pyshp-2.1.0\n",
            "Collecting contextily==1.0rc1\n",
            "  Downloading https://files.pythonhosted.org/packages/1c/1c/534e6000b49003d0dc5261298e46853e405d5e27bae8232d1f3d4745b3de/contextily-1.0rc1-py3-none-any.whl\n",
            "Requirement already satisfied: cartopy in /usr/local/lib/python3.6/dist-packages (from contextily==1.0rc1) (0.17.0)\n",
            "Requirement already satisfied: rasterio in /usr/local/lib/python3.6/dist-packages (from contextily==1.0rc1) (1.0.24)\n",
            "Requirement already satisfied: pillow in /usr/local/lib/python3.6/dist-packages (from contextily==1.0rc1) (4.3.0)\n",
            "Requirement already satisfied: requests in /usr/local/lib/python3.6/dist-packages (from contextily==1.0rc1) (2.21.0)\n",
            "Collecting mercantile (from contextily==1.0rc1)\n",
            "  Downloading https://files.pythonhosted.org/packages/f8/d4/f79431df2ea725206e845d06980a2eb7f4264aa6a7cc749647dd6233d6fa/mercantile-1.1.1-py3-none-any.whl\n",
            "Requirement already satisfied: geopy in /usr/local/lib/python3.6/dist-packages (from contextily==1.0rc1) (1.17.0)\n",
            "Requirement already satisfied: pytest in /usr/local/lib/python3.6/dist-packages (from contextily==1.0rc1) (3.6.4)\n",
            "Requirement already satisfied: pandas in /usr/local/lib/python3.6/dist-packages (from contextily==1.0rc1) (0.24.2)\n",
            "Requirement already satisfied: pyshp>=1.1.4 in /usr/local/lib/python3.6/dist-packages (from cartopy->contextily==1.0rc1) (2.1.0)\n",
            "Requirement already satisfied: six>=1.3.0 in /usr/local/lib/python3.6/dist-packages (from cartopy->contextily==1.0rc1) (1.12.0)\n",
            "Requirement already satisfied: shapely>=1.5.6 in /usr/local/lib/python3.6/dist-packages (from cartopy->contextily==1.0rc1) (1.6.4.post2)\n",
            "Requirement already satisfied: numpy>=1.10 in /usr/local/lib/python3.6/dist-packages (from cartopy->contextily==1.0rc1) (1.16.4)\n",
            "Requirement already satisfied: setuptools>=0.7.2 in /usr/local/lib/python3.6/dist-packages (from cartopy->contextily==1.0rc1) (41.0.1)\n",
            "Requirement already satisfied: affine in /usr/local/lib/python3.6/dist-packages (from rasterio->contextily==1.0rc1) (2.2.2)\n",
            "Requirement already satisfied: click<8,>=4.0 in /usr/local/lib/python3.6/dist-packages (from rasterio->contextily==1.0rc1) (7.0)\n",
            "Requirement already satisfied: cligj>=0.5 in /usr/local/lib/python3.6/dist-packages (from rasterio->contextily==1.0rc1) (0.5.0)\n",
            "Requirement already satisfied: attrs in /usr/local/lib/python3.6/dist-packages (from rasterio->contextily==1.0rc1) (19.1.0)\n",
            "Requirement already satisfied: click-plugins in /usr/local/lib/python3.6/dist-packages (from rasterio->contextily==1.0rc1) (1.1.1)\n",
            "Requirement already satisfied: snuggs>=1.4.1 in /usr/local/lib/python3.6/dist-packages (from rasterio->contextily==1.0rc1) (1.4.6)\n",
            "Requirement already satisfied: olefile in /usr/local/lib/python3.6/dist-packages (from pillow->contextily==1.0rc1) (0.46)\n",
            "Requirement already satisfied: idna<2.9,>=2.5 in /usr/local/lib/python3.6/dist-packages (from requests->contextily==1.0rc1) (2.8)\n",
            "Requirement already satisfied: chardet<3.1.0,>=3.0.2 in /usr/local/lib/python3.6/dist-packages (from requests->contextily==1.0rc1) (3.0.4)\n",
            "Requirement already satisfied: urllib3<1.25,>=1.21.1 in /usr/local/lib/python3.6/dist-packages (from requests->contextily==1.0rc1) (1.24.3)\n",
            "Requirement already satisfied: certifi>=2017.4.17 in /usr/local/lib/python3.6/dist-packages (from requests->contextily==1.0rc1) (2019.6.16)\n",
            "Requirement already satisfied: geographiclib<2,>=1.49 in /usr/local/lib/python3.6/dist-packages (from geopy->contextily==1.0rc1) (1.49)\n",
            "Requirement already satisfied: more-itertools>=4.0.0 in /usr/local/lib/python3.6/dist-packages (from pytest->contextily==1.0rc1) (7.2.0)\n",
            "Requirement already satisfied: atomicwrites>=1.0 in /usr/local/lib/python3.6/dist-packages (from pytest->contextily==1.0rc1) (1.3.0)\n",
            "Requirement already satisfied: py>=1.5.0 in /usr/local/lib/python3.6/dist-packages (from pytest->contextily==1.0rc1) (1.8.0)\n",
            "Requirement already satisfied: pluggy<0.8,>=0.5 in /usr/local/lib/python3.6/dist-packages (from pytest->contextily==1.0rc1) (0.7.1)\n",
            "Requirement already satisfied: pytz>=2011k in /usr/local/lib/python3.6/dist-packages (from pandas->contextily==1.0rc1) (2018.9)\n",
            "Requirement already satisfied: python-dateutil>=2.5.0 in /usr/local/lib/python3.6/dist-packages (from pandas->contextily==1.0rc1) (2.5.3)\n",
            "Requirement already satisfied: pyparsing>=2.1.6 in /usr/local/lib/python3.6/dist-packages (from snuggs>=1.4.1->rasterio->contextily==1.0rc1) (2.4.2)\n",
            "Installing collected packages: mercantile, contextily\n",
            "Successfully installed contextily-1.0rc1 mercantile-1.1.1\n",
            "Go to this URL in a browser: https://accounts.google.com/o/oauth2/auth?client_id=947318989803-6bn6qk8qdgf4n4g3pfee6491hc0brc4i.apps.googleusercontent.com&redirect_uri=urn%3Aietf%3Awg%3Aoauth%3A2.0%3Aoob&scope=email%20https%3A%2F%2Fwww.googleapis.com%2Fauth%2Fdocs.test%20https%3A%2F%2Fwww.googleapis.com%2Fauth%2Fdrive%20https%3A%2F%2Fwww.googleapis.com%2Fauth%2Fdrive.photos.readonly%20https%3A%2F%2Fwww.googleapis.com%2Fauth%2Fpeopleapi.readonly&response_type=code\n",
            "\n",
            "Enter your authorization code:\n",
            "··········\n",
            "Mounted at /content/drive\n",
            "/content/drive/My Drive/BWSI-Remote-Sensing/Final_exercise\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "tcoz8ygszAnC",
        "colab_type": "text"
      },
      "source": [
        ""
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "G5APaCHIbojw",
        "colab_type": "code",
        "outputId": "25d24384-05de-4f2e-f9a1-d323bba73f73",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 118
        }
      },
      "source": [
        "!ls game_grid_export"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "facilities\t      flooding.shx\t       track.cpg\n",
            "facilities_stats.csv  game_grid_all_stats.cpg  track.dbf\n",
            "flooding.cpg\t      game_grid_all_stats.dbf  track.prj\n",
            "flooding.dbf\t      game_grid_all_stats.prj  track.qpj\n",
            "flooding.prj\t      game_grid_all_stats.shp  track.shp\n",
            "flooding.shp\t      game_grid_all_stats.shx  track.shx\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "Nv4vusci5bV4",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "flooding = gpd.read_file('game_grid_export/flooding.shp')"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "g1V4pacmb6op",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "grid = gpd.read_file('game_grid_export/game_grid_all_stats.shp').to_crs(epsg = 3857)\n",
        "grid['centroid'] = grid['geometry'].centroid"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "yeGgnGXju24h",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "to_m = lambda km: km * 1000"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "Rb7An9RmIVmE",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "from shapely.geometry import MultiPolygon\n",
        "import shapely.ops as shap\n",
        "#Creates a MultiPolygon using the flood shapely file and the buffer size in meters\n",
        "def flcreate(fl, b, i): \n",
        "  flooding = fl.copy().to_crs(grid.crs)\n",
        "  flooding['geometry'] = flooding['geometry'].buffer(b)\n",
        "  if i == 0:\n",
        "    flood = flooding['geometry'].unary_union\n",
        "  else:\n",
        "    l = [MultiPolygon([flooding.loc[i, 'geometry'], \\\n",
        "                       flooding.loc[i + 1, 'geometry']]).convex_hull \\\n",
        "         for i in range(len(flooding) - 1)]\n",
        "    flood = shap.cascaded_union(MultiPolygon(l))\n",
        "    \n",
        "  return flood\n",
        "    "
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "t9-uSpIMiEAq",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "#Outputs a list of facilities in the specified distance from the flood\n",
        "def uprox(mindist, maxdist, fdf_idx, flood):\n",
        "  l = []\n",
        "  for i in range(len(fdf_idx)):\n",
        "    dist = fdf_idx.loc[i, 'geometry'].distance(flood)\n",
        "    if dist <= maxdist and dist >= mindist: #TODO: mindist <= dist <= maxdist check\n",
        "      if fdf_idx.loc[i, 'NAME'] not in l:\n",
        "        l.append(fdf_idx.loc[i, 'NAME'])\n",
        "  return l"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "dKYYPli621ve",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "def wcalc_weight(subset, t, cost_ref, threshold):\n",
        "  s = 0\n",
        "  for cell in subset:\n",
        "    if grid.loc[cell, 'Health'] <= threshold:\n",
        "      s += grid.loc[cell, 'population']\n",
        "  w = s / cost_ref[t]\n",
        "  return w"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "1nLqj-Vx0rtQ",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "#freq_chart = [0 for _ in range(game_grid.shape[0])]\n",
        "covered = set()\n",
        "\n",
        "def wset_cover(budget, facdict, cost_ref, threshold): #facdict = hosp_to_cells, EMS_to_cells, shelt_to_cells, fire_to_cells\n",
        "  to_ret = set()\n",
        "  while budget > 0 and len(facdict.items()) != 0:\n",
        "    max_fac, max_subset = max(facdict.items(), key = lambda tup: \\\n",
        "                               wcalc_weight(tup[1], tup[0][1], cost_ref, threshold))\n",
        "#     if wcalc_weight(max_subset, max_fac[1], cost_ref, threshold) == 0:\n",
        "#       max_fac, max_subset = max(facdict.items(), key = lambda tup: \\\n",
        "#                                calc_weight(tup[1], tup[0][1], cost_ref))\n",
        "    \n",
        "#Set cover could be wrong because of overlapping areas\n",
        "#    for cell in max_subset:\n",
        "#      freq_chart[cell] += 1\n",
        "#       if freq_chart[cell] == 2:\n",
        "#         done.add(cell)\n",
        "    \n",
        "    budget -= cost_ref[max_fac[1]]\n",
        "    to_ret.add(max_fac)\n",
        "    \n",
        "    for key in facdict:\n",
        "      facdict[key] -= max_subset\n",
        "    del facdict[max_fac]\n",
        "    \n",
        "  return to_ret"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "k0Mm1sux1Y4c",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "#costs is a list of costs in the order of large hospitals, small hospitals\n",
        "#large EMS, small EMS, shelters, and fire stations\n",
        "def wuoptfac (costs, budget, fac_info, threshold = 6):\n",
        "  large_hos_cost = costs[0]\n",
        "  small_hos_cost = costs[1]\n",
        "  large_ems_cost = costs[2]\n",
        "  small_ems_cost = costs[3]\n",
        "  shelter_cost = costs[4]\n",
        "  fstation_cost = costs[5]\n",
        "\n",
        "  cost_ref = {'smallhosp': small_hos_cost, 'largehosp': large_hos_cost, \\\n",
        "              'smallEMS': small_ems_cost, 'largeEMS': large_ems_cost, \\\n",
        "              'shelter': shelter_cost, 'fstation': fstation_cost}\n",
        "  fac_to_cells = {}\n",
        "  for i, row in fac_info.iterrows():\n",
        "    fac_name = row[\"NAME\"]\n",
        "    f = ''\n",
        "    if row['type'] == \"small hospital\":\n",
        "      f = 'smallhosp'\n",
        "    elif row['type'] == \"large hospital\":\n",
        "      f = 'largehosp'\n",
        "    elif row['type'] == \"small EMS\":\n",
        "      f = 'smallEMS'\n",
        "    elif row['type'] == \"large EMS\":\n",
        "      f = 'largeEMS'\n",
        "    elif row['type'] == \"shelter\":\n",
        "      f = 'shelter'\n",
        "    elif row['type'] == \"fire station\":\n",
        "      f = 'fstation'\n",
        "    else:\n",
        "      return 'no'\n",
        "    if (fac_name, f) not in fac_to_cells:\n",
        "      fac_to_cells[(fac_name, f)] = {i}\n",
        "    else:\n",
        "      fac_to_cells[(fac_name, f)].add(i)\n",
        "  \n",
        "  return wset_cover(budget, fac_to_cells, cost_ref, threshold)\n",
        "    "
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "uR_pqXL47kY_",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "def wuabracadabra(costs, budget, mindist, maxdist, buffer, fdf_idx, fdf_info, alpha_fdf, fl = flooding, graph_type = 0, threshold = 6):\n",
        "  l = uprox(mindist, maxdist, fdf_idx, flcreate(fl, buffer, graph_type))\n",
        "\n",
        "  x = wuoptfac(costs, budget, fdf_info.loc[fdf_info['NAME'].isin(l)], threshold)\n",
        "  l = [i[0] for i in x]\n",
        "  \n",
        "  return alpha_fdf.loc[l]"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "8VWWvL3fE3oY",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "#Input the current grid, the flooding shapely, the costs, budget, minimum distance, maximum distance\n",
        "#buffer, graph_type, and threshold for health\n",
        "def main(grid, flooding, costs, budget, mindist, maxdist, buffer, graph_type, threshold):\n",
        "  #Hospitals\n",
        "  hospital_grid = gpd.read_file('game_grid_export/facilities/hospitals.shp').to_crs(grid.crs)\n",
        "  hospital_grid.loc[hospital_grid.loc[:, 'BEDS'] >= 200, 'radius'] = 25\n",
        "  hospital_grid.loc[hospital_grid.loc[:, 'BEDS'] < 200, 'radius'] = 15\n",
        "  hospital_grid.loc[hospital_grid.loc[:, 'BEDS'] >= 200, 'type'] = \"large hospital\"\n",
        "  hospital_grid.loc[hospital_grid.loc[:, 'BEDS'] < 200, 'type'] = \"small hospital\"\n",
        "  hospital_grid['health'] = 2\n",
        "\n",
        "  big_hospitals=hospital_grid[hospital_grid['BEDS']>=200].to_crs(grid.crs)\n",
        "  small_hospitals=hospital_grid[hospital_grid['BEDS']<200].to_crs(grid.crs)\n",
        "  hospitals_buffered = hospital_grid.copy().to_crs(grid.crs)\n",
        "  hospitals_buffered.loc[big_hospitals.index, 'geometry'] = hospitals_buffered.loc[big_hospitals.index,'geometry'].buffer(to_m(25))\n",
        "  hospitals_buffered.loc[small_hospitals.index, 'geometry'] = hospitals_buffered.loc[small_hospitals.index,'geometry'].buffer(to_m(15))\n",
        "  hospitals_buffered['centroid'] = hospitals_buffered['geometry'].centroid\n",
        "\n",
        "  hosp_info = gpd.sjoin(grid, hospitals_buffered, how = 'inner', op = 'intersects')\n",
        "  hosp_infopop = hosp_info.groupby('NAME').sum()\n",
        "\n",
        "  alpha_hosp = hospitals_buffered.sort_values('NAME').set_index('NAME')\n",
        "  alpha_hosp['population'] = hosp_infopop['population'].values\n",
        "  hosp_idx = alpha_hosp.reset_index()\n",
        "\n",
        "  #EMS\n",
        "  EMS_grid = gpd.read_file('game_grid_export/facilities/EMS.shp').to_crs(grid.crs)\n",
        "  EMS_grid.loc[EMS_grid.loc[:, 'TOTAL_VEHI'] >= 5, 'radius'] = 15\n",
        "  EMS_grid.loc[EMS_grid.loc[:, 'TOTAL_VEHI'] < 5, 'radius'] = 10\n",
        "  EMS_grid.loc[EMS_grid.loc[:, 'TOTAL_VEHI'] >= 5, 'type'] = \"large EMS\"\n",
        "  EMS_grid.loc[EMS_grid.loc[:, 'TOTAL_VEHI'] < 5, 'type'] = \"small EMS\"\n",
        "  EMS_grid['health'] = 1\n",
        "\n",
        "  big_EMS=EMS_grid[EMS_grid['TOTAL_VEHI']>=5].to_crs(grid.crs)\n",
        "  small_EMS=EMS_grid[EMS_grid['TOTAL_VEHI']<5].to_crs(grid.crs)\n",
        "  EMS_buffered = EMS_grid.copy().to_crs(grid.crs)\n",
        "  EMS_buffered.loc[big_EMS.index, 'geometry'] = EMS_buffered.loc[big_EMS.index,'geometry'].buffer(to_m(15))\n",
        "  EMS_buffered.loc[small_EMS.index, 'geometry'] = EMS_buffered.loc[small_EMS.index,'geometry'].buffer(to_m(10))\n",
        "  EMS_buffered['centroid'] = EMS_buffered['geometry'].centroid\n",
        "\n",
        "  EMS_info = gpd.sjoin(grid, EMS_buffered, how = 'inner', op = 'intersects')\n",
        "  EMS_infopop = EMS_info.groupby('NAME').sum()\n",
        "\n",
        "  alpha_EMS = EMS_buffered.sort_values('NAME').set_index('NAME')\n",
        "  alpha_EMS = alpha_EMS.loc[~alpha_EMS.index.duplicated(keep='first')]\n",
        "  alpha_EMS.loc[EMS_infopop.index,'population'] = EMS_infopop['population'].values\n",
        "  alpha_EMS = alpha_EMS.dropna(subset=['population'])\n",
        "  EMS_idx = alpha_EMS.reset_index()\n",
        "\n",
        "  #Shelters\n",
        "  shelt_grid = gpd.read_file('game_grid_export/facilities/shelters.shp').to_crs(grid.crs)\n",
        "  shelt_grid['radius'] = 3\n",
        "  shelt_grid['type'] = \"shelter\"\n",
        "  shelt_grid['health'] = 1\n",
        "  shelt_buffered = shelt_grid.copy().to_crs(grid.crs)\n",
        "  shelt_buffered['geometry'] = shelt_buffered['geometry'].buffer(to_m(3))\n",
        "  shelt_buffered['centroid'] = shelt_buffered['geometry'].centroid\n",
        "\n",
        "  shelt_info = gpd.sjoin(grid, shelt_buffered, how = 'inner', op = 'intersects')\n",
        "  shelt_infopop = shelt_info.groupby('NAME').sum()\n",
        "\n",
        "  alpha_shelt = shelt_buffered.sort_values('NAME').set_index('NAME')\n",
        "  alpha_shelt = alpha_shelt.loc[~alpha_shelt.index.duplicated(keep='first')]\n",
        "  alpha_shelt.loc[shelt_infopop.index,'population'] = shelt_infopop['population'].values\n",
        "  alpha_shelt = alpha_shelt.dropna(subset=['population'])\n",
        "  shelt_idx = alpha_shelt.reset_index()\n",
        "\n",
        "  #Fire Stations\n",
        "  fire_grid = gpd.read_file('game_grid_export/facilities/fire_stations.shp').to_crs(grid.crs)\n",
        "  fire_grid['radius'] = 15\n",
        "  fire_grid['type'] = \"fire station\"\n",
        "  fire_grid['health'] = 1\n",
        "  fire_buffered = fire_grid.copy().to_crs(grid.crs)\n",
        "  fire_buffered['geometry'] = fire_buffered['geometry'].buffer(to_m(15))\n",
        "  fire_buffered['centroid'] = fire_buffered['geometry'].centroid\n",
        "\n",
        "  fire_info = gpd.sjoin(grid, fire_buffered, how = 'inner', op = 'intersects')\n",
        "  fire_infopop = fire_info.groupby('NAME').sum()\n",
        "\n",
        "  alpha_fire = fire_buffered.sort_values('NAME').set_index('NAME')\n",
        "  alpha_fire = alpha_fire.loc[~alpha_fire.index.duplicated(keep='first')]\n",
        "  alpha_fire.loc[fire_infopop.index,'population'] = fire_infopop['population'].values\n",
        "  alpha_fire = alpha_fire.dropna(subset=['population'])\n",
        "  fire_idx = alpha_fire.reset_index()\n",
        "\n",
        "  #Congregation\n",
        "  col = ['NAME', 'ID', 'geometry', 'radius', 'population', 'centroid', 'type', 'health']\n",
        "  hospfac = gpd.GeoDataFrame(columns = col)\n",
        "  EMSfac = gpd.GeoDataFrame(columns = col)\n",
        "  sheltfac = gpd.GeoDataFrame(columns = col)\n",
        "  firefac = gpd.GeoDataFrame(columns = col)\n",
        "\n",
        "\n",
        "  for i in col:\n",
        "    hospfac[i] = hosp_idx[i]\n",
        "    EMSfac[i] = EMS_idx[i]\n",
        "    sheltfac[i] = shelt_idx[i]\n",
        "    firefac[i] = fire_idx[i]\n",
        "\n",
        "  facdf = gpd.GeoDataFrame(columns = col)\n",
        "  fdf = facdf.append([hospfac, EMSfac, sheltfac, firefac])\n",
        "  alpha_fdf = fdf.sort_values('NAME').set_index('NAME')\n",
        "  fdf_idx = alpha_fdf.reset_index()\n",
        "  \n",
        "  fdf_info = gpd.sjoin(grid, fdf, how = 'inner', op = 'intersects')\n",
        "  \n",
        "  #wuabracadabra\n",
        "  return wuabracadabra(costs, budget, mindist, maxdist, buffer, fdf_idx, fdf_info, alpha_fdf, flooding, graph_type, threshold)\n",
        "\n",
        "  "
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "080v8dHKIExg",
        "colab_type": "code",
        "outputId": "11454fee-3f64-43f3-864e-535a5fbcf80d",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 1000
        }
      },
      "source": [
        "main(grid, flooding, [100, 50, 50, 25, 10, 15], 2000, 0, 1000, 500, 0, 6)"
      ],
      "execution_count": 40,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "/usr/local/lib/python3.6/dist-packages/geopandas/tools/sjoin.py:56: UserWarning: CRS of frames being joined does not match!({'init': 'epsg:3857', 'no_defs': True} != None)\n",
            "  '(%s != %s)' % (left_df.crs, right_df.crs))\n"
          ],
          "name": "stderr"
        },
        {
          "output_type": "execute_result",
          "data": {
            "text/html": [
              "<div>\n",
              "<style scoped>\n",
              "    .dataframe tbody tr th:only-of-type {\n",
              "        vertical-align: middle;\n",
              "    }\n",
              "\n",
              "    .dataframe tbody tr th {\n",
              "        vertical-align: top;\n",
              "    }\n",
              "\n",
              "    .dataframe thead th {\n",
              "        text-align: right;\n",
              "    }\n",
              "</style>\n",
              "<table border=\"1\" class=\"dataframe\">\n",
              "  <thead>\n",
              "    <tr style=\"text-align: right;\">\n",
              "      <th></th>\n",
              "      <th>ID</th>\n",
              "      <th>geometry</th>\n",
              "      <th>radius</th>\n",
              "      <th>population</th>\n",
              "      <th>centroid</th>\n",
              "      <th>type</th>\n",
              "      <th>health</th>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>NAME</th>\n",
              "      <th></th>\n",
              "      <th></th>\n",
              "      <th></th>\n",
              "      <th></th>\n",
              "      <th></th>\n",
              "      <th></th>\n",
              "      <th></th>\n",
              "    </tr>\n",
              "  </thead>\n",
              "  <tbody>\n",
              "    <tr>\n",
              "      <th>CHELSEA FIRE DEPARTMENT - HEADQUARTERS</th>\n",
              "      <td>10229411</td>\n",
              "      <td>POLYGON ((-7892352.722121374 5220227.991812181...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>712263.0</td>\n",
              "      <td>POINT (-7907352.722121376 5220227.99181218)</td>\n",
              "      <td>large EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>CHELSEA FIRE DEPARTMENT - HEADQUARTERS</th>\n",
              "      <td>10229411</td>\n",
              "      <td>POLYGON ((-7892352.722021262 5220227.991812181...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>712263.0</td>\n",
              "      <td>POINT (-7907352.722021261 5220227.99181218)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>REVERE FIRE DEPARTMENT HEADQUARTERS</th>\n",
              "      <td>10229190</td>\n",
              "      <td>POLYGON ((-7890069.723878649 5222672.037382862...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>602534.0</td>\n",
              "      <td>POINT (-7905069.723878649 5222672.037382863)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>REVERE FIRE DEPARTMENT HEADQUARTERS</th>\n",
              "      <td>10229190</td>\n",
              "      <td>POLYGON ((-7890069.723878649 5222672.037382862...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>602534.0</td>\n",
              "      <td>POINT (-7905069.723878649 5222672.037382863)</td>\n",
              "      <td>large EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>TUFTS MEDICAL CENTER</th>\n",
              "      <td>0001002111</td>\n",
              "      <td>POLYGON ((-7885808.293722065 5213537.72773091,...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1220580.0</td>\n",
              "      <td>POINT (-7910808.293722068 5213537.727730909)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>NORTH SHORE MEDICAL CENTER  - SALEM CAMPUS</th>\n",
              "      <td>0004501970</td>\n",
              "      <td>POLYGON ((-7868308.86978271 5237849.891362702,...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>447453.0</td>\n",
              "      <td>POINT (-7893308.869782709 5237849.891362701)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>CATALDO AMBULANCE SERVICE - REVERE BASE</th>\n",
              "      <td>10168057</td>\n",
              "      <td>POLYGON ((-7888851.126197943 5223863.39224383,...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>550189.0</td>\n",
              "      <td>POINT (-7903851.126197945 5223863.392243832)</td>\n",
              "      <td>large EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>REVERE FIRE DEPARTMENT ENGINE 5</th>\n",
              "      <td>10229188</td>\n",
              "      <td>POLYGON ((-7887695.070530999 5223984.511717383...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>502535.0</td>\n",
              "      <td>POINT (-7902695.070531 5223984.511717384)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>REVERE FIRE DEPARTMENT ENGINE 5</th>\n",
              "      <td>10229188</td>\n",
              "      <td>POLYGON ((-7892695.07063111 5223984.511717383,...</td>\n",
              "      <td>10.0</td>\n",
              "      <td>247013.0</td>\n",
              "      <td>POINT (-7902695.070631106 5223984.511717381)</td>\n",
              "      <td>small EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>VA BOSTON HEALTHCARE SYSTEM - JAMAICA PLAIN</th>\n",
              "      <td>0008502130</td>\n",
              "      <td>POLYGON ((-7890962.38616192 5210089.015941914,...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1197581.0</td>\n",
              "      <td>POINT (-7915962.386161918 5210089.015941912)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>ALDEN MILLS POINT OF PINES FIRE STATION - TRAINING FACILITY</th>\n",
              "      <td>10229189</td>\n",
              "      <td>POLYGON ((-7885144.77675957 5227098.8816414, -...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>374033.0</td>\n",
              "      <td>POINT (-7900144.77675957 5227098.881641401)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>ALDEN MILLS POINT OF PINES FIRE STATION - TRAINING FACILITY</th>\n",
              "      <td>10229189</td>\n",
              "      <td>POLYGON ((-7885144.776659458 5227098.8816414, ...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>374033.0</td>\n",
              "      <td>POINT (-7900144.776659461 5227098.881641401)</td>\n",
              "      <td>large EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>CHELSEA FIRE DEPARTMENT - ENGINE 3/LADDER 2</th>\n",
              "      <td>10229410</td>\n",
              "      <td>POLYGON ((-7896192.590091074 5220969.776106819...</td>\n",
              "      <td>10.0</td>\n",
              "      <td>348841.0</td>\n",
              "      <td>POINT (-7906192.590091074 5220969.776106817)</td>\n",
              "      <td>small EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>CHELSEA FIRE DEPARTMENT - ENGINE 3/LADDER 2</th>\n",
              "      <td>10229410</td>\n",
              "      <td>POLYGON ((-7891192.590091074 5220969.776106819...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>682214.0</td>\n",
              "      <td>POINT (-7906192.590091076 5220969.776106821)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>BOSTON FIRE DEPARTMENT DIVISION 1 DISTRICT 1 - ENGINE 56 LADDER 21</th>\n",
              "      <td>10229133</td>\n",
              "      <td>POLYGON ((-7889119.413036386 5219289.357784536...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>631975.0</td>\n",
              "      <td>POINT (-7904119.413036382 5219289.357784539)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>BOSTON FIRE DEPARTMENT DIVISION 1 DISTRICT 1 - ENGINE 56 LADDER 21</th>\n",
              "      <td>10229133</td>\n",
              "      <td>POLYGON ((-7894119.413036386 5219289.357648994...</td>\n",
              "      <td>10.0</td>\n",
              "      <td>290019.0</td>\n",
              "      <td>POINT (-7904119.413036385 5219289.357648992)</td>\n",
              "      <td>small EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>BETH ISRAEL DEACONESS MEDICAL CENTER - WEST CAMPUS</th>\n",
              "      <td>0001202215</td>\n",
              "      <td>POLYGON ((-7890839.934681447 5211730.385925083...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1200722.0</td>\n",
              "      <td>POINT (-7915839.934681444 5211730.385925081)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>CHILDREN'S HOSPITAL BOSTON</th>\n",
              "      <td>0007502115</td>\n",
              "      <td>POLYGON ((-7890528.240149499 5211624.968173864...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1201179.0</td>\n",
              "      <td>POINT (-7915528.240149499 5211624.968173862)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>REVERE FIRE DEPARTMENT ENGINE 1</th>\n",
              "      <td>10229186</td>\n",
              "      <td>POLYGON ((-7888118.510858756 5222329.798603789...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>546260.0</td>\n",
              "      <td>POINT (-7903118.510858758 5222329.798603791)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>REVERE FIRE DEPARTMENT ENGINE 1</th>\n",
              "      <td>10229186</td>\n",
              "      <td>POLYGON ((-7888118.510758645 5222329.798603789...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>546260.0</td>\n",
              "      <td>POINT (-7903118.510758645 5222329.798603791)</td>\n",
              "      <td>large EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>BETH ISRAEL DEACONESS MEDICAL CENTER - EAST CAMPUS</th>\n",
              "      <td>0000502215</td>\n",
              "      <td>POLYGON ((-7890383.524845394 5211971.34478531,...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1208502.0</td>\n",
              "      <td>POINT (-7915383.524845392 5211971.344785311)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>MOUNT AUBURN HOSPITAL</th>\n",
              "      <td>0002602138</td>\n",
              "      <td>POLYGON ((-7893634.053907593 5217183.583733298...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1203351.0</td>\n",
              "      <td>POINT (-7918634.053907591 5217183.583733297)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>BOSTON MEDICAL CENTER - EAST NEWTON CAMPUS</th>\n",
              "      <td>0001302118</td>\n",
              "      <td>POLYGON ((-7886543.002404683 5211594.849167342...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1216527.0</td>\n",
              "      <td>POINT (-7911543.002404679 5211594.849167342)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>CHELSEA FIRE DEPARTMENT - ENGINE 1</th>\n",
              "      <td>10229517</td>\n",
              "      <td>POLYGON ((-7897521.778174163 5221845.20015608,...</td>\n",
              "      <td>10.0</td>\n",
              "      <td>384541.0</td>\n",
              "      <td>POINT (-7907521.778174162 5221845.20015608)</td>\n",
              "      <td>small EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>CHELSEA FIRE DEPARTMENT - ENGINE 1</th>\n",
              "      <td>10229517</td>\n",
              "      <td>POLYGON ((-7892521.778174163 5221845.20015608,...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>701969.0</td>\n",
              "      <td>POINT (-7907521.778174164 5221845.200156081)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>MASSACHUSETTS GENERAL HOSPITAL</th>\n",
              "      <td>0007702114</td>\n",
              "      <td>POLYGON ((-7886398.287000467 5215450.87379117,...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1225365.0</td>\n",
              "      <td>POINT (-7911398.287000466 5215450.873791171)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>REVERE FIRE DEPARTMENT ENGINE 1</th>\n",
              "      <td>10229186</td>\n",
              "      <td>POLYGON ((-7888118.510858756 5222329.798603789...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>546260.0</td>\n",
              "      <td>POINT (-7903118.510858758 5222329.798603791)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>REVERE FIRE DEPARTMENT ENGINE 1</th>\n",
              "      <td>10229186</td>\n",
              "      <td>POLYGON ((-7888118.510758645 5222329.798603789...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>546260.0</td>\n",
              "      <td>POINT (-7903118.510758645 5222329.798603791)</td>\n",
              "      <td>large EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>ST. ELIZABETH'S MEDICAL CENTER</th>\n",
              "      <td>0002402135</td>\n",
              "      <td>POLYGON ((-7895214.790691317 5213356.977944344...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1187782.0</td>\n",
              "      <td>POINT (-7920214.790691316 5213356.977944341)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>SPAULDING REHABILITATION HOSPITAL</th>\n",
              "      <td>0010602129</td>\n",
              "      <td>POLYGON ((-7894144.684665511 5217847.985340781...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>758826.0</td>\n",
              "      <td>POINT (-7909144.684665511 5217847.985340782)</td>\n",
              "      <td>small hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>NAHANT FIRE DEPARTMENT</th>\n",
              "      <td>10218917</td>\n",
              "      <td>POLYGON ((-7886145.056201712 5224665.877675011...</td>\n",
              "      <td>10.0</td>\n",
              "      <td>128332.0</td>\n",
              "      <td>POINT (-7896145.056201709 5224665.87767501)</td>\n",
              "      <td>small EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>NAHANT FIRE DEPARTMENT</th>\n",
              "      <td>10218917</td>\n",
              "      <td>POLYGON ((-7881145.056201712 5224665.877675011...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>288549.0</td>\n",
              "      <td>POINT (-7896145.056201712 5224665.877675011)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>CAMBRIDGE FIRE DEPARTMENT - EAST CAMBRIDGE</th>\n",
              "      <td>10229218</td>\n",
              "      <td>POLYGON ((-7902449.745274012 5216744.257061245...</td>\n",
              "      <td>10.0</td>\n",
              "      <td>590987.0</td>\n",
              "      <td>POINT (-7912449.745274013 5216744.257061245)</td>\n",
              "      <td>small EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>CAMBRIDGE FIRE DEPARTMENT - EAST CAMBRIDGE</th>\n",
              "      <td>10229218</td>\n",
              "      <td>POLYGON ((-7897449.745274012 5216744.257061245...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>828514.0</td>\n",
              "      <td>POINT (-7912449.745274012 5216744.257061246)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>BRIGHAM AND WOMEN'S HOSPITAL</th>\n",
              "      <td>0007302115</td>\n",
              "      <td>POLYGON ((-7890572.767896872 5211429.195734219...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1204412.0</td>\n",
              "      <td>POINT (-7915572.767896872 5211429.195734218)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>REVERE FIRE DEPARTMENT ENGINE 5</th>\n",
              "      <td>10229188</td>\n",
              "      <td>POLYGON ((-7887695.070530999 5223984.511717383...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>502535.0</td>\n",
              "      <td>POINT (-7902695.070531 5223984.511717384)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>REVERE FIRE DEPARTMENT ENGINE 5</th>\n",
              "      <td>10229188</td>\n",
              "      <td>POLYGON ((-7892695.07063111 5223984.511717383,...</td>\n",
              "      <td>10.0</td>\n",
              "      <td>247013.0</td>\n",
              "      <td>POINT (-7902695.070631106 5223984.511717381)</td>\n",
              "      <td>small EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>HALLMARK HEALTH SYSTEM - MELROSE-WAKEFIELD HOSPITAL CAMPUS</th>\n",
              "      <td>0000202176</td>\n",
              "      <td>POLYGON ((-7885463.203304501 5230211.669397064...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1122502.0</td>\n",
              "      <td>POINT (-7910463.203304499 5230211.669397063)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>ALDEN MILLS POINT OF PINES FIRE STATION - TRAINING FACILITY</th>\n",
              "      <td>10229189</td>\n",
              "      <td>POLYGON ((-7885144.77675957 5227098.8816414, -...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>374033.0</td>\n",
              "      <td>POINT (-7900144.77675957 5227098.881641401)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>ALDEN MILLS POINT OF PINES FIRE STATION - TRAINING FACILITY</th>\n",
              "      <td>10229189</td>\n",
              "      <td>POLYGON ((-7885144.776659458 5227098.8816414, ...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>374033.0</td>\n",
              "      <td>POINT (-7900144.776659461 5227098.881641401)</td>\n",
              "      <td>large EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>BOSTON MEDICAL CENTER - MENINO CAMPUS</th>\n",
              "      <td>0001102118</td>\n",
              "      <td>POLYGON ((-7886832.43301289 5211248.486325764,...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1217093.0</td>\n",
              "      <td>POINT (-7911832.433012891 5211248.486325763)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>CAMBRIDGE HEALTH ALLIANCE - CAMBRIDGE CAMPUS</th>\n",
              "      <td>0006902139</td>\n",
              "      <td>POLYGON ((-7890272.205326794 5217243.857549599...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1222459.0</td>\n",
              "      <td>POINT (-7915272.205326792 5217243.857549597)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>REVERE FIRE DEPARTMENT HEADQUARTERS</th>\n",
              "      <td>10229190</td>\n",
              "      <td>POLYGON ((-7890069.723878649 5222672.037382862...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>602534.0</td>\n",
              "      <td>POINT (-7905069.723878649 5222672.037382863)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>REVERE FIRE DEPARTMENT HEADQUARTERS</th>\n",
              "      <td>10229190</td>\n",
              "      <td>POLYGON ((-7890069.723878649 5222672.037382862...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>602534.0</td>\n",
              "      <td>POINT (-7905069.723878649 5222672.037382863)</td>\n",
              "      <td>large EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>HALLMARK HEALTH SYSTEM - LAWRENCE MEMORIAL HOSPITAL CAMPUS</th>\n",
              "      <td>0003302155</td>\n",
              "      <td>POLYGON ((-7891084.837642399 5224992.248453597...</td>\n",
              "      <td>25.0</td>\n",
              "      <td>1213487.0</td>\n",
              "      <td>POINT (-7916084.837642399 5224992.248453596)</td>\n",
              "      <td>large hospital</td>\n",
              "      <td>2</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>NAHANT FIRE DEPARTMENT</th>\n",
              "      <td>10218917</td>\n",
              "      <td>POLYGON ((-7886145.056201712 5224665.877675011...</td>\n",
              "      <td>10.0</td>\n",
              "      <td>128332.0</td>\n",
              "      <td>POINT (-7896145.056201709 5224665.87767501)</td>\n",
              "      <td>small EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>NAHANT FIRE DEPARTMENT</th>\n",
              "      <td>10218917</td>\n",
              "      <td>POLYGON ((-7881145.056201712 5224665.877675011...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>288549.0</td>\n",
              "      <td>POINT (-7896145.056201712 5224665.877675011)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>CHELSEA FIRE DEPARTMENT - HEADQUARTERS</th>\n",
              "      <td>10229411</td>\n",
              "      <td>POLYGON ((-7892352.722121374 5220227.991812181...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>712263.0</td>\n",
              "      <td>POINT (-7907352.722121376 5220227.99181218)</td>\n",
              "      <td>large EMS</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "    <tr>\n",
              "      <th>CHELSEA FIRE DEPARTMENT - HEADQUARTERS</th>\n",
              "      <td>10229411</td>\n",
              "      <td>POLYGON ((-7892352.722021262 5220227.991812181...</td>\n",
              "      <td>15.0</td>\n",
              "      <td>712263.0</td>\n",
              "      <td>POINT (-7907352.722021261 5220227.99181218)</td>\n",
              "      <td>fire station</td>\n",
              "      <td>1</td>\n",
              "    </tr>\n",
              "  </tbody>\n",
              "</table>\n",
              "</div>"
            ],
            "text/plain": [
              "                                                            ID  ... health\n",
              "NAME                                                            ...       \n",
              "CHELSEA FIRE DEPARTMENT - HEADQUARTERS                10229411  ...      1\n",
              "CHELSEA FIRE DEPARTMENT - HEADQUARTERS                10229411  ...      1\n",
              "REVERE FIRE DEPARTMENT HEADQUARTERS                   10229190  ...      1\n",
              "REVERE FIRE DEPARTMENT HEADQUARTERS                   10229190  ...      1\n",
              "TUFTS MEDICAL CENTER                                0001002111  ...      2\n",
              "NORTH SHORE MEDICAL CENTER  - SALEM CAMPUS          0004501970  ...      2\n",
              "CATALDO AMBULANCE SERVICE - REVERE BASE               10168057  ...      1\n",
              "REVERE FIRE DEPARTMENT ENGINE 5                       10229188  ...      1\n",
              "REVERE FIRE DEPARTMENT ENGINE 5                       10229188  ...      1\n",
              "VA BOSTON HEALTHCARE SYSTEM - JAMAICA PLAIN         0008502130  ...      2\n",
              "ALDEN MILLS POINT OF PINES FIRE STATION - TRAIN...    10229189  ...      1\n",
              "ALDEN MILLS POINT OF PINES FIRE STATION - TRAIN...    10229189  ...      1\n",
              "CHELSEA FIRE DEPARTMENT - ENGINE 3/LADDER 2           10229410  ...      1\n",
              "CHELSEA FIRE DEPARTMENT - ENGINE 3/LADDER 2           10229410  ...      1\n",
              "BOSTON FIRE DEPARTMENT DIVISION 1 DISTRICT 1 - ...    10229133  ...      1\n",
              "BOSTON FIRE DEPARTMENT DIVISION 1 DISTRICT 1 - ...    10229133  ...      1\n",
              "BETH ISRAEL DEACONESS MEDICAL CENTER - WEST CAMPUS  0001202215  ...      2\n",
              "CHILDREN'S HOSPITAL BOSTON                          0007502115  ...      2\n",
              "REVERE FIRE DEPARTMENT ENGINE 1                       10229186  ...      1\n",
              "REVERE FIRE DEPARTMENT ENGINE 1                       10229186  ...      1\n",
              "BETH ISRAEL DEACONESS MEDICAL CENTER - EAST CAMPUS  0000502215  ...      2\n",
              "MOUNT AUBURN HOSPITAL                               0002602138  ...      2\n",
              "BOSTON MEDICAL CENTER - EAST NEWTON CAMPUS          0001302118  ...      2\n",
              "CHELSEA FIRE DEPARTMENT - ENGINE 1                    10229517  ...      1\n",
              "CHELSEA FIRE DEPARTMENT - ENGINE 1                    10229517  ...      1\n",
              "MASSACHUSETTS GENERAL HOSPITAL                      0007702114  ...      2\n",
              "REVERE FIRE DEPARTMENT ENGINE 1                       10229186  ...      1\n",
              "REVERE FIRE DEPARTMENT ENGINE 1                       10229186  ...      1\n",
              "ST. ELIZABETH'S MEDICAL CENTER                      0002402135  ...      2\n",
              "SPAULDING REHABILITATION HOSPITAL                   0010602129  ...      2\n",
              "NAHANT FIRE DEPARTMENT                                10218917  ...      1\n",
              "NAHANT FIRE DEPARTMENT                                10218917  ...      1\n",
              "CAMBRIDGE FIRE DEPARTMENT - EAST CAMBRIDGE            10229218  ...      1\n",
              "CAMBRIDGE FIRE DEPARTMENT - EAST CAMBRIDGE            10229218  ...      1\n",
              "BRIGHAM AND WOMEN'S HOSPITAL                        0007302115  ...      2\n",
              "REVERE FIRE DEPARTMENT ENGINE 5                       10229188  ...      1\n",
              "REVERE FIRE DEPARTMENT ENGINE 5                       10229188  ...      1\n",
              "HALLMARK HEALTH SYSTEM - MELROSE-WAKEFIELD HOSP...  0000202176  ...      2\n",
              "ALDEN MILLS POINT OF PINES FIRE STATION - TRAIN...    10229189  ...      1\n",
              "ALDEN MILLS POINT OF PINES FIRE STATION - TRAIN...    10229189  ...      1\n",
              "BOSTON MEDICAL CENTER - MENINO CAMPUS               0001102118  ...      2\n",
              "CAMBRIDGE HEALTH ALLIANCE - CAMBRIDGE CAMPUS        0006902139  ...      2\n",
              "REVERE FIRE DEPARTMENT HEADQUARTERS                   10229190  ...      1\n",
              "REVERE FIRE DEPARTMENT HEADQUARTERS                   10229190  ...      1\n",
              "HALLMARK HEALTH SYSTEM - LAWRENCE MEMORIAL HOSP...  0003302155  ...      2\n",
              "NAHANT FIRE DEPARTMENT                                10218917  ...      1\n",
              "NAHANT FIRE DEPARTMENT                                10218917  ...      1\n",
              "CHELSEA FIRE DEPARTMENT - HEADQUARTERS                10229411  ...      1\n",
              "CHELSEA FIRE DEPARTMENT - HEADQUARTERS                10229411  ...      1\n",
              "\n",
              "[49 rows x 7 columns]"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 40
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "mzXToFKmKp1f",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        ""
      ],
      "execution_count": 0,
      "outputs": []
    }
  ]
}