{
 "cells": [
  {
   "cell_type": "raw",
   "metadata": {
    "raw_mimetype": "text/restructuredtext"
    },
   "source": [
    "Conservative Regridding\n",
    "=======================\n",
    "\n",
    "This tutorial will explain to you the concept of conservative regridding and \n",
    "how emiproc uses geopandas to perform this operation.\n",
    "\n",
    "This tutorial is there solely for educational purposes.\n",
    "If you use emiproc, you can use the\n",
    ":func:`emiproc.regrid.remap_inventory`\n",
    "function to perform conservative regridding on your inventories.\n",
    "\n",
    "A good additional source of information on regridding is this article in German [uba_arcgis_2016]_ by the Umweltbundesamt (UBA), which describes a similar approach using ArcGIS."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import geopandas as gpd\n",
    "from emiproc.grids import RegularGrid\n",
    "from shapely.geometry import Polygon"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
       "&lt;html&gt;\n",
       "&lt;head&gt;\n",
       "    \n",
       "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
       "    \n",
       "        &lt;script&gt;\n",
       "            L_NO_TOUCH = false;\n",
       "            L_DISABLE_3D = false;\n",
       "        &lt;/script&gt;\n",
       "    \n",
       "    &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;\n",
       "    &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://code.jquery.com/jquery-3.7.1.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
       "    \n",
       "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
       "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
       "            &lt;style&gt;\n",
       "                #map_ac1bc564eaa2bc4934dda2a686b270d5 {\n",
       "                    position: relative;\n",
       "                    width: 100.0%;\n",
       "                    height: 100.0%;\n",
       "                    left: 0.0%;\n",
       "                    top: 0.0%;\n",
       "                }\n",
       "                .leaflet-container { font-size: 1rem; }\n",
       "            &lt;/style&gt;\n",
       "        \n",
       "&lt;/head&gt;\n",
       "&lt;body&gt;\n",
       "    \n",
       "    \n",
       "            &lt;div class=&quot;folium-map&quot; id=&quot;map_ac1bc564eaa2bc4934dda2a686b270d5&quot; &gt;&lt;/div&gt;\n",
       "        \n",
       "&lt;/body&gt;\n",
       "&lt;script&gt;\n",
       "    \n",
       "    \n",
       "            var map_ac1bc564eaa2bc4934dda2a686b270d5 = L.map(\n",
       "                &quot;map_ac1bc564eaa2bc4934dda2a686b270d5&quot;,\n",
       "                {\n",
       "                    center: [1.0, 1.5],\n",
       "                    crs: L.CRS.Simple,\n",
       "                    zoom: 10,\n",
       "                    zoomControl: true,\n",
       "                    preferCanvas: false,\n",
       "                }\n",
       "            );\n",
       "            L.control.scale().addTo(map_ac1bc564eaa2bc4934dda2a686b270d5);\n",
       "\n",
       "            \n",
       "\n",
       "        \n",
       "    \n",
       "            map_ac1bc564eaa2bc4934dda2a686b270d5.fitBounds(\n",
       "                [[0.0, 0.0], [2.0, 3.0]],\n",
       "                {}\n",
       "            );\n",
       "        \n",
       "    \n",
       "        function geo_json_5ac5f55d8f1e537c3b8ab417a5eab169_styler(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.5, &quot;weight&quot;: 2};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_5ac5f55d8f1e537c3b8ab417a5eab169_highlighter(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.75};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_5ac5f55d8f1e537c3b8ab417a5eab169_pointToLayer(feature, latlng) {\n",
       "            var opts = {&quot;bubblingMouseEvents&quot;: true, &quot;color&quot;: &quot;#3388ff&quot;, &quot;dashArray&quot;: null, &quot;dashOffset&quot;: null, &quot;fill&quot;: true, &quot;fillColor&quot;: &quot;#3388ff&quot;, &quot;fillOpacity&quot;: 0.2, &quot;fillRule&quot;: &quot;evenodd&quot;, &quot;lineCap&quot;: &quot;round&quot;, &quot;lineJoin&quot;: &quot;round&quot;, &quot;opacity&quot;: 1.0, &quot;radius&quot;: 2, &quot;stroke&quot;: true, &quot;weight&quot;: 3};\n",
       "            \n",
       "            let style = geo_json_5ac5f55d8f1e537c3b8ab417a5eab169_styler(feature)\n",
       "            Object.assign(opts, style)\n",
       "            \n",
       "            return new L.CircleMarker(latlng, opts)\n",
       "        }\n",
       "\n",
       "        function geo_json_5ac5f55d8f1e537c3b8ab417a5eab169_onEachFeature(feature, layer) {\n",
       "            layer.on({\n",
       "                mouseout: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                            geo_json_5ac5f55d8f1e537c3b8ab417a5eab169.resetStyle(e.target);\n",
       "                    }\n",
       "                },\n",
       "                mouseover: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                        const highlightStyle = geo_json_5ac5f55d8f1e537c3b8ab417a5eab169_highlighter(e.target.feature)\n",
       "                        e.target.setStyle(highlightStyle);\n",
       "                    }\n",
       "                },\n",
       "            });\n",
       "        };\n",
       "        var geo_json_5ac5f55d8f1e537c3b8ab417a5eab169 = L.geoJson(null, {\n",
       "                onEachFeature: geo_json_5ac5f55d8f1e537c3b8ab417a5eab169_onEachFeature,\n",
       "            \n",
       "                style: geo_json_5ac5f55d8f1e537c3b8ab417a5eab169_styler,\n",
       "                pointToLayer: geo_json_5ac5f55d8f1e537c3b8ab417a5eab169_pointToLayer,\n",
       "        });\n",
       "\n",
       "        function geo_json_5ac5f55d8f1e537c3b8ab417a5eab169_add (data) {\n",
       "            geo_json_5ac5f55d8f1e537c3b8ab417a5eab169\n",
       "                .addData(data);\n",
       "        }\n",
       "            geo_json_5ac5f55d8f1e537c3b8ab417a5eab169_add({&quot;bbox&quot;: [0.0, 0.0, 3.0, 2.0], &quot;features&quot;: [{&quot;bbox&quot;: [0.0, 0.0, 1.0, 1.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[0.0, 0.0], [0.0, 1.0], [1.0, 1.0], [1.0, 0.0], [0.0, 0.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;0&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [0.0, 1.0, 1.0, 2.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[0.0, 1.0], [0.0, 2.0], [1.0, 2.0], [1.0, 1.0], [0.0, 1.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;1&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [1.0, 0.0, 2.0, 1.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[1.0, 0.0], [1.0, 1.0], [2.0, 1.0], [2.0, 0.0], [1.0, 0.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;2&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [1.0, 1.0, 2.0, 2.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[1.0, 1.0], [1.0, 2.0], [2.0, 2.0], [2.0, 1.0], [1.0, 1.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;3&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [2.0, 0.0, 3.0, 1.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[2.0, 0.0], [2.0, 1.0], [3.0, 1.0], [3.0, 0.0], [2.0, 0.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;4&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [2.0, 1.0, 3.0, 2.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[2.0, 1.0], [2.0, 2.0], [3.0, 2.0], [3.0, 1.0], [2.0, 1.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;5&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}], &quot;type&quot;: &quot;FeatureCollection&quot;});\n",
       "\n",
       "        \n",
       "    \n",
       "            geo_json_5ac5f55d8f1e537c3b8ab417a5eab169.addTo(map_ac1bc564eaa2bc4934dda2a686b270d5);\n",
       "        \n",
       "&lt;/script&gt;\n",
       "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x1fe137833b0>"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "\n",
    "\n",
    "# Create some toy data \n",
    "grid = RegularGrid(nx=3, ny=2, dx=1, dy=1, xmin=0, ymin=0, crs=None)\n",
    "grid_serie = grid.gdf.geometry\n",
    "grid_serie.explore()\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This is a regular grid with 6 squared cells.\n",
    "We now want to create another grid in this to see how we can to the remapping."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
       "&lt;html&gt;\n",
       "&lt;head&gt;\n",
       "    \n",
       "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
       "    \n",
       "        &lt;script&gt;\n",
       "            L_NO_TOUCH = false;\n",
       "            L_DISABLE_3D = false;\n",
       "        &lt;/script&gt;\n",
       "    \n",
       "    &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;\n",
       "    &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://code.jquery.com/jquery-3.7.1.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
       "    \n",
       "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
       "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
       "            &lt;style&gt;\n",
       "                #map_a86ae684d62bdecf9775f150ef6a6858 {\n",
       "                    position: relative;\n",
       "                    width: 100.0%;\n",
       "                    height: 100.0%;\n",
       "                    left: 0.0%;\n",
       "                    top: 0.0%;\n",
       "                }\n",
       "                .leaflet-container { font-size: 1rem; }\n",
       "            &lt;/style&gt;\n",
       "        \n",
       "&lt;/head&gt;\n",
       "&lt;body&gt;\n",
       "    \n",
       "    \n",
       "            &lt;div class=&quot;folium-map&quot; id=&quot;map_a86ae684d62bdecf9775f150ef6a6858&quot; &gt;&lt;/div&gt;\n",
       "        \n",
       "&lt;/body&gt;\n",
       "&lt;script&gt;\n",
       "    \n",
       "    \n",
       "            var map_a86ae684d62bdecf9775f150ef6a6858 = L.map(\n",
       "                &quot;map_a86ae684d62bdecf9775f150ef6a6858&quot;,\n",
       "                {\n",
       "                    center: [1.0, 1.5],\n",
       "                    crs: L.CRS.Simple,\n",
       "                    zoom: 10,\n",
       "                    zoomControl: true,\n",
       "                    preferCanvas: false,\n",
       "                }\n",
       "            );\n",
       "            L.control.scale().addTo(map_a86ae684d62bdecf9775f150ef6a6858);\n",
       "\n",
       "            \n",
       "\n",
       "        \n",
       "    \n",
       "            map_a86ae684d62bdecf9775f150ef6a6858.fitBounds(\n",
       "                [[0.0, -1.0], [2.0, 4.0]],\n",
       "                {}\n",
       "            );\n",
       "        \n",
       "    \n",
       "        function geo_json_2192712679ccd98738b737569b4e51ef_styler(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.5, &quot;weight&quot;: 2};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_2192712679ccd98738b737569b4e51ef_highlighter(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.75};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_2192712679ccd98738b737569b4e51ef_pointToLayer(feature, latlng) {\n",
       "            var opts = {&quot;bubblingMouseEvents&quot;: true, &quot;color&quot;: &quot;#3388ff&quot;, &quot;dashArray&quot;: null, &quot;dashOffset&quot;: null, &quot;fill&quot;: true, &quot;fillColor&quot;: &quot;#3388ff&quot;, &quot;fillOpacity&quot;: 0.2, &quot;fillRule&quot;: &quot;evenodd&quot;, &quot;lineCap&quot;: &quot;round&quot;, &quot;lineJoin&quot;: &quot;round&quot;, &quot;opacity&quot;: 1.0, &quot;radius&quot;: 2, &quot;stroke&quot;: true, &quot;weight&quot;: 3};\n",
       "            \n",
       "            let style = geo_json_2192712679ccd98738b737569b4e51ef_styler(feature)\n",
       "            Object.assign(opts, style)\n",
       "            \n",
       "            return new L.CircleMarker(latlng, opts)\n",
       "        }\n",
       "\n",
       "        function geo_json_2192712679ccd98738b737569b4e51ef_onEachFeature(feature, layer) {\n",
       "            layer.on({\n",
       "                mouseout: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                            geo_json_2192712679ccd98738b737569b4e51ef.resetStyle(e.target);\n",
       "                    }\n",
       "                },\n",
       "                mouseover: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                        const highlightStyle = geo_json_2192712679ccd98738b737569b4e51ef_highlighter(e.target.feature)\n",
       "                        e.target.setStyle(highlightStyle);\n",
       "                    }\n",
       "                },\n",
       "            });\n",
       "        };\n",
       "        var geo_json_2192712679ccd98738b737569b4e51ef = L.geoJson(null, {\n",
       "                onEachFeature: geo_json_2192712679ccd98738b737569b4e51ef_onEachFeature,\n",
       "            \n",
       "                style: geo_json_2192712679ccd98738b737569b4e51ef_styler,\n",
       "                pointToLayer: geo_json_2192712679ccd98738b737569b4e51ef_pointToLayer,\n",
       "        });\n",
       "\n",
       "        function geo_json_2192712679ccd98738b737569b4e51ef_add (data) {\n",
       "            geo_json_2192712679ccd98738b737569b4e51ef\n",
       "                .addData(data);\n",
       "        }\n",
       "            geo_json_2192712679ccd98738b737569b4e51ef_add({&quot;bbox&quot;: [-1.0, 0.0, 4.0, 2.0], &quot;features&quot;: [{&quot;bbox&quot;: [-1.0, 0.0, 1.5, 2.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[-1.0, 1.0], [1.5, 0.0], [1.5, 2.0], [-1.0, 1.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;0&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [1.5, 0.0, 4.0, 2.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[1.5, 2.0], [1.5, 0.0], [3.0, 0.0], [4.0, 1.0], [3.0, 2.0], [1.5, 2.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;1&quot;, &quot;properties&quot;: {}, &quot;type&quot;: &quot;Feature&quot;}], &quot;type&quot;: &quot;FeatureCollection&quot;});\n",
       "\n",
       "        \n",
       "    \n",
       "            geo_json_2192712679ccd98738b737569b4e51ef.addTo(map_a86ae684d62bdecf9775f150ef6a6858);\n",
       "        \n",
       "&lt;/script&gt;\n",
       "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x1fe14c86db0>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# We \n",
    "triangle = Polygon([(-1, 1), (1.5, 0), (1.5, 2)])\n",
    "# We put another polygon on the side \n",
    "polygon = Polygon([(1.5, 2 ), (1.5, 0), (3, 0), (4, 1), (3, 2)])\n",
    "serie = gpd.GeoSeries([triangle, polygon])\n",
    "serie.explore()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Spatial join \n",
    "\n",
    "The first step is to genereate find which shapes will intersect shapes from the other grid.\n",
    "\n",
    "This can be performed using the `geopandas.overlay` function.\n",
    "It supports only geodataframes."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "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>source_index</th>\n",
       "      <th>target_index</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((0.00000 1.00000, 1.00000 1.00000, 1....</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((0.00000 1.40000, 1.00000 1.80000, 1....</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((1.00000 1.00000, 1.50000 1.00000, 1....</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 1.00000, 2.00000 0.00000, 1....</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((1.00000 1.80000, 1.50000 2.00000, 1....</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 2.00000, 2.00000 1.00000, 1....</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>4</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 1.00000, 3.00000 1.00000, 3....</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>5</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 2.00000, 3.00000 2.00000, 3....</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   source_index  target_index  \\\n",
       "0             0             0   \n",
       "1             1             0   \n",
       "2             2             0   \n",
       "3             2             1   \n",
       "4             3             0   \n",
       "5             3             1   \n",
       "6             4             1   \n",
       "7             5             1   \n",
       "\n",
       "                                            geometry  \n",
       "0  POLYGON ((0.00000 1.00000, 1.00000 1.00000, 1....  \n",
       "1  POLYGON ((0.00000 1.40000, 1.00000 1.80000, 1....  \n",
       "2  POLYGON ((1.00000 1.00000, 1.50000 1.00000, 1....  \n",
       "3  POLYGON ((2.00000 1.00000, 2.00000 0.00000, 1....  \n",
       "4  POLYGON ((1.00000 1.80000, 1.50000 2.00000, 1....  \n",
       "5  POLYGON ((2.00000 2.00000, 2.00000 1.00000, 1....  \n",
       "6  POLYGON ((2.00000 1.00000, 3.00000 1.00000, 3....  \n",
       "7  POLYGON ((2.00000 2.00000, 3.00000 2.00000, 3....  "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "grid_gdf = grid.gdf\n",
    "grid_gdf['source_index'] = grid_gdf.index\n",
    "\n",
    "gdf_out = gpd.GeoDataFrame(geometry=serie)\n",
    "gdf_out['target_index'] = gdf_out.index\n",
    "\n",
    "gdf_overlayed = gpd.overlay(grid_gdf, gdf_out, how='intersection')\n",
    "gdf_overlayed"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we have the intersections between the two grids. Each line will be a mapping between a cell in the source grid and a cell in the target grid.\n",
    "\n",
    "Cell index of the input grid is `source_index`. The column `target_index` is the index of the cell in the target grid.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<Axes: >"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAF2CAYAAABNisPlAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAroklEQVR4nO3dfXCU9b3//1dSknhYN2AFNyEKUhBowYZbcR0xSIQJtZ54Q6M4DlGwNeIcQaw3eMqB0E5S7M9Ah4M6RQ05dg5KTw8R6yGBHMHqkEBNTQC/kRsNirlZwGASzT18fn942LomgVwhyXVl9/mY+Uzdaz/XJ+/rmu21L965NhsmyQgAAMDBwu0uAAAA4EIILAAAwPEILAAAwPEILAAAwPEILAAAwPEILAAAwPEILAAAwPEILAAAwPEG2F1ATxk2bJjq6+vtLgMAAFjgdrtVWVl5wXlBEViGDRumiooKu8sAAADdEBcXd8HQEhSB5VxnJS4uji4LAAD9hNvtVkVFRZfeu4MisJxTX19PYAEAIAhx0y0AAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8S4Hl6aef1r59+1RXVyefz6etW7dqzJgxF9xv3rx5KisrU2Njo/bv36+5c+e2m5Oenq7Kyko1NDRo586dGj16tJXSAABAELMUWBISErRhwwZdf/31mj17tiIiIrRjxw4NHDiw0328Xq82b96sl19+WZMmTVJubq5yc3M1fvx4/5wnn3xSjz76qNLS0jR9+nR9/fXXys/PV1RUVPePDAAABBXT3TFkyBBjjDEzZszodM5rr71m3nzzzYBthYWF5oUXXvA/rqysNI8//rj/cXR0tGlsbDR33313l+pwu93GGGPcbne3j4XBYDAYDEbfDivv3xd1D8ugQYMkSTU1NZ3O8Xq9KigoCNiWn58vr9crSRo5cqRiY2MD5tTV1Wnv3r3+Od8VGRkpt9sdMAAEp7DICA0YepndZQCwWbe//DAsLEzr1q3Te++9pw8//LDTeTExMfL5fAHbfD6fYmJi/M+f29bZnO9avny5Vq1a1d3SAfQDEVdeoehbpunSGRP1Pdc/2V1Ov3IkIcfuEhCEvhd7xNaf3+3AsmHDBk2YMEE33nhjT9bTJZmZmcrKyvI/Pvf11AD6t7DICLmun6DoW6bpkjHD7S4HgIN0K7CsX79eP/3pT3XTTTddMChUV1fL4/EEbPN4PKqurvY//91t5x6XlJR0uGZLS4taWlq6UzoAB6KbAuBCLN/Dsn79et1xxx2aNWuWjh07dsH5hYWFSkxMDNg2e/ZsFRYWSpLKy8tVVVUVMMftdmv69On+OQCCT1hkhC69aZKGrf6Frvr/HtWgJC9hBUCnLHVYNmzYoHvvvVfJycmqr6/3d05qa2vV1NQkScrJyVFFRYWeeeYZSdLvf/97vfPOO1q2bJneeust3XPPPZo6dap+8Ytf+Nddt26dfvWrX+nIkSMqLy/Xr3/9a1VWVio3N7eHDhOAU9BNAdAdlgLL4sWLJUnvvPNOwPb7779fOTnf3OQ1fPhwnT171v9cYWGh7r33Xv3mN79RRkaGjhw5ottvvz3gRt1nn31WLpdLf/jDHzR48GC99957SkpKUnNzc7cPDIBzcG8KgIsVpm8+39yvud1u1dXVKTo6WvX19XaXA+D/0E2xB58SQm/ojU8JWXn/7vanhACgI3RTAPQGAguAHkE3BUBvIrAA6Da6KQD6CoEFgGV0UwD0NQILgC6hmwLATgQWAOdFNwWAExBYALRDNwWA0xBYAPjRTQHgVAQWIMTRTQHQHxBYgBBFNwVAf0JgAUII3RQA/RWBBQgBdFMA9HcEFiBI0U0BEEwILECQoZsCIBgRWIAgQDcFQLAjsAD9GN0UAKGCwAL0M3RTAIQiAgvQT9BNARDKCCyAg9FNAYBvEFgAB6KbAgCBCCyAQ9BNAYDOEVgAm9FNAYALI7AANggPj5Bryg/lvu16uikA0AUEFqAPuVweDYu7Tp7YSfr0Bp8ixlxld0kA0C8QWIBeFh4eoSs812pY3DQNGny13eUAQL9EYAF6ybe7KRER3JsCABcj3OoOM2bM0LZt21RRUSFjjJKTk887Pzs7W8aYduPgwYP+OStXrmz3fFlZmfWjAWwWHh6hmNjJmjz1IV3nXaorh99AWAGAHmC5w+JyuVRaWqpXXnlFW7duveD8JUuW6Omnn/7HDxwwQKWlpfrTn/4UMO/gwYO65ZZb/I/b2tqslgbYhm4KAPQuy4ElLy9PeXl5XZ5fV1enuro6/+Pk5GRddtllys7ODpjX1tYmn89ntRzANtybAgB9p8/vYVm0aJEKCgr02WefBWy/5pprVFFRoaamJhUWFmr58uU6fvx4h2tERkYqKirK/9jtdvdqzcC30U0BgL7Xp4ElNjZWc+fO1b333huwfe/evbr//vt16NAhxcbGauXKlXr33Xc1YcIEffXVV+3WWb58uVatWtVHVQN0UwDAbn0aWFJTU/Xll18qNzc3YPu3f8V04MAB7d27V59++qlSUlL0yiuvtFsnMzNTWVlZ/sdut1sVFRW9VjdCF90UAHCGPg0sCxcu1KuvvqrW1tbzzqutrdXhw4c1evToDp9vaWlRS0tLb5QI0E0BAAfqs8CSkJCga665Ri+//PIF57pcLo0aNUqvvvpqH1QGfINuCgA4V7c+1vztzsfIkSMVHx+vmpoaHT9+XBkZGYqLi1NqamrAfosWLVJRUZE+/PDDdmv+7ne/05tvvqlPP/1Uw4YNU3p6us6cOaPNmzd345CArqObAgD9g+XAMnXqVO3evdv/eO3atZKkTZs26YEHHlBsbKyGDw/8Mrfo6GjdddddWrJkSYdrXnnlldq8ebMuv/xynTx5Uu+9956uv/56nTp1ymp5QJfQTQGA/iVMkrG7iIvldrtVV1en6Oho1dfX210OHMpp3ZSj044pIp4vP0TPO5KQY3cJCELfiz3S42taef/mu4QQ9OimAED/R2BBUHJaNwUAcHEILAgqdFMAIDgRWNDv0U0BgOBHYEG/RTcFAEIHgQX9Ct0UAAhNBBb0C3RTACC0EVjgWHRTAADnEFjgOHRTAADfRWCBI9BNAQCcD4EFtqKbAgDoCgIL+hzdFACAVQQW9Bm6KQCA7iKwoFfRTQEA9AQCC3oF3RQAQE8isKDH0E0BAPQWAgsuGt0UAEBvI7CgW+imAAD6EoEFltBNAQDYgcCCC6KbAgCwG4EFnaKbAgBwCgILAtBNAQA4EYEFkuimAACcjcASwuimAAD6CwJLCKKbAgDobwgsIYJuCgCgPyOwBDm6KQCAYBBudYcZM2Zo27ZtqqiokDFGycnJ552fkJAgY0y74fF4AuYtXrxY5eXlamxsVFFRkaZNm2a1NPyf8PAIxcRO1uSpD+k671JdOfwGwgoAoF+z3GFxuVwqLS3VK6+8oq1bt3Z5vzFjxqiurs7/+MSJE/7/TklJUVZWltLS0rR3714tXbpU+fn5Gjt2rE6ePGm1xJBFNwUAEKwsB5a8vDzl5eVZ/kEnTpxQbW1th88tW7ZMGzdu1KZNmyRJaWlpuvXWW7Vw4UKtWbPG8s8KJdybAgAIBX12D0tJSYmioqJ08OBBrVq1Snv27JEkRUREaMqUKcrMzPTPNcaooKBAXq+3w7UiIyMVFRXlf+x2u3u3eAeimwIACCW9Hliqqqr00EMP6f3331dUVJQefPBB7d69W9OnT9cHH3ygIUOGaMCAAfL5fAH7+Xw+jRs3rsM1ly9frlWrVvV26Y5DNwUAEKp6PbAcPnxYhw8f9j8uLCzUqFGj9Nhjj2nBggXdWjMzM1NZWVn+x263WxUVFRddq1PRTQEAhDpbPta8b98+3XjjjZKkU6dOqa2trd2nhjwej6qrqzvcv6WlRS0tLb1ep53opgAA8A+2BJaJEyeqqqpKktTa2qri4mIlJibqjTfekCSFhYUpMTFR//7v/25HebaimwIAQHvd+ljz6NGj/Y9Hjhyp+Ph41dTU6Pjx48rIyFBcXJxSU1MlSUuWLFF5ebk+/PBDXXLJJXrwwQc1a9YszZkzx79GVlaWcnJy9P7772vfvn1aunSpXC6XsrOze+AQnY9uCgAA52c5sEydOlW7d+/2P167dq0kadOmTXrggQcUGxur4cOH+5+PjIzUc889p7i4ODU0NGj//v265ZZbAtbYsmWLhg4dqtWrVysmJkYlJSVKSkoK+FstwYhuCgAAXRMmydhdxMVyu92qq6tTdHS06uvr7S7nvOim4Jyj044pIv4qu8tAEDqSkGN3CQhC34s90uNrWnn/5ruE+gjdFAAAuo/A0ovopgAA0DMILL2AbgoAAD2LwNJD6KYAANB7CCwXiW4KAAC9j8DSDXRTAADoWwQWC+imAABgDwJLF3hiJyku7jq6KQAA2ITA0gXjfniXwsO/Z3cZAACErHC7CwAAALgQAgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8y4FlxowZ2rZtmyoqKmSMUXJy8nnn33HHHdqxY4dOnDih2tpa7dmzR3PmzAmYs3LlShljAkZZWZnV0gAAQJCyHFhcLpdKS0v1yCOPdGn+TTfdpJ07d+onP/mJpkyZol27dunNN9/UxIkTA+YdPHhQMTEx/nHjjTdaLQ0AAASpAVZ3yMvLU15eXpfnP/bYYwGP//Vf/1XJycm67bbbVFJS4t/e1tYmn89ntRwAABAC+vwelrCwMLndbtXU1ARsv+aaa1RRUaGPP/5Yf/zjH3XVVVd1ukZkZKTcbnfAAAAAwavPA8svf/lLXXrppdqyZYt/2969e3X//fcrKSlJDz/8sEaOHKl3331Xl156aYdrLF++XHV1df5RUVHRV+UDAAAb9GlgmT9/vlauXKmUlBSdPHnSvz0vL0//9V//pQMHDmjHjh36yU9+osGDByslJaXDdTIzMxUdHe0fcXFxfXUIAADABpbvYemuu+++Wy+99JJ+9rOf6X//93/PO7e2tlaHDx/W6NGjO3y+paVFLS0tvVEmAABwoD7psNxzzz3Kzs7W/Pnz9T//8z8XnO9yuTRq1ChVVVX1QXUAAMDpuvWx5vj4eMXHx0uSRo4cqfj4eP9NshkZGcrJyfHPnz9/vv7jP/5Djz/+uPbu3SuPxyOPx6Po6Gj/nN/97ne66aabNGLECHm9Xm3dulVnzpzR5s2bL/b4AABAELAcWKZOnaqSkhL/R5LXrl2rkpISrV69WpIUGxur4cOH++f/4he/UEREhJ5//nlVV1f7x+9//3v/nCuvvFKbN2/WoUOHtGXLFn3xxRe6/vrrderUqYs8PAAAEAws38PyzjvvKCwsrNPnH3jggYDHN9988wXXnD9/vtUyAABACOG7hAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgOMRWAAAgONZDiwzZszQtm3bVFFRIWOMkpOTL7hPQkKCiouL1dTUpCNHjig1NbXdnMWLF6u8vFyNjY0qKirStGnTrJYGAACClOXA4nK5VFpaqkceeaRL86+++mq99dZb2rVrlyZOnKh169bppZde0pw5c/xzUlJSlJWVpfT0dE2ePFmlpaXKz8/X0KFDrZYHAACC0ACrO+Tl5SkvL6/L89PS0lReXq5f/vKXkqSPPvpIN954ox577DHt2LFDkrRs2TJt3LhRmzZt8u9z6623auHChVqzZo3VEgEAQJCxHFis8nq9KigoCNiWn5+vdevWSZIiIiI0ZcoUZWZm+p83xqigoEBer7fDNSMjIxUVFeV/7Ha7e77wbzlzpllnz3K7D3pY6xmdbWy2u4p+Ifyfoi48CUBQ6/XAEhMTI5/PF7DN5/Np0KBBuuSSS3TZZZdpwIABHc4ZN25ch2suX75cq1at6q2S24l850OFhxFY0LMiLxuiqAGX2V0GgtCKH06wuwQEpSO2/vR++S6cmZmp6Oho/4iLi7O7JAAA0It6vcNSXV0tj8cTsM3j8ai2tlZNTU06deqU2traOpxTXV3d4ZotLS1qaWnptZoBAICz9HqHpbCwUImJiQHbZs+ercLCQklSa2uriouLA+aEhYUpMTHRPwcAAIS2bn2sOT4+XvHx8ZKkkSNHKj4+XldddZUkKSMjQzk5Of75L774on7wgx9ozZo1Gjt2rB5++GGlpKRo7dq1/jlZWVn6+c9/rgULFmjcuHF64YUX5HK5lJ2dfbHHBwAAgoDlXwlNnTpVu3fv9j8+Fzw2bdqkBx54QLGxsRo+fLj/+WPHjunWW2/V2rVrtWTJEn3++ed68MEH/R9plqQtW7Zo6NChWr16tWJiYlRSUqKkpCSdOHHiIg4NAAAEizBJxu4iLpbb7VZdXZ2io6NVX1/f4+vP0p18Sgg9rvz2IYoaO8buMhCEbp//rt0lIAhlxG/t8TWtvH/zLgwAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByPwAIAAByvW4Fl8eLFKi8vV2Njo4qKijRt2rRO5+7atUvGmHbjL3/5i39OdnZ2u+e3b9/endIAAEAQGmB1h5SUFGVlZSktLU179+7V0qVLlZ+fr7Fjx+rkyZPt5t95552KjIz0P7788stVWlqqP/3pTwHztm/frgceeMD/uLm52WppAAAgSFnusCxbtkwbN27Upk2bVFZWprS0NDU0NGjhwoUdzj99+rR8Pp9/zJ49Ww0NDe0CS3Nzc8C8L7/8slsHBAAAgo+lwBIREaEpU6aooKDAv80Yo4KCAnm93i6tsWjRIr322mtqaGgI2D5z5kz5fD599NFHev755/X973/fSmkAACCIWfqV0JAhQzRgwAD5fL6A7T6fT+PGjbvg/tOmTdO1116rRYsWBWzPy8vTf//3f6u8vFyjRo1SRkaGtm/fLq/Xq7Nnz7ZbJzIyUlFRUf7HbrfbymEAAIB+xvI9LBdj0aJF2r9/v/72t78FbH/99df9/33w4EHt379fn3zyiWbOnKm333673TrLly/XqlWrertcAADgEJZ+JXTq1Cm1tbXJ4/EEbPd4PKqurj7vvgMHDtQ999yjl19++YI/p7y8XCdPntTo0aM7fD4zM1PR0dH+ERcX1/WDAAAA/Y6lwNLa2qri4mIlJib6t4WFhSkxMVGFhYXn3fdnP/uZoqKi9Mc//vGCPycuLk6XX365qqqqOny+paVF9fX1AQMAAAQvy58SysrK0s9//nMtWLBA48aN0wsvvCCXy6Xs7GxJUk5OjjIyMtrtt2jRIuXm5qqmpiZgu8vl0rPPPqvp06drxIgRmjVrlt544w0dPXpU+fn53TwsAAAQTCzfw7JlyxYNHTpUq1evVkxMjEpKSpSUlKQTJ05IkoYPH97uRtkxY8ZoxowZmj17drv1zpw5ox//+MdKTU3V4MGDVVlZqR07dmjFihVqaWnp5mEBAIBgEibJ2F3ExXK73aqrq1N0dHSv/Hpolu5UeBjfYoCeVX77EEWNHWN3GQhCt89/1+4SEIQy4rf2+JpW3r95FwYAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYOmCD/Suqs1xnTVn7S4FAICQNMDuAvqD0zqp0zqpCEVpmBmhOI3UwDC33WUBABAyCCwWtKpZn+qwPtVhXWaGKk4/0BWKU3gYjSoAAHoTgaWb6LoAANB3CCwXia4LAAC9j8DSg+i6AADQOwgsvYCuCwAAPYvA0svougAAcPEILH2ErgsAAN1HYLEBXRcAAKwhsNiIrgsAAF3TrXfGxYsXq7y8XI2NjSoqKtK0adM6nZuamipjTMBobGxsNy89PV2VlZVqaGjQzp07NXr06O6U1m+d1kkd1F69q7d0xOxXg6m3uyQAABzDcmBJSUlRVlaW0tPTNXnyZJWWlio/P19Dhw7tdJ/a2lrFxMT4x4gRIwKef/LJJ/Xoo48qLS1N06dP19dff638/HxFRUVZP6J+7lzXZY/yVWze4TuMAABQNwLLsmXLtHHjRm3atEllZWVKS0tTQ0ODFi5c2Ok+xhj5fD7/OHHiRMDzS5cu1W9+8xtt27ZNBw4c0IIFCzRs2DDdfvvtlg8omNB1AQDgG5YCS0REhKZMmaKCggL/NmOMCgoK5PV6O93v0ksv1bFjx/TZZ58pNzdXP/rRj/zPjRw5UrGxsQFr1tXVae/evZ2uGRkZKbfbHTCCGV0XAECosxRYhgwZogEDBsjn8wVs9/l8iomJ6XCfQ4cOaeHChUpOTtZ9992n8PBw7dmzR3FxcZLk38/KmsuXL1ddXZ1/VFRUWDmMfo2uCwAgFPX6x1GKior06quvqrS0VH/9619155136uTJk3rooYe6vWZmZqaio6P941z4CSV0XQAAocTSx5pPnTqltrY2eTyegO0ej0fV1dVdWqOtrU0ffPCB/1NA5/b77hoej0clJSUdrtHS0qKWlhYrpQc1/q4LACDYWeqwtLa2qri4WImJif5tYWFhSkxMVGFhYdd+YHi4rr32WlVVVUmSysvLVVVVFbCm2+3W9OnTu7wmvkHXBQAQrCz/4bisrCzl5OTo/fff1759+7R06VK5XC5lZ2dLknJyclRRUaFnnnlGkrRixQoVFRXp6NGjGjx4sJ544gmNGDFCL730kn/NdevW6Ve/+pWOHDmi8vJy/frXv1ZlZaVyc3N75ihDEF0XAEAwsRxYtmzZoqFDh2r16tWKiYlRSUmJkpKS/B9VHj58uM6e/ce/6i+77DJt3LhRMTExOn36tIqLi3XDDTeorKzMP+fZZ5+Vy+XSH/7wBw0ePFjvvfeekpKS1Nzc3AOHGNr4a7oAgGAQJsnYXcTFcrvdqqurU3R0tOrr+dTMhUQoSsNE18Vu5bcPUdTYMXaXgSB0+/x37S4BQSgjfmuPr2nl/ZvvEgpBdF0AAP0NgSXEca8LAKA/ILBAEl0XAICzEVjQDl0XAIDTEFjQKbouAACnILCgS+i6AADsRGCBJXRdAAB2ILCg2+i6AAD6CoEFF42uCwCgtxFY0KPougAAegOBBb2CrgsAoCcRWNDr6LoAAC4WgQV9hq4LAKC7CCywBV0XAIAVBBbYiq4LAKArCCxwDLouAIDOEFjgOHRdAADfRWCBo9F1AQBIBBb0E3RdACC0EVjQ79B1AYDQQ2BBv0XXBQBCB4EFQYGuCwAENwILggpdFwAITgQWBC26LgAQPAgsCHp0XQCg/yOwIKTQdQGA/onAgpBE1wUA+pduXZ0XL16s8vJyNTY2qqioSNOmTet07oMPPqi//vWvqqmpUU1NjXbu3NlufnZ2towxAWP79u3dKQ2w7LRO6qD26l29pSNmvxpMvd0lAQC+w3JgSUlJUVZWltLT0zV58mSVlpYqPz9fQ4cO7XD+zJkztXnzZt18883yer06fvy4duzYoWHDhgXM2759u2JiYvxj/vz53TsioJvOdV32KF/F5h1Vm+M6a87aXRYAQN0ILMuWLdPGjRu1adMmlZWVKS0tTQ0NDVq4cGGH8++77z698MILKi0t1aFDh/Tggw8qPDxciYmJAfOam5vl8/n848svv+zWAQE9ga4LADiLpcASERGhKVOmqKCgwL/NGKOCggJ5vd4urTFw4EBFRESopqYmYPvMmTPl8/n00Ucf6fnnn9f3v//9TteIjIyU2+0OGEBvoOsCAM5gKbAMGTJEAwYMkM/nC9ju8/kUExPTpTXWrFmjysrKgNCTl5enBQsWKDExUU899ZQSEhK0fft2hYd3XN7y5ctVV1fnHxUVFVYOA+gWui4AYJ8+/ZTQU089pXvuuUczZ85Uc3Ozf/vrr7/u/++DBw9q//79+uSTTzRz5ky9/fbb7dbJzMxUVlaW/7Hb7Sa0oM/wCSMA6HuWAsupU6fU1tYmj8cTsN3j8ai6uvq8+z7++ON6+umndcstt+jAgQPnnVteXq6TJ09q9OjRHQaWlpYWtbS0WCkd6BX8XRcA6BuW/knY2tqq4uLigBtmw8LClJiYqMLCwk73e+KJJ7RixQolJSWpuLj4gj8nLi5Ol19+uaqqqqyUB9iGe10AoHdZ/pVQVlaWcnJy9P7772vfvn1aunSpXC6XsrOzJUk5OTmqqKjQM888I0l68skntXr1at177706duyYvzvz1Vdf6euvv5bL5dLKlSv15z//WdXV1Ro1apSeffZZHT16VPn5+T14qEDfoOsCAD3PcmDZsmWLhg4dqtWrVysmJkYlJSVKSkrSiRMnJEnDhw/X2bP/+Jflww8/rKioKP35z38OWGfVqlVKT0/XmTNn9OMf/1ipqakaPHiwKisrtWPHDq1YsYJf+6Bf414XAOg5YZKM3UVcLLfbrbq6OkVHR6u+nk9uwLkiFKVh+qbr4rtjpKLGjrG7JASh2+e/a3cJCEIZ8Vt7fE0r7998lxDQh77ddRm49we67OzXih5zrcK+x/8VAeB8uEoCNmmo/EQNb3yi6oGXavCEaRo88XpFfb/jr7gAgFBHYAFsdqbhK32xb5e+2LdbA0eM0mUTvXRdAOA7uCICjmHU8OlRNXx6lK4LAHwHgQVwILouABCIqx/gaHRdAEAisAD9Bl0XAKGMKx3Q79B1ARB6CCxAP0bXBUCo4KoGBAW6LgCCG4EFCDJ0XQAEI65gQNCi6wIgeBBYgBBA1wVAf8fVCggpdF0A9E8EFiBE0XUB0J9wZQJCHl0XAM5HYAHgR9cFgFNxFQLQAbouAJyFwALgvOi6AHACrjgAuoiuCwD7EFgAWNau6xLvVfRYui4Aeg9XFwAXga4LgL5BYAHQI+i6AOhNXEkA9DC6LgB6HoEFQK+h6wKgp3DVANAH6LoAuDjh3dlp8eLFKi8vV2Njo4qKijRt2rTzzp83b57KysrU2Nio/fv3a+7cue3mpKenq7KyUg0NDdq5c6dGjx7dndIAONy5rsvHf/itjm1+XrX/7wOZM212lwXA4SwHlpSUFGVlZSk9PV2TJ09WaWmp8vPzNXRox/9S8nq92rx5s15++WVNmjRJubm5ys3N1fjx4/1znnzyST366KNKS0vT9OnT9fXXXys/P19RUVHdPzIADvdN16Vi26s6vGG1fG+/qeaak3YXBcChwiQZKzsUFRXpb3/7m/7lX/7lmwXCwnT8+HGtX79ea9asaTf/tddek8vl0m233ebfVlhYqJKSEj388MOSpMrKSj333HN67rnnJEnR0dHy+Xy6//779frrr1+wJrfbrbq6OkVHR6u+vt7K4QBwlDDudekBt89/1+4SEIQy4rf2+JpW3r8tdVgiIiI0ZcoUFRQU+LcZY1RQUCCv19vhPl6vN2C+JOXn5/vnjxw5UrGxsQFz6urqtHfv3k7XjIyMlNvtDhgAggFdFwAds/TPlyFDhmjAgAHy+XwB230+n8aNG9fhPjExMR3Oj4mJ8T9/bltnc75r+fLlWrVqlZXSAfQz3/6E0fcu+SedaWqwu6R+4//91u4KgJ7XrZtu7ZaZmano6Gj/iIuLs7skAL3GEFYAWAssp06dUltbmzweT8B2j8ej6urqDveprq4+7/xz/2tlzZaWFtXX1wcMAAAQvCwFltbWVhUXFysxMdG/LSwsTImJiSosLOxwn8LCwoD5kjR79mz//PLyclVVVQXMcbvdmj59eqdrAgCA0GOsjJSUFNPY2GgWLFhgxo0bZ1588UVTU1NjrrjiCiPJ5OTkmIyMDP98r9drWlpazLJly8zYsWPNypUrTXNzsxk/frx/zpNPPmlqamrMbbfdZiZMmGC2bt1qPv74YxMVFdWlmtxutzHGGLfbbelYGAwGg8Fg2Dcsvn9b/wGPPPKIOXbsmGlqajJFRUXmuuuu8z+3a9cuk52dHTB/3rx55qOPPjJNTU3mwIEDZu7cue3WTE9PN1VVVaaxsdHs3LnTXHPNNb11wAwGg8FgMBwwrLx/W/47LE7E32EBAKD/6bW/wwIAAGAHAgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8S9/W7HRut9vuEgAAQBdZed8OisBy7oArKipsrgQAAFjldrsv+IfjguIv3UrSsGHDeuWv3LrdblVUVCguLo6/onsBnKuu41x1HefKGs5X13Guuq43z5Xb7VZlZeUF5wVFh0VSlw72YtTX1/OC7iLOVddxrrqOc2UN56vrOFdd1xvnqqvrcdMtAABwPAILAABwPALLBTQ3N2vVqlVqbm62uxTH41x1Heeq6zhX1nC+uo5z1XVOOFdBc9MtAAAIXnRYAACA4xFYAACA4xFYAACA4xFYAACA4xFYJC1evFjl5eVqbGxUUVGRpk2bdt758+bNU1lZmRobG7V//37NnTu3jyq1n5VzlZqaKmNMwGhsbOzDau0zY8YMbdu2TRUVFTLGKDk5+YL7JCQkqLi4WE1NTTpy5IhSU1P7oFL7WT1XCQkJ7V5Xxhh5PJ4+qtg+Tz/9tPbt26e6ujr5fD5t3bpVY8aMueB+oXjN6s65CtVrVlpamkpLS1VbW6va2lrt2bNHSUlJ593HjtdUyAeWlJQUZWVlKT09XZMnT1Zpaany8/M1dOjQDud7vV5t3rxZL7/8siZNmqTc3Fzl5uZq/PjxfVx537N6riSptrZWMTEx/jFixIg+rNg+LpdLpaWleuSRR7o0/+qrr9Zbb72lXbt2aeLEiVq3bp1eeuklzZkzp5crtZ/Vc3XOmDFjAl5bJ06c6KUKnSMhIUEbNmzQ9ddfr9mzZysiIkI7duzQwIEDO90nVK9Z3TlXUmhesz7//HM9/fTTmjJliqZOnaq3335bb7zxhn70ox91ON/O15QJ5VFUVGTWr1/vfxwWFmY+//xz89RTT3U4/7XXXjNvvvlmwLbCwkLzwgsv2H4sTjtXqamp5vTp07bXbfcwxpjk5OTzzvntb39rDhw4ELBt8+bNZvv27bbX77RzlZCQYIwxZtCgQbbXa/cYMmSIMcaYGTNmdDonlK9ZVs8V16x/jC+++MIsXLiww+fsek2FdIclIiJCU6ZMUUFBgX+bMUYFBQXyer0d7uP1egPmS1J+fn6n84NFd86VJF166aU6duyYPvvsM+Xm5naa2ENdqL6uLkZJSYkqKyu1Y8cO3XDDDXaXY4tBgwZJkmpqajqdw2vrG105VxLXrPDwcN19991yuVwqLCzscI5dr6mQDixDhgzRgAED5PP5Arb7fD7FxMR0uE9MTIyl+cGiO+fq0KFDWrhwoZKTk3XfffcpPDxce/bsUVxcXF+U3K909roaNGiQLrnkEpuqcqaqqio99NBDuuuuu3TXXXfp+PHj2r17tyZNmmR3aX0qLCxM69at03vvvacPP/yw03mhes36tq6eq1C+Zk2YMEH19fVqbm7Wiy++qDvuuENlZWUdzrXrNRU039YM5ykqKlJRUZH/8Z49e1RWVqaHHnpI//Zv/2ZjZejPDh8+rMOHD/sfFxYWatSoUXrssce0YMECGyvrWxs2bNCECRN044032l2K43X1XIXyNevQoUOaOHGiBg0apHnz5iknJ0cJCQmdhhY7hHSH5dSpU2pra2v36QKPx6Pq6uoO96murrY0P1h051x9V1tbmz744AONHj26N0rs1zp7XdXW1qqpqcmmqvqPffv2hdTrav369frpT3+qm2++WRUVFeedG6rXrHOsnKvvCqVrVmtrqz7++GP9/e9/1zPPPKPS0lItWbKkw7l2vaZCOrC0traquLhYiYmJ/m1hYWFKTEzs9Hd3hYWFAfMlafbs2Z3ODxbdOVffFR4ermuvvVZVVVW9VWa/Faqvq54yceLEkHldrV+/XnfccYdmzZqlY8eOXXB+KL+2rJ6r7wrla1Z4eLiioqI6fM7O15TtdyPbOVJSUkxjY6NZsGCBGTdunHnxxRdNTU2NueKKK4wkk5OTYzIyMvzzvV6vaWlpMcuWLTNjx441K1euNM3NzWb8+PG2H4vTztWKFSvM7NmzzciRI82kSZPMf/7nf5qGhgbzwx/+0PZj6e3hcrlMfHy8iY+PN8YYs3TpUhMfH2+uuuoqI8lkZGSYnJwc//yrr77afPXVV2bNmjVm7Nix5uGHHzatra1mzpw5th+L087VkiVLzD//8z+bUaNGmfHjx5u1a9eatrY2M2vWLNuPpbfHhg0bzOnTp81NN91kPB6Pf1xyySX+OVyzun+uQvWalZGRYWbMmGFGjBhhJkyYYDIyMsyZM2fMLbfc4rTXlP0ny+7xyCOPmGPHjpmmpiZTVFRkrrvuOv9zu3btMtnZ2QHz582bZz766CPT1NRkDhw4YObOnWv7MTjxXGVlZfnnVlVVmb/85S9m4sSJth9DX4xzH739rnPnJzs72+zatavdPn//+99NU1OTOXr0qElNTbX9OJx4rp544glz5MgR09DQYE6dOmXefvttM3PmTNuPoy9GZ779WuGa1f1zFarXrJdeesmUl5ebpqYm4/P5zM6dO/1hxUmvqbD/+w8AAADHCul7WAAAQP9AYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI73/wOfqreD0jfkGgAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAF2CAYAAABNisPlAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAA9hAAAPYQGoP6dpAAArQ0lEQVR4nO3dfXBU9b3H8c8GkjiuG2Ih3Q1RkIJICzY8FtcRoUQYqPVGhcaH6RAFWxGugnh9wFsvhHZCsWOgw6BO1YbUzkXp7SVivRDIFa0OCWjaBHCioA2KeVjAYBLNM/zuH162rkkgJyQ5Z3ffr5nf6J79nV++58zO7ofvnt11STICAABwsBi7CwAAADgfAgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHA8AgsAAHC8gXYX0FuGDh2qhoYGu8sAAAAWeDweVVVVnXdeRASWoUOHqrKy0u4yAABAD6SkpJw3tEREYDnbWUlJSaHLAgBAmPB4PKqsrOzWa3dEBJazGhoaCCwAAEQgLroFAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOZymwPPbYY9q/f7/q6+sVCAS0bds2jR49+rz7zZ8/X+Xl5WpqatKBAwc0d+7cDnOysrJUVVWlxsZG7d69W6NGjbJSGgAAiGCWAsv06dO1adMmXXPNNZo1a5ZiY2O1a9cuXXzxxV3u4/f7tWXLFr3wwguaMGGC8vPzlZ+fr7FjxwbnPPLII3rggQe0ePFiTZ06VV9++aUKCgoUHx/f8yMDAAARxfR0DBkyxBhjzLRp07qc89JLL5lXX301ZFtRUZF55plngrerqqrMQw89FLydkJBgmpqazG233datOjwejzHGGI/H0+NjYTAYDAaD0b/Dyuv3BV3DMmjQIElSbW1tl3P8fr8KCwtDthUUFMjv90uSRowYoeTk5JA59fX12rdvX3DON8XFxcnj8YQMAJEpRgN0kbru4gKIDj3+8UOXy6UNGzbo7bff1nvvvdflPJ/Pp0AgELItEAjI5/MF7z+7ras537Ry5UqtXr26p6UDCANuJegyfUc+DdPr1eV2lwNEvQHJR2z9+z0OLJs2bdK4ceN03XXX9WY93bJ27Vrl5OQEb5/9eWoA4S1GA+TVZUrRd5ToGmx3OQAcpEeBZePGjfrxj3+s66+//rxBoaamRl6vN2Sb1+tVTU1N8P5vbjt7u7S0tNM1W1tb1dra2pPSATjQ17spsa44u8sB4ECWr2HZuHGjbrnlFs2cOVNHjx497/yioiKlpaWFbJs1a5aKiookSRUVFaqurg6Z4/F4NHXq1OAcAJEnRgOUrOGarB/K75qty12jCCsAumSpw7Jp0ybdeeedSk9PV0NDQ7BzUldXp+bmZklSXl6eKisr9fjjj0uSfvvb3+rNN9/UihUr9Nprr+n222/X5MmT9fOf/zy47oYNG/SLX/xCR44cUUVFhX75y1+qqqpK+fn5vXSYAJyCbgqAnrAUWJYsWSJJevPNN0O233XXXcrLy5MkDRs2TGfOnAneV1RUpDvvvFO/+tWvlJ2drSNHjujmm28OuVD3ySeflNvt1u9+9zslJibq7bff1pw5c9TS0tLjAwPgHFybAuBCufTV55vDmsfjUX19vRISEtTQ0GB3OQD+X291UwqqynqxKgA90RefErLy+t3jTwkBQGfopgDoCwQWAL2Ca1MA9CUCC4Aeo5sCoL8QWABYRjcFQH8jsADoFropAOxEYAFwTnRTADgBgQVAB3RTADgNgQVAEN0UAE5FYAGiHN0UAOGAwAJEKbopAMIJgQWIInRTAIQrAgsQBeimAAh3BBYgQtFNARBJCCxAhKGbAiASEViACEA3BUCkI7AAYYxuCoBoQWABwgzdFADRiMAChAm6KQCiGYEFcDC6KQDwFQIL4EB0UwAgFIEFcAi6KQDQNQILYDO6KQBwfgQWwAYxGqAkDdXlGkU3BQC6gcAC9CO3EpSiEUrWcH2pBsIKAHQTgQXoY/+8NmWEEl1D/nmHsa8mAAg3BBagj3y9m8K1KQBwYWKs7jBt2jRt375dlZWVMsYoPT39nPNzc3NljOkwDh06FJyzatWqDveXl5dbPxrAZjEaoGQN12TNkN81W8NcVxJWAKAXWO6wuN1ulZWV6fe//722bdt23vnLli3TY4899s8/OHCgysrK9Kc//Slk3qFDh3TDDTcEb7e3t1stDbAN3RQA6FuWA8vOnTu1c+fObs+vr69XfX198HZ6erouvfRS5ebmhsxrb29XIBCwWg5gmy6vTQEA9Lp+v4Zl0aJFKiws1CeffBKy/corr1RlZaWam5tVVFSklStX6tixY52uERcXp/j4+OBtj8fTpzUDX0c3BQD6X78GluTkZM2dO1d33nlnyPZ9+/bprrvu0gcffKDk5GStWrVKb731lsaNG6cvvviiwzorV67U6tWr+6lqgG4KANitXwNLZmamPv/8c+Xn54ds//pbTAcPHtS+ffv08ccfKyMjQ7///e87rLN27Vrl5OQEb3s8HlVWVvZZ3YhedFMAwBn6NbAsXLhQL774otra2s45r66uTocPH9aoUaM6vb+1tVWtra19USJANwUAHKjfAsv06dN15ZVX6oUXXjjvXLfbrZEjR+rFF1/sh8qAr9BNAQDn6tHHmr/e+RgxYoRSU1NVW1urY8eOKTs7WykpKcrMzAzZb9GiRSouLtZ7773XYc3f/OY3evXVV/Xxxx9r6NChysrK0unTp7Vly5YeHBLQfXRTACA8WA4skydP1htvvBG8vX79eknS5s2bdffddys5OVnDhg0L2SchIUHz5s3TsmXLOl3zsssu05YtWzR48GCdOHFCb7/9tq655hqdPHnSanlAt9BNAYDw4lIE/KKJx+NRfX29EhIS1NDQYHc5cCindVM+N5/x44fdVFBVZncJQNQbkHyk19e08vrNbwkh4tFNAYDwR2BBRHJaNwUAcGEILIgodFMAIDIRWBD26KYAQOQjsCBs0U0BgOhBYEFYoZsCANGJwIKwQDcFAKIbgQWORTcFAHAWgQWOQzcFAPBNBBY4At0UAMC5EFhgK7opAIDuILCg39FNAQBYRWBBv6GbAgDoKQIL+hTdFABAbyCwoE/QTQEA9CYCC3oN3RQAQF8hsOCC0U0BAPQ1Agt6hG4KAKA/EVhgCd0UAIAdCCw4L7opAAC7EVjQJbopAACnILAgBN0UAIATEVggiW4KAMDZCCxRjG4KACBcEFiiEN0UAEC4IbBECbopAIBwRmCJcHRTAACRIMbqDtOmTdP27dtVWVkpY4zS09PPOX/69OkyxnQYXq83ZN6SJUtUUVGhpqYmFRcXa8qUKVZLw/+L0QAla7gma4b8rtka5rqSsAIACGuWA4vb7VZZWZmWLl1qab/Ro0fL5/MFx/Hjx4P3ZWRkKCcnR1lZWZo4caLKyspUUFCgpKQkq+VFNbcSNFqpmqYbNdY1hbd+AAARw/JbQjt37tTOnTst/6Hjx4+rrq6u0/tWrFih5557Tps3b5YkLV68WDfeeKMWLlyodevWWf5b0YRrUwAA0cByh6WnSktLVVVVpV27dunaa68Nbo+NjdWkSZNUWFgY3GaMUWFhofx+f6drxcXFyePxhIxoQzcFABBN+jywVFdX695779W8efM0b948HTt2TG+88YYmTJggSRoyZIgGDhyoQCAQsl8gEJDP5+t0zZUrV6q+vj44Kisr+/owHIFrUwAA0arPPyV0+PBhHT58OHi7qKhII0eO1IMPPqgFCxb0aM21a9cqJycneNvj8UR0aOGTPgCAaGfLx5r379+v6667TpJ08uRJtbe3d/jUkNfrVU1NTaf7t7a2qrW1tc/rtBPXpgAA8E/9dg3L140fP17V1dWSpLa2NpWUlCgtLS14v8vlUlpamoqKiuwoz1ZcmwIAQEeWOyxut1ujRo0K3h4xYoRSU1NVW1urY8eOKTs7WykpKcrMzJQkLVu2TBUVFXrvvfd00UUX6Z577tHMmTM1e/bs4Bo5OTnKy8vTu+++q/3792v58uVyu93Kzc3thUN0PropAACcm+XAMnnyZL3xxhvB2+vXr5ckbd68WXfffbeSk5M1bNiw4P1xcXF66qmnlJKSosbGRh04cEA33HBDyBpbt25VUlKS1qxZI5/Pp9LSUs2ZMyfku1oiEdemAADQPS5Jxu4iLpTH41F9fb0SEhLU0NBgdznnRDcFZ31uPlOia7DdZYSFgqoyu0sAot6A5CO9vqaV129+S6if0E0BAKDnCCx9iG4KAAC9g8DSB+imAADQuwgsvYRuCgAAfYfAcoHopgAA0PcILD1ANwUAgP5FYLGAbgoAAPYgsHRDsoYpRd+hmwIAgE0ILN3wXU1WjMuWn10CAACy6ccPAQAArCCwAAAAxyOwAAAAxyOwAAAAxyOwAAAAxyOwAAAAxyOwAAAAxyOwAAAAxyOwAAAAxyOwAAAAxyOwAAAAxyOwAAAAxyOwAAAAxyOwAAAAxyOwAAAAxyOwAAAAxyOwAAAAx7McWKZNm6bt27ersrJSxhilp6efc/4tt9yiXbt26fjx46qrq9PevXs1e/bskDmrVq2SMSZklJeXWy0NAABEKMuBxe12q6ysTEuXLu3W/Ouvv167d+/Wj370I02aNEl79uzRq6++qvHjx4fMO3TokHw+X3Bcd911VksDAAARaqDVHXbu3KmdO3d2e/6DDz4Ycvvf//3flZ6erptuukmlpaXB7e3t7QoEAlbLAQAAUaDfr2FxuVzyeDyqra0N2X7llVeqsrJSH330kf74xz/q8ssv73KNuLg4eTyekAEAACJXvweWf/u3f9Mll1yirVu3Brft27dPd911l+bMmaP77rtPI0aM0FtvvaVLLrmk0zVWrlyp+vr64KisrOyv8gEAgA36NbDccccdWrVqlTIyMnTixIng9p07d+q//uu/dPDgQe3atUs/+tGPlJiYqIyMjE7XWbt2rRISEoIjJSWlvw4BAADYwPI1LD1122236fnnn9dPfvIT/e///u8559bV1enw4cMaNWpUp/e3traqtbW1L8oEAAAO1C8dlttvv125ubm644479D//8z/nne92uzVy5EhVV1f3Q3UAAMDpevSx5tTUVKWmpkqSRowYodTU1OBFstnZ2crLywvOv+OOO/SHP/xBDz30kPbt2yev1yuv16uEhITgnN/85je6/vrrNXz4cPn9fm3btk2nT5/Wli1bLvT4AABABLAcWCZPnqzS0tLgR5LXr1+v0tJSrVmzRpKUnJysYcOGBef//Oc/V2xsrJ5++mnV1NQEx29/+9vgnMsuu0xbtmzRBx98oK1bt+qzzz7TNddco5MnT17g4QEAgEjgkmTsLuJCeTwe1dfXKyEhQQ0NDb2+/kzdqhgXv2KA3vW5+UyJrsF2lxEWCqrK7C4BiHoDko/0+ppWXr95FQYAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5HYAEAAI5nObBMmzZN27dvV2VlpYwxSk9PP+8+06dPV0lJiZqbm3XkyBFlZmZ2mLNkyRJVVFSoqalJxcXFmjJlitXSAABAhLIcWNxut8rKyrR06dJuzb/iiiv02muvac+ePRo/frw2bNig559/XrNnzw7OycjIUE5OjrKysjRx4kSVlZWpoKBASUlJVssDAAARyCXJ9HRnY4xuvvlmvfLKK13O+fWvf60bb7xRV199dXDbli1blJiYqLlz50qSiouL9c477+j+++//qiiXS8eOHdPGjRu1bt2689bh8XhUX1+vhIQENTQ09PRwujRTtyrGxbtn6F2fm8+U6BpsdxlhoaCqzO4SgKg3IPlIr69p5fV7YK//9W/w+/0qLCwM2VZQUKANGzZIkmJjYzVp0iStXbs2eL8xRoWFhfL7/Z2uGRcXp/j4+OBtj8fT+4V/zWm164xx9enfQPQxMmo3bXaXAQBhoc8Di8/nUyAQCNkWCAQ0aNAgXXTRRbr00ks1cODATueMGTOm0zVXrlyp1atX91XJHQzQQDos6HUu49JAV6zdZQBAWAjLV+G1a9cqISEhOFJSUuwuCQAA9KE+77DU1NTI6/WGbPN6vaqrq1Nzc7NOnjyp9vb2TufU1NR0umZra6taW1v7rGYAAOAsfd5hKSoqUlpaWsi2WbNmqaioSJLU1tamkpKSkDkul0tpaWnBOQAAILr16GPNqampSk1NlSSNGDFCqampuvzyyyVJ2dnZysvLC85/9tln9Z3vfEfr1q3TVVddpfvuu08ZGRlav359cE5OTo5+9rOfacGCBRozZoyeeeYZud1u5ebmXujxAQCACGD5LaHJkyfrjTfeCN4+Gzw2b96su+++W8nJyRo2bFjw/qNHj+rGG2/U+vXrtWzZMn366ae65557tGvXruCcrVu3KikpSWvWrJHP51NpaanmzJmj48ePX8ChAQCASHFB38PiFHwPC8IR38PSfXwPC2A/u7+HhVdhAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeD0KLEuWLFFFRYWamppUXFysKVOmdDl3z549MsZ0GH/5y1+Cc3Jzczvcv2PHjp6UBgAAItBAqztkZGQoJydHixcv1r59+7R8+XIVFBToqquu0okTJzrMv/XWWxUXFxe8PXjwYJWVlelPf/pTyLwdO3bo7rvvDt5uaWmxWhoAAIhQljssK1as0HPPPafNmzervLxcixcvVmNjoxYuXNjp/FOnTikQCATHrFmz1NjY2CGwtLS0hMz7/PPPe3RAAAAg8lgKLLGxsZo0aZIKCwuD24wxKiwslN/v79YaixYt0ksvvaTGxsaQ7TNmzFAgEND777+vp59+Wt/61reslAYAACKYpbeEhgwZooEDByoQCIRsDwQCGjNmzHn3nzJliq6++motWrQoZPvOnTv13//936qoqNDIkSOVnZ2tHTt2yO/368yZMx3WiYuLU3x8fPC2x+OxchgAACDMWL6G5UIsWrRIBw4c0DvvvBOy/eWXXw7+/6FDh3TgwAH94x//0IwZM/T66693WGflypVavXp1X5cLAAAcwtJbQidPnlR7e7u8Xm/Idq/Xq5qamnPue/HFF+v222/XCy+8cN6/U1FRoRMnTmjUqFGd3r927VolJCQER0pKSvcPAgAAhB1LgaWtrU0lJSVKS0sLbnO5XEpLS1NRUdE59/3JT36i+Ph4/fGPfzzv30lJSdHgwYNVXV3d6f2tra1qaGgIGQAAIHJZ/pRQTk6Ofvazn2nBggUaM2aMnnnmGbndbuXm5kqS8vLylJ2d3WG/RYsWKT8/X7W1tSHb3W63nnzySU2dOlXDhw/XzJkz9corr+jDDz9UQUFBDw8LAABEEsvXsGzdulVJSUlas2aNfD6fSktLNWfOHB0/flySNGzYsA4Xyo4ePVrTpk3TrFmzOqx3+vRpff/731dmZqYSExNVVVWlXbt26YknnlBra2sPDwsAAEQSlyRjdxEXyuPxqL6+XgkJCX3y9tBM3aoYF79igN71uflMia7BdpcRFgqqyuwuAYh6A5KP9PqaVl6/eRUGAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2ABAACOR2Dphr/rLdWYYzpjzthdCgAAUWmg3QWEg1M6oVM6oVjFa6gZrhSN0MUuj91lAQAQNQgsFrSpRR/rsD7WYV1qkpSi7+jbSlGMi0YVAAB9icDSQ3RdAADoPwSWC0TXBQCAvkdg6UV0XQAA6BsElj5A1wUAgN5FYOljdF0AALhwBJZ+QtcFAICeI7DYgK4LAADWEFhsRNcFAIDu6dEr45IlS1RRUaGmpiYVFxdrypQpXc7NzMyUMSZkNDU1dZiXlZWlqqoqNTY2avfu3Ro1alRPSgtbp3RCh7RPb+k1HTEH1Gga7C4JAADHsBxYMjIylJOTo6ysLE2cOFFlZWUqKChQUlJSl/vU1dXJ5/MFx/Dhw0Puf+SRR/TAAw9o8eLFmjp1qr788ksVFBQoPj7e+hGFubNdl70qUIl5k98wAgBAPQgsK1as0HPPPafNmzervLxcixcvVmNjoxYuXNjlPsYYBQKB4Dh+/HjI/cuXL9evfvUrbd++XQcPHtSCBQs0dOhQ3XzzzZYPKJLQdQEA4CuWAktsbKwmTZqkwsLC4DZjjAoLC+X3+7vc75JLLtHRo0f1ySefKD8/X9/73veC940YMULJyckha9bX12vfvn1drhkXFyePxxMyIhldFwBAtLMUWIYMGaKBAwcqEAiEbA8EAvL5fJ3u88EHH2jhwoVKT0/XT3/6U8XExGjv3r1KSUmRpOB+VtZcuXKl6uvrg6OystLKYYQ1ui4AgGjU5x9HKS4u1osvvqiysjL99a9/1a233qoTJ07o3nvv7fGaa9euVUJCQnCcDT/RhK4LACCaWPpY88mTJ9Xe3i6v1xuy3ev1qqampltrtLe36+9//3vwU0Bn9/vmGl6vV6WlpZ2u0draqtbWViulRzS+1wUAEOksdVja2tpUUlKitLS04DaXy6W0tDQVFRV17w/GxOjqq69WdXW1JKmiokLV1dUha3o8Hk2dOrXba+IrdF0AAJHK8hfH5eTkKC8vT++++67279+v5cuXy+12Kzc3V5KUl5enyspKPf7445KkJ554QsXFxfrwww+VmJiohx9+WMOHD9fzzz8fXHPDhg36xS9+oSNHjqiiokK//OUvVVVVpfz8/N45yihE1wUAEEksB5atW7cqKSlJa9askc/nU2lpqebMmRP8qPKwYcN05sw//1V/6aWX6rnnnpPP59OpU6dUUlKia6+9VuXl5cE5Tz75pNxut373u98pMTFRb7/9tubMmaOWlpZeOMToxrfpAgAigUuSsbuIC+XxeFRfX6+EhAQ1NPCpmfOJVbyGiq6L3T43nynRNdjuMsJCQVWZ3SUAUW9A8pFeX9PK6ze/JRSF6LoAAMINgSXKca0LACAcEFggia4LAMDZCCzogK4LAMBpCCzoEl0XAIBTEFjQLXRdAAB2IrDAErouAAA7EFjQY3RdAAD9hcCCC0bXBQDQ1wgs6FV0XQAAfYHAgj5B1wUA0JsILOhzdF0AABeKwIJ+Q9cFANBTBBbYgq4LAMAKAgtsRdcFANAdBBY4Bl0XAEBXCCxwHLouAIBvIrDA0ei6AAAkAgvCBF0XAIhuBBaEHbouABB9CCwIW3RdACB6EFgQEei6AEBkI7AgotB1AYDIRGBBxKLrAgCRg8CCiEfXBQDCH4EFUYWuCwCEJwILohJdFwAILz16dl6yZIkqKirU1NSk4uJiTZkypcu599xzj/7617+qtrZWtbW12r17d4f5ubm5MsaEjB07dvSkNMCyUzqhQ9qnt/SajpgDajQNdpcEAPgGy4ElIyNDOTk5ysrK0sSJE1VWVqaCggIlJSV1On/GjBnasmWLfvjDH8rv9+vYsWPatWuXhg4dGjJvx44d8vl8wXHHHXf07IiAHjrbddmrApWYN1VjjumMOWN3WQAASS5JxsoOxcXFeuedd3T//fd/tYDLpWPHjmnjxo1at27defePiYnRqVOn9K//+q968cUXJX3VYUlMTNQtt9xi/QgkeTwe1dfXKyEhQQ0N/OsYvSdW8RqqvrnW5XPzmRJdg3t1zUhVUFVmdwlA1BuQfKTX17Ty+m2pwxIbG6tJkyapsLAwuM0Yo8LCQvn9/m6tcfHFFys2Nla1tbUh22fMmKFAIKD3339fTz/9tL71rW91uUZcXJw8Hk/IAPoCXRcAcAZLgWXIkCEaOHCgAoFAyPZAICCfz9etNdatW6eqqqqQ0LNz504tWLBAaWlpevTRRzV9+nTt2LFDMTGdl7dy5UrV19cHR2VlpZXDAHqEa10AwD79+imhRx99VLfffrtmzJihlpaW4PaXX345+P+HDh3SgQMH9I9//EMzZszQ66+/3mGdtWvXKicnJ3jb4/EQWtBv+IQRAPQ/S4Hl5MmTam9vl9frDdnu9XpVU1Nzzn0feughPfbYY7rhhht08ODBc86tqKjQiRMnNGrUqE4DS2trq1pbW62UDvQJvtcFAPqHpX8StrW1qaSkRGlpacFtLpdLaWlpKioq6nK/hx9+WE888YTmzJmjkpKS8/6dlJQUDR48WNXV1VbKA2zDtS4A0LcsvyWUk5OjvLw8vfvuu9q/f7+WL18ut9ut3NxcSVJeXp4qKyv1+OOPS5IeeeQRrVmzRnfeeaeOHj0a7M588cUX+vLLL+V2u7Vq1Sr9+c9/Vk1NjUaOHKknn3xSH374oQoKCnrxUIH+QdcFAHqf5cCydetWJSUlac2aNfL5fCotLdWcOXN0/PhxSdKwYcN05sw//2V53333KT4+Xn/+859D1lm9erWysrJ0+vRpff/731dmZqYSExNVVVWlXbt26YknnuBtH4Q1rnUBgN5j+XtYnIjvYUG4+Pr3urSqle9h6Sa+hwWwn93fw8JvCQH96Otdl0Ql6TJD1wUAuoPAAtjkc53Q51zrAgDdQmABbMa1LgBwfgQWwEH4hBEAdI7AAjgQXRcACEVgARyOrgsAEFiAsEHXBUA0I7AAYYiuC4BoQ2ABwhhdFwDRgsACRAi6LgAiGYEFiDB0XQBEIgILEMHougCIFAQWIArQdQEQ7ggsQJSh6wIgHBFYgChF1wVAOCGwAKDrAsDxCCwAgui6AHAqAguATtF1AeAkBBYA50TXBYATEFgAdBtdFwB2IbAAsIyuC4D+RmABcEHougDoDwQWAL2CrguAvkRgAdDr6LoA6G0EFgB9hq4LgN5CYAHQL+i6ALgQPfpnzpIlS1RRUaGmpiYVFxdrypQp55w/f/58lZeXq6mpSQcOHNDcuXM7zMnKylJVVZUaGxu1e/dujRo1qielAXC4s12XvSpQiXlTNeaYzpgzdpcFwOEsB5aMjAzl5OQoKytLEydOVFlZmQoKCpSUlNTpfL/fry1btuiFF17QhAkTlJ+fr/z8fI0dOzY455FHHtEDDzygxYsXa+rUqfryyy9VUFCg+Pj4nh8ZAMc7pRM6pH16S6/piDmgRtNgd0kAHMolyVjZobi4WO+8847uv//+rxZwuXTs2DFt3LhR69at6zD/pZdektvt1k033RTcVlRUpNLSUt13332SpKqqKj311FN66qmnJEkJCQkKBAK666679PLLL5+3Jo/Ho/r6eiUkJKihgSc8IJxdqo7XuhRUldlcFYAByUd6fU0rr9+WOiyxsbGaNGmSCgsLg9uMMSosLJTf7+90H7/fHzJfkgoKCoLzR4wYoeTk5JA59fX12rdvX5drxsXFyePxhAwAkYGuC4DOWLrodsiQIRo4cKACgUDI9kAgoDFjxnS6j8/n63S+z+cL3n92W1dzvmnlypVavXq1ldIBhJmvf8LoouQ4tanV7pIA2CgsP1u4du1aJSQkBEdKSordJQHoQ4QVAJYCy8mTJ9Xe3i6v1xuy3ev1qqamptN9ampqzjn/7H+trNna2qqGhoaQAQAAIpelwNLW1qaSkhKlpaUFt7lcLqWlpamoqKjTfYqKikLmS9KsWbOC8ysqKlRdXR0yx+PxaOrUqV2uCQAAoo+xMjIyMkxTU5NZsGCBGTNmjHn22WdNbW2t+fa3v20kmby8PJOdnR2c7/f7TWtrq1mxYoW56qqrzKpVq0xLS4sZO3ZscM4jjzxiamtrzU033WTGjRtntm3bZj766CMTHx/frZo8Ho8xxhiPx2PpWBgMBoPBYNg3LL5+W/8DS5cuNUePHjXNzc2muLjY/OAHPwjet2fPHpObmxsyf/78+eb99983zc3N5uDBg2bu3Lkd1szKyjLV1dWmqanJ7N6921x55ZV9dcAMBoPBYDAcMKy8flv+HhYn4ntYAAAIP332PSwAAAB2ILAAAADHI7AAAADHI7AAAADHI7AAAADHI7AAAADHI7AAAADHs/RrzU7n8XjsLgEAAHSTldftiAgsZw+4srLS5koAAIBVHo/nvF8cFxHfdCtJQ4cO7ZNvufV4PKqsrFRKSgrfonsenKvu41x1H+fKGs5X93Guuq8vz5XH41FVVdV550VEh0VStw72QjQ0NPCA7ibOVfdxrrqPc2UN56v7OFfd1xfnqrvrcdEtAABwPAILAABwPALLebS0tGj16tVqaWmxuxTH41x1H+eq+zhX1nC+uo9z1X1OOFcRc9EtAACIXHRYAACA4xFYAACA4xFYAACA4xFYAACA4xFYJC1ZskQVFRVqampScXGxpkyZcs758+fPV3l5uZqamnTgwAHNnTu3nyq1n5VzlZmZKWNMyGhqaurHau0zbdo0bd++XZWVlTLGKD09/bz7TJ8+XSUlJWpubtaRI0eUmZnZD5Xaz+q5mj59eofHlTFGXq+3nyq2z2OPPab9+/ervr5egUBA27Zt0+jRo8+7XzQ+Z/XkXEXrc9bixYtVVlamuro61dXVae/evZozZ84597HjMRX1gSUjI0M5OTnKysrSxIkTVVZWpoKCAiUlJXU63+/3a8uWLXrhhRc0YcIE5efnKz8/X2PHju3nyvuf1XMlSXV1dfL5fMExfPjwfqzYPm63W2VlZVq6dGm35l9xxRV67bXXtGfPHo0fP14bNmzQ888/r9mzZ/dxpfazeq7OGj16dMhj6/jx431UoXNMnz5dmzZt0jXXXKNZs2YpNjZWu3bt0sUXX9zlPtH6nNWTcyVF53PWp59+qscee0yTJk3S5MmT9frrr+uVV17R9773vU7n2/mYMtE8iouLzcaNG4O3XS6X+fTTT82jjz7a6fyXXnrJvPrqqyHbioqKzDPPPGP7sTjtXGVmZppTp07ZXrfdwxhj0tPTzznn17/+tTl48GDIti1btpgdO3bYXr/TztX06dONMcYMGjTI9nrtHkOGDDHGGDNt2rQu50Tzc5bVc8Vz1j/HZ599ZhYuXNjpfXY9pqK6wxIbG6tJkyapsLAwuM0Yo8LCQvn9/k738fv9IfMlqaCgoMv5kaIn50qSLrnkEh09elSffPKJ8vPzu0zs0S5aH1cXorS0VFVVVdq1a5euvfZau8uxxaBBgyRJtbW1Xc7hsfWV7pwrieesmJgY3XbbbXK73SoqKup0jl2PqagOLEOGDNHAgQMVCARCtgcCAfl8vk738fl8luZHip6cqw8++EALFy5Uenq6fvrTnyomJkZ79+5VSkpKf5QcVrp6XA0aNEgXXXSRTVU5U3V1te69917NmzdP8+bN07Fjx/TGG29owoQJdpfWr1wulzZs2KC3335b7733XpfzovU56+u6e66i+Tlr3LhxamhoUEtLi5599lndcsstKi8v73SuXY+piPm1ZjhPcXGxiouLg7f37t2r8vJy3XvvvfqP//gPGytDODt8+LAOHz4cvF1UVKSRI0fqwQcf1IIFC2ysrH9t2rRJ48aN03XXXWd3KY7X3XMVzc9ZH3zwgcaPH69BgwZp/vz5ysvL0/Tp07sMLXaI6g7LyZMn1d7e3uHTBV6vVzU1NZ3uU1NTY2l+pOjJufqm9vZ2/f3vf9eoUaP6osSw1tXjqq6uTs3NzTZVFT72798fVY+rjRs36sc//rF++MMfqrKy8pxzo/U56ywr5+qbouk5q62tTR999JH+9re/6fHHH1dZWZmWLVvW6Vy7HlNRHVja2tpUUlKitLS04DaXy6W0tLQu37srKioKmS9Js2bN6nJ+pOjJufqmmJgYXX311aquru6rMsNWtD6uesv48eOj5nG1ceNG3XLLLZo5c6aOHj163vnR/Niyeq6+KZqfs2JiYhQfH9/pfXY+pmy/GtnOkZGRYZqamsyCBQvMmDFjzLPPPmtqa2vNt7/9bSPJ5OXlmezs7OB8v99vWltbzYoVK8xVV11lVq1aZVpaWszYsWNtPxannasnnnjCzJo1y4wYMcJMmDDB/Od//qdpbGw03/3ud20/lr4ebrfbpKammtTUVGOMMcuXLzepqanm8ssvN5JMdna2ycvLC86/4oorzBdffGHWrVtnrrrqKnPfffeZtrY2M3v2bNuPxWnnatmyZeZf/uVfzMiRI83YsWPN+vXrTXt7u5k5c6btx9LXY9OmTebUqVPm+uuvN16vNzguuuii4Byes3p+rqL1OSs7O9tMmzbNDB8+3IwbN85kZ2eb06dPmxtuuMFpjyn7T5bdY+nSpebo0aOmubnZFBcXmx/84AfB+/bs2WNyc3ND5s+fP9+8//77prm52Rw8eNDMnTvX9mNw4rnKyckJzq2urjZ/+ctfzPjx420/hv4YZz96+01nz09ubq7Zs2dPh33+9re/mebmZvPhhx+azMxM24/Diefq4YcfNkeOHDGNjY3m5MmT5vXXXzczZsyw/Tj6Y3Tl648VnrN6fq6i9Tnr+eefNxUVFaa5udkEAgGze/fuYFhx0mPK9f//AwAA4FhRfQ0LAAAIDwQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeAQWAADgeP8HfdZavpxlpW0AAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gdf_overlayed.plot('source_index')\n",
    "gdf_overlayed.plot('target_index')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Weight calculation\n",
    "\n",
    "The next step is to calculate the weights for each intersection.\n",
    "\n",
    "For that we apply a method that compares the area of the polygon that falls in the target polygon.\n",
    "\n",
    "To calculate the area of each intersection, we need to create the intersection polygon."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "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>source_index</th>\n",
       "      <th>target_index</th>\n",
       "      <th>geometry</th>\n",
       "      <th>area_overlay</th>\n",
       "      <th>area_source</th>\n",
       "      <th>weights</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((0.00000 1.00000, 1.00000 1.00000, 1....</td>\n",
       "      <td>0.60</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.60</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((0.00000 1.40000, 1.00000 1.80000, 1....</td>\n",
       "      <td>0.60</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.60</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((1.00000 1.00000, 1.50000 1.00000, 1....</td>\n",
       "      <td>0.45</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.45</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 1.00000, 2.00000 0.00000, 1....</td>\n",
       "      <td>0.50</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.50</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((1.00000 1.80000, 1.50000 2.00000, 1....</td>\n",
       "      <td>0.45</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.45</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 2.00000, 2.00000 1.00000, 1....</td>\n",
       "      <td>0.50</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.50</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>4</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 1.00000, 3.00000 1.00000, 3....</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>5</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 2.00000, 3.00000 2.00000, 3....</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.00</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   source_index  target_index  \\\n",
       "0             0             0   \n",
       "1             1             0   \n",
       "2             2             0   \n",
       "3             2             1   \n",
       "4             3             0   \n",
       "5             3             1   \n",
       "6             4             1   \n",
       "7             5             1   \n",
       "\n",
       "                                            geometry  area_overlay  \\\n",
       "0  POLYGON ((0.00000 1.00000, 1.00000 1.00000, 1....          0.60   \n",
       "1  POLYGON ((0.00000 1.40000, 1.00000 1.80000, 1....          0.60   \n",
       "2  POLYGON ((1.00000 1.00000, 1.50000 1.00000, 1....          0.45   \n",
       "3  POLYGON ((2.00000 1.00000, 2.00000 0.00000, 1....          0.50   \n",
       "4  POLYGON ((1.00000 1.80000, 1.50000 2.00000, 1....          0.45   \n",
       "5  POLYGON ((2.00000 2.00000, 2.00000 1.00000, 1....          0.50   \n",
       "6  POLYGON ((2.00000 1.00000, 3.00000 1.00000, 3....          1.00   \n",
       "7  POLYGON ((2.00000 2.00000, 3.00000 2.00000, 3....          1.00   \n",
       "\n",
       "   area_source  weights  \n",
       "0          1.0     0.60  \n",
       "1          1.0     0.60  \n",
       "2          1.0     0.45  \n",
       "3          1.0     0.50  \n",
       "4          1.0     0.45  \n",
       "5          1.0     0.50  \n",
       "6          1.0     1.00  \n",
       "7          1.0     1.00  "
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_overlayed['area_overlay'] = gdf_overlayed.area\n",
    "gdf_overlayed['area_source'] = gdf_overlayed['source_index'].map(grid_gdf.area)\n",
    "gdf_overlayed['weights'] = gdf_overlayed['area_overlay'] / gdf_overlayed['area_source']\n",
    "gdf_overlayed"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The weight express how much of the source cell is in the target cell.\n",
    "\n",
    "If you sum the weights for all intersections of a source cell, you will get 1."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
       "&lt;html&gt;\n",
       "&lt;head&gt;\n",
       "    \n",
       "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
       "    \n",
       "        &lt;script&gt;\n",
       "            L_NO_TOUCH = false;\n",
       "            L_DISABLE_3D = false;\n",
       "        &lt;/script&gt;\n",
       "    \n",
       "    &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;\n",
       "    &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://code.jquery.com/jquery-3.7.1.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
       "    \n",
       "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
       "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
       "            &lt;style&gt;\n",
       "                #map_b39b247aec3d6b4a73e76904f4ffdf5f {\n",
       "                    position: relative;\n",
       "                    width: 100.0%;\n",
       "                    height: 100.0%;\n",
       "                    left: 0.0%;\n",
       "                    top: 0.0%;\n",
       "                }\n",
       "                .leaflet-container { font-size: 1rem; }\n",
       "            &lt;/style&gt;\n",
       "        \n",
       "    \n",
       "                    &lt;style&gt;\n",
       "                        .foliumtooltip {\n",
       "                            \n",
       "                        }\n",
       "                       .foliumtooltip table{\n",
       "                            margin: auto;\n",
       "                        }\n",
       "                        .foliumtooltip tr{\n",
       "                            text-align: left;\n",
       "                        }\n",
       "                        .foliumtooltip th{\n",
       "                            padding: 2px; padding-right: 8px;\n",
       "                        }\n",
       "                    &lt;/style&gt;\n",
       "            \n",
       "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/d3/3.5.5/d3.min.js&quot;&gt;&lt;/script&gt;\n",
       "&lt;/head&gt;\n",
       "&lt;body&gt;\n",
       "    \n",
       "    \n",
       "            &lt;div class=&quot;folium-map&quot; id=&quot;map_b39b247aec3d6b4a73e76904f4ffdf5f&quot; &gt;&lt;/div&gt;\n",
       "        \n",
       "&lt;/body&gt;\n",
       "&lt;script&gt;\n",
       "    \n",
       "    \n",
       "            var map_b39b247aec3d6b4a73e76904f4ffdf5f = L.map(\n",
       "                &quot;map_b39b247aec3d6b4a73e76904f4ffdf5f&quot;,\n",
       "                {\n",
       "                    center: [1.0, 1.5],\n",
       "                    crs: L.CRS.Simple,\n",
       "                    zoom: 10,\n",
       "                    zoomControl: true,\n",
       "                    preferCanvas: false,\n",
       "                }\n",
       "            );\n",
       "            L.control.scale().addTo(map_b39b247aec3d6b4a73e76904f4ffdf5f);\n",
       "\n",
       "            \n",
       "\n",
       "        \n",
       "    \n",
       "            map_b39b247aec3d6b4a73e76904f4ffdf5f.fitBounds(\n",
       "                [[0.0, 0.0], [2.0, 3.0]],\n",
       "                {}\n",
       "            );\n",
       "        \n",
       "    \n",
       "        function geo_json_a10a3d13e99b93f69f5e0e34d601a4bc_styler(feature) {\n",
       "            switch(feature.id) {\n",
       "                case &quot;0&quot;: case &quot;1&quot;: \n",
       "                    return {&quot;color&quot;: &quot;#38588c&quot;, &quot;fillColor&quot;: &quot;#38588c&quot;, &quot;fillOpacity&quot;: 0.5, &quot;weight&quot;: 2};\n",
       "                case &quot;2&quot;: case &quot;4&quot;: \n",
       "                    return {&quot;color&quot;: &quot;#440154&quot;, &quot;fillColor&quot;: &quot;#440154&quot;, &quot;fillOpacity&quot;: 0.5, &quot;weight&quot;: 2};\n",
       "                case &quot;3&quot;: case &quot;5&quot;: \n",
       "                    return {&quot;color&quot;: &quot;#482173&quot;, &quot;fillColor&quot;: &quot;#482173&quot;, &quot;fillOpacity&quot;: 0.5, &quot;weight&quot;: 2};\n",
       "                default:\n",
       "                    return {&quot;color&quot;: &quot;#fde725&quot;, &quot;fillColor&quot;: &quot;#fde725&quot;, &quot;fillOpacity&quot;: 0.5, &quot;weight&quot;: 2};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_a10a3d13e99b93f69f5e0e34d601a4bc_highlighter(feature) {\n",
       "            switch(feature.id) {\n",
       "                default:\n",
       "                    return {&quot;fillOpacity&quot;: 0.75};\n",
       "            }\n",
       "        }\n",
       "        function geo_json_a10a3d13e99b93f69f5e0e34d601a4bc_pointToLayer(feature, latlng) {\n",
       "            var opts = {&quot;bubblingMouseEvents&quot;: true, &quot;color&quot;: &quot;#3388ff&quot;, &quot;dashArray&quot;: null, &quot;dashOffset&quot;: null, &quot;fill&quot;: true, &quot;fillColor&quot;: &quot;#3388ff&quot;, &quot;fillOpacity&quot;: 0.2, &quot;fillRule&quot;: &quot;evenodd&quot;, &quot;lineCap&quot;: &quot;round&quot;, &quot;lineJoin&quot;: &quot;round&quot;, &quot;opacity&quot;: 1.0, &quot;radius&quot;: 2, &quot;stroke&quot;: true, &quot;weight&quot;: 3};\n",
       "            \n",
       "            let style = geo_json_a10a3d13e99b93f69f5e0e34d601a4bc_styler(feature)\n",
       "            Object.assign(opts, style)\n",
       "            \n",
       "            return new L.CircleMarker(latlng, opts)\n",
       "        }\n",
       "\n",
       "        function geo_json_a10a3d13e99b93f69f5e0e34d601a4bc_onEachFeature(feature, layer) {\n",
       "            layer.on({\n",
       "                mouseout: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                            geo_json_a10a3d13e99b93f69f5e0e34d601a4bc.resetStyle(e.target);\n",
       "                    }\n",
       "                },\n",
       "                mouseover: function(e) {\n",
       "                    if(typeof e.target.setStyle === &quot;function&quot;){\n",
       "                        const highlightStyle = geo_json_a10a3d13e99b93f69f5e0e34d601a4bc_highlighter(e.target.feature)\n",
       "                        e.target.setStyle(highlightStyle);\n",
       "                    }\n",
       "                },\n",
       "            });\n",
       "        };\n",
       "        var geo_json_a10a3d13e99b93f69f5e0e34d601a4bc = L.geoJson(null, {\n",
       "                onEachFeature: geo_json_a10a3d13e99b93f69f5e0e34d601a4bc_onEachFeature,\n",
       "            \n",
       "                style: geo_json_a10a3d13e99b93f69f5e0e34d601a4bc_styler,\n",
       "                pointToLayer: geo_json_a10a3d13e99b93f69f5e0e34d601a4bc_pointToLayer,\n",
       "        });\n",
       "\n",
       "        function geo_json_a10a3d13e99b93f69f5e0e34d601a4bc_add (data) {\n",
       "            geo_json_a10a3d13e99b93f69f5e0e34d601a4bc\n",
       "                .addData(data);\n",
       "        }\n",
       "            geo_json_a10a3d13e99b93f69f5e0e34d601a4bc_add({&quot;bbox&quot;: [0.0, 0.0, 3.0, 2.0], &quot;features&quot;: [{&quot;bbox&quot;: [0.0, 0.2, 1.0, 1.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[0.0, 1.0], [1.0, 1.0], [1.0, 0.2], [0.0, 0.6], [0.0, 1.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;0&quot;, &quot;properties&quot;: {&quot;__folium_color&quot;: &quot;#38588c&quot;, &quot;area_overlay&quot;: 0.6000000000000001, &quot;area_source&quot;: 1.0, &quot;source_index&quot;: 0, &quot;target_index&quot;: 0, &quot;weights&quot;: 0.6000000000000001}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [0.0, 1.0, 1.0, 1.8], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[0.0, 1.4], [1.0, 1.8], [1.0, 1.0], [0.0, 1.0], [0.0, 1.4]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;1&quot;, &quot;properties&quot;: {&quot;__folium_color&quot;: &quot;#38588c&quot;, &quot;area_overlay&quot;: 0.6, &quot;area_source&quot;: 1.0, &quot;source_index&quot;: 1, &quot;target_index&quot;: 0, &quot;weights&quot;: 0.6}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [1.0, 0.0, 1.5, 1.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[1.0, 1.0], [1.5, 1.0], [1.5, 0.0], [1.0, 0.2], [1.0, 1.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;2&quot;, &quot;properties&quot;: {&quot;__folium_color&quot;: &quot;#440154&quot;, &quot;area_overlay&quot;: 0.45, &quot;area_source&quot;: 1.0, &quot;source_index&quot;: 2, &quot;target_index&quot;: 0, &quot;weights&quot;: 0.45}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [1.5, 0.0, 2.0, 1.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[2.0, 1.0], [2.0, 0.0], [1.5, 0.0], [1.5, 1.0], [2.0, 1.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;3&quot;, &quot;properties&quot;: {&quot;__folium_color&quot;: &quot;#482173&quot;, &quot;area_overlay&quot;: 0.5, &quot;area_source&quot;: 1.0, &quot;source_index&quot;: 2, &quot;target_index&quot;: 1, &quot;weights&quot;: 0.5}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [1.0, 1.0, 1.5, 2.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[1.0, 1.8], [1.5, 2.0], [1.5, 1.0], [1.0, 1.0], [1.0, 1.8]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;4&quot;, &quot;properties&quot;: {&quot;__folium_color&quot;: &quot;#440154&quot;, &quot;area_overlay&quot;: 0.45, &quot;area_source&quot;: 1.0, &quot;source_index&quot;: 3, &quot;target_index&quot;: 0, &quot;weights&quot;: 0.45}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [1.5, 1.0, 2.0, 2.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[2.0, 2.0], [2.0, 1.0], [1.5, 1.0], [1.5, 2.0], [2.0, 2.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;5&quot;, &quot;properties&quot;: {&quot;__folium_color&quot;: &quot;#482173&quot;, &quot;area_overlay&quot;: 0.5, &quot;area_source&quot;: 1.0, &quot;source_index&quot;: 3, &quot;target_index&quot;: 1, &quot;weights&quot;: 0.5}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [2.0, 0.0, 3.0, 1.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[2.0, 1.0], [3.0, 1.0], [3.0, 0.0], [2.0, 0.0], [2.0, 1.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;6&quot;, &quot;properties&quot;: {&quot;__folium_color&quot;: &quot;#fde725&quot;, &quot;area_overlay&quot;: 1.0, &quot;area_source&quot;: 1.0, &quot;source_index&quot;: 4, &quot;target_index&quot;: 1, &quot;weights&quot;: 1.0}, &quot;type&quot;: &quot;Feature&quot;}, {&quot;bbox&quot;: [2.0, 1.0, 3.0, 2.0], &quot;geometry&quot;: {&quot;coordinates&quot;: [[[2.0, 2.0], [3.0, 2.0], [3.0, 1.0], [2.0, 1.0], [2.0, 2.0]]], &quot;type&quot;: &quot;Polygon&quot;}, &quot;id&quot;: &quot;7&quot;, &quot;properties&quot;: {&quot;__folium_color&quot;: &quot;#fde725&quot;, &quot;area_overlay&quot;: 1.0, &quot;area_source&quot;: 1.0, &quot;source_index&quot;: 5, &quot;target_index&quot;: 1, &quot;weights&quot;: 1.0}, &quot;type&quot;: &quot;Feature&quot;}], &quot;type&quot;: &quot;FeatureCollection&quot;});\n",
       "\n",
       "        \n",
       "    \n",
       "    geo_json_a10a3d13e99b93f69f5e0e34d601a4bc.bindTooltip(\n",
       "    function(layer){\n",
       "    let div = L.DomUtil.create(&#x27;div&#x27;);\n",
       "    \n",
       "    let handleObject = feature=&gt;typeof(feature)==&#x27;object&#x27; ? JSON.stringify(feature) : feature;\n",
       "    let fields = [&quot;source_index&quot;, &quot;target_index&quot;, &quot;area_overlay&quot;, &quot;area_source&quot;, &quot;weights&quot;];\n",
       "    let aliases = [&quot;source_index&quot;, &quot;target_index&quot;, &quot;area_overlay&quot;, &quot;area_source&quot;, &quot;weights&quot;];\n",
       "    let table = &#x27;&lt;table&gt;&#x27; +\n",
       "        String(\n",
       "        fields.map(\n",
       "        (v,i)=&gt;\n",
       "        `&lt;tr&gt;\n",
       "            &lt;th&gt;${aliases[i]}&lt;/th&gt;\n",
       "            \n",
       "            &lt;td&gt;${handleObject(layer.feature.properties[v])}&lt;/td&gt;\n",
       "        &lt;/tr&gt;`).join(&#x27;&#x27;))\n",
       "    +&#x27;&lt;/table&gt;&#x27;;\n",
       "    div.innerHTML=table;\n",
       "    \n",
       "    return div\n",
       "    }\n",
       "    ,{&quot;className&quot;: &quot;foliumtooltip&quot;, &quot;sticky&quot;: true});\n",
       "                     \n",
       "    \n",
       "            geo_json_a10a3d13e99b93f69f5e0e34d601a4bc.addTo(map_b39b247aec3d6b4a73e76904f4ffdf5f);\n",
       "        \n",
       "    \n",
       "    var color_map_d2f09e21a56bd360e2497366a340a351 = {};\n",
       "\n",
       "    \n",
       "    color_map_d2f09e21a56bd360e2497366a340a351.color = d3.scale.threshold()\n",
       "              .domain([0.45, 0.45110220440881765, 0.4522044088176353, 0.45330661322645294, 0.4544088176352706, 0.45551102204408817, 0.4566132264529058, 0.45771543086172345, 0.4588176352705411, 0.45991983967935873, 0.4610220440881764, 0.462124248496994, 0.46322645290581166, 0.46432865731462925, 0.4654308617234469, 0.46653306613226453, 0.46763527054108217, 0.4687374749498998, 0.46983967935871745, 0.4709418837675351, 0.47204408817635274, 0.4731462925851704, 0.47424849699398797, 0.4753507014028056, 0.47645290581162325, 0.4775551102204409, 0.47865731462925853, 0.4797595190380762, 0.4808617234468938, 0.48196392785571146, 0.48306613226452905, 0.4841683366733467, 0.48527054108216433, 0.486372745490982, 0.4874749498997996, 0.48857715430861726, 0.4896793587174349, 0.49078156312625254, 0.4918837675350701, 0.49298597194388777, 0.4940881763527054, 0.49519038076152305, 0.4962925851703407, 0.49739478957915834, 0.498496993987976, 0.4995991983967936, 0.5007014028056113, 0.5018036072144288, 0.5029058116232465, 0.5040080160320641, 0.5051102204408818, 0.5062124248496994, 0.507314629258517, 0.5084168336673347, 0.5095190380761523, 0.51062124248497, 0.5117234468937876, 0.5128256513026053, 0.5139278557114229, 0.5150300601202404, 0.5161322645290581, 0.5172344689378757, 0.5183366733466934, 0.519438877755511, 0.5205410821643287, 0.5216432865731463, 0.522745490981964, 0.5238476953907816, 0.5249498997995992, 0.5260521042084169, 0.5271543086172344, 0.5282565130260521, 0.5293587174348697, 0.5304609218436874, 0.531563126252505, 0.5326653306613227, 0.5337675350701403, 0.5348697394789579, 0.5359719438877756, 0.5370741482965932, 0.5381763527054109, 0.5392785571142285, 0.5403807615230461, 0.5414829659318637, 0.5425851703406814, 0.543687374749499, 0.5447895791583166, 0.5458917835671343, 0.5469939879759519, 0.5480961923847696, 0.5491983967935872, 0.5503006012024049, 0.5514028056112225, 0.5525050100200402, 0.5536072144288577, 0.5547094188376753, 0.555811623246493, 0.5569138276553106, 0.5580160320641283, 0.5591182364729459, 0.5602204408817636, 0.5613226452905812, 0.5624248496993988, 0.5635270541082165, 0.564629258517034, 0.5657314629258517, 0.5668336673346693, 0.567935871743487, 0.5690380761523046, 0.5701402805611222, 0.5712424849699399, 0.5723446893787575, 0.5734468937875752, 0.5745490981963928, 0.5756513026052105, 0.576753507014028, 0.5778557114228458, 0.5789579158316633, 0.5800601202404809, 0.5811623246492986, 0.5822645290581162, 0.5833667334669339, 0.5844689378757515, 0.5855711422845692, 0.5866733466933868, 0.5877755511022045, 0.5888777555110221, 0.5899799599198396, 0.5910821643286573, 0.5921843687374749, 0.5932865731462926, 0.5943887775551102, 0.5954909819639279, 0.5965931863727455, 0.5976953907815632, 0.5987975951903808, 0.5998997995991984, 0.6010020040080161, 0.6021042084168337, 0.6032064128256514, 0.6043086172344689, 0.6054108216432866, 0.6065130260521042, 0.6076152304609219, 0.6087174348697395, 0.6098196392785571, 0.6109218436873748, 0.6120240480961924, 0.6131262525050101, 0.6142284569138277, 0.6153306613226452, 0.6164328657314629, 0.6175350701402806, 0.6186372745490982, 0.6197394789579158, 0.6208416833667335, 0.6219438877755511, 0.6230460921843688, 0.6241482965931864, 0.625250501002004, 0.6263527054108217, 0.6274549098196394, 0.628557114228457, 0.6296593186372745, 0.6307615230460922, 0.6318637274549098, 0.6329659318637275, 0.6340681362725451, 0.6351703406813627, 0.6362725450901804, 0.6373747494989981, 0.6384769539078157, 0.6395791583166333, 0.640681362725451, 0.6417835671342685, 0.6428857715430862, 0.6439879759519038, 0.6450901803607214, 0.6461923847695391, 0.6472945891783567, 0.6483967935871744, 0.649498997995992, 0.6506012024048097, 0.6517034068136273, 0.652805611222445, 0.6539078156312625, 0.6550100200400801, 0.6561122244488978, 0.6572144288577154, 0.6583166332665331, 0.6594188376753507, 0.6605210420841684, 0.661623246492986, 0.6627254509018037, 0.6638276553106213, 0.6649298597194389, 0.6660320641282566, 0.6671342685370741, 0.6682364729458918, 0.6693386773547094, 0.6704408817635271, 0.6715430861723447, 0.6726452905811624, 0.67374749498998, 0.6748496993987976, 0.6759519038076153, 0.6770541082164329, 0.6781563126252506, 0.6792585170340681, 0.6803607214428857, 0.6814629258517034, 0.6825651302605211, 0.6836673346693387, 0.6847695390781563, 0.685871743486974, 0.6869739478957916, 0.6880761523046093, 0.6891783567134269, 0.6902805611222445, 0.6913827655310621, 0.6924849699398798, 0.6935871743486974, 0.694689378757515, 0.6957915831663327, 0.6968937875751504, 0.697995991983968, 0.6990981963927856, 0.7002004008016032, 0.7013026052104209, 0.7024048096192386, 0.7035070140280562, 0.7046092184368737, 0.7057114228456914, 0.706813627254509, 0.7079158316633267, 0.7090180360721443, 0.7101202404809619, 0.7112224448897796, 0.7123246492985973, 0.7134268537074149, 0.7145290581162325, 0.7156312625250502, 0.7167334669338679, 0.7178356713426854, 0.718937875751503, 0.7200400801603206, 0.7211422845691383, 0.722244488977956, 0.7233466933867736, 0.7244488977955912, 0.7255511022044088, 0.7266533066132265, 0.7277555110220442, 0.7288577154308618, 0.7299599198396793, 0.7310621242484969, 0.7321643286573147, 0.7332665330661323, 0.7343687374749499, 0.7354709418837676, 0.7365731462925852, 0.7376753507014029, 0.7387775551102205, 0.7398797595190381, 0.7409819639278558, 0.7420841683366733, 0.743186372745491, 0.7442885771543086, 0.7453907815631262, 0.7464929859719439, 0.7475951903807616, 0.7486973947895792, 0.7497995991983968, 0.7509018036072144, 0.7520040080160322, 0.7531062124248498, 0.7542084168336673, 0.755310621242485, 0.7564128256513026, 0.7575150300601203, 0.7586172344689379, 0.7597194388777555, 0.7608216432865732, 0.7619238476953908, 0.7630260521042085, 0.7641282565130261, 0.7652304609218437, 0.7663326653306614, 0.767434869739479, 0.7685370741482966, 0.7696392785571142, 0.7707414829659319, 0.7718436873747496, 0.7729458917835672, 0.7740480961923848, 0.7751503006012024, 0.7762525050100201, 0.7773547094188378, 0.7784569138276554, 0.779559118236473, 0.7806613226452905, 0.7817635270541082, 0.7828657314629259, 0.7839679358717435, 0.7850701402805611, 0.7861723446893787, 0.7872745490981965, 0.7883767535070141, 0.7894789579158317, 0.7905811623246494, 0.791683366733467, 0.7927855711422847, 0.7938877755511022, 0.7949899799599198, 0.7960921843687375, 0.7971943887775551, 0.7982965931863728, 0.7993987975951904, 0.800501002004008, 0.8016032064128257, 0.8027054108216434, 0.803807615230461, 0.8049098196392785, 0.8060120240480961, 0.8071142284569139, 0.8082164328657315, 0.8093186372745491, 0.8104208416833668, 0.8115230460921844, 0.8126252505010021, 0.8137274549098197, 0.8148296593186373, 0.815931863727455, 0.8170340681362726, 0.8181362725450902, 0.8192384769539078, 0.8203406813627254, 0.8214428857715431, 0.8225450901803608, 0.8236472945891784, 0.824749498997996, 0.8258517034068137, 0.8269539078156314, 0.828056112224449, 0.8291583166332666, 0.8302605210420843, 0.8313627254509018, 0.8324649298597195, 0.8335671342685371, 0.8346693386773547, 0.8357715430861724, 0.83687374749499, 0.8379759519038077, 0.8390781563126253, 0.8401803607214429, 0.8412825651302606, 0.8423847695390783, 0.8434869739478958, 0.8445891783567134, 0.8456913827655311, 0.8467935871743488, 0.8478957915831664, 0.848997995991984, 0.8501002004008016, 0.8512024048096193, 0.852304609218437, 0.8534068136272546, 0.8545090180360722, 0.8556112224448897, 0.8567134268537074, 0.8578156312625251, 0.8589178356713427, 0.8600200400801603, 0.861122244488978, 0.8622244488977957, 0.8633266533066133, 0.8644288577154309, 0.8655310621242486, 0.8666332665330662, 0.8677354709418839, 0.8688376753507014, 0.869939879759519, 0.8710420841683367, 0.8721442885771543, 0.873246492985972, 0.8743486973947896, 0.8754509018036072, 0.8765531062124249, 0.8776553106212426, 0.8787575150300602, 0.8798597194388778, 0.8809619238476954, 0.8820641282565131, 0.8831663326653307, 0.8842685370741483, 0.885370741482966, 0.8864729458917836, 0.8875751503006013, 0.8886773547094189, 0.8897795591182365, 0.8908817635270542, 0.8919839679358718, 0.8930861723446895, 0.894188376753507, 0.8952905811623246, 0.8963927855711423, 0.89749498997996, 0.8985971943887776, 0.8996993987975952, 0.9008016032064129, 0.9019038076152306, 0.9030060120240482, 0.9041082164328658, 0.9052104208416833, 0.906312625250501, 0.9074148296593187, 0.9085170340681363, 0.9096192384769539, 0.9107214428857715, 0.9118236472945892, 0.9129258517034069, 0.9140280561122245, 0.9151302605210421, 0.9162324649298598, 0.9173346693386775, 0.918436873747495, 0.9195390781563126, 0.9206412825651303, 0.9217434869739479, 0.9228456913827656, 0.9239478957915832, 0.9250501002004008, 0.9261523046092185, 0.9272545090180361, 0.9283567134268538, 0.9294589178356714, 0.930561122244489, 0.9316633266533066, 0.9327655310621243, 0.9338677354709419, 0.9349699398797595, 0.9360721442885772, 0.9371743486973949, 0.9382765531062125, 0.9393787575150301, 0.9404809619238478, 0.9415831663326654, 0.9426853707414831, 0.9437875751503007, 0.9448897795591182, 0.9459919839679359, 0.9470941883767535, 0.9481963927855712, 0.9492985971943888, 0.9504008016032064, 0.9515030060120242, 0.9526052104208418, 0.9537074148296594, 0.954809619238477, 0.9559118236472945, 0.9570140280561124, 0.9581162324649299, 0.9592184368737475, 0.9603206412825651, 0.9614228456913829, 0.9625250501002005, 0.9636272545090181, 0.9647294589178357, 0.9658316633266535, 0.9669338677354711, 0.9680360721442887, 0.9691382765531062, 0.9702404809619238, 0.9713426853707416, 0.9724448897795592, 0.9735470941883768, 0.9746492985971944, 0.975751503006012, 0.9768537074148298, 0.9779559118236474, 0.979058116232465, 0.9801603206412826, 0.9812625250501004, 0.982364729458918, 0.9834669338677355, 0.9845691382765531, 0.9856713426853707, 0.9867735470941885, 0.9878757515030061, 0.9889779559118237, 0.9900801603206413, 0.9911823647294589, 0.9922845691382767, 0.9933867735470943, 0.9944889779559118, 0.9955911823647294, 0.996693386773547, 0.9977955911823648, 0.9988977955911824, 1.0])\n",
       "              .range([&#x27;#440154ff&#x27;, &#x27;#440155ff&#x27;, &#x27;#440256ff&#x27;, &#x27;#440356ff&#x27;, &#x27;#450457ff&#x27;, &#x27;#450458ff&#x27;, &#x27;#450559ff&#x27;, &#x27;#450659ff&#x27;, &#x27;#46075aff&#x27;, &#x27;#46075bff&#x27;, &#x27;#46085cff&#x27;, &#x27;#46095cff&#x27;, &#x27;#460a5dff&#x27;, &#x27;#460a5eff&#x27;, &#x27;#460b5eff&#x27;, &#x27;#460c5fff&#x27;, &#x27;#470d60ff&#x27;, &#x27;#470d61ff&#x27;, &#x27;#470e61ff&#x27;, &#x27;#470f62ff&#x27;, &#x27;#471063ff&#x27;, &#x27;#471064ff&#x27;, &#x27;#471164ff&#x27;, &#x27;#471265ff&#x27;, &#x27;#471365ff&#x27;, &#x27;#481366ff&#x27;, &#x27;#481467ff&#x27;, &#x27;#481568ff&#x27;, &#x27;#481668ff&#x27;, &#x27;#481669ff&#x27;, &#x27;#481769ff&#x27;, &#x27;#48176aff&#x27;, &#x27;#48186bff&#x27;, &#x27;#48196cff&#x27;, &#x27;#481a6cff&#x27;, &#x27;#481a6dff&#x27;, &#x27;#481b6dff&#x27;, &#x27;#481c6eff&#x27;, &#x27;#481c6eff&#x27;, &#x27;#481d6fff&#x27;, &#x27;#481d6fff&#x27;, &#x27;#481f70ff&#x27;, &#x27;#481f70ff&#x27;, &#x27;#482071ff&#x27;, &#x27;#482072ff&#x27;, &#x27;#482173ff&#x27;, &#x27;#482273ff&#x27;, &#x27;#482374ff&#x27;, &#x27;#482374ff&#x27;, &#x27;#482475ff&#x27;, &#x27;#482476ff&#x27;, &#x27;#482576ff&#x27;, &#x27;#482577ff&#x27;, &#x27;#482677ff&#x27;, &#x27;#482778ff&#x27;, &#x27;#482878ff&#x27;, &#x27;#482879ff&#x27;, &#x27;#482979ff&#x27;, &#x27;#47297aff&#x27;, &#x27;#472a7aff&#x27;, &#x27;#472b7aff&#x27;, &#x27;#472c7aff&#x27;, &#x27;#472c7bff&#x27;, &#x27;#472d7bff&#x27;, &#x27;#472d7cff&#x27;, &#x27;#472e7cff&#x27;, &#x27;#472e7dff&#x27;, &#x27;#472f7dff&#x27;, &#x27;#462f7eff&#x27;, &#x27;#46307eff&#x27;, &#x27;#46317eff&#x27;, &#x27;#46327eff&#x27;, &#x27;#46327fff&#x27;, &#x27;#46337fff&#x27;, &#x27;#463480ff&#x27;, &#x27;#453480ff&#x27;, &#x27;#453581ff&#x27;, &#x27;#453581ff&#x27;, &#x27;#453681ff&#x27;, &#x27;#453781ff&#x27;, &#x27;#453882ff&#x27;, &#x27;#443882ff&#x27;, &#x27;#443983ff&#x27;, &#x27;#443983ff&#x27;, &#x27;#443a83ff&#x27;, &#x27;#443a83ff&#x27;, &#x27;#443b84ff&#x27;, &#x27;#433c84ff&#x27;, &#x27;#433d84ff&#x27;, &#x27;#433d85ff&#x27;, &#x27;#433e85ff&#x27;, &#x27;#423e85ff&#x27;, &#x27;#423f85ff&#x27;, &#x27;#423f86ff&#x27;, &#x27;#424086ff&#x27;, &#x27;#424086ff&#x27;, &#x27;#424186ff&#x27;, &#x27;#414187ff&#x27;, &#x27;#414287ff&#x27;, &#x27;#414387ff&#x27;, &#x27;#414487ff&#x27;, &#x27;#404488ff&#x27;, &#x27;#404588ff&#x27;, &#x27;#404588ff&#x27;, &#x27;#404688ff&#x27;, &#x27;#3f4688ff&#x27;, &#x27;#3f4788ff&#x27;, &#x27;#3f4789ff&#x27;, &#x27;#3f4889ff&#x27;, &#x27;#3e4889ff&#x27;, &#x27;#3e4989ff&#x27;, &#x27;#3e4a89ff&#x27;, &#x27;#3e4a89ff&#x27;, &#x27;#3e4b8aff&#x27;, &#x27;#3d4c8aff&#x27;, &#x27;#3d4d8aff&#x27;, &#x27;#3d4d8aff&#x27;, &#x27;#3d4e8aff&#x27;, &#x27;#3c4e8aff&#x27;, &#x27;#3c4f8aff&#x27;, &#x27;#3c4f8aff&#x27;, &#x27;#3c508bff&#x27;, &#x27;#3b508bff&#x27;, &#x27;#3b518bff&#x27;, &#x27;#3b518bff&#x27;, &#x27;#3b528bff&#x27;, &#x27;#3a528bff&#x27;, &#x27;#3a538bff&#x27;, &#x27;#3a538bff&#x27;, &#x27;#3a548cff&#x27;, &#x27;#39548cff&#x27;, &#x27;#39558cff&#x27;, &#x27;#39558cff&#x27;, &#x27;#39568cff&#x27;, &#x27;#38578cff&#x27;, &#x27;#38588cff&#x27;, &#x27;#38588cff&#x27;, &#x27;#38598cff&#x27;, &#x27;#37598cff&#x27;, &#x27;#375a8cff&#x27;, &#x27;#375a8dff&#x27;, &#x27;#375b8dff&#x27;, &#x27;#365b8dff&#x27;, &#x27;#365c8dff&#x27;, &#x27;#365c8dff&#x27;, &#x27;#365d8dff&#x27;, &#x27;#355d8dff&#x27;, &#x27;#355e8dff&#x27;, &#x27;#355f8dff&#x27;, &#x27;#355f8dff&#x27;, &#x27;#34608dff&#x27;, &#x27;#34608dff&#x27;, &#x27;#34618dff&#x27;, &#x27;#34618dff&#x27;, &#x27;#33628dff&#x27;, &#x27;#33628dff&#x27;, &#x27;#33638dff&#x27;, &#x27;#32638dff&#x27;, &#x27;#32648eff&#x27;, &#x27;#32648eff&#x27;, &#x27;#32658eff&#x27;, &#x27;#31658eff&#x27;, &#x27;#31668eff&#x27;, &#x27;#31668eff&#x27;, &#x27;#31678eff&#x27;, &#x27;#31678eff&#x27;, &#x27;#31688eff&#x27;, &#x27;#30688eff&#x27;, &#x27;#30698eff&#x27;, &#x27;#30698eff&#x27;, &#x27;#306a8eff&#x27;, &#x27;#2f6a8eff&#x27;, &#x27;#2f6b8eff&#x27;, &#x27;#2f6b8eff&#x27;, &#x27;#2f6c8eff&#x27;, &#x27;#2e6c8eff&#x27;, &#x27;#2e6d8eff&#x27;, &#x27;#2e6d8eff&#x27;, &#x27;#2e6e8eff&#x27;, &#x27;#2e6e8eff&#x27;, &#x27;#2e6f8eff&#x27;, &#x27;#2d6f8eff&#x27;, &#x27;#2d708eff&#x27;, &#x27;#2d708eff&#x27;, &#x27;#2d718eff&#x27;, &#x27;#2c718eff&#x27;, &#x27;#2c718eff&#x27;, &#x27;#2c728eff&#x27;, &#x27;#2c728eff&#x27;, &#x27;#2c738eff&#x27;, &#x27;#2c738eff&#x27;, &#x27;#2b748eff&#x27;, &#x27;#2b748eff&#x27;, &#x27;#2b758eff&#x27;, &#x27;#2b758eff&#x27;, &#x27;#2a768eff&#x27;, &#x27;#2a768eff&#x27;, &#x27;#2a778eff&#x27;, &#x27;#2a778eff&#x27;, &#x27;#2a788eff&#x27;, &#x27;#29788eff&#x27;, &#x27;#29798eff&#x27;, &#x27;#29798eff&#x27;, &#x27;#297a8eff&#x27;, &#x27;#297a8eff&#x27;, &#x27;#297b8eff&#x27;, &#x27;#287b8eff&#x27;, &#x27;#287c8eff&#x27;, &#x27;#287c8eff&#x27;, &#x27;#287d8eff&#x27;, &#x27;#277d8eff&#x27;, &#x27;#277e8eff&#x27;, &#x27;#277e8eff&#x27;, &#x27;#277f8eff&#x27;, &#x27;#277f8eff&#x27;, &#x27;#27808eff&#x27;, &#x27;#26808eff&#x27;, &#x27;#26818eff&#x27;, &#x27;#26818eff&#x27;, &#x27;#26828eff&#x27;, &#x27;#26828eff&#x27;, &#x27;#26828eff&#x27;, &#x27;#25828eff&#x27;, &#x27;#25838eff&#x27;, &#x27;#25838eff&#x27;, &#x27;#25848eff&#x27;, &#x27;#25858eff&#x27;, &#x27;#25858eff&#x27;, &#x27;#24868eff&#x27;, &#x27;#24868eff&#x27;, &#x27;#24878eff&#x27;, &#x27;#24878eff&#x27;, &#x27;#23888eff&#x27;, &#x27;#23888eff&#x27;, &#x27;#23898eff&#x27;, &#x27;#23898eff&#x27;, &#x27;#238a8dff&#x27;, &#x27;#238a8dff&#x27;, &#x27;#228b8dff&#x27;, &#x27;#228b8dff&#x27;, &#x27;#228c8dff&#x27;, &#x27;#228c8dff&#x27;, &#x27;#228d8dff&#x27;, &#x27;#218d8dff&#x27;, &#x27;#218e8dff&#x27;, &#x27;#218e8dff&#x27;, &#x27;#218f8dff&#x27;, &#x27;#218f8dff&#x27;, &#x27;#21908dff&#x27;, &#x27;#21908dff&#x27;, &#x27;#21918cff&#x27;, &#x27;#20918cff&#x27;, &#x27;#20928cff&#x27;, &#x27;#20928cff&#x27;, &#x27;#20928cff&#x27;, &#x27;#20928cff&#x27;, &#x27;#20938cff&#x27;, &#x27;#1f938cff&#x27;, &#x27;#1f948cff&#x27;, &#x27;#1f948cff&#x27;, &#x27;#1f958bff&#x27;, &#x27;#1f958bff&#x27;, &#x27;#1f968bff&#x27;, &#x27;#1f968bff&#x27;, &#x27;#1f978bff&#x27;, &#x27;#1f988bff&#x27;, &#x27;#1f988bff&#x27;, &#x27;#1f998bff&#x27;, &#x27;#1f998aff&#x27;, &#x27;#1f9a8aff&#x27;, &#x27;#1f9a8aff&#x27;, &#x27;#1e9b8aff&#x27;, &#x27;#1e9b8aff&#x27;, &#x27;#1e9c8aff&#x27;, &#x27;#1e9c89ff&#x27;, &#x27;#1e9d89ff&#x27;, &#x27;#1e9d89ff&#x27;, &#x27;#1e9e89ff&#x27;, &#x27;#1f9e89ff&#x27;, &#x27;#1f9f88ff&#x27;, &#x27;#1f9f88ff&#x27;, &#x27;#1fa088ff&#x27;, &#x27;#1fa088ff&#x27;, &#x27;#1fa188ff&#x27;, &#x27;#1fa188ff&#x27;, &#x27;#1fa187ff&#x27;, &#x27;#1fa187ff&#x27;, &#x27;#1fa287ff&#x27;, &#x27;#1fa287ff&#x27;, &#x27;#1fa386ff&#x27;, &#x27;#20a386ff&#x27;, &#x27;#20a486ff&#x27;, &#x27;#20a486ff&#x27;, &#x27;#20a585ff&#x27;, &#x27;#21a585ff&#x27;, &#x27;#21a685ff&#x27;, &#x27;#21a685ff&#x27;, &#x27;#21a785ff&#x27;, &#x27;#22a785ff&#x27;, &#x27;#22a884ff&#x27;, &#x27;#22a884ff&#x27;, &#x27;#22a983ff&#x27;, &#x27;#23a983ff&#x27;, &#x27;#23aa83ff&#x27;, &#x27;#24ab83ff&#x27;, &#x27;#25ab82ff&#x27;, &#x27;#25ac82ff&#x27;, &#x27;#25ac82ff&#x27;, &#x27;#25ad82ff&#x27;, &#x27;#26ad81ff&#x27;, &#x27;#26ad81ff&#x27;, &#x27;#27ad81ff&#x27;, &#x27;#27ae81ff&#x27;, &#x27;#28ae80ff&#x27;, &#x27;#28af80ff&#x27;, &#x27;#29af7fff&#x27;, &#x27;#29b07fff&#x27;, &#x27;#2ab07fff&#x27;, &#x27;#2bb17eff&#x27;, &#x27;#2cb17eff&#x27;, &#x27;#2cb27dff&#x27;, &#x27;#2db27dff&#x27;, &#x27;#2db37cff&#x27;, &#x27;#2eb37cff&#x27;, &#x27;#2eb47cff&#x27;, &#x27;#2fb47cff&#x27;, &#x27;#30b57bff&#x27;, &#x27;#31b57bff&#x27;, &#x27;#31b67aff&#x27;, &#x27;#32b67aff&#x27;, &#x27;#33b679ff&#x27;, &#x27;#34b679ff&#x27;, &#x27;#34b779ff&#x27;, &#x27;#35b779ff&#x27;, &#x27;#36b878ff&#x27;, &#x27;#37b878ff&#x27;, &#x27;#37b977ff&#x27;, &#x27;#38b977ff&#x27;, &#x27;#39ba76ff&#x27;, &#x27;#3aba76ff&#x27;, &#x27;#3abb75ff&#x27;, &#x27;#3bbb75ff&#x27;, &#x27;#3cbc74ff&#x27;, &#x27;#3dbc74ff&#x27;, &#x27;#3ebc73ff&#x27;, &#x27;#3fbd73ff&#x27;, &#x27;#40bd72ff&#x27;, &#x27;#40be72ff&#x27;, &#x27;#41be71ff&#x27;, &#x27;#42bf71ff&#x27;, &#x27;#43bf70ff&#x27;, &#x27;#45c070ff&#x27;, &#x27;#46c06fff&#x27;, &#x27;#47c16fff&#x27;, &#x27;#48c16eff&#x27;, &#x27;#49c16eff&#x27;, &#x27;#4ac16dff&#x27;, &#x27;#4bc26cff&#x27;, &#x27;#4cc26cff&#x27;, &#x27;#4dc36bff&#x27;, &#x27;#4ec36bff&#x27;, &#x27;#4fc46aff&#x27;, &#x27;#50c46aff&#x27;, &#x27;#51c569ff&#x27;, &#x27;#52c569ff&#x27;, &#x27;#53c568ff&#x27;, &#x27;#54c568ff&#x27;, &#x27;#55c667ff&#x27;, &#x27;#56c667ff&#x27;, &#x27;#57c766ff&#x27;, &#x27;#58c765ff&#x27;, &#x27;#59c864ff&#x27;, &#x27;#5ac864ff&#x27;, &#x27;#5bc863ff&#x27;, &#x27;#5cc863ff&#x27;, &#x27;#5dc962ff&#x27;, &#x27;#5ec962ff&#x27;, &#x27;#5fca61ff&#x27;, &#x27;#60ca60ff&#x27;, &#x27;#62cb5fff&#x27;, &#x27;#63cb5fff&#x27;, &#x27;#64cb5eff&#x27;, &#x27;#65cc5dff&#x27;, &#x27;#66cc5cff&#x27;, &#x27;#67cd5cff&#x27;, &#x27;#68cd5bff&#x27;, &#x27;#6acd5bff&#x27;, &#x27;#6bcd5aff&#x27;, &#x27;#6cce59ff&#x27;, &#x27;#6ece58ff&#x27;, &#x27;#6fcf58ff&#x27;, &#x27;#70cf57ff&#x27;, &#x27;#71d057ff&#x27;, &#x27;#72d056ff&#x27;, &#x27;#74d055ff&#x27;, &#x27;#75d054ff&#x27;, &#x27;#76d153ff&#x27;, &#x27;#77d153ff&#x27;, &#x27;#78d152ff&#x27;, &#x27;#7ad151ff&#x27;, &#x27;#7bd250ff&#x27;, &#x27;#7cd250ff&#x27;, &#x27;#7dd34fff&#x27;, &#x27;#7fd34eff&#x27;, &#x27;#80d34dff&#x27;, &#x27;#81d34dff&#x27;, &#x27;#82d44cff&#x27;, &#x27;#84d44bff&#x27;, &#x27;#85d54aff&#x27;, &#x27;#86d549ff&#x27;, &#x27;#88d548ff&#x27;, &#x27;#89d548ff&#x27;, &#x27;#8ad647ff&#x27;, &#x27;#8bd646ff&#x27;, &#x27;#8dd645ff&#x27;, &#x27;#8ed645ff&#x27;, &#x27;#8fd744ff&#x27;, &#x27;#90d743ff&#x27;, &#x27;#92d742ff&#x27;, &#x27;#93d741ff&#x27;, &#x27;#94d840ff&#x27;, &#x27;#96d83fff&#x27;, &#x27;#97d83eff&#x27;, &#x27;#99d93dff&#x27;, &#x27;#9ad93cff&#x27;, &#x27;#9bd93cff&#x27;, &#x27;#9dd93bff&#x27;, &#x27;#9eda3aff&#x27;, &#x27;#9fda39ff&#x27;, &#x27;#a1da38ff&#x27;, &#x27;#a2da37ff&#x27;, &#x27;#a3db36ff&#x27;, &#x27;#a4db36ff&#x27;, &#x27;#a6db35ff&#x27;, &#x27;#a8db34ff&#x27;, &#x27;#a9dc33ff&#x27;, &#x27;#aadc32ff&#x27;, &#x27;#abdc31ff&#x27;, &#x27;#addc30ff&#x27;, &#x27;#aedd2fff&#x27;, &#x27;#b0dd2fff&#x27;, &#x27;#b1dd2eff&#x27;, &#x27;#b2dd2dff&#x27;, &#x27;#b3de2cff&#x27;, &#x27;#b5de2bff&#x27;, &#x27;#b6de2aff&#x27;, &#x27;#b8de29ff&#x27;, &#x27;#b9de28ff&#x27;, &#x27;#bade28ff&#x27;, &#x27;#bcdf27ff&#x27;, &#x27;#bddf26ff&#x27;, &#x27;#bfdf25ff&#x27;, &#x27;#c0df25ff&#x27;, &#x27;#c1df24ff&#x27;, &#x27;#c2df23ff&#x27;, &#x27;#c4e022ff&#x27;, &#x27;#c5e021ff&#x27;, &#x27;#c7e020ff&#x27;, &#x27;#c8e020ff&#x27;, &#x27;#c9e11fff&#x27;, &#x27;#cbe11eff&#x27;, &#x27;#cce11dff&#x27;, &#x27;#cee11dff&#x27;, &#x27;#cfe11cff&#x27;, &#x27;#d1e11bff&#x27;, &#x27;#d2e21bff&#x27;, &#x27;#d3e21aff&#x27;, &#x27;#d4e21aff&#x27;, &#x27;#d6e219ff&#x27;, &#x27;#d7e219ff&#x27;, &#x27;#d9e319ff&#x27;, &#x27;#dae319ff&#x27;, &#x27;#dbe318ff&#x27;, &#x27;#dde318ff&#x27;, &#x27;#dee318ff&#x27;, &#x27;#dfe318ff&#x27;, &#x27;#e0e418ff&#x27;, &#x27;#e2e418ff&#x27;, &#x27;#e3e418ff&#x27;, &#x27;#e5e418ff&#x27;, &#x27;#e6e419ff&#x27;, &#x27;#e7e419ff&#x27;, &#x27;#e8e519ff&#x27;, &#x27;#eae519ff&#x27;, &#x27;#ebe51aff&#x27;, &#x27;#ece51aff&#x27;, &#x27;#ede51bff&#x27;, &#x27;#efe51bff&#x27;, &#x27;#f0e51cff&#x27;, &#x27;#f1e51dff&#x27;, &#x27;#f3e61dff&#x27;, &#x27;#f4e61eff&#x27;, &#x27;#f5e61eff&#x27;, &#x27;#f6e61fff&#x27;, &#x27;#f7e620ff&#x27;, &#x27;#f8e621ff&#x27;, &#x27;#fae722ff&#x27;, &#x27;#fbe723ff&#x27;, &#x27;#fce724ff&#x27;, &#x27;#fde725ff&#x27;]);\n",
       "    \n",
       "\n",
       "    color_map_d2f09e21a56bd360e2497366a340a351.x = d3.scale.linear()\n",
       "              .domain([0.45, 1.0])\n",
       "              .range([0, 450 - 50]);\n",
       "\n",
       "    color_map_d2f09e21a56bd360e2497366a340a351.legend = L.control({position: &#x27;topright&#x27;});\n",
       "    color_map_d2f09e21a56bd360e2497366a340a351.legend.onAdd = function (map) {var div = L.DomUtil.create(&#x27;div&#x27;, &#x27;legend&#x27;); return div};\n",
       "    color_map_d2f09e21a56bd360e2497366a340a351.legend.addTo(map_b39b247aec3d6b4a73e76904f4ffdf5f);\n",
       "\n",
       "    color_map_d2f09e21a56bd360e2497366a340a351.xAxis = d3.svg.axis()\n",
       "        .scale(color_map_d2f09e21a56bd360e2497366a340a351.x)\n",
       "        .orient(&quot;top&quot;)\n",
       "        .tickSize(1)\n",
       "        .tickValues([0.45, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, 0.506078431372549, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, 0.562156862745098, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, 0.6182352941176471, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, 0.6743137254901961, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, 0.7303921568627452, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, 0.7864705882352943, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, 0.8425490196078431, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, 0.8986274509803922, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, 0.9547058823529413, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;, &#x27;&#x27;]);\n",
       "\n",
       "    color_map_d2f09e21a56bd360e2497366a340a351.svg = d3.select(&quot;.legend.leaflet-control&quot;).append(&quot;svg&quot;)\n",
       "        .attr(&quot;id&quot;, &#x27;legend&#x27;)\n",
       "        .attr(&quot;width&quot;, 450)\n",
       "        .attr(&quot;height&quot;, 40);\n",
       "\n",
       "    color_map_d2f09e21a56bd360e2497366a340a351.g = color_map_d2f09e21a56bd360e2497366a340a351.svg.append(&quot;g&quot;)\n",
       "        .attr(&quot;class&quot;, &quot;key&quot;)\n",
       "        .attr(&quot;transform&quot;, &quot;translate(25,16)&quot;);\n",
       "\n",
       "    color_map_d2f09e21a56bd360e2497366a340a351.g.selectAll(&quot;rect&quot;)\n",
       "        .data(color_map_d2f09e21a56bd360e2497366a340a351.color.range().map(function(d, i) {\n",
       "          return {\n",
       "            x0: i ? color_map_d2f09e21a56bd360e2497366a340a351.x(color_map_d2f09e21a56bd360e2497366a340a351.color.domain()[i - 1]) : color_map_d2f09e21a56bd360e2497366a340a351.x.range()[0],\n",
       "            x1: i &lt; color_map_d2f09e21a56bd360e2497366a340a351.color.domain().length ? color_map_d2f09e21a56bd360e2497366a340a351.x(color_map_d2f09e21a56bd360e2497366a340a351.color.domain()[i]) : color_map_d2f09e21a56bd360e2497366a340a351.x.range()[1],\n",
       "            z: d\n",
       "          };\n",
       "        }))\n",
       "      .enter().append(&quot;rect&quot;)\n",
       "        .attr(&quot;height&quot;, 40 - 30)\n",
       "        .attr(&quot;x&quot;, function(d) { return d.x0; })\n",
       "        .attr(&quot;width&quot;, function(d) { return d.x1 - d.x0; })\n",
       "        .style(&quot;fill&quot;, function(d) { return d.z; });\n",
       "\n",
       "    color_map_d2f09e21a56bd360e2497366a340a351.g.call(color_map_d2f09e21a56bd360e2497366a340a351.xAxis).append(&quot;text&quot;)\n",
       "        .attr(&quot;class&quot;, &quot;caption&quot;)\n",
       "        .attr(&quot;y&quot;, 21)\n",
       "        .text(&quot;weights&quot;);\n",
       "&lt;/script&gt;\n",
       "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x1fe14a7baa0>"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_overlayed.explore('weights')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Remapping data\n",
    "\n",
    "The good thing of working with weights, is that you need simply to calculate them once and then you can apply them to any data.\n",
    "\n",
    "For this example, we will assign some values to the source grid and remap them to the target grid."
   ]
  },
  {
   "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>geometry</th>\n",
       "      <th>source_index</th>\n",
       "      <th>value1</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((0.00000 0.00000, 0.00000 1.00000, 1....</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>POLYGON ((0.00000 1.00000, 0.00000 2.00000, 1....</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>POLYGON ((1.00000 0.00000, 1.00000 1.00000, 2....</td>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>POLYGON ((1.00000 1.00000, 1.00000 2.00000, 2....</td>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>POLYGON ((2.00000 0.00000, 2.00000 1.00000, 3....</td>\n",
       "      <td>4</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>POLYGON ((2.00000 1.00000, 2.00000 2.00000, 3....</td>\n",
       "      <td>5</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry  source_index  value1\n",
       "0  POLYGON ((0.00000 0.00000, 0.00000 1.00000, 1....             0       1\n",
       "1  POLYGON ((0.00000 1.00000, 0.00000 2.00000, 1....             1       1\n",
       "2  POLYGON ((1.00000 0.00000, 1.00000 1.00000, 2....             2       1\n",
       "3  POLYGON ((1.00000 1.00000, 1.00000 2.00000, 2....             3       1\n",
       "4  POLYGON ((2.00000 0.00000, 2.00000 1.00000, 3....             4       1\n",
       "5  POLYGON ((2.00000 1.00000, 2.00000 2.00000, 3....             5       1"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Assign values to the grid\n",
    "import numpy as np\n",
    "\n",
    "grid_gdf['value1'] = 1\n",
    "\n",
    "grid_gdf"
   ]
  },
  {
   "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>source_index</th>\n",
       "      <th>target_index</th>\n",
       "      <th>geometry</th>\n",
       "      <th>area_overlay</th>\n",
       "      <th>area_source</th>\n",
       "      <th>weights</th>\n",
       "      <th>value1</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((0.00000 1.00000, 1.00000 1.00000, 1....</td>\n",
       "      <td>0.60</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.60</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((0.00000 1.40000, 1.00000 1.80000, 1....</td>\n",
       "      <td>0.60</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.60</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((1.00000 1.00000, 1.50000 1.00000, 1....</td>\n",
       "      <td>0.45</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.45</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 1.00000, 2.00000 0.00000, 1....</td>\n",
       "      <td>0.50</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.50</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((1.00000 1.80000, 1.50000 2.00000, 1....</td>\n",
       "      <td>0.45</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.45</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 2.00000, 2.00000 1.00000, 1....</td>\n",
       "      <td>0.50</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.50</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>4</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 1.00000, 3.00000 1.00000, 3....</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>5</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 2.00000, 3.00000 2.00000, 3....</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   source_index  target_index  \\\n",
       "0             0             0   \n",
       "1             1             0   \n",
       "2             2             0   \n",
       "3             2             1   \n",
       "4             3             0   \n",
       "5             3             1   \n",
       "6             4             1   \n",
       "7             5             1   \n",
       "\n",
       "                                            geometry  area_overlay  \\\n",
       "0  POLYGON ((0.00000 1.00000, 1.00000 1.00000, 1....          0.60   \n",
       "1  POLYGON ((0.00000 1.40000, 1.00000 1.80000, 1....          0.60   \n",
       "2  POLYGON ((1.00000 1.00000, 1.50000 1.00000, 1....          0.45   \n",
       "3  POLYGON ((2.00000 1.00000, 2.00000 0.00000, 1....          0.50   \n",
       "4  POLYGON ((1.00000 1.80000, 1.50000 2.00000, 1....          0.45   \n",
       "5  POLYGON ((2.00000 2.00000, 2.00000 1.00000, 1....          0.50   \n",
       "6  POLYGON ((2.00000 1.00000, 3.00000 1.00000, 3....          1.00   \n",
       "7  POLYGON ((2.00000 2.00000, 3.00000 2.00000, 3....          1.00   \n",
       "\n",
       "   area_source  weights  value1  \n",
       "0          1.0     0.60       1  \n",
       "1          1.0     0.60       1  \n",
       "2          1.0     0.45       1  \n",
       "3          1.0     0.50       1  \n",
       "4          1.0     0.45       1  \n",
       "5          1.0     0.50       1  \n",
       "6          1.0     1.00       1  \n",
       "7          1.0     1.00       1  "
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "col = 'value1'\n",
    "# Expand the input data to multiply it with the weights \n",
    "gdf_overlayed[col] = grid_gdf.loc[gdf_overlayed['source_index'], col].values\n",
    "gdf_overlayed\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "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>source_index</th>\n",
       "      <th>target_index</th>\n",
       "      <th>geometry</th>\n",
       "      <th>area_overlay</th>\n",
       "      <th>area_source</th>\n",
       "      <th>weights</th>\n",
       "      <th>value1</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((0.00000 1.00000, 1.00000 1.00000, 1....</td>\n",
       "      <td>0.60</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.60</td>\n",
       "      <td>0.60</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((0.00000 1.40000, 1.00000 1.80000, 1....</td>\n",
       "      <td>0.60</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.60</td>\n",
       "      <td>0.60</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((1.00000 1.00000, 1.50000 1.00000, 1....</td>\n",
       "      <td>0.45</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.45</td>\n",
       "      <td>0.45</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 1.00000, 2.00000 0.00000, 1....</td>\n",
       "      <td>0.50</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.50</td>\n",
       "      <td>0.50</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>POLYGON ((1.00000 1.80000, 1.50000 2.00000, 1....</td>\n",
       "      <td>0.45</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.45</td>\n",
       "      <td>0.45</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 2.00000, 2.00000 1.00000, 1....</td>\n",
       "      <td>0.50</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.50</td>\n",
       "      <td>0.50</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>4</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 1.00000, 3.00000 1.00000, 3....</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>5</td>\n",
       "      <td>1</td>\n",
       "      <td>POLYGON ((2.00000 2.00000, 3.00000 2.00000, 3....</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1.00</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   source_index  target_index  \\\n",
       "0             0             0   \n",
       "1             1             0   \n",
       "2             2             0   \n",
       "3             2             1   \n",
       "4             3             0   \n",
       "5             3             1   \n",
       "6             4             1   \n",
       "7             5             1   \n",
       "\n",
       "                                            geometry  area_overlay  \\\n",
       "0  POLYGON ((0.00000 1.00000, 1.00000 1.00000, 1....          0.60   \n",
       "1  POLYGON ((0.00000 1.40000, 1.00000 1.80000, 1....          0.60   \n",
       "2  POLYGON ((1.00000 1.00000, 1.50000 1.00000, 1....          0.45   \n",
       "3  POLYGON ((2.00000 1.00000, 2.00000 0.00000, 1....          0.50   \n",
       "4  POLYGON ((1.00000 1.80000, 1.50000 2.00000, 1....          0.45   \n",
       "5  POLYGON ((2.00000 2.00000, 2.00000 1.00000, 1....          0.50   \n",
       "6  POLYGON ((2.00000 1.00000, 3.00000 1.00000, 3....          1.00   \n",
       "7  POLYGON ((2.00000 2.00000, 3.00000 2.00000, 3....          1.00   \n",
       "\n",
       "   area_source  weights  value1  \n",
       "0          1.0     0.60    0.60  \n",
       "1          1.0     0.60    0.60  \n",
       "2          1.0     0.45    0.45  \n",
       "3          1.0     0.50    0.50  \n",
       "4          1.0     0.45    0.45  \n",
       "5          1.0     0.50    0.50  \n",
       "6          1.0     1.00    1.00  \n",
       "7          1.0     1.00    1.00  "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Multiply the values by the weights\n",
    "gdf_overlayed[col] *= gdf_overlayed['weights'] \n",
    "gdf_overlayed"
   ]
  },
  {
   "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>geometry</th>\n",
       "      <th>target_index</th>\n",
       "      <th>value1</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-1.00000 1.00000, 1.50000 0.00000, 1...</td>\n",
       "      <td>0</td>\n",
       "      <td>2.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>POLYGON ((1.50000 2.00000, 1.50000 0.00000, 3....</td>\n",
       "      <td>1</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry  target_index  value1\n",
       "0  POLYGON ((-1.00000 1.00000, 1.50000 0.00000, 1...             0     2.1\n",
       "1  POLYGON ((1.50000 2.00000, 1.50000 0.00000, 3....             1     3.0"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Sum the values belonging to the same target index\n",
    "gdf_out[col] = gdf_overlayed.groupby('target_index')[col].sum().values\n",
    "gdf_out"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we can try to have other values and see the results."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "grid_gdf['value2'] = np.arange(len(grid_gdf))\n",
    "grid_gdf['value3'] = np.arange(len(grid_gdf))**2\n",
    "\n",
    "for col in ['value2', 'value3']:\n",
    "    gdf_overlayed[col] = grid_gdf.loc[gdf_overlayed['source_index'], col].values\n",
    "    gdf_overlayed[col] *= gdf_overlayed['weights'] \n",
    "    gdf_out[col] = gdf_overlayed.groupby('target_index')[col].sum().values"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "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>geometry</th>\n",
       "      <th>target_index</th>\n",
       "      <th>value1</th>\n",
       "      <th>value2</th>\n",
       "      <th>value3</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>POLYGON ((-1.00000 1.00000, 1.50000 0.00000, 1...</td>\n",
       "      <td>0</td>\n",
       "      <td>2.1</td>\n",
       "      <td>2.85</td>\n",
       "      <td>6.45</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>POLYGON ((1.50000 2.00000, 1.50000 0.00000, 3....</td>\n",
       "      <td>1</td>\n",
       "      <td>3.0</td>\n",
       "      <td>11.50</td>\n",
       "      <td>47.50</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                            geometry  target_index  value1  \\\n",
       "0  POLYGON ((-1.00000 1.00000, 1.50000 0.00000, 1...             0     2.1   \n",
       "1  POLYGON ((1.50000 2.00000, 1.50000 0.00000, 3....             1     3.0   \n",
       "\n",
       "   value2  value3  \n",
       "0    2.85    6.45  \n",
       "1   11.50   47.50  "
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_out"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.12.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}