资讯专栏INFORMATION COLUMN

如何用Python完成百度与高德地图转换

89542767 / 479人阅读

  众所周知,Python的一个使用场景还是比较多的,在工作当中,也会涉及到多方面的一些事情。那么,今天小编写这篇文章的一个主要目的,给大家来介绍关于如何用Python完成百度与搞得地图转换,下面就给大家详细介绍下。


  一、地理编码与逆编码


  地理编码与逆编码表示的是地名地址与地理坐标(经纬度)互相转换的过程。其中,将地址信息映射为地理坐标的过程称之为地理编码;将地理坐标转换为地址信息的过程称之为逆地理编码。(ps:猜猜我在哪)

01.png

  地理编码能对一些只有名称或地址的数据进行空间化,从而可以开展空间分析和制图,对于数据分析人员十分重要。


  大量地图厂商都提供了相关的API,可以直接利用这些API进行转化。国外的如谷歌、esri、osm等,国内的有百度、高德、腾讯等公司,python提供了专业的geopy包集成了这些API调用。


  地理处理包将多带带研究,本文直接使用高德和百度的开放API进行地理编码。使用开放API前需要注册应用获取KEY


  提示:以下是本篇文章正文内容,下面案例可供参考


  二、高德地图地理编码


  详细信息可参考官方开发文档,可以看到,请求参数必填的有key和address,其他的根据需要设置

