最近看了一份人民大学的报告,《中国城市政商关系排行榜2020》,https://new.qq.com/omn/20201230/20201230A0F3MY00.html。
讲的是中国各个城市的政商关系健康指数,决定拿这份报告中的数据,来做一份html的专题图。
效果如下(有审图号,地图不违规):
一、数据获取
做这份专题图,需要两个数据,一个是各省的边界数据,一个是各省的政商关系健康指数,后者从报告中摘取就行。
前者的参考腾讯地图web api:https://lbs.qq.com/service/webService/webServiceGuide/webServiceDistrict(因为百度没有开放获取省边界数据的接口)。
获取数据python3脚本:
import requests
import time
import shapely
# 获取省code和省name列表
def getAllProvince(key):
url = 'http://apis.map.qq.com/ws/district/v1/list?key='+key
reponse = requests.get(url=url)
reponse.encoding = 'utf-8'
data = reponse.json()
provincelist = []
for r in data['result'][0]:
provincelist.append(r['id']+'\t'+r['name'])
return provincelist
# 获取省围栏
def getProvincePolygon(key,provinceCode):
url = 'https://apis.map.qq.com/ws/district/v1/search?&keyword='+provinceCode+'&key='+key+'&get_polygon=2&max_offset=3000'
print(url)
reponse = requests.get(url=url)
reponse.encoding = 'utf-8'
data = reponse.json()
print(data)
path = data['result'][0][0]['polygon']
polygonlist = []
# 对响应结果进行差分解压,lng lat,lng lat,lng lat|lng lat……格式
for p in path:
print(p)
ringlist = []
pointnum = int(len(p)/2)
for i in range(0,pointnum):
ringlist.append(str(p[i*2])+' '+str(p[i*2+1]))
ringlist.append(ringlist[0])
polygonlist.append('POLYGON(('+','.join(ringlist)+'))')
return polygonlist
key = '你的key'
healthlist = {'北京':86.33,'上海':81.84,'天津':62.73,'海南':51.43,'浙江':49.42,'山东':48.91,
'广东':47.47,'江苏':45,'重庆':44.64,'福建':40.39,'贵州':38.4,'四川':36.74,'安徽':36.52,'广西':34.95,'江西':33.77,
'湖北':31.56,'宁夏':28.82,'湖南':28.03,'辽宁':27.7,'山西':25.73,'内蒙古':25.27,'陕西':23.95,'西藏':23.94,'甘肃':22.56,'新疆':21.99,
'青海':21.69,'河北':21.39,'吉林':21.16,'黑龙江':21.08,'河南':20.49,'云南':19.17}
f = open(r'province.txt','a',encoding='utf-8')
f.write('\t'.join(['code','name','health','polygon'])+'\n')
province_list = getAllProvince(key)
for p in province_list:
code,name = p.split('\t')
if healthlist.__contains__(name):
health = healthlist[name]
else:
health = 0.0
time.sleep(1)
polygonlist = getProvincePolygon(key,code)
for pl in polygonlist:
print(pl)
f.write('\t'.join([code,name,str(health),pl])+'\n')
f.close()
二、坐标系转换
因为省边界数据是gcj02坐标系的,而百度底图是bd09坐标系的,所以需要将上面获取到的数据转换成bd09坐标系。