geovectorslib.utils

Complexity

Total Complexity

3

Size/Duplication

Total Lines	40
Duplicated Lines	0 %

Importance

Changes

0

3 Functions

Rating	Name	Duplication	Size	Complexity
A	wrap90deg()	0	9	1
A	wrap180deg()	0	14	1
A	wrap360deg()	0	9	1

import numpy as np


def wrap90deg(deg: np.ndarray) -> np.ndarray:
    """
    Make sure that degrees (e.g. lats) is within -90;90 range.

    Ref: https://www.movable-type.co.uk/scripts/latlong-vincenty.html
    """
    mask = (deg < -90) | (deg > 90)
    deg[mask] = np.abs((deg[mask] % 360 + 270) % 360 - 180) - 90
    return deg


def wrap180deg(deg: np.ndarray) -> np.ndarray:
    """
    Make sure that degrees (e.g. lon) is within +/-(0-180) range.

    Ref: https://www.movable-type.co.uk/scripts/latlong-vincenty.html
    """
    # experimental value to fix https://github.com/omdv/geovectors/issues/2
    eps = 1e-6

    mask = (deg <= -180) | (deg >= 180)
    deg[mask] = (deg[mask] + 540) % 360 - 180

    deg = np.where(deg > 180 - eps, -180, deg)
    return deg


def wrap360deg(deg: np.ndarray) -> np.ndarray:
    """
    Make sure that degrees (e.g. heading) is within 360 range.

    Ref: https://www.movable-type.co.uk/scripts/latlong-vincenty.html
    """
    mask = (deg < 0) | (deg > 360)
    deg[mask] = (deg[mask] % 360 + 360) % 360
    return deg