02.png

  具体代码如下:


  def getGDCor(key,addr):
  baseUrl='http://restapi.amap.com/v3/geocode/geo?'
  params={'key':key,##应用key
  'address':addr,
  'city':u'武汉'##指定城市,限制范围
  }
  url=baseUrl+urllib.parse.urlencode(params)
  req=urllib.request.Request(url)
  content=urllib.request.urlopen(req).read()
  jsonData=json.loads(content)
  lon,lat='',''
  if jsonData['status']=='1':
  try:
  corr=jsonData['geocodes'][0]['location']
  lon,lat=corr.split(',')[0],corr.split(',')[1]
  except:
  lon,lat='0','0'
  else:
  print('error')
  return(lon,lat)


  三、百度地图地理编码


  百度地图与高德类似,同样需要注册应用获取key,详细内容参考官网文档


  def getBDCor(ak,addr):
  #sleep(random.random())
  baseUrl='http://api.map.baidu.com/geocoding/v3/?'
  params={'address':addr,
  'city':u'武汉市',
  'output':'json',
  'ak':ak,##应用key
  'callback':'showLocation'
  }
  url=baseUrl+urllib.parse.urlencode(params)
  req=urllib.request.urlopen(url).read().decode()#'utf-8'针对中文需要设施编码
  cont1=req.replace("showLocation&&showLocation(",'')
  cont=cont1.replace(")",'')
  jsonData=json.loads(cont)
  lon,lat=0,0
  if jsonData['status']==0:
  try:
  lon=jsonData['result']['location']['lng']
  lat=jsonData['result']['location']['lat']
  except:
  print(addr)
  return(lon,lat)


  四、坐标转换和空间化


  高德或百度获取的地理位置都是经过加密的,高德地图位置使用的是火星坐标,百度地图是bd-09坐标,与WGS84坐标系间有偏差,通常需要进行转换。


  火星坐标(GCJ-02):国测局坐标,由WGS-84加密而成,国内必须至少使用GCJ-02坐标系,或者使用在GCJ-02加密后再进行加密的坐标系,如百度坐标系。高德和Google在国内都是使用GCJ-02坐标系,可以说,GCJ-02是国内最广泛使用的坐标系;


  百度坐标系是在GCJ-02坐标系的基础上再次加密偏移后形成的坐标系。


  各坐标间可以通过坐标转换互转,通常有三参数或七参数法,网上可找到相关的转换参数,贴一下转换函数(不知道哪位大神写的)


  import math
  x_pi=3.14159265358979324*3000.0/180.0
  pi=3.1415926535897932384626#π
  a=6378245.0#长半轴
  ee=0.00669342162296594323#扁率
  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)/(a/sqrtmagic*math.cos(radlat)*pi)
  mglat=lat+dlat
  mglng=lng+dlng
  return[mglng,mglat]
  def gcj02_to_wgs84(lng,lat):
  """
  GCJ02(火星坐标系)转GPS84
  :param lng:火星坐标系的经度
  :param lat:火星坐标系纬度
  :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)/(a/sqrtmagic*math.cos(radlat)*pi)
  mglat=lat+dlat
  mglng=lng+dlng
  return[lng*2-mglng,lat*2-mglat]
  def bd09_to_wgs84(bd_lon,bd_lat):
  lon,lat=bd09_to_gcj02(bd_lon,bd_lat)
  return gcj02_to_wgs84(lon,lat)
  def wgs84_to_bd09(lon,lat):
  lon,lat=wgs84_to_gcj02(lon,lat)
  return gcj02_to_bd09(lon,lat)
  def _transformlat(lng,lat):
  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.0
  ret+=(20.0*math.sin(lat*pi)+40.0*
  math.sin(lat/3.0*pi))*2.0/3.0
  ret+=(160.0*math.sin(lat/12.0*pi)+320*
  math.sin(lat*pi/30.0))*2.0/3.0
  return ret
  def _transformlng(lng,lat):
  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.0
  ret+=(20.0*math.sin(lng*pi)+40.0*
  math.sin(lng/3.0*pi))*2.0/3.0
  ret+=(150.0*math.sin(lng/12.0*pi)+300.0*
  math.sin(lng/30.0*pi))*2.0/3.0
  return ret
  def out_of_china(lng,lat):
  """
  判断是否在国内,不在国内不做偏移
  :param lng:
  :param lat:
  :return:
  """
  return not(lng>73.66 and lng<135.05 and lat>3.86 and lat<53.55)

  到此为止,这篇文章就给大家介绍到这里了,希望可以给大家带来更多的帮助。

文章版权归作者所有,未经允许请勿转载,若此文章存在违规行为,您可以联系管理员删除。

转载请注明本文地址:https://www.ucloud.cn/yun/128237.html

相关文章

  • 百度高德,Google地图定位偏移以及坐标系转换

    摘要:高德和在国内都是使用坐标系或在此基础上面加密为直接的坐标系,可以说,是国内最广泛使用的坐标系百度坐标系,百度坐标系是在坐标系的基础上再次加密偏移后形成的坐标系,只适用于百度地图。 概述 一:在进行地图开发过程中,我们一般能接触到以下三种类型的地图坐标系: 1.WGS-84原始坐标系,一般用国际GPS纪录仪记录下来的经纬度,通过GPS定位拿到的原始经纬度,Google和高德地图定位的的经...

    nidaye 评论0 收藏0
  • 前端小知识--地图坐标转换

    摘要:实际中我们可能会用到不同的地图,那么就对应到不同坐标系的转换,比如说,你有一份的数据服务,你要展现在百度或者高德地图上,这时候你就需要转换了。 地图坐标转换 LBS,基于位置的服务(Location Based Service),近年来已经无处不在,尤其是我们前端,相信或多或少都有接触一些地图API服务,比如高德、百度啊、谷歌啊~但是用的时候可能看到下面这些字眼:比如BD09、火星坐标...

    liangdas 评论0 收藏0
  • 高德地图坐标百度地图坐标相互转换

    摘要:来自侵删百度坐标转高德传入经度纬度高德坐标转百度传入经度纬度来自:https://blog.csdn.net/wang415229224ye/article/details/79207791 侵删//百度坐标转高德(传入经度、纬度) function bd_decrypt(bd_lng, bd_lat) { var X_PI = Math.PI * 3000.0 / 180.0; ...

    genedna 评论0 收藏0
  • gcoord: 转换WGS84、GCJ02、BD09坐标,解决百度地图高德地图坐标系不统一的问题

    摘要:做过地图相关开发的同学肯定会遇到这样一个问题同样的经纬度坐标,在百度地图和高德地图上位置不一样。解决方案百度地图以及高德地图都提供了一些方法来转换不同坐标系下的坐标,但是它们都需要进行网络请求,性能很差。 做过地图相关开发的同学肯定会遇到这样一个问题:同样的经纬度坐标,在百度地图和高德地图上位置不一样。showImg(https://segmentfault.com/img/remot...

    yimo 评论0 收藏0
  • JS中高德百度地图坐标互相转换

    //百度坐标转高德(传入经度、纬度) function bd_decrypt(bd_lng, bd_lat) { var X_PI = Math.PI * 3000.0 / 180.0; var x = bd_lng - 0.0065; var y = bd_lat - 0.006; var z = Math.sqrt(x * x + y * y) - 0.0000...

    dreamGong 评论0 收藏0

发表评论

0条评论

最新活动
阅读需要支付1元查看
<