大地坐标转为大地空间直角坐标 jupyter

首先需要导入所需的库，包括 math 库和 numpy 库。然后定义计算大地空间直角坐标的函数。

import math
import numpy as np

def geodetic_to_rectangular(latitude, longitude, height):
    # WGS84椭球体参数
    a = 6378137.0
    f = 1.0 / 298.257223563
    b = a * (1.0 - f)
    e2 = 1.0 - (b / a) ** 2
    n = a / math.sqrt(1 - e2 * math.sin(latitude) ** 2)

    x = (n + height) * math.cos(latitude) * math.cos(longitude)
    y = (n + height) * math.cos(latitude) * math.sin(longitude)
    z = (n * (1 - e2) + height) * math.sin(latitude)

    return np.array([x, y, z])

其中 latitude 为纬度，longitude 为经度，height 为高度。函数返回一个 numpy 数组，包含三个元素，分别为大地空间直角坐标系下的 x、y、z 坐标。

使用示例：

latitude = math.radians(40)  # 纬度为40度
longitude = math.radians(116)  # 经度为116度
height = 1000  # 高度为1000米

rectangular = geodetic_to_rectangular(latitude, longitude, height)
print(rectangular)  # 输出：[ 4.61048316e+06  3.83202015e+06  4.66247253e+06]