天地图JS API制作专题图

最近看了一份人民大学的报告，《中国城市政商关系排行榜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坐标系的，而天地图底图是CGCS2000坐标系的，CGC