Read UBC Gravity Observations

Read a UBC gravity observations file

from pyvista import examples

# sphinx_gallery_thumbnail_number = 1
import PVGeo

Download sample data files and keep track of names:

url = 'https://github.com/OpenGeoVis/PVGeo/raw/main/tests/data/Craig-Chile/LdM_grav_obs.grv'
grav_file, _ = examples.downloads._retrieve_file(url, 'LdM_grav_obs.grv')
grav = PVGeo.ubc.GravObsReader().apply(grav_file)
grav
HeaderData Arrays
PolyDataInformation
N Cells191
N Points191
X Bounds3.565e+05, 3.702e+05
Y Bounds6.000e+06, 6.015e+06
Z Bounds2.155e+03, 2.865e+03
N Arrays2
NameFieldTypeN CompMinMax
GravPointsfloat641-1.863e+019.917e+00
ErrPointsfloat6415.000e-025.000e-02


grav.plot(render_points_as_spheres=True, point_size=10)
read grv obs

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

Gallery generated by Sphinx-Gallery