Skip to content

stemflow.utils.sphere.coordinate_transform

lonlat_cartesian_3D_transformer

Transformer between longitude,latitude and 3d dimension (x,y,z).

Source code in stemflow/utils/sphere/coordinate_transform.py 
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
class lonlat_cartesian_3D_transformer:
    """Transformer between longitude,latitude and 3d dimension (x,y,z)."""

    def __init__(self) -> None:
        pass

    def transform(lng: np.ndarray, lat: np.ndarray, radius: float = 6371.0) -> Tuple[np.ndarray, np.ndarray]:
        """Transform lng, lat to x,y,z

        Args:
            lng (np.ndarray): lng
            lat (np.ndarray): lat
            radius (float, optional): radius of earth in km. Defaults to 6371.

        Returns:
            Tuple[np.ndarray, np.ndarray]: x,y,z
        """

        # Convert latitude and longitude from degrees to radians
        lat_rad = np.radians(lat)
        lng_rad = np.radians(lng)

        # Calculate Cartesian coordinates
        x = radius * np.cos(lat_rad) * np.cos(lng_rad)
        y = radius * np.cos(lat_rad) * np.sin(lng_rad)
        z = radius * np.sin(lat_rad)

        return x, y, z

    def inverse_transform(
        x: np.ndarray, y: np.ndarray, z: np.ndarray, r: float = None
    ) -> Tuple[np.ndarray, np.ndarray]:
        """transform x,y,z to lon, lat

        Args:
            x (np.ndarray): x
            y (np.ndarray): y
            z (np.ndarray): z
            r (float, optional): Radius of your spherical coordinate. If not given, calculate from x,y,z. Defaults to None.

        Returns:
            Tuple[np.ndarray, np.ndarray]: longitude, latitude
        """
        if r is None:
            r = np.sqrt(x**2 + y**2 + z**2)
        latitude = np.degrees(np.arcsin(z / r))
        longitude = np.degrees(np.arctan2(y, x))
        return longitude, latitude

inverse_transform(x, y, z, r=None)

transform x,y,z to lon, lat

Parameters:

  • x (ndarray) –

    x

  • y (ndarray) –

    y

  • z (ndarray) –

    z

  • r (float, default: None ) –

    Radius of your spherical coordinate. If not given, calculate from x,y,z. Defaults to None.

Returns:

  • Tuple[ndarray, ndarray] –

    Tuple[np.ndarray, np.ndarray]: longitude, latitude

Source code in stemflow/utils/sphere/coordinate_transform.py 
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
def inverse_transform(
    x: np.ndarray, y: np.ndarray, z: np.ndarray, r: float = None
) -> Tuple[np.ndarray, np.ndarray]:
    """transform x,y,z to lon, lat

    Args:
        x (np.ndarray): x
        y (np.ndarray): y
        z (np.ndarray): z
        r (float, optional): Radius of your spherical coordinate. If not given, calculate from x,y,z. Defaults to None.

    Returns:
        Tuple[np.ndarray, np.ndarray]: longitude, latitude
    """
    if r is None:
        r = np.sqrt(x**2 + y**2 + z**2)
    latitude = np.degrees(np.arcsin(z / r))
    longitude = np.degrees(np.arctan2(y, x))
    return longitude, latitude

transform(lng, lat, radius=6371.0)

Transform lng, lat to x,y,z

Parameters:

  • lng (ndarray) –

    lng

  • lat (ndarray) –

    lat

  • radius (float, default: 6371.0 ) –

    radius of earth in km. Defaults to 6371.

Returns:

  • Tuple[ndarray, ndarray] –

    Tuple[np.ndarray, np.ndarray]: x,y,z

Source code in stemflow/utils/sphere/coordinate_transform.py 
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
def transform(lng: np.ndarray, lat: np.ndarray, radius: float = 6371.0) -> Tuple[np.ndarray, np.ndarray]:
    """Transform lng, lat to x,y,z

    Args:
        lng (np.ndarray): lng
        lat (np.ndarray): lat
        radius (float, optional): radius of earth in km. Defaults to 6371.

    Returns:
        Tuple[np.ndarray, np.ndarray]: x,y,z
    """

    # Convert latitude and longitude from degrees to radians
    lat_rad = np.radians(lat)
    lng_rad = np.radians(lng)

    # Calculate Cartesian coordinates
    x = radius * np.cos(lat_rad) * np.cos(lng_rad)
    y = radius * np.cos(lat_rad) * np.sin(lng_rad)
    z = radius * np.sin(lat_rad)

    return x, y, z

continuous_interpolation_3D_plotting(p1, p2, radius=6371.0)

interpolate 10 points on earth surface between the given two points. For plotting.

Parameters:

  • p1 (ndarray) –

    p1

  • p2 (ndarray) –

    p2

  • radius (float, default: 6371.0 ) –

    radius of earth in km. Defaults to 6371.0.

Returns:

  • Tuple[ndarray, ndarray, ndarray] –

    Tuple[np.ndarray, np.ndarray, np.ndarray]: 10 x, 10 y, 10 z

Source code in stemflow/utils/sphere/coordinate_transform.py 
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
def continuous_interpolation_3D_plotting(
    p1: np.ndarray, p2: np.ndarray, radius: float = 6371.0
) -> Tuple[np.ndarray, np.ndarray, np.ndarray]:
    """interpolate 10 points on earth surface between the given two points. For plotting.

    Args:
        p1 (np.ndarray): p1
        p2 (np.ndarray): p2
        radius (float, optional): radius of earth in km. Defaults to 6371.0.

    Returns:
        Tuple[np.ndarray, np.ndarray, np.ndarray]: 10 x, 10 y, 10 z
    """
    v1 = np.array([p1[0], p1[1], p1[2]])
    v2 = np.array([p2[0], p2[1], p2[2]])

    x_ = []
    y_ = []
    z_ = []
    for bins_ in np.linspace(0, 1, 10):
        v3 = v1 * bins_ + v2 * (1 - bins_)
        v3 = v3 * (radius / np.linalg.norm(v3))

        x_.append(v3[0])
        y_.append(v3[1])
        z_.append(v3[2])

    return np.array(x_), np.array(y_), np.array(z_)

get_midpoint_3D(p1, p2, radius=6371.0)

Get the mid-point of three QPoint_3D objet (vector)

Parameters:

  • p1 (QPoint_3D) –

    p1

  • p2 (QPoint_3D) –

    p2

  • radius (float, default: 6371.0 ) –

    radius of earth in km. Defaults to 6371.0.

Returns:

Source code in stemflow/utils/sphere/coordinate_transform.py 
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
def get_midpoint_3D(p1: QPoint_3D, p2: QPoint_3D, radius: float = 6371.0) -> QPoint_3D:
    """Get the mid-point of three QPoint_3D objet (vector)

    Args:
        p1 (QPoint_3D): p1
        p2 (QPoint_3D): p2
        radius (float, optional): radius of earth in km. Defaults to 6371.0.

    Returns:
        QPoint_3D: mid-point.
    """
    v1 = np.array([p1.x, p1.y, p1.z])
    v2 = np.array([p2.x, p2.y, p2.z])

    v3 = v1 + v2
    v3 = v3 * (radius / np.linalg.norm(v3))

    p3 = QPoint_3D(None, v3[0], v3[1], v3[2])

    return p3