百度坐标转WGS84坐标程序

# -*- coding: utf-8 -*-
import json
import urllib
import math

x_pi = 3.14159265358979324 * 3000.0 / 180.0
pi = 3.1415926535897932384626  # π
a = 6378245.0  # 长半轴
ee = 0.00669342162296594323  # 扁率


class Geocoding:
    def __init__（self api_key）:
        self.api_key = api_key

    def geocode（self address）:
        “““
        利用高德geocoding服务解析地址获取位置坐标
        :param address:需要解析的地址
        :return:
        “““
        geocoding = {‘s‘: ‘rsv3‘
                     ‘key‘: self.api_key
                     ‘city‘: ‘全国‘
                     ‘address‘: address}
        geocoding = urllib.urlencode（geocoding）
        ret = urllib.urlopen（“%s?%s“ % （“http://restapi.amap.com/v3/geocode/geo“ geocoding））

        if ret.getcode（） == 200:
            res = ret.read（）
            json_obj = json.loads（res）
            if json_obj[‘status‘] == ‘1‘ and int（json_obj[‘count‘]） >= 1:
                geocodes = json_obj[‘geocodes‘][0]
                lng = float（geocodes.get（‘location‘）.split（‘‘）[0]）
                lat = float（geocodes.get（‘location‘）.split（‘‘）[1]）
                return [lng lat]
            else:
                return None
        else:
            return None


def gcj02_to_bd09（lng lat）:
    “““
    火星坐标系（GCJ-02）转百度坐标系（BD-09）
    谷歌、高德——>百度
    :param lng:火星坐标经度
    :param lat:火星坐标纬度
    :return:
    “““
    z = math.sqrt（lng * lng + lat * lat） + 0.00002 * math.sin（lat * x_pi）
    theta = math.atan2（lat lng） + 0.000003 * math.cos（lng * x_pi）
    bd_lng = z * math.cos（theta） + 0.0065
    bd_lat = z * math.sin（theta） + 0.006
    return [bd_lng bd_lat]


def bd09_to_gcj02（bd_lon bd_lat）:
    “““
    百度坐标系（BD-09）转火星坐标系（GCJ-02）
    百度——>谷歌、高德
    :param bd_lat:百度坐标纬度
    :param bd_lon:百度坐标经度
    :return:转换后的坐标列表形式
    “““
    x = bd_lon - 0.0065
    y = bd_lat - 0.006
    z = math.sqrt（x * x + y * y） - 0.00002 * math.sin（y * x_pi）
    theta = math.atan2（y x） - 0.000003 * math.cos（x * x_pi）
    gg_lng = z * math.cos（theta）
    gg_lat = z * math.sin（theta）
    return [gg_lng gg_lat]


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 = _transformlat（lng - 105.0 lat - 35.0）
    dlng = _transformlng（lng - 105.0 lat - 35.0）
    radlat = lat / 180.0 * pi
    magic = math.sin（radlat）
    magic = 1 - ee * magic * magic
    sqrtmagic = math.sqrt（magic）
    dlat = （dlat * 180.0） / （（a * （1 - ee）） / （magic * sqrtmagic） * pi）
    dlng = （dlng * 180.0） /