notebooks/2_processen_bressen.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "from pathlib import Path\n",
    "\n",
    "from bressen import BasisGegevens, Bressen, read_files"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Paden naar gegevens\n",
    "Gegevens zijn aangeleverd in verschillende directories:\n",
    "- 01_AGV_Dijktrajecten_Legger: dijktrajecten van waterschap Amstel Gooi en Vecht\n",
    "- 02_ARK_Kering: kering van het Amsterdam-Rijnkanaal\n",
    "- 03_Doorbraaklocaties: de locaties van de doorbraken\n",
    "- 04_Watervlakken: alle beschikbare watervlakken\n",
    "\n",
    "BASISGEGEVENS_GPKG en BRESSEN_GPKG verwijzen naar de basisgegevens (keringen, peil- en watervlakken) en het bressen-bestand"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "DATA_DIR = Path(r\"d:\\projecten\\D2405.waternet.bressen\\01.gegevens\")\n",
    "BASISGEGEVENS_GPKG = DATA_DIR / \"basisgegevens.gpkg\"\n",
    "BRESSEN_GPKG = DATA_DIR / \"bressen.gpkg\""
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Inlezen\n",
    "Inlezen van basisgegevens en bressen"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Inlezen basisgegevens vanuit een GeoPackage\n",
    "basisgegevens = BasisGegevens.from_gpkg(BASISGEGEVENS_GPKG)\n",
    "\n",
    "# Inlezen bressen vanuit meerdere bestanden naar een GeoDataFrame\n",
    "bressen_gdf = read_files(\n",
    "    [\n",
    "        DATA_DIR.joinpath(\"03_Doorbraaklocaties\", \"240122_breslocaties_gelijk_groter_1000.shp\"),\n",
    "        DATA_DIR.joinpath(\"03_Doorbraaklocaties\", \"240212_breslocaties_kleiner_1000_angle.shp\"),\n",
    "        DATA_DIR.joinpath(\"03_Doorbraaklocaties\", \"240212_breslocaties_Waterkeringen_Blaeu_ARK.shp\"),\n",
    "    ]\n",
    ")\n",
    "\n",
    "# Vullen van de bressenklasse met basisgegevens en de bressen GeoDataFrame\n",
    "bressen = Bressen.from_gdf(bressen_gdf=bressen_gdf, basisgegevens=basisgegevens)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Verwerken bressen\n",
    "We hanteren de volgende instellingen:\n",
    "\n",
    "- `offset`: offset van de bres t.o.v. kering in meter\n",
    "- `breedte`: breedte van de bres haaks op de kering in meter\n",
    "- `lengte`: lengte van de bres langs de kering in meter\n",
    "\n",
    "\n",
    "### Complete set"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "bressen.bereken_vlakken(offset=50, lengte=200, breedte=40)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Een of enkele bressen \n",
    "Wanneer je 1 of meerdere bressen afzonderlijk wilt processen met andere instellingen, kun je individuele bressen opvragen en als dusdanig processen."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "bres = bressen.get_bres(fid=1329)\n",
    "bres.bereken_vlak(offset=40, lengte=150, breedte=20, basisgegevens=bressen.basisgegevens)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Bekijken resultaat\n",
    "Resultaten kunnen worden opgevraagd als geodataframe. Hier een uitleg van de verschillend resultaat-velden:\n",
    "\n",
    "### Gekozen instellingen\n",
    "Velden met de gekozen instellingen voor het berekenen van het bresvlak:\n",
    "* `bres_offset`: de gekozen offset t.o.v. de kering\n",
    "* `bres_breedte`: de breedte van de bres haaks op de kering\n",
    "* `bres_lengte`: de lengte van de bres in langsrichting van de kering\n",
    "\n",
    "### Onderbouwing van het resultaat\n",
    "Onder bouwing voor het wel, of niet, bepalen van het bresvlak:\n",
    "* `bres_fout`: wanneer dit veld niet leeg is, staat er uitleg waarom het algoritme geen bres kan berekenen. Alle hieronder genoemde velden zijn dan leeg.\n",
    "* `bres_offset_locatie`: de, door het algoritme, gekozen kant van de bres t.o.v. de kering. Wanneer er een watervlak naast de kering ligt, valt de keuze aan de ándere kant van de kering: `naast watervlak`. Zonder watervlak, maar met 2 peilgebieden met verschillende peilen, valt de keuze aan de kant van het `laagste peilgebied`. Is er maar aan 1 kant en peilgebied, dan valt de keuze aan de kant van het bestaande peilgebied, dus `binnen beheergebied`.\n",
    "\n",
    "\n",
    "### Validatie van de resultaten\n",
    "Als `True` in een van deze kolommen, dan is het goed deze bres te bekijken; deze is wellicht niet bruikbaar:\n",
    "* `afwijking_oppervlak`: `True` wanneer `bres_lengte` * `bres_breedte` meer dan 2% afwijkt van het geometrisch oppervlak van het bresvlak\n",
    "* `snijdt_kering`: `True` wanneer het berekende bresvlak een kering snijdt\n",
    "* `snijdt_watervlak`: `True` wanneer het berekende bresvlak een watervlak snijdt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "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>bres_offset</th>\n",
       "      <th>bres_offset_locatie</th>\n",
       "      <th>bres_breedte</th>\n",
       "      <th>bres_lengte</th>\n",
       "      <th>bres_fout</th>\n",
       "      <th>geometry</th>\n",
       "      <th>afwijking_oppervlak</th>\n",
       "      <th>snijdt_kering</th>\n",
       "      <th>snijdt_watervlak</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>fid</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>50</td>\n",
       "      <td>binnen beheergebied</td>\n",
       "      <td>40</td>\n",
       "      <td>200</td>\n",
       "      <td>None</td>\n",
       "      <td>POLYGON ((129537.703 463435.395, 129539.941 46...</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>50</td>\n",
       "      <td>naast watervlak</td>\n",
       "      <td>40</td>\n",
       "      <td>200</td>\n",
       "      <td>None</td>\n",
       "      <td>POLYGON ((129862.151 463153.456, 129862.272 46...</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>50</td>\n",
       "      <td>naast watervlak</td>\n",
       "      <td>40</td>\n",
       "      <td>200</td>\n",
       "      <td>None</td>\n",
       "      <td>POLYGON ((130152.052 462768.863, 130152.345 46...</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>50</td>\n",
       "      <td>binnen beheergebied</td>\n",
       "      <td>40</td>\n",
       "      <td>200</td>\n",
       "      <td>None</td>\n",
       "      <td>POLYGON ((130600.571 462599.507, 130601.947 46...</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>50</td>\n",
       "      <td>naast watervlak</td>\n",
       "      <td>40</td>\n",
       "      <td>200</td>\n",
       "      <td>None</td>\n",
       "      <td>POLYGON ((130690.266 462167.927, 130690.223 46...</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1545</th>\n",
       "      <td>50</td>\n",
       "      <td>naast watervlak</td>\n",
       "      <td>40</td>\n",
       "      <td>200</td>\n",
       "      <td>None</td>\n",
       "      <td>POLYGON ((130039.882 477605.796, 130079.882 47...</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "      <td>True</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1546</th>\n",
       "      <td>50</td>\n",
       "      <td>naast watervlak</td>\n",
       "      <td>40</td>\n",
       "      <td>200</td>\n",
       "      <td>None</td>\n",
       "      <td>POLYGON ((130138.607 478106.108, 130178.607 47...</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1547</th>\n",
       "      <td>50</td>\n",
       "      <td>naast watervlak</td>\n",
       "      <td>40</td>\n",
       "      <td>200</td>\n",
       "      <td>None</td>\n",
       "      <td>POLYGON ((130136.903 478606.173, 130176.903 47...</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1548</th>\n",
       "      <td>50</td>\n",
       "      <td>naast watervlak</td>\n",
       "      <td>40</td>\n",
       "      <td>200</td>\n",
       "      <td>None</td>\n",
       "      <td>POLYGON ((130135.007 479106.193, 130175.007 47...</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1549</th>\n",
       "      <td>50</td>\n",
       "      <td>naast watervlak</td>\n",
       "      <td>40</td>\n",
       "      <td>200</td>\n",
       "      <td>None</td>\n",
       "      <td>POLYGON ((130128.394 479447.669, 130128.387 47...</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "      <td>False</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1549 rows × 9 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "      bres_offset  bres_offset_locatie  bres_breedte  bres_lengte bres_fout  \\\n",
       "fid                                                                           \n",
       "1              50  binnen beheergebied            40          200      None   \n",
       "2              50      naast watervlak            40          200      None   \n",
       "3              50      naast watervlak            40          200      None   \n",
       "4              50  binnen beheergebied            40          200      None   \n",
       "5              50      naast watervlak            40          200      None   \n",
       "...           ...                  ...           ...          ...       ...   \n",
       "1545           50      naast watervlak            40          200      None   \n",
       "1546           50      naast watervlak            40          200      None   \n",
       "1547           50      naast watervlak            40          200      None   \n",
       "1548           50      naast watervlak            40          200      None   \n",
       "1549           50      naast watervlak            40          200      None   \n",
       "\n",
       "                                               geometry afwijking_oppervlak  \\\n",
       "fid                                                                           \n",
       "1     POLYGON ((129537.703 463435.395, 129539.941 46...               False   \n",
       "2     POLYGON ((129862.151 463153.456, 129862.272 46...               False   \n",
       "3     POLYGON ((130152.052 462768.863, 130152.345 46...               False   \n",
       "4     POLYGON ((130600.571 462599.507, 130601.947 46...               False   \n",
       "5     POLYGON ((130690.266 462167.927, 130690.223 46...               False   \n",
       "...                                                 ...                 ...   \n",
       "1545  POLYGON ((130039.882 477605.796, 130079.882 47...               False   \n",
       "1546  POLYGON ((130138.607 478106.108, 130178.607 47...               False   \n",
       "1547  POLYGON ((130136.903 478606.173, 130176.903 47...               False   \n",
       "1548  POLYGON ((130135.007 479106.193, 130175.007 47...               False   \n",
       "1549  POLYGON ((130128.394 479447.669, 130128.387 47...               False   \n",
       "\n",
       "     snijdt_kering snijdt_watervlak  \n",
       "fid                                  \n",
       "1            False             True  \n",
       "2            False            False  \n",
       "3            False            False  \n",
       "4            False            False  \n",
       "5            False            False  \n",
       "...            ...              ...  \n",
       "1545         False             True  \n",
       "1546         False            False  \n",
       "1547         False            False  \n",
       "1548         False            False  \n",
       "1549         False            False  \n",
       "\n",
       "[1549 rows x 9 columns]"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf = bressen.to_gdf(eigenschappen=False)\n",
    "\n",
    "gdf"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Wegschrijven resultaat\n",
    "\n",
    "In de laatste stap schrijven we de resultaten weg naar een GeoPackage. Deze bevat drie lagen\n",
    "- `bressen` met de originele punten samengevoegd in 1 laag (punt-feature)\n",
    "- `bresvlakken` met de bresvlakken hierboven berekend (polygoon-feature)\n",
    "- `fouten` met alle vlakken die niet kunnen worden berekend (geen geometrie)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "bressen.to_gpkg(DATA_DIR / \"resultaat.gpkg\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.12.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}