05_Procesamiento_Escenario_Futuro_RCP8.5.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# ESCENARIO PRESENTE (2010)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# PAQUETES"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import rasterio\n",
    "import geopandas as gpd\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "from rasterio.mask import mask\n",
    "from shapely.geometry import box\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# CARGAR DATOS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "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>ET_ID</th>\n",
       "      <th>COUNTRY</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>Malta</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3101</th>\n",
       "      <td>9369</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((26.92771 70.58708, 25.63446 70.69412...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3102</th>\n",
       "      <td>9370</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.20775 70.47181, 26.92771 70.58708...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3103</th>\n",
       "      <td>9371</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((29.47384 70.34844, 28.20775 70.47181...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3104</th>\n",
       "      <td>9372</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((30.72534 70.21713, 29.47384 70.34844...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3105</th>\n",
       "      <td>9485</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.60702 70.90647, 27.30206 71.02446...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3106 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      ET_ID COUNTRY                                           geometry\n",
       "0       249   Spain  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...\n",
       "1       280   Malta  POLYGON ((14.68659 35.53849, 14.14122 35.56449...\n",
       "2       363   Spain  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...\n",
       "3       412  Greece  POLYGON ((23.95653 35.05738, 23.41935 35.13773...\n",
       "4       413  Greece  POLYGON ((24.49271 34.97391, 23.95653 35.05738...\n",
       "...     ...     ...                                                ...\n",
       "3101   9369  Norway  POLYGON ((26.92771 70.58708, 25.63446 70.69412...\n",
       "3102   9370  Norway  POLYGON ((28.20775 70.47181, 26.92771 70.58708...\n",
       "3103   9371  Norway  POLYGON ((29.47384 70.34844, 28.20775 70.47181...\n",
       "3104   9372  Norway  POLYGON ((30.72534 70.21713, 29.47384 70.34844...\n",
       "3105   9485  Norway  POLYGON ((28.60702 70.90647, 27.30206 71.02446...\n",
       "\n",
       "[3106 rows x 3 columns]"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#Cargar el grid\n",
    "Grid_50 = gpd.read_file(r'Data/Grid50/Grid_50Km_Europe_WGS84.shp')\n",
    "Grid_50"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Abrir el raster y leer los datos fuera del bucle\n",
    "src = rasterio.open(r'Data/GLOBIO/RCP8_5/Globio4_landuse_10sec_2050_SSP5RCP85_cropint_pasint.tif')\n",
    "\n",
    "# Obtener la transformación y el sistema de referencia del raster\n",
    "raster_transform = src.transform\n",
    "raster_crs = src.crs\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Reproyectar la malla al sistema de coordenadas del raster si es necesario\n",
    "if Grid_50.crs != raster_crs:\n",
    "    Grid_50 = Grid_50.to_crs(raster_crs)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# PROCESAMIENTO DATOS DE HABITAT"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Con los datos descargados de GLOBIO para los usos del suelo del año 2010 (periodo coincidente con los datos de temperatura) se procede a extraer la información asociada a cada categoría de hábitat definida por GLOBIO que sea de interés para el modelo de distribución de especies. Para el SDM se ha establecido 4 tipologías de habitas, las cuales están constituidas por distintas categorías de GLOBIO.\n",
    "\n",
    "- Habitat 1 - Zonas acuáticas\n",
    "- Habitat 2 - Vegetación Arbórea Perenne\n",
    "- Habitat 3 - Vegetación Arbórea Caduca\n",
    "- Habitat 4 - Matorrales y Pastizales\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Habitat 1 - Zonas acuáticas\n",
    "\n",
    "El habitat 1, corresponderá con ambientes acuaticos, por lo que se incluyen las siguiente tipologias disponible en GLOBIO:\n",
    " - 160: Cubierta arbórea inundada con agua dulce o salobre\n",
    " - 170: Cubierta arbórea inundada con agua salina\n",
    " - 180: Cobertura de arbustos o herbácea inundada con agua dulce/salina/salobre\n",
    " - 210: Cuerpos de agua"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "def extraccion_calculo_area(code):\n",
    "    # Crear una columna para almacenar el área del código específico por celda\n",
    "    column_name = f'area_{code}'\n",
    "    Grid_50[column_name] = 0.0\n",
    "    \n",
    "    # Iterar sobre cada celda en la malla\n",
    "    for idx, cell in Grid_50.iterrows():\n",
    "        # Crear una máscara para la celda actual\n",
    "        geom = [cell.geometry]\n",
    "        out_image, out_transform = mask(src, geom, crop=True, filled=False)\n",
    "        \n",
    "        # Extraer el código específico\n",
    "        code_mask = (out_image == code).astype(np.uint8)\n",
    "        \n",
    "        # Calcular el área del código específico en la celda actual (en m²)\n",
    "        area_code = np.sum(code_mask) * (raster_transform[0] ** 2)  # Área en m²\n",
    "        \n",
    "        # Almacenar el resultado en la celda correspondiente\n",
    "        Grid_50.at[idx, column_name] = area_code"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "extraccion_calculo_area(160)\n",
    "extraccion_calculo_area(170)\n",
    "extraccion_calculo_area(180)\n",
    "extraccion_calculo_area(210)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "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>ET_ID</th>\n",
       "      <th>COUNTRY</th>\n",
       "      <th>geometry</th>\n",
       "      <th>area_160</th>\n",
       "      <th>area_170</th>\n",
       "      <th>area_180</th>\n",
       "      <th>area_210</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.087106</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>Malta</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.224738</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.207083</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000015</td>\n",
       "      <td>0.113789</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.123380</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3101</th>\n",
       "      <td>9369</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((26.92771 70.58708, 25.63446 70.69412...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.004144</td>\n",
       "      <td>0.441667</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3102</th>\n",
       "      <td>9370</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.20775 70.47181, 26.92771 70.58708...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.010941</td>\n",
       "      <td>0.180910</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3103</th>\n",
       "      <td>9371</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((29.47384 70.34844, 28.20775 70.47181...</td>\n",
       "      <td>0.000008</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.006057</td>\n",
       "      <td>0.130077</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3104</th>\n",
       "      <td>9372</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((30.72534 70.21713, 29.47384 70.34844...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.003503</td>\n",
       "      <td>0.214846</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3105</th>\n",
       "      <td>9485</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.60702 70.90647, 27.30206 71.02446...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.002099</td>\n",
       "      <td>0.528511</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3106 rows × 7 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      ET_ID COUNTRY                                           geometry  \\\n",
       "0       249   Spain  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...   \n",
       "1       280   Malta  POLYGON ((14.68659 35.53849, 14.14122 35.56449...   \n",
       "2       363   Spain  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...   \n",
       "3       412  Greece  POLYGON ((23.95653 35.05738, 23.41935 35.13773...   \n",
       "4       413  Greece  POLYGON ((24.49271 34.97391, 23.95653 35.05738...   \n",
       "...     ...     ...                                                ...   \n",
       "3101   9369  Norway  POLYGON ((26.92771 70.58708, 25.63446 70.69412...   \n",
       "3102   9370  Norway  POLYGON ((28.20775 70.47181, 26.92771 70.58708...   \n",
       "3103   9371  Norway  POLYGON ((29.47384 70.34844, 28.20775 70.47181...   \n",
       "3104   9372  Norway  POLYGON ((30.72534 70.21713, 29.47384 70.34844...   \n",
       "3105   9485  Norway  POLYGON ((28.60702 70.90647, 27.30206 71.02446...   \n",
       "\n",
       "      area_160  area_170  area_180  area_210  \n",
       "0     0.000000       0.0  0.000000  0.087106  \n",
       "1     0.000000       0.0  0.000000  0.224738  \n",
       "2     0.000000       0.0  0.000000  0.207083  \n",
       "3     0.000000       0.0  0.000015  0.113789  \n",
       "4     0.000000       0.0  0.000000  0.123380  \n",
       "...        ...       ...       ...       ...  \n",
       "3101  0.000000       0.0  0.004144  0.441667  \n",
       "3102  0.000000       0.0  0.010941  0.180910  \n",
       "3103  0.000008       0.0  0.006057  0.130077  \n",
       "3104  0.000000       0.0  0.003503  0.214846  \n",
       "3105  0.000000       0.0  0.002099  0.528511  \n",
       "\n",
       "[3106 rows x 7 columns]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Grid_50"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Seleccionar las columnas de especies (las que comienzan con 'sp_')\n",
    "area_columns = [col for col in Grid_50.columns if col.startswith('area_')]\n",
    "#Realizar el sumatorio de las distintas columnas con el area de cada tipología\n",
    "Grid_50['Area_H_1'] = Grid_50[area_columns].sum(axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "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>ET_ID</th>\n",
       "      <th>COUNTRY</th>\n",
       "      <th>geometry</th>\n",
       "      <th>area_160</th>\n",
       "      <th>area_170</th>\n",
       "      <th>area_180</th>\n",
       "      <th>area_210</th>\n",
       "      <th>Area_H_1</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.087106</td>\n",
       "      <td>0.087106</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>Malta</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.224738</td>\n",
       "      <td>0.224738</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.207083</td>\n",
       "      <td>0.207083</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000015</td>\n",
       "      <td>0.113789</td>\n",
       "      <td>0.113804</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.123380</td>\n",
       "      <td>0.123380</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   ET_ID COUNTRY                                           geometry  area_160  \\\n",
       "0    249   Spain  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...       0.0   \n",
       "1    280   Malta  POLYGON ((14.68659 35.53849, 14.14122 35.56449...       0.0   \n",
       "2    363   Spain  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...       0.0   \n",
       "3    412  Greece  POLYGON ((23.95653 35.05738, 23.41935 35.13773...       0.0   \n",
       "4    413  Greece  POLYGON ((24.49271 34.97391, 23.95653 35.05738...       0.0   \n",
       "\n",
       "   area_170  area_180  area_210  Area_H_1  \n",
       "0       0.0  0.000000  0.087106  0.087106  \n",
       "1       0.0  0.000000  0.224738  0.224738  \n",
       "2       0.0  0.000000  0.207083  0.207083  \n",
       "3       0.0  0.000015  0.113789  0.113804  \n",
       "4       0.0  0.000000  0.123380  0.123380  "
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Grid_50.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Guardar el resultado en un nuevo shapefile\n",
    "#Grid_50.to_file('Data/GLOBIO/RCP8_5/Grid_H_1.shp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['ET_ID', 'COUNTRY', 'geometry', 'area_160', 'area_170', 'area_180',\n",
       "       'area_210', 'Area_H_1'],\n",
       "      dtype='object')"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Grid_50.columns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "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>ET_ID</th>\n",
       "      <th>COUNTRY</th>\n",
       "      <th>geometry</th>\n",
       "      <th>Area_H_1</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>0.087106</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>Malta</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>0.224738</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>0.207083</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "      <td>0.113804</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "      <td>0.123380</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3101</th>\n",
       "      <td>9369</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((26.92771 70.58708, 25.63446 70.69412...</td>\n",
       "      <td>0.445810</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3102</th>\n",
       "      <td>9370</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.20775 70.47181, 26.92771 70.58708...</td>\n",
       "      <td>0.191852</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3103</th>\n",
       "      <td>9371</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((29.47384 70.34844, 28.20775 70.47181...</td>\n",
       "      <td>0.136142</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3104</th>\n",
       "      <td>9372</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((30.72534 70.21713, 29.47384 70.34844...</td>\n",
       "      <td>0.218349</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3105</th>\n",
       "      <td>9485</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.60702 70.90647, 27.30206 71.02446...</td>\n",
       "      <td>0.530610</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3106 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      ET_ID COUNTRY                                           geometry  \\\n",
       "0       249   Spain  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...   \n",
       "1       280   Malta  POLYGON ((14.68659 35.53849, 14.14122 35.56449...   \n",
       "2       363   Spain  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...   \n",
       "3       412  Greece  POLYGON ((23.95653 35.05738, 23.41935 35.13773...   \n",
       "4       413  Greece  POLYGON ((24.49271 34.97391, 23.95653 35.05738...   \n",
       "...     ...     ...                                                ...   \n",
       "3101   9369  Norway  POLYGON ((26.92771 70.58708, 25.63446 70.69412...   \n",
       "3102   9370  Norway  POLYGON ((28.20775 70.47181, 26.92771 70.58708...   \n",
       "3103   9371  Norway  POLYGON ((29.47384 70.34844, 28.20775 70.47181...   \n",
       "3104   9372  Norway  POLYGON ((30.72534 70.21713, 29.47384 70.34844...   \n",
       "3105   9485  Norway  POLYGON ((28.60702 70.90647, 27.30206 71.02446...   \n",
       "\n",
       "      Area_H_1  \n",
       "0     0.087106  \n",
       "1     0.224738  \n",
       "2     0.207083  \n",
       "3     0.113804  \n",
       "4     0.123380  \n",
       "...        ...  \n",
       "3101  0.445810  \n",
       "3102  0.191852  \n",
       "3103  0.136142  \n",
       "3104  0.218349  \n",
       "3105  0.530610  \n",
       "\n",
       "[3106 rows x 4 columns]"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "columns_select = ['ET_ID', 'COUNTRY', 'geometry',\n",
    "                  'Area_H_1'\n",
    "                  ]\n",
    "\n",
    "Grid_Habitats = Grid_50[columns_select]\n",
    "Grid_Habitats"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Habitat 2 - Vegetación Arbórea Perenne\n",
    "\n",
    "El habitat 2, corresponde con vegetación arbórea perenne, por lo que se incluyen las siguiente tipologias disponible en GLOBIO:\n",
    "\n",
    " - 50: Cubierta arbórea de hoja ancha perenne cerrada a abierta (>15%)\n",
    " - 70: Cubierta arbórea de hoja acicular perenne cerrada a abierta (>15%)\n",
    " - 71: Cubierta arbórea de hoja acicular perenne cerrada (>40%)\n",
    " - 72: Cubierta arbórea de hoja acicular perenne abierta (15-40%)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Cargar el grid nuevamente\n",
    "Grid_50 = gpd.read_file(r'Data/Grid50/Grid_50Km_Europe_WGS84.shp')\n",
    "# Reproyectar la malla al sistema de coordenadas del raster si es necesario\n",
    "if Grid_50.crs != raster_crs:\n",
    "    Grid_50 = Grid_50.to_crs(raster_crs)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "extraccion_calculo_area(50)\n",
    "extraccion_calculo_area(70)\n",
    "extraccion_calculo_area(71)\n",
    "extraccion_calculo_area(72)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "area_columns = [col for col in Grid_50.columns if col.startswith('area_')]\n",
    "\n",
    "Grid_50['Area_H_2'] = Grid_50[area_columns].sum(axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "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>ET_ID</th>\n",
       "      <th>COUNTRY</th>\n",
       "      <th>geometry</th>\n",
       "      <th>area_50</th>\n",
       "      <th>area_70</th>\n",
       "      <th>area_71</th>\n",
       "      <th>area_72</th>\n",
       "      <th>Area_H_2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000069</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000069</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>Malta</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000201</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000201</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.003881</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.003881</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.001119</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.001119</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   ET_ID COUNTRY                                           geometry  area_50  \\\n",
       "0    249   Spain  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...      0.0   \n",
       "1    280   Malta  POLYGON ((14.68659 35.53849, 14.14122 35.56449...      0.0   \n",
       "2    363   Spain  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...      0.0   \n",
       "3    412  Greece  POLYGON ((23.95653 35.05738, 23.41935 35.13773...      0.0   \n",
       "4    413  Greece  POLYGON ((24.49271 34.97391, 23.95653 35.05738...      0.0   \n",
       "\n",
       "    area_70  area_71  area_72  Area_H_2  \n",
       "0  0.000069      0.0      0.0  0.000069  \n",
       "1  0.000000      0.0      0.0  0.000000  \n",
       "2  0.000201      0.0      0.0  0.000201  \n",
       "3  0.003881      0.0      0.0  0.003881  \n",
       "4  0.001119      0.0      0.0  0.001119  "
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Grid_50.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Guardar el resultado en un nuevo shapefile\n",
    "#Grid_50.to_file('Data/GLOBIO/RCP8_5/Grid_H_2.shp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "Grid_50_s = Grid_50[['ET_ID', 'Area_H_2']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>ET_ID</th>\n",
       "      <th>COUNTRY</th>\n",
       "      <th>geometry</th>\n",
       "      <th>Area_H_1</th>\n",
       "      <th>Area_H_2</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>0.087106</td>\n",
       "      <td>0.000069</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>Malta</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>0.224738</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>0.207083</td>\n",
       "      <td>0.000201</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "      <td>0.113804</td>\n",
       "      <td>0.003881</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "      <td>0.123380</td>\n",
       "      <td>0.001119</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3101</th>\n",
       "      <td>9369</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((26.92771 70.58708, 25.63446 70.69412...</td>\n",
       "      <td>0.445810</td>\n",
       "      <td>0.003835</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3102</th>\n",
       "      <td>9370</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.20775 70.47181, 26.92771 70.58708...</td>\n",
       "      <td>0.191852</td>\n",
       "      <td>0.017670</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3103</th>\n",
       "      <td>9371</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((29.47384 70.34844, 28.20775 70.47181...</td>\n",
       "      <td>0.136142</td>\n",
       "      <td>0.007948</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3104</th>\n",
       "      <td>9372</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((30.72534 70.21713, 29.47384 70.34844...</td>\n",
       "      <td>0.218349</td>\n",
       "      <td>0.003403</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3105</th>\n",
       "      <td>9485</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.60702 70.90647, 27.30206 71.02446...</td>\n",
       "      <td>0.530610</td>\n",
       "      <td>0.002461</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3106 rows × 5 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      ET_ID COUNTRY                                           geometry  \\\n",
       "0       249   Spain  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...   \n",
       "1       280   Malta  POLYGON ((14.68659 35.53849, 14.14122 35.56449...   \n",
       "2       363   Spain  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...   \n",
       "3       412  Greece  POLYGON ((23.95653 35.05738, 23.41935 35.13773...   \n",
       "4       413  Greece  POLYGON ((24.49271 34.97391, 23.95653 35.05738...   \n",
       "...     ...     ...                                                ...   \n",
       "3101   9369  Norway  POLYGON ((26.92771 70.58708, 25.63446 70.69412...   \n",
       "3102   9370  Norway  POLYGON ((28.20775 70.47181, 26.92771 70.58708...   \n",
       "3103   9371  Norway  POLYGON ((29.47384 70.34844, 28.20775 70.47181...   \n",
       "3104   9372  Norway  POLYGON ((30.72534 70.21713, 29.47384 70.34844...   \n",
       "3105   9485  Norway  POLYGON ((28.60702 70.90647, 27.30206 71.02446...   \n",
       "\n",
       "      Area_H_1  Area_H_2  \n",
       "0     0.087106  0.000069  \n",
       "1     0.224738  0.000000  \n",
       "2     0.207083  0.000201  \n",
       "3     0.113804  0.003881  \n",
       "4     0.123380  0.001119  \n",
       "...        ...       ...  \n",
       "3101  0.445810  0.003835  \n",
       "3102  0.191852  0.017670  \n",
       "3103  0.136142  0.007948  \n",
       "3104  0.218349  0.003403  \n",
       "3105  0.530610  0.002461  \n",
       "\n",
       "[3106 rows x 5 columns]"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#  Unir los DataFrames usando el campo ET_ID como clave\n",
    "Grid_Habitats = pd.merge(Grid_Habitats, Grid_50_s, on='ET_ID', how='left')\n",
    "Grid_Habitats"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Habitat 3 - Vegetación Arbórea Caduca\n",
    "\n",
    "El habitat 3, corresponde con vegetación arbórea caduca, por lo que se incluyen las siguiente tipologias disponible en GLOBIO:\n",
    "\n",
    " - 60: Cubierta arbórea de hoja ancha caduca cerrada a abierta (>15%)\n",
    " - 61: Cubierta arbórea de hoja ancha caduca cerrada (>40%)\n",
    " - 62: Cubierta arbórea de hoja ancha caduca cerrada (>40%)\n",
    " - 80: Cubierta arbórea de hoja acicular caduca cerrada a abierta (>15%)\n",
    " - 81: Cubierta arbórea de hoja acicular caduca cerrada (>40%)\n",
    " - 82: Cubierta arbórea de hoja acicular caduca abierta (15-40%)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>ET_ID</th>\n",
       "      <th>COUNTRY</th>\n",
       "      <th>geometry</th>\n",
       "      <th>area_60</th>\n",
       "      <th>area_61</th>\n",
       "      <th>area_62</th>\n",
       "      <th>area_80</th>\n",
       "      <th>area_81</th>\n",
       "      <th>area_82</th>\n",
       "      <th>Area_H_3</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>0.000293</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000293</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>Malta</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>0.000046</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000046</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "      <td>0.000139</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000139</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   ET_ID COUNTRY                                           geometry   area_60  \\\n",
       "0    249   Spain  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...  0.000293   \n",
       "1    280   Malta  POLYGON ((14.68659 35.53849, 14.14122 35.56449...  0.000000   \n",
       "2    363   Spain  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...  0.000046   \n",
       "3    412  Greece  POLYGON ((23.95653 35.05738, 23.41935 35.13773...  0.000139   \n",
       "4    413  Greece  POLYGON ((24.49271 34.97391, 23.95653 35.05738...  0.000000   \n",
       "\n",
       "   area_61  area_62  area_80  area_81  area_82  Area_H_3  \n",
       "0      0.0      0.0      0.0      0.0      0.0  0.000293  \n",
       "1      0.0      0.0      0.0      0.0      0.0  0.000000  \n",
       "2      0.0      0.0      0.0      0.0      0.0  0.000046  \n",
       "3      0.0      0.0      0.0      0.0      0.0  0.000139  \n",
       "4      0.0      0.0      0.0      0.0      0.0  0.000000  "
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#Cargar el grid nuevamente\n",
    "Grid_50 = gpd.read_file(r'Data/Grid50/Grid_50Km_Europe_WGS84.shp')\n",
    "# Reproyectar la malla al sistema de coordenadas del raster si es necesario\n",
    "if Grid_50.crs != raster_crs:\n",
    "    Grid_50 = Grid_50.to_crs(raster_crs)\n",
    "\n",
    "extraccion_calculo_area(60)\n",
    "extraccion_calculo_area(61)\n",
    "extraccion_calculo_area(62)\n",
    "extraccion_calculo_area(80)\n",
    "extraccion_calculo_area(81)\n",
    "extraccion_calculo_area(82)\n",
    "\n",
    "area_columns = [col for col in Grid_50.columns if col.startswith('area_')]\n",
    "Grid_50['Area_H_3'] = Grid_50[area_columns].sum(axis=1)\n",
    "\n",
    "Grid_50.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Guardar el resultado en un nuevo shapefile\n",
    "#Grid_50.to_file('Data/GLOBIO/CP8_5/Grid_H_3.shp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "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>ET_ID</th>\n",
       "      <th>COUNTRY</th>\n",
       "      <th>geometry</th>\n",
       "      <th>Area_H_1</th>\n",
       "      <th>Area_H_2</th>\n",
       "      <th>Area_H_3</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>0.087106</td>\n",
       "      <td>0.000069</td>\n",
       "      <td>0.000293</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>Malta</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>0.224738</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>0.207083</td>\n",
       "      <td>0.000201</td>\n",
       "      <td>0.000046</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "      <td>0.113804</td>\n",
       "      <td>0.003881</td>\n",
       "      <td>0.000139</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "      <td>0.123380</td>\n",
       "      <td>0.001119</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3101</th>\n",
       "      <td>9369</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((26.92771 70.58708, 25.63446 70.69412...</td>\n",
       "      <td>0.445810</td>\n",
       "      <td>0.003835</td>\n",
       "      <td>0.004298</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3102</th>\n",
       "      <td>9370</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.20775 70.47181, 26.92771 70.58708...</td>\n",
       "      <td>0.191852</td>\n",
       "      <td>0.017670</td>\n",
       "      <td>0.016150</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3103</th>\n",
       "      <td>9371</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((29.47384 70.34844, 28.20775 70.47181...</td>\n",
       "      <td>0.136142</td>\n",
       "      <td>0.007948</td>\n",
       "      <td>0.010316</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3104</th>\n",
       "      <td>9372</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((30.72534 70.21713, 29.47384 70.34844...</td>\n",
       "      <td>0.218349</td>\n",
       "      <td>0.003403</td>\n",
       "      <td>0.003356</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3105</th>\n",
       "      <td>9485</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.60702 70.90647, 27.30206 71.02446...</td>\n",
       "      <td>0.530610</td>\n",
       "      <td>0.002461</td>\n",
       "      <td>0.000571</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3106 rows × 6 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      ET_ID COUNTRY                                           geometry  \\\n",
       "0       249   Spain  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...   \n",
       "1       280   Malta  POLYGON ((14.68659 35.53849, 14.14122 35.56449...   \n",
       "2       363   Spain  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...   \n",
       "3       412  Greece  POLYGON ((23.95653 35.05738, 23.41935 35.13773...   \n",
       "4       413  Greece  POLYGON ((24.49271 34.97391, 23.95653 35.05738...   \n",
       "...     ...     ...                                                ...   \n",
       "3101   9369  Norway  POLYGON ((26.92771 70.58708, 25.63446 70.69412...   \n",
       "3102   9370  Norway  POLYGON ((28.20775 70.47181, 26.92771 70.58708...   \n",
       "3103   9371  Norway  POLYGON ((29.47384 70.34844, 28.20775 70.47181...   \n",
       "3104   9372  Norway  POLYGON ((30.72534 70.21713, 29.47384 70.34844...   \n",
       "3105   9485  Norway  POLYGON ((28.60702 70.90647, 27.30206 71.02446...   \n",
       "\n",
       "      Area_H_1  Area_H_2  Area_H_3  \n",
       "0     0.087106  0.000069  0.000293  \n",
       "1     0.224738  0.000000  0.000000  \n",
       "2     0.207083  0.000201  0.000046  \n",
       "3     0.113804  0.003881  0.000139  \n",
       "4     0.123380  0.001119  0.000000  \n",
       "...        ...       ...       ...  \n",
       "3101  0.445810  0.003835  0.004298  \n",
       "3102  0.191852  0.017670  0.016150  \n",
       "3103  0.136142  0.007948  0.010316  \n",
       "3104  0.218349  0.003403  0.003356  \n",
       "3105  0.530610  0.002461  0.000571  \n",
       "\n",
       "[3106 rows x 6 columns]"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Grid_50_s = Grid_50[['ET_ID', 'Area_H_3']]\n",
    "Grid_Habitats = pd.merge(Grid_Habitats, Grid_50_s, on='ET_ID', how='left')\n",
    "Grid_Habitats"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Habitat 4 - Matorrales y Pastizales\n",
    "\n",
    "El habitat 4, corresponde con matorrales y pastizales, por lo que se incluyen las siguiente tipologias disponible en GLOBIO:\n",
    "\n",
    " - 2: Tierras de cultivo\n",
    " - 3: Pastizales en general\n",
    " - 4: Pastizales utilizados para la cría de ganado y la agricultura extensiva en zonas de pastoreo\n",
    " - 120: Matorrales\n",
    " - 121: Matorrales perennes\n",
    " - 122: Matorrales caducifolios\n",
    " - 130: Pastizales\n",
    " - 150: Vegetación dispersa (árboles, arbustos, cobertura herbácea) (<15%)\n",
    " - 152: Arbustos dispersos (<15%)\n",
    " - 153: Cobertura herbácea dispersa (<15%)\n",
    " - 110: Mosaico de cobertura herbácea (>50%) / árboles y arbustos (<50%)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "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>ET_ID</th>\n",
       "      <th>COUNTRY</th>\n",
       "      <th>geometry</th>\n",
       "      <th>area_2</th>\n",
       "      <th>area_3</th>\n",
       "      <th>area_4</th>\n",
       "      <th>area_120</th>\n",
       "      <th>area_121</th>\n",
       "      <th>area_122</th>\n",
       "      <th>area_130</th>\n",
       "      <th>area_150</th>\n",
       "      <th>area_152</th>\n",
       "      <th>area_153</th>\n",
       "      <th>area_110</th>\n",
       "      <th>Area_H_4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.038179</td>\n",
       "      <td>0.027685</td>\n",
       "      <td>0.002870</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000556</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000162</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.069452</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>Malta</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000039</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000008</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000054</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000100</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.012623</td>\n",
       "      <td>0.008225</td>\n",
       "      <td>0.000856</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000054</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000062</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.021821</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.047346</td>\n",
       "      <td>0.006281</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.004660</td>\n",
       "      <td>0.000571</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.058858</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.039931</td>\n",
       "      <td>0.003789</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000085</td>\n",
       "      <td>0.014514</td>\n",
       "      <td>0.001073</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000201</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.059591</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   ET_ID COUNTRY                                           geometry  area_2  \\\n",
       "0    249   Spain  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...     0.0   \n",
       "1    280   Malta  POLYGON ((14.68659 35.53849, 14.14122 35.56449...     0.0   \n",
       "2    363   Spain  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...     0.0   \n",
       "3    412  Greece  POLYGON ((23.95653 35.05738, 23.41935 35.13773...     0.0   \n",
       "4    413  Greece  POLYGON ((24.49271 34.97391, 23.95653 35.05738...     0.0   \n",
       "\n",
       "     area_3    area_4  area_120  area_121  area_122  area_130  area_150  \\\n",
       "0  0.038179  0.027685  0.002870       0.0  0.000000  0.000000  0.000556   \n",
       "1  0.000000  0.000000  0.000039       0.0  0.000008  0.000000  0.000054   \n",
       "2  0.012623  0.008225  0.000856       0.0  0.000000  0.000000  0.000054   \n",
       "3  0.000000  0.047346  0.006281       0.0  0.000000  0.004660  0.000571   \n",
       "4  0.000000  0.039931  0.003789       0.0  0.000085  0.014514  0.001073   \n",
       "\n",
       "   area_152  area_153  area_110  Area_H_4  \n",
       "0       0.0  0.000162       0.0  0.069452  \n",
       "1       0.0  0.000000       0.0  0.000100  \n",
       "2       0.0  0.000062       0.0  0.021821  \n",
       "3       0.0  0.000000       0.0  0.058858  \n",
       "4       0.0  0.000201       0.0  0.059591  "
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#Cargar el grid nuevamente\n",
    "Grid_50 = gpd.read_file(r'Data/Grid50/Grid_50Km_Europe_WGS84.shp')\n",
    "# Reproyectar la malla al sistema de coordenadas del raster si es necesario\n",
    "if Grid_50.crs != raster_crs:\n",
    "    Grid_50 = Grid_50.to_crs(raster_crs)\n",
    "\n",
    "extraccion_calculo_area(2)\n",
    "extraccion_calculo_area(3)\n",
    "extraccion_calculo_area(4)\n",
    "extraccion_calculo_area(120)\n",
    "extraccion_calculo_area(121)\n",
    "extraccion_calculo_area(122)\n",
    "extraccion_calculo_area(130)\n",
    "extraccion_calculo_area(150)\n",
    "extraccion_calculo_area(152)\n",
    "extraccion_calculo_area(153)\n",
    "extraccion_calculo_area(110)\n",
    "\n",
    "area_columns = [col for col in Grid_50.columns if col.startswith('area_')]\n",
    "Grid_50['Area_H_4'] = Grid_50[area_columns].sum(axis=1)\n",
    "Grid_50.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Guardar el resultado en un nuevo shapefile\n",
    "#Grid_50.to_file('Data/GLOBIO/RCP8_5/Grid_H_4.shp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "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>ET_ID</th>\n",
       "      <th>COUNTRY</th>\n",
       "      <th>geometry</th>\n",
       "      <th>Area_H_1</th>\n",
       "      <th>Area_H_2</th>\n",
       "      <th>Area_H_3</th>\n",
       "      <th>Area_H_4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>0.087106</td>\n",
       "      <td>0.000069</td>\n",
       "      <td>0.000293</td>\n",
       "      <td>0.069452</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>Malta</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>0.224738</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000100</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>0.207083</td>\n",
       "      <td>0.000201</td>\n",
       "      <td>0.000046</td>\n",
       "      <td>0.021821</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "      <td>0.113804</td>\n",
       "      <td>0.003881</td>\n",
       "      <td>0.000139</td>\n",
       "      <td>0.058858</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "      <td>0.123380</td>\n",
       "      <td>0.001119</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.059591</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3101</th>\n",
       "      <td>9369</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((26.92771 70.58708, 25.63446 70.69412...</td>\n",
       "      <td>0.445810</td>\n",
       "      <td>0.003835</td>\n",
       "      <td>0.004298</td>\n",
       "      <td>0.114128</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3102</th>\n",
       "      <td>9370</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.20775 70.47181, 26.92771 70.58708...</td>\n",
       "      <td>0.191852</td>\n",
       "      <td>0.017670</td>\n",
       "      <td>0.016150</td>\n",
       "      <td>0.271250</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3103</th>\n",
       "      <td>9371</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((29.47384 70.34844, 28.20775 70.47181...</td>\n",
       "      <td>0.136142</td>\n",
       "      <td>0.007948</td>\n",
       "      <td>0.010316</td>\n",
       "      <td>0.285023</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3104</th>\n",
       "      <td>9372</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((30.72534 70.21713, 29.47384 70.34844...</td>\n",
       "      <td>0.218349</td>\n",
       "      <td>0.003403</td>\n",
       "      <td>0.003356</td>\n",
       "      <td>0.235787</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3105</th>\n",
       "      <td>9485</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.60702 70.90647, 27.30206 71.02446...</td>\n",
       "      <td>0.530610</td>\n",
       "      <td>0.002461</td>\n",
       "      <td>0.000571</td>\n",
       "      <td>0.052685</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3106 rows × 7 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      ET_ID COUNTRY                                           geometry  \\\n",
       "0       249   Spain  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...   \n",
       "1       280   Malta  POLYGON ((14.68659 35.53849, 14.14122 35.56449...   \n",
       "2       363   Spain  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...   \n",
       "3       412  Greece  POLYGON ((23.95653 35.05738, 23.41935 35.13773...   \n",
       "4       413  Greece  POLYGON ((24.49271 34.97391, 23.95653 35.05738...   \n",
       "...     ...     ...                                                ...   \n",
       "3101   9369  Norway  POLYGON ((26.92771 70.58708, 25.63446 70.69412...   \n",
       "3102   9370  Norway  POLYGON ((28.20775 70.47181, 26.92771 70.58708...   \n",
       "3103   9371  Norway  POLYGON ((29.47384 70.34844, 28.20775 70.47181...   \n",
       "3104   9372  Norway  POLYGON ((30.72534 70.21713, 29.47384 70.34844...   \n",
       "3105   9485  Norway  POLYGON ((28.60702 70.90647, 27.30206 71.02446...   \n",
       "\n",
       "      Area_H_1  Area_H_2  Area_H_3  Area_H_4  \n",
       "0     0.087106  0.000069  0.000293  0.069452  \n",
       "1     0.224738  0.000000  0.000000  0.000100  \n",
       "2     0.207083  0.000201  0.000046  0.021821  \n",
       "3     0.113804  0.003881  0.000139  0.058858  \n",
       "4     0.123380  0.001119  0.000000  0.059591  \n",
       "...        ...       ...       ...       ...  \n",
       "3101  0.445810  0.003835  0.004298  0.114128  \n",
       "3102  0.191852  0.017670  0.016150  0.271250  \n",
       "3103  0.136142  0.007948  0.010316  0.285023  \n",
       "3104  0.218349  0.003403  0.003356  0.235787  \n",
       "3105  0.530610  0.002461  0.000571  0.052685  \n",
       "\n",
       "[3106 rows x 7 columns]"
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Guardar el resultado en un nuevo shapefile\n",
    "Grid_50_s = Grid_50[['ET_ID', 'Area_H_4']]\n",
    "\n",
    "#  Unir los DataFrames usando el campo ET_ID como clave\n",
    "Grid_Habitats = pd.merge(Grid_Habitats, Grid_50_s, on='ET_ID', how='left')\n",
    "Grid_Habitats"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Guardar el resultado en un nuevo shapefile\n",
    "#Grid_Habitats.to_file('Data/GLOBIO/RCP8_5/Grid_Habitats.shp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Cerrar el dataset raster\n",
    "src.close()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# % DE ÁREA DISPONIBLE POR CUADRÍCULA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\molic\\AppData\\Local\\Temp\\ipykernel_11008\\850118656.py:2: UserWarning: Geometry is in a geographic CRS. Results from 'area' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.\n",
      "\n",
      "  Grid_Habitats['Area_grid'] = Grid_Habitats.geometry.area #en m2\n"
     ]
    },
    {
     "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>ET_ID</th>\n",
       "      <th>COUNTRY</th>\n",
       "      <th>geometry</th>\n",
       "      <th>Area_H_1</th>\n",
       "      <th>Area_H_2</th>\n",
       "      <th>Area_H_3</th>\n",
       "      <th>Area_H_4</th>\n",
       "      <th>Area_grid</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>0.087106</td>\n",
       "      <td>0.000069</td>\n",
       "      <td>0.000293</td>\n",
       "      <td>0.069452</td>\n",
       "      <td>0.246970</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>Malta</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>0.224738</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000100</td>\n",
       "      <td>0.249205</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>0.207083</td>\n",
       "      <td>0.000201</td>\n",
       "      <td>0.000046</td>\n",
       "      <td>0.021821</td>\n",
       "      <td>0.248088</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "      <td>0.113804</td>\n",
       "      <td>0.003881</td>\n",
       "      <td>0.000139</td>\n",
       "      <td>0.058858</td>\n",
       "      <td>0.247813</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "      <td>0.123380</td>\n",
       "      <td>0.001119</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.059591</td>\n",
       "      <td>0.247566</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3101</th>\n",
       "      <td>9369</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((26.92771 70.58708, 25.63446 70.69412...</td>\n",
       "      <td>0.445810</td>\n",
       "      <td>0.003835</td>\n",
       "      <td>0.004298</td>\n",
       "      <td>0.114128</td>\n",
       "      <td>0.612105</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3102</th>\n",
       "      <td>9370</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.20775 70.47181, 26.92771 70.58708...</td>\n",
       "      <td>0.191852</td>\n",
       "      <td>0.017670</td>\n",
       "      <td>0.016150</td>\n",
       "      <td>0.271250</td>\n",
       "      <td>0.608670</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3103</th>\n",
       "      <td>9371</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((29.47384 70.34844, 28.20775 70.47181...</td>\n",
       "      <td>0.136142</td>\n",
       "      <td>0.007948</td>\n",
       "      <td>0.010316</td>\n",
       "      <td>0.285023</td>\n",
       "      <td>0.605030</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3104</th>\n",
       "      <td>9372</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((30.72534 70.21713, 29.47384 70.34844...</td>\n",
       "      <td>0.218349</td>\n",
       "      <td>0.003403</td>\n",
       "      <td>0.003356</td>\n",
       "      <td>0.235787</td>\n",
       "      <td>0.601196</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3105</th>\n",
       "      <td>9485</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.60702 70.90647, 27.30206 71.02446...</td>\n",
       "      <td>0.530610</td>\n",
       "      <td>0.002461</td>\n",
       "      <td>0.000571</td>\n",
       "      <td>0.052685</td>\n",
       "      <td>0.622198</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3106 rows × 8 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      ET_ID COUNTRY                                           geometry  \\\n",
       "0       249   Spain  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...   \n",
       "1       280   Malta  POLYGON ((14.68659 35.53849, 14.14122 35.56449...   \n",
       "2       363   Spain  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...   \n",
       "3       412  Greece  POLYGON ((23.95653 35.05738, 23.41935 35.13773...   \n",
       "4       413  Greece  POLYGON ((24.49271 34.97391, 23.95653 35.05738...   \n",
       "...     ...     ...                                                ...   \n",
       "3101   9369  Norway  POLYGON ((26.92771 70.58708, 25.63446 70.69412...   \n",
       "3102   9370  Norway  POLYGON ((28.20775 70.47181, 26.92771 70.58708...   \n",
       "3103   9371  Norway  POLYGON ((29.47384 70.34844, 28.20775 70.47181...   \n",
       "3104   9372  Norway  POLYGON ((30.72534 70.21713, 29.47384 70.34844...   \n",
       "3105   9485  Norway  POLYGON ((28.60702 70.90647, 27.30206 71.02446...   \n",
       "\n",
       "      Area_H_1  Area_H_2  Area_H_3  Area_H_4  Area_grid  \n",
       "0     0.087106  0.000069  0.000293  0.069452   0.246970  \n",
       "1     0.224738  0.000000  0.000000  0.000100   0.249205  \n",
       "2     0.207083  0.000201  0.000046  0.021821   0.248088  \n",
       "3     0.113804  0.003881  0.000139  0.058858   0.247813  \n",
       "4     0.123380  0.001119  0.000000  0.059591   0.247566  \n",
       "...        ...       ...       ...       ...        ...  \n",
       "3101  0.445810  0.003835  0.004298  0.114128   0.612105  \n",
       "3102  0.191852  0.017670  0.016150  0.271250   0.608670  \n",
       "3103  0.136142  0.007948  0.010316  0.285023   0.605030  \n",
       "3104  0.218349  0.003403  0.003356  0.235787   0.601196  \n",
       "3105  0.530610  0.002461  0.000571  0.052685   0.622198  \n",
       "\n",
       "[3106 rows x 8 columns]"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Calculo del area para cada cuadrícula en km²\n",
    "Grid_Habitats['Area_grid'] = Grid_Habitats.geometry.area #en m2\n",
    "Grid_Habitats"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "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>ET_ID</th>\n",
       "      <th>COUNTRY</th>\n",
       "      <th>geometry</th>\n",
       "      <th>Area_H_1</th>\n",
       "      <th>Area_H_2</th>\n",
       "      <th>Area_H_3</th>\n",
       "      <th>Area_H_4</th>\n",
       "      <th>Area_grid</th>\n",
       "      <th>%_H_1</th>\n",
       "      <th>%_H_2</th>\n",
       "      <th>%_H_3</th>\n",
       "      <th>%_H_4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>0.087106</td>\n",
       "      <td>0.000069</td>\n",
       "      <td>0.000293</td>\n",
       "      <td>0.069452</td>\n",
       "      <td>0.246970</td>\n",
       "      <td>35.270045</td>\n",
       "      <td>0.028119</td>\n",
       "      <td>0.118723</td>\n",
       "      <td>28.121682</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>Malta</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>0.224738</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000100</td>\n",
       "      <td>0.249205</td>\n",
       "      <td>90.181696</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.040251</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>Spain</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>0.207083</td>\n",
       "      <td>0.000201</td>\n",
       "      <td>0.000046</td>\n",
       "      <td>0.021821</td>\n",
       "      <td>0.248088</td>\n",
       "      <td>83.471585</td>\n",
       "      <td>0.080865</td>\n",
       "      <td>0.018661</td>\n",
       "      <td>8.795650</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "      <td>0.113804</td>\n",
       "      <td>0.003881</td>\n",
       "      <td>0.000139</td>\n",
       "      <td>0.058858</td>\n",
       "      <td>0.247813</td>\n",
       "      <td>45.923425</td>\n",
       "      <td>1.566173</td>\n",
       "      <td>0.056046</td>\n",
       "      <td>23.751026</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>Greece</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "      <td>0.123380</td>\n",
       "      <td>0.001119</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.059591</td>\n",
       "      <td>0.247566</td>\n",
       "      <td>49.837149</td>\n",
       "      <td>0.451932</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>24.070813</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3101</th>\n",
       "      <td>9369</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((26.92771 70.58708, 25.63446 70.69412...</td>\n",
       "      <td>0.445810</td>\n",
       "      <td>0.003835</td>\n",
       "      <td>0.004298</td>\n",
       "      <td>0.114128</td>\n",
       "      <td>0.612105</td>\n",
       "      <td>72.832347</td>\n",
       "      <td>0.626507</td>\n",
       "      <td>0.702141</td>\n",
       "      <td>18.645192</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3102</th>\n",
       "      <td>9370</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.20775 70.47181, 26.92771 70.58708...</td>\n",
       "      <td>0.191852</td>\n",
       "      <td>0.017670</td>\n",
       "      <td>0.016150</td>\n",
       "      <td>0.271250</td>\n",
       "      <td>0.608670</td>\n",
       "      <td>31.519831</td>\n",
       "      <td>2.903009</td>\n",
       "      <td>2.653274</td>\n",
       "      <td>44.564356</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3103</th>\n",
       "      <td>9371</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((29.47384 70.34844, 28.20775 70.47181...</td>\n",
       "      <td>0.136142</td>\n",
       "      <td>0.007948</td>\n",
       "      <td>0.010316</td>\n",
       "      <td>0.285023</td>\n",
       "      <td>0.605030</td>\n",
       "      <td>22.501694</td>\n",
       "      <td>1.313577</td>\n",
       "      <td>1.705099</td>\n",
       "      <td>47.108936</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3104</th>\n",
       "      <td>9372</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((30.72534 70.21713, 29.47384 70.34844...</td>\n",
       "      <td>0.218349</td>\n",
       "      <td>0.003403</td>\n",
       "      <td>0.003356</td>\n",
       "      <td>0.235787</td>\n",
       "      <td>0.601196</td>\n",
       "      <td>36.319039</td>\n",
       "      <td>0.566001</td>\n",
       "      <td>0.558300</td>\n",
       "      <td>39.219634</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3105</th>\n",
       "      <td>9485</td>\n",
       "      <td>Norway</td>\n",
       "      <td>POLYGON ((28.60702 70.90647, 27.30206 71.02446...</td>\n",
       "      <td>0.530610</td>\n",
       "      <td>0.002461</td>\n",
       "      <td>0.000571</td>\n",
       "      <td>0.052685</td>\n",
       "      <td>0.622198</td>\n",
       "      <td>85.279899</td>\n",
       "      <td>0.395601</td>\n",
       "      <td>0.091769</td>\n",
       "      <td>8.467596</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3106 rows × 12 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      ET_ID COUNTRY                                           geometry  \\\n",
       "0       249   Spain  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...   \n",
       "1       280   Malta  POLYGON ((14.68659 35.53849, 14.14122 35.56449...   \n",
       "2       363   Spain  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...   \n",
       "3       412  Greece  POLYGON ((23.95653 35.05738, 23.41935 35.13773...   \n",
       "4       413  Greece  POLYGON ((24.49271 34.97391, 23.95653 35.05738...   \n",
       "...     ...     ...                                                ...   \n",
       "3101   9369  Norway  POLYGON ((26.92771 70.58708, 25.63446 70.69412...   \n",
       "3102   9370  Norway  POLYGON ((28.20775 70.47181, 26.92771 70.58708...   \n",
       "3103   9371  Norway  POLYGON ((29.47384 70.34844, 28.20775 70.47181...   \n",
       "3104   9372  Norway  POLYGON ((30.72534 70.21713, 29.47384 70.34844...   \n",
       "3105   9485  Norway  POLYGON ((28.60702 70.90647, 27.30206 71.02446...   \n",
       "\n",
       "      Area_H_1  Area_H_2  Area_H_3  Area_H_4  Area_grid      %_H_1     %_H_2  \\\n",
       "0     0.087106  0.000069  0.000293  0.069452   0.246970  35.270045  0.028119   \n",
       "1     0.224738  0.000000  0.000000  0.000100   0.249205  90.181696  0.000000   \n",
       "2     0.207083  0.000201  0.000046  0.021821   0.248088  83.471585  0.080865   \n",
       "3     0.113804  0.003881  0.000139  0.058858   0.247813  45.923425  1.566173   \n",
       "4     0.123380  0.001119  0.000000  0.059591   0.247566  49.837149  0.451932   \n",
       "...        ...       ...       ...       ...        ...        ...       ...   \n",
       "3101  0.445810  0.003835  0.004298  0.114128   0.612105  72.832347  0.626507   \n",
       "3102  0.191852  0.017670  0.016150  0.271250   0.608670  31.519831  2.903009   \n",
       "3103  0.136142  0.007948  0.010316  0.285023   0.605030  22.501694  1.313577   \n",
       "3104  0.218349  0.003403  0.003356  0.235787   0.601196  36.319039  0.566001   \n",
       "3105  0.530610  0.002461  0.000571  0.052685   0.622198  85.279899  0.395601   \n",
       "\n",
       "         %_H_3      %_H_4  \n",
       "0     0.118723  28.121682  \n",
       "1     0.000000   0.040251  \n",
       "2     0.018661   8.795650  \n",
       "3     0.056046  23.751026  \n",
       "4     0.000000  24.070813  \n",
       "...        ...        ...  \n",
       "3101  0.702141  18.645192  \n",
       "3102  2.653274  44.564356  \n",
       "3103  1.705099  47.108936  \n",
       "3104  0.558300  39.219634  \n",
       "3105  0.091769   8.467596  \n",
       "\n",
       "[3106 rows x 12 columns]"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Calculo porcentaje de habitat\n",
    "Grid_Habitats['%_H_1'] = (Grid_Habitats['Area_H_1']*100)/Grid_Habitats['Area_grid']\n",
    "Grid_Habitats['%_H_2'] = (Grid_Habitats['Area_H_2']*100)/Grid_Habitats['Area_grid']\n",
    "Grid_Habitats['%_H_3'] = (Grid_Habitats['Area_H_3']*100)/Grid_Habitats['Area_grid']\n",
    "Grid_Habitats['%_H_4'] = (Grid_Habitats['Area_H_4']*100)/Grid_Habitats['Area_grid']\n",
    "\n",
    "Grid_Habitats"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "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>ET_ID</th>\n",
       "      <th>geometry</th>\n",
       "      <th>%_H_1</th>\n",
       "      <th>%_H_2</th>\n",
       "      <th>%_H_3</th>\n",
       "      <th>%_H_4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>35.270045</td>\n",
       "      <td>0.028119</td>\n",
       "      <td>0.118723</td>\n",
       "      <td>28.121682</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>90.181696</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.040251</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>83.471585</td>\n",
       "      <td>0.080865</td>\n",
       "      <td>0.018661</td>\n",
       "      <td>8.795650</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "      <td>45.923425</td>\n",
       "      <td>1.566173</td>\n",
       "      <td>0.056046</td>\n",
       "      <td>23.751026</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "      <td>49.837149</td>\n",
       "      <td>0.451932</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>24.070813</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3101</th>\n",
       "      <td>9369</td>\n",
       "      <td>POLYGON ((26.92771 70.58708, 25.63446 70.69412...</td>\n",
       "      <td>72.832347</td>\n",
       "      <td>0.626507</td>\n",
       "      <td>0.702141</td>\n",
       "      <td>18.645192</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3102</th>\n",
       "      <td>9370</td>\n",
       "      <td>POLYGON ((28.20775 70.47181, 26.92771 70.58708...</td>\n",
       "      <td>31.519831</td>\n",
       "      <td>2.903009</td>\n",
       "      <td>2.653274</td>\n",
       "      <td>44.564356</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3103</th>\n",
       "      <td>9371</td>\n",
       "      <td>POLYGON ((29.47384 70.34844, 28.20775 70.47181...</td>\n",
       "      <td>22.501694</td>\n",
       "      <td>1.313577</td>\n",
       "      <td>1.705099</td>\n",
       "      <td>47.108936</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3104</th>\n",
       "      <td>9372</td>\n",
       "      <td>POLYGON ((30.72534 70.21713, 29.47384 70.34844...</td>\n",
       "      <td>36.319039</td>\n",
       "      <td>0.566001</td>\n",
       "      <td>0.558300</td>\n",
       "      <td>39.219634</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3105</th>\n",
       "      <td>9485</td>\n",
       "      <td>POLYGON ((28.60702 70.90647, 27.30206 71.02446...</td>\n",
       "      <td>85.279899</td>\n",
       "      <td>0.395601</td>\n",
       "      <td>0.091769</td>\n",
       "      <td>8.467596</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3106 rows × 6 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      ET_ID                                           geometry      %_H_1  \\\n",
       "0       249  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...  35.270045   \n",
       "1       280  POLYGON ((14.68659 35.53849, 14.14122 35.56449...  90.181696   \n",
       "2       363  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...  83.471585   \n",
       "3       412  POLYGON ((23.95653 35.05738, 23.41935 35.13773...  45.923425   \n",
       "4       413  POLYGON ((24.49271 34.97391, 23.95653 35.05738...  49.837149   \n",
       "...     ...                                                ...        ...   \n",
       "3101   9369  POLYGON ((26.92771 70.58708, 25.63446 70.69412...  72.832347   \n",
       "3102   9370  POLYGON ((28.20775 70.47181, 26.92771 70.58708...  31.519831   \n",
       "3103   9371  POLYGON ((29.47384 70.34844, 28.20775 70.47181...  22.501694   \n",
       "3104   9372  POLYGON ((30.72534 70.21713, 29.47384 70.34844...  36.319039   \n",
       "3105   9485  POLYGON ((28.60702 70.90647, 27.30206 71.02446...  85.279899   \n",
       "\n",
       "         %_H_2     %_H_3      %_H_4  \n",
       "0     0.028119  0.118723  28.121682  \n",
       "1     0.000000  0.000000   0.040251  \n",
       "2     0.080865  0.018661   8.795650  \n",
       "3     1.566173  0.056046  23.751026  \n",
       "4     0.451932  0.000000  24.070813  \n",
       "...        ...       ...        ...  \n",
       "3101  0.626507  0.702141  18.645192  \n",
       "3102  2.903009  2.653274  44.564356  \n",
       "3103  1.313577  1.705099  47.108936  \n",
       "3104  0.566001  0.558300  39.219634  \n",
       "3105  0.395601  0.091769   8.467596  \n",
       "\n",
       "[3106 rows x 6 columns]"
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Grid_Habitats_Variables = Grid_Habitats[['ET_ID', 'geometry', '%_H_1', '%_H_2', '%_H_3', '%_H_4']]\n",
    "Grid_Habitats_Variables"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Guardar el resultado en un nuevo shapefile\n",
    "#Grid_Habitats_Variables.to_file('Data/GLOBIO/RCP8_5/Grid_Habitats_Variables.shp')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# UNION CON VARIABLES CLIMÁTICAS"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Se tienen en cuenta solo las variables seleccionadas con el Jupyter Notebook 06_Analisis_variables_independientes.ipynb"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "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>ET_ID</th>\n",
       "      <th>%_H_1</th>\n",
       "      <th>%_H_2</th>\n",
       "      <th>%_H_3</th>\n",
       "      <th>%_H_4</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>35.270045</td>\n",
       "      <td>0.028119</td>\n",
       "      <td>0.118723</td>\n",
       "      <td>28.121682</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>90.181696</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.040251</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>83.471585</td>\n",
       "      <td>0.080865</td>\n",
       "      <td>0.018661</td>\n",
       "      <td>8.795650</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>412</td>\n",
       "      <td>45.923425</td>\n",
       "      <td>1.566173</td>\n",
       "      <td>0.056046</td>\n",
       "      <td>23.751026</td>\n",
       "      <td>POLYGON ((23.95653 35.05738, 23.41935 35.13773...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>413</td>\n",
       "      <td>49.837149</td>\n",
       "      <td>0.451932</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>24.070813</td>\n",
       "      <td>POLYGON ((24.49271 34.97391, 23.95653 35.05738...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3101</th>\n",
       "      <td>9369</td>\n",
       "      <td>72.832347</td>\n",
       "      <td>0.626507</td>\n",
       "      <td>0.702141</td>\n",
       "      <td>18.645192</td>\n",
       "      <td>POLYGON ((26.92771 70.58708, 25.63446 70.69412...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3102</th>\n",
       "      <td>9370</td>\n",
       "      <td>31.519831</td>\n",
       "      <td>2.903009</td>\n",
       "      <td>2.653274</td>\n",
       "      <td>44.564356</td>\n",
       "      <td>POLYGON ((28.20775 70.47181, 26.92771 70.58708...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3103</th>\n",
       "      <td>9371</td>\n",
       "      <td>22.501694</td>\n",
       "      <td>1.313577</td>\n",
       "      <td>1.705099</td>\n",
       "      <td>47.108936</td>\n",
       "      <td>POLYGON ((29.47384 70.34844, 28.20775 70.47181...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3104</th>\n",
       "      <td>9372</td>\n",
       "      <td>36.319039</td>\n",
       "      <td>0.566001</td>\n",
       "      <td>0.558300</td>\n",
       "      <td>39.219634</td>\n",
       "      <td>POLYGON ((30.72534 70.21713, 29.47384 70.34844...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3105</th>\n",
       "      <td>9485</td>\n",
       "      <td>85.279899</td>\n",
       "      <td>0.395601</td>\n",
       "      <td>0.091769</td>\n",
       "      <td>8.467596</td>\n",
       "      <td>POLYGON ((28.60702 70.90647, 27.30206 71.02446...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3106 rows × 6 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      ET_ID      %_H_1     %_H_2     %_H_3      %_H_4  \\\n",
       "0       249  35.270045  0.028119  0.118723  28.121682   \n",
       "1       280  90.181696  0.000000  0.000000   0.040251   \n",
       "2       363  83.471585  0.080865  0.018661   8.795650   \n",
       "3       412  45.923425  1.566173  0.056046  23.751026   \n",
       "4       413  49.837149  0.451932  0.000000  24.070813   \n",
       "...     ...        ...       ...       ...        ...   \n",
       "3101   9369  72.832347  0.626507  0.702141  18.645192   \n",
       "3102   9370  31.519831  2.903009  2.653274  44.564356   \n",
       "3103   9371  22.501694  1.313577  1.705099  47.108936   \n",
       "3104   9372  36.319039  0.566001  0.558300  39.219634   \n",
       "3105   9485  85.279899  0.395601  0.091769   8.467596   \n",
       "\n",
       "                                               geometry  \n",
       "0     POLYGON ((-2.18083 34.85145, -2.71695 34.77839...  \n",
       "1     POLYGON ((14.68659 35.53849, 14.14122 35.56449...  \n",
       "2     POLYGON ((-2.79499 35.22702, -3.33331 35.15029...  \n",
       "3     POLYGON ((23.95653 35.05738, 23.41935 35.13773...  \n",
       "4     POLYGON ((24.49271 34.97391, 23.95653 35.05738...  \n",
       "...                                                 ...  \n",
       "3101  POLYGON ((26.92771 70.58708, 25.63446 70.69412...  \n",
       "3102  POLYGON ((28.20775 70.47181, 26.92771 70.58708...  \n",
       "3103  POLYGON ((29.47384 70.34844, 28.20775 70.47181...  \n",
       "3104  POLYGON ((30.72534 70.21713, 29.47384 70.34844...  \n",
       "3105  POLYGON ((28.60702 70.90647, 27.30206 71.02446...  \n",
       "\n",
       "[3106 rows x 6 columns]"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#Cargar el grid\n",
    "#Grid_Habitats_Variables = gpd.read_file(r'Data/GLOBIO/RCP8_5/Grid_Habitats_Variables.shp')\n",
    "Grid_Habitats_Variables"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['C:/Users/molic/Desktop/TFM_V1/Chelsa_bio/ESCENARIOS_FUTUROS/Semejante_RCP_8_5/Estadistica_Zona\\\\Bio_01.shp',\n",
       " 'C:/Users/molic/Desktop/TFM_V1/Chelsa_bio/ESCENARIOS_FUTUROS/Semejante_RCP_8_5/Estadistica_Zona\\\\Bio_02.shp',\n",
       " 'C:/Users/molic/Desktop/TFM_V1/Chelsa_bio/ESCENARIOS_FUTUROS/Semejante_RCP_8_5/Estadistica_Zona\\\\Bio_03.shp',\n",
       " 'C:/Users/molic/Desktop/TFM_V1/Chelsa_bio/ESCENARIOS_FUTUROS/Semejante_RCP_8_5/Estadistica_Zona\\\\Bio_04.shp',\n",
       " 'C:/Users/molic/Desktop/TFM_V1/Chelsa_bio/ESCENARIOS_FUTUROS/Semejante_RCP_8_5/Estadistica_Zona\\\\Bio_05.shp',\n",
       " 'C:/Users/molic/Desktop/TFM_V1/Chelsa_bio/ESCENARIOS_FUTUROS/Semejante_RCP_8_5/Estadistica_Zona\\\\Bio_08.shp',\n",
       " 'C:/Users/molic/Desktop/TFM_V1/Chelsa_bio/ESCENARIOS_FUTUROS/Semejante_RCP_8_5/Estadistica_Zona\\\\Bio_12.shp',\n",
       " 'C:/Users/molic/Desktop/TFM_V1/Chelsa_bio/ESCENARIOS_FUTUROS/Semejante_RCP_8_5/Estadistica_Zona\\\\Bio_15.shp']"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import glob\n",
    "import os\n",
    "\n",
    "directorio = 'Data/CHELSA/ESCENARIOS_FUTUROS/Semejante_RCP_8_5/Estadistica_Zona' \n",
    "archivos_shp = glob.glob(directorio + '/Bio_*.shp') \n",
    "archivos_shp"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      ET_ID COUNTRY   x   y   Bio_12_mea  \\\n",
      "0       249   Spain  -2  35   285.961681   \n",
      "1       280   Malta  14  36   436.340764   \n",
      "2       363   Spain  -3  35   305.008734   \n",
      "3       412  Greece  24  35   706.663425   \n",
      "4       413  Greece  24  35   646.812883   \n",
      "...     ...     ...  ..  ..          ...   \n",
      "3101   9369  Norway  26  71   939.118965   \n",
      "3102   9370  Norway  28  71  1070.334074   \n",
      "3103   9371  Norway  29  71  1016.507321   \n",
      "3104   9372  Norway  30  70   910.152189   \n",
      "3105   9485  Norway  28  71   841.002412   \n",
      "\n",
      "                                               geometry  \n",
      "0     POLYGON ((-2.18083 34.85145, -2.71695 34.77839...  \n",
      "1     POLYGON ((14.68659 35.53849, 14.14122 35.56449...  \n",
      "2     POLYGON ((-2.79499 35.22702, -3.33331 35.15029...  \n",
      "3     POLYGON ((23.95653 35.05738, 23.41935 35.13773...  \n",
      "4     POLYGON ((24.49271 34.97391, 23.95653 35.05738...  \n",
      "...                                                 ...  \n",
      "3101  POLYGON ((26.92771 70.58708, 25.63446 70.69412...  \n",
      "3102  POLYGON ((28.20775 70.47181, 26.92771 70.58708...  \n",
      "3103  POLYGON ((29.47384 70.34844, 28.20775 70.47181...  \n",
      "3104  POLYGON ((30.72534 70.21713, 29.47384 70.34844...  \n",
      "3105  POLYGON ((28.60702 70.90647, 27.30206 71.02446...  \n",
      "\n",
      "[3106 rows x 6 columns]\n"
     ]
    }
   ],
   "source": [
    "# Cargar cada archivo shapefile y asignarlo a una variable con el nombre del archivo\n",
    "for shp in archivos_shp:\n",
    "    # Obtener el nombre base del archivo sin la extensión\n",
    "    file_name = os.path.splitext(os.path.basename(shp))[0]\n",
    "    \n",
    "    # Cargar el shapefile en un GeoDataFrame\n",
    "    gdf = gpd.read_file(shp)\n",
    "    \n",
    "    # Asignar el GeoDataFrame a una variable con el nombre del archivo\n",
    "    globals()[file_name] = gdf\n",
    "\n",
    "print(Bio_12)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Guardar el resultado en un nuevo shapefile\n",
    "Bio_01 = Bio_01[['ET_ID', 'Bio_01_mea']]\n",
    "Bio_02 = Bio_02[['ET_ID', 'Bio_02_mea']]\n",
    "Bio_03 = Bio_03[['ET_ID', 'Bio_03_mea']]\n",
    "Bio_04 = Bio_04[['ET_ID', 'Bio_04_mea']]\n",
    "Bio_05 = Bio_05[['ET_ID', 'Bio_05_mea']]\n",
    "Bio_08 = Bio_08[['ET_ID', 'Bio_08_mea']]\n",
    "Bio_12 = Bio_12[['ET_ID', 'Bio_12_mea']]\n",
    "Bio_15 = Bio_15[['ET_ID', 'Bio_15_mea']]\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_names = ['Bio_01', 'Bio_02', 'Bio_03', 'Bio_04', 'Bio_05', 'Bio_08', 'Bio_12', 'Bio_15']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "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>ET_ID</th>\n",
       "      <th>%_H_1</th>\n",
       "      <th>%_H_2</th>\n",
       "      <th>%_H_3</th>\n",
       "      <th>%_H_4</th>\n",
       "      <th>geometry</th>\n",
       "      <th>Bio_01_mea</th>\n",
       "      <th>Bio_02_mea</th>\n",
       "      <th>Bio_03_mea</th>\n",
       "      <th>Bio_04_mea</th>\n",
       "      <th>Bio_05_mea</th>\n",
       "      <th>Bio_08_mea</th>\n",
       "      <th>Bio_12_mea</th>\n",
       "      <th>Bio_15_mea</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>35.270045</td>\n",
       "      <td>0.028119</td>\n",
       "      <td>0.118723</td>\n",
       "      <td>28.121682</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>21.897371</td>\n",
       "      <td>9.207563</td>\n",
       "      <td>0.363398</td>\n",
       "      <td>491.788839</td>\n",
       "      <td>33.930236</td>\n",
       "      <td>16.690062</td>\n",
       "      <td>285.961681</td>\n",
       "      <td>54.397891</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>90.181696</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.040251</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>22.521384</td>\n",
       "      <td>1.715520</td>\n",
       "      <td>0.110632</td>\n",
       "      <td>475.048984</td>\n",
       "      <td>30.947298</td>\n",
       "      <td>19.187448</td>\n",
       "      <td>436.340764</td>\n",
       "      <td>85.713708</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   ET_ID      %_H_1     %_H_2     %_H_3      %_H_4  \\\n",
       "0    249  35.270045  0.028119  0.118723  28.121682   \n",
       "1    280  90.181696  0.000000  0.000000   0.040251   \n",
       "\n",
       "                                            geometry  Bio_01_mea  Bio_02_mea  \\\n",
       "0  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...   21.897371    9.207563   \n",
       "1  POLYGON ((14.68659 35.53849, 14.14122 35.56449...   22.521384    1.715520   \n",
       "\n",
       "   Bio_03_mea  Bio_04_mea  Bio_05_mea  Bio_08_mea  Bio_12_mea  Bio_15_mea  \n",
       "0    0.363398  491.788839   33.930236   16.690062  285.961681   54.397891  \n",
       "1    0.110632  475.048984   30.947298   19.187448  436.340764   85.713708  "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "for df_name in df_names:\n",
    "    df = globals()[df_name] \n",
    "    Grid_Habitats_Variables = pd.merge(Grid_Habitats_Variables, df, on='ET_ID', how='left')\n",
    "\n",
    "Grid_Habitats_Variables.head(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "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>ET_ID</th>\n",
       "      <th>%_H_1</th>\n",
       "      <th>%_H_2</th>\n",
       "      <th>%_H_3</th>\n",
       "      <th>%_H_4</th>\n",
       "      <th>geometry</th>\n",
       "      <th>Bio_01_mea</th>\n",
       "      <th>Bio_02_mea</th>\n",
       "      <th>Bio_03_mea</th>\n",
       "      <th>Bio_04_mea</th>\n",
       "      <th>Bio_05_mea</th>\n",
       "      <th>Bio_08_mea</th>\n",
       "      <th>Bio_12_mea</th>\n",
       "      <th>Bio_15_mea</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>35.2700</td>\n",
       "      <td>0.0281</td>\n",
       "      <td>0.1187</td>\n",
       "      <td>28.1217</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>21.8974</td>\n",
       "      <td>9.2076</td>\n",
       "      <td>0.3634</td>\n",
       "      <td>491.7888</td>\n",
       "      <td>33.9302</td>\n",
       "      <td>16.6901</td>\n",
       "      <td>285.9617</td>\n",
       "      <td>54.3979</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>90.1817</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.0403</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>22.5214</td>\n",
       "      <td>1.7155</td>\n",
       "      <td>0.1106</td>\n",
       "      <td>475.0490</td>\n",
       "      <td>30.9473</td>\n",
       "      <td>19.1874</td>\n",
       "      <td>436.3408</td>\n",
       "      <td>85.7137</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   ET_ID    %_H_1   %_H_2   %_H_3    %_H_4  \\\n",
       "0    249  35.2700  0.0281  0.1187  28.1217   \n",
       "1    280  90.1817  0.0000  0.0000   0.0403   \n",
       "\n",
       "                                            geometry  Bio_01_mea  Bio_02_mea  \\\n",
       "0  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...     21.8974      9.2076   \n",
       "1  POLYGON ((14.68659 35.53849, 14.14122 35.56449...     22.5214      1.7155   \n",
       "\n",
       "   Bio_03_mea  Bio_04_mea  Bio_05_mea  Bio_08_mea  Bio_12_mea  Bio_15_mea  \n",
       "0      0.3634    491.7888     33.9302     16.6901    285.9617     54.3979  \n",
       "1      0.1106    475.0490     30.9473     19.1874    436.3408     85.7137  "
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Redondear todas las columnas float64 a 4 decimales\n",
    "Grid_Habitats_Variables = Grid_Habitats_Variables.round(4)\n",
    "\n",
    "# Asegurarse de que no hay duplicados\n",
    "Grid_Habitats_Variables = Grid_Habitats_Variables.drop_duplicates(subset=['ET_ID'])\n",
    "Grid_Habitats_Variables.head(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "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>ET_ID</th>\n",
       "      <th>%_H_1</th>\n",
       "      <th>%_H_2</th>\n",
       "      <th>%_H_3</th>\n",
       "      <th>%_H_4</th>\n",
       "      <th>geometry</th>\n",
       "      <th>Bio_01_mea</th>\n",
       "      <th>Bio_02_mea</th>\n",
       "      <th>Bio_03_mea</th>\n",
       "      <th>Bio_04_mea</th>\n",
       "      <th>Bio_05_mea</th>\n",
       "      <th>Bio_08_mea</th>\n",
       "      <th>Bio_12_mea</th>\n",
       "      <th>Bio_15_mea</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>249</td>\n",
       "      <td>35.2700</td>\n",
       "      <td>0.0281</td>\n",
       "      <td>0.1187</td>\n",
       "      <td>28.1217</td>\n",
       "      <td>POLYGON ((-2.18083 34.85145, -2.71695 34.77839...</td>\n",
       "      <td>21.8974</td>\n",
       "      <td>9.2076</td>\n",
       "      <td>0.3634</td>\n",
       "      <td>491.7888</td>\n",
       "      <td>33.9302</td>\n",
       "      <td>16.6901</td>\n",
       "      <td>285.9617</td>\n",
       "      <td>54.3979</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>280</td>\n",
       "      <td>90.1817</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>0.0403</td>\n",
       "      <td>POLYGON ((14.68659 35.53849, 14.14122 35.56449...</td>\n",
       "      <td>22.5214</td>\n",
       "      <td>1.7155</td>\n",
       "      <td>0.1106</td>\n",
       "      <td>475.0490</td>\n",
       "      <td>30.9473</td>\n",
       "      <td>19.1874</td>\n",
       "      <td>436.3408</td>\n",
       "      <td>85.7137</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>363</td>\n",
       "      <td>83.4716</td>\n",
       "      <td>0.0809</td>\n",
       "      <td>0.0187</td>\n",
       "      <td>8.7956</td>\n",
       "      <td>POLYGON ((-2.79499 35.22702, -3.33331 35.15029...</td>\n",
       "      <td>21.4139</td>\n",
       "      <td>4.2293</td>\n",
       "      <td>0.2429</td>\n",
       "      <td>395.9235</td>\n",
       "      <td>30.0194</td>\n",
       "      <td>17.0383</td>\n",
       "      <td>305.0087</td>\n",
       "      <td>60.2356</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   ET_ID    %_H_1   %_H_2   %_H_3    %_H_4  \\\n",
       "0    249  35.2700  0.0281  0.1187  28.1217   \n",
       "1    280  90.1817  0.0000  0.0000   0.0403   \n",
       "2    363  83.4716  0.0809  0.0187   8.7956   \n",
       "\n",
       "                                            geometry  Bio_01_mea  Bio_02_mea  \\\n",
       "0  POLYGON ((-2.18083 34.85145, -2.71695 34.77839...     21.8974      9.2076   \n",
       "1  POLYGON ((14.68659 35.53849, 14.14122 35.56449...     22.5214      1.7155   \n",
       "2  POLYGON ((-2.79499 35.22702, -3.33331 35.15029...     21.4139      4.2293   \n",
       "\n",
       "   Bio_03_mea  Bio_04_mea  Bio_05_mea  Bio_08_mea  Bio_12_mea  Bio_15_mea  \n",
       "0      0.3634    491.7888     33.9302     16.6901    285.9617     54.3979  \n",
       "1      0.1106    475.0490     30.9473     19.1874    436.3408     85.7137  \n",
       "2      0.2429    395.9235     30.0194     17.0383    305.0087     60.2356  "
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Redondear todas las columnas float64 a 4 decimales\n",
    "Grid_Habitats_Variables = Grid_Habitats_Variables.round(4)\n",
    "Grid_Habitats_Variables.head(3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Exportar el resultado\n",
    "Grid_Habitats_Variables.to_csv('Data/Variables/RCP8_5/Grid_Variables_RCP8_5.csv', index=False, decimal='.')\n",
    "Grid_Habitats_Variables.to_file('Data/Variables/RCP8_5/Grid_Variables_RCP8_5.shp')\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "anura",
   "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.12.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}