![](\"../images/tutorial56_cover.png\")
"
]
},
{
"cell_type": "markdown",
"id": "other-portable",
"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": "classical-virus",
"metadata": {
"ExecuteTime": {
"end_time": "2022-12-10T20:31:04.923848Z",
"start_time": "2022-12-10T20:31:02.213634Z"
}
},
"outputs": [],
"source": [
"import warnings\n",
"warnings.filterwarnings(\"ignore\")\n",
"\n",
"# import gemgis as gg\n",
"\n",
"from segysak.segy import segy_loader\n",
"\n",
"import numpy as np\n",
"\n",
"from matplotlib.colors import ListedColormap\n",
"\n",
"import matplotlib as mpl\n",
"import matplotlib.pyplot as plt\n",
"\n",
"import pyvista as pv\n",
"\n",
"from shapely.geometry import LineString"
]
},
{
"cell_type": "markdown",
"id": "desperate-listing",
"metadata": {},
"source": [
"## Conflict between GeoPandas and Segysak\n",
"\n",
"Importing ```gemgis``` and its dependencies together with ```segysak``` will result in an ```ImportError```. Therefore, please comment the ```segysak``` import, import ```gemgis```, download the data, restart the notebook, comment the gemgis import and uncomment the ```segysak``` import to continue running the notebook. "
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "accomplished-signal",
"metadata": {
"ExecuteTime": {
"end_time": "2022-12-10T20:31:04.939848Z",
"start_time": "2022-12-10T20:31:04.926844Z"
}
},
"outputs": [],
"source": [
"file_path ='data/56_displaying_seismic_data_in_pyvista/'\n",
"# gg.download_gemgis_data.download_tutorial_data(filename=\"56_displaying_seismic_data_in_pyvista.zip\", dirpath=file_path)"
]
},
{
"cell_type": "markdown",
"id": "endless-president",
"metadata": {},
"source": [
"## Loading Seismic Data\n",
"\n",
"The used seismic data can be freely downloaded from https://www.nlog.nl/en/scan-2d-seismic-data and https://www.nlog.nl/en/map-2d-seismic-data. It will be opened using the ```segysak``` package (https://segysak.readthedocs.io/en/latest/).\n",
"\n",
"We are loading the recently acquired SCAN Seismic Line 29 as ```xarray``` object. The remaining SCAN Line 20 as well as the Seismic Lines shot within the Framework of the Californie Geothermal project will be loaded directly as vtk file for the final plotting. All data have been prepared the same way as shown here in this notebook."
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "secondary-revision",
"metadata": {
"ExecuteTime": {
"end_time": "2022-12-10T20:31:16.665209Z",
"start_time": "2022-12-10T20:31:04.942849Z"
}
},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "ab108ca1e6ff4d78a959a6c478b1f2c4",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"HBox(children=(HTML(value=''), FloatProgress(value=0.0, max=19393.0), HTML(value='')))"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"Loading as 2D\n"
]
},
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "18d73615b82a450f99f5ed04e1cfe09c",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"HBox(children=(HTML(value='Converting SEGY'), FloatProgress(value=0.0, max=19393.0), HTML(value='')))"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n"
]
},
{
"data": {
"text/html": [
"<xarray.Dataset>\n",
"Dimensions: (cdp: 19393, twt: 5001)\n",
"Coordinates:\n",
" * cdp (cdp) uint16 100 101 102 103 104 ... 19488 19489 19490 19491 19492\n",
" * twt (twt) float64 0.0 2.0 4.0 6.0 ... 9.996e+03 9.998e+03 1e+04\n",
" cdp_x (cdp) float32 1.697e+05 1.697e+05 1.697e+05 ... 2.164e+05 2.164e+05\n",
" cdp_y (cdp) float32 3.753e+05 3.753e+05 3.753e+05 ... 3.874e+05 3.874e+05\n",
"Data variables:\n",
" data (cdp, twt) float32 0.0 0.03689 -0.05879 -0.3133 ... 0.0 0.0 0.0 0.0\n",
"Attributes: (12/13)\n",
" ns: None\n",
" sample_rate: 2.0\n",
" text: C01 CLIENT: EBN AREA: SCAN AREAS, NETHE...\n",
" measurement_system: m\n",
" d3_domain: None\n",
" epsg: None\n",
" ... ...\n",
" corner_points_xy: None\n",
" source_file: L2EBN2020ASCAN029_PreSTM_final_full_AGC.sgy\n",
" srd: None\n",
" datatype: None\n",
" percentiles: [-4.339022828226371, -3.723384707669676, -0.46175039...\n",
" coord_scalar: -100.0| \n", " | \n", " | data | \n", "cdp_x | \n", "cdp_y | \n", "
|---|---|---|---|---|
| cdp | \n", "twt | \n", "\n", " | \n", " | \n", " |
| 100 | \n", "0.0 | \n", "0.000000 | \n", "169653.796875 | \n", "375267.0000 | \n", "
| 2.0 | \n", "0.036886 | \n", "169653.796875 | \n", "375267.0000 | \n", "|
| 4.0 | \n", "-0.058794 | \n", "169653.796875 | \n", "375267.0000 | \n", "|
| 6.0 | \n", "-0.313349 | \n", "169653.796875 | \n", "375267.0000 | \n", "|
| 8.0 | \n", "-0.186058 | \n", "169653.796875 | \n", "375267.0000 | \n", "|
| ... | \n", "... | \n", "... | \n", "... | \n", "... | \n", "
| 19492 | \n", "9992.0 | \n", "0.000000 | \n", "216394.500000 | \n", "387366.1875 | \n", "
| 9994.0 | \n", "0.000000 | \n", "216394.500000 | \n", "387366.1875 | \n", "|
| 9996.0 | \n", "0.000000 | \n", "216394.500000 | \n", "387366.1875 | \n", "|
| 9998.0 | \n", "0.000000 | \n", "216394.500000 | \n", "387366.1875 | \n", "|
| 10000.0 | \n", "0.000000 | \n", "216394.500000 | \n", "387366.1875 | \n", "
96984393 rows × 3 columns
\n", "| \n", " | \n", " | cdp_x | \n", "cdp_y | \n", "
|---|---|---|---|
| cdp | \n", "twt | \n", "\n", " | \n", " |
| 12000 | \n", "0.0 | \n", "198476.09375 | \n", "381918.71875 | \n", "
| 2.0 | \n", "198476.09375 | \n", "381918.71875 | \n", "|
| 4.0 | \n", "198476.09375 | \n", "381918.71875 | \n", "|
| 6.0 | \n", "198476.09375 | \n", "381918.71875 | \n", "|
| 8.0 | \n", "198476.09375 | \n", "381918.71875 | \n", "|
| ... | \n", "... | \n", "... | \n", "... | \n", "
| 19492 | \n", "9992.0 | \n", "216394.50000 | \n", "387366.18750 | \n", "
| 9994.0 | \n", "216394.50000 | \n", "387366.18750 | \n", "|
| 9996.0 | \n", "216394.50000 | \n", "387366.18750 | \n", "|
| 9998.0 | \n", "216394.50000 | \n", "387366.18750 | \n", "|
| 10000.0 | \n", "216394.50000 | \n", "387366.18750 | \n", "
37472493 rows × 2 columns
\n", "| Header | Data Arrays | ||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|
| Header | Data Arrays | ||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
\n",
"
| \n",
"
|