mpas_tools.transects.intersects

mpas_tools.transects.intersects(a1, a2, b1, b2)[source]

Based on https://stackoverflow.com/a/26669130/7728169 Determine if the great circle arc from a1 to a2 intersects that from b1 to b2.

Parameters
  • a1 (mpas_tools.transects.Vector) – Cartesian coordinates of the end point of a great circle arc. The types of the attributes x, y, and z must either be numpy.arrays of identical size for all 4 vectors (in which case intersections are found element-wise), or scalars for at least one of either a1 and a2 or b1 and b2.

  • a2 (mpas_tools.transects.Vector) – Cartesian coordinates of the other end point of a great circle arc.

  • b1 (mpas_tools.transects.Vector) – Cartesian coordinates of an end point of a second great circle arc.

  • b2 (mpas_tools.transects.Vector) – Cartesian coordinates of the other end point of the second great circle arc.

Returns

intersect (numpy.ndarray) – A boolean array of the same size as a1 and a2 or b1 and b2, whichever is greater, indicating if the particular pair of arcs intersects