mpas_tools.transects.intersection

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

Based on https://stackoverflow.com/a/26669130/7728169 Find the intersection point as a unit vector between great circle arc from a1 to a2 and from b1 to b2. The arcs should have already have been found to intersect by calling intersects()

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

points (mpas_tools.transects.Vector) – An array of Cartesian points on the unit sphere indicating where the arcs intersect