查看原文:【数据seminar】数据资源 | 干货满满!利用Python获取中国行政区边界经纬度坐标
启研学社由知名学者担任学术顾问,由高校师生与企研数据科学团队联合组建的,是以大数据资源及相关技术助力中国学术、智库与行业研究为宗旨的研究组织。团队当前的主要目标是挖掘行政、经济与社会大数据资源在经济学学术、智库与相关行业研究领域中的应用价值,以学术研究为标准开展大数据治理研究,努力探索大数据分析技术融入中国经济社会研究的可行进路。
目录
Part 2 利用Python从高德开放平台获取行政区边界经纬度
某君手中有一份某市企业的经纬度坐标数据,希望统计出该市从南至北,从东到西,每平方公里上企业数量(下文简称为:企业数量地理分布)。
处理思路如下:要求得该市的企业数量地理分布,首先,需要获取该市行政区边界经纬度坐标;然后,将企业与行政区边界的经纬度坐标换算为地理上的距离;最后,统计出该市企业数量的地理分布。
本文解决了统计某市企业数量地理分布的第一步:获取行政区边界经纬度坐标(之后会推送后续的步骤,敬请关注!)。首先介绍几种获取中国行政区边界经纬度的方法,然后基于高德开放平台[1],提供2份已实现功能的Python代码,供读者参考。
注:本文适合具有Python基础并了解HTTP协议的读者。若无相关知识储备亦无妨,本文已写好相关注释,以便读者能够将代码运行起来。
Part 1 获取中国行政区边界经纬度坐标的几种方法
本节介绍4种获取中国行政区边界经纬度坐标的方法。
方法1
从“国家基础地理信息中心”运行并维护的“全国地理信息资源目录服务系统”中下载全国基础地理数据[2]。
注:原始矢量数据,不是最终地图。该数据整体现势性为2015年,现势性,是指地图对地理信息现状的反映程度。
方法2
通过“国家基础地理信息中心”承办的“国家地理信息公共服务平台天地图2021版[3]”提供的API接口进行查询。
方法3
通过“高德开放平台”提供的API接口进行查询。本文使用该方法。
方法4
通过“百度地图开放平台[4]”提供的API接口进行查询。
上述方法优劣如下 :
获取方法 | 优势 | 劣势 |
---|---|---|
方法1:下载矢量数据 | 全量下载,获取方便 | 数据整体现势性为2015年 |
方法2:天地图平台 | 更新及时 | 需要申请权限 |
方法3:高德开放平台 | 更新及时 | 需要申请权限 |
方法4:百度地图开放平台 | 更新及时 | 需要申请权限 |
注:上述优劣比较是笔者使用中的体会,仅供参考。
Part 2 利用Python从高德开放平台获取行政区边界经纬度
本节提供了两份python代码供读者参考。第一份代码,基于高德开放平台提供的案例,利用该案例,我们无需申请KEY,亦可获取行政区边界经纬度。第二份代码,基于高德开放平台提供的API接口,需申请KEY。
读者若要运行本节的代码,请先安装好python环境和requests模块。
1、基于高德提供的模板
高德开放平台已经为广大开发人员提供了一个获取行政区边界的实用案例[5]。利用该案例我们可以非常方便地获取行政区边界的经纬度。
代码的基本思路如下:
-
python模拟浏览器发起请求
-
解析网站返回的响应数据
-
保存我们需要的数据
python代码如下:
import requests
# 这里的 AK 是高德样例中内置的 key. 不需要再申请key.
AK = '6e79f6d236e295632f21b385e363b6e8'
def get_boundary_from_gaode_demo(level: str, keyword: str) -> dict:
"""
功能: 获取某省市区的边界 gis 坐标
level: 行政级别 district 区、 city 市、 province 省
keywords: 某个具体的省, 市, 区
return: 返回网页响应.
"""
# 模拟浏览器
headers = {
'Referer': 'https://lbs.amap.com/',
'Accept': 'application/json, text/javascript, */*; q=0.01',
'User-Agent': 'Mozilla/5.0 (Windows NT 10.0; Win64; x64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/87.0.4280.88 Safari/537.36',
'Accept-Language': 'zh-CN,zh;q=0.9',
'accept-encoding': 'gzip, deflate, br'
}
# 设置参数
params = (
# 获取边界不需要返回下级行政区
('subdistrict', '0'),
# 返回行政区边界坐标组等具体信息
('extensions', 'all'),
# 行政级别 district、city、province
('level', level),
# AK, 这里是高德样例提供的 key
('key', AK),
('s', 'rsv3'),
# 返回结果为json格式
('output', 'json'),
# 检索的行政区划名称
('keywords', keyword),
('platform', 'JS'),
('logversion', '2.0'),
('appname', 'https://lbs.amap.com/demo/javascript-api/example/district-search/draw-district-boundaries'),
('sdkversion', '1.4.15'),
)
# 发起请求,获得响应.
response = requests.get('https://restapi.amap.com/v3/config/district', headers=headers, params=params)
# 返回响应结果.
res = response.json()
boundary = res['districts'][0]['polyline']
return boundary
if __name__ == '__main__':
# 参数1 : 行政级别 district 区、 city 市、 province 省
level = 'district'
# 参数2 : 待查询边界的行政区
keyword = '杭州市'
# 调用函数 get_boundary_from_gaode_demo() 获取边界经纬度
boundary = get_boundary_from_gaode_demo(level, keyword)
# 查看边界有多少个经纬度点组成
print(len(boundary))
2、基于高德开发平台的web API接口
第一份代码, 直接爬取高德提供的样例中的数据,因此不需要申请key, 但高德平台可能会更改这个key,导致代码失效。
通过查阅高德的开发文档,我们发现,高德平台提供了获取行政区边界的经纬度的web api接口文档[6],但这个方式需要向高德平台申请key。
python代码如下:
注意: 代码中需要补充填写 key 值。
import requests
def get_boundary_from_api(keywords: str) -> dict:
"""
获取某个行政区划地理信息
keywords: 具体的行政区划, 例如 '杭州市' 等
"""
# 向高德开发平台申请的key
key = 'xxxxxxxxxxxxxxxxxxxxxxx'
# 发起请求的url
url = r'https://restapi.amap.com/v3/config/district'
# 设置参数
params = (
# key参数, 是我们申请的key
("key", key),
# keywords参数, 是我们需要查询的行政区划
("keywords", keywords),
# extensions=='all' , 则响应中会返回边界的经纬度.
("extensions", 'all')
)
# 发起请求,获取响应
response = requests.get(url, params=params)
# 解析响应,返回数据
res = response.json()
boundary = res["districts"][0]["polyline"]
return boundary
if __name__ == '__main__':
# 调用函数获取数据
boundary = get_boundary_from_api('杭州市')
# 查看边界由多少经纬度点组成.
print(len(boundary))
3、小结
本节提供了从高德开放平台获取行政区边界经纬度的两份python代码, 它们各有优劣:
-
基于样例的方法,不需要向高德申请key,拿来就能用,但可能不稳定。
-
基于接口的方法,需要申请key,但更为稳定。
Part 3 总结和预告
本文分享了获取中国行政区边界经纬度坐标的四种方法,并解决了统计企业数量地理分布的第一步。之后,我们将陆续推出,经纬度坐标距离换算等相关推文,逐步实现“企业数量地理分布”的计算。敬请期待!
友情链接
[1]高德开放平台: https://lbs.amap.com/
[2]1:100万全国基础地理数据库: https://www.webmap.cn/commres.do?method=result100W
[3]国家地理信息公共服务平台: https://www.tianditu.gov.cn/
[4]百度地图开放平台: https://lbsyun.baidu.com/
[5]高德开放平台的实用案例: https://lbs.amap.com/demo/javascript-api/example/district-search/draw-district-boundaries
[6]高德web api接口文档: https://lbs.amap.com/api/webservice/guide/api/district