Group of Surfaces

In [1]:

# Surface Group

# Surface Group

In [2]:

import numpy as np
import pandas as pd
from volumetricspy import Surface, SurfaceGroup
import matplotlib.pyplot as plt
import seaborn as sns

import numpy as np import pandas as pd from volumetricspy import Surface, SurfaceGroup import matplotlib.pyplot as plt import seaborn as sns

In [3]:

s1 = Surface(name='anticline')
s1.from_z_map('anticline',crs=3117)

sg1 = SurfaceGroup(name='Test')

sg1.add_surface(s1)

s1 = Surface(name='anticline') s1.from_z_map('anticline',crs=3117) sg1 = SurfaceGroup(name='Test') sg1.add_surface(s1)

In [4]:

sg1.surfaces['anticline'].contourf()

sg1.surfaces['anticline'].contourf()

Out[4]:

<matplotlib.contour.QuadContourSet at 0x7fc490a43490>

In [5]:

sg1.create_parallel_surfaces('anticline',[-180,-250,-320])

sg1.create_parallel_surfaces('anticline',[-180,-250,-320])

In [6]:

sg1.surfaces.keys()

sg1.surfaces.keys()

Out[6]:

dict_keys(['anticline', 'anticline_-180.0', 'anticline_-250.0', 'anticline_-320.0'])

In [7]:

sg1.surfaces['anticline_-180.0'].contourf()

sg1.surfaces['anticline_-180.0'].contourf()

Out[7]:

<matplotlib.contour.QuadContourSet at 0x7fc490acb4c0>

In [8]:

rv, area = sg1.get_volume_bounds(
    top_surface = 'anticline',
    bottom_surface = 'anticline_-180.0',
    method ='contours',
    zmin = -5400,zmax=-4900
)
rv

rv, area = sg1.get_volume_bounds( top_surface = 'anticline', bottom_surface = 'anticline_-180.0', method ='contours', zmin = -5400,zmax=-4900 ) rv

Out[8]:

800015.9260285863

In [9]:

fig, ax = plt.subplots(1,1,figsize=(10,10))
sns.lineplot(data=area.reset_index(),x='index',y='area_top',ax=ax)
sns.lineplot(data=area.reset_index(),x='index',y='area_bottom',ax=ax)

fig, ax = plt.subplots(1,1,figsize=(10,10)) sns.lineplot(data=area.reset_index(),x='index',y='area_top',ax=ax) sns.lineplot(data=area.reset_index(),x='index',y='area_bottom',ax=ax)

Out[9]:

<AxesSubplot:xlabel='index', ylabel='area_top'>

In [10]:

surv1_grid = sg1.surfaces['anticline'].make_grid(dz=-100,nz=5)

surv1_grid = sg1.surfaces['anticline'].make_grid(dz=-100,nz=5)

In [11]:

#surv1_vtk.plot()

#surv1_vtk.plot()

In [12]:

surv_grid_vtk = surv1_grid.unstructured_grid_vtk()

surv_grid_vtk = surv1_grid.unstructured_grid_vtk()

In [17]:

surv_grid_vtk.save('test.vtk')

surv_grid_vtk.save('test.vtk')

In [28]:

surv_grid_vtk

surv_grid_vtk

Out[28]:

Header	Data Arrays
UnstructuredGrid Information N Cells 71868 N Points 574944 X Bounds 1.102e+06, 1.108e+06 Y Bounds 1.012e+06, 1.020e+06 Z Bounds -6.058e+03, -4.708e+03 N Arrays 1	Name Field Type N Comp Min Max Elevation Points float32 1 -6.058e+03 -4.708e+03

In [32]:

surv1_grid.nz

surv1_grid.nz

Out[32]:

4

In [ ]: