![](\"../images/tutorial55_cover.png\")
"
]
},
{
"cell_type": "markdown",
"id": "technical-times",
"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": "latter-premium",
"metadata": {
"ExecuteTime": {
"end_time": "2022-04-06T06:09:12.146812Z",
"start_time": "2022-04-06T06:09:09.434392Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"WARNING (theano.configdefaults): g++ not available, if using conda: `conda install m2w64-toolchain`\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 warnings\n",
"warnings.filterwarnings(\"ignore\")\n",
"import gemgis as gg\n",
"import geopandas as gpd\n",
"import rasterio\n",
"import gempy as gp\n",
"import pyvista as pv\n",
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"import pandas as pd"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "welsh-ireland",
"metadata": {
"ExecuteTime": {
"end_time": "2022-04-06T06:09:12.240352Z",
"start_time": "2022-04-06T06:09:12.148797Z"
}
},
"outputs": [],
"source": [
"file_path ='data/55_extracting_well_tops_from_pyvista_meshes/'\n",
"gg.download_gemgis_data.download_tutorial_data(filename=\"55_extracting_well_tops_from_pyvista_meshes.zip\", dirpath=file_path)"
]
},
{
"cell_type": "markdown",
"id": "capable-dayton",
"metadata": {},
"source": [
"## Creating GemPy Model\n",
"\n",
"### Loading Topography"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "visible-leisure",
"metadata": {
"ExecuteTime": {
"end_time": "2022-04-06T06:09:12.288364Z",
"start_time": "2022-04-06T06:09:12.242346Z"
}
},
"outputs": [],
"source": [
"topo_raster = rasterio.open(file_path + 'raster1.tif')"
]
},
{
"cell_type": "markdown",
"id": "naval-pittsburgh",
"metadata": {},
"source": [
"### Creating Interfaces"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "recovered-flashing",
"metadata": {
"ExecuteTime": {
"end_time": "2022-04-06T06:09:12.383385Z",
"start_time": "2022-04-06T06:09:12.290349Z"
}
},
"outputs": [
{
"data": {
"text/html": [
"| \n", " | formation | \n", "geometry | \n", "X | \n", "Y | \n", "Z | \n", "
|---|---|---|---|---|---|
| 0 | \n", "Sand1 | \n", "POINT (0.256 264.862) | \n", "0.26 | \n", "264.86 | \n", "353.97 | \n", "
| 1 | \n", "Sand1 | \n", "POINT (10.593 276.734) | \n", "10.59 | \n", "276.73 | \n", "359.04 | \n", "
| 2 | \n", "Sand1 | \n", "POINT (17.135 289.090) | \n", "17.13 | \n", "289.09 | \n", "364.28 | \n", "
| 3 | \n", "Sand1 | \n", "POINT (19.150 293.313) | \n", "19.15 | \n", "293.31 | \n", "364.99 | \n", "
| 4 | \n", "Sand1 | \n", "POINT (27.795 310.572) | \n", "27.80 | \n", "310.57 | \n", "372.81 | \n", "
| \n", " | formation | \n", "dip | \n", "azimuth | \n", "geometry | \n", "X | \n", "Y | \n", "Z | \n", "polarity | \n", "
|---|---|---|---|---|---|---|---|---|
| 0 | \n", "Ton | \n", "30.50 | \n", "180.00 | \n", "POINT (96.471 451.564) | \n", "96.47 | \n", "451.56 | \n", "440.59 | \n", "1 | \n", "
| 1 | \n", "Ton | \n", "30.50 | \n", "180.00 | \n", "POINT (172.761 661.877) | \n", "172.76 | \n", "661.88 | \n", "556.38 | \n", "1 | \n", "
| 2 | \n", "Ton | \n", "30.50 | \n", "180.00 | \n", "POINT (383.074 957.758) | \n", "383.07 | \n", "957.76 | \n", "729.02 | \n", "1 | \n", "
| 3 | \n", "Ton | \n", "30.50 | \n", "180.00 | \n", "POINT (592.356 722.702) | \n", "592.36 | \n", "722.70 | \n", "601.55 | \n", "1 | \n", "
| 4 | \n", "Ton | \n", "30.50 | \n", "180.00 | \n", "POINT (766.586 348.469) | \n", "766.59 | \n", "348.47 | \n", "378.63 | \n", "1 | \n", "
| 5 | \n", "Ton | \n", "30.50 | \n", "180.00 | \n", "POINT (843.907 167.023) | \n", "843.91 | \n", "167.02 | \n", "282.61 | \n", "1 | \n", "
| 6 | \n", "Ton | \n", "30.50 | \n", "180.00 | \n", "POINT (941.846 428.883) | \n", "941.85 | \n", "428.88 | \n", "423.45 | \n", "1 | \n", "
| 7 | \n", "Ton | \n", "30.50 | \n", "180.00 | \n", "POINT (22.142 299.553) | \n", "22.14 | \n", "299.55 | \n", "368.05 | \n", "1 | \n", "
| Header | Data Arrays | ||||||||||||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| \n", " | Surface | \n", "Z | \n", "
|---|---|---|
| 0 | \n", "Sand1 | \n", "465.48 | \n", "
| 1 | \n", "Ton | \n", "408.51 | \n", "
| Header | Data Arrays | ||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| \n", " | Surface | \n", "Z | \n", "
|---|---|---|
| 0 | \n", "Sand1 | \n", "497.21 | \n", "
| 1 | \n", "Ton | \n", "377.33 | \n", "
| 2 | \n", "Sand1 | \n", "357.01 | \n", "
| 3 | \n", "Ton | \n", "313.96 | \n", "