![](\"../images/tutorial12_cover.png\")
"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Data Preparation\n",
"\n",
"The data used for GemGIS is obtained from [OpenDataNRW](#https://www.opengeodata.nrw.de/produkte/geologie/geologie/GK/ISGK100/GK100analog/). It will be used under Datenlizenz Deutschland – Namensnennung – Version 2.0 (https://www.govdata.de/dl-de/by-2-0) with © Geowissenschaftliche Daten: Analoges Kartenwerk der Geologischen Karte von Nordrhein-Westfalen 1:100.000 (2020).\n",
"\n",
"### Geological Maps\n",
"For the visualization of cross sections, the geological maps ``Münster`` and ``Rheine`` were taken as shown below.\n",
"\n",
"![](\"../images/geological_maps.png\")
\n",
"\n",
"### Geological Cross Section\n",
"From these maps, the geological cross sections were extracted with an image editing software, tilted so that they are completely horizontal and cropped to the extent of the cross section. The axes were also cut off.\n",
"\n",
"![](\"../images/profile1.jpg\")
\n",
"![](\"../images/profile2.jpg\")
\n",
"\n",
"### Traces of Geological Cross Sections\n",
"\n",
"In order to locate the cross sections in space, the traces of the profiles need to be digitized and the vertical ranges of the cross sections need to be provided as attributes of the shape file. This can be done for example in QGIS.\n",
"\n",
"![](\"../images/traces.PNG\")
\n",
"\n",
"![](\"../images/table.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-17T11:08:49.176819Z",
"start_time": "2021-03-17T11:08:46.913311Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"C:\\Users\\ale93371\\Anaconda3\\envs\\pygeomechanical\\lib\\site-packages\\numpy\\_distributor_init.py:30: UserWarning: loaded more than 1 DLL from .libs:\n",
"C:\\Users\\ale93371\\Anaconda3\\envs\\pygeomechanical\\lib\\site-packages\\numpy\\.libs\\libopenblas.FB5AE2TYXYH2IJRDKGDGQ3XBKLKTF43H.gfortran-win_amd64.dll\n",
"C:\\Users\\ale93371\\Anaconda3\\envs\\pygeomechanical\\lib\\site-packages\\numpy\\.libs\\libopenblas64__v0.3.23-246-g3d31191b-gcc_10_3_0.dll\n",
" warnings.warn(\"loaded more than 1 DLL from .libs:\"\n"
]
}
],
"source": [
"import gemgis as gg\n",
"\n",
"file_path ='data/12_visualizing_cross_sections_in_pyvista/'"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"ExecuteTime": {
"end_time": "2021-03-17T11:08:49.288513Z",
"start_time": "2021-03-17T11:08:49.226871Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Downloading file '12_visualizing_cross_sections_in_pyvista.zip' from 'https://rwth-aachen.sciebo.de/s/AfXRsZywYDbUF34/download?path=%2F/12_visualizing_cross_sections_in_pyvista.zip' to 'C:\\Users\\ale93371\\Documents\\gemgis1\\docs\\getting_started\\tutorial\\data\\12_visualizing_cross_sections_in_pyvista'.\n"
]
}
],
"source": [
"gg.download_gemgis_data.download_tutorial_data(filename=\"12_visualizing_cross_sections_in_pyvista.zip\", dirpath=file_path)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Loading Data\n",
"\n",
"The traces were saved as shape file and are now loaded as GeoDataFrame using GeoPandas."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"ExecuteTime": {
"end_time": "2021-01-02T08:04:56.164520Z",
"start_time": "2021-01-02T08:04:51.877307Z"
}
},
"outputs": [
{
"data": {
"text/html": [
"| \n", " | id | \n", "zmax | \n", "zmin | \n", "name | \n", "geometry | \n", "
|---|---|---|---|---|---|
| 0 | \n", "NaN | \n", "500 | \n", "-6000 | \n", "Muenster | \n", "LINESTRING (32386148.890 5763304.720, 32393549... | \n", "
| 1 | \n", "NaN | \n", "500 | \n", "-2000 | \n", "Rheine | \n", "LINESTRING (32402275.136 5761828.501, 32431165... | \n", "
| PolyData | Information |
|---|---|
| N Cells | 20 |
| N Points | 22 |
| N Strips | 0 |
| X Bounds | 3.239e+07, 3.242e+07 |
| Y Bounds | 5.717e+06, 5.763e+06 |
| Z Bounds | -6.000e+03, 5.000e+02 |
| N Arrays | 0 |