百度坐标转wgs84_【Python与GIS】用Python将火星，百度坐标转WGS84坐标

GCJ02&＃xff1a;火星坐标系&＃xff0c;是由中国国家测绘局制订的地理信息系统的坐标系统。由WGS84坐标系经加密后的坐标系&＃xff0c;加密方法并非线性加密&＃xff0c;每个地区偏移的都不一样。BD09&＃xff1a;为百度坐标系&＃xff0c;在GCJ02坐标系基础上再次加密。其中bd09ll表示百度经纬度坐标&＃xff0c;bd09mc表示百度墨卡托米制坐标WGS84&＃xff1a;国际坐标系&＃xff0c;为一种大地坐标系&＃xff0c;也是目前广泛使用的GPS全球卫星定位系统使用的坐标系。

以下demo实现火星坐标系转WGS84坐标、百度坐标系转WGS84坐标&＃xff0c;此反算适用于精度要求不高时。

import mathimport json#火星坐标系转成wgs84坐标系def gcj02towgs84(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]#百度坐标系转成wgs84坐标系def baidutowgs84(lng,lat): x &＃61; lng - 0.0065; y &＃61; lat - 0.006; z &＃61; math.sqrt(x*x &＃43; y*y) - 0.00002 *math.sin(y*PI) theta &＃61; math.atan2(y,x) - 0.000003*math.cos(x*PI) gg_lng &＃61; z*math.cos(theta) gg_lat &＃61; z*math.sin(theta) return gcj02towgs84(gg_lng,gg_lat)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(abs(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 ret;def 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(abs(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 ret;if __name__ &＃61;&＃61; &＃39;__main__&＃39;:    global a global ee global PI a &＃61; 6378245.0 ee &＃61; 0.00669342162296594323; PI &＃61; math.pi  huoxing_lon_ &＃61; 113 huoxing_lat &＃61; 32 [wgs84_lon,wgs84_lat] &＃61; gcj02towgs84(huoxing_lon_,huoxing_lat) print([wgs84_lon,wgs84_lat])