segysak.SeisGeom

class segysak.SeisGeom(xarray_obj)
__init__(xarray_obj)

Methods

__init__(xarray_obj)

calc_corner_points()

Calculate the corner points of the geometry or end points of a 2D line.

fill_cdpna()

Fills NaN cdp locations by fitting known cdp x and y values to the local grid using a planar surface relationshipt.

get_affine_transform()

Calculate the forward iline/xline -> cdp_x, cdp_y Affine transform for Matplotlib using corner point geometry.

get_dead_trace_map([scan, zeros_as_nan])

Scan the vertical axis of a volume to find traces that are all NaN and return an DataArray which maps the all dead traces.

get_measurement_system()

Return measurement_system if present, else None

grid_rotation()

Calculate the rotation of the grid using the sum under the curve method.

interp_line(cdpx, cdpy[, extra, ...])

Select data at x and y coordinates

is_2d()

Returns True if the dataset is 2D peformant else False

is_2dgath()

Returns True if the dataset is 2D peformant and has offset or angle else False

is_3d()

Returns True if the dataset is 3D peformant else False

is_3dgath()

Returns True if the dataset is 3D peformant and has offset or angle else False

is_depth()

Check if seisnc volume is in depth

is_empty()

Check if empty

is_twt()

Check if seisnc volume is in twt

plot_bounds([ax])

Plot survey bbounding box to a new or existing axis

subsample_dims(**dim_kwargs)

Return a dictionary of subsampled dims suitable for xarray.interp.

surface_from_points(points, attr[, left, ...])

Sample a 2D point set with an attribute (like Z) to the seisnc geometry using interpolation.

xysel(cdp_x, cdp_y[, method, sample_dim_name])

Select data at x and y coordinates

zeros_like()

Create a new dataset with the same attributes and coordinates and dimensions but with data filled by zeros.

Attributes

humanbytes

Prints Human Friendly size of Dataset to return the bytes as an int use xarray.Dataset.nbytes