Read SGeMS Grid

Read SGeMS Grid file formats.

# sphinx_gallery_thumbnail_number = 2
import os

import pooch

from PVGeo.gslib import SGeMSGridReader

# grid_url = 'http://www.trainingimages.org/uploads/3/4/7/0/34703305/a_wlreferencecat.zip'
url = "https://raw.githubusercontent.com/pyvista/vtk-data/master/Data/A_WLreferenceCAT.sgems"
file_path = pooch.retrieve(url=url, known_hash=None)

grid = SGeMSGridReader().apply(file_path)
grid

Header	Data Arrays
ImageData Information N Cells 78000 N Points 157122 X Bounds 0.000e+00, 2.600e+02 Y Bounds 0.000e+00, 3.000e+02 Z Bounds 0.000e+00, 1.000e+00 Dimensions 261, 301, 2 Spacing 1.000e+00, 1.000e+00, 1.000e+00 N Arrays 1	Name Field Type N Comp Min Max var Cells int64 1 0.000e+00 2.000e+00

warped = grid.cell_data_to_point_data().warp_by_scalar(scale_factor=5)
warped.plot()

Static Scene

Interactive Scene

# grid_url = 'http://www.trainingimages.org/uploads/3/4/7/0/34703305/maules_creek_3d.zip'
url = "https://raw.githubusercontent.com/pyvista/vtk-data/master/Data/Maules_Creek_3D.SGEMS.zip"
file_paths = pooch.retrieve(url=url, known_hash=None, processor=pooch.Unzip())
file_path = [f for f in file_paths if os.path.basename(f) == "Maules_Creek_3D.SGEMS"][0]
file_path

'/home/runner/.cache/pooch/db952eed0c325d1a10fc03d53f46f9a9-Maules_Creek_3D.SGEMS.zip.unzip/Maules_Creek_3D.SGEMS'

grid = SGeMSGridReader().apply(file_path)
grid

Header	Data Arrays
ImageData Information N Cells 5440000 N Points 5551821 X Bounds 0.000e+00, 3.400e+02 Y Bounds 0.000e+00, 2.000e+02 Z Bounds 0.000e+00, 8.000e+01 Dimensions 341, 201, 81 Spacing 1.000e+00, 1.000e+00, 1.000e+00 N Arrays 1	Name Field Type N Comp Min Max facies Cells int64 1 0.000e+00 1.000e+00

grid.plot(categories=True)

Static Scene

Interactive Scene

# grid_url = 'http://www.trainingimages.org/uploads/3/4/7/0/34703305/ti_horizons_continuous.zip'
url = "https://raw.githubusercontent.com/pyvista/vtk-data/master/Data/TI_horizons_continuous.SGEMS.zip"
file_paths = pooch.retrieve(url=url, known_hash=None, processor=pooch.Unzip())
file_path = [f for f in file_paths if os.path.basename(f) == "TI_horizons_continuous.SGEMS"][0]
file_path

'/home/runner/.cache/pooch/cb332143f44792c75fe40091d3efa739-TI_horizons_continuous.SGEMS.zip.unzip/TI_horizons_continuous.SGEMS'

grid = SGeMSGridReader().apply(file_path)
grid.threshold([-4, 1.06]).plot(clim=grid.get_data_range())

Static Scene

Interactive Scene

# grid_url = 'http://www.trainingimages.org/uploads/3/4/7/0/34703305/ti.zip'
url = "https://raw.githubusercontent.com/pyvista/vtk-data/master/Data/ti.sgems.zip"
file_paths = pooch.retrieve(url=url, known_hash=None, processor=pooch.Unzip())
file_path = [f for f in file_paths if os.path.basename(f) == "ti.sgems"][0]
file_path

grid = SGeMSGridReader().apply(file_path)
grid.plot(scalars="photo", cpos="xy", cmap="bone")

Static Scene

Interactive Scene

grid.plot(scalars="seismic", cpos="xy")

Static Scene

Interactive Scene