make_rida.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Creating municipality RIDA layers\n",
    "\n",
    "Assuming that the data we received from Septima by WeTransfer (2 folders with fgb files, plus one separate fgb file) is in the folder `data/RIDA` (from Dropbox [here](https://www.dropbox.com/scl/fo/38gkcvyfex0atdbqe3hfn/h?rlkey=z9xwq2mh6ga03c301b7f66wqx&dl=0))\n",
    "\n",
    "This script imports each of the fgb layers we need and cuts them into muncipality-sized bites; then filters merges them to create our evaluation layers.\n",
    "\n",
    "Run from main repo folder!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "# import packages\n",
    "import geopandas as gpd\n",
    "import pandas as pd\n",
    "import os\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "# read in kommuner data and make subfolders\n",
    "muni = gpd.read_file(\"data/raw/municipality_boundaries/muni_boundary.gpkg\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "# create subfolders for each municipality\n",
    "for _, row in muni.iterrows():\n",
    "    os.makedirs(f\"data/raw/municipality_data/{row.kommunekode}\", exist_ok=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "# dictionary: key is fgb file name (origin file), value is gpkg file name (destination file)\n",
    "rida_dict = {\n",
    "    \"facilit_indkoeb.fgb\": \"indkoeb.gpkg\",\n",
    "    \"facilit_overnatning.fgb\": \"overnatning.gpkg\",\n",
    "    \"facilit_rasteplads.fgb\": \"rasteplads.gpkg\",\n",
    "    'facilit_service.fgb': \"infoservice.gpkg\",\n",
    "    'land_anvendelse.fgb': \"anvendelse.gpkg\", # new! formerly 5 layers (byomrade, dyrket areal, naturareal, skovinddeling, teknisk areal)\n",
    "    'land_beskyttnatur_flade.fgb': \"beskyttet_natur.gpkg\",\n",
    "    'land_beskyttnatur_linje.fgb': \"beskyttet_vandloeb.gpkg\",\n",
    "    'land_besoegs.fgb': \"besoeg.gpkg\",\n",
    "    \"land_fortid.fgb\": \"fortidsminder.gpkg\",\n",
    "    'land_frednatpark.fgb': \"frednatpark.gpkg\", # new!\n",
    "    'land_landemaerke.fgb': \"landemaerker.gpkg\",\n",
    "    'land_udflugt.fgb': \"udflugt.gpkg\",\n",
    "    'land_vaerdifuld.fgb': \"vaerdifuld.gpkg\", # new!\n",
    "    'facilit_overnatning.fgb': \"overnatning.gpkg\",    \n",
    "}\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "reading land_beskyttnatur_linje.fgb\n",
      "reading land_besoegs.fgb\n",
      "reading land_fortid.fgb\n",
      "reading land_frednatpark.fgb\n",
      "reading land_landemaerke.fgb\n",
      "reading land_udflugt.fgb\n",
      "reading land_vaerdifuld.fgb\n",
      "reading facilit_overnatning.fgb\n",
      "CPU times: user 7min 30s, sys: 15.2 s, total: 7min 45s\n",
      "Wall time: 7min 53s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "# read in fgb file; \n",
    "# for each of the kommuner, clip fgb file to its extent; \n",
    "# then save in corresponding folder\n",
    "for k, v in rida_dict.items():\n",
    "    print(f\"reading {k}\")\n",
    "    gdf = gpd.read_file(f\"data/municipality_data/{k}\")\n",
    "    assert gdf.crs == muni.crs\n",
    "    for _, row in muni.iterrows():\n",
    "        gdf_curr = gdf.copy()\n",
    "        gdf_curr = gdf_curr.clip(row.geometry)\n",
    "        gdf_curr = gdf_curr.reset_index(drop=True)\n",
    "        if not gdf_curr.empty:\n",
    "            gdf_curr.to_file(f\"data/raw/municipality_data/{row.kommunekode}/{v}\")    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "def merge_gdfs(gdf_list):\n",
    "    # make sure all gdfs are the same crs\n",
    "    assert len(set([gdf.crs for gdf in gdf_list])) == 1\n",
    "    # make sure we have the same geometries everywhere\n",
    "    assert len(set([t for gdf in gdf_list for t in gdf.type])) == 1\n",
    "    # concatenate with pandas\n",
    "    gdf_main = pd.concat([gdf[[\"gruppe\", \"type\", \"geometry\"]].copy() for gdf in gdf_list])\n",
    "    gdf_main = gdf_main.reset_index(drop=True)\n",
    "    return gdf_main"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "munis_remove = [\n",
    "    \"0161\", \"0240\", \"0163\", \"0411\", # no besoeg\n",
    "    \"0269\", \"0159\", # no overnatning\n",
    "    \"0183\", \"0147\", \"0190\", \"0187\", \"0825\", # no fortidsminder\n",
    "    \"0157\", \"0175\" # no udflugt\n",
    "]\n",
    "\n",
    "munis = list(muni.kommunekode)\n",
    "for m in munis_remove:\n",
    "    munis.remove(m)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "start kode 0370\n",
      "\t data read in\n",
      "\t 0370 done\n",
      "start kode 0330\n",
      "\t data read in\n",
      "\t 0330 done\n",
      "start kode 0766\n",
      "\t data read in\n",
      "\t 0766 done\n",
      "start kode 0510\n",
      "\t data read in\n",
      "\t 0510 done\n",
      "start kode 0550\n",
      "\t data read in\n",
      "\t 0550 done\n",
      "start kode 0751\n",
      "\t data read in\n",
      "\t 0751 done\n",
      "start kode 0165\n",
      "\t data read in\n",
      "\t 0165 done\n",
      "start kode 0621\n",
      "\t data read in\n",
      "\t 0621 done\n",
      "start kode 0607\n",
      "\t data read in\n",
      "\t 0607 done\n",
      "start kode 0259\n",
      "\t data read in\n",
      "\t 0259 done\n",
      "start kode 0630\n",
      "\t data read in\n",
      "\t 0630 done\n",
      "start kode 0787\n",
      "\t data read in\n",
      "\t 0787 done\n",
      "start kode 0326\n",
      "\t data read in\n",
      "\t 0326 done\n",
      "start kode 0253\n",
      "\t data read in\n",
      "\t 0253 done\n",
      "start kode 0101\n",
      "\t data read in\n",
      "\t 0101 done\n",
      "start kode 0760\n",
      "\t data read in\n",
      "\t 0760 done\n",
      "start kode 0820\n",
      "\t data read in\n",
      "\t 0820 done\n",
      "start kode 0791\n",
      "\t data read in\n",
      "\t 0791 done\n",
      "start kode 0846\n",
      "\t data read in\n",
      "\t 0846 done\n",
      "start kode 0840\n",
      "\t data read in\n",
      "\t 0840 done\n",
      "start kode 0153\n",
      "\t data read in\n",
      "\t 0153 done\n",
      "start kode 0400\n",
      "\t data read in\n",
      "\t 0400 done\n",
      "start kode 0185\n",
      "\t data read in\n",
      "\t 0185 done\n",
      "start kode 0155\n",
      "\t data read in\n",
      "\t 0155 done\n",
      "start kode 0740\n",
      "\t data read in\n",
      "\t 0740 done\n",
      "start kode 0561\n",
      "\t data read in\n",
      "\t 0561 done\n",
      "start kode 0410\n",
      "\t data read in\n",
      "\t 0410 done\n",
      "start kode 0480\n",
      "\t data read in\n",
      "\t 0480 done\n",
      "start kode 0756\n",
      "\t data read in\n",
      "\t 0756 done\n",
      "start kode 0615\n",
      "\t data read in\n",
      "\t 0615 done\n",
      "start kode 0746\n",
      "\t data read in\n",
      "\t 0746 done\n",
      "start kode 0173\n",
      "\t data read in\n",
      "\t 0173 done\n",
      "start kode 0223\n",
      "\t data read in\n",
      "\t 0223 done\n",
      "start kode 0575\n",
      "\t data read in\n",
      "\t 0575 done\n",
      "start kode 0265\n",
      "\t data read in\n",
      "\t 0265 done\n",
      "start kode 0706\n",
      "\t data read in\n",
      "\t 0706 done\n",
      "start kode 0169\n",
      "\t data read in\n",
      "\t 0169 done\n",
      "start kode 0707\n",
      "\t data read in\n",
      "\t 0707 done\n",
      "start kode 0580\n",
      "\t data read in\n",
      "\t 0580 done\n",
      "start kode 0727\n",
      "\t data read in\n",
      "\t 0727 done\n",
      "start kode 0492\n",
      "\t data read in\n",
      "\t 0492 done\n",
      "start kode 0482\n",
      "\t data read in\n",
      "\t 0482 done\n",
      "start kode 0440\n",
      "\t data read in\n",
      "\t 0440 done\n",
      "start kode 0461\n",
      "\t data read in\n",
      "\t 0461 done\n",
      "start kode 0306\n",
      "\t data read in\n",
      "\t 0306 done\n",
      "start kode 0665\n",
      "\t data read in\n",
      "\t 0665 done\n",
      "start kode 0260\n",
      "\t data read in\n",
      "\t 0260 done\n",
      "start kode 0151\n",
      "\t data read in\n",
      "\t 0151 done\n",
      "start kode 0360\n",
      "\t data read in\n",
      "\t 0360 done\n",
      "start kode 0741\n",
      "\t data read in\n",
      "\t 0741 done\n",
      "start kode 0773\n",
      "\t data read in\n",
      "\t 0773 done\n",
      "start kode 0813\n",
      "\t data read in\n",
      "\t 0813 done\n",
      "start kode 0450\n",
      "\t data read in\n",
      "\t 0450 done\n",
      "start kode 0376\n",
      "\t data read in\n",
      "\t 0376 done\n",
      "start kode 0540\n",
      "\t data read in\n",
      "\t 0540 done\n",
      "start kode 0657\n",
      "\t data read in\n",
      "\t 0657 done\n",
      "start kode 0219\n",
      "\t data read in\n",
      "\t 0219 done\n",
      "start kode 0810\n",
      "\t data read in\n",
      "\t 0810 done\n",
      "start kode 0563\n",
      "\t data read in\n",
      "\t 0563 done\n",
      "start kode 0336\n",
      "\t data read in\n",
      "\t 0336 done\n",
      "start kode 0270\n",
      "\t data read in\n",
      "\t 0270 done\n",
      "start kode 0167\n",
      "\t data read in\n",
      "\t 0167 done\n",
      "start kode 0340\n",
      "\t data read in\n",
      "\t 0340 done\n",
      "start kode 0661\n",
      "\t data read in\n",
      "\t 0661 done\n",
      "start kode 0201\n",
      "\t data read in\n",
      "\t 0201 done\n",
      "start kode 0316\n",
      "\t data read in\n",
      "\t 0316 done\n",
      "start kode 0573\n",
      "\t data read in\n",
      "\t 0573 done\n",
      "start kode 0217\n",
      "\t data read in\n",
      "\t 0217 done\n",
      "start kode 0860\n",
      "\t data read in\n",
      "\t 0860 done\n",
      "start kode 0479\n",
      "\t data read in\n",
      "\t 0479 done\n",
      "start kode 0730\n",
      "\t data read in\n",
      "\t 0730 done\n",
      "start kode 0710\n",
      "\t data read in\n",
      "\t 0710 done\n",
      "start kode 0530\n",
      "\t data read in\n",
      "\t 0530 done\n",
      "start kode 0420\n",
      "\t data read in\n",
      "\t 0420 done\n",
      "start kode 0320\n",
      "\t data read in\n",
      "\t 0320 done\n",
      "start kode 0350\n",
      "\t data read in\n",
      "\t 0350 done\n",
      "start kode 0849\n",
      "\t data read in\n",
      "\t 0849 done\n",
      "start kode 0430\n",
      "\t data read in\n",
      "\t 0430 done\n",
      "start kode 0329\n",
      "\t data read in\n",
      "\t 0329 done\n",
      "start kode 0671\n",
      "\t data read in\n",
      "\t 0671 done\n",
      "start kode 0230\n",
      "\t data read in\n",
      "\t 0230 done\n",
      "start kode 0250\n",
      "\t data read in\n",
      "\t 0250 done\n",
      "start kode 0210\n",
      "\t data read in\n",
      "\t 0210 done\n",
      "start kode 0851\n",
      "\t data read in\n",
      "\t 0851 done\n",
      "start kode 0779\n",
      "\t data read in\n",
      "\t 0779 done\n",
      "start kode 0390\n",
      "\t data read in\n",
      "\t 0390 done\n"
     ]
    }
   ],
   "source": [
    "# evaluation layers:\n",
    "# read in all fgb-to-gpkg layers;\n",
    "# drop not needed columns;\n",
    "# explode;\n",
    "# drop not needed geometries;\n",
    "# filter and merge;\n",
    "# save in corresponding kommune subfolder\n",
    " \n",
    "for kode in munis:\n",
    "    \n",
    "    print(f\"start kode {kode}\")\n",
    "\n",
    "    anve = gpd.read_file(f\"data/raw/municipality_data/{kode}/anvendelse.gpkg\")\n",
    "    anve = anve[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    anve = anve.explode(index_parts=False)\n",
    "    anve = anve[anve.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    besk = gpd.read_file(f\"data/raw/municipality_data/{kode}/beskyttet_natur.gpkg\")\n",
    "    besk = besk[[\"type\", \"geometry\"]]\n",
    "    besk[\"gruppe\"] = \"beskyttet_natur\"\n",
    "    besk = besk.explode(index_parts=False)\n",
    "    besk = besk[besk.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    beso = gpd.read_file(f\"data/raw/municipality_data/{kode}/besoeg.gpkg\")\n",
    "    beso = beso[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    beso = beso.explode(index_parts=False)\n",
    "    beso = beso[beso.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fort = gpd.read_file(f\"data/raw/municipality_data/{kode}/fortidsminder.gpkg\")\n",
    "    fort = fort[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    fort = fort.explode(index_parts=False)\n",
    "    fort = fort[fort.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fred = gpd.read_file(f\"data/raw/municipality_data/{kode}/frednatpark.gpkg\")\n",
    "    fred = fred[[\"type\", \"geometry\"]]\n",
    "    fred[\"gruppe\"] = \"frednatpark\"\n",
    "    fred = fred.explode(index_parts=False)\n",
    "    fred = fred[fred.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    indk = gpd.read_file(f\"data/raw/municipality_data/{kode}/indkoeb.gpkg\")\n",
    "    indk = indk[[\"type\", \"geometry\"]]\n",
    "    indk[\"gruppe\"] = \"indkoeb\"\n",
    "    indk = indk.explode(index_parts=False)\n",
    "    indk = indk[indk.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    info = gpd.read_file(f\"data/raw/municipality_data/{kode}/infoservice.gpkg\")\n",
    "    info = info[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    info = info.explode(index_parts=False)\n",
    "    info = info[info.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    land = gpd.read_file(f\"data/raw/municipality_data/{kode}/landemaerker.gpkg\")\n",
    "    land = land[[\"type\", \"geometry\"]]\n",
    "    land[\"gruppe\"] = \"landemaerker\"\n",
    "    land = land.explode(index_parts=False)\n",
    "    land = land[land.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    over = gpd.read_file(f\"data/raw/municipality_data/{kode}/overnatning.gpkg\")\n",
    "    over = over[[\"type\", \"geometry\"]]\n",
    "    over[\"gruppe\"] = \"overnatning\"\n",
    "    over = over.explode(index_parts=False)\n",
    "    over = over[over.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    rast = gpd.read_file(f\"data/raw/municipality_data/{kode}/rasteplads.gpkg\")\n",
    "    rast = rast[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    rast = rast.explode(index_parts=False)\n",
    "    rast = rast[rast.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    udfl = gpd.read_file(f\"data/raw/municipality_data/{kode}/udflugt.gpkg\")\n",
    "    udfl = udfl[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    udfl = udfl.explode(index_parts=False)\n",
    "    udfl = udfl[udfl.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    vaer = gpd.read_file(f\"data/raw/municipality_data/{kode}/vaerdifuld.gpkg\")\n",
    "    vaer = vaer[[\"type\", \"geometry\"]]\n",
    "    vaer[\"gruppe\"] = \"vaerdifuld\"\n",
    "    vaer = vaer.explode(index_parts=False)\n",
    "    vaer = vaer[vaer.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    print(\"\\t data read in\")\n",
    "\n",
    "    nature = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Naturtyper\") & -anve[\"type\"].isin([\"Råstofområde\", \"Gartneri\"])].copy(),\n",
    "            anve[anve[\"gruppe\"]==\"Skovinddeling\"].copy(),\n",
    "            besk.copy(),\n",
    "            vaer[vaer[\"type\"]==\"Naturbeskyttelsesområde\"].copy(),\n",
    "            fred.copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(nature) > 0:\n",
    "        nature.to_file(f\"data/raw/municipality_data/{kode}/nature.gpkg\", index = False)\n",
    "        \n",
    "    culture = merge_gdfs(\n",
    "        [\n",
    "            anve[anve[\"type\"].isin(['Bykerne'])].copy(),\n",
    "            vaer[vaer[\"type\"].isin(\n",
    "                [\n",
    "                    'Kulturhistorisk bevaringsværdi', \n",
    "                    'Værdifuldt kulturmiljø', \n",
    "                    'Bevaringsværdig landskab'\n",
    "                ]\n",
    "            )].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(culture) > 0:\n",
    "        culture.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/culture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    sommerhus = anve[\n",
    "        anve[\"type\"].isin(\n",
    "            ['Sommerhusområde', 'Sommerhusområde skov']\n",
    "            )\n",
    "        ].copy().reset_index(drop=True)\n",
    "    if len(sommerhus) > 0:\n",
    "        sommerhus.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/sommerhus.gpkg\", \n",
    "            index = False)\n",
    "\n",
    "    agri = anve[\n",
    "        anve[\"gruppe\"]==\"Dyrkede arealer\"].copy().reset_index(\n",
    "            drop=True)\n",
    "\n",
    "    if len(agri) > 0:\n",
    "        agri.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/agriculture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    verify = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Teknisk areal\") & (-anve[\"type\"].isin([\"Sportsanlæg\"]))].copy(),\n",
    "            anve[(anve[\"gruppe\"]==\"Bymæssig anvendelse\") & (anve[\"type\"].isin([\"Høj bebyggelse\", \"Erhverv\"]))].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(verify) > 0:\n",
    "        verify.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/verify.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "\n",
    "    pois = merge_gdfs(\n",
    "        [beso, fort, udfl, land]\n",
    "    )\n",
    "\n",
    "    if len(pois) > 0:\n",
    "        pois.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/pois.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    faci = merge_gdfs(\n",
    "        [\n",
    "            info[-info[\"type\"].isin([\"Turistkontor\"])].copy(),\n",
    "            rast\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(faci) > 0:\n",
    "        faci.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/facilities.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    service = merge_gdfs(\n",
    "        [\n",
    "            indk, \n",
    "            over, \n",
    "            info[info[\"type\"].isin([\"Turistkontor\"])]\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(service) > 0:\n",
    "        service.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/service.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    print(f\"\\t {kode} done\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* ~~0269: no over~~\n",
    "* ~~\"0183\", \"0147\", \"0190\", \"0825\": no fortidsminder~~\n",
    "* ~~\"0157\", \"0175\": no udflugt~~\n",
    "* ~~0161, 0163: no besoeg, no fortidsminder~~\n",
    "* ~~0240: no besoeg, no overnatning~~\n",
    "* ~~0411: no besoeg, no indkoeb, no vaerdifuld, no bad (!)~~\n",
    "* ~~0187: no fortidsminder, no overnatning~~\n",
    "* 0159: no over, no udflugt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "start kode 0269\n",
      "\t data read in\n",
      "\t 0269 done\n"
     ]
    }
   ],
   "source": [
    "# SPECIAL CASES: NO OVERNATNING\n",
    "\n",
    "for kode in [\"0269\"]:\n",
    "    \n",
    "    print(f\"start kode {kode}\")\n",
    "\n",
    "    anve = gpd.read_file(f\"data/raw/municipality_data/{kode}/anvendelse.gpkg\")\n",
    "    anve = anve[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    anve = anve.explode(index_parts=False)\n",
    "    anve = anve[anve.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    besk = gpd.read_file(f\"data/raw/municipality_data/{kode}/beskyttet_natur.gpkg\")\n",
    "    besk = besk[[\"type\", \"geometry\"]]\n",
    "    besk[\"gruppe\"] = \"beskyttet_natur\"\n",
    "    besk = besk.explode(index_parts=False)\n",
    "    besk = besk[besk.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    beso = gpd.read_file(f\"data/raw/municipality_data/{kode}/besoeg.gpkg\")\n",
    "    beso = beso[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    beso = beso.explode(index_parts=False)\n",
    "    beso = beso[beso.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fort = gpd.read_file(f\"data/raw/municipality_data/{kode}/fortidsminder.gpkg\")\n",
    "    fort = fort[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    fort = fort.explode(index_parts=False)\n",
    "    fort = fort[fort.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fred = gpd.read_file(f\"data/raw/municipality_data/{kode}/frednatpark.gpkg\")\n",
    "    fred = fred[[\"type\", \"geometry\"]]\n",
    "    fred[\"gruppe\"] = \"frednatpark\"\n",
    "    fred = fred.explode(index_parts=False)\n",
    "    fred = fred[fred.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    indk = gpd.read_file(f\"data/raw/municipality_data/{kode}/indkoeb.gpkg\")\n",
    "    indk = indk[[\"type\", \"geometry\"]]\n",
    "    indk[\"gruppe\"] = \"indkoeb\"\n",
    "    indk = indk.explode(index_parts=False)\n",
    "    indk = indk[indk.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    info = gpd.read_file(f\"data/raw/municipality_data/{kode}/infoservice.gpkg\")\n",
    "    info = info[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    info = info.explode(index_parts=False)\n",
    "    info = info[info.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    land = gpd.read_file(f\"data/raw/municipality_data/{kode}/landemaerker.gpkg\")\n",
    "    land = land[[\"type\", \"geometry\"]]\n",
    "    land[\"gruppe\"] = \"landemaerker\"\n",
    "    land = land.explode(index_parts=False)\n",
    "    land = land[land.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # over = gpd.read_file(f\"data/raw/municipality_data/{kode}/overnatning.gpkg\")\n",
    "    # over = over[[\"type\", \"geometry\"]]\n",
    "    # over[\"gruppe\"] = \"overnatning\"\n",
    "    # over = over.explode(index_parts=False)\n",
    "    # over = over[over.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    rast = gpd.read_file(f\"data/raw/municipality_data/{kode}/rasteplads.gpkg\")\n",
    "    rast = rast[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    rast = rast.explode(index_parts=False)\n",
    "    rast = rast[rast.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    udfl = gpd.read_file(f\"data/raw/municipality_data/{kode}/udflugt.gpkg\")\n",
    "    udfl = udfl[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    udfl = udfl.explode(index_parts=False)\n",
    "    udfl = udfl[udfl.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    vaer = gpd.read_file(f\"data/raw/municipality_data/{kode}/vaerdifuld.gpkg\")\n",
    "    vaer = vaer[[\"type\", \"geometry\"]]\n",
    "    vaer[\"gruppe\"] = \"vaerdifuld\"\n",
    "    vaer = vaer.explode(index_parts=False)\n",
    "    vaer = vaer[vaer.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    print(\"\\t data read in\")\n",
    "\n",
    "    nature = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Naturtyper\") & -anve[\"type\"].isin([\"Råstofområde\", \"Gartneri\"])].copy(),\n",
    "            anve[anve[\"gruppe\"]==\"Skovinddeling\"].copy(),\n",
    "            besk.copy(),\n",
    "            vaer[vaer[\"type\"]==\"Naturbeskyttelsesområde\"].copy(),\n",
    "            fred.copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(nature) > 0:\n",
    "        nature.to_file(f\"data/raw/municipality_data/{kode}/nature.gpkg\", index = False)\n",
    "    \n",
    "        \n",
    "    culture = merge_gdfs(\n",
    "        [\n",
    "            anve[anve[\"type\"].isin(['Bykerne'])].copy(),\n",
    "            vaer[vaer[\"type\"].isin(\n",
    "                [\n",
    "                    'Kulturhistorisk bevaringsværdi', \n",
    "                    'Værdifuldt kulturmiljø', \n",
    "                    'Bevaringsværdig landskab'\n",
    "                ]\n",
    "            )].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(culture) > 0:\n",
    "        culture.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/culture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    sommerhus = anve[\n",
    "        anve[\"type\"].isin(\n",
    "            ['Sommerhusområde', 'Sommerhusområde skov']\n",
    "            )\n",
    "        ].copy().reset_index(drop=True)\n",
    "    if len(sommerhus) > 0:\n",
    "        sommerhus.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/sommerhus.gpkg\", \n",
    "            index = False)\n",
    "\n",
    "    agri = anve[\n",
    "        anve[\"gruppe\"]==\"Dyrkede arealer\"].copy().reset_index(\n",
    "            drop=True)\n",
    "\n",
    "    if len(agri) > 0:\n",
    "        agri.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/agriculture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    verify = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Teknisk areal\") & (-anve[\"type\"].isin([\"Sportsanlæg\"]))].copy(),\n",
    "            anve[(anve[\"gruppe\"]==\"Bymæssig anvendelse\") & (anve[\"type\"].isin([\"Høj bebyggelse\", \"Erhverv\"]))].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(verify) > 0:\n",
    "        verify.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/verify.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "\n",
    "    pois = merge_gdfs(\n",
    "        [beso, fort, udfl, land]\n",
    "    )\n",
    "\n",
    "    if len(pois) > 0:\n",
    "        pois.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/pois.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    faci = merge_gdfs(\n",
    "        [\n",
    "            info[-info[\"type\"].isin([\"Turistkontor\"])].copy(),\n",
    "            rast,\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(faci) > 0:\n",
    "        faci.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/facilities.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    service = merge_gdfs(\n",
    "        [\n",
    "            indk, \n",
    "            #over, \n",
    "            info[info[\"type\"].isin([\"Turistkontor\"])]]\n",
    "    )\n",
    "\n",
    "    if len(service) > 0:\n",
    "        service.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/service.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    print(f\"\\t {kode} done\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "start kode 0183\n",
      "\t data read in\n",
      "\t 0183 done\n",
      "start kode 0147\n",
      "\t data read in\n",
      "\t 0147 done\n",
      "start kode 0190\n",
      "\t data read in\n",
      "\t 0190 done\n",
      "start kode 0825\n",
      "\t data read in\n",
      "\t 0825 done\n"
     ]
    }
   ],
   "source": [
    "# SPECIAL CASES: NO FORTIDSMINDER\n",
    "\n",
    "for kode in [\"0183\", \"0147\", \"0190\", \"0825\"]:\n",
    "    \n",
    "    print(f\"start kode {kode}\")\n",
    "\n",
    "    anve = gpd.read_file(f\"data/raw/municipality_data/{kode}/anvendelse.gpkg\")\n",
    "    anve = anve[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    anve = anve.explode(index_parts=False)\n",
    "    anve = anve[anve.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    besk = gpd.read_file(f\"data/raw/municipality_data/{kode}/beskyttet_natur.gpkg\")\n",
    "    besk = besk[[\"type\", \"geometry\"]]\n",
    "    besk[\"gruppe\"] = \"beskyttet_natur\"\n",
    "    besk = besk.explode(index_parts=False)\n",
    "    besk = besk[besk.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    beso = gpd.read_file(f\"data/raw/municipality_data/{kode}/besoeg.gpkg\")\n",
    "    beso = beso[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    beso = beso.explode(index_parts=False)\n",
    "    beso = beso[beso.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # fort = gpd.read_file(f\"data/raw/municipality_data/{kode}/fortidsminder.gpkg\")\n",
    "    # fort = fort[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    # fort = fort.explode(index_parts=False)\n",
    "    # fort = fort[fort.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fred = gpd.read_file(f\"data/raw/municipality_data/{kode}/frednatpark.gpkg\")\n",
    "    fred = fred[[\"type\", \"geometry\"]]\n",
    "    fred[\"gruppe\"] = \"frednatpark\"\n",
    "    fred = fred.explode(index_parts=False)\n",
    "    fred = fred[fred.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    indk = gpd.read_file(f\"data/raw/municipality_data/{kode}/indkoeb.gpkg\")\n",
    "    indk = indk[[\"type\", \"geometry\"]]\n",
    "    indk[\"gruppe\"] = \"indkoeb\"\n",
    "    indk = indk.explode(index_parts=False)\n",
    "    indk = indk[indk.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    info = gpd.read_file(f\"data/raw/municipality_data/{kode}/infoservice.gpkg\")\n",
    "    info = info[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    info = info.explode(index_parts=False)\n",
    "    info = info[info.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    land = gpd.read_file(f\"data/raw/municipality_data/{kode}/landemaerker.gpkg\")\n",
    "    land = land[[\"type\", \"geometry\"]]\n",
    "    land[\"gruppe\"] = \"landemaerker\"\n",
    "    land = land.explode(index_parts=False)\n",
    "    land = land[land.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    over = gpd.read_file(f\"data/raw/municipality_data/{kode}/overnatning.gpkg\")\n",
    "    over = over[[\"type\", \"geometry\"]]\n",
    "    over[\"gruppe\"] = \"overnatning\"\n",
    "    over = over.explode(index_parts=False)\n",
    "    over = over[over.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    rast = gpd.read_file(f\"data/raw/municipality_data/{kode}/rasteplads.gpkg\")\n",
    "    rast = rast[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    rast = rast.explode(index_parts=False)\n",
    "    rast = rast[rast.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    udfl = gpd.read_file(f\"data/raw/municipality_data/{kode}/udflugt.gpkg\")\n",
    "    udfl = udfl[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    udfl = udfl.explode(index_parts=False)\n",
    "    udfl = udfl[udfl.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    vaer = gpd.read_file(f\"data/raw/municipality_data/{kode}/vaerdifuld.gpkg\")\n",
    "    vaer = vaer[[\"type\", \"geometry\"]]\n",
    "    vaer[\"gruppe\"] = \"vaerdifuld\"\n",
    "    vaer = vaer.explode(index_parts=False)\n",
    "    vaer = vaer[vaer.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    print(\"\\t data read in\")\n",
    "\n",
    "    nature = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Naturtyper\") & -anve[\"type\"].isin([\"Råstofområde\", \"Gartneri\"])].copy(),\n",
    "            anve[anve[\"gruppe\"]==\"Skovinddeling\"].copy(),\n",
    "            besk.copy(),\n",
    "            vaer[vaer[\"type\"]==\"Naturbeskyttelsesområde\"].copy(),\n",
    "            fred.copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(nature) > 0:\n",
    "        nature.to_file(f\"data/raw/municipality_data/{kode}/nature.gpkg\", index = False)\n",
    "    \n",
    "        \n",
    "    culture = merge_gdfs(\n",
    "        [\n",
    "            anve[anve[\"type\"].isin(['Bykerne'])].copy(),\n",
    "            vaer[vaer[\"type\"].isin(\n",
    "                [\n",
    "                    'Kulturhistorisk bevaringsværdi', \n",
    "                    'Værdifuldt kulturmiljø', \n",
    "                    'Bevaringsværdig landskab'\n",
    "                ]\n",
    "            )].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(culture) > 0:\n",
    "        culture.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/culture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    sommerhus = anve[\n",
    "        anve[\"type\"].isin(\n",
    "            ['Sommerhusområde', 'Sommerhusområde skov']\n",
    "            )\n",
    "        ].copy().reset_index(drop=True)\n",
    "    if len(sommerhus) > 0:\n",
    "        sommerhus.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/sommerhus.gpkg\", \n",
    "            index = False)\n",
    "\n",
    "    agri = anve[\n",
    "        anve[\"gruppe\"]==\"Dyrkede arealer\"].copy().reset_index(\n",
    "            drop=True)\n",
    "\n",
    "    if len(agri) > 0:\n",
    "        agri.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/agriculture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    verify = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Teknisk areal\") & (-anve[\"type\"].isin([\"Sportsanlæg\"]))].copy(),\n",
    "            anve[(anve[\"gruppe\"]==\"Bymæssig anvendelse\") & (anve[\"type\"].isin([\"Høj bebyggelse\", \"Erhverv\"]))].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(verify) > 0:\n",
    "        verify.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/verify.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "\n",
    "    pois = merge_gdfs(\n",
    "        [\n",
    "            beso, \n",
    "            #fort, \n",
    "            udfl, \n",
    "            land]\n",
    "    )\n",
    "\n",
    "    if len(pois) > 0:\n",
    "        pois.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/pois.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    faci = merge_gdfs(\n",
    "        [\n",
    "            info[-info[\"type\"].isin([\"Turistkontor\"])].copy(),\n",
    "            rast\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(faci) > 0:\n",
    "        faci.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/facilities.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    service = merge_gdfs(\n",
    "        [indk, over, info[info[\"type\"].isin([\"Turistkontor\"])]]\n",
    "    )\n",
    "\n",
    "    if len(service) > 0:\n",
    "        service.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/service.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    print(f\"\\t {kode} done\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "start kode 0157\n",
      "\t data read in\n",
      "\t 0157 done\n",
      "start kode 0175\n",
      "\t data read in\n",
      "\t 0175 done\n"
     ]
    }
   ],
   "source": [
    "# SPECIAL CASES: NO UDFLUGT\n",
    "\n",
    "for kode in [\"0157\", \"0175\"]:\n",
    "    \n",
    "    print(f\"start kode {kode}\")\n",
    "\n",
    "    anve = gpd.read_file(f\"data/raw/municipality_data/{kode}/anvendelse.gpkg\")\n",
    "    anve = anve[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    anve = anve.explode(index_parts=False)\n",
    "    anve = anve[anve.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    besk = gpd.read_file(f\"data/raw/municipality_data/{kode}/beskyttet_natur.gpkg\")\n",
    "    besk = besk[[\"type\", \"geometry\"]]\n",
    "    besk[\"gruppe\"] = \"beskyttet_natur\"\n",
    "    besk = besk.explode(index_parts=False)\n",
    "    besk = besk[besk.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    beso = gpd.read_file(f\"data/raw/municipality_data/{kode}/besoeg.gpkg\")\n",
    "    beso = beso[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    beso = beso.explode(index_parts=False)\n",
    "    beso = beso[beso.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fort = gpd.read_file(f\"data/raw/municipality_data/{kode}/fortidsminder.gpkg\")\n",
    "    fort = fort[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    fort = fort.explode(index_parts=False)\n",
    "    fort = fort[fort.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fred = gpd.read_file(f\"data/raw/municipality_data/{kode}/frednatpark.gpkg\")\n",
    "    fred = fred[[\"type\", \"geometry\"]]\n",
    "    fred[\"gruppe\"] = \"frednatpark\"\n",
    "    fred = fred.explode(index_parts=False)\n",
    "    fred = fred[fred.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    indk = gpd.read_file(f\"data/raw/municipality_data/{kode}/indkoeb.gpkg\")\n",
    "    indk = indk[[\"type\", \"geometry\"]]\n",
    "    indk[\"gruppe\"] = \"indkoeb\"\n",
    "    indk = indk.explode(index_parts=False)\n",
    "    indk = indk[indk.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    info = gpd.read_file(f\"data/raw/municipality_data/{kode}/infoservice.gpkg\")\n",
    "    info = info[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    info = info.explode(index_parts=False)\n",
    "    info = info[info.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    land = gpd.read_file(f\"data/raw/municipality_data/{kode}/landemaerker.gpkg\")\n",
    "    land = land[[\"type\", \"geometry\"]]\n",
    "    land[\"gruppe\"] = \"landemaerker\"\n",
    "    land = land.explode(index_parts=False)\n",
    "    land = land[land.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    over = gpd.read_file(f\"data/raw/municipality_data/{kode}/overnatning.gpkg\")\n",
    "    over = over[[\"type\", \"geometry\"]]\n",
    "    over[\"gruppe\"] = \"overnatning\"\n",
    "    over = over.explode(index_parts=False)\n",
    "    over = over[over.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    rast = gpd.read_file(f\"data/raw/municipality_data/{kode}/rasteplads.gpkg\")\n",
    "    rast = rast[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    rast = rast.explode(index_parts=False)\n",
    "    rast = rast[rast.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # udfl = gpd.read_file(f\"data/raw/municipality_data/{kode}/udflugt.gpkg\")\n",
    "    # udfl = udfl[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    # udfl = udfl.explode(index_parts=False)\n",
    "    # udfl = udfl[udfl.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    vaer = gpd.read_file(f\"data/raw/municipality_data/{kode}/vaerdifuld.gpkg\")\n",
    "    vaer = vaer[[\"type\", \"geometry\"]]\n",
    "    vaer[\"gruppe\"] = \"vaerdifuld\"\n",
    "    vaer = vaer.explode(index_parts=False)\n",
    "    vaer = vaer[vaer.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    print(\"\\t data read in\")\n",
    "\n",
    "    nature = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Naturtyper\") & -anve[\"type\"].isin([\"Råstofområde\", \"Gartneri\"])].copy(),\n",
    "            anve[anve[\"gruppe\"]==\"Skovinddeling\"].copy(),\n",
    "            besk.copy(),\n",
    "            vaer[vaer[\"type\"]==\"Naturbeskyttelsesområde\"].copy(),\n",
    "            fred.copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(nature) > 0:\n",
    "        nature.to_file(f\"data/raw/municipality_data/{kode}/nature.gpkg\", index = False)\n",
    "    \n",
    "        \n",
    "    culture = merge_gdfs(\n",
    "        [\n",
    "            anve[anve[\"type\"].isin(['Bykerne'])].copy(),\n",
    "            vaer[vaer[\"type\"].isin(\n",
    "                [\n",
    "                    'Kulturhistorisk bevaringsværdi', \n",
    "                    'Værdifuldt kulturmiljø', \n",
    "                    'Bevaringsværdig landskab'\n",
    "                ]\n",
    "            )].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(culture) > 0:\n",
    "        culture.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/culture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    sommerhus = anve[\n",
    "        anve[\"type\"].isin(\n",
    "            ['Sommerhusområde', 'Sommerhusområde skov']\n",
    "            )\n",
    "        ].copy().reset_index(drop=True)\n",
    "    if len(sommerhus) > 0:\n",
    "        sommerhus.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/sommerhus.gpkg\", \n",
    "            index = False)\n",
    "\n",
    "    agri = anve[\n",
    "        anve[\"gruppe\"]==\"Dyrkede arealer\"].copy().reset_index(\n",
    "            drop=True)\n",
    "\n",
    "    if len(agri) > 0:\n",
    "        agri.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/agriculture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    verify = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Teknisk areal\") & (-anve[\"type\"].isin([\"Sportsanlæg\"]))].copy(),\n",
    "            anve[(anve[\"gruppe\"]==\"Bymæssig anvendelse\") & (anve[\"type\"].isin([\"Høj bebyggelse\", \"Erhverv\"]))].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(verify) > 0:\n",
    "        verify.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/verify.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "\n",
    "    pois = merge_gdfs(\n",
    "        [\n",
    "            beso, \n",
    "            fort, \n",
    "            # udfl, \n",
    "            land\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(pois) > 0:\n",
    "        pois.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/pois.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    faci = merge_gdfs(\n",
    "        [\n",
    "            info[-info[\"type\"].isin([\"Turistkontor\"])].copy(),\n",
    "            rast,\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(faci) > 0:\n",
    "        faci.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/facilities.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    service = merge_gdfs(\n",
    "        [indk, over, info[info[\"type\"].isin([\"Turistkontor\"])]]\n",
    "    )\n",
    "\n",
    "    if len(service) > 0:\n",
    "        service.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/service.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    print(f\"\\t {kode} done\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "start kode 0161\n",
      "\t data read in\n",
      "\t 0161 done\n",
      "start kode 0163\n",
      "\t data read in\n",
      "\t 0163 done\n"
     ]
    }
   ],
   "source": [
    "# SPECIAL CASES: NO BESOEG, NO FORTIDSMINDER\n",
    "\n",
    "for kode in [\"0161\", \"0163\"]:\n",
    "    \n",
    "    print(f\"start kode {kode}\")\n",
    "\n",
    "    anve = gpd.read_file(f\"data/raw/municipality_data/{kode}/anvendelse.gpkg\")\n",
    "    anve = anve[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    anve = anve.explode(index_parts=False)\n",
    "    anve = anve[anve.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    besk = gpd.read_file(f\"data/raw/municipality_data/{kode}/beskyttet_natur.gpkg\")\n",
    "    besk = besk[[\"type\", \"geometry\"]]\n",
    "    besk[\"gruppe\"] = \"beskyttet_natur\"\n",
    "    besk = besk.explode(index_parts=False)\n",
    "    besk = besk[besk.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # beso = gpd.read_file(f\"data/raw/municipality_data/{kode}/besoeg.gpkg\")\n",
    "    # beso = beso[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    # beso = beso.explode(index_parts=False)\n",
    "    # beso = beso[beso.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # fort = gpd.read_file(f\"data/raw/municipality_data/{kode}/fortidsminder.gpkg\")\n",
    "    # fort = fort[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    # fort = fort.explode(index_parts=False)\n",
    "    # fort = fort[fort.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fred = gpd.read_file(f\"data/raw/municipality_data/{kode}/frednatpark.gpkg\")\n",
    "    fred = fred[[\"type\", \"geometry\"]]\n",
    "    fred[\"gruppe\"] = \"frednatpark\"\n",
    "    fred = fred.explode(index_parts=False)\n",
    "    fred = fred[fred.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    indk = gpd.read_file(f\"data/raw/municipality_data/{kode}/indkoeb.gpkg\")\n",
    "    indk = indk[[\"type\", \"geometry\"]]\n",
    "    indk[\"gruppe\"] = \"indkoeb\"\n",
    "    indk = indk.explode(index_parts=False)\n",
    "    indk = indk[indk.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    info = gpd.read_file(f\"data/raw/municipality_data/{kode}/infoservice.gpkg\")\n",
    "    info = info[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    info = info.explode(index_parts=False)\n",
    "    info = info[info.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    land = gpd.read_file(f\"data/raw/municipality_data/{kode}/landemaerker.gpkg\")\n",
    "    land = land[[\"type\", \"geometry\"]]\n",
    "    land[\"gruppe\"] = \"landemaerker\"\n",
    "    land = land.explode(index_parts=False)\n",
    "    land = land[land.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    over = gpd.read_file(f\"data/raw/municipality_data/{kode}/overnatning.gpkg\")\n",
    "    over = over[[\"type\", \"geometry\"]]\n",
    "    over[\"gruppe\"] = \"overnatning\"\n",
    "    over = over.explode(index_parts=False)\n",
    "    over = over[over.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    rast = gpd.read_file(f\"data/raw/municipality_data/{kode}/rasteplads.gpkg\")\n",
    "    rast = rast[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    rast = rast.explode(index_parts=False)\n",
    "    rast = rast[rast.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    udfl = gpd.read_file(f\"data/raw/municipality_data/{kode}/udflugt.gpkg\")\n",
    "    udfl = udfl[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    udfl = udfl.explode(index_parts=False)\n",
    "    udfl = udfl[udfl.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    vaer = gpd.read_file(f\"data/raw/municipality_data/{kode}/vaerdifuld.gpkg\")\n",
    "    vaer = vaer[[\"type\", \"geometry\"]]\n",
    "    vaer[\"gruppe\"] = \"vaerdifuld\"\n",
    "    vaer = vaer.explode(index_parts=False)\n",
    "    vaer = vaer[vaer.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    print(\"\\t data read in\")\n",
    "\n",
    "    nature = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Naturtyper\") & -anve[\"type\"].isin([\"Råstofområde\", \"Gartneri\"])].copy(),\n",
    "            anve[anve[\"gruppe\"]==\"Skovinddeling\"].copy(),\n",
    "            besk.copy(),\n",
    "            vaer[vaer[\"type\"]==\"Naturbeskyttelsesområde\"].copy(),\n",
    "            fred.copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(nature) > 0:\n",
    "        nature.to_file(f\"data/raw/municipality_data/{kode}/nature.gpkg\", index = False)\n",
    "    \n",
    "        \n",
    "    culture = merge_gdfs(\n",
    "        [\n",
    "            anve[anve[\"type\"].isin(['Bykerne'])].copy(),\n",
    "            vaer[vaer[\"type\"].isin(\n",
    "                [\n",
    "                    'Kulturhistorisk bevaringsværdi', \n",
    "                    'Værdifuldt kulturmiljø', \n",
    "                    'Bevaringsværdig landskab'\n",
    "                ]\n",
    "            )].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(culture) > 0:\n",
    "        culture.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/culture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    sommerhus = anve[\n",
    "        anve[\"type\"].isin(\n",
    "            ['Sommerhusområde', 'Sommerhusområde skov']\n",
    "            )\n",
    "        ].copy().reset_index(drop=True)\n",
    "    if len(sommerhus) > 0:\n",
    "        sommerhus.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/sommerhus.gpkg\", \n",
    "            index = False)\n",
    "\n",
    "    agri = anve[\n",
    "        anve[\"gruppe\"]==\"Dyrkede arealer\"].copy().reset_index(\n",
    "            drop=True)\n",
    "\n",
    "    if len(agri) > 0:\n",
    "        agri.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/agriculture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    verify = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Teknisk areal\") & (-anve[\"type\"].isin([\"Sportsanlæg\"]))].copy(),\n",
    "            anve[(anve[\"gruppe\"]==\"Bymæssig anvendelse\") & (anve[\"type\"].isin([\"Høj bebyggelse\", \"Erhverv\"]))].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(verify) > 0:\n",
    "        verify.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/verify.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "\n",
    "    pois = merge_gdfs(\n",
    "        [\n",
    "            # beso, \n",
    "            # fort, \n",
    "            udfl, \n",
    "            land\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(pois) > 0:\n",
    "        pois.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/pois.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    faci = merge_gdfs(\n",
    "        [\n",
    "            info[-info[\"type\"].isin([\"Turistkontor\"])].copy(),\n",
    "            rast,\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(faci) > 0:\n",
    "        faci.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/facilities.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    service = merge_gdfs(\n",
    "        [indk, over, info[info[\"type\"].isin([\"Turistkontor\"])]]\n",
    "    )\n",
    "\n",
    "    if len(service) > 0:\n",
    "        service.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/service.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    print(f\"\\t {kode} done\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "start kode 0240\n",
      "\t data read in\n",
      "\t 0240 done\n"
     ]
    }
   ],
   "source": [
    "# SPECIAL CASES: 0240: no besoeg, no overnatning\n",
    "\n",
    "for kode in [\"0240\"]:\n",
    "    \n",
    "    print(f\"start kode {kode}\")\n",
    "\n",
    "    anve = gpd.read_file(f\"data/raw/municipality_data/{kode}/anvendelse.gpkg\")\n",
    "    anve = anve[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    anve = anve.explode(index_parts=False)\n",
    "    anve = anve[anve.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    besk = gpd.read_file(f\"data/raw/municipality_data/{kode}/beskyttet_natur.gpkg\")\n",
    "    besk = besk[[\"type\", \"geometry\"]]\n",
    "    besk[\"gruppe\"] = \"beskyttet_natur\"\n",
    "    besk = besk.explode(index_parts=False)\n",
    "    besk = besk[besk.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # beso = gpd.read_file(f\"data/raw/municipality_data/{kode}/besoeg.gpkg\")\n",
    "    # beso = beso[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    # beso = beso.explode(index_parts=False)\n",
    "    # beso = beso[beso.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fort = gpd.read_file(f\"data/raw/municipality_data/{kode}/fortidsminder.gpkg\")\n",
    "    fort = fort[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    fort = fort.explode(index_parts=False)\n",
    "    fort = fort[fort.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fred = gpd.read_file(f\"data/raw/municipality_data/{kode}/frednatpark.gpkg\")\n",
    "    fred = fred[[\"type\", \"geometry\"]]\n",
    "    fred[\"gruppe\"] = \"frednatpark\"\n",
    "    fred = fred.explode(index_parts=False)\n",
    "    fred = fred[fred.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    indk = gpd.read_file(f\"data/raw/municipality_data/{kode}/indkoeb.gpkg\")\n",
    "    indk = indk[[\"type\", \"geometry\"]]\n",
    "    indk[\"gruppe\"] = \"indkoeb\"\n",
    "    indk = indk.explode(index_parts=False)\n",
    "    indk = indk[indk.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    info = gpd.read_file(f\"data/raw/municipality_data/{kode}/infoservice.gpkg\")\n",
    "    info = info[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    info = info.explode(index_parts=False)\n",
    "    info = info[info.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    land = gpd.read_file(f\"data/raw/municipality_data/{kode}/landemaerker.gpkg\")\n",
    "    land = land[[\"type\", \"geometry\"]]\n",
    "    land[\"gruppe\"] = \"landemaerker\"\n",
    "    land = land.explode(index_parts=False)\n",
    "    land = land[land.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # over = gpd.read_file(f\"data/raw/municipality_data/{kode}/overnatning.gpkg\")\n",
    "    # over = over[[\"type\", \"geometry\"]]\n",
    "    # over[\"gruppe\"] = \"overnatning\"\n",
    "    # over = over.explode(index_parts=False)\n",
    "    # over = over[over.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    rast = gpd.read_file(f\"data/raw/municipality_data/{kode}/rasteplads.gpkg\")\n",
    "    rast = rast[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    rast = rast.explode(index_parts=False)\n",
    "    rast = rast[rast.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    udfl = gpd.read_file(f\"data/raw/municipality_data/{kode}/udflugt.gpkg\")\n",
    "    udfl = udfl[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    udfl = udfl.explode(index_parts=False)\n",
    "    udfl = udfl[udfl.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    vaer = gpd.read_file(f\"data/raw/municipality_data/{kode}/vaerdifuld.gpkg\")\n",
    "    vaer = vaer[[\"type\", \"geometry\"]]\n",
    "    vaer[\"gruppe\"] = \"vaerdifuld\"\n",
    "    vaer = vaer.explode(index_parts=False)\n",
    "    vaer = vaer[vaer.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    print(\"\\t data read in\")\n",
    "\n",
    "    nature = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Naturtyper\") & -anve[\"type\"].isin([\"Råstofområde\", \"Gartneri\"])].copy(),\n",
    "            anve[anve[\"gruppe\"]==\"Skovinddeling\"].copy(),\n",
    "            besk.copy(),\n",
    "            vaer[vaer[\"type\"]==\"Naturbeskyttelsesområde\"].copy(),\n",
    "            fred.copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(nature) > 0:\n",
    "        nature.to_file(f\"data/raw/municipality_data/{kode}/nature.gpkg\", index = False)\n",
    "    \n",
    "        \n",
    "    culture = merge_gdfs(\n",
    "        [\n",
    "            anve[anve[\"type\"].isin(['Bykerne'])].copy(),\n",
    "            vaer[vaer[\"type\"].isin(\n",
    "                [\n",
    "                    'Kulturhistorisk bevaringsværdi', \n",
    "                    'Værdifuldt kulturmiljø', \n",
    "                    'Bevaringsværdig landskab'\n",
    "                ]\n",
    "            )].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(culture) > 0:\n",
    "        culture.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/culture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    sommerhus = anve[\n",
    "        anve[\"type\"].isin(\n",
    "            ['Sommerhusområde', 'Sommerhusområde skov']\n",
    "            )\n",
    "        ].copy().reset_index(drop=True)\n",
    "    if len(sommerhus) > 0:\n",
    "        sommerhus.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/sommerhus.gpkg\", \n",
    "            index = False)\n",
    "\n",
    "    agri = anve[\n",
    "        anve[\"gruppe\"]==\"Dyrkede arealer\"].copy().reset_index(\n",
    "            drop=True)\n",
    "\n",
    "    if len(agri) > 0:\n",
    "        agri.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/agriculture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    verify = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Teknisk areal\") & (-anve[\"type\"].isin([\"Sportsanlæg\"]))].copy(),\n",
    "            anve[(anve[\"gruppe\"]==\"Bymæssig anvendelse\") & (anve[\"type\"].isin([\"Høj bebyggelse\", \"Erhverv\"]))].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(verify) > 0:\n",
    "        verify.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/verify.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "\n",
    "    pois = merge_gdfs(\n",
    "        [\n",
    "            # beso, \n",
    "            fort, \n",
    "            udfl, \n",
    "            land\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(pois) > 0:\n",
    "        pois.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/pois.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    faci = merge_gdfs(\n",
    "        [\n",
    "            info[-info[\"type\"].isin([\"Turistkontor\"])].copy(),\n",
    "            rast,\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(faci) > 0:\n",
    "        faci.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/facilities.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    service = merge_gdfs(\n",
    "        [\n",
    "            indk, \n",
    "            # over, \n",
    "            info[info[\"type\"].isin([\"Turistkontor\"])]]\n",
    "    )\n",
    "\n",
    "    if len(service) > 0:\n",
    "        service.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/service.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    print(f\"\\t {kode} done\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "start kode 0411\n",
      "\t data read in\n",
      "\t 0411 done\n"
     ]
    }
   ],
   "source": [
    "# SPECIAL CASES: 0411: no besoeg, no indkoeb, no vaerdifuld\n",
    "\n",
    "for kode in [\"0411\"]:\n",
    "    \n",
    "    print(f\"start kode {kode}\")\n",
    "\n",
    "    anve = gpd.read_file(f\"data/raw/municipality_data/{kode}/anvendelse.gpkg\")\n",
    "    anve = anve[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    anve = anve.explode(index_parts=False)\n",
    "    anve = anve[anve.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    besk = gpd.read_file(f\"data/raw/municipality_data/{kode}/beskyttet_natur.gpkg\")\n",
    "    besk = besk[[\"type\", \"geometry\"]]\n",
    "    besk[\"gruppe\"] = \"beskyttet_natur\"\n",
    "    besk = besk.explode(index_parts=False)\n",
    "    besk = besk[besk.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # beso = gpd.read_file(f\"data/raw/municipality_data/{kode}/besoeg.gpkg\")\n",
    "    # beso = beso[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    # beso = beso.explode(index_parts=False)\n",
    "    # beso = beso[beso.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fort = gpd.read_file(f\"data/raw/municipality_data/{kode}/fortidsminder.gpkg\")\n",
    "    fort = fort[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    fort = fort.explode(index_parts=False)\n",
    "    fort = fort[fort.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fred = gpd.read_file(f\"data/raw/municipality_data/{kode}/frednatpark.gpkg\")\n",
    "    fred = fred[[\"type\", \"geometry\"]]\n",
    "    fred[\"gruppe\"] = \"frednatpark\"\n",
    "    fred = fred.explode(index_parts=False)\n",
    "    fred = fred[fred.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # indk = gpd.read_file(f\"data/raw/municipality_data/{kode}/indkoeb.gpkg\")\n",
    "    # indk = indk[[\"type\", \"geometry\"]]\n",
    "    # indk[\"gruppe\"] = \"indkoeb\"\n",
    "    # indk = indk.explode(index_parts=False)\n",
    "    # indk = indk[indk.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    info = gpd.read_file(f\"data/raw/municipality_data/{kode}/infoservice.gpkg\")\n",
    "    info = info[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    info = info.explode(index_parts=False)\n",
    "    info = info[info.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    land = gpd.read_file(f\"data/raw/municipality_data/{kode}/landemaerker.gpkg\")\n",
    "    land = land[[\"type\", \"geometry\"]]\n",
    "    land[\"gruppe\"] = \"landemaerker\"\n",
    "    land = land.explode(index_parts=False)\n",
    "    land = land[land.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    over = gpd.read_file(f\"data/raw/municipality_data/{kode}/overnatning.gpkg\")\n",
    "    over = over[[\"type\", \"geometry\"]]\n",
    "    over[\"gruppe\"] = \"overnatning\"\n",
    "    over = over.explode(index_parts=False)\n",
    "    over = over[over.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    rast = gpd.read_file(f\"data/raw/municipality_data/{kode}/rasteplads.gpkg\")\n",
    "    rast = rast[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    rast = rast.explode(index_parts=False)\n",
    "    rast = rast[rast.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    udfl = gpd.read_file(f\"data/raw/municipality_data/{kode}/udflugt.gpkg\")\n",
    "    udfl = udfl[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    udfl = udfl.explode(index_parts=False)\n",
    "    udfl = udfl[udfl.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # vaer = gpd.read_file(f\"data/raw/municipality_data/{kode}/vaerdifuld.gpkg\")\n",
    "    # vaer = vaer[[\"type\", \"geometry\"]]\n",
    "    # vaer[\"gruppe\"] = \"vaerdifuld\"\n",
    "    # vaer = vaer.explode(index_parts=False)\n",
    "    # vaer = vaer[vaer.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    print(\"\\t data read in\")\n",
    "\n",
    "    nature = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Naturtyper\") & -anve[\"type\"].isin([\"Råstofområde\", \"Gartneri\"])].copy(),\n",
    "            anve[anve[\"gruppe\"]==\"Skovinddeling\"].copy(),\n",
    "            besk.copy(),\n",
    "            # vaer[vaer[\"type\"]==\"Naturbeskyttelsesområde\"].copy(),\n",
    "            fred.copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(nature) > 0:\n",
    "        nature.to_file(f\"data/raw/municipality_data/{kode}/nature.gpkg\", index = False)\n",
    "    \n",
    "        \n",
    "    culture = anve[anve[\"type\"].isin(['Bykerne'])].copy().reset_index(drop=True)\n",
    "\n",
    "    if len(culture) > 0:\n",
    "        culture.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/culture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    sommerhus = anve[\n",
    "        anve[\"type\"].isin(\n",
    "            ['Sommerhusområde', 'Sommerhusområde skov']\n",
    "            )\n",
    "        ].copy().reset_index(drop=True)\n",
    "    if len(sommerhus) > 0:\n",
    "        sommerhus.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/sommerhus.gpkg\", \n",
    "            index = False)\n",
    "\n",
    "    agri = anve[\n",
    "        anve[\"gruppe\"]==\"Dyrkede arealer\"].copy().reset_index(\n",
    "            drop=True)\n",
    "\n",
    "    if len(agri) > 0:\n",
    "        agri.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/agriculture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    # verify = merge_gdfs(\n",
    "    #     [\n",
    "    #         anve[(anve[\"gruppe\"]==\"Teknisk areal\") & (-anve[\"type\"].isin([\"Sportsanlæg\"]))].copy(),\n",
    "    #         anve[(anve[\"gruppe\"]==\"Bymæssig anvendelse\") & (anve[\"type\"].isin([\"Høj bebyggelse\", \"Erhverv\"]))].copy()\n",
    "    #     ]\n",
    "    # )\n",
    "\n",
    "    # if len(verify) > 0:\n",
    "    #     verify.to_file(\n",
    "    #         f\"data/raw/municipality_data/{kode}/verify.gpkg\", \n",
    "    #         index = False)\n",
    "        \n",
    "\n",
    "    pois = merge_gdfs(\n",
    "        [\n",
    "            # beso, \n",
    "            fort, \n",
    "            udfl, \n",
    "            land\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(pois) > 0:\n",
    "        pois.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/pois.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    faci = merge_gdfs(\n",
    "        [\n",
    "            info[-info[\"type\"].isin([\"Turistkontor\"])].copy(),\n",
    "            rast,\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(faci) > 0:\n",
    "        faci.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/facilities.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    service = merge_gdfs(\n",
    "        [indk, over, info[info[\"type\"].isin([\"Turistkontor\"])]]\n",
    "    )\n",
    "\n",
    "    if len(service) > 0:\n",
    "        service.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/service.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    print(f\"\\t {kode} done\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "start kode 0187\n",
      "\t data read in\n",
      "\t 0187 done\n"
     ]
    }
   ],
   "source": [
    "# SPECIAL CASES: 0187: no fortidsminder, no overnatning\n",
    "\n",
    "for kode in [\"0187\"]:\n",
    "    \n",
    "    print(f\"start kode {kode}\")\n",
    "\n",
    "    anve = gpd.read_file(f\"data/raw/municipality_data/{kode}/anvendelse.gpkg\")\n",
    "    anve = anve[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    anve = anve.explode(index_parts=False)\n",
    "    anve = anve[anve.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    besk = gpd.read_file(f\"data/raw/municipality_data/{kode}/beskyttet_natur.gpkg\")\n",
    "    besk = besk[[\"type\", \"geometry\"]]\n",
    "    besk[\"gruppe\"] = \"beskyttet_natur\"\n",
    "    besk = besk.explode(index_parts=False)\n",
    "    besk = besk[besk.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    beso = gpd.read_file(f\"data/raw/municipality_data/{kode}/besoeg.gpkg\")\n",
    "    beso = beso[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    beso = beso.explode(index_parts=False)\n",
    "    beso = beso[beso.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # fort = gpd.read_file(f\"data/raw/municipality_data/{kode}/fortidsminder.gpkg\")\n",
    "    # fort = fort[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    # fort = fort.explode(index_parts=False)\n",
    "    # fort = fort[fort.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fred = gpd.read_file(f\"data/raw/municipality_data/{kode}/frednatpark.gpkg\")\n",
    "    fred = fred[[\"type\", \"geometry\"]]\n",
    "    fred[\"gruppe\"] = \"frednatpark\"\n",
    "    fred = fred.explode(index_parts=False)\n",
    "    fred = fred[fred.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    indk = gpd.read_file(f\"data/raw/municipality_data/{kode}/indkoeb.gpkg\")\n",
    "    indk = indk[[\"type\", \"geometry\"]]\n",
    "    indk[\"gruppe\"] = \"indkoeb\"\n",
    "    indk = indk.explode(index_parts=False)\n",
    "    indk = indk[indk.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    info = gpd.read_file(f\"data/raw/municipality_data/{kode}/infoservice.gpkg\")\n",
    "    info = info[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    info = info.explode(index_parts=False)\n",
    "    info = info[info.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    land = gpd.read_file(f\"data/raw/municipality_data/{kode}/landemaerker.gpkg\")\n",
    "    land = land[[\"type\", \"geometry\"]]\n",
    "    land[\"gruppe\"] = \"landemaerker\"\n",
    "    land = land.explode(index_parts=False)\n",
    "    land = land[land.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # over = gpd.read_file(f\"data/raw/municipality_data/{kode}/overnatning.gpkg\")\n",
    "    # over = over[[\"type\", \"geometry\"]]\n",
    "    # over[\"gruppe\"] = \"overnatning\"\n",
    "    # over = over.explode(index_parts=False)\n",
    "    # over = over[over.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    rast = gpd.read_file(f\"data/raw/municipality_data/{kode}/rasteplads.gpkg\")\n",
    "    rast = rast[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    rast = rast.explode(index_parts=False)\n",
    "    rast = rast[rast.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    udfl = gpd.read_file(f\"data/raw/municipality_data/{kode}/udflugt.gpkg\")\n",
    "    udfl = udfl[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    udfl = udfl.explode(index_parts=False)\n",
    "    udfl = udfl[udfl.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    vaer = gpd.read_file(f\"data/raw/municipality_data/{kode}/vaerdifuld.gpkg\")\n",
    "    vaer = vaer[[\"type\", \"geometry\"]]\n",
    "    vaer[\"gruppe\"] = \"vaerdifuld\"\n",
    "    vaer = vaer.explode(index_parts=False)\n",
    "    vaer = vaer[vaer.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    print(\"\\t data read in\")\n",
    "\n",
    "    nature = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Naturtyper\") & -anve[\"type\"].isin([\"Råstofområde\", \"Gartneri\"])].copy(),\n",
    "            anve[anve[\"gruppe\"]==\"Skovinddeling\"].copy(),\n",
    "            besk.copy(),\n",
    "            vaer[vaer[\"type\"]==\"Naturbeskyttelsesområde\"].copy(),\n",
    "            fred.copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(nature) > 0:\n",
    "        nature.to_file(f\"data/raw/municipality_data/{kode}/nature.gpkg\", index = False)\n",
    "    \n",
    "        \n",
    "    culture = merge_gdfs(\n",
    "        [\n",
    "            anve[anve[\"type\"].isin(['Bykerne'])].copy(),\n",
    "            vaer[vaer[\"type\"].isin(\n",
    "                [\n",
    "                    'Kulturhistorisk bevaringsværdi', \n",
    "                    'Værdifuldt kulturmiljø', \n",
    "                    'Bevaringsværdig landskab'\n",
    "                ]\n",
    "            )].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(culture) > 0:\n",
    "        culture.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/culture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    sommerhus = anve[\n",
    "        anve[\"type\"].isin(\n",
    "            ['Sommerhusområde', 'Sommerhusområde skov']\n",
    "            )\n",
    "        ].copy().reset_index(drop=True)\n",
    "    if len(sommerhus) > 0:\n",
    "        sommerhus.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/sommerhus.gpkg\", \n",
    "            index = False)\n",
    "\n",
    "    agri = anve[\n",
    "        anve[\"gruppe\"]==\"Dyrkede arealer\"].copy().reset_index(\n",
    "            drop=True)\n",
    "\n",
    "    if len(agri) > 0:\n",
    "        agri.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/agriculture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    verify = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Teknisk areal\") & (-anve[\"type\"].isin([\"Sportsanlæg\"]))].copy(),\n",
    "            anve[(anve[\"gruppe\"]==\"Bymæssig anvendelse\") & (anve[\"type\"].isin([\"Høj bebyggelse\", \"Erhverv\"]))].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(verify) > 0:\n",
    "        verify.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/verify.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "\n",
    "    pois = merge_gdfs(\n",
    "        [\n",
    "            beso, \n",
    "            # fort, \n",
    "            udfl, \n",
    "            land]\n",
    "    )\n",
    "\n",
    "    if len(pois) > 0:\n",
    "        pois.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/pois.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    faci = merge_gdfs(\n",
    "        [\n",
    "            info[-info[\"type\"].isin([\"Turistkontor\"])].copy(),\n",
    "            rast,\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(faci) > 0:\n",
    "        faci.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/facilities.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    service = merge_gdfs(\n",
    "        [\n",
    "            indk, \n",
    "            # over, \n",
    "            info[info[\"type\"].isin([\"Turistkontor\"])]]\n",
    "    )\n",
    "\n",
    "    if len(service) > 0:\n",
    "        service.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/service.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    print(f\"\\t {kode} done\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "start kode 0159\n",
      "\t data read in\n",
      "\t 0159 done\n"
     ]
    }
   ],
   "source": [
    "# SPECIAL CASES: 0159: no over, no udflugt\n",
    "\n",
    "for kode in [\"0159\"]:\n",
    "    \n",
    "    print(f\"start kode {kode}\")\n",
    "\n",
    "    anve = gpd.read_file(f\"data/raw/municipality_data/{kode}/anvendelse.gpkg\")\n",
    "    anve = anve[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    anve = anve.explode(index_parts=False)\n",
    "    anve = anve[anve.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    besk = gpd.read_file(f\"data/raw/municipality_data/{kode}/beskyttet_natur.gpkg\")\n",
    "    besk = besk[[\"type\", \"geometry\"]]\n",
    "    besk[\"gruppe\"] = \"beskyttet_natur\"\n",
    "    besk = besk.explode(index_parts=False)\n",
    "    besk = besk[besk.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    beso = gpd.read_file(f\"data/raw/municipality_data/{kode}/besoeg.gpkg\")\n",
    "    beso = beso[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    beso = beso.explode(index_parts=False)\n",
    "    beso = beso[beso.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fort = gpd.read_file(f\"data/raw/municipality_data/{kode}/fortidsminder.gpkg\")\n",
    "    fort = fort[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    fort = fort.explode(index_parts=False)\n",
    "    fort = fort[fort.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    fred = gpd.read_file(f\"data/raw/municipality_data/{kode}/frednatpark.gpkg\")\n",
    "    fred = fred[[\"type\", \"geometry\"]]\n",
    "    fred[\"gruppe\"] = \"frednatpark\"\n",
    "    fred = fred.explode(index_parts=False)\n",
    "    fred = fred[fred.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    indk = gpd.read_file(f\"data/raw/municipality_data/{kode}/indkoeb.gpkg\")\n",
    "    indk = indk[[\"type\", \"geometry\"]]\n",
    "    indk[\"gruppe\"] = \"indkoeb\"\n",
    "    indk = indk.explode(index_parts=False)\n",
    "    indk = indk[indk.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    info = gpd.read_file(f\"data/raw/municipality_data/{kode}/infoservice.gpkg\")\n",
    "    info = info[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    info = info.explode(index_parts=False)\n",
    "    info = info[info.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    land = gpd.read_file(f\"data/raw/municipality_data/{kode}/landemaerker.gpkg\")\n",
    "    land = land[[\"type\", \"geometry\"]]\n",
    "    land[\"gruppe\"] = \"landemaerker\"\n",
    "    land = land.explode(index_parts=False)\n",
    "    land = land[land.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # over = gpd.read_file(f\"data/raw/municipality_data/{kode}/overnatning.gpkg\")\n",
    "    # over = over[[\"type\", \"geometry\"]]\n",
    "    # over[\"gruppe\"] = \"overnatning\"\n",
    "    # over = over.explode(index_parts=False)\n",
    "    # over = over[over.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    rast = gpd.read_file(f\"data/raw/municipality_data/{kode}/rasteplads.gpkg\")\n",
    "    rast = rast[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    rast = rast.explode(index_parts=False)\n",
    "    rast = rast[rast.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    # udfl = gpd.read_file(f\"data/raw/municipality_data/{kode}/udflugt.gpkg\")\n",
    "    # udfl = udfl[[\"gruppe\", \"type\", \"geometry\"]]\n",
    "    # udfl = udfl.explode(index_parts=False)\n",
    "    # udfl = udfl[udfl.geometry.type == \"Point\"].copy().reset_index(drop=True)\n",
    "\n",
    "    vaer = gpd.read_file(f\"data/raw/municipality_data/{kode}/vaerdifuld.gpkg\")\n",
    "    vaer = vaer[[\"type\", \"geometry\"]]\n",
    "    vaer[\"gruppe\"] = \"vaerdifuld\"\n",
    "    vaer = vaer.explode(index_parts=False)\n",
    "    vaer = vaer[vaer.geometry.type == \"Polygon\"].copy().reset_index(drop=True)\n",
    "\n",
    "    print(\"\\t data read in\")\n",
    "\n",
    "    nature = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Naturtyper\") & -anve[\"type\"].isin([\"Råstofområde\", \"Gartneri\"])].copy(),\n",
    "            anve[anve[\"gruppe\"]==\"Skovinddeling\"].copy(),\n",
    "            besk.copy(),\n",
    "            vaer[vaer[\"type\"]==\"Naturbeskyttelsesområde\"].copy(),\n",
    "            fred.copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(nature) > 0:\n",
    "        nature.to_file(f\"data/raw/municipality_data/{kode}/nature.gpkg\", index = False)\n",
    "    \n",
    "        \n",
    "    culture = merge_gdfs(\n",
    "        [\n",
    "            anve[anve[\"type\"].isin(['Bykerne'])].copy(),\n",
    "            vaer[vaer[\"type\"].isin(\n",
    "                [\n",
    "                    'Kulturhistorisk bevaringsværdi', \n",
    "                    'Værdifuldt kulturmiljø', \n",
    "                    'Bevaringsværdig landskab'\n",
    "                ]\n",
    "            )].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(culture) > 0:\n",
    "        culture.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/culture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    sommerhus = anve[\n",
    "        anve[\"type\"].isin(\n",
    "            ['Sommerhusområde', 'Sommerhusområde skov']\n",
    "            )\n",
    "        ].copy().reset_index(drop=True)\n",
    "    if len(sommerhus) > 0:\n",
    "        sommerhus.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/sommerhus.gpkg\", \n",
    "            index = False)\n",
    "\n",
    "    agri = anve[\n",
    "        anve[\"gruppe\"]==\"Dyrkede arealer\"].copy().reset_index(\n",
    "            drop=True)\n",
    "\n",
    "    if len(agri) > 0:\n",
    "        agri.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/agriculture.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    verify = merge_gdfs(\n",
    "        [\n",
    "            anve[(anve[\"gruppe\"]==\"Teknisk areal\") & (-anve[\"type\"].isin([\"Sportsanlæg\"]))].copy(),\n",
    "            anve[(anve[\"gruppe\"]==\"Bymæssig anvendelse\") & (anve[\"type\"].isin([\"Høj bebyggelse\", \"Erhverv\"]))].copy()\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(verify) > 0:\n",
    "        verify.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/verify.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "\n",
    "    pois = merge_gdfs(\n",
    "        [\n",
    "            beso, \n",
    "            fort, \n",
    "            # udfl, \n",
    "            land\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(pois) > 0:\n",
    "        pois.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/pois.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    faci = merge_gdfs(\n",
    "        [\n",
    "            info[-info[\"type\"].isin([\"Turistkontor\"])].copy(),\n",
    "            rast,\n",
    "        ]\n",
    "    )\n",
    "\n",
    "    if len(faci) > 0:\n",
    "        faci.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/facilities.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    service = merge_gdfs(\n",
    "        [\n",
    "            indk, \n",
    "            #over, \n",
    "            info[info[\"type\"].isin([\"Turistkontor\"])]]\n",
    "    )\n",
    "\n",
    "    if len(service) > 0:\n",
    "        service.to_file(\n",
    "            f\"data/raw/municipality_data/{kode}/service.gpkg\", \n",
    "            index = False)\n",
    "        \n",
    "    print(f\"\\t {kode} done\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "knupu",
   "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.11.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}