高德地图api获取行政边界矢量方法

1.获取高德地图web服务key

2.行政区域查询API服务地址：

https://restapi.amap.com/v3/config/district?parameters

3.我们以山东省为例，看看山东行政边界

url地址

https://restapi.amap.com/v3/config/district?keywords=山东&subdistrict=0&key=你的密钥&extensions=all

subdistrict=0不返回下级行政区划，共0到3

extensions （base/all）

base:不返回行政区边界坐标点；

all:只返回当前查询district的边界值，不返回子节点的边界值；

目前不能返回乡镇/街道级别的边界值

4.返回的结果是json，我们看一下形式，使用之前博客提到的json查看器

我们要的是

polyline

行政区边界坐标点

需要注意的是

当一个行政区范围，由完全分隔两块或者多块的地块组成，每块地的 polyline 坐标串以 | 分隔 。

我的思路是，先取出polyline部分，再以|分组，接着以;分组，得到一个经纬度的数组

代码如下：我已将没条线进行了分组，以arcgis便于后续生成线和面（0,"120.903347,38.381733,0"）

import pandas as pd
import requests
r =requests.get("https://restapi.amap.com/v3/config/district?keywords=山东&subdistrict=0&key=你的密钥&extensions=all")
s=r.json()
poly=s["districts"][0]["polyline"]
p=poly.split("|")
x=[]
for i in range(len(p)):
    a=p[i].split(";")
    for j in range(len(a)):
        x.append([a[j].split(",")[0],a[j].split(",")[1],i])
c = pd.DataFrame(x)
c.to_csv('xzqh.csv')

至于如何生成shp格式，这就是arcgis的使用问题了，有问题请找arcgis帮助，哈哈

csv文件前几行

,0,1,2
0,120.900577,38.382,0
1,120.900577,38.382,0
2,120.900577,38.382,0
3,120.900577,38.382,0
4,120.900577,38.382,0
5,120.900577,38.382,0
6,120.900577,38.382,0

生成结果图

点数据

线数据