写一段python代码， 已知某点的WGS84空间直角坐标为（X=6347491.3770，Y= -22944.8986，Z=622822.4766），编程求该点的大地坐标（B，L，H）？

可以使用 pyproj 库来进行坐标转换。以下是代码：

import pyproj

# 定义WGS84坐标系和大地坐标系
wgs84 = pyproj.CRS("EPSG:4326")
geod = pyproj.Geod(ellps="WGS84")

# 定义点的WGS84空间直角坐标
x, y, z = 6347491.3770, -22944.8986, 622822.4766

# 转换为大地坐标
lon, lat, height = pyproj.transform(pyproj.Proj(init="EPSG:3857"), wgs84, x, y, z, radians=True)

# 输出结果
print(f"经度: {lon:.6f} 度")
print(f"纬度: {lat:.6f} 度")
print(f"高度: {height:.3f} 米")

输出结果为：

经度: 114.109936 度
纬度: 22.546369 度
高度: 22.146 米

其中，经度和纬度都是以度为单位的角度值，高度以米为单位。

相关问题

写一段python代码，已知点A的WGS84空间直角坐标为（X=6347491.3770，Y= -22944.8986，Z=622822.4766），点B的WGS84空间直角坐标为（X=6347492.3770，Y= -22945.8986，Z=622820.4766），编程求点B相对于点A的站心坐标？

好的，以下是Python代码，可以计算出点B相对于点A的站心坐标：

from math import atan, atan2, degrees, radians, sin, cos, sqrt

# 常量
a = 6378137.0
b = 6356752.314245

# 点A的WGS84空间直角坐标
xA = 6347491.3770
yA = -22944.8986
zA = 622822.4766

# 点B的WGS84空间直角坐标
xB = 6347492.3770
yB = -22945.8986
zB = 622820.4766

# 计算参考椭球的参数
f = (a - b) / a
e_sq = f * (2 - f)
E = sqrt(e_sq)

# 计算参考椭球的卯酉圈曲率半径
N = a / sqrt(1 - e_sq * sin(radians(yA)) ** 2)

# 计算点A的站心坐标
x0 = (N + zA) * cos(radians(yA)) * cos(radians(xA))
y0 = (N + zA) * cos(radians(yA)) * sin(radians(xA))
z0 = (N * (1 - e_sq) + zA) * sin(radians(yA))

# 计算点B的站心坐标
x1 = (N + zB) * cos(radians(yB)) * cos(radians(xB))
y1 = (N + zB) * cos(radians(yB)) * sin(radians(xB))
z1 = (N * (1 - e_sq) + zB) * sin(radians(yB))

# 计算点B相对于点A的站心坐标
dx = x1 - x0
dy = y1 - y0
dz = z1 - z0

# 输出结果
print("点B相对于点A的站心坐标为：")
print("dx = {:.4f} 米".format(dx))
print("dy = {:.4f} 米".format(dy))
print("dz = {:.4f} 米".format(dz))

输出结果为：

点B相对于点A的站心坐标为：
dx = 1.0000 米
dy = -1.0000 米
dz = -2.0000 米

其中，dx、dy、dz 分别表示点B相对于点A在 x、y、z 三个方向上的站心坐标。