经纬度换算

经纬度换算是一个常见的地理坐标转换问题，通常涉及将经纬度（Lat, Lon）转换为地理坐标（如地图上的坐标、距离、方向等），或者将地理坐标转换为经纬度。

一、经纬度的定义

• 纬度（Latitude）：从赤道（0°）向北和向南延伸的角度，范围是 -90° 到 +90°。
• 经度（Longitude）：从本初子午线（0°）向东和向西延伸的角度，范围是 -180° 到 +180°。

二、经纬度换算的常见用途

1. 将经纬度转换为距离（如地球表面距离）
2. 将经纬度转换为地图上的坐标（如 Google Maps）
3. 将经纬度转换为方向（如方位角）
4. 将经纬度转换为坐标（如笛卡尔坐标系）

三、经纬度换算公式（地球半径为 6371 km）

1. 将经纬度转换为笛卡尔坐标（x, y, z）

公式如下：

$$ x = R \cdot \cos(\phi) \cdot \cos(\lambda) $$ $$ y = R \cdot \cos(\phi) \cdot \sin(\lambda) $$ $$ z = R \cdot \sin(\phi) $$

其中：

• $ R = 6371 $ km（地球平均半径）
• $ \phi $：纬度（弧度）
• $ \lambda $：经度（弧度）

2. 将笛卡尔坐标转换为经纬度

公式如下：

$$ \phi = \arctan\left( \frac{z}{\sqrt{x^2 + y^2}} \right) $$ $$ \lambda = \arctan\left( \frac{y}{x} \right) $$

四、使用工具或编程实现经纬度换算

1. Python（使用 geopy 库）

from geopy import geoposition
from geopy.distance import distance

# 将经纬度转换为地理坐标
lat, lon = 35.6895, 139.6917  # 东京
pos = geoposition.Geoposition(lat, lon)
print(pos)

# 将地理坐标转换为笛卡尔坐标
x, y, z = pos.to_cartesian()
print(f"笛卡尔坐标: ({x}, {y}, {z})")

2. 在线工具

• https://www.latlong.net/
• https://www.movable-type.co.uk/scripts/latlong.html

五、经纬度换算示例

示例 1：将经纬度转换为笛卡尔坐标

• 纬度：35.6895°N
• 经度：139.6917°E
• 地球半径：6371 km

$$ x = 6371 \cdot \cos(35.6895^\circ) \cdot \cos(139.6917^\circ) $$ $$ y = 6371 \cdot \cos(35.6895^\circ) \cdot \sin(139.6917^\circ) $$ $$ z = 6371 \cdot \sin(35.6895^\circ) $$

计算结果（近似）：

• $ x \approx 18200000 $
• $ y \approx 19500000 $
• $ z \approx 12000000 $

六、经纬度换算的注意事项

• 纬度和经度的单位是 度（°），转换为弧度时需乘以 $ \pi/180 $
• 地球是椭球体，实际换算时需考虑地球的扁率（如 WGS84）
• 精度要求：一般用于地图、导航时精度在 100 米以内即可

七、总结

如需进一步帮助，可以告诉我你具体想做什么（如：转换、计算距离、方向等），我可以为你提供更详细的指导。