写一段python代码, 已知某点的WGS84空间直角坐标为(X=6347491.3770,Y= -22944.8986,Z=622822.4766),编程求该点的大地坐标(B,L,H)?
时间: 2024-02-27 09:51:50 浏览: 65
PY3经典三种坐标转换.txt
可以使用 pyproj 库来进行坐标转换。以下是代码:
```python
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 米
```
其中,经度和纬度都是以度为单位的角度值,高度以米为单位。
阅读全文