利用Python|Arcgis实现交通等时圈|可达性
交通等时圈
第一次写,希望能帮到各位吧,先上成品图:深圳福田区各加油站小汽车出行可达性
Arcgis栅格数据准备
用福田区行政地图生成渔网点,利用arcgis中的添加XY坐标,生成渔网点的坐标(我的图层是设置为WGS-84坐标)
之后将加油站坐标打在地图上,我设置了2000m缓冲区,在与栅格点裁剪,生成每个站点2km范围内的坐标点。注意,生成渔网点坐标的目的是固定渔网的经纬度点,之后要做去重处理。
高德API获得加油站至各点的出行时间
由于我想得到比较精细的出行图,所以生成了50m*50m的渔网图,每个站点2公里周边大概4千多个点,利用api得到各个站点的出行时间。
先将84坐标转为高德坐标,包括加油站的坐标,destination是目标点的经纬度,key填入自己的密匙。获取出行时间代码如下,自己可根据需求调整。
def get_time(coordinate):
api_addr="https://restapi.amap.com/v3/direction/driving?origin="+coordinate+"&destination=113.994191,22.533795&strategy=10 &extensions = base&output=json&key="
req=requests.get(api_addr)
content=req.content
sjson=json.loads(content