Terrain variables analysis.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 235,
   "id": "aa736333",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "id": "105603f1",
   "metadata": {},
   "outputs": [],
   "source": [
    "area = pd.read_csv(\"parameters/area_terrain_statistics.csv\").dropna().drop(\"Unnamed: 0\", axis = 1).drop([''])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 149,
   "id": "eca5d165",
   "metadata": {
    "scrolled": true
   },
   "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>directory</th>\n",
       "      <th>terrain_variable</th>\n",
       "      <th>mean_value</th>\n",
       "      <th>variance_biased</th>\n",
       "      <th>variance_unbiased</th>\n",
       "      <th>standard_deviation_biased</th>\n",
       "      <th>standard_deviation_unbiased</th>\n",
       "      <th>max_value</th>\n",
       "      <th>min_value</th>\n",
       "      <th>median_value</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>dem</td>\n",
       "      <td>2.019486e+02</td>\n",
       "      <td>1.707965e+04</td>\n",
       "      <td>1.708020e+04</td>\n",
       "      <td>1.306891e+02</td>\n",
       "      <td>1.306912e+02</td>\n",
       "      <td>993.000000</td>\n",
       "      <td>8.200000e+01</td>\n",
       "      <td>1.570000e+02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>distance_stream_downstream</td>\n",
       "      <td>1.390115e+04</td>\n",
       "      <td>1.176266e+08</td>\n",
       "      <td>1.176311e+08</td>\n",
       "      <td>1.084558e+04</td>\n",
       "      <td>1.084579e+04</td>\n",
       "      <td>41508.400000</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>1.111416e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>EQIL_area</td>\n",
       "      <td>6.864827e-02</td>\n",
       "      <td>5.129056e-01</td>\n",
       "      <td>5.129221e-01</td>\n",
       "      <td>7.161743e-01</td>\n",
       "      <td>7.161858e-01</td>\n",
       "      <td>10.000000</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>ground_failure_estimate_cropped</td>\n",
       "      <td>4.470170e-04</td>\n",
       "      <td>3.800000e-07</td>\n",
       "      <td>3.800000e-07</td>\n",
       "      <td>6.160860e-04</td>\n",
       "      <td>6.160960e-04</td>\n",
       "      <td>0.015452</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>5.625600e-04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>local_relief</td>\n",
       "      <td>3.330574e+00</td>\n",
       "      <td>8.717205e+00</td>\n",
       "      <td>8.717494e+00</td>\n",
       "      <td>2.952491e+00</td>\n",
       "      <td>2.952540e+00</td>\n",
       "      <td>9.000000</td>\n",
       "      <td>1.000000e+00</td>\n",
       "      <td>1.000000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>profile_curvature</td>\n",
       "      <td>1.112838e-01</td>\n",
       "      <td>3.514898e+03</td>\n",
       "      <td>3.515013e+03</td>\n",
       "      <td>5.928657e+01</td>\n",
       "      <td>5.928754e+01</td>\n",
       "      <td>6686.076700</td>\n",
       "      <td>-6.033268e+03</td>\n",
       "      <td>-2.750000e-06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>slope</td>\n",
       "      <td>3.509013e+00</td>\n",
       "      <td>1.867733e+01</td>\n",
       "      <td>1.867794e+01</td>\n",
       "      <td>4.321728e+00</td>\n",
       "      <td>4.321798e+00</td>\n",
       "      <td>30.360462</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>1.857782e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>TWI_raster</td>\n",
       "      <td>-1.301636e+01</td>\n",
       "      <td>4.404442e+00</td>\n",
       "      <td>4.404584e+00</td>\n",
       "      <td>2.098676e+00</td>\n",
       "      <td>2.098710e+00</td>\n",
       "      <td>-4.960811</td>\n",
       "      <td>-1.782958e+01</td>\n",
       "      <td>-1.303805e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>vector_ruggedness_measure</td>\n",
       "      <td>1.722950e-03</td>\n",
       "      <td>1.210000e-05</td>\n",
       "      <td>1.210000e-05</td>\n",
       "      <td>3.480999e-03</td>\n",
       "      <td>3.481055e-03</td>\n",
       "      <td>0.047728</td>\n",
       "      <td>9.270000e-08</td>\n",
       "      <td>3.751010e-04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>dem</td>\n",
       "      <td>7.744444e+02</td>\n",
       "      <td>3.574136e+03</td>\n",
       "      <td>3.784379e+03</td>\n",
       "      <td>5.978408e+01</td>\n",
       "      <td>6.151731e+01</td>\n",
       "      <td>876.000000</td>\n",
       "      <td>6.560000e+02</td>\n",
       "      <td>7.745000e+02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>distance_stream_downstream</td>\n",
       "      <td>1.526995e+04</td>\n",
       "      <td>2.368485e+07</td>\n",
       "      <td>2.507808e+07</td>\n",
       "      <td>4.866709e+03</td>\n",
       "      <td>5.007802e+03</td>\n",
       "      <td>26276.666000</td>\n",
       "      <td>1.265530e+04</td>\n",
       "      <td>1.326696e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>EQIL_area</td>\n",
       "      <td>8.447705e-01</td>\n",
       "      <td>3.675252e+00</td>\n",
       "      <td>3.891443e+00</td>\n",
       "      <td>1.917095e+00</td>\n",
       "      <td>1.972674e+00</td>\n",
       "      <td>6.082763</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>local_relief</td>\n",
       "      <td>2.777778e+00</td>\n",
       "      <td>1.728395e-01</td>\n",
       "      <td>1.830065e-01</td>\n",
       "      <td>4.157397e-01</td>\n",
       "      <td>4.277926e-01</td>\n",
       "      <td>3.000000</td>\n",
       "      <td>2.000000e+00</td>\n",
       "      <td>3.000000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>pga_contour_raster_cropped</td>\n",
       "      <td>2.724381e-01</td>\n",
       "      <td>9.890000e-05</td>\n",
       "      <td>1.046980e-04</td>\n",
       "      <td>9.943904e-03</td>\n",
       "      <td>1.023219e-02</td>\n",
       "      <td>0.299562</td>\n",
       "      <td>2.600000e-01</td>\n",
       "      <td>2.693489e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>profile_curvature</td>\n",
       "      <td>6.400000e-07</td>\n",
       "      <td>3.510000e-13</td>\n",
       "      <td>3.710000e-13</td>\n",
       "      <td>5.920000e-07</td>\n",
       "      <td>6.090000e-07</td>\n",
       "      <td>0.000002</td>\n",
       "      <td>8.770000e-08</td>\n",
       "      <td>4.490000e-07</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>slope</td>\n",
       "      <td>2.442523e+01</td>\n",
       "      <td>6.276500e+00</td>\n",
       "      <td>6.645706e+00</td>\n",
       "      <td>2.505294e+00</td>\n",
       "      <td>2.577927e+00</td>\n",
       "      <td>28.023220</td>\n",
       "      <td>1.912525e+01</td>\n",
       "      <td>2.492959e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>TWI_raster</td>\n",
       "      <td>-1.570290e+01</td>\n",
       "      <td>4.878731e-01</td>\n",
       "      <td>5.165715e-01</td>\n",
       "      <td>6.984791e-01</td>\n",
       "      <td>7.187291e-01</td>\n",
       "      <td>-14.334189</td>\n",
       "      <td>-1.710681e+01</td>\n",
       "      <td>-1.575450e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>vector_ruggedness_measure</td>\n",
       "      <td>6.596914e-03</td>\n",
       "      <td>1.110000e-05</td>\n",
       "      <td>1.180000e-05</td>\n",
       "      <td>3.336264e-03</td>\n",
       "      <td>3.432988e-03</td>\n",
       "      <td>0.014293</td>\n",
       "      <td>2.736294e-03</td>\n",
       "      <td>5.140999e-03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>dem</td>\n",
       "      <td>7.205886e+02</td>\n",
       "      <td>1.523901e+04</td>\n",
       "      <td>1.528739e+04</td>\n",
       "      <td>1.234464e+02</td>\n",
       "      <td>1.236422e+02</td>\n",
       "      <td>981.000000</td>\n",
       "      <td>4.040000e+02</td>\n",
       "      <td>7.205000e+02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>distance_stream_downstream</td>\n",
       "      <td>2.262184e+04</td>\n",
       "      <td>8.892256e+06</td>\n",
       "      <td>8.920485e+06</td>\n",
       "      <td>2.981989e+03</td>\n",
       "      <td>2.986718e+03</td>\n",
       "      <td>27561.291000</td>\n",
       "      <td>1.263372e+04</td>\n",
       "      <td>2.140696e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>ground_failure_estimate_cropped</td>\n",
       "      <td>2.261018e-03</td>\n",
       "      <td>3.850000e-06</td>\n",
       "      <td>3.860000e-06</td>\n",
       "      <td>1.962442e-03</td>\n",
       "      <td>1.965554e-03</td>\n",
       "      <td>0.011195</td>\n",
       "      <td>6.240580e-04</td>\n",
       "      <td>1.263709e-03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>local_relief</td>\n",
       "      <td>4.183544e+00</td>\n",
       "      <td>3.181501e+00</td>\n",
       "      <td>3.191601e+00</td>\n",
       "      <td>1.783676e+00</td>\n",
       "      <td>1.786505e+00</td>\n",
       "      <td>9.000000</td>\n",
       "      <td>2.000000e+00</td>\n",
       "      <td>5.000000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>pga_contour_raster_cropped</td>\n",
       "      <td>2.906823e-01</td>\n",
       "      <td>1.627930e-04</td>\n",
       "      <td>1.633090e-04</td>\n",
       "      <td>1.275902e-02</td>\n",
       "      <td>1.277926e-02</td>\n",
       "      <td>0.319561</td>\n",
       "      <td>2.727869e-01</td>\n",
       "      <td>2.898045e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>profile_curvature</td>\n",
       "      <td>1.010000e-05</td>\n",
       "      <td>1.130000e-08</td>\n",
       "      <td>1.130000e-08</td>\n",
       "      <td>1.060680e-04</td>\n",
       "      <td>1.062360e-04</td>\n",
       "      <td>0.001765</td>\n",
       "      <td>-4.190000e-05</td>\n",
       "      <td>2.420000e-09</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>slope</td>\n",
       "      <td>1.430640e+01</td>\n",
       "      <td>1.631690e+01</td>\n",
       "      <td>1.636870e+01</td>\n",
       "      <td>4.039418e+00</td>\n",
       "      <td>4.045825e+00</td>\n",
       "      <td>25.352772</td>\n",
       "      <td>2.168196e+00</td>\n",
       "      <td>1.460791e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>TWI_raster</td>\n",
       "      <td>-1.510587e+01</td>\n",
       "      <td>2.425306e+00</td>\n",
       "      <td>2.433005e+00</td>\n",
       "      <td>1.557339e+00</td>\n",
       "      <td>1.559809e+00</td>\n",
       "      <td>-10.349417</td>\n",
       "      <td>-1.752993e+01</td>\n",
       "      <td>-1.543937e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>vector_ruggedness_measure</td>\n",
       "      <td>1.002198e-02</td>\n",
       "      <td>3.750000e-05</td>\n",
       "      <td>3.760000e-05</td>\n",
       "      <td>6.124666e-03</td>\n",
       "      <td>6.134380e-03</td>\n",
       "      <td>0.034200</td>\n",
       "      <td>1.212827e-03</td>\n",
       "      <td>8.357128e-03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>dem</td>\n",
       "      <td>1.664955e+02</td>\n",
       "      <td>1.002100e+04</td>\n",
       "      <td>1.002115e+04</td>\n",
       "      <td>1.001049e+02</td>\n",
       "      <td>1.001057e+02</td>\n",
       "      <td>993.000000</td>\n",
       "      <td>8.200000e+01</td>\n",
       "      <td>1.360000e+02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>distance_stream_downstream</td>\n",
       "      <td>1.569840e+04</td>\n",
       "      <td>1.554591e+08</td>\n",
       "      <td>1.554627e+08</td>\n",
       "      <td>1.246832e+04</td>\n",
       "      <td>1.246847e+04</td>\n",
       "      <td>52587.348000</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>1.239146e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>EQIL_area</td>\n",
       "      <td>3.194292e-02</td>\n",
       "      <td>2.398340e-01</td>\n",
       "      <td>2.398376e-01</td>\n",
       "      <td>4.897285e-01</td>\n",
       "      <td>4.897322e-01</td>\n",
       "      <td>10.000000</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>ground_failure_estimate_cropped</td>\n",
       "      <td>3.079860e-04</td>\n",
       "      <td>2.280000e-07</td>\n",
       "      <td>2.280000e-07</td>\n",
       "      <td>4.776440e-04</td>\n",
       "      <td>4.776480e-04</td>\n",
       "      <td>0.015452</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>0.000000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>35</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>local_relief</td>\n",
       "      <td>2.462099e+00</td>\n",
       "      <td>6.698878e+00</td>\n",
       "      <td>6.698989e+00</td>\n",
       "      <td>2.588219e+00</td>\n",
       "      <td>2.588240e+00</td>\n",
       "      <td>9.000000</td>\n",
       "      <td>1.000000e+00</td>\n",
       "      <td>1.000000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>37</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>pga_contour_raster_cropped</td>\n",
       "      <td>1.325230e-01</td>\n",
       "      <td>4.558253e-03</td>\n",
       "      <td>4.558322e-03</td>\n",
       "      <td>6.751483e-02</td>\n",
       "      <td>6.751534e-02</td>\n",
       "      <td>0.336039</td>\n",
       "      <td>4.299609e-02</td>\n",
       "      <td>1.162887e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>39</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>profile_curvature</td>\n",
       "      <td>-5.978387e-02</td>\n",
       "      <td>2.318865e+03</td>\n",
       "      <td>2.318902e+03</td>\n",
       "      <td>4.815460e+01</td>\n",
       "      <td>4.815498e+01</td>\n",
       "      <td>6686.076700</td>\n",
       "      <td>-6.033268e+03</td>\n",
       "      <td>-5.760000e-06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>slope</td>\n",
       "      <td>2.494904e+00</td>\n",
       "      <td>1.190447e+01</td>\n",
       "      <td>1.190466e+01</td>\n",
       "      <td>3.450286e+00</td>\n",
       "      <td>3.450312e+00</td>\n",
       "      <td>30.360462</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>1.084164e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>TWI_raster</td>\n",
       "      <td>-1.286445e+01</td>\n",
       "      <td>3.705568e+00</td>\n",
       "      <td>3.705623e+00</td>\n",
       "      <td>1.924985e+00</td>\n",
       "      <td>1.925000e+00</td>\n",
       "      <td>-4.960811</td>\n",
       "      <td>-1.782958e+01</td>\n",
       "      <td>-1.283767e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>vector_ruggedness_measure</td>\n",
       "      <td>1.053972e-03</td>\n",
       "      <td>6.980000e-06</td>\n",
       "      <td>6.980000e-06</td>\n",
       "      <td>2.641462e-03</td>\n",
       "      <td>2.641482e-03</td>\n",
       "      <td>0.047728</td>\n",
       "      <td>-2.110000e-08</td>\n",
       "      <td>1.072590e-04</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                 directory                 terrain_variable    mean_value  \\\n",
       "0                 PGA_zone                              dem  2.019486e+02   \n",
       "1                 PGA_zone       distance_stream_downstream  1.390115e+04   \n",
       "2                 PGA_zone                        EQIL_area  6.864827e-02   \n",
       "3                 PGA_zone  ground_failure_estimate_cropped  4.470170e-04   \n",
       "4                 PGA_zone                     local_relief  3.330574e+00   \n",
       "5                 PGA_zone                profile_curvature  1.112838e-01   \n",
       "6                 PGA_zone                            slope  3.509013e+00   \n",
       "7                 PGA_zone                       TWI_raster -1.301636e+01   \n",
       "8                 PGA_zone        vector_ruggedness_measure  1.722950e-03   \n",
       "9      ground_failure_zone                              dem  7.744444e+02   \n",
       "10     ground_failure_zone       distance_stream_downstream  1.526995e+04   \n",
       "11     ground_failure_zone                        EQIL_area  8.447705e-01   \n",
       "12     ground_failure_zone                     local_relief  2.777778e+00   \n",
       "13     ground_failure_zone       pga_contour_raster_cropped  2.724381e-01   \n",
       "14     ground_failure_zone                profile_curvature  6.400000e-07   \n",
       "15     ground_failure_zone                            slope  2.442523e+01   \n",
       "16     ground_failure_zone                       TWI_raster -1.570290e+01   \n",
       "17     ground_failure_zone        vector_ruggedness_measure  6.596914e-03   \n",
       "18  EQIL_distribution_area                              dem  7.205886e+02   \n",
       "19  EQIL_distribution_area       distance_stream_downstream  2.262184e+04   \n",
       "20  EQIL_distribution_area  ground_failure_estimate_cropped  2.261018e-03   \n",
       "21  EQIL_distribution_area                     local_relief  4.183544e+00   \n",
       "22  EQIL_distribution_area       pga_contour_raster_cropped  2.906823e-01   \n",
       "23  EQIL_distribution_area                profile_curvature  1.010000e-05   \n",
       "24  EQIL_distribution_area                            slope  1.430640e+01   \n",
       "25  EQIL_distribution_area                       TWI_raster -1.510587e+01   \n",
       "26  EQIL_distribution_area        vector_ruggedness_measure  1.002198e-02   \n",
       "27             binary_maps                              dem  1.664955e+02   \n",
       "29             binary_maps       distance_stream_downstream  1.569840e+04   \n",
       "31             binary_maps                        EQIL_area  3.194292e-02   \n",
       "33             binary_maps  ground_failure_estimate_cropped  3.079860e-04   \n",
       "35             binary_maps                     local_relief  2.462099e+00   \n",
       "37             binary_maps       pga_contour_raster_cropped  1.325230e-01   \n",
       "39             binary_maps                profile_curvature -5.978387e-02   \n",
       "41             binary_maps                            slope  2.494904e+00   \n",
       "43             binary_maps                       TWI_raster -1.286445e+01   \n",
       "45             binary_maps        vector_ruggedness_measure  1.053972e-03   \n",
       "\n",
       "    variance_biased  variance_unbiased  standard_deviation_biased  \\\n",
       "0      1.707965e+04       1.708020e+04               1.306891e+02   \n",
       "1      1.176266e+08       1.176311e+08               1.084558e+04   \n",
       "2      5.129056e-01       5.129221e-01               7.161743e-01   \n",
       "3      3.800000e-07       3.800000e-07               6.160860e-04   \n",
       "4      8.717205e+00       8.717494e+00               2.952491e+00   \n",
       "5      3.514898e+03       3.515013e+03               5.928657e+01   \n",
       "6      1.867733e+01       1.867794e+01               4.321728e+00   \n",
       "7      4.404442e+00       4.404584e+00               2.098676e+00   \n",
       "8      1.210000e-05       1.210000e-05               3.480999e-03   \n",
       "9      3.574136e+03       3.784379e+03               5.978408e+01   \n",
       "10     2.368485e+07       2.507808e+07               4.866709e+03   \n",
       "11     3.675252e+00       3.891443e+00               1.917095e+00   \n",
       "12     1.728395e-01       1.830065e-01               4.157397e-01   \n",
       "13     9.890000e-05       1.046980e-04               9.943904e-03   \n",
       "14     3.510000e-13       3.710000e-13               5.920000e-07   \n",
       "15     6.276500e+00       6.645706e+00               2.505294e+00   \n",
       "16     4.878731e-01       5.165715e-01               6.984791e-01   \n",
       "17     1.110000e-05       1.180000e-05               3.336264e-03   \n",
       "18     1.523901e+04       1.528739e+04               1.234464e+02   \n",
       "19     8.892256e+06       8.920485e+06               2.981989e+03   \n",
       "20     3.850000e-06       3.860000e-06               1.962442e-03   \n",
       "21     3.181501e+00       3.191601e+00               1.783676e+00   \n",
       "22     1.627930e-04       1.633090e-04               1.275902e-02   \n",
       "23     1.130000e-08       1.130000e-08               1.060680e-04   \n",
       "24     1.631690e+01       1.636870e+01               4.039418e+00   \n",
       "25     2.425306e+00       2.433005e+00               1.557339e+00   \n",
       "26     3.750000e-05       3.760000e-05               6.124666e-03   \n",
       "27     1.002100e+04       1.002115e+04               1.001049e+02   \n",
       "29     1.554591e+08       1.554627e+08               1.246832e+04   \n",
       "31     2.398340e-01       2.398376e-01               4.897285e-01   \n",
       "33     2.280000e-07       2.280000e-07               4.776440e-04   \n",
       "35     6.698878e+00       6.698989e+00               2.588219e+00   \n",
       "37     4.558253e-03       4.558322e-03               6.751483e-02   \n",
       "39     2.318865e+03       2.318902e+03               4.815460e+01   \n",
       "41     1.190447e+01       1.190466e+01               3.450286e+00   \n",
       "43     3.705568e+00       3.705623e+00               1.924985e+00   \n",
       "45     6.980000e-06       6.980000e-06               2.641462e-03   \n",
       "\n",
       "    standard_deviation_unbiased     max_value     min_value  median_value  \n",
       "0                  1.306912e+02    993.000000  8.200000e+01  1.570000e+02  \n",
       "1                  1.084579e+04  41508.400000  0.000000e+00  1.111416e+04  \n",
       "2                  7.161858e-01     10.000000  0.000000e+00  0.000000e+00  \n",
       "3                  6.160960e-04      0.015452  0.000000e+00  5.625600e-04  \n",
       "4                  2.952540e+00      9.000000  1.000000e+00  1.000000e+00  \n",
       "5                  5.928754e+01   6686.076700 -6.033268e+03 -2.750000e-06  \n",
       "6                  4.321798e+00     30.360462  0.000000e+00  1.857782e+00  \n",
       "7                  2.098710e+00     -4.960811 -1.782958e+01 -1.303805e+01  \n",
       "8                  3.481055e-03      0.047728  9.270000e-08  3.751010e-04  \n",
       "9                  6.151731e+01    876.000000  6.560000e+02  7.745000e+02  \n",
       "10                 5.007802e+03  26276.666000  1.265530e+04  1.326696e+04  \n",
       "11                 1.972674e+00      6.082763  0.000000e+00  0.000000e+00  \n",
       "12                 4.277926e-01      3.000000  2.000000e+00  3.000000e+00  \n",
       "13                 1.023219e-02      0.299562  2.600000e-01  2.693489e-01  \n",
       "14                 6.090000e-07      0.000002  8.770000e-08  4.490000e-07  \n",
       "15                 2.577927e+00     28.023220  1.912525e+01  2.492959e+01  \n",
       "16                 7.187291e-01    -14.334189 -1.710681e+01 -1.575450e+01  \n",
       "17                 3.432988e-03      0.014293  2.736294e-03  5.140999e-03  \n",
       "18                 1.236422e+02    981.000000  4.040000e+02  7.205000e+02  \n",
       "19                 2.986718e+03  27561.291000  1.263372e+04  2.140696e+04  \n",
       "20                 1.965554e-03      0.011195  6.240580e-04  1.263709e-03  \n",
       "21                 1.786505e+00      9.000000  2.000000e+00  5.000000e+00  \n",
       "22                 1.277926e-02      0.319561  2.727869e-01  2.898045e-01  \n",
       "23                 1.062360e-04      0.001765 -4.190000e-05  2.420000e-09  \n",
       "24                 4.045825e+00     25.352772  2.168196e+00  1.460791e+01  \n",
       "25                 1.559809e+00    -10.349417 -1.752993e+01 -1.543937e+01  \n",
       "26                 6.134380e-03      0.034200  1.212827e-03  8.357128e-03  \n",
       "27                 1.001057e+02    993.000000  8.200000e+01  1.360000e+02  \n",
       "29                 1.246847e+04  52587.348000  0.000000e+00  1.239146e+04  \n",
       "31                 4.897322e-01     10.000000  0.000000e+00  0.000000e+00  \n",
       "33                 4.776480e-04      0.015452  0.000000e+00  0.000000e+00  \n",
       "35                 2.588240e+00      9.000000  1.000000e+00  1.000000e+00  \n",
       "37                 6.751534e-02      0.336039  4.299609e-02  1.162887e-01  \n",
       "39                 4.815498e+01   6686.076700 -6.033268e+03 -5.760000e-06  \n",
       "41                 3.450312e+00     30.360462  0.000000e+00  1.084164e+00  \n",
       "43                 1.925000e+00     -4.960811 -1.782958e+01 -1.283767e+01  \n",
       "45                 2.641482e-03      0.047728 -2.110000e-08  1.072590e-04  "
      ]
     },
     "execution_count": 149,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "area"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 162,
   "id": "19b99daf",
   "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>directory</th>\n",
       "      <th>terrain_variable</th>\n",
       "      <th>mean_value</th>\n",
       "      <th>variance_biased</th>\n",
       "      <th>variance_unbiased</th>\n",
       "      <th>standard_deviation_biased</th>\n",
       "      <th>standard_deviation_unbiased</th>\n",
       "      <th>max_value</th>\n",
       "      <th>min_value</th>\n",
       "      <th>median_value</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>distance_stream_downstream</td>\n",
       "      <td>13901.147480</td>\n",
       "      <td>1.176266e+08</td>\n",
       "      <td>1.176311e+08</td>\n",
       "      <td>10845.578680</td>\n",
       "      <td>10845.787950</td>\n",
       "      <td>41508.400000</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>11114.157000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>local_relief</td>\n",
       "      <td>3.330574</td>\n",
       "      <td>8.717205e+00</td>\n",
       "      <td>8.717494e+00</td>\n",
       "      <td>2.952491</td>\n",
       "      <td>2.952540</td>\n",
       "      <td>9.000000</td>\n",
       "      <td>1.000000e+00</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>profile_curvature</td>\n",
       "      <td>0.111284</td>\n",
       "      <td>3.514898e+03</td>\n",
       "      <td>3.515013e+03</td>\n",
       "      <td>59.286571</td>\n",
       "      <td>59.287544</td>\n",
       "      <td>6686.076700</td>\n",
       "      <td>-6.033268e+03</td>\n",
       "      <td>-0.000003</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>slope</td>\n",
       "      <td>3.509013</td>\n",
       "      <td>1.867733e+01</td>\n",
       "      <td>1.867794e+01</td>\n",
       "      <td>4.321728</td>\n",
       "      <td>4.321798</td>\n",
       "      <td>30.360462</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>1.857782</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>TWI_raster</td>\n",
       "      <td>-13.016357</td>\n",
       "      <td>4.404442e+00</td>\n",
       "      <td>4.404584e+00</td>\n",
       "      <td>2.098676</td>\n",
       "      <td>2.098710</td>\n",
       "      <td>-4.960811</td>\n",
       "      <td>-1.782958e+01</td>\n",
       "      <td>-13.038050</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PGA_zone</td>\n",
       "      <td>vector_ruggedness_measure</td>\n",
       "      <td>0.001723</td>\n",
       "      <td>1.210000e-05</td>\n",
       "      <td>1.210000e-05</td>\n",
       "      <td>0.003481</td>\n",
       "      <td>0.003481</td>\n",
       "      <td>0.047728</td>\n",
       "      <td>9.270000e-08</td>\n",
       "      <td>0.000375</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  directory            terrain_variable    mean_value  variance_biased  \\\n",
       "1  PGA_zone  distance_stream_downstream  13901.147480     1.176266e+08   \n",
       "4  PGA_zone                local_relief      3.330574     8.717205e+00   \n",
       "5  PGA_zone           profile_curvature      0.111284     3.514898e+03   \n",
       "6  PGA_zone                       slope      3.509013     1.867733e+01   \n",
       "7  PGA_zone                  TWI_raster    -13.016357     4.404442e+00   \n",
       "8  PGA_zone   vector_ruggedness_measure      0.001723     1.210000e-05   \n",
       "\n",
       "   variance_unbiased  standard_deviation_biased  standard_deviation_unbiased  \\\n",
       "1       1.176311e+08               10845.578680                 10845.787950   \n",
       "4       8.717494e+00                   2.952491                     2.952540   \n",
       "5       3.515013e+03                  59.286571                    59.287544   \n",
       "6       1.867794e+01                   4.321728                     4.321798   \n",
       "7       4.404584e+00                   2.098676                     2.098710   \n",
       "8       1.210000e-05                   0.003481                     0.003481   \n",
       "\n",
       "      max_value     min_value  median_value  \n",
       "1  41508.400000  0.000000e+00  11114.157000  \n",
       "4      9.000000  1.000000e+00      1.000000  \n",
       "5   6686.076700 -6.033268e+03     -0.000003  \n",
       "6     30.360462  0.000000e+00      1.857782  \n",
       "7     -4.960811 -1.782958e+01    -13.038050  \n",
       "8      0.047728  9.270000e-08      0.000375  "
      ]
     },
     "execution_count": 162,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "PGA_zone = area[area['directory'] == 'PGA_zone']\n",
    "PGA_zone = PGA_zone[PGA_zone['terrain_variable'] != 'dem']\n",
    "PGA_zone = PGA_zone[PGA_zone['terrain_variable'] != 'EQIL_area']\n",
    "PGA_zone = PGA_zone[PGA_zone['terrain_variable'] != 'ground_failure_estimate_cropped']\n",
    "\n",
    "PGA_zone"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 166,
   "id": "e4c7c390",
   "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>directory</th>\n",
       "      <th>terrain_variable</th>\n",
       "      <th>mean_value</th>\n",
       "      <th>variance_biased</th>\n",
       "      <th>variance_unbiased</th>\n",
       "      <th>standard_deviation_biased</th>\n",
       "      <th>standard_deviation_unbiased</th>\n",
       "      <th>max_value</th>\n",
       "      <th>min_value</th>\n",
       "      <th>median_value</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>distance_stream_downstream</td>\n",
       "      <td>1.526995e+04</td>\n",
       "      <td>2.368485e+07</td>\n",
       "      <td>2.507808e+07</td>\n",
       "      <td>4.866709e+03</td>\n",
       "      <td>5.007802e+03</td>\n",
       "      <td>26276.666000</td>\n",
       "      <td>1.265530e+04</td>\n",
       "      <td>1.326696e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>local_relief</td>\n",
       "      <td>2.777778e+00</td>\n",
       "      <td>1.728395e-01</td>\n",
       "      <td>1.830065e-01</td>\n",
       "      <td>4.157397e-01</td>\n",
       "      <td>4.277926e-01</td>\n",
       "      <td>3.000000</td>\n",
       "      <td>2.000000e+00</td>\n",
       "      <td>3.000000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>profile_curvature</td>\n",
       "      <td>6.400000e-07</td>\n",
       "      <td>3.510000e-13</td>\n",
       "      <td>3.710000e-13</td>\n",
       "      <td>5.920000e-07</td>\n",
       "      <td>6.090000e-07</td>\n",
       "      <td>0.000002</td>\n",
       "      <td>8.770000e-08</td>\n",
       "      <td>4.490000e-07</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>slope</td>\n",
       "      <td>2.442523e+01</td>\n",
       "      <td>6.276500e+00</td>\n",
       "      <td>6.645706e+00</td>\n",
       "      <td>2.505294e+00</td>\n",
       "      <td>2.577927e+00</td>\n",
       "      <td>28.023220</td>\n",
       "      <td>1.912525e+01</td>\n",
       "      <td>2.492959e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>TWI_raster</td>\n",
       "      <td>-1.570290e+01</td>\n",
       "      <td>4.878731e-01</td>\n",
       "      <td>5.165715e-01</td>\n",
       "      <td>6.984791e-01</td>\n",
       "      <td>7.187291e-01</td>\n",
       "      <td>-14.334189</td>\n",
       "      <td>-1.710681e+01</td>\n",
       "      <td>-1.575450e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>ground_failure_zone</td>\n",
       "      <td>vector_ruggedness_measure</td>\n",
       "      <td>6.596914e-03</td>\n",
       "      <td>1.110000e-05</td>\n",
       "      <td>1.180000e-05</td>\n",
       "      <td>3.336264e-03</td>\n",
       "      <td>3.432988e-03</td>\n",
       "      <td>0.014293</td>\n",
       "      <td>2.736294e-03</td>\n",
       "      <td>5.140999e-03</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "              directory            terrain_variable    mean_value  \\\n",
       "10  ground_failure_zone  distance_stream_downstream  1.526995e+04   \n",
       "12  ground_failure_zone                local_relief  2.777778e+00   \n",
       "14  ground_failure_zone           profile_curvature  6.400000e-07   \n",
       "15  ground_failure_zone                       slope  2.442523e+01   \n",
       "16  ground_failure_zone                  TWI_raster -1.570290e+01   \n",
       "17  ground_failure_zone   vector_ruggedness_measure  6.596914e-03   \n",
       "\n",
       "    variance_biased  variance_unbiased  standard_deviation_biased  \\\n",
       "10     2.368485e+07       2.507808e+07               4.866709e+03   \n",
       "12     1.728395e-01       1.830065e-01               4.157397e-01   \n",
       "14     3.510000e-13       3.710000e-13               5.920000e-07   \n",
       "15     6.276500e+00       6.645706e+00               2.505294e+00   \n",
       "16     4.878731e-01       5.165715e-01               6.984791e-01   \n",
       "17     1.110000e-05       1.180000e-05               3.336264e-03   \n",
       "\n",
       "    standard_deviation_unbiased     max_value     min_value  median_value  \n",
       "10                 5.007802e+03  26276.666000  1.265530e+04  1.326696e+04  \n",
       "12                 4.277926e-01      3.000000  2.000000e+00  3.000000e+00  \n",
       "14                 6.090000e-07      0.000002  8.770000e-08  4.490000e-07  \n",
       "15                 2.577927e+00     28.023220  1.912525e+01  2.492959e+01  \n",
       "16                 7.187291e-01    -14.334189 -1.710681e+01 -1.575450e+01  \n",
       "17                 3.432988e-03      0.014293  2.736294e-03  5.140999e-03  "
      ]
     },
     "execution_count": 166,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ground_failure_zone = area[area['directory'] == 'ground_failure_zone']\n",
    "ground_failure_zone = ground_failure_zone[ground_failure_zone['terrain_variable'] != 'dem']\n",
    "ground_failure_zone = ground_failure_zone[ground_failure_zone['terrain_variable'] != 'pga_contour_raster_cropped']\n",
    "ground_failure_zone = ground_failure_zone[ground_failure_zone['terrain_variable'] != 'EQIL_area']\n",
    "ground_failure_zone"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 167,
   "id": "22f4af5e",
   "metadata": {
    "scrolled": true
   },
   "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>directory</th>\n",
       "      <th>terrain_variable</th>\n",
       "      <th>mean_value</th>\n",
       "      <th>variance_biased</th>\n",
       "      <th>variance_unbiased</th>\n",
       "      <th>standard_deviation_biased</th>\n",
       "      <th>standard_deviation_unbiased</th>\n",
       "      <th>max_value</th>\n",
       "      <th>min_value</th>\n",
       "      <th>median_value</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>distance_stream_downstream</td>\n",
       "      <td>22621.837090</td>\n",
       "      <td>8.892256e+06</td>\n",
       "      <td>8.920485e+06</td>\n",
       "      <td>2981.988599</td>\n",
       "      <td>2986.718163</td>\n",
       "      <td>27561.291000</td>\n",
       "      <td>12633.724000</td>\n",
       "      <td>2.140696e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>local_relief</td>\n",
       "      <td>4.183544</td>\n",
       "      <td>3.181501e+00</td>\n",
       "      <td>3.191601e+00</td>\n",
       "      <td>1.783676</td>\n",
       "      <td>1.786505</td>\n",
       "      <td>9.000000</td>\n",
       "      <td>2.000000</td>\n",
       "      <td>5.000000e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>profile_curvature</td>\n",
       "      <td>0.000010</td>\n",
       "      <td>1.130000e-08</td>\n",
       "      <td>1.130000e-08</td>\n",
       "      <td>0.000106</td>\n",
       "      <td>0.000106</td>\n",
       "      <td>0.001765</td>\n",
       "      <td>-0.000042</td>\n",
       "      <td>2.420000e-09</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>slope</td>\n",
       "      <td>14.306402</td>\n",
       "      <td>1.631690e+01</td>\n",
       "      <td>1.636870e+01</td>\n",
       "      <td>4.039418</td>\n",
       "      <td>4.045825</td>\n",
       "      <td>25.352772</td>\n",
       "      <td>2.168196</td>\n",
       "      <td>1.460791e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>TWI_raster</td>\n",
       "      <td>-15.105869</td>\n",
       "      <td>2.425306e+00</td>\n",
       "      <td>2.433005e+00</td>\n",
       "      <td>1.557339</td>\n",
       "      <td>1.559809</td>\n",
       "      <td>-10.349417</td>\n",
       "      <td>-17.529930</td>\n",
       "      <td>-1.543937e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>EQIL_distribution_area</td>\n",
       "      <td>vector_ruggedness_measure</td>\n",
       "      <td>0.010022</td>\n",
       "      <td>3.750000e-05</td>\n",
       "      <td>3.760000e-05</td>\n",
       "      <td>0.006125</td>\n",
       "      <td>0.006134</td>\n",
       "      <td>0.034200</td>\n",
       "      <td>0.001213</td>\n",
       "      <td>8.357128e-03</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                 directory            terrain_variable    mean_value  \\\n",
       "19  EQIL_distribution_area  distance_stream_downstream  22621.837090   \n",
       "21  EQIL_distribution_area                local_relief      4.183544   \n",
       "23  EQIL_distribution_area           profile_curvature      0.000010   \n",
       "24  EQIL_distribution_area                       slope     14.306402   \n",
       "25  EQIL_distribution_area                  TWI_raster    -15.105869   \n",
       "26  EQIL_distribution_area   vector_ruggedness_measure      0.010022   \n",
       "\n",
       "    variance_biased  variance_unbiased  standard_deviation_biased  \\\n",
       "19     8.892256e+06       8.920485e+06                2981.988599   \n",
       "21     3.181501e+00       3.191601e+00                   1.783676   \n",
       "23     1.130000e-08       1.130000e-08                   0.000106   \n",
       "24     1.631690e+01       1.636870e+01                   4.039418   \n",
       "25     2.425306e+00       2.433005e+00                   1.557339   \n",
       "26     3.750000e-05       3.760000e-05                   0.006125   \n",
       "\n",
       "    standard_deviation_unbiased     max_value     min_value  median_value  \n",
       "19                  2986.718163  27561.291000  12633.724000  2.140696e+04  \n",
       "21                     1.786505      9.000000      2.000000  5.000000e+00  \n",
       "23                     0.000106      0.001765     -0.000042  2.420000e-09  \n",
       "24                     4.045825     25.352772      2.168196  1.460791e+01  \n",
       "25                     1.559809    -10.349417    -17.529930 -1.543937e+01  \n",
       "26                     0.006134      0.034200      0.001213  8.357128e-03  "
      ]
     },
     "execution_count": 167,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "EQIL_distribution_area = area[area['directory'] == 'EQIL_distribution_area']\n",
    "EQIL_distribution_area = EQIL_distribution_area[EQIL_distribution_area['terrain_variable'] != 'dem']\n",
    "EQIL_distribution_area = EQIL_distribution_area[EQIL_distribution_area['terrain_variable'] != 'pga_contour_raster_cropped']\n",
    "EQIL_distribution_area = EQIL_distribution_area[EQIL_distribution_area['terrain_variable'] != 'ground_failure_estimate_cropped']\n",
    "EQIL_distribution_area"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 168,
   "id": "5cd93788",
   "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>directory</th>\n",
       "      <th>terrain_variable</th>\n",
       "      <th>mean_value</th>\n",
       "      <th>variance_biased</th>\n",
       "      <th>variance_unbiased</th>\n",
       "      <th>standard_deviation_biased</th>\n",
       "      <th>standard_deviation_unbiased</th>\n",
       "      <th>max_value</th>\n",
       "      <th>min_value</th>\n",
       "      <th>median_value</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>distance_stream_downstream</td>\n",
       "      <td>15698.401450</td>\n",
       "      <td>1.554591e+08</td>\n",
       "      <td>1.554627e+08</td>\n",
       "      <td>12468.323720</td>\n",
       "      <td>12468.469430</td>\n",
       "      <td>52587.348000</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>12391.461000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>35</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>local_relief</td>\n",
       "      <td>2.462099</td>\n",
       "      <td>6.698878e+00</td>\n",
       "      <td>6.698989e+00</td>\n",
       "      <td>2.588219</td>\n",
       "      <td>2.588240</td>\n",
       "      <td>9.000000</td>\n",
       "      <td>1.000000e+00</td>\n",
       "      <td>1.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>39</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>profile_curvature</td>\n",
       "      <td>-0.059784</td>\n",
       "      <td>2.318865e+03</td>\n",
       "      <td>2.318902e+03</td>\n",
       "      <td>48.154596</td>\n",
       "      <td>48.154982</td>\n",
       "      <td>6686.076700</td>\n",
       "      <td>-6.033268e+03</td>\n",
       "      <td>-0.000006</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>slope</td>\n",
       "      <td>2.494904</td>\n",
       "      <td>1.190447e+01</td>\n",
       "      <td>1.190466e+01</td>\n",
       "      <td>3.450286</td>\n",
       "      <td>3.450312</td>\n",
       "      <td>30.360462</td>\n",
       "      <td>0.000000e+00</td>\n",
       "      <td>1.084164</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>TWI_raster</td>\n",
       "      <td>-12.864450</td>\n",
       "      <td>3.705568e+00</td>\n",
       "      <td>3.705623e+00</td>\n",
       "      <td>1.924985</td>\n",
       "      <td>1.925000</td>\n",
       "      <td>-4.960811</td>\n",
       "      <td>-1.782958e+01</td>\n",
       "      <td>-12.837674</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>binary_maps</td>\n",
       "      <td>vector_ruggedness_measure</td>\n",
       "      <td>0.001054</td>\n",
       "      <td>6.980000e-06</td>\n",
       "      <td>6.980000e-06</td>\n",
       "      <td>0.002641</td>\n",
       "      <td>0.002641</td>\n",
       "      <td>0.047728</td>\n",
       "      <td>-2.110000e-08</td>\n",
       "      <td>0.000107</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      directory            terrain_variable    mean_value  variance_biased  \\\n",
       "29  binary_maps  distance_stream_downstream  15698.401450     1.554591e+08   \n",
       "35  binary_maps                local_relief      2.462099     6.698878e+00   \n",
       "39  binary_maps           profile_curvature     -0.059784     2.318865e+03   \n",
       "41  binary_maps                       slope      2.494904     1.190447e+01   \n",
       "43  binary_maps                  TWI_raster    -12.864450     3.705568e+00   \n",
       "45  binary_maps   vector_ruggedness_measure      0.001054     6.980000e-06   \n",
       "\n",
       "    variance_unbiased  standard_deviation_biased  standard_deviation_unbiased  \\\n",
       "29       1.554627e+08               12468.323720                 12468.469430   \n",
       "35       6.698989e+00                   2.588219                     2.588240   \n",
       "39       2.318902e+03                  48.154596                    48.154982   \n",
       "41       1.190466e+01                   3.450286                     3.450312   \n",
       "43       3.705623e+00                   1.924985                     1.925000   \n",
       "45       6.980000e-06                   0.002641                     0.002641   \n",
       "\n",
       "       max_value     min_value  median_value  \n",
       "29  52587.348000  0.000000e+00  12391.461000  \n",
       "35      9.000000  1.000000e+00      1.000000  \n",
       "39   6686.076700 -6.033268e+03     -0.000006  \n",
       "41     30.360462  0.000000e+00      1.084164  \n",
       "43     -4.960811 -1.782958e+01    -12.837674  \n",
       "45      0.047728 -2.110000e-08      0.000107  "
      ]
     },
     "execution_count": 168,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "binary_maps = area[area['directory'] == 'binary_maps']\n",
    "binary_maps = binary_maps[binary_maps['terrain_variable'] != 'dem']\n",
    "binary_maps = binary_maps[binary_maps['terrain_variable'] != 'EQIL_area']\n",
    "binary_maps = binary_maps[binary_maps['terrain_variable'] != 'ground_failure_estimate_cropped']\n",
    "binary_maps = binary_maps[binary_maps['terrain_variable'] != 'pga_contour_raster_cropped']\n",
    "\n",
    "binary_maps"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "760bc14b",
   "metadata": {},
   "source": [
    "We can discard dem, EQIL area, PGA\n",
    "We are interested in distance_stream_downstream, local_relief, profile_curvature, slope, TWI, VRM"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 293,
   "id": "e22486b2",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(5    0.111284\n",
       " Name: mean_value, dtype: float64,\n",
       " 14    6.400000e-07\n",
       " Name: mean_value, dtype: float64,\n",
       " 23    0.00001\n",
       " Name: mean_value, dtype: float64,\n",
       " 39   -0.059784\n",
       " Name: mean_value, dtype: float64)"
      ]
     },
     "execution_count": 293,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "PGA_zone_profile_curvature = PGA_zone[PGA_zone['terrain_variable'] == 'profile_curvature'].mean_value\n",
    "ground_failure_zone_profile_curvature = ground_failure_zone[ground_failure_zone['terrain_variable'] == 'profile_curvature'].mean_value\n",
    "EQIL_distribution_area_profile_curvature = EQIL_distribution_area[EQIL_distribution_area['terrain_variable'] == 'profile_curvature'].mean_value\n",
    "binary_maps_profile_curvature = binary_maps[binary_maps['terrain_variable'] == 'profile_curvature'].mean_value\n",
    "PGA_zone_profile_curvature, ground_failure_zone_profile_curvature, EQIL_distribution_area_profile_curvature, binary_maps_profile_curvature"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "41dd4fd4",
   "metadata": {},
   "source": [
    "Profile curvature is a measure describing the concavity/convexity of slope along the vertical direction. A positive value indicates concave slope, while a negative value indicated a convex one. Having a concave surface can increase slope instability by increasing the subsurface drainage that can cause high water pressure (Pierson 1980). By analyzing the profile curvature mean we can see on average the slope stability. Given the results above, we can see tha the PGA_zone, ground_failure_Zone and EQIL_distribution area are concave, although the latter two are extremely small compared to the first one. On the other hand, the binary map, which is the entire area, seems to be slightly convex"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 299,
   "id": "22f04da0",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(8    9.270000e-08\n",
       " Name: min_value, dtype: float64,\n",
       " 17    0.002736\n",
       " Name: min_value, dtype: float64,\n",
       " 26    0.001213\n",
       " Name: min_value, dtype: float64,\n",
       " 45   -2.110000e-08\n",
       " Name: min_value, dtype: float64)"
      ]
     },
     "execution_count": 299,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "PGA_zone_VRM2 = PGA_zone[PGA_zone['terrain_variable'] == 'vector_ruggedness_measure'].min_value\n",
    "ground_failure_zone_VRM2  = ground_failure_zone[ground_failure_zone['terrain_variable'] == 'vector_ruggedness_measure'].min_value\n",
    "EQIL_distribution_area_VRM2  = EQIL_distribution_area[EQIL_distribution_area['terrain_variable'] == 'vector_ruggedness_measure'].min_value\n",
    "binary_maps_VRM2  = binary_maps[binary_maps['terrain_variable'] == 'vector_ruggedness_measure'].min_value\n",
    "\n",
    "PGA_zone_VRM2 , ground_failure_zone_VRM2 , EQIL_distribution_area_VRM2 , binary_maps_VRM2 "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d5d8d177",
   "metadata": {},
   "source": [
    "Tanyas et al. (2019) showed that the majority of EQIL are initiated at low VRM values. Furthemore, the lower the VRM value is then the chances of a landslide occuring is higher. In this case, if one area has drastically lower VRM than the others then it means that area is most likely to have an earthquake occur. Thus, we are looking for the lowest value possible in all areas/maps. As observed, the binary map is the only one with a negative value, although the others are positive but still very close to 0. The discrepancy is not too considerable, but it shows the ara most likely to have a landslide occur"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 301,
   "id": "122fa147",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(4    9.0\n",
       " Name: max_value, dtype: float64,\n",
       " 12    3.0\n",
       " Name: max_value, dtype: float64,\n",
       " 21    9.0\n",
       " Name: max_value, dtype: float64,\n",
       " 35    9.0\n",
       " Name: max_value, dtype: float64)"
      ]
     },
     "execution_count": 301,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "PGA_zone_local = PGA_zone[PGA_zone['terrain_variable'] == 'local_relief'].max_value\n",
    "ground_failure_zone_local  = ground_failure_zone[ground_failure_zone['terrain_variable'] == 'local_relief'].max_value\n",
    "EQIL_distribution_area_local  = EQIL_distribution_area[EQIL_distribution_area['terrain_variable'] == 'local_relief'].max_value\n",
    "binary_maps_local  = binary_maps[binary_maps['terrain_variable'] == 'local_relief'].max_value\n",
    "\n",
    "PGA_zone_local , ground_failure_zone_local , EQIL_distribution_area_local , binary_maps_local "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "eea9930c",
   "metadata": {},
   "source": [
    "Local relied is the maximum difference in height in a local neighborhood of each pixel and can be related to slope instability caused by tectonic uplift. Thus, the maximum value in those maps indicates the maximum difference in height and a higher slope instablity. As it can be seen, all  maps have the same value of 9.0, with exception of ground_failure_zone_local which is quite lower than the others at 3.0 . "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 303,
   "id": "0f35c4dd",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(7    2.09871\n",
       " Name: standard_deviation_unbiased, dtype: float64,\n",
       " 16    0.718729\n",
       " Name: standard_deviation_unbiased, dtype: float64,\n",
       " 25    1.559809\n",
       " Name: standard_deviation_unbiased, dtype: float64,\n",
       " 43    1.925\n",
       " Name: standard_deviation_unbiased, dtype: float64)"
      ]
     },
     "execution_count": 303,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "PGA_zone_TWI = PGA_zone[PGA_zone['terrain_variable'] == 'TWI_raster'].standard_deviation_unbiased\n",
    "ground_failure_zone_TWI  = ground_failure_zone[ground_failure_zone['terrain_variable'] == 'TWI_raster'].standard_deviation_unbiased\n",
    "EQIL_distribution_area_TWI  = EQIL_distribution_area[EQIL_distribution_area['terrain_variable'] == 'TWI_raster'].standard_deviation_unbiased\n",
    "binary_maps_TWI  = binary_maps[binary_maps['terrain_variable'] == 'TWI_raster'].standard_deviation_unbiased\n",
    "\n",
    "PGA_zone_TWI , ground_failure_zone_TWI , EQIL_distribution_area_TWI , binary_maps_TWI"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "17cf1e3a",
   "metadata": {},
   "source": [
    "TWI (Moore et al. 1991) is a proxy for potential soil wetness used to estimate the spatial variability of wetness within a landscape (e.g., Nowicki Jessee et al. 2018). It can take a role in slope stability by changing the pore water pressure. By looking at the standard deviation of the TWI of our maps, we can see the values are quite high in three of the cases, which means the values are quite scattered around the means in great distances."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 305,
   "id": "079b637d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(6    1.857782\n",
       " Name: median_value, dtype: float64,\n",
       " 15    24.929588\n",
       " Name: median_value, dtype: float64,\n",
       " 24    14.60791\n",
       " Name: median_value, dtype: float64,\n",
       " 41    1.084164\n",
       " Name: median_value, dtype: float64)"
      ]
     },
     "execution_count": 305,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "PGA_zone_slope = PGA_zone[PGA_zone['terrain_variable'] == 'slope'].median_value\n",
    "ground_failure_zone_slope  = ground_failure_zone[ground_failure_zone['terrain_variable'] == 'slope'].median_value\n",
    "EQIL_distribution_area_slope  = EQIL_distribution_area[EQIL_distribution_area['terrain_variable'] == 'slope'].median_value\n",
    "binary_maps_slope  = binary_maps[binary_maps['terrain_variable'] == 'slope'].median_value\n",
    "\n",
    "PGA_zone_slope , ground_failure_zone_slope , EQIL_distribution_area_slope , binary_maps_slope"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ea9c6cdc",
   "metadata": {},
   "source": [
    "Slope is a factor controlling the normal and shear stresses, which take a role in slope stability. Budimir et al. (2015) stated that in all those studies slope was found as a significant variable. There is quite a discrepancy of the median value between the four maps. The PGA_zone and binary_maps seems to have quite similar values between 1 and 2, whereas the EQIL_distribution_area has higher values at around the 14 mark, followed by the ground_failure_zone with the highest value of almost 25. That makes sense for the ground_failure_zone since the slope is a very important parameter and a higher slope indicated higher chances of a landslide."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 307,
   "id": "f3c641db",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(1    13901.14748\n",
       " Name: mean_value, dtype: float64,\n",
       " 10    15269.94775\n",
       " Name: mean_value, dtype: float64,\n",
       " 19    22621.83709\n",
       " Name: mean_value, dtype: float64,\n",
       " 29    15698.40145\n",
       " Name: mean_value, dtype: float64)"
      ]
     },
     "execution_count": 307,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "PGA_zone_distance = PGA_zone[PGA_zone['terrain_variable'] == 'distance_stream_downstream'].mean_value\n",
    "ground_failure_zone_distance  = ground_failure_zone[ground_failure_zone['terrain_variable'] == 'distance_stream_downstream'].mean_value\n",
    "EQIL_distribution_area_distance = EQIL_distribution_area[EQIL_distribution_area['terrain_variable'] == 'distance_stream_downstream'].mean_value\n",
    "binary_maps_distance  = binary_maps[binary_maps['terrain_variable'] == 'distance_stream_downstream'].mean_value\n",
    "\n",
    "PGA_zone_distance , ground_failure_zone_distance , EQIL_distribution_area_distance , binary_maps_distance"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "cc4639fd",
   "metadata": {},
   "source": [
    "Distance to stream is proxy related to fluvial undercutting (e.g., Kritikos et al. 2015) that cause high rates of shear stress as a result of loss of lateral support (Korup 2004). Tanyaş et al. (2017) showed that the majority of EQIL are initiated close to river channels and the frequency of observed landslides gradually decreases while going far away from channels. The lower the distance_stream_downstream, the lower the distance to rive channels which results in higher chances of occuring landslides.The PGA_zone has the lowest distance_stream_downstream, closely followed by he ground_failure_zone and binary_maps which have very similar. Those have considerably higher chances of occuring landslidse compared to the EQIL_distribution_area as they are closer to river channels."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.8.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}