python 将XYZ坐标系转换为WGS84坐标系

要将XYZ坐标系转换为WGS84坐标系，需要进行以下步骤：

1. 将XYZ坐标系转换为大地坐标系（经度、纬度、高度）。可以使用以下公式：

- 经度：lon = atan2(y, x)
- 纬度：lat = atan2(z, sqrt(x^2 + y^2))
- 高度：h = sqrt(x^2 + y^2 + z^2) - a

其中，a为地球半径。

2. 将大地坐标系转换为WGS84坐标系。可以使用以下公式：

- a = 6378137.0
- f = 1 / 298.257223563
- b = a * (1 - f)
- e = sqrt((a^2 - b^2) / a^2)
- N = a / sqrt(1 - e^2 * sin(lat)^2)
- X = (N + h) * cos(lat) * cos(lon)
- Y = (N + h) * cos(lat) * sin(lon)
- Z = (N * (1 - e^2) + h) * sin(lat)

其中，a为地球长半径，f为扁率，b为地球短半径，e为第一偏心率，N为卯酉圈半径，h为高度，(X, Y, Z)为XYZ坐标。

下面是Python代码实现：

import math

def xyz_to_wgs84(x, y, z):
    a = 6378137.0
    f = 1 / 298.257223563
    b = a * (1 - f)
    e = math.sqrt((a**2 - b**2) / a**2)
    lon = math.atan2(y, x)
    lat = math.atan2(z, math.sqrt(x**2 + y**2))
    N = a / math.sqrt(1 - e**2 * math.sin(lat)**2)
    h = math.sqrt(x**2 + y**2 + z**2) - a
    X = (N + h) * math.cos(lat) * math.cos(lon)
    Y = (N + h) * math.cos(lat) * math.sin(lon)
    Z = (N * (1 - e**2) + h) * math.sin(lat)
    return (X, Y, Z)