python获取城市区域边界坐标
生活随笔
收集整理的這篇文章主要介紹了
python获取城市区域边界坐标
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
# -*- coding:utf-8 -*-
# 第一行必須有,否則報中文字符非ascii碼錯誤
import urllib.request
from urllib.parse import quote
import numpy as np
import json
import pandas as pd
from pandas import Series, DataFrame#TODO
# 高德上申請的key
key = '9f99fc570ccaf6abc209780433d9f4c1'
#搜素的城市名(全名)
addr_name = '福州'url = 'http://restapi.amap.com/v3/config/district?'
def getlnglat(address):uri = url + 'keywords=' + quote(address) + '&key=' + key + '&subdistrict=1' + '&extensions=all'# 訪問鏈接后,api會回傳給一個json格式的數據temp = urllib.request.urlopen(uri)temp = json.loads(temp.read())# polyline是坐標,name是區域的名字Data = temp["districts"][0]['polyline']name = temp["districts"][0]['name']# polyline數據是一整個純文本數據,不同的地理塊按照|分,塊里面的地理信息按照;分,橫縱坐標按照,分,因此要對文本進行三次處理Data_Div1 = Data.split('|') # 對結果進行第一次切割,按照|符號len_Div1 = len(Data_Div1) # 求得第一次切割長度num = 0len_Div2 = 0 # 求得第二次切割長度,也即整個數據的總長度while num < len_Div1:len_Div2 += len(Data_Div1[num].split(';'))num += 1num = 0num_base = 0output = np.zeros((len_Div2, 5)).astype(np.float) # 循環2次,分割;與,while num < len_Div1:temp = Data_Div1[num].split(';')len_temp = len(temp)num_temp = 0while num_temp < len_temp:output[num_temp + num_base, :2] = np.array(temp[num_temp].split(',')) # 得到橫縱坐標output[num_temp + num_base, 2] = num_temp + 1 # 得到橫縱坐標的連接順序output[num_temp + num_base, 3] = num + 1 # 得到塊的序號num_temp += 1num_base += len_tempnum += 1output = DataFrame(output, columns=['經度', '緯度', '連接順序', '塊', '名稱'])output['名稱'] = namereturn outputdef getSubName(address): # 獲取搜索區域的名稱,部分區域例如鼓樓重名太多,因此返回城市代碼,將城市代碼作為參數給上述函數uri = url + 'keywords=' + quote(address) + '&key=' + key + '&subdistrict=1' + '&extensions=all'temp = urllib.request.urlopen(uri)temp = json.loads(temp.read())list0 = temp['districts'][0]['districts']num_Qu = 0output = []while num_Qu < len(list0):output.append(list0[num_Qu]['adcode'])num_Qu += 1return outputnum = 0
ad = getSubName(addr_name) # 得到福州下屬區域的城市代碼
add = getlnglat(addr_name) # 得到福州整個的邊界數據
while num < len(ad):add = pd.concat([add, getlnglat(ad[num].encode("utf-8"))]) # 得到福州下屬的全部區域的邊界數據num += 1
add.to_csv('{0}.csv'.format(addr_name), encoding='gbk')
需要在TODO下面的地方修改高德地圖申請的用于web服務的密鑰key,addr_name修改為需要的城市名,代碼執行后,會在當前目錄下生成cvs文件,里面的即為區域邊界坐標數據。
參考文章忘了,如有冒犯,請見諒。
總結
以上是生活随笔為你收集整理的python获取城市区域边界坐标的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 在类中用class时数据是共有还是私有_
- 下一篇: python编程软件排行榜_PYPL 9