结合奥维地图与百度地图api进行路径规划开发测试（python、HTML）

最新推荐文章于 2021-03-16 10:53:14 发布

原创

最新推荐文章于 2021-03-16 10:53:14 发布 · 3.2k 阅读

6 ·

CC 4.0 BY-SA版权

文章标签：

#python #html #经验分享

本文详细介绍如何使用Python进行地图坐标系统之间的批量转换，包括WGS84、GCJ02和BD09坐标系，并演示如何将转换后的经纬度点导入百度地图，形成路径。

第一步：奥维地图如何得到规划路径上轨迹点的经纬度（参见我之前的博客）。

第二步：各种坐标系转换，选择对应转换坐标系以及输入输出文件路径后可以直接调用（python批量转换，每次运行前清空输出文件result_ouput.txt）。源代码可参见：

https://blog.youkuaiyun.com/weixin_42032107/article/details/97497327

#!/usr/bin/env python
#-*- coding: utf-8 -*-
# @File  : ×ø±ê×ª»».py
# @Author: hukelin
# @Date  : 2019/7/18
# @Desc  :hello world!
 
'''
百度坐标转换
'''
import json
import urllib
import math
from urllib.request import urlopen, quote
import string
 
x_pi = 3.14159265358979324 * 3000.0 / 180.0
pi = 3.1415926535897932384626  # ¦Ð
a = 6378245.0  # ³¤°ëÖá
ee = 0.00669342162296594323  # Æ«ÐÄÂÊÆ½·½
coordinate = []
lng = []
lat = []
converted_lng = []
converted_lat = []
 
 
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

最低0.47元/天解锁文章