assetmaster_modelcreator.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",
    "\n",
    "%load_ext autoreload\n",
    "%autoreload 2\n",
    "import warnings; warnings.filterwarnings(\"ignore\")\n",
    "import pandas as pd\n",
    "import geopandas as gpd\n",
    "%matplotlib inline\n",
    "import h5py\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from matplotlib.colors import LogNorm\n",
    "import os\n",
    "import folium"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read  CHILE GEM SARA Model (csv format)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "path = '/Users/pittore/Documents/PROJECTS/BMBF/RIESGOS/data/GEM_sara_chile/chl-l3-exposure'\n",
    "infile = 'chl-l3-exposure_total_mod_tax.csv'\n",
    "sara_data = pd.read_csv(os.path.join(path,infile))\n",
    "\n",
    "#clean 'name' and use it as ID\n",
    "sara_data['name'] = sara_data.Region.str.normalize('NFKD').str.encode('ascii', errors='ignore').str.decode('utf-8')\n",
    "#save all names in an array\n",
    "sara_names = sara_data['name'].unique()\n",
    "\n",
    "#group data based on 'name'\n",
    "datagrp = sara_data.groupby('name')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 124,
   "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",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <td>5709</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>CR-PC-LWAL-H1-3</td>\n",
       "      <td>78.0</td>\n",
       "      <td>15.6</td>\n",
       "      <td>360000</td>\n",
       "      <td>308.7</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5710</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>CR-LWAL-DNO-H1-3</td>\n",
       "      <td>35.6</td>\n",
       "      <td>8.9</td>\n",
       "      <td>288000</td>\n",
       "      <td>140.8</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5711</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>CR-LWAL-DNO-H4-7</td>\n",
       "      <td>1.7</td>\n",
       "      <td>0.1</td>\n",
       "      <td>1080000</td>\n",
       "      <td>6.7</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5712</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>CR-LWAL-DUC-H1-3</td>\n",
       "      <td>6.7</td>\n",
       "      <td>1.6</td>\n",
       "      <td>336000</td>\n",
       "      <td>26.4</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5713</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>CR-LWAL-DUC-H4-7</td>\n",
       "      <td>1.4</td>\n",
       "      <td>0.1</td>\n",
       "      <td>1260000</td>\n",
       "      <td>5.5</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5714</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>CR-LWAL-DUC-H8-19</td>\n",
       "      <td>1.2</td>\n",
       "      <td>0.0</td>\n",
       "      <td>4032000</td>\n",
       "      <td>4.7</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5715</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>ER-ETR-H1-2</td>\n",
       "      <td>0.6</td>\n",
       "      <td>0.5</td>\n",
       "      <td>43750</td>\n",
       "      <td>2.4</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5716</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>MCF-DNO-H1</td>\n",
       "      <td>3.8</td>\n",
       "      <td>3.7</td>\n",
       "      <td>35000</td>\n",
       "      <td>13.3</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5717</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>MCF-DNO-H1-3</td>\n",
       "      <td>129.4</td>\n",
       "      <td>86.2</td>\n",
       "      <td>94500</td>\n",
       "      <td>508.8</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5718</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>MCF-DUC-H1-3</td>\n",
       "      <td>88.0</td>\n",
       "      <td>22.0</td>\n",
       "      <td>288000</td>\n",
       "      <td>346.5</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5719</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>MR-DNO-H1-3</td>\n",
       "      <td>114.7</td>\n",
       "      <td>23.0</td>\n",
       "      <td>360000</td>\n",
       "      <td>447.4</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5720</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>MR-DUC-H1-3</td>\n",
       "      <td>27.2</td>\n",
       "      <td>5.5</td>\n",
       "      <td>420000</td>\n",
       "      <td>104.4</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5721</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>MUR-ADO-H1-2</td>\n",
       "      <td>2.4</td>\n",
       "      <td>1.9</td>\n",
       "      <td>43750</td>\n",
       "      <td>9.5</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5722</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>MUR-STDRE-H1-2</td>\n",
       "      <td>2.0</td>\n",
       "      <td>1.7</td>\n",
       "      <td>43750</td>\n",
       "      <td>7.8</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5723</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>MUR-H1-3</td>\n",
       "      <td>8.2</td>\n",
       "      <td>5.5</td>\n",
       "      <td>52500</td>\n",
       "      <td>32.4</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5724</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>UNK</td>\n",
       "      <td>359.8</td>\n",
       "      <td>359.8</td>\n",
       "      <td>35000</td>\n",
       "      <td>1411.3</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5725</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>W-WLI-H1</td>\n",
       "      <td>146.9</td>\n",
       "      <td>146.9</td>\n",
       "      <td>72000</td>\n",
       "      <td>567.9</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5726</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>W-WLI-H1-3</td>\n",
       "      <td>3530.9</td>\n",
       "      <td>2353.9</td>\n",
       "      <td>108000</td>\n",
       "      <td>13692.3</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5727</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>W-WS-H1</td>\n",
       "      <td>169.3</td>\n",
       "      <td>169.3</td>\n",
       "      <td>35000</td>\n",
       "      <td>657.1</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5728</td>\n",
       "      <td>AREA # 11201</td>\n",
       "      <td>Aisén</td>\n",
       "      <td>W-WS-H1-2</td>\n",
       "      <td>1073.1</td>\n",
       "      <td>858.5</td>\n",
       "      <td>43750</td>\n",
       "      <td>4058.8</td>\n",
       "      <td>Aisen</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                id Region           Taxonomy  Dwellings  Buildings  \\\n",
       "5709  AREA # 11201  Aisén    CR-PC-LWAL-H1-3       78.0       15.6   \n",
       "5710  AREA # 11201  Aisén   CR-LWAL-DNO-H1-3       35.6        8.9   \n",
       "5711  AREA # 11201  Aisén   CR-LWAL-DNO-H4-7        1.7        0.1   \n",
       "5712  AREA # 11201  Aisén   CR-LWAL-DUC-H1-3        6.7        1.6   \n",
       "5713  AREA # 11201  Aisén   CR-LWAL-DUC-H4-7        1.4        0.1   \n",
       "5714  AREA # 11201  Aisén  CR-LWAL-DUC-H8-19        1.2        0.0   \n",
       "5715  AREA # 11201  Aisén        ER-ETR-H1-2        0.6        0.5   \n",
       "5716  AREA # 11201  Aisén         MCF-DNO-H1        3.8        3.7   \n",
       "5717  AREA # 11201  Aisén       MCF-DNO-H1-3      129.4       86.2   \n",
       "5718  AREA # 11201  Aisén       MCF-DUC-H1-3       88.0       22.0   \n",
       "5719  AREA # 11201  Aisén        MR-DNO-H1-3      114.7       23.0   \n",
       "5720  AREA # 11201  Aisén        MR-DUC-H1-3       27.2        5.5   \n",
       "5721  AREA # 11201  Aisén       MUR-ADO-H1-2        2.4        1.9   \n",
       "5722  AREA # 11201  Aisén     MUR-STDRE-H1-2        2.0        1.7   \n",
       "5723  AREA # 11201  Aisén           MUR-H1-3        8.2        5.5   \n",
       "5724  AREA # 11201  Aisén                UNK      359.8      359.8   \n",
       "5725  AREA # 11201  Aisén           W-WLI-H1      146.9      146.9   \n",
       "5726  AREA # 11201  Aisén         W-WLI-H1-3     3530.9     2353.9   \n",
       "5727  AREA # 11201  Aisén            W-WS-H1      169.3      169.3   \n",
       "5728  AREA # 11201  Aisén          W-WS-H1-2     1073.1      858.5   \n",
       "\n",
       "      Repl-cost-USD-bdg  Population   name  \n",
       "5709             360000       308.7  Aisen  \n",
       "5710             288000       140.8  Aisen  \n",
       "5711            1080000         6.7  Aisen  \n",
       "5712             336000        26.4  Aisen  \n",
       "5713            1260000         5.5  Aisen  \n",
       "5714            4032000         4.7  Aisen  \n",
       "5715              43750         2.4  Aisen  \n",
       "5716              35000        13.3  Aisen  \n",
       "5717              94500       508.8  Aisen  \n",
       "5718             288000       346.5  Aisen  \n",
       "5719             360000       447.4  Aisen  \n",
       "5720             420000       104.4  Aisen  \n",
       "5721              43750         9.5  Aisen  \n",
       "5722              43750         7.8  Aisen  \n",
       "5723              52500        32.4  Aisen  \n",
       "5724              35000      1411.3  Aisen  \n",
       "5725              72000       567.9  Aisen  \n",
       "5726             108000     13692.3  Aisen  \n",
       "5727              35000       657.1  Aisen  \n",
       "5728              43750      4058.8  Aisen  "
      ]
     },
     "execution_count": 124,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "a = datagrp.get_group(\"Aisen\")\n",
    "a"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 145,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "tax_lut = {\"W-WS-H1\":\"W-WS-H1-2\",\n",
    "            \"W-WLI-H1\":\"W-WLI-H1-3\",\n",
    "            \"MCF-DNO-H1\":\"MCF-DNO-H1-3\",\n",
    "            \"ER-ETR-H1\":\"ER-ETR-H1-2\",\n",
    "            \"MUR-H1\":\"MUR-H1-3\",\n",
    "            \"MUR-ADO-H1\":\"MUR-ADO-H1-2\",\n",
    "            \"MCF-DNO-H1\":\"MCF-DNO-H1-3\"}\n",
    "\n",
    "#tax_props = {\"W-WS-H1\":{\"Dwellings\": 1,\"Repl-cost-USD-bdg\":35000,},\n",
    "#            \"W-WLI-H1\":{}}\n",
    "\n",
    "def updategroup(grp):\n",
    "    group_tax = set(grp.Taxonomy)\n",
    "    droplist = []\n",
    "    if 1:\n",
    "        for i,item in grp.iterrows():\n",
    "            if item.Taxonomy in tax_lut.keys():\n",
    "                if tax_lut[item.Taxonomy] in group_tax:\n",
    "                    #print (tax_lut[item.Taxonomy])\n",
    "                    tg=grp[grp.Taxonomy==tax_lut[item.Taxonomy]]\n",
    "                    #print(tg.index)\n",
    "                    grp.loc[tg.index,\"Dwellings\"] = tg.Dwellings + item.Dwellings\n",
    "                    grp.loc[tg.index,\"Buildings\"] = tg.Buildings + item.Buildings\n",
    "                    grp.loc[tg.index,\"Population\"] = tg.Population + item.Population\n",
    "                    #put index in droplist\n",
    "                    droplist.append(i)\n",
    "                else:\n",
    "                    #just change the name of the taxonomy...\n",
    "                    grp.loc[i,\"Taxonomy\"]=tax_lut[item.Taxonomy]\n",
    "\n",
    "    grp.drop(droplist,inplace=True)\n",
    "    return(grp)\n",
    "\n",
    "#testg = datagrp.get_group(\"Aisen\")\n",
    "#updategroup(testg)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read Chile admin_3 areas from GADM"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "302"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "chile_map = '/Users/pittore/Documents/PROJECTS/BMBF/RIESGOS/data/GADM_chile/gadm36_CHL_3.json'\n",
    "chile_map_gpd = gpd.read_file(chile_map)\n",
    "chile_map_gpd['name'] = chile_map_gpd.NAME_3.str.normalize('NFKD').str.encode('ascii', errors='ignore').str.decode('utf-8')\n",
    "len(chile_map_gpd)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Associate exposure model for each admin area and change geometry type to MultiPolygon"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 146,
   "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>geometry</th>\n",
       "      <th>expo</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <td>0</td>\n",
       "      <td>Aisen</td>\n",
       "      <td>CHL.1.1.1_1</td>\n",
       "      <td>(POLYGON ((-73.49375152587885 -44.965415954589...</td>\n",
       "      <td>{\"id\":{\"5709\":\"AREA # 11201\",\"5710\":\"AREA # 11...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>1</td>\n",
       "      <td>Cisnes</td>\n",
       "      <td>CHL.1.1.2_1</td>\n",
       "      <td>(POLYGON ((-73.31124877929682 -45.111248016357...</td>\n",
       "      <td>{\"id\":{\"5729\":\"AREA # 11202\",\"5730\":\"AREA # 11...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>2</td>\n",
       "      <td>Guaitecas</td>\n",
       "      <td>CHL.1.1.3_1</td>\n",
       "      <td>(POLYGON ((-73.49153137207031 -43.711528778076...</td>\n",
       "      <td>{\"id\":{\"5747\":\"AREA # 11203\",\"5748\":\"AREA # 11...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>3</td>\n",
       "      <td>Cochrane</td>\n",
       "      <td>CHL.1.2.1_1</td>\n",
       "      <td>(POLYGON ((-72.42082214355469 -47.743259429931...</td>\n",
       "      <td>{\"id\":{\"5757\":\"AREA # 11301\",\"5758\":\"AREA # 11...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>4</td>\n",
       "      <td>O'Higgins</td>\n",
       "      <td>CHL.1.2.2_1</td>\n",
       "      <td>(POLYGON ((-73.17006683349609 -49.224048614501...</td>\n",
       "      <td>{\"id\":{\"5775\":\"AREA # 11302\",\"5776\":\"AREA # 11...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        name          gid                                           geometry  \\\n",
       "0      Aisen  CHL.1.1.1_1  (POLYGON ((-73.49375152587885 -44.965415954589...   \n",
       "1     Cisnes  CHL.1.1.2_1  (POLYGON ((-73.31124877929682 -45.111248016357...   \n",
       "2  Guaitecas  CHL.1.1.3_1  (POLYGON ((-73.49153137207031 -43.711528778076...   \n",
       "3   Cochrane  CHL.1.2.1_1  (POLYGON ((-72.42082214355469 -47.743259429931...   \n",
       "4  O'Higgins  CHL.1.2.2_1  (POLYGON ((-73.17006683349609 -49.224048614501...   \n",
       "\n",
       "                                                expo  \n",
       "0  {\"id\":{\"5709\":\"AREA # 11201\",\"5710\":\"AREA # 11...  \n",
       "1  {\"id\":{\"5729\":\"AREA # 11202\",\"5730\":\"AREA # 11...  \n",
       "2  {\"id\":{\"5747\":\"AREA # 11203\",\"5748\":\"AREA # 11...  \n",
       "3  {\"id\":{\"5757\":\"AREA # 11301\",\"5758\":\"AREA # 11...  \n",
       "4  {\"id\":{\"5775\":\"AREA # 11302\",\"5776\":\"AREA # 11...  "
      ]
     },
     "execution_count": 146,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from shapely.geometry import MultiPolygon\n",
    "\n",
    "outmodel = chile_map_gpd[['name','GID_3','geometry']]\n",
    "outmodel['expo'] = \"\"\n",
    "#outmodel.head()\n",
    "outmodel = outmodel.set_index('name')\n",
    "\n",
    "for gcname,gc in outmodel.iterrows():\n",
    "    if (gcname in sara_names):\n",
    "            testgrp = updategroup(datagrp.get_group(gcname))\n",
    "            #add initial damage state - by default Null Damage (D0)\n",
    "            testgrp['Damage'] = 'D0'\n",
    "            outmodel.at[gcname,'expo'] = testgrp.to_json()\n",
    "    #convert polygons into multipolygons to avoid errors in file output\n",
    "    if (gc.geometry.geom_type == 'Polygon'):\n",
    "        outmodel.at[gcname,'geometry'] = MultiPolygon([gc.geometry])\n",
    "\n",
    "outmodel.columns = ['gid', 'geometry', 'expo']\n",
    "outmodel = outmodel.reset_index()\n",
    "outmodel.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Drop geocells with empty exposure "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 147,
   "metadata": {},
   "outputs": [],
   "source": [
    "blacklist  = outmodel[outmodel.expo ==''].index\n",
    "outmodel1 = outmodel.drop(blacklist)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 155,
   "metadata": {},
   "outputs": [],
   "source": [
    "#pd.read_json(outmodel1.iloc[2].expo)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Save to GEOPACKAGE "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 156,
   "metadata": {},
   "outputs": [],
   "source": [
    "#ioschema = gpd.io.file.infer_schema(outmodel)\n",
    "#ioschema['geometry']='MultiPolygon'\n",
    "#ioschema\n",
    "outpath = \"/Users/pittore/Documents/workspace/RIESGOS/GFZRIESGOS/updated_schemas/assetmaster\"\n",
    "outfile =os.path.join(outpath,\"SARA_v1.0_data.gpkg\")\n",
    "outmodel1.to_file(outfile, driver=\"GPKG\")#,schema = ioschema)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Display the boundary of the model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "map_osa = folium.Map(location=[-32,-71],tiles='Stamen Terrain',zoom_start=10)\n",
    "map_osa.choropleth(geo_data=outmodel1,fill_opacity=0.3)\n",
    "#map_osa.choropleth(geo_data=chile_map,fill_opacity=0.3)\n",
    "#folium.Marker([-32, -71], popup='test').add_to(map_osa)\n",
    "map_osa"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#map_osa = folium.Map(location=[df.iloc[event_index].latitude,df.iloc[event_index].longitude])\n",
    "#map_osm.create_map(path='osm.html')\n",
    "#map_osa.save('testfolium.html')\n",
    "#map_osa"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## UPDATING THE FRAGILITY MODEL"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 171,
   "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>taxonomy</th>\n",
       "      <th>D1_mean</th>\n",
       "      <th>D1_stddev</th>\n",
       "      <th>D2_mean</th>\n",
       "      <th>D2_stddev</th>\n",
       "      <th>D3_mean</th>\n",
       "      <th>D3_stddev</th>\n",
       "      <th>D4_mean</th>\n",
       "      <th>D4_stddev</th>\n",
       "      <th>id</th>\n",
       "      <th>imt</th>\n",
       "      <th>imu</th>\n",
       "      <th>im_min</th>\n",
       "      <th>im_max</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <td>0</td>\n",
       "      <td>MUR-H1</td>\n",
       "      <td>-1.418</td>\n",
       "      <td>0.310</td>\n",
       "      <td>-0.709</td>\n",
       "      <td>0.328</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>1</td>\n",
       "      <td>ER-ETR-H1-2</td>\n",
       "      <td>-1.418</td>\n",
       "      <td>0.310</td>\n",
       "      <td>-0.709</td>\n",
       "      <td>0.328</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>2</td>\n",
       "      <td>MUR-H1-3</td>\n",
       "      <td>-1.418</td>\n",
       "      <td>0.310</td>\n",
       "      <td>-0.709</td>\n",
       "      <td>0.328</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>3</td>\n",
       "      <td>MUR-ADO-H1-2</td>\n",
       "      <td>-1.782</td>\n",
       "      <td>0.371</td>\n",
       "      <td>-0.894</td>\n",
       "      <td>0.326</td>\n",
       "      <td>-0.644</td>\n",
       "      <td>0.328</td>\n",
       "      <td>-0.335</td>\n",
       "      <td>0.346</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>4</td>\n",
       "      <td>MUR-STRUB-H1-2</td>\n",
       "      <td>-1.502</td>\n",
       "      <td>0.295</td>\n",
       "      <td>-0.718</td>\n",
       "      <td>0.310</td>\n",
       "      <td>-0.492</td>\n",
       "      <td>0.311</td>\n",
       "      <td>-0.215</td>\n",
       "      <td>0.322</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5</td>\n",
       "      <td>MR-DUC-H1-3</td>\n",
       "      <td>-0.532</td>\n",
       "      <td>0.407</td>\n",
       "      <td>0.251</td>\n",
       "      <td>0.394</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.377</td>\n",
       "      <td>0.578</td>\n",
       "      <td>0.382</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>6</td>\n",
       "      <td>MCF-DUC-H1-3</td>\n",
       "      <td>-0.532</td>\n",
       "      <td>0.407</td>\n",
       "      <td>0.251</td>\n",
       "      <td>0.394</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.377</td>\n",
       "      <td>0.578</td>\n",
       "      <td>0.382</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>7</td>\n",
       "      <td>MCF-DNO-H1-3</td>\n",
       "      <td>-0.804</td>\n",
       "      <td>0.307</td>\n",
       "      <td>-0.046</td>\n",
       "      <td>0.320</td>\n",
       "      <td>0.127</td>\n",
       "      <td>0.314</td>\n",
       "      <td>0.386</td>\n",
       "      <td>0.367</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>8</td>\n",
       "      <td>MR-DNO-H1-3</td>\n",
       "      <td>-0.804</td>\n",
       "      <td>0.307</td>\n",
       "      <td>-0.046</td>\n",
       "      <td>0.320</td>\n",
       "      <td>0.127</td>\n",
       "      <td>0.314</td>\n",
       "      <td>0.386</td>\n",
       "      <td>0.367</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>9</td>\n",
       "      <td>CR-LFM-DUC-H1-3</td>\n",
       "      <td>0.088</td>\n",
       "      <td>0.348</td>\n",
       "      <td>0.919</td>\n",
       "      <td>0.398</td>\n",
       "      <td>1.179</td>\n",
       "      <td>0.381</td>\n",
       "      <td>1.616</td>\n",
       "      <td>0.507</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_03</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>10</td>\n",
       "      <td>CR-LFLS-DUC-H1-3</td>\n",
       "      <td>0.088</td>\n",
       "      <td>0.348</td>\n",
       "      <td>0.919</td>\n",
       "      <td>0.398</td>\n",
       "      <td>1.179</td>\n",
       "      <td>0.381</td>\n",
       "      <td>1.616</td>\n",
       "      <td>0.507</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_03</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>11</td>\n",
       "      <td>CR-LFM-DUC-H4-7</td>\n",
       "      <td>-1.257</td>\n",
       "      <td>0.327</td>\n",
       "      <td>-0.419</td>\n",
       "      <td>0.350</td>\n",
       "      <td>-0.058</td>\n",
       "      <td>0.404</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_01</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>12</td>\n",
       "      <td>CR-LFLS-DUC-H4-7</td>\n",
       "      <td>-1.257</td>\n",
       "      <td>0.327</td>\n",
       "      <td>-0.419</td>\n",
       "      <td>0.350</td>\n",
       "      <td>-0.058</td>\n",
       "      <td>0.404</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_01</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>13</td>\n",
       "      <td>CR-LDUAL-DUC-H4-7</td>\n",
       "      <td>-1.257</td>\n",
       "      <td>0.327</td>\n",
       "      <td>-0.419</td>\n",
       "      <td>0.350</td>\n",
       "      <td>-0.058</td>\n",
       "      <td>0.404</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_01</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>14</td>\n",
       "      <td>CR-LFLS-DNO-H1-3</td>\n",
       "      <td>-0.785</td>\n",
       "      <td>0.333</td>\n",
       "      <td>-0.006</td>\n",
       "      <td>0.319</td>\n",
       "      <td>-0.256</td>\n",
       "      <td>0.372</td>\n",
       "      <td>-0.643</td>\n",
       "      <td>0.437</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_03</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>15</td>\n",
       "      <td>CR-LFM-DNO-H1-3</td>\n",
       "      <td>-0.785</td>\n",
       "      <td>0.333</td>\n",
       "      <td>-0.006</td>\n",
       "      <td>0.319</td>\n",
       "      <td>-0.256</td>\n",
       "      <td>0.372</td>\n",
       "      <td>-0.643</td>\n",
       "      <td>0.437</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_03</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>16</td>\n",
       "      <td>CR-LFM-DNO-H4-7</td>\n",
       "      <td>-1.841</td>\n",
       "      <td>0.292</td>\n",
       "      <td>-1.068</td>\n",
       "      <td>0.403</td>\n",
       "      <td>-0.717</td>\n",
       "      <td>0.474</td>\n",
       "      <td>-0.195</td>\n",
       "      <td>0.591</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_01</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>17</td>\n",
       "      <td>CR-LFM-DNO-SOS-H1-3</td>\n",
       "      <td>0.820</td>\n",
       "      <td>0.317</td>\n",
       "      <td>0.054</td>\n",
       "      <td>0.334</td>\n",
       "      <td>0.205</td>\n",
       "      <td>0.383</td>\n",
       "      <td>0.582</td>\n",
       "      <td>0.469</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>18</td>\n",
       "      <td>CR-LFM-DNO-SOS-H4-7</td>\n",
       "      <td>-1.715</td>\n",
       "      <td>0.342</td>\n",
       "      <td>-1.187</td>\n",
       "      <td>0.423</td>\n",
       "      <td>-0.977</td>\n",
       "      <td>0.480</td>\n",
       "      <td>-0.638</td>\n",
       "      <td>0.566</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_01</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>19</td>\n",
       "      <td>CR-LFINF-DUC-H1-3</td>\n",
       "      <td>-0.401</td>\n",
       "      <td>0.410</td>\n",
       "      <td>0.482</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.694</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>20</td>\n",
       "      <td>CR-LWAL-DUC-H1-3</td>\n",
       "      <td>-0.401</td>\n",
       "      <td>0.410</td>\n",
       "      <td>0.482</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.694</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>21</td>\n",
       "      <td>CR-PC-LWAL-H1-3</td>\n",
       "      <td>-0.401</td>\n",
       "      <td>0.410</td>\n",
       "      <td>0.482</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.694</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>22</td>\n",
       "      <td>CR-LFINF-DUC-H4-7</td>\n",
       "      <td>-2.023</td>\n",
       "      <td>0.465</td>\n",
       "      <td>-0.827</td>\n",
       "      <td>0.589</td>\n",
       "      <td>-0.307</td>\n",
       "      <td>0.617</td>\n",
       "      <td>0.288</td>\n",
       "      <td>0.636</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_01</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>23</td>\n",
       "      <td>CR-LWAL-DUC-H4-7</td>\n",
       "      <td>-2.023</td>\n",
       "      <td>0.465</td>\n",
       "      <td>-0.827</td>\n",
       "      <td>0.589</td>\n",
       "      <td>-0.307</td>\n",
       "      <td>0.617</td>\n",
       "      <td>0.288</td>\n",
       "      <td>0.636</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_01</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>24</td>\n",
       "      <td>CR-LFINF-DNO-H1-3</td>\n",
       "      <td>-1.126</td>\n",
       "      <td>0.363</td>\n",
       "      <td>-0.330</td>\n",
       "      <td>0.350</td>\n",
       "      <td>-0.119</td>\n",
       "      <td>0.331</td>\n",
       "      <td>0.162</td>\n",
       "      <td>0.400</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>25</td>\n",
       "      <td>CR-LWAL-DNO-H1-3</td>\n",
       "      <td>-1.126</td>\n",
       "      <td>0.363</td>\n",
       "      <td>-0.330</td>\n",
       "      <td>0.350</td>\n",
       "      <td>-0.119</td>\n",
       "      <td>0.331</td>\n",
       "      <td>0.162</td>\n",
       "      <td>0.400</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>pga</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>26</td>\n",
       "      <td>CR-LWAL-DNO-H4-7</td>\n",
       "      <td>-2.401</td>\n",
       "      <td>0.362</td>\n",
       "      <td>-1.400</td>\n",
       "      <td>0.430</td>\n",
       "      <td>-0.981</td>\n",
       "      <td>0.444</td>\n",
       "      <td>-0.402</td>\n",
       "      <td>0.548</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_01</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>27</td>\n",
       "      <td>CR-LDUAL-DUC-H8-19</td>\n",
       "      <td>-0.816</td>\n",
       "      <td>0.294</td>\n",
       "      <td>-0.260</td>\n",
       "      <td>0.361</td>\n",
       "      <td>-0.012</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.371</td>\n",
       "      <td>0.490</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_01</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>28</td>\n",
       "      <td>CR-LWAL-DUC-H8-19</td>\n",
       "      <td>-0.816</td>\n",
       "      <td>0.294</td>\n",
       "      <td>-0.260</td>\n",
       "      <td>0.361</td>\n",
       "      <td>-0.012</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.371</td>\n",
       "      <td>0.490</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_01</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>29</td>\n",
       "      <td>W-WLI-H1-3</td>\n",
       "      <td>0.107</td>\n",
       "      <td>0.403</td>\n",
       "      <td>0.853</td>\n",
       "      <td>0.450</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_03</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>30</td>\n",
       "      <td>W-WHE-H1-3</td>\n",
       "      <td>0.107</td>\n",
       "      <td>0.403</td>\n",
       "      <td>0.853</td>\n",
       "      <td>0.450</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_03</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>31</td>\n",
       "      <td>W-WBB-H1</td>\n",
       "      <td>0.107</td>\n",
       "      <td>0.403</td>\n",
       "      <td>0.853</td>\n",
       "      <td>0.450</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_03</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>32</td>\n",
       "      <td>W-WS-H1-2</td>\n",
       "      <td>0.107</td>\n",
       "      <td>0.403</td>\n",
       "      <td>0.853</td>\n",
       "      <td>0.450</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_03</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>33</td>\n",
       "      <td>UNK</td>\n",
       "      <td>-0.877</td>\n",
       "      <td>0.360</td>\n",
       "      <td>-0.143</td>\n",
       "      <td>0.401</td>\n",
       "      <td>0.151</td>\n",
       "      <td>0.362</td>\n",
       "      <td>0.489</td>\n",
       "      <td>0.378</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>SA_03</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>34</td>\n",
       "      <td>S-LFM-H4-7</td>\n",
       "      <td>0.227</td>\n",
       "      <td>0.121</td>\n",
       "      <td>0.310</td>\n",
       "      <td>0.147</td>\n",
       "      <td>0.347</td>\n",
       "      <td>0.144</td>\n",
       "      <td>0.390</td>\n",
       "      <td>0.162</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>PGA</td>\n",
       "      <td>g</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "               taxonomy  D1_mean  D1_stddev  D2_mean  D2_stddev  D3_mean  \\\n",
       "0                MUR-H1   -1.418      0.310   -0.709      0.328   -0.496   \n",
       "1           ER-ETR-H1-2   -1.418      0.310   -0.709      0.328   -0.496   \n",
       "2              MUR-H1-3   -1.418      0.310   -0.709      0.328   -0.496   \n",
       "3          MUR-ADO-H1-2   -1.782      0.371   -0.894      0.326   -0.644   \n",
       "4        MUR-STRUB-H1-2   -1.502      0.295   -0.718      0.310   -0.492   \n",
       "5           MR-DUC-H1-3   -0.532      0.407    0.251      0.394    0.364   \n",
       "6          MCF-DUC-H1-3   -0.532      0.407    0.251      0.394    0.364   \n",
       "7          MCF-DNO-H1-3   -0.804      0.307   -0.046      0.320    0.127   \n",
       "8           MR-DNO-H1-3   -0.804      0.307   -0.046      0.320    0.127   \n",
       "9       CR-LFM-DUC-H1-3    0.088      0.348    0.919      0.398    1.179   \n",
       "10     CR-LFLS-DUC-H1-3    0.088      0.348    0.919      0.398    1.179   \n",
       "11      CR-LFM-DUC-H4-7   -1.257      0.327   -0.419      0.350   -0.058   \n",
       "12     CR-LFLS-DUC-H4-7   -1.257      0.327   -0.419      0.350   -0.058   \n",
       "13    CR-LDUAL-DUC-H4-7   -1.257      0.327   -0.419      0.350   -0.058   \n",
       "14     CR-LFLS-DNO-H1-3   -0.785      0.333   -0.006      0.319   -0.256   \n",
       "15      CR-LFM-DNO-H1-3   -0.785      0.333   -0.006      0.319   -0.256   \n",
       "16      CR-LFM-DNO-H4-7   -1.841      0.292   -1.068      0.403   -0.717   \n",
       "17  CR-LFM-DNO-SOS-H1-3    0.820      0.317    0.054      0.334    0.205   \n",
       "18  CR-LFM-DNO-SOS-H4-7   -1.715      0.342   -1.187      0.423   -0.977   \n",
       "19    CR-LFINF-DUC-H1-3   -0.401      0.410    0.482      0.391    0.694   \n",
       "20     CR-LWAL-DUC-H1-3   -0.401      0.410    0.482      0.391    0.694   \n",
       "21      CR-PC-LWAL-H1-3   -0.401      0.410    0.482      0.391    0.694   \n",
       "22    CR-LFINF-DUC-H4-7   -2.023      0.465   -0.827      0.589   -0.307   \n",
       "23     CR-LWAL-DUC-H4-7   -2.023      0.465   -0.827      0.589   -0.307   \n",
       "24    CR-LFINF-DNO-H1-3   -1.126      0.363   -0.330      0.350   -0.119   \n",
       "25     CR-LWAL-DNO-H1-3   -1.126      0.363   -0.330      0.350   -0.119   \n",
       "26     CR-LWAL-DNO-H4-7   -2.401      0.362   -1.400      0.430   -0.981   \n",
       "27   CR-LDUAL-DUC-H8-19   -0.816      0.294   -0.260      0.361   -0.012   \n",
       "28    CR-LWAL-DUC-H8-19   -0.816      0.294   -0.260      0.361   -0.012   \n",
       "29           W-WLI-H1-3    0.107      0.403    0.853      0.450    1.144   \n",
       "30           W-WHE-H1-3    0.107      0.403    0.853      0.450    1.144   \n",
       "31             W-WBB-H1    0.107      0.403    0.853      0.450    1.144   \n",
       "32            W-WS-H1-2    0.107      0.403    0.853      0.450    1.144   \n",
       "33                  UNK   -0.877      0.360   -0.143      0.401    0.151   \n",
       "34           S-LFM-H4-7    0.227      0.121    0.310      0.147    0.347   \n",
       "\n",
       "    D3_stddev  D4_mean  D4_stddev         id    imt imu  im_min  im_max  \n",
       "0       0.322   -0.231      0.317  unique_id    pga   g     0.0     1.0  \n",
       "1       0.322   -0.231      0.317  unique_id    pga   g     0.0     1.0  \n",
       "2       0.322   -0.231      0.317  unique_id    pga   g     0.0     1.0  \n",
       "3       0.328   -0.335      0.346  unique_id    pga   g     0.0     1.0  \n",
       "4       0.311   -0.215      0.322  unique_id    pga   g     0.0     1.0  \n",
       "5       0.377    0.578      0.382  unique_id    pga   g     0.0     1.0  \n",
       "6       0.377    0.578      0.382  unique_id    pga   g     0.0     1.0  \n",
       "7       0.314    0.386      0.367  unique_id    pga   g     0.0     1.0  \n",
       "8       0.314    0.386      0.367  unique_id    pga   g     0.0     1.0  \n",
       "9       0.381    1.616      0.507  unique_id  SA_03   g     0.0     3.0  \n",
       "10      0.381    1.616      0.507  unique_id  SA_03   g     0.0     3.0  \n",
       "11      0.404    0.436      0.525  unique_id  SA_01   g     0.0     1.2  \n",
       "12      0.404    0.436      0.525  unique_id  SA_01   g     0.0     1.2  \n",
       "13      0.404    0.436      0.525  unique_id  SA_01   g     0.0     1.2  \n",
       "14      0.372   -0.643      0.437  unique_id  SA_03   g     0.0     3.0  \n",
       "15      0.372   -0.643      0.437  unique_id  SA_03   g     0.0     3.0  \n",
       "16      0.474   -0.195      0.591  unique_id  SA_01   g     0.0     1.2  \n",
       "17      0.383    0.582      0.469  unique_id    pga   g     0.0     1.0  \n",
       "18      0.480   -0.638      0.566  unique_id  SA_01   g     0.0     1.2  \n",
       "19      0.364    0.973      0.406  unique_id    pga   g     0.0     1.0  \n",
       "20      0.364    0.973      0.406  unique_id    pga   g     0.0     1.0  \n",
       "21      0.364    0.973      0.406  unique_id    pga   g     0.0     1.0  \n",
       "22      0.617    0.288      0.636  unique_id  SA_01   g     0.0     1.2  \n",
       "23      0.617    0.288      0.636  unique_id  SA_01   g     0.0     1.2  \n",
       "24      0.331    0.162      0.400  unique_id    pga   g     0.0     1.0  \n",
       "25      0.331    0.162      0.400  unique_id    pga   g     0.0     1.0  \n",
       "26      0.444   -0.402      0.548  unique_id  SA_01   g     0.0     1.2  \n",
       "27      0.391    0.371      0.490  unique_id  SA_01   g     0.0     1.2  \n",
       "28      0.391    0.371      0.490  unique_id  SA_01   g     0.0     1.2  \n",
       "29      0.449    1.594      0.573  unique_id  SA_03   g     0.0     3.0  \n",
       "30      0.449    1.594      0.573  unique_id  SA_03   g     0.0     3.0  \n",
       "31      0.449    1.594      0.573  unique_id  SA_03   g     0.0     3.0  \n",
       "32      0.449    1.594      0.573  unique_id  SA_03   g     0.0     3.0  \n",
       "33      0.362    0.489      0.378  unique_id  SA_03   g     0.0     3.0  \n",
       "34      0.144    0.390      0.162  unique_id    PGA   g     0.0     1.0  "
      ]
     },
     "execution_count": 171,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import json\n",
    "#tstfile = \"/Users/pittore/PowerFolders/riesgos/CODES/WEB_services_and_software/assetmaster/schemas/SARA_v1.0/SARA_v1.0_data.gpkg\"\n",
    "#tst_df = gpd.read_file(outfile,driver=\"GPKG\")\n",
    "#pd.read_json(tst_df.iloc[0].expo)\n",
    "\n",
    "tstfile0=\"/Users/pittore/Documents/workspace/RIESGOS/GFZRIESGOS/updated_schemas/modelprop/SARA_v1.0.json\"\n",
    "jsondata = json.load(open(tstfile0,'r'))\n",
    "df1 = pd.DataFrame(jsondata[\"data\"])\n",
    "df1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 173,
   "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>taxonomy</th>\n",
       "      <th>D1_mean</th>\n",
       "      <th>D1_stddev</th>\n",
       "      <th>D2_mean</th>\n",
       "      <th>D2_stddev</th>\n",
       "      <th>D3_mean</th>\n",
       "      <th>D3_stddev</th>\n",
       "      <th>D4_mean</th>\n",
       "      <th>D4_stddev</th>\n",
       "      <th>id</th>\n",
       "      <th>...</th>\n",
       "      <th>D1_3_mean</th>\n",
       "      <th>D1_3_stddev</th>\n",
       "      <th>D1_4_mean</th>\n",
       "      <th>D1_4_stddev</th>\n",
       "      <th>D2_3_mean</th>\n",
       "      <th>D2_3_stddev</th>\n",
       "      <th>D2_4_mean</th>\n",
       "      <th>D2_4_stddev</th>\n",
       "      <th>D3_4_mean</th>\n",
       "      <th>D3_4_stddev</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <td>0</td>\n",
       "      <td>MUR-H1</td>\n",
       "      <td>-1.418</td>\n",
       "      <td>0.310</td>\n",
       "      <td>-0.709</td>\n",
       "      <td>0.328</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>1</td>\n",
       "      <td>ER_ETR_H1</td>\n",
       "      <td>-1.418</td>\n",
       "      <td>0.310</td>\n",
       "      <td>-0.709</td>\n",
       "      <td>0.328</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>2</td>\n",
       "      <td>ER_ETR_H1_2</td>\n",
       "      <td>-1.418</td>\n",
       "      <td>0.310</td>\n",
       "      <td>-0.709</td>\n",
       "      <td>0.328</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>3</td>\n",
       "      <td>MUR-H1_3</td>\n",
       "      <td>-1.418</td>\n",
       "      <td>0.310</td>\n",
       "      <td>-0.709</td>\n",
       "      <td>0.328</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>-0.496</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "      <td>-0.231</td>\n",
       "      <td>0.317</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>4</td>\n",
       "      <td>MUR_ADO_H1_2</td>\n",
       "      <td>-1.782</td>\n",
       "      <td>0.371</td>\n",
       "      <td>-0.894</td>\n",
       "      <td>0.326</td>\n",
       "      <td>-0.644</td>\n",
       "      <td>0.328</td>\n",
       "      <td>-0.335</td>\n",
       "      <td>0.346</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.644</td>\n",
       "      <td>0.328</td>\n",
       "      <td>-0.335</td>\n",
       "      <td>0.346</td>\n",
       "      <td>-0.644</td>\n",
       "      <td>0.328</td>\n",
       "      <td>-0.335</td>\n",
       "      <td>0.346</td>\n",
       "      <td>-0.335</td>\n",
       "      <td>0.346</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>5</td>\n",
       "      <td>MUR-STRUB-H1-2</td>\n",
       "      <td>-1.502</td>\n",
       "      <td>0.295</td>\n",
       "      <td>-0.718</td>\n",
       "      <td>0.310</td>\n",
       "      <td>-0.492</td>\n",
       "      <td>0.311</td>\n",
       "      <td>-0.215</td>\n",
       "      <td>0.322</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.492</td>\n",
       "      <td>0.311</td>\n",
       "      <td>-0.215</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.492</td>\n",
       "      <td>0.311</td>\n",
       "      <td>-0.215</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.215</td>\n",
       "      <td>0.322</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>6</td>\n",
       "      <td>MUR-STDRE-H1-2</td>\n",
       "      <td>-1.502</td>\n",
       "      <td>0.295</td>\n",
       "      <td>-0.718</td>\n",
       "      <td>0.310</td>\n",
       "      <td>-0.492</td>\n",
       "      <td>0.311</td>\n",
       "      <td>-0.215</td>\n",
       "      <td>0.322</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.492</td>\n",
       "      <td>0.311</td>\n",
       "      <td>-0.215</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.492</td>\n",
       "      <td>0.311</td>\n",
       "      <td>-0.215</td>\n",
       "      <td>0.322</td>\n",
       "      <td>-0.215</td>\n",
       "      <td>0.322</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>7</td>\n",
       "      <td>MR_DUC_H1_3</td>\n",
       "      <td>-0.532</td>\n",
       "      <td>0.407</td>\n",
       "      <td>0.251</td>\n",
       "      <td>0.394</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.377</td>\n",
       "      <td>0.578</td>\n",
       "      <td>0.382</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.377</td>\n",
       "      <td>0.578</td>\n",
       "      <td>0.382</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.377</td>\n",
       "      <td>0.578</td>\n",
       "      <td>0.382</td>\n",
       "      <td>0.578</td>\n",
       "      <td>0.382</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>8</td>\n",
       "      <td>MCF_DUC_H1_3</td>\n",
       "      <td>-0.532</td>\n",
       "      <td>0.407</td>\n",
       "      <td>0.251</td>\n",
       "      <td>0.394</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.377</td>\n",
       "      <td>0.578</td>\n",
       "      <td>0.382</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.377</td>\n",
       "      <td>0.578</td>\n",
       "      <td>0.382</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.377</td>\n",
       "      <td>0.578</td>\n",
       "      <td>0.382</td>\n",
       "      <td>0.578</td>\n",
       "      <td>0.382</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>9</td>\n",
       "      <td>MCF-DNO-H1-3</td>\n",
       "      <td>-0.804</td>\n",
       "      <td>0.307</td>\n",
       "      <td>-0.046</td>\n",
       "      <td>0.320</td>\n",
       "      <td>0.127</td>\n",
       "      <td>0.314</td>\n",
       "      <td>0.386</td>\n",
       "      <td>0.367</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>0.127</td>\n",
       "      <td>0.314</td>\n",
       "      <td>0.386</td>\n",
       "      <td>0.367</td>\n",
       "      <td>0.127</td>\n",
       "      <td>0.314</td>\n",
       "      <td>0.386</td>\n",
       "      <td>0.367</td>\n",
       "      <td>0.386</td>\n",
       "      <td>0.367</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>10</td>\n",
       "      <td>MR-DNO-H1-3</td>\n",
       "      <td>-0.804</td>\n",
       "      <td>0.307</td>\n",
       "      <td>-0.046</td>\n",
       "      <td>0.320</td>\n",
       "      <td>0.127</td>\n",
       "      <td>0.314</td>\n",
       "      <td>0.386</td>\n",
       "      <td>0.367</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>0.127</td>\n",
       "      <td>0.314</td>\n",
       "      <td>0.386</td>\n",
       "      <td>0.367</td>\n",
       "      <td>0.127</td>\n",
       "      <td>0.314</td>\n",
       "      <td>0.386</td>\n",
       "      <td>0.367</td>\n",
       "      <td>0.386</td>\n",
       "      <td>0.367</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>11</td>\n",
       "      <td>CR-LFM-DUC-H1-3</td>\n",
       "      <td>0.088</td>\n",
       "      <td>0.348</td>\n",
       "      <td>0.919</td>\n",
       "      <td>0.398</td>\n",
       "      <td>1.179</td>\n",
       "      <td>0.381</td>\n",
       "      <td>1.616</td>\n",
       "      <td>0.507</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>1.179</td>\n",
       "      <td>0.381</td>\n",
       "      <td>1.616</td>\n",
       "      <td>0.507</td>\n",
       "      <td>1.179</td>\n",
       "      <td>0.381</td>\n",
       "      <td>1.616</td>\n",
       "      <td>0.507</td>\n",
       "      <td>1.616</td>\n",
       "      <td>0.507</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>12</td>\n",
       "      <td>CR-LFLS-DUC-H1-3</td>\n",
       "      <td>0.088</td>\n",
       "      <td>0.348</td>\n",
       "      <td>0.919</td>\n",
       "      <td>0.398</td>\n",
       "      <td>1.179</td>\n",
       "      <td>0.381</td>\n",
       "      <td>1.616</td>\n",
       "      <td>0.507</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>1.179</td>\n",
       "      <td>0.381</td>\n",
       "      <td>1.616</td>\n",
       "      <td>0.507</td>\n",
       "      <td>1.179</td>\n",
       "      <td>0.381</td>\n",
       "      <td>1.616</td>\n",
       "      <td>0.507</td>\n",
       "      <td>1.616</td>\n",
       "      <td>0.507</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>13</td>\n",
       "      <td>CR-LFM-DUC-H4-7</td>\n",
       "      <td>-1.257</td>\n",
       "      <td>0.327</td>\n",
       "      <td>-0.419</td>\n",
       "      <td>0.350</td>\n",
       "      <td>-0.058</td>\n",
       "      <td>0.404</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.058</td>\n",
       "      <td>0.404</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "      <td>-0.058</td>\n",
       "      <td>0.404</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>14</td>\n",
       "      <td>CR-LFLS-DUC-H4-7</td>\n",
       "      <td>-1.257</td>\n",
       "      <td>0.327</td>\n",
       "      <td>-0.419</td>\n",
       "      <td>0.350</td>\n",
       "      <td>-0.058</td>\n",
       "      <td>0.404</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.058</td>\n",
       "      <td>0.404</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "      <td>-0.058</td>\n",
       "      <td>0.404</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>15</td>\n",
       "      <td>CR-LDUAL-DUC-H4-7</td>\n",
       "      <td>-1.257</td>\n",
       "      <td>0.327</td>\n",
       "      <td>-0.419</td>\n",
       "      <td>0.350</td>\n",
       "      <td>-0.058</td>\n",
       "      <td>0.404</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.058</td>\n",
       "      <td>0.404</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "      <td>-0.058</td>\n",
       "      <td>0.404</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "      <td>0.436</td>\n",
       "      <td>0.525</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>16</td>\n",
       "      <td>CR-LFLS-DNO-H1-3</td>\n",
       "      <td>-0.785</td>\n",
       "      <td>0.333</td>\n",
       "      <td>-0.006</td>\n",
       "      <td>0.319</td>\n",
       "      <td>-0.256</td>\n",
       "      <td>0.372</td>\n",
       "      <td>-0.643</td>\n",
       "      <td>0.437</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.256</td>\n",
       "      <td>0.372</td>\n",
       "      <td>-0.643</td>\n",
       "      <td>0.437</td>\n",
       "      <td>-0.256</td>\n",
       "      <td>0.372</td>\n",
       "      <td>-0.643</td>\n",
       "      <td>0.437</td>\n",
       "      <td>-0.643</td>\n",
       "      <td>0.437</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>17</td>\n",
       "      <td>CR-LFM-DNO-H1-3</td>\n",
       "      <td>-0.785</td>\n",
       "      <td>0.333</td>\n",
       "      <td>-0.006</td>\n",
       "      <td>0.319</td>\n",
       "      <td>-0.256</td>\n",
       "      <td>0.372</td>\n",
       "      <td>-0.643</td>\n",
       "      <td>0.437</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.256</td>\n",
       "      <td>0.372</td>\n",
       "      <td>-0.643</td>\n",
       "      <td>0.437</td>\n",
       "      <td>-0.256</td>\n",
       "      <td>0.372</td>\n",
       "      <td>-0.643</td>\n",
       "      <td>0.437</td>\n",
       "      <td>-0.643</td>\n",
       "      <td>0.437</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>18</td>\n",
       "      <td>CR-LFM-DNO-H4-7</td>\n",
       "      <td>-1.841</td>\n",
       "      <td>0.292</td>\n",
       "      <td>-1.068</td>\n",
       "      <td>0.403</td>\n",
       "      <td>-0.717</td>\n",
       "      <td>0.474</td>\n",
       "      <td>-0.195</td>\n",
       "      <td>0.591</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.717</td>\n",
       "      <td>0.474</td>\n",
       "      <td>-0.195</td>\n",
       "      <td>0.591</td>\n",
       "      <td>-0.717</td>\n",
       "      <td>0.474</td>\n",
       "      <td>-0.195</td>\n",
       "      <td>0.591</td>\n",
       "      <td>-0.195</td>\n",
       "      <td>0.591</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>19</td>\n",
       "      <td>CR-LFM-DNO-SOS-H1-3</td>\n",
       "      <td>0.820</td>\n",
       "      <td>0.317</td>\n",
       "      <td>0.054</td>\n",
       "      <td>0.334</td>\n",
       "      <td>0.205</td>\n",
       "      <td>0.383</td>\n",
       "      <td>0.582</td>\n",
       "      <td>0.469</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>0.205</td>\n",
       "      <td>0.383</td>\n",
       "      <td>0.582</td>\n",
       "      <td>0.469</td>\n",
       "      <td>0.205</td>\n",
       "      <td>0.383</td>\n",
       "      <td>0.582</td>\n",
       "      <td>0.469</td>\n",
       "      <td>0.582</td>\n",
       "      <td>0.469</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>20</td>\n",
       "      <td>CR-LFM-DNO-SOS-H4-7</td>\n",
       "      <td>-1.715</td>\n",
       "      <td>0.342</td>\n",
       "      <td>-1.187</td>\n",
       "      <td>0.423</td>\n",
       "      <td>-0.977</td>\n",
       "      <td>0.480</td>\n",
       "      <td>-0.638</td>\n",
       "      <td>0.566</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.977</td>\n",
       "      <td>0.480</td>\n",
       "      <td>-0.638</td>\n",
       "      <td>0.566</td>\n",
       "      <td>-0.977</td>\n",
       "      <td>0.480</td>\n",
       "      <td>-0.638</td>\n",
       "      <td>0.566</td>\n",
       "      <td>-0.638</td>\n",
       "      <td>0.566</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>21</td>\n",
       "      <td>CR-LFINF-DUC-H1-3</td>\n",
       "      <td>-0.401</td>\n",
       "      <td>0.410</td>\n",
       "      <td>0.482</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.694</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>0.694</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "      <td>0.694</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>22</td>\n",
       "      <td>CR-LWAL-DUC-H1-3</td>\n",
       "      <td>-0.401</td>\n",
       "      <td>0.410</td>\n",
       "      <td>0.482</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.694</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>0.694</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "      <td>0.694</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>23</td>\n",
       "      <td>CR-PC-LWAL-H1-3</td>\n",
       "      <td>-0.401</td>\n",
       "      <td>0.410</td>\n",
       "      <td>0.482</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.694</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>0.694</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "      <td>0.694</td>\n",
       "      <td>0.364</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "      <td>0.973</td>\n",
       "      <td>0.406</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>24</td>\n",
       "      <td>CR-LFINF-DUC-H4-7</td>\n",
       "      <td>-2.023</td>\n",
       "      <td>0.465</td>\n",
       "      <td>-0.827</td>\n",
       "      <td>0.589</td>\n",
       "      <td>-0.307</td>\n",
       "      <td>0.617</td>\n",
       "      <td>0.288</td>\n",
       "      <td>0.636</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.307</td>\n",
       "      <td>0.617</td>\n",
       "      <td>0.288</td>\n",
       "      <td>0.636</td>\n",
       "      <td>-0.307</td>\n",
       "      <td>0.617</td>\n",
       "      <td>0.288</td>\n",
       "      <td>0.636</td>\n",
       "      <td>0.288</td>\n",
       "      <td>0.636</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>25</td>\n",
       "      <td>CR-LWAL-DUC-H4-7</td>\n",
       "      <td>-2.023</td>\n",
       "      <td>0.465</td>\n",
       "      <td>-0.827</td>\n",
       "      <td>0.589</td>\n",
       "      <td>-0.307</td>\n",
       "      <td>0.617</td>\n",
       "      <td>0.288</td>\n",
       "      <td>0.636</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.307</td>\n",
       "      <td>0.617</td>\n",
       "      <td>0.288</td>\n",
       "      <td>0.636</td>\n",
       "      <td>-0.307</td>\n",
       "      <td>0.617</td>\n",
       "      <td>0.288</td>\n",
       "      <td>0.636</td>\n",
       "      <td>0.288</td>\n",
       "      <td>0.636</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>26</td>\n",
       "      <td>CR-LFINF-DNO-H1-3</td>\n",
       "      <td>-1.126</td>\n",
       "      <td>0.363</td>\n",
       "      <td>-0.330</td>\n",
       "      <td>0.350</td>\n",
       "      <td>-0.119</td>\n",
       "      <td>0.331</td>\n",
       "      <td>0.162</td>\n",
       "      <td>0.400</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.119</td>\n",
       "      <td>0.331</td>\n",
       "      <td>0.162</td>\n",
       "      <td>0.400</td>\n",
       "      <td>-0.119</td>\n",
       "      <td>0.331</td>\n",
       "      <td>0.162</td>\n",
       "      <td>0.400</td>\n",
       "      <td>0.162</td>\n",
       "      <td>0.400</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>27</td>\n",
       "      <td>CR-LWAL-DNO-H1-3</td>\n",
       "      <td>-1.126</td>\n",
       "      <td>0.363</td>\n",
       "      <td>-0.330</td>\n",
       "      <td>0.350</td>\n",
       "      <td>-0.119</td>\n",
       "      <td>0.331</td>\n",
       "      <td>0.162</td>\n",
       "      <td>0.400</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.119</td>\n",
       "      <td>0.331</td>\n",
       "      <td>0.162</td>\n",
       "      <td>0.400</td>\n",
       "      <td>-0.119</td>\n",
       "      <td>0.331</td>\n",
       "      <td>0.162</td>\n",
       "      <td>0.400</td>\n",
       "      <td>0.162</td>\n",
       "      <td>0.400</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>28</td>\n",
       "      <td>CR-LWAL-DNO-H4-7</td>\n",
       "      <td>-2.401</td>\n",
       "      <td>0.362</td>\n",
       "      <td>-1.400</td>\n",
       "      <td>0.430</td>\n",
       "      <td>-0.981</td>\n",
       "      <td>0.444</td>\n",
       "      <td>-0.402</td>\n",
       "      <td>0.548</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.981</td>\n",
       "      <td>0.444</td>\n",
       "      <td>-0.402</td>\n",
       "      <td>0.548</td>\n",
       "      <td>-0.981</td>\n",
       "      <td>0.444</td>\n",
       "      <td>-0.402</td>\n",
       "      <td>0.548</td>\n",
       "      <td>-0.402</td>\n",
       "      <td>0.548</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>29</td>\n",
       "      <td>CR-LDUAL-DUC-H8-19</td>\n",
       "      <td>-0.816</td>\n",
       "      <td>0.294</td>\n",
       "      <td>-0.260</td>\n",
       "      <td>0.361</td>\n",
       "      <td>-0.012</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.371</td>\n",
       "      <td>0.490</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.012</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.371</td>\n",
       "      <td>0.490</td>\n",
       "      <td>-0.012</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.371</td>\n",
       "      <td>0.490</td>\n",
       "      <td>0.371</td>\n",
       "      <td>0.490</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>30</td>\n",
       "      <td>CR-LWAL-DUC-H8-19</td>\n",
       "      <td>-0.816</td>\n",
       "      <td>0.294</td>\n",
       "      <td>-0.260</td>\n",
       "      <td>0.361</td>\n",
       "      <td>-0.012</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.371</td>\n",
       "      <td>0.490</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.012</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.371</td>\n",
       "      <td>0.490</td>\n",
       "      <td>-0.012</td>\n",
       "      <td>0.391</td>\n",
       "      <td>0.371</td>\n",
       "      <td>0.490</td>\n",
       "      <td>0.371</td>\n",
       "      <td>0.490</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>31</td>\n",
       "      <td>W-WLI-H1</td>\n",
       "      <td>0.107</td>\n",
       "      <td>0.403</td>\n",
       "      <td>0.853</td>\n",
       "      <td>0.450</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>32</td>\n",
       "      <td>W-WLI-H1-3</td>\n",
       "      <td>0.107</td>\n",
       "      <td>0.403</td>\n",
       "      <td>0.853</td>\n",
       "      <td>0.450</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>33</td>\n",
       "      <td>W-WHE-H1-3</td>\n",
       "      <td>0.107</td>\n",
       "      <td>0.403</td>\n",
       "      <td>0.853</td>\n",
       "      <td>0.450</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>34</td>\n",
       "      <td>W-WBB-H1</td>\n",
       "      <td>0.107</td>\n",
       "      <td>0.403</td>\n",
       "      <td>0.853</td>\n",
       "      <td>0.450</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>35</td>\n",
       "      <td>W-WS-H1</td>\n",
       "      <td>0.107</td>\n",
       "      <td>0.403</td>\n",
       "      <td>0.853</td>\n",
       "      <td>0.450</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>36</td>\n",
       "      <td>W-WS-H1-2</td>\n",
       "      <td>0.107</td>\n",
       "      <td>0.403</td>\n",
       "      <td>0.853</td>\n",
       "      <td>0.450</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>1.144</td>\n",
       "      <td>0.449</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "      <td>1.594</td>\n",
       "      <td>0.573</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>37</td>\n",
       "      <td>UNK</td>\n",
       "      <td>-0.877</td>\n",
       "      <td>0.360</td>\n",
       "      <td>-0.143</td>\n",
       "      <td>0.401</td>\n",
       "      <td>0.151</td>\n",
       "      <td>0.362</td>\n",
       "      <td>0.489</td>\n",
       "      <td>0.378</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>0.151</td>\n",
       "      <td>0.362</td>\n",
       "      <td>0.489</td>\n",
       "      <td>0.378</td>\n",
       "      <td>0.151</td>\n",
       "      <td>0.362</td>\n",
       "      <td>0.489</td>\n",
       "      <td>0.378</td>\n",
       "      <td>0.489</td>\n",
       "      <td>0.378</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <td>38</td>\n",
       "      <td>S-LFM-H4-7</td>\n",
       "      <td>0.227</td>\n",
       "      <td>0.121</td>\n",
       "      <td>0.310</td>\n",
       "      <td>0.147</td>\n",
       "      <td>0.347</td>\n",
       "      <td>0.144</td>\n",
       "      <td>0.390</td>\n",
       "      <td>0.162</td>\n",
       "      <td>unique_id</td>\n",
       "      <td>...</td>\n",
       "      <td>0.347</td>\n",
       "      <td>0.144</td>\n",
       "      <td>0.390</td>\n",
       "      <td>0.162</td>\n",
       "      <td>0.347</td>\n",
       "      <td>0.144</td>\n",
       "      <td>0.390</td>\n",
       "      <td>0.162</td>\n",
       "      <td>0.390</td>\n",
       "      <td>0.162</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>39 rows × 26 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "               taxonomy  D1_mean  D1_stddev  D2_mean  D2_stddev  D3_mean  \\\n",
       "0                MUR-H1   -1.418      0.310   -0.709      0.328   -0.496   \n",
       "1             ER_ETR_H1   -1.418      0.310   -0.709      0.328   -0.496   \n",
       "2           ER_ETR_H1_2   -1.418      0.310   -0.709      0.328   -0.496   \n",
       "3              MUR-H1_3   -1.418      0.310   -0.709      0.328   -0.496   \n",
       "4          MUR_ADO_H1_2   -1.782      0.371   -0.894      0.326   -0.644   \n",
       "5        MUR-STRUB-H1-2   -1.502      0.295   -0.718      0.310   -0.492   \n",
       "6        MUR-STDRE-H1-2   -1.502      0.295   -0.718      0.310   -0.492   \n",
       "7           MR_DUC_H1_3   -0.532      0.407    0.251      0.394    0.364   \n",
       "8          MCF_DUC_H1_3   -0.532      0.407    0.251      0.394    0.364   \n",
       "9          MCF-DNO-H1-3   -0.804      0.307   -0.046      0.320    0.127   \n",
       "10          MR-DNO-H1-3   -0.804      0.307   -0.046      0.320    0.127   \n",
       "11      CR-LFM-DUC-H1-3    0.088      0.348    0.919      0.398    1.179   \n",
       "12     CR-LFLS-DUC-H1-3    0.088      0.348    0.919      0.398    1.179   \n",
       "13      CR-LFM-DUC-H4-7   -1.257      0.327   -0.419      0.350   -0.058   \n",
       "14     CR-LFLS-DUC-H4-7   -1.257      0.327   -0.419      0.350   -0.058   \n",
       "15    CR-LDUAL-DUC-H4-7   -1.257      0.327   -0.419      0.350   -0.058   \n",
       "16     CR-LFLS-DNO-H1-3   -0.785      0.333   -0.006      0.319   -0.256   \n",
       "17      CR-LFM-DNO-H1-3   -0.785      0.333   -0.006      0.319   -0.256   \n",
       "18      CR-LFM-DNO-H4-7   -1.841      0.292   -1.068      0.403   -0.717   \n",
       "19  CR-LFM-DNO-SOS-H1-3    0.820      0.317    0.054      0.334    0.205   \n",
       "20  CR-LFM-DNO-SOS-H4-7   -1.715      0.342   -1.187      0.423   -0.977   \n",
       "21    CR-LFINF-DUC-H1-3   -0.401      0.410    0.482      0.391    0.694   \n",
       "22     CR-LWAL-DUC-H1-3   -0.401      0.410    0.482      0.391    0.694   \n",
       "23      CR-PC-LWAL-H1-3   -0.401      0.410    0.482      0.391    0.694   \n",
       "24    CR-LFINF-DUC-H4-7   -2.023      0.465   -0.827      0.589   -0.307   \n",
       "25     CR-LWAL-DUC-H4-7   -2.023      0.465   -0.827      0.589   -0.307   \n",
       "26    CR-LFINF-DNO-H1-3   -1.126      0.363   -0.330      0.350   -0.119   \n",
       "27     CR-LWAL-DNO-H1-3   -1.126      0.363   -0.330      0.350   -0.119   \n",
       "28     CR-LWAL-DNO-H4-7   -2.401      0.362   -1.400      0.430   -0.981   \n",
       "29   CR-LDUAL-DUC-H8-19   -0.816      0.294   -0.260      0.361   -0.012   \n",
       "30    CR-LWAL-DUC-H8-19   -0.816      0.294   -0.260      0.361   -0.012   \n",
       "31             W-WLI-H1    0.107      0.403    0.853      0.450    1.144   \n",
       "32           W-WLI-H1-3    0.107      0.403    0.853      0.450    1.144   \n",
       "33           W-WHE-H1-3    0.107      0.403    0.853      0.450    1.144   \n",
       "34             W-WBB-H1    0.107      0.403    0.853      0.450    1.144   \n",
       "35              W-WS-H1    0.107      0.403    0.853      0.450    1.144   \n",
       "36            W-WS-H1-2    0.107      0.403    0.853      0.450    1.144   \n",
       "37                  UNK   -0.877      0.360   -0.143      0.401    0.151   \n",
       "38           S-LFM-H4-7    0.227      0.121    0.310      0.147    0.347   \n",
       "\n",
       "    D3_stddev  D4_mean  D4_stddev         id  ... D1_3_mean D1_3_stddev  \\\n",
       "0       0.322   -0.231      0.317  unique_id  ...    -0.496       0.322   \n",
       "1       0.322   -0.231      0.317  unique_id  ...    -0.496       0.322   \n",
       "2       0.322   -0.231      0.317  unique_id  ...    -0.496       0.322   \n",
       "3       0.322   -0.231      0.317  unique_id  ...    -0.496       0.322   \n",
       "4       0.328   -0.335      0.346  unique_id  ...    -0.644       0.328   \n",
       "5       0.311   -0.215      0.322  unique_id  ...    -0.492       0.311   \n",
       "6       0.311   -0.215      0.322  unique_id  ...    -0.492       0.311   \n",
       "7       0.377    0.578      0.382  unique_id  ...     0.364       0.377   \n",
       "8       0.377    0.578      0.382  unique_id  ...     0.364       0.377   \n",
       "9       0.314    0.386      0.367  unique_id  ...     0.127       0.314   \n",
       "10      0.314    0.386      0.367  unique_id  ...     0.127       0.314   \n",
       "11      0.381    1.616      0.507  unique_id  ...     1.179       0.381   \n",
       "12      0.381    1.616      0.507  unique_id  ...     1.179       0.381   \n",
       "13      0.404    0.436      0.525  unique_id  ...    -0.058       0.404   \n",
       "14      0.404    0.436      0.525  unique_id  ...    -0.058       0.404   \n",
       "15      0.404    0.436      0.525  unique_id  ...    -0.058       0.404   \n",
       "16      0.372   -0.643      0.437  unique_id  ...    -0.256       0.372   \n",
       "17      0.372   -0.643      0.437  unique_id  ...    -0.256       0.372   \n",
       "18      0.474   -0.195      0.591  unique_id  ...    -0.717       0.474   \n",
       "19      0.383    0.582      0.469  unique_id  ...     0.205       0.383   \n",
       "20      0.480   -0.638      0.566  unique_id  ...    -0.977       0.480   \n",
       "21      0.364    0.973      0.406  unique_id  ...     0.694       0.364   \n",
       "22      0.364    0.973      0.406  unique_id  ...     0.694       0.364   \n",
       "23      0.364    0.973      0.406  unique_id  ...     0.694       0.364   \n",
       "24      0.617    0.288      0.636  unique_id  ...    -0.307       0.617   \n",
       "25      0.617    0.288      0.636  unique_id  ...    -0.307       0.617   \n",
       "26      0.331    0.162      0.400  unique_id  ...    -0.119       0.331   \n",
       "27      0.331    0.162      0.400  unique_id  ...    -0.119       0.331   \n",
       "28      0.444   -0.402      0.548  unique_id  ...    -0.981       0.444   \n",
       "29      0.391    0.371      0.490  unique_id  ...    -0.012       0.391   \n",
       "30      0.391    0.371      0.490  unique_id  ...    -0.012       0.391   \n",
       "31      0.449    1.594      0.573  unique_id  ...     1.144       0.449   \n",
       "32      0.449    1.594      0.573  unique_id  ...     1.144       0.449   \n",
       "33      0.449    1.594      0.573  unique_id  ...     1.144       0.449   \n",
       "34      0.449    1.594      0.573  unique_id  ...     1.144       0.449   \n",
       "35      0.449    1.594      0.573  unique_id  ...     1.144       0.449   \n",
       "36      0.449    1.594      0.573  unique_id  ...     1.144       0.449   \n",
       "37      0.362    0.489      0.378  unique_id  ...     0.151       0.362   \n",
       "38      0.144    0.390      0.162  unique_id  ...     0.347       0.144   \n",
       "\n",
       "    D1_4_mean  D1_4_stddev  D2_3_mean  D2_3_stddev  D2_4_mean  D2_4_stddev  \\\n",
       "0      -0.231        0.317     -0.496        0.322     -0.231        0.317   \n",
       "1      -0.231        0.317     -0.496        0.322     -0.231        0.317   \n",
       "2      -0.231        0.317     -0.496        0.322     -0.231        0.317   \n",
       "3      -0.231        0.317     -0.496        0.322     -0.231        0.317   \n",
       "4      -0.335        0.346     -0.644        0.328     -0.335        0.346   \n",
       "5      -0.215        0.322     -0.492        0.311     -0.215        0.322   \n",
       "6      -0.215        0.322     -0.492        0.311     -0.215        0.322   \n",
       "7       0.578        0.382      0.364        0.377      0.578        0.382   \n",
       "8       0.578        0.382      0.364        0.377      0.578        0.382   \n",
       "9       0.386        0.367      0.127        0.314      0.386        0.367   \n",
       "10      0.386        0.367      0.127        0.314      0.386        0.367   \n",
       "11      1.616        0.507      1.179        0.381      1.616        0.507   \n",
       "12      1.616        0.507      1.179        0.381      1.616        0.507   \n",
       "13      0.436        0.525     -0.058        0.404      0.436        0.525   \n",
       "14      0.436        0.525     -0.058        0.404      0.436        0.525   \n",
       "15      0.436        0.525     -0.058        0.404      0.436        0.525   \n",
       "16     -0.643        0.437     -0.256        0.372     -0.643        0.437   \n",
       "17     -0.643        0.437     -0.256        0.372     -0.643        0.437   \n",
       "18     -0.195        0.591     -0.717        0.474     -0.195        0.591   \n",
       "19      0.582        0.469      0.205        0.383      0.582        0.469   \n",
       "20     -0.638        0.566     -0.977        0.480     -0.638        0.566   \n",
       "21      0.973        0.406      0.694        0.364      0.973        0.406   \n",
       "22      0.973        0.406      0.694        0.364      0.973        0.406   \n",
       "23      0.973        0.406      0.694        0.364      0.973        0.406   \n",
       "24      0.288        0.636     -0.307        0.617      0.288        0.636   \n",
       "25      0.288        0.636     -0.307        0.617      0.288        0.636   \n",
       "26      0.162        0.400     -0.119        0.331      0.162        0.400   \n",
       "27      0.162        0.400     -0.119        0.331      0.162        0.400   \n",
       "28     -0.402        0.548     -0.981        0.444     -0.402        0.548   \n",
       "29      0.371        0.490     -0.012        0.391      0.371        0.490   \n",
       "30      0.371        0.490     -0.012        0.391      0.371        0.490   \n",
       "31      1.594        0.573      1.144        0.449      1.594        0.573   \n",
       "32      1.594        0.573      1.144        0.449      1.594        0.573   \n",
       "33      1.594        0.573      1.144        0.449      1.594        0.573   \n",
       "34      1.594        0.573      1.144        0.449      1.594        0.573   \n",
       "35      1.594        0.573      1.144        0.449      1.594        0.573   \n",
       "36      1.594        0.573      1.144        0.449      1.594        0.573   \n",
       "37      0.489        0.378      0.151        0.362      0.489        0.378   \n",
       "38      0.390        0.162      0.347        0.144      0.390        0.162   \n",
       "\n",
       "    D3_4_mean  D3_4_stddev  \n",
       "0      -0.231        0.317  \n",
       "1      -0.231        0.317  \n",
       "2      -0.231        0.317  \n",
       "3      -0.231        0.317  \n",
       "4      -0.335        0.346  \n",
       "5      -0.215        0.322  \n",
       "6      -0.215        0.322  \n",
       "7       0.578        0.382  \n",
       "8       0.578        0.382  \n",
       "9       0.386        0.367  \n",
       "10      0.386        0.367  \n",
       "11      1.616        0.507  \n",
       "12      1.616        0.507  \n",
       "13      0.436        0.525  \n",
       "14      0.436        0.525  \n",
       "15      0.436        0.525  \n",
       "16     -0.643        0.437  \n",
       "17     -0.643        0.437  \n",
       "18     -0.195        0.591  \n",
       "19      0.582        0.469  \n",
       "20     -0.638        0.566  \n",
       "21      0.973        0.406  \n",
       "22      0.973        0.406  \n",
       "23      0.973        0.406  \n",
       "24      0.288        0.636  \n",
       "25      0.288        0.636  \n",
       "26      0.162        0.400  \n",
       "27      0.162        0.400  \n",
       "28     -0.402        0.548  \n",
       "29      0.371        0.490  \n",
       "30      0.371        0.490  \n",
       "31      1.594        0.573  \n",
       "32      1.594        0.573  \n",
       "33      1.594        0.573  \n",
       "34      1.594        0.573  \n",
       "35      1.594        0.573  \n",
       "36      1.594        0.573  \n",
       "37      0.489        0.378  \n",
       "38      0.390        0.162  \n",
       "\n",
       "[39 rows x 26 columns]"
      ]
     },
     "execution_count": 173,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "\n",
    "#tstfile=\"/Users/pittore/Documents/workspace/RIESGOS/GFZRIESGOS/updated_schemas/modelprop/SARA_v1.0_modelprop.json\"\n",
    "#pd.read_json(tstfile)\n",
    "#jsondata = json.load(open(tstfile,'r'))\n",
    "#df = pd.DataFrame(jsondata[\"data\"])\n",
    "df[\"D1_2_mean\"] = df[\"D2_mean\"]\n",
    "df[\"D1_2_stddev\"]= df[\"D2_stddev\"]\n",
    "\n",
    "df[\"D1_3_mean\"] = df[\"D3_mean\"]\n",
    "df[\"D1_3_stddev\"]= df[\"D3_stddev\"]\n",
    "\n",
    "df[\"D1_4_mean\"] = df[\"D4_mean\"]\n",
    "df[\"D1_4_stddev\"]= df[\"D4_stddev\"]\n",
    "\n",
    "df[\"D2_3_mean\"] = df[\"D3_mean\"]\n",
    "df[\"D2_3_stddev\"]= df[\"D3_stddev\"]\n",
    "\n",
    "df[\"D2_4_mean\"] = df[\"D4_mean\"]\n",
    "df[\"D2_4_stddev\"]= df[\"D4_stddev\"]\n",
    "\n",
    "df[\"D3_4_mean\"] = df[\"D4_mean\"]\n",
    "df[\"D3_4_stddev\"]= df[\"D4_stddev\"]\n",
    "\n",
    "df\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 174,
   "metadata": {},
   "outputs": [],
   "source": [
    "datastr = df.to_json(orient='records')\n",
    "jsondata[\"data\"] = eval(datastr)\n",
    "\n",
    "tstfile2=\"/Users/pittore/Documents/workspace/RIESGOS/GFZRIESGOS/updated_schemas/modelprop/SARA_v1.0_modelprop.json\"\n",
    "\n",
    "json.dump(jsondata,open(tstfile2,'w'))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}