python3实现GPS经纬度坐标（WGS84）国测局火星坐标（GCJ02）百度坐标（BD09）相互转换

  • A+
所属分类：编程语言

在进行坐标转换前先来了解一下目前的的坐标体系分类：

一是GPS坐标，也即WGS-84坐标是一个国际的标准，一般卫星导航，原始的GPS设备中的数据都是采用这一坐标系。国外的Google地图、OSM等采用的都是这一坐标。

二是国测局坐标，国测局坐标GCJ-02坐标也叫火星坐标，是国家测绘局为了国家安全在原始坐标的基础上进行偏移得到的坐标，基本国内的电子地图、导航设备都是采用的这一坐标系，如：高德、搜搜、51地图MapABC地图，谷歌中国地图也是。

三是百度坐标，百度坐标BD-09坐标是百度公司出于商业保护在国测局坐标基础上进行的二次加密。

另外，还有一些其他公司也进行了二次加密形成了自己的坐标如：图吧坐标、搜狗坐标等。

知道了上述坐标的基本情况后下面具体怎么进行相互转换，直接上代码：

  1. # -*- coding: utf-8 -*-
  2. import json
  3. import requests
  4. import math
  6. x_pi = 3.14159265358979324 * 3000.0 / 180.0
  7. pi = 3.1415926535897932384626  # π
  8. a = 6378245.0  # 长半轴
  9. ee = 0.00669342162296594323  # 偏心率平方
  12. class Geocoding:
  13.     def __init__(self, api_key):
  14.         self.api_key = api_key
  16.     def geocode(self, address):
  17.         """
  18.         利用高德geocoding服务解析地址获取位置坐标
  19.         :param address:需要解析的地址
  20.         :return:
  21.         """
  22.         geocoding = {'s': 'rsv3',
  23.                      'key': self.api_key,
  24.                      'city': '全国',
  25.                      'address': address}
  26.         # geocoding = urllib.urlencode(geocoding)
  27.         # ret = urllib.urlopen("http://restapi.amap.com/v3/geocode/geo{}".format(geocoding))
  28.         url = "http://restapi.amap.com/v3/geocode/geo?"
  29.         ret = requests.get(url, params=geocoding)
  30.         if ret.status_code == 200:
  31.             # res = ret.json()
  32.             # json_obj = json.loads(res)
  33.             json_obj = ret.json()
  34.             if json_obj['status'] == '1' and int(json_obj['count']) >= 1:
  35.                 geocodes = json_obj['geocodes'][0]
  36.                 lng = float(geocodes.get('location').split(',')[0])
  37.                 lat = float(geocodes.get('location').split(',')[1])
  38.                 return [lng, lat]
  39.             else:
  40.                 return None
  41.         else:
  42.             return None
  45. def gcj02_to_bd09(lng, lat):
  46.     """
  47.     火星坐标系(GCJ-02)转百度坐标系(BD-09)
  48.     谷歌、高德——>百度
  49.     :param lng:火星坐标经度
  50.     :param lat:火星坐标纬度
  51.     :return:
  52.     """
  53.     z = math.sqrt(lng * lng + lat * lat) + 0.00002 * math.sin(lat * x_pi)
  54.     theta = math.atan2(lat, lng) + 0.000003 * math.cos(lng * x_pi)
  55.     bd_lng = z * math.cos(theta) + 0.0065
  56.     bd_lat = z * math.sin(theta) + 0.006
  57.     return [bd_lng, bd_lat]
  60. def bd09_to_gcj02(bd_lon, bd_lat):
  61.     """
  62.     百度坐标系(BD-09)转火星坐标系(GCJ-02)
  63.     百度——>谷歌、高德
  64.     :param bd_lat:百度坐标纬度
  65.     :param bd_lon:百度坐标经度
  66.     :return:转换后的坐标列表形式
  67.     """
  68.     x = bd_lon - 0.0065
  69.     y = bd_lat - 0.006
  70.     z = math.sqrt(x * x + y * y) - 0.00002 * math.sin(y * x_pi)
  71.     theta = math.atan2(y, x) - 0.000003 * math.cos(x * x_pi)
  72.     gg_lng = z * math.cos(theta)
  73.     gg_lat = z * math.sin(theta)
  74.     return [gg_lng, gg_lat]
  77. def wgs84_to_gcj02(lng, lat):
  78.     """
  79.     WGS84转GCJ02(火星坐标系)
  80.     :param lng:WGS84坐标系的经度
  81.     :param lat:WGS84坐标系的纬度
  82.     :return:
  83.     """
  84.     if out_of_china(lng, lat):  # 判断是否在国内
  85.         return [lng, lat]
  86.     dlat = _transformlat(lng - 105.0, lat - 35.0)
  87.     dlng = _transformlng(lng - 105.0, lat - 35.0)
  88.     radlat = lat / 180.0 * pi
  89.     magic = math.sin(radlat)
  90.     magic = 1 - ee * magic * magic
  91.     sqrtmagic = math.sqrt(magic)
  92.     dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
  93.     dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi)
  94.     mglat = lat + dlat
  95.     mglng = lng + dlng
  96.     return [mglng, mglat]
  99. def gcj02_to_wgs84(lng, lat):
  100.     """
  101.     GCJ02(火星坐标系)转GPS84
  102.     :param lng:火星坐标系的经度
  103.     :param lat:火星坐标系纬度
  104.     :return:
  105.     """
  106.     if out_of_china(lng, lat):
  107.         return [lng, lat]
  108.     dlat = _transformlat(lng - 105.0, lat - 35.0)
  109.     dlng = _transformlng(lng - 105.0, lat - 35.0)
  110.     radlat = lat / 180.0 * pi
  111.     magic = math.sin(radlat)
  112.     magic = 1 - ee * magic * magic
  113.     sqrtmagic = math.sqrt(magic)
  114.     dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
  115.     dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi)
  116.     mglat = lat + dlat
  117.     mglng = lng + dlng
  118.     return [lng * 2 - mglng, lat * 2 - mglat]
  121. def bd09_to_wgs84(bd_lon, bd_lat):
  122.     lon, lat = bd09_to_gcj02(bd_lon, bd_lat)
  123.     return gcj02_to_wgs84(lon, lat)
  126. def wgs84_to_bd09(lon, lat):
  127.     lon, lat = wgs84_to_gcj02(lon, lat)
  128.     return gcj02_to_bd09(lon, lat)
  131. def _transformlat(lng, lat):
  132.     ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + \
  133.           0.1 * lng * lat + 0.2 * math.sqrt(math.fabs(lng))
  134.     ret += (20.0 * math.sin(6.0 * lng * pi) + 20.0 *
  135.             math.sin(2.0 * lng * pi)) * 2.0 / 3.0
  136.     ret += (20.0 * math.sin(lat * pi) + 40.0 *
  137.             math.sin(lat / 3.0 * pi)) * 2.0 / 3.0
  138.     ret += (160.0 * math.sin(lat / 12.0 * pi) + 320 *
  139.             math.sin(lat * pi / 30.0)) * 2.0 / 3.0
  140.     return ret
  143. def _transformlng(lng, lat):
  144.     ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + \
  145.           0.1 * lng * lat + 0.1 * math.sqrt(math.fabs(lng))
  146.     ret += (20.0 * math.sin(6.0 * lng * pi) + 20.0 *
  147.             math.sin(2.0 * lng * pi)) * 2.0 / 3.0
  148.     ret += (20.0 * math.sin(lng * pi) + 40.0 *
  149.             math.sin(lng / 3.0 * pi)) * 2.0 / 3.0
  150.     ret += (150.0 * math.sin(lng / 12.0 * pi) + 300.0 *
  151.             math.sin(lng / 30.0 * pi)) * 2.0 / 3.0
  152.     return ret
  155. def out_of_china(lng, lat):
  156.     """
  157.     判断是否在国内，不在国内不做偏移
  158.     :param lng:
  159.     :param lat:
  160.     :return:
  161.     """
  162.     return not (lng > 73.66 and lng < 135.05 and lat > 3.86 and lat < 53.55)
  165. if __name__ == '__main__':
  166.     lng = 116.382997
  167.     lat = 39.915156
  168.     result1 = gcj02_to_bd09(lng, lat)
  169.     result2 = bd09_to_gcj02(lng, lat)
  170.     result3 = wgs84_to_gcj02(lng, lat)
  171.     result4 = gcj02_to_wgs84(lng, lat)
  172.     result5 = bd09_to_wgs84(lng, lat)
  173.     result6 = wgs84_to_bd09(lng, lat)
  175.     g = Geocoding('apikey')  # 这里填写你的高德api的key
  176.     result7 = g.geocode('广东省深圳市南山区')
  177. print(result1, result2, result3, result4, result5, result6, result7)
weinxin
我的微信公众号
爱真理，得永生！　　　　　　　　　　爱在灵灵久博客，网罗天下，福利大家！

发表评论

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen: