mpas_tools.transects.subdivide_great_circle

mpas_tools.transects.subdivide_great_circle(x, y, z, maxRes, earthRadius)[source]

Subdivide each segment of the transect so the horizontal resolution approximately matches the requested resolution

Uses a formula for interpolating unit vectors on the sphere from https://en.wikipedia.org/wiki/Slerp

Parameters:

  • x (numpy.ndarray) – The Cartesian x coordinate of a transect, where the number of segments is len(x) - 1. x, y and z are of the same length.

  • y (numpy.ndarray) – The Cartesian y coordinate of the transect

  • z (numpy.ndarray) – The Cartesian z coordinate of the transect

  • maxRes (float) – The maximum allowed spacing in m after subdivision

  • earthRadius (float) – The radius of the Earth in m

Returns:

  • xOut (numpy.ndarray) – The Cartesian x values of the transect subdivided into segments with segment length at most maxRes. All the points in x, y and z are guaranteed to be included.

  • yOut (numpy.ndarray) – The Cartesian y values of the subdivided transect

  • zOut (numpy.ndarray) – The Cartesian y values of the subdivided transect

  • dIn (numpy.ndarray) – The distance along the transect before subdivision

  • dOut (numpy.ndarray) – The distance along the transect after subdivision