TRAITEMENTS_BDMOS_1999_2020.ipynb

{
 "cells": [
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "##### Evaluer les pertes de terres naturelles et agricoles en 20 ans de croissance de l'agglomération nantaise\n",
    "Traitement de la BDMOS de Loire-Atlantique, millésimes 1999, 2004, 2009, 2012, 2016 et 2020. \n",
    "\n",
    "*Source : https://data.loire-atlantique.fr/explore/dataset/224400028_modes-occupation-sols-loire-atlantique/information/?disjunctive.commune&disjunctive.epci&disjunctive.scot&disjunctive.insee_dep&disjunctive.delegation&disjunctive.ocs_1_1949&disjunctive.ocs_1_1999&disjunctive.ocs_2_1999&disjunctive.ocs_3_1999&disjunctive.ocs_1_2004&disjunctive.ocs_2_2004&disjunctive.ocs_3_2004&disjunctive.ocs_1_2009&disjunctive.ocs_2_2009&disjunctive.ocs_3_2009&disjunctive.ocs_2_2012&disjunctive.ocs_3_2012&disjunctive.ocs_1_2016&disjunctive.ocs_2_2016&disjunctive.ocs_3_2016&disjunctive.ocs_1_2020&disjunctive.ocs_2_2020&disjunctive.ocs_3_2020&sort=objectid*\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import geopandas as gpd\n",
    "from math import pi\n",
    "import numpy as np\n",
    "import fiona\n",
    "import matplotlib as mpl\n",
    "import matplotlib.pyplot as plt\n",
    "pd.set_option(\"display.max_columns\", None)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Chargement du fichier geojson de la BDMOS : https://data.loire-atlantique.fr/explore/dataset/224400028_modes-occupation-sols-loire-atlantique/information/?disjunctive.commune&disjunctive.epci&disjunctive.scot&disjunctive.insee_dep&disjunctive.delegation&disjunctive.ocs_1_1949&disjunctive.ocs_1_1999&disjunctive.ocs_2_1999&disjunctive.ocs_3_1999&disjunctive.ocs_1_2004&disjunctive.ocs_2_2004&disjunctive.ocs_3_2004&disjunctive.ocs_1_2009&disjunctive.ocs_2_2009&disjunctive.ocs_3_2009&disjunctive.ocs_2_2012&disjunctive.ocs_3_2012&disjunctive.ocs_1_2016&disjunctive.ocs_2_2016&disjunctive.ocs_3_2016&disjunctive.ocs_1_2020&disjunctive.ocs_2_2020&disjunctive.ocs_3_2020&sort=objectid \n",
    "mos=gpd.read_file(\"MOS-loire-atlantique_1949_2020.geojson\")\n",
    "\n",
    "# Communes de l'aire d'attraction de Nantes :\n",
    "aav=[\"44009\",\"44018\",\"44020\",\"44024\",\"44026\",\"44035\",\"44047\",\"44074\",\"44094\",\"44101\",\"44109\",\"44114\",\"44120\",\"44143\",\"44150\",\"44162\",\"44166\",\"44171\",\"44172\",\"44190\",\"44194\",\"44198\",\"44204\",\"44215\",\"44001\",\"44002\",\"44005\",\"44014\",\"44015\",\"44016\",\"44019\",\"44023\",\"44025\",\"44027\",\"44028\",\"44029\",\"44032\",\"44033\",\"44037\",\"44039\",\"44041\",\"44043\",\"44045\",\"44048\",\"44056\",\"44061\",\"44062\",\"44063\",\"44064\",\"44066\",\"44070\",\"44071\",\"44073\",\"44076\",\"44077\",\"44079\",\"44080\",\"44082\",\"44083\",\"44084\",\"44088\",\"44089\",\"44090\",\"44091\",\"44100\",\"44102\",\"44107\",\"44108\",\"44110\",\"44111\",\"44113\",\"44115\",\"44117\",\"44119\",\"44122\",\"44127\",\"44130\",\"44133\",\"44137\",\"44138\",\"44139\",\"44140\",\"44141\",\"44142\",\"44145\",\"44149\",\"44155\",\"44156\",\"44158\",\"44159\",\"44164\",\"44165\",\"44169\",\"44173\",\"44174\",\"44178\",\"44179\",\"44186\",\"44188\",\"44195\",\"44201\",\"44202\",\"44203\",\"44205\",\"44207\",\"44208\",\"44209\",\"44212\",\"44214\",\"44216\",\"44217\",\"44220\",\"44221\",\"44223\",\"44224\"]\n",
    "# Communes de Nantes Métropole :\n",
    "nm=['44172', '44120', '44009', '44024', '44020', '44204', '44047','44166', '44114', '44094', '44109', '44171', '44194', '44150','44074', '44035', '44101', '44026', '44143', '44018', '44198','44162', '44215', '44190']\n",
    "# Communes de l'aire d'attraction exterieures à Nantes Métropole :\n",
    "hors_nm=[\"44001\",\"44002\",\"44005\",\"44014\",\"44015\",\"44016\",\"44019\",\"44023\",\"44025\",\"44027\",\"44028\",\"44029\",\"44032\",\"44033\",\"44037\",\"44039\",\"44041\",\"44043\",\"44045\",\"44048\",\"44056\",\"44061\",\"44062\",\"44063\",\"44064\",\"44066\",\"44070\",\"44071\",\"44073\",\"44076\",\"44077\",\"44079\",\"44080\",\"44082\",\"44083\",\"44084\",\"44088\",\"44089\",\"44090\",\"44091\",\"44100\",\"44102\",\"44107\",\"44108\",\"44110\",\"44111\",\"44113\",\"44115\",\"44117\",\"44119\",\"44122\",\"44127\",\"44130\",\"44133\",\"44137\",\"44138\",\"44139\",\"44140\",\"44141\",\"44142\",\"44145\",\"44149\",\"44155\",\"44156\",\"44158\",\"44159\",\"44164\",\"44165\",\"44169\",\"44173\",\"44174\",\"44178\",\"44179\",\"44186\",\"44188\",\"44195\",\"44201\",\"44202\",\"44203\",\"44205\",\"44207\",\"44208\",\"44209\",\"44212\",\"44214\",\"44216\",\"44217\",\"44220\",\"44221\",\"44223\",\"44224\"]\n",
    "\n",
    "# Ajout d'un variable qui précise le secteur gégraphique de chaque parcelle (Communes de Nantes Métropole / Autres communes de l'aire d'attraction)\n",
    "mos=mos[mos['insee_com'].isin(aav)]\n",
    "mos.loc[mos['insee_com'].astype('str').isin(nm), \"SECTEUR\"] = \"NM\"\n",
    "mos.loc[mos['insee_com'].astype('str').isin(hors_nm), \"SECTEUR\"] = \"HORS_NM\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "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>ocs_code_2_2009</th>\n",
       "      <th>ocs_2_2012</th>\n",
       "      <th>siren_epci</th>\n",
       "      <th>scot</th>\n",
       "      <th>ocs_code_3_2012</th>\n",
       "      <th>ocs_code_3_1999</th>\n",
       "      <th>ocs_1_2004</th>\n",
       "      <th>ocs_code_1_2012</th>\n",
       "      <th>ocs_1_2020</th>\n",
       "      <th>delegation</th>\n",
       "      <th>ocs_1_1999</th>\n",
       "      <th>ocs_code_1_2004</th>\n",
       "      <th>epci</th>\n",
       "      <th>insee_com</th>\n",
       "      <th>ocs_3_2009</th>\n",
       "      <th>ocs_3_2004</th>\n",
       "      <th>ocs_3_2020</th>\n",
       "      <th>ocs_2_2009</th>\n",
       "      <th>surface_ha</th>\n",
       "      <th>ocs_3_2016</th>\n",
       "      <th>ocs_code_1_1999</th>\n",
       "      <th>ocs_code_1_2009</th>\n",
       "      <th>ocs_code_3_2004</th>\n",
       "      <th>commune</th>\n",
       "      <th>insee_dep</th>\n",
       "      <th>ocs_2_2004</th>\n",
       "      <th>ocs_code_2_2012</th>\n",
       "      <th>ocs_1_2009</th>\n",
       "      <th>ocs_code_3_2016</th>\n",
       "      <th>ocs_code_1_2020</th>\n",
       "      <th>ocs_code_3_2009</th>\n",
       "      <th>ocs_1_2016</th>\n",
       "      <th>objectid</th>\n",
       "      <th>ocs_code_1_2016</th>\n",
       "      <th>ocs_3_1999</th>\n",
       "      <th>ocs_code_2_2016</th>\n",
       "      <th>ocs_code_2_1999</th>\n",
       "      <th>ocs_2_1999</th>\n",
       "      <th>ocs_2_2016</th>\n",
       "      <th>ocs_code_2_2004</th>\n",
       "      <th>ocs_code_1_1949</th>\n",
       "      <th>ocs_1_1949</th>\n",
       "      <th>ocs_2_2020</th>\n",
       "      <th>ocs_3_2012</th>\n",
       "      <th>ocs_code_3_2020</th>\n",
       "      <th>ocs_1_2012</th>\n",
       "      <th>ocs_code_2_2020</th>\n",
       "      <th>geometry</th>\n",
       "      <th>SECTEUR</th>\n",
       "      <th>area</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>10100</th>\n",
       "      <td>21</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>244400537</td>\n",
       "      <td>SCOT du Pays de Châteaubriant</td>\n",
       "      <td>231</td>\n",
       "      <td>231</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>2</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>Délégation Châteaubriant</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>2</td>\n",
       "      <td>CC de Nozay</td>\n",
       "      <td>44138</td>\n",
       "      <td>Cultures annuelles et terres labourées</td>\n",
       "      <td>Cultures annuelles et terres labourées</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Terres arables</td>\n",
       "      <td>4.478066</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>211</td>\n",
       "      <td>Puceul</td>\n",
       "      <td>44</td>\n",
       "      <td>Terres arables</td>\n",
       "      <td>23</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>231</td>\n",
       "      <td>2</td>\n",
       "      <td>211</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>10530</td>\n",
       "      <td>2</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>23</td>\n",
       "      <td>23</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>21</td>\n",
       "      <td>2</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>231</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>23</td>\n",
       "      <td>POLYGON ((355433.859 6725027.737, 355364.408 6...</td>\n",
       "      <td>HORS_NM</td>\n",
       "      <td>44780.659751</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>369875</th>\n",
       "      <td>21</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>200072734</td>\n",
       "      <td>SCOT de la métropole Nantes - Saint-Nazaire</td>\n",
       "      <td>231</td>\n",
       "      <td>211</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>2</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>Délégation Saint-Nazaire</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>2</td>\n",
       "      <td>CC Estuaire et Sillon</td>\n",
       "      <td>44025</td>\n",
       "      <td>Cultures annuelles et terres labourées</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Terres arables</td>\n",
       "      <td>4.006541</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>231</td>\n",
       "      <td>Campbon</td>\n",
       "      <td>44</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>23</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>231</td>\n",
       "      <td>2</td>\n",
       "      <td>211</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>397902</td>\n",
       "      <td>2</td>\n",
       "      <td>Cultures annuelles et terres labourées</td>\n",
       "      <td>23</td>\n",
       "      <td>21</td>\n",
       "      <td>Terres arables</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>23</td>\n",
       "      <td>2</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>231</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>23</td>\n",
       "      <td>POLYGON ((322334.256 6714600.310, 322514.540 6...</td>\n",
       "      <td>HORS_NM</td>\n",
       "      <td>40065.409911</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>380484</th>\n",
       "      <td>21</td>\n",
       "      <td>Terres arables</td>\n",
       "      <td>244400552</td>\n",
       "      <td>SCOT du Pays d'Ancenis</td>\n",
       "      <td>211</td>\n",
       "      <td>231</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>2</td>\n",
       "      <td>Forêts et milieux semi naturels</td>\n",
       "      <td>Délégation Ancenis</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>2</td>\n",
       "      <td>CC du Pays d'Ancenis</td>\n",
       "      <td>44077</td>\n",
       "      <td>Cultures annuelles et terres labourées</td>\n",
       "      <td>Cultures annuelles et terres labourées</td>\n",
       "      <td>Peuplements d'arbres en alignement</td>\n",
       "      <td>Terres arables</td>\n",
       "      <td>0.000009</td>\n",
       "      <td>Cultures annuelles et terres labourées</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>211</td>\n",
       "      <td>Joué-sur-Erdre</td>\n",
       "      <td>44</td>\n",
       "      <td>Terres arables</td>\n",
       "      <td>21</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>211</td>\n",
       "      <td>3</td>\n",
       "      <td>211</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>417795</td>\n",
       "      <td>2</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>21</td>\n",
       "      <td>23</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Terres arables</td>\n",
       "      <td>21</td>\n",
       "      <td>2</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>Cultures annuelles et terres labourées</td>\n",
       "      <td>312</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>31</td>\n",
       "      <td>POLYGON ((369616.953 6722384.551, 369537.719 6...</td>\n",
       "      <td>HORS_NM</td>\n",
       "      <td>0.094406</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>246996</th>\n",
       "      <td>11</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>200067866</td>\n",
       "      <td>SCOT du Pays du Vignoble nantais</td>\n",
       "      <td>112</td>\n",
       "      <td>112</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>1</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>Délégation vignoble</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>1</td>\n",
       "      <td>CC Sèvre et Loire</td>\n",
       "      <td>44108</td>\n",
       "      <td>Hameau rural</td>\n",
       "      <td>Hameau rural</td>\n",
       "      <td>Hameau rural</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>0.000005</td>\n",
       "      <td>Hameau rural</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>112</td>\n",
       "      <td>Mouzillon</td>\n",
       "      <td>44</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>11</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>112</td>\n",
       "      <td>1</td>\n",
       "      <td>112</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>267504</td>\n",
       "      <td>1</td>\n",
       "      <td>Hameau rural</td>\n",
       "      <td>11</td>\n",
       "      <td>11</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>11</td>\n",
       "      <td>2</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>Hameau rural</td>\n",
       "      <td>112</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>11</td>\n",
       "      <td>POLYGON ((377978.341 6676903.254, 377970.136 6...</td>\n",
       "      <td>HORS_NM</td>\n",
       "      <td>0.050482</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>283578</th>\n",
       "      <td>11</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>244400404</td>\n",
       "      <td>SCOT de la métropole Nantes - Saint-Nazaire</td>\n",
       "      <td>114</td>\n",
       "      <td>222</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>1</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>Délégation Nantes</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>2</td>\n",
       "      <td>Nantes Métropole</td>\n",
       "      <td>44198</td>\n",
       "      <td>Habitat pavillonnaire</td>\n",
       "      <td>Maraîchage, horticulture</td>\n",
       "      <td>Habitat pavillonnaire</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>0.308112</td>\n",
       "      <td>Habitat pavillonnaire</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>222</td>\n",
       "      <td>Les Sorinières</td>\n",
       "      <td>44</td>\n",
       "      <td>Cultures permanentes et maraîchage</td>\n",
       "      <td>11</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>114</td>\n",
       "      <td>1</td>\n",
       "      <td>114</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>304086</td>\n",
       "      <td>1</td>\n",
       "      <td>Maraîchage, horticulture</td>\n",
       "      <td>11</td>\n",
       "      <td>22</td>\n",
       "      <td>Cultures permanentes et maraîchage</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>22</td>\n",
       "      <td>2</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>Habitat pavillonnaire</td>\n",
       "      <td>114</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>11</td>\n",
       "      <td>POLYGON ((357754.613 6681606.781, 357751.447 6...</td>\n",
       "      <td>NM</td>\n",
       "      <td>3081.119196</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14195</th>\n",
       "      <td>23</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>244400503</td>\n",
       "      <td>SCOT de la métropole Nantes - Saint-Nazaire</td>\n",
       "      <td>231</td>\n",
       "      <td>231</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>2</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>Délégation Châteaubriant</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>2</td>\n",
       "      <td>CC d'Erdre et Gesvres</td>\n",
       "      <td>44073</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>0.898499</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>231</td>\n",
       "      <td>Héric</td>\n",
       "      <td>44</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>23</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>231</td>\n",
       "      <td>2</td>\n",
       "      <td>231</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>14625</td>\n",
       "      <td>2</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>23</td>\n",
       "      <td>23</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>23</td>\n",
       "      <td>2</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>231</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>23</td>\n",
       "      <td>POLYGON ((354865.204 6713235.322, 354844.574 6...</td>\n",
       "      <td>HORS_NM</td>\n",
       "      <td>8984.990985</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>384012</th>\n",
       "      <td>11</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>244400404</td>\n",
       "      <td>SCOT de la métropole Nantes - Saint-Nazaire</td>\n",
       "      <td>114</td>\n",
       "      <td>114</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>1</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>Délégation Nantes</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>1</td>\n",
       "      <td>Nantes Métropole</td>\n",
       "      <td>44047</td>\n",
       "      <td>Habitat pavillonnaire</td>\n",
       "      <td>Habitat pavillonnaire</td>\n",
       "      <td>Habitat pavillonnaire</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>2537.327937</td>\n",
       "      <td>Habitat pavillonnaire</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>114</td>\n",
       "      <td>Couëron</td>\n",
       "      <td>44</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>11</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>114</td>\n",
       "      <td>1</td>\n",
       "      <td>114</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>421323</td>\n",
       "      <td>1</td>\n",
       "      <td>Habitat pavillonnaire</td>\n",
       "      <td>11</td>\n",
       "      <td>11</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>11</td>\n",
       "      <td>3</td>\n",
       "      <td>Forêts et milieux semi naturels</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>Habitat pavillonnaire</td>\n",
       "      <td>114</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>11</td>\n",
       "      <td>POLYGON ((343613.007 6690399.642, 343604.392 6...</td>\n",
       "      <td>NM</td>\n",
       "      <td>2537.327937</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6673</th>\n",
       "      <td>23</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>244400537</td>\n",
       "      <td>SCOT du Pays de Châteaubriant</td>\n",
       "      <td>231</td>\n",
       "      <td>231</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>2</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>Délégation Châteaubriant</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>2</td>\n",
       "      <td>CC de Nozay</td>\n",
       "      <td>44001</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>0.195641</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>231</td>\n",
       "      <td>Abbaretz</td>\n",
       "      <td>44</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>23</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>231</td>\n",
       "      <td>2</td>\n",
       "      <td>231</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>6676</td>\n",
       "      <td>2</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>23</td>\n",
       "      <td>23</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>23</td>\n",
       "      <td>2</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>Prairies</td>\n",
       "      <td>231</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>23</td>\n",
       "      <td>POLYGON ((360677.249 6730906.045, 360680.982 6...</td>\n",
       "      <td>HORS_NM</td>\n",
       "      <td>1956.407206</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>102491</th>\n",
       "      <td>31</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>244400404</td>\n",
       "      <td>SCOT de la métropole Nantes - Saint-Nazaire</td>\n",
       "      <td>311</td>\n",
       "      <td>311</td>\n",
       "      <td>Forêts et milieux semi naturels</td>\n",
       "      <td>3</td>\n",
       "      <td>Forêts et milieux semi naturels</td>\n",
       "      <td>Délégation Nantes</td>\n",
       "      <td>Forêts et milieux semi naturels</td>\n",
       "      <td>3</td>\n",
       "      <td>Nantes Métropole</td>\n",
       "      <td>44018</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>0.090154</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>3</td>\n",
       "      <td>3</td>\n",
       "      <td>311</td>\n",
       "      <td>Bouaye</td>\n",
       "      <td>44</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>31</td>\n",
       "      <td>Forêts et milieux semi naturels</td>\n",
       "      <td>311</td>\n",
       "      <td>3</td>\n",
       "      <td>311</td>\n",
       "      <td>Forêts et milieux semi naturels</td>\n",
       "      <td>105136</td>\n",
       "      <td>3</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>31</td>\n",
       "      <td>31</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>31</td>\n",
       "      <td>3</td>\n",
       "      <td>Forêts et milieux semi naturels</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>Bois et forêts</td>\n",
       "      <td>311</td>\n",
       "      <td>Forêts et milieux semi naturels</td>\n",
       "      <td>31</td>\n",
       "      <td>POLYGON ((344121.222 6682704.427, 344121.469 6...</td>\n",
       "      <td>NM</td>\n",
       "      <td>901.543248</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>213891</th>\n",
       "      <td>11</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>200067635</td>\n",
       "      <td>SCOT du Pays du Vignoble nantais</td>\n",
       "      <td>112</td>\n",
       "      <td>112</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>1</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>Délégation vignoble</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>1</td>\n",
       "      <td>CA Clisson Sèvre et Maine Agglo</td>\n",
       "      <td>44165</td>\n",
       "      <td>Hameau rural</td>\n",
       "      <td>Hameau rural</td>\n",
       "      <td>Hameau rural</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>1.152150</td>\n",
       "      <td>Hameau rural</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>112</td>\n",
       "      <td>Saint-Hilaire-de-Clisson</td>\n",
       "      <td>44</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>11</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>112</td>\n",
       "      <td>1</td>\n",
       "      <td>112</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>234399</td>\n",
       "      <td>1</td>\n",
       "      <td>Hameau rural</td>\n",
       "      <td>11</td>\n",
       "      <td>11</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>11</td>\n",
       "      <td>2</td>\n",
       "      <td>Territoires agricoles</td>\n",
       "      <td>Zones Urbanisées liées à l'habitat</td>\n",
       "      <td>Hameau rural</td>\n",
       "      <td>112</td>\n",
       "      <td>Territoires artificialisés</td>\n",
       "      <td>11</td>\n",
       "      <td>MULTIPOLYGON (((372491.875 6668534.827, 372485...</td>\n",
       "      <td>HORS_NM</td>\n",
       "      <td>11521.504251</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        ocs_code_2_2009                          ocs_2_2012 siren_epci  \\\n",
       "10100                21                            Prairies  244400537   \n",
       "369875               21                            Prairies  200072734   \n",
       "380484               21                      Terres arables  244400552   \n",
       "246996               11  Zones Urbanisées liées à l'habitat  200067866   \n",
       "283578               11  Zones Urbanisées liées à l'habitat  244400404   \n",
       "14195                23                            Prairies  244400503   \n",
       "384012               11  Zones Urbanisées liées à l'habitat  244400404   \n",
       "6673                 23                            Prairies  244400537   \n",
       "102491               31                      Bois et forêts  244400404   \n",
       "213891               11  Zones Urbanisées liées à l'habitat  200067635   \n",
       "\n",
       "                                               scot ocs_code_3_2012  \\\n",
       "10100                 SCOT du Pays de Châteaubriant             231   \n",
       "369875  SCOT de la métropole Nantes - Saint-Nazaire             231   \n",
       "380484                       SCOT du Pays d'Ancenis             211   \n",
       "246996             SCOT du Pays du Vignoble nantais             112   \n",
       "283578  SCOT de la métropole Nantes - Saint-Nazaire             114   \n",
       "14195   SCOT de la métropole Nantes - Saint-Nazaire             231   \n",
       "384012  SCOT de la métropole Nantes - Saint-Nazaire             114   \n",
       "6673                  SCOT du Pays de Châteaubriant             231   \n",
       "102491  SCOT de la métropole Nantes - Saint-Nazaire             311   \n",
       "213891             SCOT du Pays du Vignoble nantais             112   \n",
       "\n",
       "       ocs_code_3_1999                       ocs_1_2004  ocs_code_1_2012  \\\n",
       "10100              231            Territoires agricoles                2   \n",
       "369875             211            Territoires agricoles                2   \n",
       "380484             231            Territoires agricoles                2   \n",
       "246996             112       Territoires artificialisés                1   \n",
       "283578             222            Territoires agricoles                1   \n",
       "14195              231            Territoires agricoles                2   \n",
       "384012             114       Territoires artificialisés                1   \n",
       "6673               231            Territoires agricoles                2   \n",
       "102491             311  Forêts et milieux semi naturels                3   \n",
       "213891             112       Territoires artificialisés                1   \n",
       "\n",
       "                             ocs_1_2020                delegation  \\\n",
       "10100             Territoires agricoles  Délégation Châteaubriant   \n",
       "369875            Territoires agricoles  Délégation Saint-Nazaire   \n",
       "380484  Forêts et milieux semi naturels        Délégation Ancenis   \n",
       "246996       Territoires artificialisés       Délégation vignoble   \n",
       "283578       Territoires artificialisés         Délégation Nantes   \n",
       "14195             Territoires agricoles  Délégation Châteaubriant   \n",
       "384012       Territoires artificialisés         Délégation Nantes   \n",
       "6673              Territoires agricoles  Délégation Châteaubriant   \n",
       "102491  Forêts et milieux semi naturels         Délégation Nantes   \n",
       "213891       Territoires artificialisés       Délégation vignoble   \n",
       "\n",
       "                             ocs_1_1999  ocs_code_1_2004  \\\n",
       "10100             Territoires agricoles                2   \n",
       "369875            Territoires agricoles                2   \n",
       "380484            Territoires agricoles                2   \n",
       "246996       Territoires artificialisés                1   \n",
       "283578            Territoires agricoles                2   \n",
       "14195             Territoires agricoles                2   \n",
       "384012       Territoires artificialisés                1   \n",
       "6673              Territoires agricoles                2   \n",
       "102491  Forêts et milieux semi naturels                3   \n",
       "213891       Territoires artificialisés                1   \n",
       "\n",
       "                                   epci insee_com  \\\n",
       "10100                       CC de Nozay     44138   \n",
       "369875            CC Estuaire et Sillon     44025   \n",
       "380484             CC du Pays d'Ancenis     44077   \n",
       "246996                CC Sèvre et Loire     44108   \n",
       "283578                 Nantes Métropole     44198   \n",
       "14195             CC d'Erdre et Gesvres     44073   \n",
       "384012                 Nantes Métropole     44047   \n",
       "6673                        CC de Nozay     44001   \n",
       "102491                 Nantes Métropole     44018   \n",
       "213891  CA Clisson Sèvre et Maine Agglo     44165   \n",
       "\n",
       "                                    ocs_3_2009  \\\n",
       "10100   Cultures annuelles et terres labourées   \n",
       "369875  Cultures annuelles et terres labourées   \n",
       "380484  Cultures annuelles et terres labourées   \n",
       "246996                            Hameau rural   \n",
       "283578                   Habitat pavillonnaire   \n",
       "14195                                 Prairies   \n",
       "384012                   Habitat pavillonnaire   \n",
       "6673                                  Prairies   \n",
       "102491                          Bois et forêts   \n",
       "213891                            Hameau rural   \n",
       "\n",
       "                                    ocs_3_2004  \\\n",
       "10100   Cultures annuelles et terres labourées   \n",
       "369875                                Prairies   \n",
       "380484  Cultures annuelles et terres labourées   \n",
       "246996                            Hameau rural   \n",
       "283578                Maraîchage, horticulture   \n",
       "14195                                 Prairies   \n",
       "384012                   Habitat pavillonnaire   \n",
       "6673                                  Prairies   \n",
       "102491                          Bois et forêts   \n",
       "213891                            Hameau rural   \n",
       "\n",
       "                                ocs_3_2020  \\\n",
       "10100                             Prairies   \n",
       "369875                            Prairies   \n",
       "380484  Peuplements d'arbres en alignement   \n",
       "246996                        Hameau rural   \n",
       "283578               Habitat pavillonnaire   \n",
       "14195                             Prairies   \n",
       "384012               Habitat pavillonnaire   \n",
       "6673                              Prairies   \n",
       "102491                      Bois et forêts   \n",
       "213891                        Hameau rural   \n",
       "\n",
       "                                ocs_2_2009   surface_ha  \\\n",
       "10100                       Terres arables     4.478066   \n",
       "369875                      Terres arables     4.006541   \n",
       "380484                      Terres arables     0.000009   \n",
       "246996  Zones Urbanisées liées à l'habitat     0.000005   \n",
       "283578  Zones Urbanisées liées à l'habitat     0.308112   \n",
       "14195                             Prairies     0.898499   \n",
       "384012  Zones Urbanisées liées à l'habitat  2537.327937   \n",
       "6673                              Prairies     0.195641   \n",
       "102491                      Bois et forêts     0.090154   \n",
       "213891  Zones Urbanisées liées à l'habitat     1.152150   \n",
       "\n",
       "                                    ocs_3_2016  ocs_code_1_1999  \\\n",
       "10100                                 Prairies                2   \n",
       "369875                                Prairies                2   \n",
       "380484  Cultures annuelles et terres labourées                2   \n",
       "246996                            Hameau rural                1   \n",
       "283578                   Habitat pavillonnaire                2   \n",
       "14195                                 Prairies                2   \n",
       "384012                   Habitat pavillonnaire                1   \n",
       "6673                                  Prairies                2   \n",
       "102491                          Bois et forêts                3   \n",
       "213891                            Hameau rural                1   \n",
       "\n",
       "        ocs_code_1_2009 ocs_code_3_2004                   commune insee_dep  \\\n",
       "10100                 2             211                    Puceul        44   \n",
       "369875                2             231                   Campbon        44   \n",
       "380484                2             211            Joué-sur-Erdre        44   \n",
       "246996                1             112                 Mouzillon        44   \n",
       "283578                1             222            Les Sorinières        44   \n",
       "14195                 2             231                     Héric        44   \n",
       "384012                1             114                   Couëron        44   \n",
       "6673                  2             231                  Abbaretz        44   \n",
       "102491                3             311                    Bouaye        44   \n",
       "213891                1             112  Saint-Hilaire-de-Clisson        44   \n",
       "\n",
       "                                ocs_2_2004  ocs_code_2_2012  \\\n",
       "10100                       Terres arables               23   \n",
       "369875                            Prairies               23   \n",
       "380484                      Terres arables               21   \n",
       "246996  Zones Urbanisées liées à l'habitat               11   \n",
       "283578  Cultures permanentes et maraîchage               11   \n",
       "14195                             Prairies               23   \n",
       "384012  Zones Urbanisées liées à l'habitat               11   \n",
       "6673                              Prairies               23   \n",
       "102491                      Bois et forêts               31   \n",
       "213891  Zones Urbanisées liées à l'habitat               11   \n",
       "\n",
       "                             ocs_1_2009 ocs_code_3_2016  ocs_code_1_2020  \\\n",
       "10100             Territoires agricoles             231                2   \n",
       "369875            Territoires agricoles             231                2   \n",
       "380484            Territoires agricoles             211                3   \n",
       "246996       Territoires artificialisés             112                1   \n",
       "283578       Territoires artificialisés             114                1   \n",
       "14195             Territoires agricoles             231                2   \n",
       "384012       Territoires artificialisés             114                1   \n",
       "6673              Territoires agricoles             231                2   \n",
       "102491  Forêts et milieux semi naturels             311                3   \n",
       "213891       Territoires artificialisés             112                1   \n",
       "\n",
       "       ocs_code_3_2009                       ocs_1_2016 objectid  \\\n",
       "10100              211            Territoires agricoles    10530   \n",
       "369875             211            Territoires agricoles   397902   \n",
       "380484             211            Territoires agricoles   417795   \n",
       "246996             112       Territoires artificialisés   267504   \n",
       "283578             114       Territoires artificialisés   304086   \n",
       "14195              231            Territoires agricoles    14625   \n",
       "384012             114       Territoires artificialisés   421323   \n",
       "6673               231            Territoires agricoles     6676   \n",
       "102491             311  Forêts et milieux semi naturels   105136   \n",
       "213891             112       Territoires artificialisés   234399   \n",
       "\n",
       "        ocs_code_1_2016                              ocs_3_1999  \\\n",
       "10100                 2                                Prairies   \n",
       "369875                2  Cultures annuelles et terres labourées   \n",
       "380484                2                                Prairies   \n",
       "246996                1                            Hameau rural   \n",
       "283578                1                Maraîchage, horticulture   \n",
       "14195                 2                                Prairies   \n",
       "384012                1                   Habitat pavillonnaire   \n",
       "6673                  2                                Prairies   \n",
       "102491                3                          Bois et forêts   \n",
       "213891                1                            Hameau rural   \n",
       "\n",
       "        ocs_code_2_2016  ocs_code_2_1999                          ocs_2_1999  \\\n",
       "10100                23               23                            Prairies   \n",
       "369875               23               21                      Terres arables   \n",
       "380484               21               23                            Prairies   \n",
       "246996               11               11  Zones Urbanisées liées à l'habitat   \n",
       "283578               11               22  Cultures permanentes et maraîchage   \n",
       "14195                23               23                            Prairies   \n",
       "384012               11               11  Zones Urbanisées liées à l'habitat   \n",
       "6673                 23               23                            Prairies   \n",
       "102491               31               31                      Bois et forêts   \n",
       "213891               11               11  Zones Urbanisées liées à l'habitat   \n",
       "\n",
       "                                ocs_2_2016  ocs_code_2_2004  ocs_code_1_1949  \\\n",
       "10100                             Prairies               21                2   \n",
       "369875                            Prairies               23                2   \n",
       "380484                      Terres arables               21                2   \n",
       "246996  Zones Urbanisées liées à l'habitat               11                2   \n",
       "283578  Zones Urbanisées liées à l'habitat               22                2   \n",
       "14195                             Prairies               23                2   \n",
       "384012  Zones Urbanisées liées à l'habitat               11                3   \n",
       "6673                              Prairies               23                2   \n",
       "102491                      Bois et forêts               31                3   \n",
       "213891  Zones Urbanisées liées à l'habitat               11                2   \n",
       "\n",
       "                             ocs_1_1949                          ocs_2_2020  \\\n",
       "10100             Territoires agricoles                            Prairies   \n",
       "369875            Territoires agricoles                            Prairies   \n",
       "380484            Territoires agricoles                      Bois et forêts   \n",
       "246996            Territoires agricoles  Zones Urbanisées liées à l'habitat   \n",
       "283578            Territoires agricoles  Zones Urbanisées liées à l'habitat   \n",
       "14195             Territoires agricoles                            Prairies   \n",
       "384012  Forêts et milieux semi naturels  Zones Urbanisées liées à l'habitat   \n",
       "6673              Territoires agricoles                            Prairies   \n",
       "102491  Forêts et milieux semi naturels                      Bois et forêts   \n",
       "213891            Territoires agricoles  Zones Urbanisées liées à l'habitat   \n",
       "\n",
       "                                    ocs_3_2012 ocs_code_3_2020  \\\n",
       "10100                                 Prairies             231   \n",
       "369875                                Prairies             231   \n",
       "380484  Cultures annuelles et terres labourées             312   \n",
       "246996                            Hameau rural             112   \n",
       "283578                   Habitat pavillonnaire             114   \n",
       "14195                                 Prairies             231   \n",
       "384012                   Habitat pavillonnaire             114   \n",
       "6673                                  Prairies             231   \n",
       "102491                          Bois et forêts             311   \n",
       "213891                            Hameau rural             112   \n",
       "\n",
       "                             ocs_1_2012  ocs_code_2_2020  \\\n",
       "10100             Territoires agricoles               23   \n",
       "369875            Territoires agricoles               23   \n",
       "380484            Territoires agricoles               31   \n",
       "246996       Territoires artificialisés               11   \n",
       "283578       Territoires artificialisés               11   \n",
       "14195             Territoires agricoles               23   \n",
       "384012       Territoires artificialisés               11   \n",
       "6673              Territoires agricoles               23   \n",
       "102491  Forêts et milieux semi naturels               31   \n",
       "213891       Territoires artificialisés               11   \n",
       "\n",
       "                                                 geometry  SECTEUR  \\\n",
       "10100   POLYGON ((355433.859 6725027.737, 355364.408 6...  HORS_NM   \n",
       "369875  POLYGON ((322334.256 6714600.310, 322514.540 6...  HORS_NM   \n",
       "380484  POLYGON ((369616.953 6722384.551, 369537.719 6...  HORS_NM   \n",
       "246996  POLYGON ((377978.341 6676903.254, 377970.136 6...  HORS_NM   \n",
       "283578  POLYGON ((357754.613 6681606.781, 357751.447 6...       NM   \n",
       "14195   POLYGON ((354865.204 6713235.322, 354844.574 6...  HORS_NM   \n",
       "384012  POLYGON ((343613.007 6690399.642, 343604.392 6...       NM   \n",
       "6673    POLYGON ((360677.249 6730906.045, 360680.982 6...  HORS_NM   \n",
       "102491  POLYGON ((344121.222 6682704.427, 344121.469 6...       NM   \n",
       "213891  MULTIPOLYGON (((372491.875 6668534.827, 372485...  HORS_NM   \n",
       "\n",
       "                area  \n",
       "10100   44780.659751  \n",
       "369875  40065.409911  \n",
       "380484      0.094406  \n",
       "246996      0.050482  \n",
       "283578   3081.119196  \n",
       "14195    8984.990985  \n",
       "384012   2537.327937  \n",
       "6673     1956.407206  \n",
       "102491    901.543248  \n",
       "213891  11521.504251  "
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Ajout d'un variable indiquant la surface, en m2, de chaque parcelle :\n",
    "mos=mos.to_crs('EPSG:2154')\n",
    "mos['area']=mos['geometry'].area\n",
    "\n",
    "# Transformation du code de niveau en chaîne de caractères : \n",
    "mos['ocs_code_3_1999']=mos['ocs_code_3_1999'].astype('str')\n",
    "mos['ocs_code_3_2004']=mos['ocs_code_3_2004'].astype('str')\n",
    "mos['ocs_code_3_2009']=mos['ocs_code_3_2009'].astype('str')\n",
    "mos['ocs_code_3_2012']=mos['ocs_code_3_2012'].astype('str')\n",
    "mos['ocs_code_3_2016']=mos['ocs_code_3_2016'].astype('str')\n",
    "mos['ocs_code_3_2020']=mos['ocs_code_3_2020'].astype('str')\n",
    "\n",
    "# Regroupements des modes d'occupations du sols en catégories :\n",
    "artif=('111','112','113','114','115','121','122','123','124','131','132','133','141','142','143','144','152','153') \n",
    "# Dans les versions 1 et 2 de ce notebook, les catégories 143 et 144 (Chantier, terrains en mutation, Terrains vacants et friches urbaines) avaient été exclus : trop de faux positifs\n",
    "habitat=('111','112','113','114','115')\n",
    "transp=['131','132','133']\n",
    "activ_commerc=('121','122')\n",
    "services_urbains=('123')\n",
    "agric=('211','221','222','223','224','231')\n",
    "forets_landes_friches_broussailles=('311','321','322')\n",
    "zones_humides=('411','412','413','414','415')\n",
    "rivieres=('511','512')\n",
    "# Attention, ici on va au-delà du périmètre habituel des \"espaces naturels, agricoles et forestiers\", en incluant les parcs et jardins urbains (catégorie 151) et les zones humides ('411','412','413','414','415')\n",
    "vert=['151','211','221','222','223','224','231','311','321','322','411','412','413','414','415','511','512']\n",
    "\n",
    "mos.sample(10)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "##### Artificialisation en Loire-Atlantique de 1999 à 2020"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Sélection, à Nantes Métropole, des espaces agricoles ou naturels en début de période devenus articiels :\n",
    "mos_NM=mos[mos['SECTEUR']=='NM']\n",
    "perte_vert_2020_NM=mos_NM[mos_NM['ocs_code_3_1999'].str.contains('|'.join(vert)) & mos_NM['ocs_code_3_2020'].str.contains('|'.join(artif))]\n",
    "perte_vert_2020_NM=perte_vert_2020_NM[['ocs_code_3_1999','ocs_3_1999','ocs_code_3_2020','ocs_3_2020','area','insee_com','commune','SECTEUR','geometry']]\n",
    "perte_vert_2020_NM=perte_vert_2020_NM[perte_vert_2020_NM['area']>=100]\n",
    "\n",
    "# Export d'un tableur \"origine destination\" des espaces artificialisés\" de Nantes Métropole:\n",
    "bilan_croise_artif_NM_2020=perte_vert_2020_NM.pivot_table(values='area', index=['ocs_3_1999','ocs_3_2020'], aggfunc='sum')\n",
    "bilan_croise_artif_NM_2020.to_csv('bilan_croise_artif_NM_1999-2020.csv', sep=\";\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entre 1999 et 2020, 87.12 km2 ont été artificialisés dans l'ensemble de l'aire d'attraction de Nantes\n"
     ]
    }
   ],
   "source": [
    "# Sélection des espaces agricoles ou naturels en début de période devenus articiels sur l'ensemble de l'aire d'attraction :\n",
    "perte_vert_2020=mos[mos['ocs_code_3_1999'].str.contains('|'.join(vert)) & mos['ocs_code_3_2020'].str.contains('|'.join(artif))]\n",
    "perte_vert_2020=perte_vert_2020[['ocs_code_3_1999','ocs_3_1999','ocs_code_3_2020','ocs_3_2020','area','insee_com','commune','SECTEUR','geometry']]\n",
    "perte_vert_2020=perte_vert_2020[perte_vert_2020['area']>=100]\n",
    "\n",
    "# Reprojection en WGS84 pour un usage web:\n",
    "perte_vert_2020=perte_vert_2020.to_crs('EPSG:4326')\n",
    "#perte_vert_2020.to_file('perte_vert_1999-2020.geojson')\n",
    "\n",
    "# Calcul de la surface totale artificialisée : \n",
    "area=(perte_vert_2020['area'].sum())/1000000\n",
    "area=area.round(decimals=2)\n",
    "area=area.astype(\"str\")\n",
    "print(\"Entre 1999 et 2020, \"+area+\" km2 ont été artificialisés dans l'ensemble de l'aire d'attraction de Nantes\")\n",
    "\n",
    "# Export d'un tableur des pertes d'espaces agricoles ou naturels par commune et par type d'espaces :\n",
    "bilan_perte_vert_2020=perte_vert_2020.pivot_table(values='area', index=['commune','ocs_3_1999'], aggfunc='sum')\n",
    "#bilan_perte_vert_2020.to_csv('bilan_perte_vert_1999-2020.csv')\n",
    "# Idem par secteurs (Métropole et reste de l'aire d'attraction) :\n",
    "bilan_perte_vert_SECTEURS_2020=perte_vert_2020.pivot_table(values='area', index=['SECTEUR','ocs_3_1999'], aggfunc='sum')\n",
    "#bilan_perte_vert_SECTEURS_2020.to_csv('SECTEURS_bilan_perte_vert_1999-2020.csv')\n",
    "\n",
    "# Export d'un tableur des gains d'espaces artificiels par commune et par usage :\n",
    "bilan_artif_2020=perte_vert_2020.pivot_table(values='area', index=['commune','ocs_3_2020'], aggfunc='sum')\n",
    "#bilan_artif_2020.to_csv('bilan_artif_1999-2020.csv')\n",
    "# Idem par secteurs (Métropole et reste de l'aire d'attraction) :\n",
    "bilan_artif_SECTEURS_2020=perte_vert_2020.pivot_table(values='area', index=['SECTEUR','ocs_3_2020'], aggfunc='sum')\n",
    "bilan_artif_SECTEURS_2020.to_csv('SECTEURS_bilan_artif_1999-2020.csv')\n",
    "\n",
    "# Export d'un tableur \"origine destination\" des espaces artificialisés\" :\n",
    "bilan_croise_artif_2020=perte_vert_2020.pivot_table(values='area', index=['ocs_3_1999','ocs_3_2020'], aggfunc='sum')\n",
    "##bilan_croise_artif_2020.to_csv('bilan_croise_artif_1999-2020.csv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Extraction des zones artificialisées pour chaque millésime :\n",
    "\n",
    "artif_1999=mos[mos['ocs_code_3_1999'].str.contains('|'.join(artif))]\n",
    "artif_1999=artif_1999\n",
    "artif_1999=artif_1999.to_crs('EPSG:4326')\n",
    "artif_1999=artif_1999[['commune','ocs_code_3_1999','ocs_3_1999','geometry','area','insee_com','SECTEUR']]\n",
    "#artif_1999.to_file('BDMOS_44_1999_artif.geojson')\n",
    "\n",
    "artif_2004=mos[mos['ocs_code_3_2004'].str.contains('|'.join(artif))]\n",
    "artif_2004=artif_2004\n",
    "artif_2004=artif_2004.to_crs('EPSG:4326')\n",
    "artif_2004=artif_2004[['commune','ocs_code_3_2004','ocs_3_2004','geometry','area','insee_com','SECTEUR']]\n",
    "#artif_2004.to_file('BDMOS_44_2004_artif.geojson')\n",
    "\n",
    "artif_2009=mos[mos['ocs_code_3_2009'].str.contains('|'.join(artif))]\n",
    "artif_2009=artif_2009\n",
    "artif_2009=artif_2009.to_crs('EPSG:4326')\n",
    "artif_2009=artif_2009[['commune','ocs_code_3_2009','ocs_3_2009','geometry','area','insee_com','SECTEUR']]\n",
    "#artif_2009.to_file('BDMOS_44_2009_artif.geojson')\n",
    "\n",
    "artif_2012=mos[mos['ocs_code_3_2012'].str.contains('|'.join(artif))]\n",
    "artif_2012=artif_2012\n",
    "artif_2012=artif_2012.to_crs('EPSG:4326')\n",
    "artif_2012=artif_2012[['commune','ocs_code_3_2012','ocs_3_2012','geometry','area','insee_com','SECTEUR']]\n",
    "#artif_2012.to_file('BDMOS_44_2012_artif.geojson')\n",
    "\n",
    "artif_2016=mos[mos['ocs_code_3_2016'].str.contains('|'.join(artif))]\n",
    "artif_2016=artif_2016\n",
    "artif_2016=artif_2016.to_crs('EPSG:4326')\n",
    "artif_2016=artif_2016[['commune','ocs_code_3_2016','ocs_3_2016','geometry','area','insee_com','SECTEUR']]\n",
    "#artif_2016.to_file('BDMOS_44_2016_artif.geojson')\n",
    "\n",
    "artif_2020=mos[mos['ocs_code_3_2020'].str.contains('|'.join(artif))]\n",
    "artif_2020=artif_2020\n",
    "artif_2020=artif_2020.to_crs('EPSG:4326')\n",
    "artif_2020=artif_2020[['commune','ocs_code_3_2020','ocs_3_2020','geometry','area','insee_com','SECTEUR']]\n",
    "#artif_2020.to_file('BDMOS_44_2020_artif.geojson')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Extraction des zones agricoles ou naturelles en 1999 :\n",
    "\n",
    "agric_1999=mos[mos['ocs_code_3_1999'].str.contains('|'.join(agric))]\n",
    "agric_1999=agric_1999.to_crs('EPSG:4326')\n",
    "#agric_1999.to_file('BDMOS_44_1999_agric.geojson')\n",
    "\n",
    "forets_landes_friches_broussailles_1999=mos[mos['ocs_code_3_1999'].str.contains('|'.join(forets_landes_friches_broussailles))]\n",
    "forets_landes_friches_broussailles_1999=forets_landes_friches_broussailles_1999.to_crs('EPSG:4326')\n",
    "##forets_landes_friches_broussailles_1999.to_file('BDMOS_44_1999_forets_landes_friches_broussailles.geojson')\n",
    "\n",
    "zones_humides_1999=mos[mos['ocs_code_3_1999'].str.contains('|'.join(zones_humides))]\n",
    "zones_humides_1999=zones_humides_1999.to_crs('EPSG:4326')\n",
    "zones_humides_1999.to_file('BDMOS_44_1999_zones_humides.geojson')\n",
    "\n",
    "vert_1999=mos[mos['ocs_code_3_1999'].str.contains('|'.join(vert))]\n",
    "vert_1999=vert_1999.to_crs('EPSG:4326')\n",
    "#vert_1999.to_file('BDMOS_44_1999_vert.geojson')\n",
    "vert_1999_plat=vert_1999.pivot_table(values='area', index=['commune','ocs_3_1999'], aggfunc='sum')\n",
    "#vert_1999_plat.to_csv(\"BDMOS_44_1999_vert.csv\")\n",
    "\n",
    "# Extraction des zones déjà artificialisées en 1999 :\n",
    "gris_1999=mos[mos['ocs_code_3_1999'].str.contains('|'.join(artif))]\n",
    "gris_1999=gris_1999.to_crs('EPSG:4326')\n",
    "#gris_1999.to_file('gris_1999.geojson')"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "##### Artificialisation entre 1999 et 2004"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entre 1999 et 2004,  29.18  km2 ont été artificialisés\n",
      "(3849, 9)\n"
     ]
    }
   ],
   "source": [
    "# Sélection des espaces agricoles ou naturels en début de période devenus articiels :\n",
    "perte_vert_2004=mos[mos['ocs_code_3_1999'].str.contains('|'.join(vert)) & mos['ocs_code_3_2004'].str.contains('|'.join(artif))]\n",
    "perte_vert_2004=perte_vert_2004[['ocs_code_3_1999','ocs_3_1999','ocs_code_3_2004','ocs_3_2004','area','insee_com','commune','SECTEUR','geometry']]\n",
    "perte_vert_2004=perte_vert_2004[perte_vert_2004['area']>=100]\n",
    "\n",
    "# Reprojection en WGS84 pour un usage web:\n",
    "perte_vert_2004=perte_vert_2004.to_crs('EPSG:4326')\n",
    "#perte_vert_2004.to_file('perte_vert_1999-2004.geojson')\n",
    "\n",
    "# Calcul de la surface totale artificialisée : \n",
    "area=(perte_vert_2004['area'].sum())/1000000\n",
    "area=area.round(decimals=2)\n",
    "print('Entre 1999 et 2004, ', area,' km2 ont été artificialisés')\n",
    "print(perte_vert_2004.shape)\n",
    "\n",
    "# Export d'un tableur des pertes d'espaces agricoles ou naturels par commune et par type d'espaces :\n",
    "bilan_perte_vert_2004=perte_vert_2004.pivot_table(values='area', index=['commune','ocs_3_1999'], aggfunc='sum')\n",
    "#bilan_perte_vert_2004.to_csv('bilan_perte_vert_1999-2004.csv')\n",
    "\n",
    "# Export d'un tableur des gains d'espaces artificiels par commune et par usage :\n",
    "bilan_artif_2004=perte_vert_2004.pivot_table(values='area', index=['commune','ocs_3_2004'], aggfunc='sum')\n",
    "#bilan_artif_2004.to_csv('bilan_artif_1999-2004.csv')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": []
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "##### Artificialisation entre 2004 et 2009"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entre 2004 et 2009,  31.72  km2 ont été artificialisés\n",
      "(6512, 9)\n"
     ]
    }
   ],
   "source": [
    "# Sélection des espaces agricoles ou naturels en début de période devenus artificiels :\n",
    "perte_vert_2009=mos[mos['ocs_code_3_2004'].str.contains('|'.join(vert)) & mos['ocs_code_3_2009'].str.contains('|'.join(artif))]\n",
    "perte_vert_2009=perte_vert_2009[['ocs_code_3_2004','ocs_3_2004','ocs_code_3_2009','ocs_3_2009','area','insee_com','commune','SECTEUR','geometry']]\n",
    "perte_vert_2009=perte_vert_2009[perte_vert_2009['area']>=100]\n",
    "\n",
    "# Reprojection en WGS84 pour un usage web:\n",
    "perte_vert_2009=perte_vert_2009.to_crs('EPSG:4326')\n",
    "#perte_vert_2009.to_file('perte_vert_2004-2009.geojson')\n",
    "\n",
    "# Calcul de la surface totale artificialisée : \n",
    "area=(perte_vert_2009['area'].sum())/1000000\n",
    "area=area.round(decimals=2)\n",
    "print('Entre 2004 et 2009, ', area,' km2 ont été artificialisés')\n",
    "print(perte_vert_2009.shape)\n",
    "\n",
    "# Export d'un tableur des pertes d'espaces agricoles ou naturels par commune et par type d'espaces :\n",
    "bilan_perte_vert_2009=perte_vert_2009.pivot_table(values='area', index=['commune','ocs_3_2004'], aggfunc='sum')\n",
    "#bilan_perte_vert_2009.to_csv('bilan_perte_vert_2004-2009.csv')\n",
    "\n",
    "# Export d'un tableur des gains d'espaces artificiels par commune et par usage :\n",
    "bilan_artif_2009=perte_vert_2009.pivot_table(values='area', index=['commune','ocs_3_2009'], aggfunc='sum')\n",
    "#bilan_artif_2009.to_csv('bilan_artif_2004-2009.csv')"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "##### Artificialisation entre 2009 et 2012"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entre 2009 et 2012,  14.15  km2 ont été artificialisés\n",
      "(3465, 9)\n"
     ]
    }
   ],
   "source": [
    "# Sélection des espaces agricoles ou naturels en début de période devenus articiels :\n",
    "perte_vert_2012=mos[mos['ocs_code_3_2009'].str.contains('|'.join(vert)) & mos['ocs_code_3_2012'].str.contains('|'.join(artif))]\n",
    "perte_vert_2012=perte_vert_2012[['ocs_code_3_2009','ocs_3_2009','ocs_code_3_2012','ocs_3_2012','area','insee_com','commune','SECTEUR','geometry']]\n",
    "perte_vert_2012=perte_vert_2012[perte_vert_2012['area']>=100]\n",
    "\n",
    "# Reprojection en WGS84 pour un usage web:\n",
    "perte_vert_2012=perte_vert_2012.to_crs('EPSG:4326')\n",
    "#perte_vert_2012.to_file('perte_vert_2009-2012.geojson')\n",
    "\n",
    "# Calcul de la surface totale artificialisée : \n",
    "area=(perte_vert_2012['area'].sum())/1000000\n",
    "area=area.round(decimals=2)\n",
    "print('Entre 2009 et 2012, ', area,' km2 ont été artificialisés')\n",
    "print(perte_vert_2012.shape)\n",
    "\n",
    "# Export d'un tableur des pertes d'espaces agricoles ou naturels par commune et par type d'espaces :\n",
    "bilan_perte_vert_2012=perte_vert_2012.pivot_table(values='area', index=['commune','ocs_3_2009'], aggfunc='sum')\n",
    "#bilan_perte_vert_2012.to_csv('bilan_perte_vert_2009-2012.csv')\n",
    "\n",
    "# Export d'un tableur des gains d'espaces artificiels par commune et par usage :\n",
    "bilan_artif_2012=perte_vert_2012.pivot_table(values='area', index=['commune','ocs_3_2012'], aggfunc='sum')\n",
    "#bilan_artif_2012.to_csv('bilan_artif_2009-2012.csv')"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "##### Artificialisation entre 2012 et 2016"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entre 2012 et 2016,  12.1  km2 ont été artificialisés\n",
      "(2716, 9)\n"
     ]
    }
   ],
   "source": [
    "# Sélection des espaces agricoles ou naturels en début de période devenus articiels :\n",
    "perte_vert_2016=mos[mos['ocs_code_3_2012'].str.contains('|'.join(vert)) & mos['ocs_code_3_2016'].str.contains('|'.join(artif))]\n",
    "perte_vert_2016=perte_vert_2016[['ocs_code_3_2012','ocs_3_2012','ocs_code_3_2016','ocs_3_2016','area','insee_com','commune','SECTEUR','geometry']]\n",
    "perte_vert_2016=perte_vert_2016[perte_vert_2016['area']>=100]\n",
    "\n",
    "# Reprojection en WGS84 pour un usage web:\n",
    "perte_vert_2016=perte_vert_2016.to_crs('EPSG:4326')\n",
    "#perte_vert_2016.to_file('perte_vert_2012-2016.geojson')\n",
    "\n",
    "# Calcul de la surface totale artificialisée : \n",
    "area=(perte_vert_2016['area'].sum())/1000000\n",
    "area=area.round(decimals=2)\n",
    "print('Entre 2012 et 2016, ', area,' km2 ont été artificialisés')\n",
    "print(perte_vert_2016.shape)\n",
    "\n",
    "# Export d'un tableur des pertes d'espaces agricoles ou naturels par commune et par type d'espaces :\n",
    "bilan_perte_vert_2016=perte_vert_2016.pivot_table(values='area', index=['commune','ocs_3_2012'], aggfunc='sum')\n",
    "#bilan_perte_vert_2016.to_csv('bilan_perte_vert_2012-2016.csv')\n",
    "\n",
    "# Export d'un tableur des gains d'espaces artificiels par commune et par usage :\n",
    "bilan_artif_2016=perte_vert_2016.pivot_table(values='area', index=['commune','ocs_3_2016'], aggfunc='sum')\n",
    "#bilan_artif_2016.to_csv('bilan_artif_2012-2016.csv')"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "##### Artificialisation entre 2016 et 2020"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entre 2016 et 2020,  8.95  km2 ont été artificialisés\n",
      "(1733, 9)\n"
     ]
    }
   ],
   "source": [
    "# Sélection des espaces agricoles ou naturels en début de période devenus articiels :\n",
    "perte_vert_2020=mos[mos['ocs_code_3_2016'].str.contains('|'.join(vert)) & mos['ocs_code_3_2020'].str.contains('|'.join(artif))]\n",
    "perte_vert_2020=perte_vert_2020[['ocs_code_3_2016','ocs_3_2016','ocs_code_3_2020','ocs_3_2020','area','insee_com','commune','SECTEUR','geometry']]\n",
    "perte_vert_2020=perte_vert_2020[perte_vert_2020['area']>=100]\n",
    "\n",
    "# Reprojection en WGS84 pour un usage web:\n",
    "perte_vert_2020=perte_vert_2020.to_crs('EPSG:4326')\n",
    "#perte_vert_2020.to_file('perte_vert_2016-2020.geojson')\n",
    "\n",
    "# Calcul de la surface totale artificialisée : \n",
    "area=(perte_vert_2020['area'].sum())/1000000\n",
    "area=area.round(decimals=2)\n",
    "print('Entre 2016 et 2020, ', area,' km2 ont été artificialisés')\n",
    "print(perte_vert_2020.shape)\n",
    "\n",
    "# Export d'un tableur des pertes d'espaces agricoles ou naturels par commune et par type d'espaces :\n",
    "bilan_perte_vert_2020=perte_vert_2020.pivot_table(values='area', index=['commune','ocs_3_2016'], aggfunc='sum')\n",
    "#bilan_perte_vert_2020.to_csv('bilan_perte_vert_2016-2020.csv')\n",
    "\n",
    "# Export d'un tableur des gains d'espaces artificiels par commune et par usage :\n",
    "bilan_artif_2020=perte_vert_2020.pivot_table(values='area', index=['commune','ocs_3_2020'], aggfunc='sum')\n",
    "#bilan_artif_2020.to_csv('bilan_artif_2016-2020.csv')"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Synthèse par secteurs (Nantes Métropole / reste de l'aire d'attraction)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "###### 1999-2004 ##########\n",
    "# Export d'un tableur des pertes d'espaces agricoles ou naturels par secteur et par type d'espaces :\n",
    "bilan_perte_vert_SECTEURS_2004=perte_vert_2004.pivot_table(values='area', index=['SECTEUR','ocs_3_1999'], aggfunc='sum')\n",
    "#bilan_perte_vert_SECTEURS_2004.to_csv('SECTEUR_bilan_perte_vert_1999-2004.csv')\n",
    "# Export d'un tableur des gains d'espaces artificiels par secteur et par usage :\n",
    "bilan_artif_SECTEURS_2004=perte_vert_2004.pivot_table(values='area', index=['SECTEUR','ocs_3_2004'], aggfunc='sum')\n",
    "#bilan_artif_SECTEURS_2004.to_csv('SECTEUR_bilan_artif_1999-2004.csv')\n",
    "\n",
    "###### 2004-2009 ##########\n",
    "# Export d'un tableur des pertes d'espaces agricoles ou naturels par secteur et par type d'espaces :\n",
    "bilan_perte_vert_SECTEURS_2009=perte_vert_2009.pivot_table(values='area', index=['SECTEUR','ocs_3_2004'], aggfunc='sum')\n",
    "#bilan_perte_vert_SECTEURS_2009.to_csv('SECTEUR_bilan_perte_vert_2004-2009.csv')\n",
    "# Export d'un tableur des gains d'espaces artificiels par secteur et par usage :\n",
    "bilan_artif_SECTEURS_2009=perte_vert_2009.pivot_table(values='area', index=['SECTEUR','ocs_3_2009'], aggfunc='sum')\n",
    "#bilan_artif_SECTEURS_2009.to_csv('SECTEUR_bilan_artif_2004-2009.csv')\n",
    "\n",
    "###### 2009-2012 ##########\n",
    "# Export d'un tableur des pertes d'espaces agricoles ou naturels par secteur et par type d'espaces :\n",
    "bilan_perte_vert_SECTEURS_2012=perte_vert_2012.pivot_table(values='area', index=['SECTEUR','ocs_3_2009'], aggfunc='sum')\n",
    "#bilan_perte_vert_SECTEURS_2012.to_csv('SECTEUR_bilan_perte_vert_2009-2012.csv')\n",
    "# Export d'un tableur des gains d'espaces artificiels par secteur et par usage :\n",
    "bilan_artif_SECTEURS_2012=perte_vert_2012.pivot_table(values='area', index=['SECTEUR','ocs_3_2012'], aggfunc='sum')\n",
    "#bilan_artif_SECTEURS_2012.to_csv('SECTEUR_bilan_artif_2009-2012.csv')\n",
    "\n",
    "###### 2012-2016 ##########\n",
    "# Export d'un tableur des pertes d'espaces agricoles ou naturels par secteur et par type d'espaces :\n",
    "bilan_perte_vert_SECTEURS_2016=perte_vert_2016.pivot_table(values='area', index=['SECTEUR','ocs_3_2012'], aggfunc='sum')\n",
    "#bilan_perte_vert_SECTEURS_2016.to_csv('SECTEUR_bilan_perte_vert_2012-2016.csv')\n",
    "# Export d'un tableur des gains d'espaces artificiels par secteur et par usage :\n",
    "bilan_artif_SECTEURS_2016=perte_vert_2016.pivot_table(values='area', index=['SECTEUR','ocs_3_2016'], aggfunc='sum')\n",
    "#bilan_artif_SECTEURS_2016.to_csv('SECTEUR_bilan_artif_2012-2016.csv')\n",
    "\n",
    "###### 2016-2020 ##########\n",
    "# Export d'un tableur des pertes d'espaces agricoles ou naturels par secteur et par type d'espaces :\n",
    "bilan_perte_vert_SECTEURS_2020=perte_vert_2020.pivot_table(values='area', index=['SECTEUR','ocs_3_2016'], aggfunc='sum')\n",
    "#bilan_perte_vert_SECTEURS_2020.to_csv('SECTEUR_bilan_perte_vert_2016-2020.csv')\n",
    "# Export d'un tableur des gains d'espaces artificiels par secteur et par usage :\n",
    "bilan_artif_SECTEURS_2020=perte_vert_2020.pivot_table(values='area', index=['SECTEUR','ocs_3_2020'], aggfunc='sum')\n",
    "#bilan_artif_SECTEURS_2020.to_csv('SECTEUR_bilan_artif_2016-2020.csv')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3.9.9 ('DATAXPLORE')",
   "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.9.9"
  },
  "orig_nbformat": 4,
  "vscode": {
   "interpreter": {
    "hash": "375f3d611b5238bc734cda3c649c29caee2849da8d5fb7814243c33fcb6d771d"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}