![](\"../images/tutorial49_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:15:39.645345Z",
"start_time": "2021-03-17T12:15:37.457090Z"
}
},
"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/49_slicing_gempy_lith_blocks_in_pyvista_with_gemgis/'"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"ExecuteTime": {
"end_time": "2021-03-17T12:15:39.961870Z",
"start_time": "2021-03-17T12:15:39.678305Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Downloading file '49_slicing_gempy_lith_blocks_in_pyvista_with_gemgis.zip' from 'https://rwth-aachen.sciebo.de/s/AfXRsZywYDbUF34/download?path=%2F49_slicing_gempy_lith_blocks_in_pyvista_with_gemgis.zip' to 'C:\\Users\\ale93371\\Documents\\gemgis\\docs\\getting_started\\tutorial\\data\\49_slicing_gempy_lith_blocks_in_pyvista_with_gemgis'.\n"
]
}
],
"source": [
"gg.download_gemgis_data.download_tutorial_data(filename=\"49_slicing_gempy_lith_blocks_in_pyvista_with_gemgis.zip\", dirpath=file_path)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Loading the lith block\n",
"\n",
"The lith block of a computed GemPy model is loaded using NumPy and reshaped to the original dimensions (resolution) of the model."
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"ExecuteTime": {
"end_time": "2021-01-22T12:58:55.924016Z",
"start_time": "2021-01-22T12:58:53.162484Z"
}
},
"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/plain": [
"array([[1., 1., 1., ..., 1., 1., 1.],\n",
" [1., 1., 1., ..., 1., 1., 1.],\n",
" [1., 1., 1., ..., 1., 1., 1.],\n",
" ...,\n",
" [3., 3., 3., ..., 2., 2., 1.],\n",
" [3., 3., 3., ..., 2., 2., 2.],\n",
" [3., 3., 3., ..., 2., 2., 2.]])"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"import pyvista as pv\n",
"import numpy as np\n",
"\n",
"lith_block = np.load(file_path + 'lith_block.npy').reshape(50,50,50)\n",
"lith_block[0]"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"All values above the topography representing air are assigned NaN values so that they can be hidden from the visualization."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"ExecuteTime": {
"end_time": "2021-01-22T12:58:55.939538Z",
"start_time": "2021-01-22T12:58:55.926018Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"array([], shape=(0, 50, 50), dtype=float64)"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"lith_block[lith_block == 1] = np.nan\n",
"lith_block[:0]"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Wrapping the object\n",
"\n",
"The array can now be wrapped using ``pv.wrap(..)`` to create a volume."
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"ExecuteTime": {
"end_time": "2021-01-22T12:58:55.955541Z",
"start_time": "2021-01-22T12:58:55.941538Z"
}
},
"outputs": [
{
"data": {
"text/html": [
"| Header | Data Arrays | ||||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Information | Blocks | ||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Information | Blocks | ||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Information | Blocks | ||||||||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Information | Blocks | ||||||||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|