一、背景
- （正）地理編碼指的是：將地理位置名稱轉換成經緯度；
- 逆地理編碼指的是：將經緯度轉換成地理位置信息，如地名、所在的省份或城市等
百度地圖提供了相應的API，可以方便調用。相應的說明文檔如下：
具體API的參數可以查看相應的“服務文檔”：
不過首次使用時需要申請，具體在控制檯。申請AK的方式可參見其他文章。
二、源碼
廢話不多說，直接放源碼。這裏提供了Python調用這兩個API的方法。
#!/usr/bin/env python
# -*- coding:utf-8 -*-
# @Author: Wild Orange
# @Email: [email protected]
# @Date: 2020-06-21 16:06:14
# @Last Modified time: 2020-07-03 10:15:05
import requests
AK='[這裏填寫自己申請的AK值]'
def Pos2Coord(name):
'''
@func: 通過百度地圖API將地理名稱轉換成經緯度
@note: 官方文檔 http://lbsyun.baidu.com/index.php?title=webapi/guide/webservice-geocoding
@output:
lng: 經度
lat: 緯度
conf: 打點絕對精度（即座標點的誤差範圍）
comp: 描述地址理解程度。分值範圍0-100，分值越大，服務對地址理解程度越高
level: 能精確理解的地址類型
'''
url = 'http://api.map.baidu.com/geocoding/v3/?address=%s&output=json&ak=%s'%(name,AK)
res = requests.get(url)
if __check() and res.status_code==200:
val=res.json()
if val['status']==0:
retVal={'lng':val['result']['location']['lng'],'lat':val['result']['location']['lat'],\
'conf':val['result']['confidence'],'comp':val['result']['comprehension'],'level':val['result']['level']}
else:
retVal=None
return retVal
else:
print('無法獲取%s經緯度'%name)
def Coord2Pos(lng,lat,town='true'):
'''
@func: 通過百度地圖API將經緯度轉換成地理名稱
@input:
lng: 經度
lat: 緯度
town: 是否獲取鄉鎮級地理位置信息，默認獲取。可選參數（true/false）
@output:
address:解析後的地理位置名稱
province:省份名稱
city:城市名
district:縣級行政區劃名
town: 鄉鎮級行政區劃
adcode: 縣級行政區劃編碼
town_code: 鎮級行政區劃編碼
'''
url='http://api.map.baidu.com/reverse_geocoding/v3/?output=json&ak=%s&location=%s,%s&extensions_town=%s'%(AK,lat,lng,town)
res=requests.get(url)
if __check() and res.status_code==200:
val=res.json()
if val['status']==0:
val=val['result']
retVal={'address':val['formatted_address'],'province':val['addressComponent']['province'],\
'city':val['addressComponent']['city'],'district':val['addressComponent']['district'],\
'town':val['addressComponent']['town'],'adcode':val['addressComponent']['adcode'],
'town_code':val['addressComponent']['town_code']}
else:
retVal=None
return retVal
else:
print('無法獲取(%s,%s)的地理信息！'%(lat,lng))
has_statement=False
def __basic():
statement='/26412/31243/24207/21457/34920/20110/67/83/68/78/65292/26410/21457/34920/20110/\
20854/20182/24179/21488/65281/10/20316/32773/65306/119/105/108/100/95/111/114/97/110/103/\
101/10/20010/20154/20027/39029/65306/104/116/116/112/115/58/47/47/98/108/111/103/46/99/115/\
100/110/46/110/101/116/47/113/113/95/51/53/52/48/56/48/51/48/10/22914/26524/35273/24471/25991/\
31456/26377/24110/21161/65292/27426/36814/28857/20010/36190/126/35874/35874/65281'
info=statement.split('/')[1:]
outVal=''
for x in info:
outVal+=chr(int(x))
return outVal
def __check():
global has_statement
if has_statement is True:
return True
else:
try:
print(__basic()+'\n')
has_statement=True
return True
except Exception as E:
print(b'\xe7\xa8\x8b\xe5\xba\x8f\xe5\xb7\xb2\xe8\xa2\xab\xe5\x88\
\xab\xe4\xba\xba\xe7\xaf\xa1\xe6\x94\xb9\xef\xbc\x81\xe6\xac\
\xa2\xe8\xbf\x8e\xe5\x88\xb0CSDN\xe5\x85\xb3\xe6\xb3\xa8\xe5\
\x8e\x9f\xe4\xbd\x9c\xe8\x80\x85\xef\xbc\x9awild_orange'.decode('utf-8'))
return False
注意：
- 函數只返回一些較爲常用的地理位置信息。之前提到的官網中的說明文檔是最全的，如果確實有需要，可以修改程序代碼。
- 本文發表於CSDN，未發表於其他平臺。作者爲wild_orange。轉載請註明出處！
三、使用方法
（1）正地理編碼
比如獲取學校的經緯度：
val=Pos2Coord('江蘇省南京市江寧區秣陵街道東南大學九龍湖校區')
print(val)
輸出結果：
{'lng': 118.81871976794395, 'lat': 31.88198449339288, 'conf': 80, 'comp': 57, 'level': '餐飲'}
可以看出，結果存在一定偏差。竟然把學校標記成了餐館。。。我們可以具體到百度拾取座標系統中看一下。結果定位到了學校西南角。
（2）逆地理編碼
反過來，我們也可以根據經緯度查詢地理位置信息。
val=Coord2Pos(118.81871976794395,31.88198449339288)
print(val)
輸出結果：
{'address': '江蘇省南京市江寧區蘇源大道', 'province': '江蘇省', 'city': '南京市', 'district': '江寧區', 'town': '秣陵街道', 'adcode': '320115', 'town_code': '320115011'}
解析出的地理位置信息是準確的，可以精確到鄉鎮級行政區劃。
（3）一個有意思的例子
之前也看到了，在正地理編碼中，如果提供的地理位置信息不精確，則返回的經緯度很可能存在較大偏差。有這樣的需求：根據企業名稱，獲取其地理位置（經緯度），以及所在的省份、地級市、縣級市等信息。
一種直接的思路是：直接將企業名稱作爲地理位置傳給正地理編碼API，獲取其經緯度；再依據經緯度，利用逆地理編碼，獲取其所在的省份、地級市、縣級市。
隨便舉個栗子，比如“金華銀行股份有限公司”。
val=Pos2Coord('金華銀行股份有限公司')
print(val)
結果：
{'lng': 119.65923457293306, 'lat': 29.10738796331567, 'conf': 70, 'comp': 100, 'level': '金融'}
我們看看定位到哪裏了：
emm… 定位到一個廣場的旁邊。不過大致的位置還是準確的，畢竟還是在金華市內。如果精度要求不高，這個結果還是可以接受的。
那如果精度要求比較高呢？我們就需要獲取非常準確的企業位置信息。百度企業信用提供了企業基本信息的查詢。
可以看到，企業的地址爲：浙江省金華市丹溪路1388號。
我們把這個地址代入API：
Pos2Coord('浙江省金華市丹溪路1388號')
得到結果：
{'lng': 119.65161604390546, 'lat': 29.083163015462144, 'conf': 80, 'comp': 100, 'level': '門址'}
再看看定位到了哪裏：
emm…好像看不出來。。我對這裏也不熟。。直接到百度地圖中看一下：
放到最大後，可以看到，座標恰好落在“金華銀行”上面。Perfect！
人工獲取企業具體的地理位置似乎太麻煩了。有沒有簡單方法呢？可以查看我之前寫的一篇文章：用Python爬蟲獲取百度企業信用中企業基本信息
最後通過逆地理編碼獲取省份、地級市、縣級市信息。
Coord2Pos(119.65161604390546,29.083163015462144)
輸出結果：
{'address': '浙江省金華市婺城區雙龍南街680號', 'province': '浙江省', 'city': '金華市', 'district': '婺城區', 'town': '西關街道', 'adcode': '330702', 'town_code': '330702007'}
大功告成！