wgs84坐标系转化为bj54坐标系_坐标系列篇（三）—地理坐标系纠偏转化验证

"接上篇WGS84(GPS)坐标要与国内的各种应用地图的坐标要对上&＃xff0c;就要进行纠偏&＃xff0c;本文就纠偏的一些代码(来自CSDN)做个验证&＃xff0c;为后续做软件插件提供质量保证。"

01

—

三种坐标

三种坐标分别为谷歌(WGS84)、火星坐标、百度坐标。

02

—

偏移代码

导入模块跟设置全局参数&＃xff1a;

import mathx_pi &＃61; 3.14159265358979324 * 3000.0 / 180.0pi &＃61; 3.1415926535897932384626 # πa &＃61; 6378245.0 # 长半轴ee &＃61; 0.00669342162296594323 # 偏心率平方

火星坐标系转百度坐标系&＃xff1a;

def gcj02_to_bd09(lng, lat): """ 火星坐标系(GCJ-02)转百度坐标系(BD-09) 谷歌、高德——>百度 :param lng:火星坐标经度 :param lat:火星坐标纬度 :return: """ z &＃61; math.sqrt(lng * lng &＃43; lat * lat) &＃43; 0.00002 * math.sin(lat * x_pi) theta &＃61; math.atan2(lat, lng) &＃43; 0.000003 * math.cos(lng * x_pi) bd_lng &＃61; z * math.cos(theta) &＃43; 0.0065 bd_lat &＃61; z * math.sin(theta) &＃43; 0.006 return [bd_lng, bd_lat]

百度坐标系转火星坐标系&＃xff1a;

def bd09_to_wgs84(bd_lon, bd_lat): lon, lat &＃61; bd09_to_gcj02(bd_lon, bd_lat) return gcj02_to_wgs84(lon, lat)

谷歌坐标转火星坐标系&＃xff1a;

def wgs84_to_gcj02(lng, lat): """ WGS84转GCJ02(火星坐标系) :param lng:WGS84坐标系的经度 :param lat:WGS84坐标系的纬度 :return: """ if out_of_china(lng, lat): # 判断是否在国内 return [lng, lat] dlat &＃61; _transformlat(lng - 105.0, lat - 35.0) dlng &＃61; _transformlng(lng - 105.0, lat - 35.0) radlat &＃61; lat / 180.0 * pi magic &＃61; math.sin(radlat) magic &＃61; 1 - ee * magic * magic sqrtmagic &＃61; math.sqrt(magic) dlat &＃61; (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi) dlng &＃61; (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi) mglat &＃61; lat &＃43; dlat mglng &＃61; lng &＃43; dlng return [mglng, mglat]

火星坐标系转谷歌坐标系&＃xff1a;

def gcj02_to_wgs84(lng, lat): """ GCJ02(火星坐标系)转GPS84 :param lng:火星坐标系的经度 :param lat:火星坐标系纬度 :return: """ if out_of_china(lng, lat): return [lng, lat] dlat &＃61; _transformlat(lng - 105.0, lat - 35.0) dlng &＃61; _transformlng(lng - 105.0, lat - 35.0) radlat &＃61; lat / 180.0 * pi magic &＃61; math.sin(radlat) magic &＃61; 1 - ee * magic * magic sqrtmagic &＃61; math.sqrt(magic) dlat &＃61; (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi) dlng &＃61; (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi) mglat &＃61; lat &＃43; dlat mglng &＃61; lng &＃43; dlng return [lng * 2 - mglng, lat * 2 - mglat]

百度坐标系转谷歌坐标系&＃xff1a;

def bd09_to_wgs84(bd_lon, bd_lat): lon, lat &＃61; bd09_to_gcj02(bd_lon, bd_lat) return gcj02_to_wgs84(lon, lat)

谷歌坐标系转百度坐标系&＃xff1a;

def wgs84_to_bd09(lon, lat): lon, lat &＃61; wgs84_to_gcj02(lon, lat) return gcj02_to_bd09(lon, lat)

用到的一些方法&＃xff1a;

def _transformlat(lng, lat): ret &＃61; -100.0 &＃43; 2.0 * lng &＃43; 3.0 * lat &＃43; 0.2 * lat * lat &＃43; \ 0.1 * lng * lat &＃43; 0.2 * math.sqrt(math.fabs(lng)) ret &＃43;&＃61; (20.0 * math.sin(6.0 * lng * pi) &＃43; 20.0 * math.sin(2.0 * lng * pi)) * 2.0 / 3.0 ret &＃43;&＃61; (20.0 * math.sin(lat * pi) &＃43; 40.0 * math.sin(lat / 3.0 * pi)) * 2.0 / 3.0 ret &＃43;&＃61; (160.0 * math.sin(lat / 12.0 * pi) &＃43; 320 * math.sin(lat * pi / 30.0)) * 2.0 / 3.0 return retdef _transformlng(lng, lat): ret &＃61; 300.0 &＃43; lng &＃43; 2.0 * lat &＃43; 0.1 * lng * lng &＃43; \ 0.1 * lng * lat &＃43; 0.1 * math.sqrt(math.fabs(lng)) ret &＃43;&＃61; (20.0 * math.sin(6.0 * lng * pi) &＃43; 20.0 * math.sin(2.0 * lng * pi)) * 2.0 / 3.0 ret &＃43;&＃61; (20.0 * math.sin(lng * pi) &＃43; 40.0 * math.sin(lng / 3.0 * pi)) * 2.0 / 3.0 ret &＃43;&＃61; (150.0 * math.sin(lng / 12.0 * pi) &＃43; 300.0 * math.sin(lng / 30.0 * pi)) * 2.0 / 3.0 return retdef out_of_china(lng, lat): """ 判断是否在国内&＃xff0c;不在国内不做偏移 :param lng: :param lat: :return: """ return not (lng > 73.66 and lng <135.05 and lat > 3.86 and lat <53.55)

03

—

验证测试

找一个点测试偏移后的值在三个坐标&＃xff0c;这里选用的坐标点是长沙橘子洲头雕塑这个地标&＃xff0c;分别在googleearth、高德地图、百度地图进行测试&＃xff1a;

步骤一&＃xff1a;在谷歌地球应用下找到橘子洲头雕塑地标&＃xff0c;并复制经纬度&＃xff0c;十进制[112.954951,28.171073]

步骤二&＃xff1a;调用前面的方法&＃xff0c;wsg84_to_gcj02&＃xff0c;得到偏移后的火星坐标经纬度[112.960436,28.167598]&＃xff1b;用高德地图开放接口的坐标拾缀找到该位置&＃xff1a;

步骤三&＃xff1a;调用前面的方法&＃xff0c;wsg84_to_bd09&＃xff0c;得到偏移后的百度坐标经纬度[112.966961,28.173417]&＃xff1b;用百度地图开放接口的坐标拾缀找到该位置&＃xff1a;

单一的坐标转换代码没问题后&＃xff0c;后续可以写一个支持txt、excel、csv格式的纠偏小程序exe。

后台回复关键词&＃xff1a;教育&＃xff0c;可查看我的个人在线课程&＃xff0c;欢迎大家来沟通、交流、共同学习&＃xff01;