service_test.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Copyright © 2021 Helmholtz Centre Potsdam GFZ German Research Centre for Geosciences, Potsdam, Germany\n",
    "#\n",
    "# Licensed under the Apache License, Version 2.0 (the \"License\"); you may not use this file except in compliance with the License. You may obtain a copy of the License at\n",
    "# \n",
    "# https://www.apache.org/licenses/LICENSE-2.0\n",
    "# \n",
    "# Unless required by applicable law or agreed to in writing, software distributed under the License is distributed on an \"AS IS\" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. See the License for the specific language governing permissions and limitations under the License.\n",
    "%load_ext autoreload\n",
    "%autoreload 2\n",
    "import warnings; warnings.filterwarnings(\"ignore\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import os\n",
    "import geopandas as gp\n",
    "from shapely.geometry import Polygon\n",
    "#import matplotlib\n",
    "#import matplotlib.pyplot as plt\n",
    "import expo_nrml as nrml\n",
    "import assetmaster"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### define a ROI"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "def get_roi(latmax, latmin, lonmax, lonmin):\n",
    "    lat_point_list = [latmin,latmax , latmax, latmin, latmin]\n",
    "    lon_point_list = [lonmin, lonmin,lonmax,lonmax,lonmin]\n",
    "    roi_geom = Polygon(zip(lon_point_list, lat_point_list))\n",
    "    crs = {'init': 'epsg:4326'}\n",
    "    roi = gp.GeoDataFrame(index=[0], crs=crs, geometry=[roi_geom])  \n",
    "    return (roi)\n",
    "\n",
    "latmax=-32.6\n",
    "latmin=-33.6\n",
    "lonmax=-71.2\n",
    "lonmin=-72.2\n",
    "\n",
    "roi = get_roi(latmax, latmin, lonmax, lonmin)\n",
    "#roi"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Instantiate the class and run the main function"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "#create an argument list\n",
    "from argparse import Namespace\n",
    "args = Namespace(assettype='res', latmax=latmax, latmin=latmin, \n",
    "                 lonmax=lonmax, lonmin=lonmin, schema='SARA_v1.0',\n",
    "                querymode='intersects')\n",
    "\n",
    "#instantiate the class and run the main function\n",
    "model = assetmaster.Main(args)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "model._check_longitude()\n",
    "model._compute_roi()\n",
    "if (model._check_schema()):\n",
    "    foldername = os.path.join(model.folder,\"schemas/{}\".format(model.schema))\n",
    "    model.path_expo_dict = foldername\n",
    "    model.path_metadatefile = foldername\n",
    "    model.path_infile = foldername\n",
    "    model.in_file = \"{}_data.gpkg\".format(model.schema)\n",
    "    model.dict_file = \"{}_prop.csv\".format(model.schema)\n",
    "    model.metadata_file = \"{}_meta.json\".format(model.schema)\n",
    "else:\n",
    "    raise Exception (\"schema {} not supported\".format(model.schema))\n",
    "\n",
    "#read the exposure model metadata (with the list of taxonomies)\n",
    "model.metadata = nrml.read_metadata(os.path.join(model.path_metadatefile,model.metadata_file))\n",
    "model.taxonomies = model.metadata['taxonomies']\n",
    "#print(self.taxonomies)\n",
    "\n",
    "#get a dataframe with the basic properties of the buildings\n",
    "model.dicts = nrml.load_expo_dicts(os.path.join(model.path_expo_dict,model.dict_file))\n",
    "btypes = model.dicts.btype"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'CR-LDUAL-DUC-H4-7',\n",
       " 'CR-LDUAL-DUC-H8-19',\n",
       " 'CR-LFINF-DNO-H1-3',\n",
       " 'CR-LFINF-DUC-H1-3',\n",
       " 'CR-LFINF-DUC-H4-7',\n",
       " 'CR-LFLS-DNO-H1-3',\n",
       " 'CR-LFLS-DUC-H1-3',\n",
       " 'CR-LFLS-DUC-H4-7',\n",
       " 'CR-LFM-DNO-H1-3',\n",
       " 'CR-LFM-DNO-H4-7',\n",
       " 'CR-LFM-DNO-SOS-H1-3',\n",
       " 'CR-LFM-DNO-SOS-H4-7',\n",
       " 'CR-LFM-DUC-H1-3',\n",
       " 'CR-LFM-DUC-H4-7',\n",
       " 'CR-LWAL-DNO-H1-3',\n",
       " 'CR-LWAL-DNO-H4-7',\n",
       " 'CR-LWAL-DUC-H1-3',\n",
       " 'CR-LWAL-DUC-H4-7',\n",
       " 'CR-LWAL-DUC-H8-19',\n",
       " 'CR-PC-LWAL-H1-3',\n",
       " 'ER-ETR-H1-2',\n",
       " 'MCF-DNO-H1-3',\n",
       " 'MCF-DUC-H1-3',\n",
       " 'MR-DNO-H1-3',\n",
       " 'MR-DUC-H1-3',\n",
       " 'MUR-ADO-H1-2',\n",
       " 'MUR-H1-3',\n",
       " 'MUR-STDRE-H1-2',\n",
       " 'MUR-STRUB-H1-2',\n",
       " 'S-LFM-H4-7',\n",
       " 'UNK',\n",
       " 'W-WBB-H1',\n",
       " 'W-WHE-H1-3',\n",
       " 'W-WLI-H1-3',\n",
       " 'W-WS-H1-2'}"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "set(btypes)\n",
    "#model.taxonomies"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'CR-LDUAL-DUC-H4-7',\n",
       " 'CR-LDUAL-DUC-H8-19',\n",
       " 'CR-LFINF-DNO-H1-3',\n",
       " 'CR-LFINF-DUC-H1-3',\n",
       " 'CR-LFINF-DUC-H4-7',\n",
       " 'CR-LFLS-DNO-H1-3',\n",
       " 'CR-LFLS-DUC-H1-3',\n",
       " 'CR-LFLS-DUC-H4-7',\n",
       " 'CR-LFM-DNO-H1-3',\n",
       " 'CR-LFM-DNO-H4-7',\n",
       " 'CR-LFM-DNO-SOS-H1-3',\n",
       " 'CR-LFM-DNO-SOS-H4-7',\n",
       " 'CR-LFM-DUC-H1-3',\n",
       " 'CR-LFM-DUC-H4-7',\n",
       " 'CR-LWAL-DNO-H1-3',\n",
       " 'CR-LWAL-DNO-H4-7',\n",
       " 'CR-LWAL-DUC-H1-3',\n",
       " 'CR-LWAL-DUC-H4-7',\n",
       " 'CR-LWAL-DUC-H8-19',\n",
       " 'CR-PC-LWAL-H1-3',\n",
       " 'ER-ETR-H1-2',\n",
       " 'MCF-DNO-H1-3',\n",
       " 'MCF-DUC-H1-3',\n",
       " 'MR-DNO-H1-3',\n",
       " 'MR-DUC-H1-3',\n",
       " 'MUR-ADO-H1-2',\n",
       " 'MUR-H1-3',\n",
       " 'MUR-STDRE-H1-2',\n",
       " 'MUR-STRUB-H1-2',\n",
       " 'S-LFM-H4-7',\n",
       " 'UNK',\n",
       " 'W-WBB-H1',\n",
       " 'W-WHE-H1-3',\n",
       " 'W-WLI-H1-3',\n",
       " 'W-WS-H1-2'}"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "set(model.taxonomies)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "set()"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "set(model.taxonomies).difference(set(btypes))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model.run()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Check the results"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "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>name</th>\n",
       "      <th>gid</th>\n",
       "      <th>expo</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <td>0</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>0 1</td>\n",
       "      <td>{\"id\":{\"0\":\"AREA # 0 1\",\"1\":\"AREA # 0 1\",\"2\":\"...</td>\n",
       "      <td>POLYGON ((-71.4765109995829 -33.01033200032151...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>0 10</td>\n",
       "      <td>{\"id\":{\"0\":\"AREA # 0 10\",\"1\":\"AREA # 0 10\",\"2\"...</td>\n",
       "      <td>POLYGON ((-71.48944699987908 -32.9985789997210...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>2</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>0 100</td>\n",
       "      <td>{\"id\":{\"0\":\"AREA # 0 100\",\"1\":\"AREA # 0 100\",\"...</td>\n",
       "      <td>(POLYGON ((-71.5503686401006 -32.9586378602489...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>3</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>0 1001</td>\n",
       "      <td>{\"id\":{\"0\":\"AREA # 0 1001\",\"1\":\"AREA # 0 1001\"...</td>\n",
       "      <td>POLYGON ((-71.51637599975254 -32.9961340002838...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>4</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>0 1005</td>\n",
       "      <td>{\"id\":{\"0\":\"AREA # 0 1005\",\"1\":\"AREA # 0 1005\"...</td>\n",
       "      <td>POLYGON ((-71.55798300021524 -33.0261810000530...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>2901</td>\n",
       "      <td>VALPARAISO</td>\n",
       "      <td>1 991</td>\n",
       "      <td>{\"id\":{\"0\":\"AREA # 1 991\",\"1\":\"AREA # 1 991\",\"...</td>\n",
       "      <td>POLYGON ((-71.46055899959811 -33.1811739997802...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>2902</td>\n",
       "      <td>VALPARAISO</td>\n",
       "      <td>1 992</td>\n",
       "      <td>{\"id\":{\"0\":\"AREA # 1 992\",\"1\":\"AREA # 1 992\",\"...</td>\n",
       "      <td>POLYGON ((-71.57578299998545 -33.1330710001265...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>2903</td>\n",
       "      <td>VALPARAISO</td>\n",
       "      <td>1 993</td>\n",
       "      <td>{\"id\":{\"0\":\"AREA # 1 993\",\"1\":\"AREA # 1 993\",\"...</td>\n",
       "      <td>POLYGON ((-71.59109300004189 -33.0657669999470...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>2904</td>\n",
       "      <td>VALPARAISO</td>\n",
       "      <td>1 995</td>\n",
       "      <td>{\"id\":{\"0\":\"AREA # 1 995\",\"1\":\"AREA # 1 995\",\"...</td>\n",
       "      <td>POLYGON ((-71.67574599987796 -33.1487259999049...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>2905</td>\n",
       "      <td>VALPARAISO</td>\n",
       "      <td>1 999</td>\n",
       "      <td>{\"id\":{\"0\":\"AREA # 1 999\",\"1\":\"AREA # 1 999\",\"...</td>\n",
       "      <td>POLYGON ((-71.64934299997941 -33.0467499998126...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2906 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "              name     gid                                               expo  \\\n",
       "0     VINA DEL MAR     0 1  {\"id\":{\"0\":\"AREA # 0 1\",\"1\":\"AREA # 0 1\",\"2\":\"...   \n",
       "1     VINA DEL MAR    0 10  {\"id\":{\"0\":\"AREA # 0 10\",\"1\":\"AREA # 0 10\",\"2\"...   \n",
       "2     VINA DEL MAR   0 100  {\"id\":{\"0\":\"AREA # 0 100\",\"1\":\"AREA # 0 100\",\"...   \n",
       "3     VINA DEL MAR  0 1001  {\"id\":{\"0\":\"AREA # 0 1001\",\"1\":\"AREA # 0 1001\"...   \n",
       "4     VINA DEL MAR  0 1005  {\"id\":{\"0\":\"AREA # 0 1005\",\"1\":\"AREA # 0 1005\"...   \n",
       "...            ...     ...                                                ...   \n",
       "2901    VALPARAISO   1 991  {\"id\":{\"0\":\"AREA # 1 991\",\"1\":\"AREA # 1 991\",\"...   \n",
       "2902    VALPARAISO   1 992  {\"id\":{\"0\":\"AREA # 1 992\",\"1\":\"AREA # 1 992\",\"...   \n",
       "2903    VALPARAISO   1 993  {\"id\":{\"0\":\"AREA # 1 993\",\"1\":\"AREA # 1 993\",\"...   \n",
       "2904    VALPARAISO   1 995  {\"id\":{\"0\":\"AREA # 1 995\",\"1\":\"AREA # 1 995\",\"...   \n",
       "2905    VALPARAISO   1 999  {\"id\":{\"0\":\"AREA # 1 999\",\"1\":\"AREA # 1 999\",\"...   \n",
       "\n",
       "                                               geometry  \n",
       "0     POLYGON ((-71.4765109995829 -33.01033200032151...  \n",
       "1     POLYGON ((-71.48944699987908 -32.9985789997210...  \n",
       "2     (POLYGON ((-71.5503686401006 -32.9586378602489...  \n",
       "3     POLYGON ((-71.51637599975254 -32.9961340002838...  \n",
       "4     POLYGON ((-71.55798300021524 -33.0261810000530...  \n",
       "...                                                 ...  \n",
       "2901  POLYGON ((-71.46055899959811 -33.1811739997802...  \n",
       "2902  POLYGON ((-71.57578299998545 -33.1330710001265...  \n",
       "2903  POLYGON ((-71.59109300004189 -33.0657669999470...  \n",
       "2904  POLYGON ((-71.67574599987796 -33.1487259999049...  \n",
       "2905  POLYGON ((-71.64934299997941 -33.0467499998126...  \n",
       "\n",
       "[2906 rows x 4 columns]"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "#check the results\n",
    "model.query_result"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read a single exposure model and convert it into a pandas structure"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "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>Region</th>\n",
       "      <th>Taxonomy</th>\n",
       "      <th>Dwellings</th>\n",
       "      <th>Buildings</th>\n",
       "      <th>Repl-cost-USD-bdg</th>\n",
       "      <th>Population</th>\n",
       "      <th>name</th>\n",
       "      <th>Damage</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <td>0</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>CR-LWAL-DNO-H1-3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>1</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>CR-LWAL-DNO-H4-7</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>2</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>CR-LWAL-DUC-H1-3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>3</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>CR-LWAL-DUC-H4-7</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>4</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>CR-LWAL-DUC-H8-19</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>CR-PC-LWAL-H1-3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>6</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>ER-ETR-H1-2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>7</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>MCF-DNO-H1-3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>8</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>MCF-DUC-H1-3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>9</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>MR-DNO-H1-3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>10</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>MR-DUC-H1-3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>11</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>MUR-ADO-H1-2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>12</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>MUR-H1-3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>13</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>MUR-STDRE-H1-2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>14</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>UNK</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>15</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>W-WLI-H1-3</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>16</td>\n",
       "      <td>AREA # 0 1</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>W-WS-H1-2</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>15</td>\n",
       "      <td>VINA DEL MAR</td>\n",
       "      <td>D0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "            id        Region           Taxonomy  Dwellings  Buildings  \\\n",
       "0   AREA # 0 1  VINA DEL MAR   CR-LWAL-DNO-H1-3          0          0   \n",
       "1   AREA # 0 1  VINA DEL MAR   CR-LWAL-DNO-H4-7          0          0   \n",
       "2   AREA # 0 1  VINA DEL MAR   CR-LWAL-DUC-H1-3          0          0   \n",
       "3   AREA # 0 1  VINA DEL MAR   CR-LWAL-DUC-H4-7          0          0   \n",
       "4   AREA # 0 1  VINA DEL MAR  CR-LWAL-DUC-H8-19          0          0   \n",
       "5   AREA # 0 1  VINA DEL MAR    CR-PC-LWAL-H1-3          0          0   \n",
       "6   AREA # 0 1  VINA DEL MAR        ER-ETR-H1-2          0          0   \n",
       "7   AREA # 0 1  VINA DEL MAR       MCF-DNO-H1-3          0          0   \n",
       "8   AREA # 0 1  VINA DEL MAR       MCF-DUC-H1-3          0          0   \n",
       "9   AREA # 0 1  VINA DEL MAR        MR-DNO-H1-3          0          0   \n",
       "10  AREA # 0 1  VINA DEL MAR        MR-DUC-H1-3          0          0   \n",
       "11  AREA # 0 1  VINA DEL MAR       MUR-ADO-H1-2          0          0   \n",
       "12  AREA # 0 1  VINA DEL MAR           MUR-H1-3          0          0   \n",
       "13  AREA # 0 1  VINA DEL MAR     MUR-STDRE-H1-2          0          0   \n",
       "14  AREA # 0 1  VINA DEL MAR                UNK          0          0   \n",
       "15  AREA # 0 1  VINA DEL MAR         W-WLI-H1-3          0          0   \n",
       "16  AREA # 0 1  VINA DEL MAR          W-WS-H1-2          0          0   \n",
       "\n",
       "    Repl-cost-USD-bdg  Population          name Damage  \n",
       "0                   0          15  VINA DEL MAR     D0  \n",
       "1                   0          15  VINA DEL MAR     D0  \n",
       "2                   0          15  VINA DEL MAR     D0  \n",
       "3                   0          15  VINA DEL MAR     D0  \n",
       "4                   0          15  VINA DEL MAR     D0  \n",
       "5                   0          15  VINA DEL MAR     D0  \n",
       "6                   0          15  VINA DEL MAR     D0  \n",
       "7                   0          15  VINA DEL MAR     D0  \n",
       "8                   0          15  VINA DEL MAR     D0  \n",
       "9                   0          15  VINA DEL MAR     D0  \n",
       "10                  0          15  VINA DEL MAR     D0  \n",
       "11                  0          15  VINA DEL MAR     D0  \n",
       "12                  0          15  VINA DEL MAR     D0  \n",
       "13                  0          15  VINA DEL MAR     D0  \n",
       "14                  0          15  VINA DEL MAR     D0  \n",
       "15                  0          15  VINA DEL MAR     D0  \n",
       "16                  0          15  VINA DEL MAR     D0  "
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd.read_json(model.query_result.iloc[0].expo)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Visualize the ROI and the selected geocells"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "ename": "ModuleNotFoundError",
     "evalue": "No module named 'folium'",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mModuleNotFoundError\u001b[0m                       Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-21-ff95877a3921>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0;32mimport\u001b[0m \u001b[0mfolium\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m      2\u001b[0m \u001b[0;32mimport\u001b[0m \u001b[0mpandas\u001b[0m \u001b[0;32mas\u001b[0m \u001b[0mpd\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      3\u001b[0m \u001b[0;32mimport\u001b[0m \u001b[0mjson\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      4\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      5\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mget_pivoted_df\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdf\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvalues\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'Buildings'\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mcolumns\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'index'\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mindex\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'Region'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mModuleNotFoundError\u001b[0m: No module named 'folium'"
     ]
    }
   ],
   "source": [
    "import folium\n",
    "import pandas as pd \n",
    "import json\n",
    "\n",
    "def get_pivoted_df(df, values='Buildings',columns='index',index='Region'):\n",
    "    test_1=df.expo.apply(lambda x: pd.DataFrame(json.loads(x)).reset_index().pivot(values=values,columns=columns,index=index).reset_index())\n",
    "    test2 = pd.concat(list(test_1),join='outer',ignore_index=True)\n",
    "    return (test2)\n",
    "\n",
    "def plot_expo(model):#,taxonomies):\n",
    "    piv_bdgs = get_pivoted_df(model,values='Buildings',columns='Taxonomy',index='Region')\n",
    "    taxonomies  = list(set(piv_bdgs.columns).difference(['Region']))\n",
    "\n",
    "    #tmp = model.copy()\n",
    "    piv_bdgs['total'] = (piv_bdgs[taxonomies].sum(axis=1))\n",
    "    #change the significant digits\n",
    "    piv_bdgs[taxonomies] = piv_bdgs[taxonomies].divide(piv_bdgs.total,axis='index')*100.\n",
    "\n",
    "    li = ['name','total']\n",
    "    li.extend(taxonomies)\n",
    "    al = ['geocell','tot buildings'] \n",
    "    al.extend(taxonomies)\n",
    "\n",
    "    m = folium.Map([-33.08,-71.6], zoom_start=9, tiles='cartodbpositron')\n",
    "    folium.GeoJson(roi).add_to(m)\n",
    "    folium.GeoJson(model.set_index('name').join(piv_bdgs.set_index('Region')).reset_index(),name=\"Exposure model\",\n",
    "                   style_function=lambda x: {\"weight\":2, 'color':'black','fillColor':'white', 'fillOpacity':0.2},\n",
    "                  highlight_function=lambda x: {'weight':3, 'color':'black'},\n",
    "                   smooth_factor=2.0,\n",
    "                  tooltip=folium.features.GeoJsonTooltip(fields=li,aliases=al,sticky=True,labels=True,localize=True)\n",
    "                  ).add_to(m)\n",
    "    return(m)\n",
    "\n",
    "#test\n",
    "plot_expo(model.query_result)\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}