geoparticle.ops.Intersect

class Intersect(geometries: Sequence[Geometry], rmax: float = 1e-05, name=None)[source]

Bases: Geometry

Pointwise intersection of multiple geometries.

Keeps points from the first geometry that are within rmax of at least one point in every other geometry (common intersection under tolerance).

Usage: - Intersect(g1, g2, g3, …, rmax=1e-5) - Intersect([g1, g2, g3, …], rmax=1e-5)

Shortest import: from geoparticle import Intersect

__init__(geometries: Sequence[Geometry], rmax: float = 1e-05, name=None)[source]

Initialize an Intersect object that computes the intersection of multiple geometries.

Parameters:
  • geometries (Sequence[Geometry]) – Geometries to intersect.

  • rmax (float, optional) – Maximum distance for intersection. Defaults to 1e-5.

  • name (str, optional) – Name of the resulting geometry. Defaults to None.

Methods

__init__(geometries[, rmax, name])

Initialize an Intersect object that computes the intersection of multiple geometries.

check_overlap([tol])

Check if there are overlapping points in the geometry.

clip(keep, *[, plane_name, plane_normal, ...])

Clip geometry by a half-space defined by a named plane or an arbitrary plane.

coord2id(x, y, z)

Find the nearest vertex index/indices and its/their coordinates.

copy([name])

Create a deep copy of the geometry object.

equal(geo)

Check if two geometries are equal based on their coordinates.

get_and_delete(ids)

Extract points by their indices and remove them from the geometry.

get_counter()

Get the current value of the class-wide counter.

intersect(geometries[, rmax, inplace, name])

Keep points from self that are within rmax of at least one point in every other geometry.

load_from(other)

Copy coordinate arrays and dimension from another Geometry into this instance.

mirror(plane_name, plane_pos[, inplace, name])

Mirror the geometry across a principal plane.

plot([ax, ms, alpha])

Plot the geometry points in 2D.

rotate(angle_deg[, axis_direction, ...])

Rotate the geometry around a principal axis or a custom axis.

set_coord(xs, ys, zs)

Set the coordinates of the geometry, broadcasting scalars if necessary.

shift([x, y, z, inplace, name])

Translate the geometry by the given offsets.

stack(axis, n_axis, dl, dimension[, ...])

Stack a planar layer (self) along axis by repeating points at spacing dl.

subtract(geo2[, rmax, inplace, name])

Return points from self that are at least rmax away from any point in geo2.

union(geometries[, inplace, name])

Concatenate this geometry with others and return a new Geometry.

Attributes

counter

flatten_coords

Get the flattened array of coordinates.

matrix_coords

Get the coordinates as a 2D array.

size

Get the number of points in the geometry.

counter = 0