mpas_tools.transects.intersection

mpas_tools.transects.intersection(a1, a2, b1, b2)

Based on https://stackoverflow.com/a/26669130/7728169 Find the intersection point 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 points of two great circle arcs. 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 the a1 and a2 or b1 and b2.

• a2 (mpas_tools.transects.Vector) – Cartesian coordinates of the end points of two great circle arcs. 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 the a1 and a2 or b1 and b2.

• b1 (mpas_tools.transects.Vector) – Cartesian coordinates of the end points of two great circle arcs. 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 the a1 and a2 or b1 and b2.

• b2 (mpas_tools.transects.Vector) – Cartesian coordinates of the end points of two great circle arcs. 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 the a1 and a2 or b1 and b2.

Returns

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