高德地图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

生成结果图

点数据

线数据

互相交流GIS数据和技术,有小钱钱的指导制图请加博主微信(w912917507)

阅读更多
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页