DownloadDataForDash.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Download Data for Dash app\n",
    "\n",
    "Needs png with transparent background and the coordinates, to be able to plug it in in the map."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib widget\n",
    "import matplotlib.pylab as plt\n",
    "import seaborn as sns\n",
    "\n",
    "import numpy as np\n",
    "\n",
    "import pandas as pd\n",
    "\n",
    "import re\n",
    "\n",
    "from tqdm import tqdm\n",
    "import glob\n",
    "import os\n",
    "\n",
    "from owslib.wms import WebMapService\n",
    "import rasterio"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Floading maps 2020\n",
    "\n",
    "For Denmark: https://gis.nst.dk/server/rest/services/ekstern\n",
    "\n",
    "The newest ones (Plantrin2) are divided in 14 zones. Although they have more scenarios (events happening in 1/20, 1/100 years and extreme) they do not include most of Denmark."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>zone</th>\n",
       "      <th>url</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Aabenraa</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Aab...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Esbjerg</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Esb...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Fredericia</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Fre...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Holstebro</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Hol...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Juelsminde</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Jue...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>Kobenhavn_NS</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Kob...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>Kobenhavn_Nord</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Koe...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>Kobenhavn_Syd</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Koe...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>Kolding</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Kol...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>Korsoer</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Kor...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>Koege</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Koe...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>Nyborg</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Nyb...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>Odense</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Ode...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>RandersFjord</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Ran...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>Skaerbaek</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Ska...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>Sydlolland</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Syd...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>Vejle</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Vej...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>Vordinborg</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Vor...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "              zone                                                url\n",
       "0         Aabenraa  https://gis.nst.dk/server/services/ekstern/Aab...\n",
       "1          Esbjerg  https://gis.nst.dk/server/services/ekstern/Esb...\n",
       "2       Fredericia  https://gis.nst.dk/server/services/ekstern/Fre...\n",
       "3        Holstebro  https://gis.nst.dk/server/services/ekstern/Hol...\n",
       "4       Juelsminde  https://gis.nst.dk/server/services/ekstern/Jue...\n",
       "5     Kobenhavn_NS  https://gis.nst.dk/server/services/ekstern/Kob...\n",
       "6   Kobenhavn_Nord  https://gis.nst.dk/server/services/ekstern/Koe...\n",
       "7    Kobenhavn_Syd  https://gis.nst.dk/server/services/ekstern/Koe...\n",
       "8          Kolding  https://gis.nst.dk/server/services/ekstern/Kol...\n",
       "9          Korsoer  https://gis.nst.dk/server/services/ekstern/Kor...\n",
       "10           Koege  https://gis.nst.dk/server/services/ekstern/Koe...\n",
       "11          Nyborg  https://gis.nst.dk/server/services/ekstern/Nyb...\n",
       "12          Odense  https://gis.nst.dk/server/services/ekstern/Ode...\n",
       "13    RandersFjord  https://gis.nst.dk/server/services/ekstern/Ran...\n",
       "14       Skaerbaek  https://gis.nst.dk/server/services/ekstern/Ska...\n",
       "15      Sydlolland  https://gis.nst.dk/server/services/ekstern/Syd...\n",
       "16           Vejle  https://gis.nst.dk/server/services/ekstern/Vej...\n",
       "17      Vordinborg  https://gis.nst.dk/server/services/ekstern/Vor..."
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_flood_maps = pd.read_csv('Plantrin_2.csv')\n",
    "df_flood_maps.head(20)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Quick fix. These two are not in WMS format. Needs to be fixed later\n",
    "df_flood_maps.drop([6,7], inplace=True);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Check what layers are available"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0                             Risikoområde ['EPSG:4326', 'EPSG:25832']\n",
      "3 Kulturarv_Seværdige_Fortidsminder_punkter ['EPSG:4326', 'EPSG:25832']\n",
      "4  Kulturarv_Seværdige_Fortidsminder_areal ['EPSG:4326', 'EPSG:25832']\n",
      "5                         Kulturarv_Kirker ['EPSG:4326', 'EPSG:25832']\n",
      "6              Kulturarv_Fredede_Bygninger ['EPSG:4326', 'EPSG:25832']\n",
      "7          Kulturarv_Fortidsminder_punkter ['EPSG:4326', 'EPSG:25832']\n",
      "8            Kulturarv_Fortidsminder_areal ['EPSG:4326', 'EPSG:25832']\n",
      "9                       Bes_sten_jorddiger ['EPSG:4326', 'EPSG:25832']\n",
      "10                       Begravelsesomraade ['EPSG:4326', 'EPSG:25832']\n",
      "12                         Vandvaerksboring ['EPSG:4326', 'EPSG:25832']\n",
      "13    Miljoeinteresser_Vandforsyningsboring ['EPSG:4326', 'EPSG:25832']\n",
      "14              Miljoeinteresser_Natura2000 ['EPSG:4326', 'EPSG:25832']\n",
      "15         Miljoeinteresser_Beskyttet_natur ['EPSG:4326', 'EPSG:25832']\n",
      "16    Miljoeinteresser_Aktiv_vandindvinding ['EPSG:4326', 'EPSG:25832']\n",
      "18             Risikovirksomheder_Holstebro ['EPSG:4326', 'EPSG:25832']\n",
      "19                  Forurenede_Virksomheder ['EPSG:4326', 'EPSG:25832']\n",
      "21                                Holstebro ['EPSG:4326', 'EPSG:25832']\n",
      "23                         Beredskab_Politi ['EPSG:4326', 'EPSG:25832']\n",
      "24                     Beredskab_Hospitaler ['EPSG:4326', 'EPSG:25832']\n",
      "25               Beredskab_Beredskabscentre ['EPSG:4326', 'EPSG:25832']\n",
      "27        Berørte_indbyggere_HOL_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "28         Berørte_indbyggere_HOL_MT50_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "29         Berørte_indbyggere_HOL_MT50_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "30        Berørte_indbyggere_HOL_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "31         Berørte_indbyggere_HOL_MT50_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "32         Berørte_indbyggere_HOL_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "34    Kritisk_Infrastruktur_Vandforsyninger ['EPSG:4326', 'EPSG:25832']\n",
      "35      Kritisk_Infrastruktur_Tranformation ['EPSG:4326', 'EPSG:25832']\n",
      "36     Kritisk_Infrastruktur_Rensningsanlæg ['EPSG:4326', 'EPSG:25832']\n",
      "37   Kritisk_Infrastruktur_Kraftvarmeværker ['EPSG:4326', 'EPSG:25832']\n",
      "38     Ejendomme_forsyninger_HOL_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "39      Ejendomme_forsyninger_HOL_MT50_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "40      Ejendomme_forsyninger_HOL_MT50_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "41     Ejendomme_forsyninger_HOL_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "42      Ejendomme_forsyninger_HOL_MT50_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "43      Ejendomme_forsyninger_HOL_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "46                    Husdyr_HOL_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "47                     Husdyr_HOL_MT50_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "48                     Husdyr_HOL_MT50_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "49                    Husdyr_HOL_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "50                     Husdyr_HOL_MT50_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "51                     Husdyr_HOL_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "53                  Afgrøder_HOL_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "54                   Afgrøder_HOL_MT50_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "55                   Afgrøder_HOL_MT50_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "56                  Afgrøder_HOL_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "57                   Afgrøder_HOL_MT50_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "58                   Afgrøder_HOL_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "60             Infrastruktur_HOL_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "61              Infrastruktur_HOL_MT50_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "62              Infrastruktur_HOL_MT50_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "63             Infrastruktur_HOL_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "64              Infrastruktur_HOL_MT50_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "65              Infrastruktur_HOL_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "67            Virksomhedstab_HOL_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "68             Virksomhedstab_HOL_MT50_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "69             Virksomhedstab_HOL_MT50_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "70            Virksomhedstab_HOL_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "71             Virksomhedstab_HOL_MT50_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "72             Virksomhedstab_HOL_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "74                IndboSkade_HOL_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "75                 IndboSkade_HOL_MT50_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "76                 IndboSkade_HOL_MT50_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "77                IndboSkade_HOL_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "78                 IndboSkade_HOL_MT50_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "79                 IndboSkade_HOL_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "81             Bygningsskade_HOL_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "82              Bygningsskade_HOL_MT50_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "83              Bygningsskade_HOL_MT50_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "84             Bygningsskade_HOL_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "85              Bygningsskade_HOL_MT50_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "86              Bygningsskade_HOL_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "89                     Dybde_HOL_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "90                      Dybde_HOL_MT50_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "91                      Dybde_HOL_MT50_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "92                     Dybde_HOL_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "93                      Dybde_HOL_MT50_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "94                      Dybde_HOL_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "96                Udbredelse_HOL_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "97                 Udbredelse_HOL_MT50_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "98                 Udbredelse_HOL_MT50_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "99                Udbredelse_HOL_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "100                 Udbredelse_HOL_MT50_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "101                 Udbredelse_HOL_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "103     Total_Økonomisk_skade_HOL_Mt100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "104      Total_Økonomisk_skade_HOL_MT50_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "105      Total_Økonomisk_skade_HOL_MT50_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "106     Total_Økonomisk_skade_HOL_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "107      Total_Økonomisk_skade_HOL_MT50_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "108      Total_Økonomisk_skade_HOL_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "110                    Risiko_HOL_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "111                     Risiko_HOL_MT50_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "112                     Risiko_HOL_MT50_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "113                    Risiko_HOL_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "114                     Risiko_HOL_MT50_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "115                     Risiko_HOL_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n"
     ]
    }
   ],
   "source": [
    "wms = WebMapService(df_flood_maps.iloc[3].url+'/WmsServer')\n",
    "\n",
    "for layer_nb in wms.contents:\n",
    "    print('%s %40s %s'%(layer_nb,wms[layer_nb].title,wms[layer_nb].crsOptions))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0                             Risikoområde ['EPSG:4326', 'EPSG:25832']\n",
      "3 Kulturarv_Seværdige_Fortidsminder_punkter ['EPSG:4326', 'EPSG:25832']\n",
      "4  Kulturarv_Seværdige_Fortidsminder_areal ['EPSG:4326', 'EPSG:25832']\n",
      "5                         Kulturarv_Kirker ['EPSG:4326', 'EPSG:25832']\n",
      "6              Kulturarv_Fredede_Bygninger ['EPSG:4326', 'EPSG:25832']\n",
      "7          Kulturarv_Fortidsminder_punkter ['EPSG:4326', 'EPSG:25832']\n",
      "8            Kulturarv_Fortidsminder_areal ['EPSG:4326', 'EPSG:25832']\n",
      "9                       Bes_sten_jorddiger ['EPSG:4326', 'EPSG:25832']\n",
      "10                       Begravelsesomraade ['EPSG:4326', 'EPSG:25832']\n",
      "12                         Vandvaerksboring ['EPSG:4326', 'EPSG:25832']\n",
      "13    Miljoeinteresser_Vandforsyningsboring ['EPSG:4326', 'EPSG:25832']\n",
      "14         Miljoeinteresser_Beskyttet_natur ['EPSG:4326', 'EPSG:25832']\n",
      "15    Miljoeinteresser_Aktiv_vandindvinding ['EPSG:4326', 'EPSG:25832']\n",
      "17        Aabenraa_forurenende_virksomheder ['EPSG:4326', 'EPSG:25832']\n",
      "18                       Risikovirksomheder ['EPSG:4326', 'EPSG:25832']\n",
      "20                                 Aabenraa ['EPSG:4326', 'EPSG:25832']\n",
      "22                         Beredskab_Politi ['EPSG:4326', 'EPSG:25832']\n",
      "23                     Beredskab_Hospitaler ['EPSG:4326', 'EPSG:25832']\n",
      "24               Beredskab_Beredskabscentre ['EPSG:4326', 'EPSG:25832']\n",
      "26        Berørte_indbyggere_AAB_MTEXT_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "27        Berørte_indbyggere_AAB_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "28        Berørte_indbyggere_AAB_MT100_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "29        Berørte_indbyggere_AAB_MTEXT_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "30        Berørte_indbyggere_AAB_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "31         Berørte_indbyggere_AAB_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "33     Ejendomme_forsyninger_AAB_MTEXT_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "34     Ejendomme_forsyninger_AAB_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "35     Ejendomme_forsyninger_AAB_MT100_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "36     Ejendomme_forsyninger_AAB_MTEXT_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "37     Ejendomme_forsyninger_AAB_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "38      Ejendomme_forsyninger_AAB_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "39    Kritisk_Infrastruktur_Vandforsyninger ['EPSG:4326', 'EPSG:25832']\n",
      "40      Kritisk_Infrastruktur_Tranformation ['EPSG:4326', 'EPSG:25832']\n",
      "41     Kritisk_Infrastruktur_Rensningsanlæg ['EPSG:4326', 'EPSG:25832']\n",
      "42   Kritisk_Infrastruktur_Kraftvarmeværker ['EPSG:4326', 'EPSG:25832']\n",
      "45                    Husdyr_AAB_MTEXT_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "46                    Husdyr_AAB_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "47                    Husdyr_AAB_MT100_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "48                    Husdyr_AAB_MTEXT_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "49                    Husdyr_AAB_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "50                     Husdyr_AAB_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "52                  Afgrøder_AAB_MTEXT_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "53                  Afgrøder_AAB_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "54                  Afgrøder_AAB_MT100_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "55                  Afgrøder_AAB_MTEXT_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "56                  Afgrøder_AAB_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "57                   Afgrøder_AAB_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "59             Infrastruktur_AAB_MTEXT_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "60             Infrastruktur_AAB_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "61             Infrastruktur_AAB_MT100_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "62             Infrastruktur_AAB_MTEXT_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "63             Infrastruktur_AAB_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "64              Infrastruktur_AAB_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "66            Virksomhedstab_AAB_MTEXT_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "67            Virksomhedstab_AAB_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "68            Virksomhedstab_AAB_MT100_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "69            Virksomhedstab_AAB_MTEXT_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "70            Virksomhedstab_AAB_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "71             Virksomhedstab_AAB_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "73                IndboSkade_AAB_MTEXT_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "74                IndboSkade_AAB_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "75                IndboSkade_AAB_MT100_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "76                IndboSkade_AAB_MTEXT_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "77                IndboSkade_AAB_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "78                 IndboSkade_AAB_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "80             Bygningsskade_AAB_MTEXT_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "81             Bygningsskade_AAB_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "82             Bygningsskade_AAB_MT100_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "83             Bygningsskade_AAB_MTEXT_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "84             Bygningsskade_AAB_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "85              Bygningsskade_AAB_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "88                     Dybde_AAB_MTEXT_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "89                     Dybde_AAB_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "90                     Dybde_AAB_MT100_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "91                     Dybde_AAB_MTEXT_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "92                     Dybde_AAB_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "93                      Dybde_AAB_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "95                Udbredelse_AAB_MTEXT_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "96                Udbredelse_AAB_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "97                Udbredelse_AAB_MT100_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "98                Udbredelse_AAB_MTEXT_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "99                Udbredelse_AAB_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "100                 Udbredelse_AAB_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "102     Total_økonomisk_skade_AAB_MTEXT_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "103     Total_økonomisk_skade_AAB_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "104     Total_økonomisk_skade_AAB_MT100_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "105     Total_økonomisk_skade_AAB_MTEXT_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "106     Total_økonomisk_skade_AAB_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "107      Total_økonomisk_skade_AAB_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "109                    Risiko_AAB_MTEXT_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "110                    Risiko_AAB_MT100_2115 ['EPSG:4326', 'EPSG:25832']\n",
      "111                    Risiko_AAB_MT100_2065 ['EPSG:4326', 'EPSG:25832']\n",
      "112                    Risiko_AAB_MTEXT_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "113                    Risiko_AAB_MT100_2019 ['EPSG:4326', 'EPSG:25832']\n",
      "114                     Risiko_AAB_MT20_2019 ['EPSG:4326', 'EPSG:25832']\n"
     ]
    }
   ],
   "source": [
    "wms = WebMapService(df_flood_maps.iloc[0].url+'/WmsServer')\n",
    "\n",
    "for layer_nb in wms.contents:\n",
    "    print('%s %40s %s'%(layer_nb,wms[layer_nb].title,wms[layer_nb].crsOptions))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We'll be interested in the *Udbredelse* layers for the diferent scenarios:\n",
    "\n",
    "    Udbredelse_AAB_MT20_2019\n",
    "    Udbredelse_AAB_MT100_2115\n",
    "    Udbredelse_AAB_MT100_2065\n",
    "    Udbredelse_AAB_MTEXT_2019\n",
    "    Udbredelse_AAB_MTEXT_2115\n",
    "    Udbredelse_AAB_MT100_2019\n",
    "    \n",
    "Since they can be in different layer numbers for the different zones, we'll check for each one"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "def check_layer_number(url,layer_name):\n",
    "    wms = WebMapService(url+'/WmsServer')\n",
    "    layer_titles = [wms[layer_nb].title for layer_nb in wms.contents]\n",
    "    layer_number = [layer_nb for layer_nb in wms.contents if re.match(layer_name,wms[layer_nb].title) ]\n",
    "    return layer_number"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "MTEXT_2115\n",
      "MTEXT_2019\n",
      "MT100_2115\n",
      "MT100_2065\n",
      "MT100_2019\n",
      "MT20_2019\n",
      "MT50_2115\n",
      "MT50_2065\n",
      "MT50_2019\n"
     ]
    }
   ],
   "source": [
    "scenarios = ['MTEXT_2115',\n",
    "             'MTEXT_2019',\n",
    "             'MT100_2115',\n",
    "             'MT100_2065',\n",
    "             'MT100_2019',\n",
    "             'MT20_2019',\n",
    "             'MT50_2115',\n",
    "             'MT50_2065',\n",
    "             'MT50_2019',\n",
    "            ]\n",
    "\n",
    "for scenario in scenarios:\n",
    "    print(scenario)\n",
    "    df_flood_maps['Udbredelse_'+scenario] = df_flood_maps.url.apply(lambda x: check_layer_number(x,'Udbredelse_\\w+_'+scenario))"
   ]
  },
  {
   "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>zone</th>\n",
       "      <th>url</th>\n",
       "      <th>Udbredelse_MTEXT_2115</th>\n",
       "      <th>Udbredelse_MTEXT_2019</th>\n",
       "      <th>Udbredelse_MT100_2115</th>\n",
       "      <th>Udbredelse_MT100_2065</th>\n",
       "      <th>Udbredelse_MT100_2019</th>\n",
       "      <th>Udbredelse_MT20_2019</th>\n",
       "      <th>Udbredelse_MT50_2115</th>\n",
       "      <th>Udbredelse_MT50_2065</th>\n",
       "      <th>Udbredelse_MT50_2019</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Aabenraa</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Aab...</td>\n",
       "      <td>[95]</td>\n",
       "      <td>[98]</td>\n",
       "      <td>[96]</td>\n",
       "      <td>[97]</td>\n",
       "      <td>[99]</td>\n",
       "      <td>[100]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Esbjerg</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Esb...</td>\n",
       "      <td>[97]</td>\n",
       "      <td>[100]</td>\n",
       "      <td>[98]</td>\n",
       "      <td>[99]</td>\n",
       "      <td>[101]</td>\n",
       "      <td>[102]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Fredericia</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Fre...</td>\n",
       "      <td>[92]</td>\n",
       "      <td>[95]</td>\n",
       "      <td>[93]</td>\n",
       "      <td>[94]</td>\n",
       "      <td>[96]</td>\n",
       "      <td>[97]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Holstebro</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Hol...</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[96]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[99]</td>\n",
       "      <td>[101]</td>\n",
       "      <td>[97]</td>\n",
       "      <td>[98]</td>\n",
       "      <td>[100]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Juelsminde</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Jue...</td>\n",
       "      <td>[95]</td>\n",
       "      <td>[98]</td>\n",
       "      <td>[96]</td>\n",
       "      <td>[97]</td>\n",
       "      <td>[99]</td>\n",
       "      <td>[100]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>Kobenhavn_NS</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Kob...</td>\n",
       "      <td>[95, 209]</td>\n",
       "      <td>[98, 212]</td>\n",
       "      <td>[96, 210]</td>\n",
       "      <td>[97, 211]</td>\n",
       "      <td>[99, 213]</td>\n",
       "      <td>[100, 214]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>Kolding</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Kol...</td>\n",
       "      <td>[104]</td>\n",
       "      <td>[108]</td>\n",
       "      <td>[105]</td>\n",
       "      <td>[106]</td>\n",
       "      <td>[109]</td>\n",
       "      <td>[110]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>Korsoer</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Kor...</td>\n",
       "      <td>[96]</td>\n",
       "      <td>[99]</td>\n",
       "      <td>[97]</td>\n",
       "      <td>[98]</td>\n",
       "      <td>[100]</td>\n",
       "      <td>[101]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>Koege</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Koe...</td>\n",
       "      <td>[95]</td>\n",
       "      <td>[98]</td>\n",
       "      <td>[96]</td>\n",
       "      <td>[97]</td>\n",
       "      <td>[99]</td>\n",
       "      <td>[100]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>Nyborg</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Nyb...</td>\n",
       "      <td>[97]</td>\n",
       "      <td>[100]</td>\n",
       "      <td>[98]</td>\n",
       "      <td>[99]</td>\n",
       "      <td>[101]</td>\n",
       "      <td>[102]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>Odense</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Ode...</td>\n",
       "      <td>[104]</td>\n",
       "      <td>[108]</td>\n",
       "      <td>[105]</td>\n",
       "      <td>[106]</td>\n",
       "      <td>[109]</td>\n",
       "      <td>[110]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>RandersFjord</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Ran...</td>\n",
       "      <td>[104]</td>\n",
       "      <td>[108]</td>\n",
       "      <td>[105]</td>\n",
       "      <td>[106]</td>\n",
       "      <td>[109]</td>\n",
       "      <td>[110]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>Skaerbaek</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Ska...</td>\n",
       "      <td>[91]</td>\n",
       "      <td>[94]</td>\n",
       "      <td>[92]</td>\n",
       "      <td>[93]</td>\n",
       "      <td>[95]</td>\n",
       "      <td>[96]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>Sydlolland</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Syd...</td>\n",
       "      <td>[114, 115]</td>\n",
       "      <td>[118, 119]</td>\n",
       "      <td>[116]</td>\n",
       "      <td>[117]</td>\n",
       "      <td>[120]</td>\n",
       "      <td>[121]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>Vejle</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Vej...</td>\n",
       "      <td>[103]</td>\n",
       "      <td>[107]</td>\n",
       "      <td>[104]</td>\n",
       "      <td>[105]</td>\n",
       "      <td>[108]</td>\n",
       "      <td>[109]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>Vordinborg</td>\n",
       "      <td>https://gis.nst.dk/server/services/ekstern/Vor...</td>\n",
       "      <td>[95]</td>\n",
       "      <td>[98]</td>\n",
       "      <td>[96]</td>\n",
       "      <td>[97]</td>\n",
       "      <td>[99]</td>\n",
       "      <td>[100]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "      <td>[]</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "            zone                                                url  \\\n",
       "0       Aabenraa  https://gis.nst.dk/server/services/ekstern/Aab...   \n",
       "1        Esbjerg  https://gis.nst.dk/server/services/ekstern/Esb...   \n",
       "2     Fredericia  https://gis.nst.dk/server/services/ekstern/Fre...   \n",
       "3      Holstebro  https://gis.nst.dk/server/services/ekstern/Hol...   \n",
       "4     Juelsminde  https://gis.nst.dk/server/services/ekstern/Jue...   \n",
       "5   Kobenhavn_NS  https://gis.nst.dk/server/services/ekstern/Kob...   \n",
       "8        Kolding  https://gis.nst.dk/server/services/ekstern/Kol...   \n",
       "9        Korsoer  https://gis.nst.dk/server/services/ekstern/Kor...   \n",
       "10         Koege  https://gis.nst.dk/server/services/ekstern/Koe...   \n",
       "11        Nyborg  https://gis.nst.dk/server/services/ekstern/Nyb...   \n",
       "12        Odense  https://gis.nst.dk/server/services/ekstern/Ode...   \n",
       "13  RandersFjord  https://gis.nst.dk/server/services/ekstern/Ran...   \n",
       "14     Skaerbaek  https://gis.nst.dk/server/services/ekstern/Ska...   \n",
       "15    Sydlolland  https://gis.nst.dk/server/services/ekstern/Syd...   \n",
       "16         Vejle  https://gis.nst.dk/server/services/ekstern/Vej...   \n",
       "17    Vordinborg  https://gis.nst.dk/server/services/ekstern/Vor...   \n",
       "\n",
       "   Udbredelse_MTEXT_2115 Udbredelse_MTEXT_2019 Udbredelse_MT100_2115  \\\n",
       "0                   [95]                  [98]                  [96]   \n",
       "1                   [97]                 [100]                  [98]   \n",
       "2                   [92]                  [95]                  [93]   \n",
       "3                     []                    []                  [96]   \n",
       "4                   [95]                  [98]                  [96]   \n",
       "5              [95, 209]             [98, 212]             [96, 210]   \n",
       "8                  [104]                 [108]                 [105]   \n",
       "9                   [96]                  [99]                  [97]   \n",
       "10                  [95]                  [98]                  [96]   \n",
       "11                  [97]                 [100]                  [98]   \n",
       "12                 [104]                 [108]                 [105]   \n",
       "13                 [104]                 [108]                 [105]   \n",
       "14                  [91]                  [94]                  [92]   \n",
       "15            [114, 115]            [118, 119]                 [116]   \n",
       "16                 [103]                 [107]                 [104]   \n",
       "17                  [95]                  [98]                  [96]   \n",
       "\n",
       "   Udbredelse_MT100_2065 Udbredelse_MT100_2019 Udbredelse_MT20_2019  \\\n",
       "0                   [97]                  [99]                [100]   \n",
       "1                   [99]                 [101]                [102]   \n",
       "2                   [94]                  [96]                 [97]   \n",
       "3                     []                  [99]                [101]   \n",
       "4                   [97]                  [99]                [100]   \n",
       "5              [97, 211]             [99, 213]           [100, 214]   \n",
       "8                  [106]                 [109]                [110]   \n",
       "9                   [98]                 [100]                [101]   \n",
       "10                  [97]                  [99]                [100]   \n",
       "11                  [99]                 [101]                [102]   \n",
       "12                 [106]                 [109]                [110]   \n",
       "13                 [106]                 [109]                [110]   \n",
       "14                  [93]                  [95]                 [96]   \n",
       "15                 [117]                 [120]                [121]   \n",
       "16                 [105]                 [108]                [109]   \n",
       "17                  [97]                  [99]                [100]   \n",
       "\n",
       "   Udbredelse_MT50_2115 Udbredelse_MT50_2065 Udbredelse_MT50_2019  \n",
       "0                    []                   []                   []  \n",
       "1                    []                   []                   []  \n",
       "2                    []                   []                   []  \n",
       "3                  [97]                 [98]                [100]  \n",
       "4                    []                   []                   []  \n",
       "5                    []                   []                   []  \n",
       "8                    []                   []                   []  \n",
       "9                    []                   []                   []  \n",
       "10                   []                   []                   []  \n",
       "11                   []                   []                   []  \n",
       "12                   []                   []                   []  \n",
       "13                   []                   []                   []  \n",
       "14                   []                   []                   []  \n",
       "15                   []                   []                   []  \n",
       "16                   []                   []                   []  \n",
       "17                   []                   []                   []  "
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_flood_maps.head(100)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Download data and save it as transparent png and get the coordinates"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "def download_image(url,layer,size,crs,dir_path):\n",
    "    \n",
    "    wms = WebMapService(url+'/WmsServer')\n",
    "    srs = wms[layer].crsOptions\n",
    "    bbox = wms[layer].boundingBoxWGS84 \n",
    "    output_name = wms[layer].title\n",
    "        \n",
    "    if crs not in srs:\n",
    "        print('CRS not available (available options:)',srs)\n",
    "        pass\n",
    "        \n",
    "    # Get the image from the server\n",
    "    img = wms.getmap(layers=[layer],\n",
    "                     styles=['default'],\n",
    "                     srs=crs,\n",
    "                     bbox=bbox,\n",
    "                     size=size,\n",
    "                     transparent=True,\n",
    "                     format='image/png8'\n",
    "                    )\n",
    "    \n",
    "    # Save the array as an image\n",
    "    print('Saved', dir_path+output_name+'.png')\n",
    "    with open(dir_path+output_name+'.png', 'wb') as out:\n",
    "        out.write(img.read())\n",
    "        \n",
    "    return dir_path+output_name+'.png',bbox"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "scenarios = ['MTEXT_2115',\n",
    "             'MTEXT_2019',\n",
    "             'MT100_2115',\n",
    "             'MT100_2065',\n",
    "             'MT100_2019',\n",
    "             'MT20_2019',\n",
    "             'MT50_2115',\n",
    "             'MT50_2065',\n",
    "             'MT50_2019',\n",
    "            ]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": true,
    "jupyter": {
     "outputs_hidden": true
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "  0%|          | 0/9 [00:00<?, ?it/s]"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_AAB_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_ESB_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_FRE_MTEXT_2115.png\n",
      "Scenario for  [] does not exist\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_JUE_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_KHS_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_KHN_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_KOL_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_KOR_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_KBU_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_NYB_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_ODE_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_RAF_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_SKB_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_SYD_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_NAK_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_VEJ_MTEXT_2115.png\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      " 11%|█         | 1/9 [00:34<04:36, 34.53s/it]"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved ForDash/Udbredelse_MTEXT_2115/Udbredelse_VOR_MTEXT_2115.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_AAB_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_ESB_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_FRE_MTEXT_2019.png\n",
      "Scenario for  [] does not exist\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_JUE_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_KHS_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_KHN_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_KOL_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_KOR_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_KBU_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_NYB_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_ODE_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_RAF_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_SKB_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_SYD_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_NAK_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_VEJ_MTEXT_2019.png\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      " 22%|██▏       | 2/9 [01:09<04:02, 34.62s/it]"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved ForDash/Udbredelse_MTEXT_2019/Udbredelse_VOR_MTEXT_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_AAB_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_ESB_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_FRE_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_HOL_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_JUE_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_KHS_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_KHN_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_KOL_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_KOR_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_KBU_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_NYB_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_ODE_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_RAF_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_SKB_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_SLO_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_VEJ_MT100_2115.png\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      " 33%|███▎      | 3/9 [01:44<03:29, 34.85s/it]"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved ForDash/Udbredelse_MT100_2115/Udbredelse_VOR_MT100_2115.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_AAB_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_ESB_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_FRE_MT100_2065.png\n",
      "Scenario for  [] does not exist\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_JUE_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_KHS_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_KHN_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_KOL_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_KOR_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_KBU_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_NYB_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_ODE_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_RAF_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_SKB_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_SLO_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_VEJ_MT100_2065.png\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      " 44%|████▍     | 4/9 [02:18<02:52, 34.46s/it]"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved ForDash/Udbredelse_MT100_2065/Udbredelse_VOR_MT100_2065.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_AAB_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_ESB_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_FRE_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_HOL_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_JUE_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_KHS_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_KHN_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_KOL_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_KOR_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_KBU_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_NYB_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_ODE_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_RAF_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_SKB_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_SLO_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_VEJ_MT100_2019.png\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      " 56%|█████▌    | 5/9 [02:52<02:17, 34.39s/it]"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved ForDash/Udbredelse_MT100_2019/Udbredelse_VOR_MT100_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_AAB_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_ESB_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_FRE_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_HOL_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_JUE_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_KHS_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_KHN_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_KOL_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_KOR_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_KBU_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_NYB_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_ODE_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_RAF_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_SKB_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_SLO_MT20_2019.png\n",
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_VEJ_MT20_2019.png\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      " 67%|██████▋   | 6/9 [03:24<01:40, 33.66s/it]"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved ForDash/Udbredelse_MT20_2019/Udbredelse_VOR_MT20_2019.png\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      " 78%|███████▊  | 7/9 [03:25<00:48, 24.02s/it]"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved ForDash/Udbredelse_MT50_2115/Udbredelse_HOL_MT50_2115.png\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      " 89%|████████▉ | 8/9 [03:27<00:17, 17.34s/it]"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved ForDash/Udbredelse_MT50_2065/Udbredelse_HOL_MT50_2065.png\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "100%|██████████| 9/9 [03:29<00:00, 23.26s/it]"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved ForDash/Udbredelse_MT50_2019/Udbredelse_HOL_MT50_2019.png\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n",
      "Scenario for  [] does not exist\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "\n"
     ]
    }
   ],
   "source": [
    "file_scenario_list = []\n",
    "file_name_list = []\n",
    "bbox_list = []\n",
    "\n",
    "for scenario in tqdm(scenarios):\n",
    "    \n",
    "    # Create path\n",
    "    dir_path = 'ForDash/Udbredelse_'+scenario+'/'\n",
    "    if not os.path.exists(dir_path):\n",
    "        os.makedirs(dir_path)\n",
    "        \n",
    "    # Create empty column for the bbox and image name to be updated during download\n",
    "    df_flood_maps['bbox_'+scenario] = 0\n",
    "    df_flood_maps['png_file_'+scenario] = 0\n",
    "    \n",
    "    for idx,row in df_flood_maps.iterrows():\n",
    "        if len(row['Udbredelse_'+scenario]) == 0:\n",
    "            print('Scenario for ',row['Udbredelse_'+scenario],'does not exist')\n",
    "            continue\n",
    "        else:\n",
    "            # there can be more than one layer (eg. cph north and south)\n",
    "            for layer in row['Udbredelse_'+scenario]:\n",
    "                file_name,bbox = download_image(row['url'],\n",
    "                                                  layer=layer,\n",
    "                                                  size=(2000,2000),\n",
    "                                                  crs='EPSG:4326',\n",
    "                                                  dir_path=dir_path)\n",
    "                file_name_list.append(file_name)\n",
    "                bbox_list.append(bbox)\n",
    "                file_scenario_list.append(scenario)\n",
    "                "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Data from 2019\n",
    "\n",
    "This is for only one scenario (1/1000 years floods) but includes all Denmark\n",
    "\n",
    "\n",
    "https://gis.nst.dk/server/services/ekstern/OSD_2016_2021_Plantrin1_Endelige_data/MapServer/WMSServer?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetCapabilities\n",
    "\n",
    " \n",
    "There are several layers in this one, but the two that covers all of Denmark for an extreme flooding is\n",
    "\n",
    " \n",
    "\n",
    "-          Fare fra hav (ekstrem hændelse)\n",
    "\n",
    "-          Fare fra vandløb (1000 år)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0                         Hovedvandoplande ['EPSG:4326', 'EPSG:25832']\n",
      "1                     Vandområdedistrekter ['EPSG:4326', 'EPSG:25832']\n",
      "2                           Kommunegrænser ['EPSG:4326', 'EPSG:25832']\n",
      "3      Potentielt Forurenende Virksomheder ['EPSG:4326', 'EPSG:25832']\n",
      "4                      Økonomisk Aktivitet ['EPSG:4326', 'EPSG:25832']\n",
      "5                                Kulturarv ['EPSG:4326', 'EPSG:25832']\n",
      "6                    Kritisk Infrastruktur ['EPSG:4326', 'EPSG:25832']\n",
      "7                            Infrastruktur ['EPSG:4326', 'EPSG:25832']\n",
      "8                                Beredskab ['EPSG:4326', 'EPSG:25832']\n",
      "9                               Befolkning ['EPSG:4326', 'EPSG:25832']\n",
      "10                          Arealanvendelse ['EPSG:4326', 'EPSG:25832']\n",
      "11               Fare fra vandløb (1000 år) ['EPSG:4326', 'EPSG:25832']\n",
      "12          Fare fra hav (ekstrem hændelse) ['EPSG:4326', 'EPSG:25832']\n",
      "13                         Samlet sårbarhed ['EPSG:4326', 'EPSG:25832']\n",
      "14                       Risiko for vandløb ['EPSG:4326', 'EPSG:25832']\n",
      "15                           Risiko for hav ['EPSG:4326', 'EPSG:25832']\n",
      "17                   Risikoområder_Aabenraa ['EPSG:4326', 'EPSG:25832']\n",
      "18                Risikoområder_Vordingborg ['EPSG:4326', 'EPSG:25832']\n",
      "19                      Risikoområder_Vejle ['EPSG:4326', 'EPSG:25832']\n",
      "20                 Risikoområder_Sydlolland ['EPSG:4326', 'EPSG:25832']\n",
      "21              Risikoområder_Randers_Fjord ['EPSG:4326', 'EPSG:25832']\n",
      "22               Risikoområder_Odense_Fjord ['EPSG:4326', 'EPSG:25832']\n",
      "23                     Risikoområder_Nyborg ['EPSG:4326', 'EPSG:25832']\n",
      "24                     Risikoområder_Korsor ['EPSG:4326', 'EPSG:25832']\n",
      "25                    Risikoområder_Kolding ['EPSG:4326', 'EPSG:25832']\n",
      "26        Risikoområder_Koge_Bugt_Kobenhavn ['EPSG:4326', 'EPSG:25832']\n",
      "27                 Risikoområder_Juelsminde ['EPSG:4326', 'EPSG:25832']\n",
      "28                  Risikoområder_Holstebro ['EPSG:4326', 'EPSG:25832']\n",
      "29                 Risikoområder_Fredericia ['EPSG:4326', 'EPSG:25832']\n",
      "30                    Risikoområder_Esbjerg ['EPSG:4326', 'EPSG:25832']\n",
      "31                  Risikoområder Forenklet ['EPSG:4326', 'EPSG:25832']\n"
     ]
    }
   ],
   "source": [
    "url_2019 = 'https://gis.nst.dk/server/services/ekstern/OSD_2016_2021_Plantrin1_Endelige_data/MapServer'\n",
    "wms_2019= WebMapService(url_2019+'/WMSServer')\n",
    "\n",
    "for layer_nb in wms_2019.contents:\n",
    "    print('%s %40s %s'%(layer_nb,wms_2019[layer_nb].title,wms_2019[layer_nb].crsOptions))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved ForDash/FloodingDK2019/Fare fra vandløb (1000 år).png\n",
      "Saved ForDash/FloodingDK2019/Fare fra hav (ekstrem hændelse).png\n"
     ]
    }
   ],
   "source": [
    "# Create path\n",
    "dir_path = 'ForDash/FloodingDK2019/'\n",
    "if not os.path.exists(dir_path):\n",
    "    os.makedirs(dir_path)\n",
    "\n",
    "file_name,bbox = download_image(url_2019,layer='11',size=(4000,4000),crs='EPSG:4326',dir_path='ForDash/FloodingDK2019/')\n",
    "file_scenario_list.append('2019')\n",
    "file_name_list.append(file_name)\n",
    "bbox_list.append(bbox)\n",
    "                                \n",
    "file_name,bbox = download_image(url_2019,layer='12',size=(4000,4000),crs='EPSG:4326',dir_path='ForDash/FloodingDK2019/')\n",
    "file_scenario_list.append('2019')\n",
    "file_name_list.append(file_name)\n",
    "bbox_list.append(bbox)                               "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Put all info in a dataframe"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "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>scenario</th>\n",
       "      <th>file</th>\n",
       "      <th>bbox</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>MTEXT_2115</td>\n",
       "      <td>ForDash/Udbredelse_MTEXT_2115/Udbredelse_KOL_M...</td>\n",
       "      <td>(9.377683, 55.47553, 9.553077, 55.510132)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>MTEXT_2115</td>\n",
       "      <td>ForDash/Udbredelse_MTEXT_2115/Udbredelse_KOR_M...</td>\n",
       "      <td>(11.090517, 55.307539, 11.217802, 55.363792)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>MTEXT_2115</td>\n",
       "      <td>ForDash/Udbredelse_MTEXT_2115/Udbredelse_KBU_M...</td>\n",
       "      <td>(12.120189, 55.410787, 12.468536, 55.639661)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>MTEXT_2115</td>\n",
       "      <td>ForDash/Udbredelse_MTEXT_2115/Udbredelse_NYB_M...</td>\n",
       "      <td>(10.780247, 55.284468, 10.857497, 55.321463)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>MTEXT_2115</td>\n",
       "      <td>ForDash/Udbredelse_MTEXT_2115/Udbredelse_ODE_M...</td>\n",
       "      <td>(10.356789, 55.390362, 10.653339, 55.57148)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>101</th>\n",
       "      <td>MT50_2115</td>\n",
       "      <td>ForDash/Udbredelse_MT50_2115/Udbredelse_HOL_MT...</td>\n",
       "      <td>(8.590139, 56.34769, 8.638287, 56.360013)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>102</th>\n",
       "      <td>MT50_2065</td>\n",
       "      <td>ForDash/Udbredelse_MT50_2065/Udbredelse_HOL_MT...</td>\n",
       "      <td>(8.590139, 56.34769, 8.638287, 56.360013)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>103</th>\n",
       "      <td>MT50_2019</td>\n",
       "      <td>ForDash/Udbredelse_MT50_2019/Udbredelse_HOL_MT...</td>\n",
       "      <td>(8.59072, 56.348062, 8.638312, 56.36004)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>104</th>\n",
       "      <td>2019</td>\n",
       "      <td>ForDash/FloodingDK2019/Fare fra vandløb (1000 ...</td>\n",
       "      <td>(8.019848, 54.445098, 15.562142, 57.653192)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>105</th>\n",
       "      <td>2019</td>\n",
       "      <td>ForDash/FloodingDK2019/Fare fra hav (ekstrem h...</td>\n",
       "      <td>(8.008285, 54.435388, 15.599735, 57.769964)</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>100 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       scenario                                               file  \\\n",
       "6    MTEXT_2115  ForDash/Udbredelse_MTEXT_2115/Udbredelse_KOL_M...   \n",
       "7    MTEXT_2115  ForDash/Udbredelse_MTEXT_2115/Udbredelse_KOR_M...   \n",
       "8    MTEXT_2115  ForDash/Udbredelse_MTEXT_2115/Udbredelse_KBU_M...   \n",
       "9    MTEXT_2115  ForDash/Udbredelse_MTEXT_2115/Udbredelse_NYB_M...   \n",
       "10   MTEXT_2115  ForDash/Udbredelse_MTEXT_2115/Udbredelse_ODE_M...   \n",
       "..          ...                                                ...   \n",
       "101   MT50_2115  ForDash/Udbredelse_MT50_2115/Udbredelse_HOL_MT...   \n",
       "102   MT50_2065  ForDash/Udbredelse_MT50_2065/Udbredelse_HOL_MT...   \n",
       "103   MT50_2019  ForDash/Udbredelse_MT50_2019/Udbredelse_HOL_MT...   \n",
       "104        2019  ForDash/FloodingDK2019/Fare fra vandløb (1000 ...   \n",
       "105        2019  ForDash/FloodingDK2019/Fare fra hav (ekstrem h...   \n",
       "\n",
       "                                             bbox  \n",
       "6       (9.377683, 55.47553, 9.553077, 55.510132)  \n",
       "7    (11.090517, 55.307539, 11.217802, 55.363792)  \n",
       "8    (12.120189, 55.410787, 12.468536, 55.639661)  \n",
       "9    (10.780247, 55.284468, 10.857497, 55.321463)  \n",
       "10    (10.356789, 55.390362, 10.653339, 55.57148)  \n",
       "..                                            ...  \n",
       "101     (8.590139, 56.34769, 8.638287, 56.360013)  \n",
       "102     (8.590139, 56.34769, 8.638287, 56.360013)  \n",
       "103      (8.59072, 56.348062, 8.638312, 56.36004)  \n",
       "104   (8.019848, 54.445098, 15.562142, 57.653192)  \n",
       "105   (8.008285, 54.435388, 15.599735, 57.769964)  \n",
       "\n",
       "[100 rows x 3 columns]"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_png = pd.DataFrame({'scenario':file_scenario_list,\n",
    "                       'file':file_name_list,\n",
    "                       'bbox':bbox_list})\n",
    "df_png.to_csv('png_file_list')\n",
    "df_png.tail(100)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Data from Sweden\n",
    "\n",
    "Maps are available here: https://gisapp.msb.se/Apps/oversvamningsportal/hemta-data.html\n",
    "\n",
    "We'll be using the 200-årsflöde (Klimatanpassat flöde för slutet av seklet) and 100-årsflöde (Klimatanpassat flöde för slutet av seklet)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Check which layers are available"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0      Översiktliga, beräknat högsta flöde ['EPSG:3006', 'EPSG:4326']\n",
      "1               Översiktliga, 100-årsflöde ['EPSG:3006', 'EPSG:4326']\n",
      "2  Tvärsektioner, översiktliga karteringar ['EPSG:3006', 'EPSG:4326']\n",
      "3         Torne älv, beräknat högsta flöde ['EPSG:3006', 'EPSG:4326']\n",
      "4                  Torne älv, 250-årsflöde ['EPSG:3006', 'EPSG:4326']\n",
      "5                  Torne älv, 100-årsflöde ['EPSG:3006', 'EPSG:4326']\n",
      "6                 Torne älv, tvärsektioner ['EPSG:3006', 'EPSG:4326']\n",
      "7 Göta älv, 1400 m³/s (Motsvarande beräknat högsta flöde) ['EPSG:3006', 'EPSG:4326']\n",
      "8 Göta älv, 1200 m³/s (Motsvarande 200-årsflöde) ['EPSG:3006', 'EPSG:4326']\n",
      "9 Göta älv, 1030 m³/s (Motsvarande 100-årsflöde) ['EPSG:3006', 'EPSG:4326']\n",
      "10                  Göta älv, tvärsektioner ['EPSG:3006', 'EPSG:4326']\n",
      "11                    Beräknat högsta flöde ['EPSG:3006', 'EPSG:4326']\n",
      "12 200-årsflöde (Klimatanpassat flöde för slutet av seklet) ['EPSG:3006', 'EPSG:4326']\n",
      "13 100-årsflöde (Klimatanpassat flöde för slutet av seklet) ['EPSG:3006', 'EPSG:4326']\n",
      "14                            Tvärsektioner ['EPSG:3006', 'EPSG:4326']\n",
      "15 Översvämningskarterade vattendrag, översikt ['EPSG:3006', 'EPSG:4326']\n"
     ]
    }
   ],
   "source": [
    "url_swe = 'https://gisapp.msb.se/arcgis/services/Oversvamningskarteringar/karteringar/MapServer'\n",
    "wms_swe= WebMapService(url_swe+'/WMSServer')\n",
    "\n",
    "for layer_nb in wms_swe.contents:\n",
    "    print('%s %40s %s'%(layer_nb,wms[layer_nb].title,wms[layer_nb].crsOptions))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create path\n",
    "dir_path = 'ForDash/FloodingSWE2020/'\n",
    "if not os.path.exists(dir_path):\n",
    "    os.makedirs(dir_path)\n",
    "\n",
    "# save info\n",
    "file_scenario_list = []\n",
    "file_name_list = []\n",
    "bbox_list = []\n",
    "\n",
    "file_name,bbox = download_image(url_swe,layer='12',size=(4000,4000),crs='EPSG:4326',dir_path='ForDash/FloodingSWE2020/')\n",
    "file_scenario_list.append('200-years')\n",
    "file_name_list.append(file_name)\n",
    "bbox_list.append(bbox)\n",
    "                                \n",
    "file_name,bbox = download_image(url_swe,layer='13',size=(4000,4000),crs='EPSG:4326',dir_path='ForDash/FloodingSWE2020/')\n",
    "file_scenario_list.append('100-years')\n",
    "file_name_list.append(file_name)\n",
    "bbox_list.append(bbox)\n",
    "\n",
    "df_png = pd.DataFrame({'scenario':file_scenario_list,\n",
    "                       'file':file_name_list,\n",
    "                       'bbox':bbox_list})\n",
    "df_png.to_csv('png_file_list_swe')\n",
    "df_png.tail(100)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved ForDash/FloodingSWE2020_EPSG3006/200-årsflöde (Klimatanpassat flöde för slutet av seklet).png\n",
      "Saved ForDash/FloodingSWE2020_EPSG3006/100-årsflöde (Klimatanpassat flöde för slutet av seklet).png\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>scenario</th>\n",
       "      <th>file</th>\n",
       "      <th>bbox</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>200-years</td>\n",
       "      <td>ForDash/FloodingSWE2020_EPSG3006/200-årsflöde ...</td>\n",
       "      <td>(10.424994, 55.472992, 23.711993, 67.43329)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>100-years</td>\n",
       "      <td>ForDash/FloodingSWE2020_EPSG3006/100-årsflöde ...</td>\n",
       "      <td>(10.255073, 55.472992, 24.032033, 68.288574)</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    scenario                                               file  \\\n",
       "0  200-years  ForDash/FloodingSWE2020_EPSG3006/200-årsflöde ...   \n",
       "1  100-years  ForDash/FloodingSWE2020_EPSG3006/100-årsflöde ...   \n",
       "\n",
       "                                           bbox  \n",
       "0   (10.424994, 55.472992, 23.711993, 67.43329)  \n",
       "1  (10.255073, 55.472992, 24.032033, 68.288574)  "
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Create path\n",
    "dir_path = 'ForDash/FloodingSWE2020_EPSG3006/'\n",
    "if not os.path.exists(dir_path):\n",
    "    os.makedirs(dir_path)\n",
    "\n",
    "# save info\n",
    "file_scenario_list = []\n",
    "file_name_list = []\n",
    "bbox_list = []\n",
    "\n",
    "file_name,bbox = download_image(url_swe,layer='12',size=(4000,4000),crs='EPSG:3006',dir_path='ForDash/FloodingSWE2020_EPSG3006/')\n",
    "file_scenario_list.append('200-years')\n",
    "file_name_list.append(file_name)\n",
    "bbox_list.append(bbox)\n",
    "                                \n",
    "file_name,bbox = download_image(url_swe,layer='13',size=(4000,4000),crs='EPSG:3006',dir_path='ForDash/FloodingSWE2020_EPSG3006/')\n",
    "file_scenario_list.append('100-years')\n",
    "file_name_list.append(file_name)\n",
    "bbox_list.append(bbox)\n",
    "\n",
    "df_png = pd.DataFrame({'scenario':file_scenario_list,\n",
    "                       'file':file_name_list,\n",
    "                       'bbox':bbox_list})\n",
    "df_png.to_csv('png_file_list_swe_epsg3006')\n",
    "df_png.tail(100)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}