经纬度换算成xy坐标

经纬度转换为平面坐标系（如笛卡尔坐标系）的XY坐标，通常使用的是地理坐标系到平面坐标系的投影方法。以下是一个常用的方法，即将经纬度转换为Web Mercator投影坐标，这是一种广泛用于在线地图服务中的投影方法。

Web Mercator投影的公式如下：

```

x = R λ

y = R ln(tan(π/4 + φ/2))

```

其中：

`R` 是地球半径，通常取平均值约为6371千米。

`λ` 是经度（以弧度为单位）。

`φ` 是纬度（以弧度为单位）。

`ln` 是自然对数。

以下是Python代码示例，将经纬度转换为Web Mercator投影的XY坐标：

```python

import math

def latlon_to_mercator(lat, lon):

R = 6378137.0 地球半径，单位：米

x = R math.radians(lon)

y = R math.log(math.tan(math.pi/4 + math.radians(lat)/2))

return (x, y)

示例使用

latitude = 39.9042 北京的纬度

longitude = 116.4074 北京的经度

xy = latlon_to_mercator(latitude, longitude)

print(f"经纬度({latitude