<em>Mac</em>Book项目 2009年学校开始实施<em>Mac</em>Book项目,所有师生配备一本<em>Mac</em>Book,并同步更新了校园无线网络。学校每周进行电脑技术更新,每月发送技术支持资料,极大改变了教学及学习方式。因此2011
2021-06-01 09:32:01
地理編碼與逆編碼表示的是地名地址與地理座標(經緯度)互相轉換的過程。其中,將地址資訊對映為地理座標的過程稱之為地理編碼;將地理座標轉換為地址資訊的過程稱之為逆地理編碼。(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, 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座標系間有偏差,通常需要進行轉換。
各座標間可以通過座標轉換互轉,通常有三引數或七引數法,網上可找到相關的轉換引數,貼一下轉換函數(不知道哪位大神寫的)
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)
到此這篇關於python使用百度或高德地圖獲取地理位置並轉換的文章就介紹到這了,更多相關python獲取地理位置並轉換內容請搜尋it145.com以前的文章或繼續瀏覽下面的相關文章希望大家以後多多支援it145.com!
相關文章
<em>Mac</em>Book项目 2009年学校开始实施<em>Mac</em>Book项目,所有师生配备一本<em>Mac</em>Book,并同步更新了校园无线网络。学校每周进行电脑技术更新,每月发送技术支持资料,极大改变了教学及学习方式。因此2011
2021-06-01 09:32:01
综合看Anker超能充系列的性价比很高,并且与不仅和iPhone12/苹果<em>Mac</em>Book很配,而且适合多设备充电需求的日常使用或差旅场景,不管是安卓还是Switch同样也能用得上它,希望这次分享能给准备购入充电器的小伙伴们有所
2021-06-01 09:31:42
除了L4WUDU与吴亦凡已经多次共事,成为了明面上的厂牌成员,吴亦凡还曾带领20XXCLUB全队参加2020年的一场音乐节,这也是20XXCLUB首次全员合照,王嗣尧Turbo、陈彦希Regi、<em>Mac</em> Ova Seas、林渝植等人全部出场。然而让
2021-06-01 09:31:34
目前应用IPFS的机构:1 谷歌<em>浏览器</em>支持IPFS分布式协议 2 万维网 (历史档案博物馆)数据库 3 火狐<em>浏览器</em>支持 IPFS分布式协议 4 EOS 等数字货币数据存储 5 美国国会图书馆,历史资料永久保存在 IPFS 6 加
2021-06-01 09:31:24
开拓者的车机是兼容苹果和<em>安卓</em>,虽然我不怎么用,但确实兼顾了我家人的很多需求:副驾的门板还配有解锁开关,有的时候老婆开车,下车的时候偶尔会忘记解锁,我在副驾驶可以自己开门:第二排设计很好,不仅配置了一个很大的
2021-06-01 09:30:48
不仅是<em>安卓</em>手机,苹果手机的降价力度也是前所未有了,iPhone12也“跳水价”了,发布价是6799元,如今已经跌至5308元,降价幅度超过1400元,最新定价确认了。iPhone12是苹果首款5G手机,同时也是全球首款5nm芯片的智能机,它
2021-06-01 09:30:45