日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程语言 > python >内容正文

python

python获取城市区域边界坐标

發布時間:2023/12/9 python 25 豆豆
生活随笔 收集整理的這篇文章主要介紹了 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获取城市区域边界坐标的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。