![](\"../images/tutorial29_cover.png\")
\n",
"\n",
"Source: Powell, D. (1995): Interpretation geologischer Strukturen durch Karten - Eine praktische Anleitung mit Aufgaben und Lösungen, page 15, figure 10 A, Springer Verlag Berlin, Heidelberg, New York, ISBN: 978-3-540-58607-4."
]
},
{
"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:24:43.556208Z",
"start_time": "2021-03-17T11:24:41.410353Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"WARNING (theano.configdefaults): g++ not available, if using conda: `conda install m2w64-toolchain`\n",
"C:\\Users\\ale93371\\Anaconda3\\envs\\test_gempy\\lib\\site-packages\\theano\\configdefaults.py:560: UserWarning: DeprecationWarning: there is no c++ compiler.This is deprecated and with Theano 0.11 a c++ compiler will be mandatory\n",
" warnings.warn(\"DeprecationWarning: there is no c++ compiler.\"\n",
"WARNING (theano.configdefaults): g++ not detected ! Theano will be unable to execute optimized C-implementations (for both CPU and GPU) and will default to Python implementations. Performance will be severely degraded. To remove this warning, set Theano flags cxx to an empty string.\n",
"WARNING (theano.tensor.blas): Using NumPy C-API based implementation for BLAS functions.\n"
]
}
],
"source": [
"import gemgis as gg\n",
"\n",
"file_path ='data/29_calculating_orientations_from_strike_lines/'"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"ExecuteTime": {
"end_time": "2021-03-17T11:24:43.936331Z",
"start_time": "2021-03-17T11:24:43.590212Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Downloading file '29_calculating_orientations_from_strike_lines.zip' from 'https://rwth-aachen.sciebo.de/s/AfXRsZywYDbUF34/download?path=%2F29_calculating_orientations_from_strike_lines.zip' to 'C:\\Users\\ale93371\\Documents\\gemgis\\docs\\getting_started\\tutorial\\data\\29_calculating_orientations_from_strike_lines'.\n"
]
}
],
"source": [
"gg.download_gemgis_data.download_tutorial_data(filename=\"29_calculating_orientations_from_strike_lines.zip\", dirpath=file_path)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Loading Data\n",
"\n",
"A shape file containing the intersection points between topographic contours and layer boundaries is loaded as GeoDataFrame."
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"ExecuteTime": {
"end_time": "2020-12-22T12:23:36.842108Z",
"start_time": "2020-12-22T12:23:34.002670Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"WARNING (theano.configdefaults): g++ not available, if using conda: `conda install m2w64-toolchain`\n",
"C:\\Users\\ale93371\\Anaconda3\\envs\\test_gempy\\lib\\site-packages\\theano\\configdefaults.py:560: UserWarning: DeprecationWarning: there is no c++ compiler.This is deprecated and with Theano 0.11 a c++ compiler will be mandatory\n",
" warnings.warn(\"DeprecationWarning: there is no c++ compiler.\"\n",
"WARNING (theano.configdefaults): g++ not detected ! Theano will be unable to execute optimized C-implementations (for both CPU and GPU) and will default to Python implementations. Performance will be severely degraded. To remove this warning, set Theano flags cxx to an empty string.\n",
"WARNING (theano.tensor.blas): Using NumPy C-API based implementation for BLAS functions.\n"
]
},
{
"data": {
"text/html": [
"| \n", " | id | \n", "Z | \n", "formation | \n", "geometry | \n", "
|---|---|---|---|---|
| 0 | \n", "2 | \n", "400 | \n", "Ton | \n", "POINT (62.513 382.832) | \n", "
| 1 | \n", "3 | \n", "500 | \n", "Ton | \n", "POINT (141.489 555.155) | \n", "
| 2 | \n", "4 | \n", "600 | \n", "Ton | \n", "POINT (191.943 720.949) | \n", "
| 3 | \n", "5 | \n", "700 | \n", "Ton | \n", "POINT (287.541 893.302) | \n", "
| 4 | \n", "5 | \n", "700 | \n", "Ton | \n", "POINT (481.908 893.302) | \n", "
| 5 | \n", "4 | \n", "600 | \n", "Ton | \n", "POINT (594.500 720.942) | \n", "
| 6 | \n", "3 | \n", "500 | \n", "Ton | \n", "POINT (684.515 555.153) | \n", "
| 7 | \n", "1 | \n", "300 | \n", "Ton | \n", "POINT (807.441 213.514) | \n", "
| 8 | \n", "1 | \n", "300 | \n", "Ton | \n", "POINT (878.392 213.516) | \n", "
| 9 | \n", "2 | \n", "400 | \n", "Ton | \n", "POINT (927.489 382.827) | \n", "
| \n", " | id | \n", "Z | \n", "formation | \n", "geometry | \n", "
|---|---|---|---|---|
| 0 | \n", "2 | \n", "400 | \n", "Ton | \n", "POINT (62.513 382.832) | \n", "
| 1 | \n", "3 | \n", "500 | \n", "Ton | \n", "POINT (141.489 555.155) | \n", "
| 2 | \n", "4 | \n", "600 | \n", "Ton | \n", "POINT (191.943 720.949) | \n", "
| 3 | \n", "5 | \n", "700 | \n", "Ton | \n", "POINT (287.541 893.302) | \n", "
| 4 | \n", "5 | \n", "700 | \n", "Ton | \n", "POINT (481.908 893.302) | \n", "
| 5 | \n", "4 | \n", "600 | \n", "Ton | \n", "POINT (594.500 720.942) | \n", "
| 6 | \n", "3 | \n", "500 | \n", "Ton | \n", "POINT (684.515 555.153) | \n", "
| 7 | \n", "1 | \n", "300 | \n", "Ton | \n", "POINT (807.441 213.514) | \n", "
| 8 | \n", "1 | \n", "300 | \n", "Ton | \n", "POINT (878.392 213.516) | \n", "
| 9 | \n", "2 | \n", "400 | \n", "Ton | \n", "POINT (927.489 382.827) | \n", "
| \n", " | index | \n", "id | \n", "Z | \n", "formation | \n", "geometry | \n", "
|---|---|---|---|---|---|
| 0 | \n", "7 | \n", "1 | \n", "300 | \n", "Ton | \n", "LINESTRING (807.441 213.514, 878.392 213.516) | \n", "
| 1 | \n", "0 | \n", "2 | \n", "400 | \n", "Ton | \n", "LINESTRING (62.513 382.832, 927.489 382.827) | \n", "
| 2 | \n", "1 | \n", "3 | \n", "500 | \n", "Ton | \n", "LINESTRING (141.489 555.155, 684.515 555.153) | \n", "
| 3 | \n", "2 | \n", "4 | \n", "600 | \n", "Ton | \n", "LINESTRING (191.943 720.949, 594.500 720.942) | \n", "
| 4 | \n", "3 | \n", "5 | \n", "700 | \n", "Ton | \n", "LINESTRING (287.541 893.302, 481.908 893.302) | \n", "
| \n", " | dip | \n", "azimuth | \n", "Z | \n", "geometry | \n", "polarity | \n", "formation | \n", "X | \n", "Y | \n", "
|---|---|---|---|---|---|---|---|---|
| 0 | \n", "30.57 | \n", "180.00 | \n", "350.00 | \n", "POINT (668.959 298.172) | \n", "1.00 | \n", "Ton | \n", "668.96 | \n", "298.17 | \n", "
| 1 | \n", "30.13 | \n", "180.00 | \n", "450.00 | \n", "POINT (454.001 468.992) | \n", "1.00 | \n", "Ton | \n", "454.00 | \n", "468.99 | \n", "
| 2 | \n", "31.10 | \n", "180.00 | \n", "550.00 | \n", "POINT (403.112 638.050) | \n", "1.00 | \n", "Ton | \n", "403.11 | \n", "638.05 | \n", "
| 3 | \n", "30.12 | \n", "180.00 | \n", "650.00 | \n", "POINT (388.973 807.124) | \n", "1.00 | \n", "Ton | \n", "388.97 | \n", "807.12 | \n", "