Read UBC Gravity Observations

Read a UBC gravity observations file

# sphinx_gallery_thumbnail_number = 1
import PVGeo
import pyvista
from pyvista import examples

Download sample data files and keep track of names:

url = 'https://github.com/OpenGeoVis/PVGeo/raw/master/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)
print(grav)

Out:

PolyData (0x7fb546c4d888)
  N Cells:      191
  N Points:     191
  X Bounds:     3.565e+05, 3.702e+05
  Y Bounds:     6.000e+06, 6.015e+06
  Z Bounds:     2.155e+03, 2.865e+03
  N Arrays:     2
grav.plot(render_points_as_spheres=True, point_size=10)
../../_images/sphx_glr_read-grv-obs_001.png

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

Gallery generated by Sphinx-Gallery