{ "cells": [ { "cell_type": "markdown", "id": "aea86869", "metadata": {}, "source": [ "# 66 Generating Voronoi Polygons\n", "\n", "This notebook illustrates the creation of Thiessen Polygons to illustrate the the data quantity of borehole data according to page 236 of the following report of the BGE: https://www.bge.de/fileadmin/user_upload/Standortsuche/Wesentliche_Unterlagen/Methodik/Phase_I_Schritt_2/rvSU-Methodik/20220328_Anlage_zu_rvSU_Konzept_Methodenbeschreibung_barrierefrei.pdf\n" ] }, { "cell_type": "markdown", "id": "448e8438", "metadata": {}, "source": [ "## Set File Paths and download Tutorial Data\n", "\n", "If you downloaded the latest `GemGIS` version from the Github repository, append the path so that the package can be imported successfully. Otherwise, it is recommended to install `GemGIS` via `pip install gemgis` and import `GemGIS` using `import gemgis as gg`. In addition, the file path to the folder where the data is being stored is set. The tutorial data is downloaded using Pooch (https://www.fatiando.org/pooch/latest/index.html) and stored in the specified folder. Use `pip install pooch` if Pooch is not installed on your system yet." ] }, { "cell_type": "code", "execution_count": 1, "id": "00408871", "metadata": { "ExecuteTime": { "end_time": "2023-08-08T11:44:12.264775Z", "start_time": "2023-08-08T11:44:09.818313Z" } }, "outputs": [], "source": [ "import warnings\n", "warnings.filterwarnings(\"ignore\")\n", "\n", "import gemgis as gg\n", "import geopandas as gpd\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "id": "9657b36d", "metadata": { "ExecuteTime": { "end_time": "2023-08-08T11:44:12.280758Z", "start_time": "2023-08-08T11:44:12.268757Z" } }, "outputs": [], "source": [ "file_path ='data/66_generating_voronoi_polygons/'\n", "# gg.download_gemgis_data.download_tutorial_data(filename=\"66_generating_voronoi_polygons.zip\", dirpath=file_path)" ] }, { "cell_type": "markdown", "id": "221929b8", "metadata": {}, "source": [ "## Loading Borehole Data\n", "\n", "The borehole data was extracted from the borehole database of the Geological Survey of NRW (https://www.gd.nrw.de/gd_archive_dabo.htm)." ] }, { "cell_type": "code", "execution_count": 3, "id": "4fdf0ae1", "metadata": { "ExecuteTime": { "end_time": "2023-08-08T11:44:12.616880Z", "start_time": "2023-08-08T11:44:12.284761Z" } }, "outputs": [ { "data": { "text/html": [ "
\n", " | Index | \n", "DABO No. | \n", "Name | \n", "Number | \n", "Depth | \n", "X | \n", "Y | \n", "Z | \n", "X_GK | \n", "Y_GK | \n", "... | \n", "Kind | \n", "Procedure | \n", "Confidenti | \n", "Record Typ | \n", "Lithlog Ve | \n", "Quality | \n", "Drilling P | \n", "Remarks | \n", "Availabili | \n", "geometry | \n", "
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
0 | \n", "GD10092 | \n", "DABO_102029 | \n", "Klingerhuf1(1982) | \n", "NaN | \n", "1118.00 | \n", "32332544.77 | \n", "5701554.56 | \n", "30.80 | \n", "2541103.04 | \n", "5700806.87 | \n", "... | \n", "Bohrung | \n", "gemischtes Bohrverfahren | \n", "offen; Bohrung mit dokumentiertem Freigabeverm... | \n", "Übertragung eines alten Archivbestandes | \n", "1 | \n", "Schichtdaten von guter Qualität; genaue strati... | \n", "30.9.1982 - 1.12.1982 | \n", "Karbon Cumulus | \n", "Original-Schichtenverzeichnis liegt vor | \n", "POINT (32332544.770 5701554.560) | \n", "
1 | \n", "GD50035 | \n", "DABO_1521 | \n", "Salzbergen2 | \n", "NaN | \n", "1875.70 | \n", "32388961.53 | \n", "5798219.66 | \n", "30.00 | \n", "2593500.00 | \n", "5799735.00 | \n", "... | \n", "Bohrung | \n", "gemischtes Bohrverfahren | \n", "vertraulich, offen nur mit Einwilligung; Eigen... | \n", "Aufnahme durch Sachbearbeiter; überarbeitet na... | \n", "1 | \n", "Schichtdaten von guter Qualität; genaue strati... | \n", "15.2.1949 - 27.5.1950 | \n", "paläontologisch untersucht; Ausführliche Unter... | \n", "Original-Schichtenverzeichnis liegt vor | \n", "POINT (32388961.530 5798219.660) | \n", "
2 | \n", "GD50073 | \n", "DABO_1564 | \n", "B.Salzbergen3 | \n", "3 | \n", "1908.90 | \n", "32390628.91 | \n", "5799824.32 | \n", "36.00 | \n", "2595100.00 | \n", "5801408.00 | \n", "... | \n", "Bohrung | \n", "Spülbohrung (Rotary-Verfahren) | \n", "vertraulich, offen nur mit Einwilligung; Eigen... | \n", "Übertragung eines alten Archivbestandes | \n", "1 | \n", "Schichtdaten von guter Qualität; genaue strati... | \n", "19.10.1950 - 25.5.1951 | \n", "siehe Original-Schichtenverzeichnis | \n", "Original-Schichtenverzeichnis liegt vor | \n", "POINT (32390628.910 5799824.320) | \n", "
3 | \n", "GD50261 | \n", "DABO_2319 | \n", "Ibbenbüren-Westfeld1 | \n", "UB1121 | \n", "1352.00 | \n", "32409115.65 | \n", "5798184.48 | \n", "118.20 | \n", "3409152.00 | \n", "5800063.00 | \n", "... | \n", "Bohrung | \n", "gemischtes Bohrverfahren | \n", "offen; Bohrung mit dokumentiertem Freigabeverm... | \n", "Übertragung eines alten Archivbestandes | \n", "1 | \n", "Schichtdaten von guter Qualität; genaue strati... | \n", "11.11.1974 - 14.1.1975 | \n", "Karbon Cumulus | \n", "Original-Schichtenverzeichnis liegt vor | \n", "POINT (32409115.650 5798184.480) | \n", "
4 | \n", "GD51096 | \n", "DABO_7331 | \n", "GildehausZ1 | \n", "NaN | \n", "1683.10 | \n", "32370464.85 | \n", "5790221.74 | \n", "40.00 | \n", "2575345.00 | \n", "5790975.00 | \n", "... | \n", "Bohrung | \n", "gemischtes Bohrverfahren | \n", "vertraulich, offen nur mit Einwilligung; Eigen... | \n", "Aufnahme durch Sachbearbeiter; überarbeitet na... | \n", "1 | \n", "Schichtdaten von guter Qualität; genaue strati... | \n", "..1956 | \n", "Ausführliche Unterlagen im Original | \n", "Original-Schichtenverzeichnis liegt vor | \n", "POINT (32370464.850 5790221.740) | \n", "
5 rows × 26 columns
\n", "\n", " | osm_id | \n", "class | \n", "type | \n", "name | \n", "address | \n", "extratags | \n", "geometry | \n", "
---|---|---|---|---|---|---|---|
0 | \n", "62761.00 | \n", "boundary | \n", "administrative | \n", "North Rhine-Westphalia, Germany | \n", "\"state\"=>\"North Rhine-Westphalia\",\"ISO3166-2-l... | \n", "\"ref:nuts\"=>\"DEA\",\"wikidata\"=>\"Q1198\",\"wikiped... | \n", "MULTIPOLYGON (((300636.227 5600913.724, 300639... | \n", "
\n", " | geometry | \n", "area | \n", "
---|---|---|
0 | \n", "POLYGON ((32473194.081 5718253.741, 32469287.8... | \n", "665698318.29 | \n", "
1 | \n", "POLYGON ((32473175.498 5773862.463, 32481354.2... | \n", "184545066.31 | \n", "
2 | \n", "POLYGON ((32481354.247 5765667.587, 32481234.9... | \n", "304148453.73 | \n", "
3 | \n", "POLYGON ((32459819.403 5765274.413, 32473175.4... | \n", "490884327.21 | \n", "
4 | \n", "POLYGON ((32317260.014 5640927.969, 32314568.7... | \n", "250801576.02 | \n", "
... | \n", "... | \n", "... | \n", "
1131 | \n", "POLYGON ((32326629.143 5727191.374, 32326669.3... | \n", "922649.29 | \n", "
1132 | \n", "POLYGON ((32326895.182 5726701.188, 32327147.6... | \n", "1407420.98 | \n", "
1133 | \n", "POLYGON ((32327531.454 5727716.570, 32327549.9... | \n", "349756.27 | \n", "
1134 | \n", "POLYGON ((32327677.535 5727156.396, 32328944.7... | \n", "1725230.07 | \n", "
1135 | \n", "POLYGON ((32327531.454 5727716.570, 32327549.9... | \n", "3340904.68 | \n", "
1135 rows × 2 columns
\n", "