一、地理编码与逆编码

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

地理编码能对一些只有名称或地址的数据进行空间化，从而可以开展空间分析和制图，对于数据分析人员十分重要。
大量地图厂商都提供了相关的API，可以直接利用这些API进行转化。国外的如谷歌、esri、osm等，国内的有百度、高德、腾讯等公司，python提供了专业的geopy包集成了这些API调用。
地理处理包将单独研究，本文直接使用高德和百度的开放API进行地理编码。使用开放API前需要注册应用获取KEY

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

二、高德地图地理编码

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

具体代码如下：

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, 0if 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 mathx_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.0065bd_lat = z * math.sin(theta) + 0.006return [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.0065y = bd_lat - 0.006z = 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, latdlat = _transformlat(lng - 105.0, lat - 35.0)dlng = _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 + dlngreturn [mglng, mglat]def gcj02_to_wgs84(lng, lat):"""GCJ02(火星坐标系)转GPS84:param lng:火星坐标系的经度:param lat:火星坐标系纬度:return:"""if out_of_china(lng, lat):return lng, latdlat = _transformlat(lng - 105.0, lat - 35.0)dlng = _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 + dlngreturn [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.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(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.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(lng, lat):"""判断是否在国内，不在国内不做偏移:param lng::param lat::return:"""return not (lng > 73.66 and lng < 135.05 and lat > 3.86 and lat < 53.55)

python利用百度/高德地图获取地理位置并转换相关推荐

ios 高德挪动地图获取经纬度_高德地图获取地理位置经纬度并将经维度转化为详细地址信息...
浏览器定位经纬度 116.396574, 39.992706 的地理编码结果: /*************************************** 由于Chrome.IOS10等已不再 ...
【精准三点定位求解汇总】利用Python或JavaScript高德地图开放平台实现精准三点定位（经纬度坐标与平面坐标转换法求解、几何绘图法求解）
[精准三点定位求解汇总]利用Python或JavaScript高德地图开放平台实现精准三点定位(经纬度坐标与平面坐标转换法求解.几何绘图法求解) 众所周知,如果已知三个点的坐标,到一个未知点的距离,则 ...
flutter利用高德如何获取地理位置信息bug处理
flutter中对于利用高德如何获取地理位置信息, 1.在高德开发平台,注册成为开发者,并获取key 2.安装依赖对于如何安装第三方插件请点击此处pub.dev amap_location 3.修 ...
android 高德地图 sh1,百度、高德地图获取发布版（Release）SHA1
一.简介: 在百度.高德地图开发中,申请key的时候,要两个版本的sha1值.一个是开发版(debug),一个是发布版(release). debug版本的sha1比较好获取,网上资料一堆,这里引用一 ...
百度高德地图全国poi数据
对于地图产品而言,某个地理位置周边的信息,称之为 POI .本文作者将简单地聊聊自己对于地图产品中POI的了解和看法. 作为一个标准的路痴,以前信奉的一句话是「地图长在嘴巴上」,到一个不熟悉的地方,多 ...
python代码示例百度云-python利用百度云接口实现车牌识别的示例
一个小需求---实现车牌识别. 目前有两个想法 1. 调云在线的接口或者使用SDK做开发(配置环境和编译第三方库很麻烦,当然使用python可以避免这些问题) 2. 自己实现车牌识别算法(复杂) 一开 ...
vue 获取当前位置高德_vue高德地图获取当前位置
一:安装依赖 npm install vue-amap --save 二:main.js中的配置 import VueAMap from 'vue-amap'; Vue.use(VueAMap); V ...
利用vue+高德地图API 实现用户的运动轨迹
利用vue+高德地图API 实现用户的运动轨迹高德地图网址:https://lbs.amap.com/api/jsapi-v2/guide/abc/prepare 任务一:实现地图显示先完成准备工 ...
Python利用PyQt5制作一个获取网络实时数据NBA数据播报GUI
现在NBA联赛也进行到半决赛了,我们怎么样才能以更快的方法获取NBA的数据呢?这里我们就自己来做一个数据播报的程序文章目录制作NBA数据爬虫捋顺思路编写代码 NBAReporter.py NB ...

python利用百度/高德地图获取地理位置并转换

文章目录

一、地理编码与逆编码

二、高德地图地理编码

三、百度地图地理编码

四、坐标转换和空间化

python利用百度/高德地图获取地理位置并转换相关推荐

最新文章

热门文章