cs2000坐标系转GCJ02的python代码
时间: 2023-10-03 20:05:54 浏览: 52
CS2000坐标系是一种大地坐标系,而GCJ02是一种转换后的坐标系。需要用到一些转换工具库才能进行转换。
以下是一个Python函数,可以将CS2000坐标系转换为GCJ02坐标系:
```python
import math
x_PI = 3.14159265358979324 * 3000.0 / 180.0
a = 6378245.0
ee = 0.00669342162296594323
def transform_lat(x, y):
ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * math.sqrt(abs(x))
ret += (20.0 * math.sin(6.0 * x * math.pi) + 20.0 * math.sin(2.0 * x * math.pi)) * 2.0 / 3.0
ret += (20.0 * math.sin(y * math.pi) + 40.0 * math.sin(y / 3.0 * math.pi)) * 2.0 / 3.0
ret += (160.0 * math.sin(y / 12.0 * math.pi) + 320 * math.sin(y * math.pi / 30.0)) * 2.0 / 3.0
return ret
def transform_lon(x, y):
ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * math.sqrt(abs(x))
ret += (20.0 * math.sin(6.0 * x * math.pi) + 20.0 * math.sin(2.0 * x * math.pi)) * 2.0 / 3.0
ret += (20.0 * math.sin(x * math.pi) + 40.0 * math.sin(x / 3.0 * math.pi)) * 2.0 / 3.0
ret += (150.0 * math.sin(x / 12.0 * math.pi) + 300.0 * math.sin(x / 30.0 * math.pi)) * 2.0 / 3.0
return ret
def transform(cs_lat, cs_lon):
dLat = transform_lat(cs_lon - 105.0, cs_lat - 35.0)
dLon = transform_lon(cs_lon - 105.0, cs_lat - 35.0)
radLat = cs_lat / 180.0 * math.pi
magic = math.sin(radLat)
magic = 1 - ee * magic * magic
sqrtMagic = math.sqrt(magic)
dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * math.pi)
dLon = (dLon * 180.0) / (a / sqrtMagic * math.cos(radLat) * math.pi)
mgLat = cs_lat + dLat
mgLon = cs_lon + dLon
return mgLat, mgLon
if __name__ == '__main__':
cs_lat = 39.908720
cs_lon = 116.397480
gcj_lat, gcj_lon = transform(cs_lat, cs_lon)
print(gcj_lat, gcj_lon)
```
使用时,只需要将`cs_lat`和`cs_lon`替换成对应的CS2000坐标系的经纬度值即可。