python2接口:WGS-84转BD坐标系

本文主要是介绍python2接口:WGS-84转BD坐标系，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

一、背景
关于python 解析 PL2303设备的原始数据得到WGS坐标系的经纬度的文章查看上一篇；

二、直接上干货，俩种方法，若有其他的请赐教
1.依赖百度API（耗网操作）

    def wgs84tobd09ByNetwork(self, lon, lat):#传入你的akqueryStr = '/geoconv/v1/?coords={},{}&from=1&to=5&ak=yourak'.format(lon, lat)# 对queryStr进行转码，safe内的保留字符不转换encodedStr = urllib.quote(queryStr, safe="/:=&?#+!$,;'@()*[]")# 在最后直接追加上yourskrawStr = encodedStr + 'fFyAglwnmgxnYI6c---------------------------'# md5计算出的sn值my_sn = hashlib.md5(urllib.quote_plus(rawStr)).hexdigest()url = 'http://api.map.baidu.com' + queryStr + "&sn=" + my_snres = requests.get(url)# get收到的内容json_str = res.content# print json_strdictData = json.loads(json_str)print(dictData)print (dictData["result"][0]["x"])print (dictData["result"][0]["y"])return dictData["result"][0]["x"], dictData["result"][0]["y"]

2.算法（误差极小）

    def transformlat(self,lng, lat):pi = 3.1415926535897932384626  # πret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + \0.1 * lng * lat + 0.2 * math.sqrt(math.fabs(lng))ret += (20.0 * math.sin(6.0 * lng * pi) + 20.0 *math.sin(2.0 * lng * pi)) * 2.0 / 3.0ret += (20.0 * math.sin(lat * pi) + 40.0 *math.sin(lat / 3.0 * pi)) * 2.0 / 3.0ret += (160.0 * math.sin(lat / 12.0 * pi) + 320 *math.sin(lat * pi / 30.0)) * 2.0 / 3.0return retdef transformlng(self,lng, lat):pi = 3.1415926535897932384626  # πret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + \0.1 * lng * lat + 0.1 * math.sqrt(math.fabs(lng))ret += (20.0 * math.sin(6.0 * lng * pi) + 20.0 *math.sin(2.0 * lng * pi)) * 2.0 / 3.0ret += (20.0 * math.sin(lng * pi) + 40.0 *math.sin(lng / 3.0 * pi)) * 2.0 / 3.0ret += (150.0 * math.sin(lng / 12.0 * pi) + 300.0 *math.sin(lng / 30.0 * pi)) * 2.0 / 3.0return retdef out_of_china(self,lng, lat):"""判断是否在国内，不在国内不做偏移:param lng::param lat::return:"""if lng < 72.004 or lng > 137.8347:return Trueif lat < 0.8293 or lat > 55.8271:return Truereturn Falsedef wgs84tobd09(self,lng, lat):x_pi = 3.14159265358979324 * 3000.0 / 180.0pi = 3.1415926535897932384626   # πa = 6378245.0  # 长半轴ee = 0.00669342162296594323      # 扁率if self.out_of_china(lng, lat):  # 判断是否在国内return lng, latdlat = self.transformlat(lng - 105.0, lat - 35.0)dlng = self.transformlng(lng - 105.0, lat - 35.0)radlat = lat / 180.0 * pimagic = math.sin(radlat)magic = 1 - ee * magic * magicsqrtmagic = math.sqrt(magic)dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * pi)mglat = lat + dlatmglng = lng + dlngz = math.sqrt(mglng * mglng + mglat * mglat) + 0.00002 * math.sin(mglat * x_pi)theta = math.atan2(mglat, mglng) + 0.000003 * math.cos(mglng * x_pi)bd_lng = z * math.cos(theta) + 0.0065bd_lat = z * math.sin(theta) + 0.006print('转换成百度坐标:%s,%s'%(bd_lng,bd_lat))return bd_lng,bd_lat

这篇关于python2接口:WGS-84转BD坐标系的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！

---