Read SGeMS Grid

Read SGeMS Grid file formats.

# sphinx_gallery_thumbnail_number = 2
from pyvista import examples

from PVGeo.gslib import SGeMSGridReader
# grid_url = 'http://www.trainingimages.org/uploads/3/4/7/0/34703305/a_wlreferencecat.zip'
filename, _ = examples.downloads._download_file('A_WLreferenceCAT.sgems')

grid = SGeMSGridReader().apply(filename)
grid
HeaderData Arrays
UniformGridInformation
N Cells78000
N Points157122
X Bounds0.000e+00, 2.600e+02
Y Bounds0.000e+00, 3.000e+02
Z Bounds0.000e+00, 1.000e+00
Dimensions261, 301, 2
Spacing1.000e+00, 1.000e+00, 1.000e+00
N Arrays1
NameFieldTypeN CompMinMax
varCellsint6410.000e+002.000e+00


warped = grid.cell_data_to_point_data().warp_by_scalar(scale_factor=5)
warped.plot()
read sgems grid
# grid_url = 'http://www.trainingimages.org/uploads/3/4/7/0/34703305/maules_creek_3d.zip'
filename, _ = examples.downloads._download_file('Maules_Creek_3D.SGEMS.zip')

grid = SGeMSGridReader().apply(filename)
grid
HeaderData Arrays
UniformGridInformation
N Cells5440000
N Points5551821
X Bounds0.000e+00, 3.400e+02
Y Bounds0.000e+00, 2.000e+02
Z Bounds0.000e+00, 8.000e+01
Dimensions341, 201, 81
Spacing1.000e+00, 1.000e+00, 1.000e+00
N Arrays1
NameFieldTypeN CompMinMax
faciesCellsint6410.000e+001.000e+00


grid.plot(categories=True)
read sgems grid
# grid_url = 'http://www.trainingimages.org/uploads/3/4/7/0/34703305/ti_horizons_continuous.zip'
filename, _ = examples.downloads._download_file('TI_horizons_continuous.SGEMS.zip')

grid = SGeMSGridReader().apply(filename)
grid.threshold([-4, 1.06]).plot(clim=grid.get_data_range())
read sgems grid
# grid_url = 'http://www.trainingimages.org/uploads/3/4/7/0/34703305/ti.zip'
filename, _ = examples.downloads._download_file('ti.sgems.zip')

grid = SGeMSGridReader().apply(filename)
grid.plot(scalars='photo', cpos='xy', cmap='bone')
read sgems grid
grid.plot(scalars='seismic', cpos='xy')
read sgems grid

Total running time of the script: ( 0 minutes 38.896 seconds)

Gallery generated by Sphinx-Gallery