![](\"../images/tutorial37_cover.png\")
"
]
},
{
"cell_type": "markdown",
"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,
"metadata": {
"ExecuteTime": {
"end_time": "2021-03-17T12:00:08.306134Z",
"start_time": "2021-03-17T12:00:06.595927Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"C:\\Users\\ale93371\\Anaconda3\\envs\\gemgis\\lib\\site-packages\\gemgis\\gemgis.py:27: UserWarning: Shapely 2.0 is installed, but because PyGEOS is also installed, GeoPandas will still use PyGEOS by default for now. To force to use and test Shapely 2.0, you have to set the environment variable USE_PYGEOS=0. You can do this before starting the Python process, or in your code before importing geopandas:\n",
"\n",
"import os\n",
"os.environ['USE_PYGEOS'] = '0'\n",
"import geopandas\n",
"\n",
"In a future release, GeoPandas will switch to using Shapely by default. If you are using PyGEOS directly (calling PyGEOS functions on geometries from GeoPandas), this will then stop working and you are encouraged to migrate from PyGEOS to Shapely 2.0 (https://shapely.readthedocs.io/en/latest/migration_pygeos.html).\n",
" import geopandas as gpd\n"
]
}
],
"source": [
"import gemgis as gg\n",
"\n",
"file_path ='data/37_delaunay_triangulation_for_isoline_maps/'"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"ExecuteTime": {
"end_time": "2021-03-17T12:00:09.023826Z",
"start_time": "2021-03-17T12:00:08.340150Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Downloading file '37_delaunay_triangulation_for_isoline_maps.zip' from 'https://rwth-aachen.sciebo.de/s/AfXRsZywYDbUF34/download?path=%2F37_delaunay_triangulation_for_isoline_maps.zip' to 'C:\\Users\\ale93371\\Documents\\gemgis\\docs\\getting_started\\tutorial\\data\\37_delaunay_triangulation_for_isoline_maps'.\n"
]
}
],
"source": [
"gg.download_gemgis_data.download_tutorial_data(filename=\"37_delaunay_triangulation_for_isoline_maps.zip\", dirpath=file_path)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Loading Data\n",
"\n",
"The data used for GemGIS is obtained from the [Geological Survey NRW](#https://www.gd.nrw.de/). It will be used under Datenlizenz Deutschland – Namensnennung – Version 2.0 (https://www.govdata.de/dl-de/by-2-0) with © Geowissenschaftliche Daten: Untergrundmodell NRW (2020)."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"ExecuteTime": {
"end_time": "2021-03-17T11:59:49.489803Z",
"start_time": "2021-03-17T11:59:48.815183Z"
}
},
"outputs": [
{
"data": {
"text/html": [
"| \n", " | OBJECTID | \n", "Z | \n", "EINHEIT | \n", "Shape_Leng | \n", "geometry | \n", "
|---|---|---|---|---|---|
| 0 | \n", "1.00 | \n", "-2450 | \n", "gg_kru_b_l_Z50m | \n", "3924.67 | \n", "LINESTRING (32403313.109 5785053.637, 32402917... | \n", "
| 1 | \n", "2.00 | \n", "-2400 | \n", "gg_kru_b_l_Z50m | \n", "26332.90 | \n", "LINESTRING (32410198.859 5781110.785, 32409807... | \n", "
| 2 | \n", "3.00 | \n", "-2350 | \n", "gg_kru_b_l_Z50m | \n", "31104.28 | \n", "LINESTRING (32409587.930 5780538.824, 32408824... | \n", "
| 3 | \n", "4.00 | \n", "-2300 | \n", "gg_kru_b_l_Z50m | \n", "35631.73 | \n", "LINESTRING (32408977.008 5779966.863, 32408808... | \n", "
| 4 | \n", "5.00 | \n", "-2250 | \n", "gg_kru_b_l_Z50m | \n", "41702.52 | \n", "LINESTRING (32407319.922 5779788.672, 32407246... | \n", "
| Header | Data Arrays | ||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| \n", " | OBJECTID | \n", "Z | \n", "EINHEIT | \n", "Shape_Leng | \n", "geometry | \n", "X | \n", "Y | \n", "
|---|---|---|---|---|---|---|---|
| 0 | \n", "1.00 | \n", "-2450.00 | \n", "gg_kru_b_l_Z50m | \n", "3924.67 | \n", "POINT (32403313.109 5785053.637) | \n", "32403313.11 | \n", "5785053.64 | \n", "
| 1 | \n", "1.00 | \n", "-2450.00 | \n", "gg_kru_b_l_Z50m | \n", "3924.67 | \n", "POINT (32402917.820 5785157.270) | \n", "32402917.82 | \n", "5785157.27 | \n", "
| 2 | \n", "1.00 | \n", "-2450.00 | \n", "gg_kru_b_l_Z50m | \n", "3924.67 | \n", "POINT (32402585.375 5785498.707) | \n", "32402585.38 | \n", "5785498.71 | \n", "
| 3 | \n", "1.00 | \n", "-2450.00 | \n", "gg_kru_b_l_Z50m | \n", "3924.67 | \n", "POINT (32402870.727 5785657.359) | \n", "32402870.73 | \n", "5785657.36 | \n", "
| 4 | \n", "1.00 | \n", "-2450.00 | \n", "gg_kru_b_l_Z50m | \n", "3924.67 | \n", "POINT (32403086.758 5786086.730) | \n", "32403086.76 | \n", "5786086.73 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
| 23004 | \n", "733.00 | \n", "400.00 | \n", "gg_kru_b_l_Z50m | \n", "153.94 | \n", "POINT (32474651.758 5704841.910) | \n", "32474651.76 | \n", "5704841.91 | \n", "
| 23005 | \n", "733.00 | \n", "400.00 | \n", "gg_kru_b_l_Z50m | \n", "153.94 | \n", "POINT (32474603.586 5704848.586) | \n", "32474603.59 | \n", "5704848.59 | \n", "
| 23006 | \n", "733.00 | \n", "400.00 | \n", "gg_kru_b_l_Z50m | \n", "153.94 | \n", "POINT (32474597.918 5704823.645) | \n", "32474597.92 | \n", "5704823.64 | \n", "
| 23007 | \n", "734.00 | \n", "400.00 | \n", "gg_kru_b_l_Z50m | \n", "94.02 | \n", "POINT (32479877.676 5703950.059) | \n", "32479877.68 | \n", "5703950.06 | \n", "
| 23008 | \n", "734.00 | \n", "400.00 | \n", "gg_kru_b_l_Z50m | \n", "94.02 | \n", "POINT (32479970.945 5703938.168) | \n", "32479970.95 | \n", "5703938.17 | \n", "
23009 rows × 7 columns
\n", "| PolyData | Information |
|---|---|
| N Cells | 23009 |
| N Points | 23009 |
| N Strips | 0 |
| X Bounds | 3.233e+07, 3.250e+07 |
| Y Bounds | 5.702e+06, 5.798e+06 |
| Z Bounds | -2.450e+03, 4.000e+02 |
| N Arrays | 0 |
| PolyData | Information |
|---|---|
| N Cells | 45536 |
| N Points | 23009 |
| N Strips | 0 |
| X Bounds | 3.233e+07, 3.250e+07 |
| Y Bounds | 5.702e+06, 5.798e+06 |
| Z Bounds | -2.450e+03, 4.000e+02 |
| N Arrays | 0 |
| Header | Data Arrays | ||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| \n", " | geometry | \n", "Z | \n", "
|---|---|---|
| 0 | \n", "LINESTRING Z (32409587.930 5780538.824 -2350.0... | \n", "-2350.00 | \n", "
| 1 | \n", "LINESTRING Z (32407304.336 5777048.086 -2050.0... | \n", "-2050.00 | \n", "
| 2 | \n", "LINESTRING Z (32408748.977 5778005.047 -2200.0... | \n", "-2200.00 | \n", "
| 3 | \n", "LINESTRING Z (32403693.547 5786613.994 -2400.0... | \n", "-2400.00 | \n", "
| 4 | \n", "LINESTRING Z (32404738.664 5782672.480 -2350.0... | \n", "-2350.00 | \n", "
| ... | \n", "... | \n", "... | \n", "
| 36332 | \n", "LINESTRING Z (32472712.875 5705828.297 350.000... | \n", "350.00 | \n", "
| 36333 | \n", "LINESTRING Z (32476371.684 5706520.105 350.000... | \n", "350.00 | \n", "
| 36334 | \n", "LINESTRING Z (32476542.145 5706398.895 350.000... | \n", "350.00 | \n", "
| 36335 | \n", "LINESTRING Z (32476538.465 5706556.277 350.000... | \n", "350.00 | \n", "
| 36336 | \n", "LINESTRING Z (32476540.793 5706593.047 350.000... | \n", "350.00 | \n", "
36337 rows × 2 columns
\n", "| \n", " | geometry | \n", "Z | \n", "
|---|---|---|
| 0 | \n", "LINESTRING Z (32331825.641 5708789.973 -200.00... | \n", "-200.00 | \n", "
| 1 | \n", "LINESTRING Z (32334315.359 5723032.766 -250.00... | \n", "-250.00 | \n", "
| 2 | \n", "LINESTRING Z (32332516.312 5722028.768 -250.00... | \n", "-250.00 | \n", "
| 3 | \n", "LINESTRING Z (32332712.750 5721717.561 -250.00... | \n", "-250.00 | \n", "
| 4 | \n", "LINESTRING Z (32332516.312 5722028.768 -250.00... | \n", "-250.00 | \n", "
| ... | \n", "... | \n", "... | \n", "
| 628 | \n", "LINESTRING Z (32480957.211 5755586.262 200.000... | \n", "200.00 | \n", "
| 629 | \n", "LINESTRING Z (32488784.838 5750893.873 250.000... | \n", "250.00 | \n", "
| 630 | \n", "LINESTRING Z (32489046.086 5750945.794 300.000... | \n", "300.00 | \n", "
| 631 | \n", "LINESTRING Z (32490632.693 5749614.501 250.000... | \n", "250.00 | \n", "
| 632 | \n", "LINESTRING Z (32492043.488 5748252.048 300.000... | \n", "300.00 | \n", "
633 rows × 2 columns
\n", "