【批量地点检索-高德百度地图地理编码对比】
批量地点检索
基于python地理位置描述转经纬度坐标
前提：
需要安装部署python环境
安装导入模块：
import requests
import json
import math
发送请求
调用函数
html = requests.get(url)
发送url获取html，调用函数
JasonData = json.loads（html.txt）
获取返回文本信息，字符串解码为字典
1. 百度地图开放平台 web服务api
百度地图返回百度坐标系，需要注册申请百度地图开放平台个人账户，每日额度30w条，申请api接口key（AK）:
https://lbsyun.baidu.com/apiconsole/key#/home
1.1 地理位置检索poi接口
可通过该功能，检索某一行政区划内（目前最细到城市级别）的地点信息。此方法搜索同百度地图主界面地点搜索；注：宽容度低，可能会有大量地点搜索失败，返回空值，但可强制约束行政区，可返回WGS84坐标。
url格式：
https://api.map.baidu.com/place/v2/search?query=ATM机&tag=银行®ion=北京&output=json&ak=您的ak
参数解释详见：
https://lbsyun.baidu.com/index.php?title=webapi/guide/webservice-placeapi
1.2 地理编码检索
地理编码服务提供将结构化地址数据（如：北京市海淀区上地十街十号）转换为对应坐标点（经纬度）功能，地理编码服务当前未推出国际化服务，解析地址仅限国内；注：此方法宽容度较高，地点数据库较大，与百度地图主界面检索不同，但有可能返回较大错误地址，不可强制约束行政区，不可返回WGS84坐标。
url格式：
https://api.map.baidu.com/geocoding/v3/?address=北京市海淀区上地十街10号&output=json&ak=您的ak&callback=showLocation
参数解释详见：
https://lbsyun.baidu.com/index.php?title=webapi/guide/webservice-geocoding
2. 高德地图开放平台 web服务api
高德地图使用国家火星坐标系，使用需要注册申请高德地图开放平台个人账户，每日额度30w条，申请api接口key（AK），https://lbs.amap.com/?ref=http%3A%2F%2Flbs.amap.com%2Fdev%2F
2.1 地点检索poi接口
地点搜索服务2.0是一类Web API接口服务；服务提供多种场景的地点搜索能力，包括关键字搜索、周边搜索、多边形区域搜索、ID搜索。可通过文本关键字搜索地点信息，文本可以是结构化地址，例如：北京市朝阳区望京阜荣街10号；也可以是POI名称，输出格式为火星坐标系，例如：首开广场；
url格式：（参数间使用‘&’符号连接）
https://restapi.amap.com/v5/place/text?parameters
参数详见：
https://lbs.amap.com/api/webservice/guide/api/newpoisearch
2.2 地理（逆）编码
结构化地址的定义： 首先，地址肯定是一串字符，内含国家、省份、城市、区县、城镇、乡村、街道、门牌号码、屋邨、大厦等建筑物名称。按照由大区域名称到小区域名称组合在一起的字符。一个有效的地址应该是独一无二的。注意：针对大陆、港、澳地区的地理编码转换时可以将国家信息选择性的忽略，但省、市、城镇等级别的地址构成是不能忽略的。暂时不支持返回台湾省的详细地址信息，输出不会出现大范围偏移现象，输出格式为火星坐标系。
结构化地址举例：北京市朝阳区阜通东大街6号
转换后经纬度：116.480881,39.989410
url格式：
https://restapi.amap.com/v3/geocode/geo?parameters
参数详见：
https://lbs.amap.com/api/webservice/guide/api/georegeo
3. 百度/火星坐标系转WGS84坐标系
详见作者：https://blog.csdn.net/yzyssg1/article/details/76120617
4. 参考代码
以下范例为使用百度地理编码接口批量获取地点坐标
# Editor: Jason Huang
# 输入地理名称检索WGS84坐标`
# 导入csv文件，有且仅有首列保存地理位置描述，去除表头`
import requests
import json
import math```
pi = 3.1415926535897932384626`
x_pi = 3.14159265358979324 * 3000.0 / 180.0`
a = 6378245.0 # 长半轴`
ee = 0.00669342162296594323 # 扁率`
def bd09togcj02(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 gcj02towgs84(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 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:
"""
if lng < 72.004 or lng > 137.8347:
return True
if lat < 0.8293 or lat > 55.8271:
return True
return False
def location_transfer(filename):
addressList = []
with open(filename, 'r') as file:
for item in file:
if item != '':
city,address = item.split(",",1)
addressList.append([city,address[:-1]])
ak = '你的key（账号管理中复制）'
savelocation = []
failedlocation = []
lat = []
lon = []
with open('输出坐标文件.csv','w') as f:
with open('搜索失败.csv','w') as f2:
f.write('Address Name'+','+'Latitude'+','+'Longitude'+'\n')
for i in addressList: # 设置搜索条目数量
url = 'http://api.map.baidu.com/geocoding/v3/?address='+i[1]+'&city='+ i[0]+ '&output=json&ak=' + ak
#
html = requests.get(url)
# print(html.text)
jsonData = json.loads(html.text)
if jsonData["status"] == 0:
if len(jsonData["results"]) != 0:
lat = (jsonData["results"][0]["location"]["lat"])
lon = (jsonData["results"][0]["location"]["lng"])
gg_lng, gg_lat = bd09togcj02(lon,lat)
wg_lon, wg_lat = gcj02towgs84(gg_lng, gg_lat)
item = str(i[1])+',' + str(lat)+','+str(lon)
f.write((item)+'\n')
else:
item2 = str(i[1]) + ',null,null'
f.write(item2+'\n')
f2.write(item2+'\n')
print(i[1] + ' failed to find address')
print("excute successfully")
return
if __name__ == '__main__':
"""
输入文件需仅保留地理位置名称于文件首列
"""
location_transfer('输入地理描述文件.csv')
