Compute Intersection Between An Object and A GridObject¶

In [1]:

import geolime as geo
from pyproj import CRS

import geolime as geo from pyproj import CRS

In [3]:

dh = geo.datasets.load("rocklea_dome/dh_hyper.geo")

dh = geo.datasets.load("rocklea_dome/dh_hyper.geo")

In [4]:

grid = geo.Voxel(
    name="NewGrid", 
    shape=[75, 150], 
    origin=[545400, 7472700, 430], 
    axis=[[50, 0, 0], [0, 50, 0]]
)

grid = geo.Voxel( name="NewGrid", shape=[75, 150], origin=[545400, 7472700, 430], axis=[[50, 0, 0], [0, 50, 0]] )

In [5]:

gs = geo.GriddedSurface(name="GridSurf", xyz=grid.coords())

gs = geo.GriddedSurface(name="GridSurf", xyz=grid.coords())

In [6]:

gs.set_property_expr(name="XY", expr="sin(X / 2) * sin (Y / 2)")

gs.set_property_expr(name="XY", expr="sin(X / 2) * sin (Y / 2)")

In [7]:

geo.compute_surface_elevation_intersection(
    obj=dh,
    grid_object=gs,
    z_prop="XY",
    new_property_name="intersection",
    ignore_warnings=False
)

geo.compute_surface_elevation_intersection( obj=dh, grid_object=gs, z_prop="XY", new_property_name="intersection", ignore_warnings=False )

In [8]:

geo.plot_2d(
    dh, 
    property="intersection", 
    agg_method="mean", 
    crs=CRS("EPSG:20350"), 
    width=650, 
    height=650, 
    mapbox_zoom=12
)

geo.plot_2d( dh, property="intersection", agg_method="mean", crs=CRS("EPSG:20350"), width=650, height=650, mapbox_zoom=12 )

---

Last update: 2023-10-23