本文主要记录一下最近的一些研究成果。
文章将分成方法部分与代码部分来写(如果有空可能会写数据可视化),为大家提供一种不同于Acrgis OD矩阵的方法来测算公共服务设施可达性。
一、概念部分
首先让我们明确一下公共服务设施、可达性与步行可达性的概念。
公共服务设施是为市民提供公共服务产品的各种公共性、服务性设施,具体包括教育、医疗卫生、文化、体育、交通、社会福利与保障、行政管理与社区服务、邮政电信和商业金融服务等设施。[1]
可达性在图论中,是指在图中从一个顶点到另一个顶点的容易程度。在无向图中,可以通过识别图的连接分量来确定所有顶点对之间的可达性。
步行可达性主要指居民从起点到终点的距离,分成时间成本与距离成本两个部分。
二、可达性的测算方法
这部分的方法主要来自《公共服务设施空间可达性及其度量方法》该文的总结 [2]
目前国内外关于可达性的测算主要分为比例法 、最近距离法 、基于机会累积的方法和基于空间相互作用的方法四个方法。
2.1 比例法
比例法主要是运用在研究范围较大的情况下,通过计算服务资源总量(设施数量、服务数量等)与服务人口总量的比值,比值越大,可达性越好。考虑到服务人口数据的获取难度(其实也不难,只要拿到人口的栅格数据然后在Acrgis里做空间匹配就可以了),同时国内也有研究表明单位面积内公共服务设施的数量越多表示设施越临近,可达性越高[3]。因此本文用用最简单的比例法及测量一定范围内该类设施的总数量来测算可达性,数量越多,可达性越好。该方法最终获取的数据为周围XX公里XX设施的总数量。
2.2 最近距离法
最近距离法主要运用在测算一些紧急设施的可达性上。顾名思义,该方法假设居民始终选择最邻近公共服务设施,距离越近,可达性越好。该方法从情理上来说是很正确的方法,毕竟周边公园再多,我也会选择时间成本与距离成本最小的那个,所以该方法最终获取的数据就一个筛选出来的点的距离与时间。
2.3 基于机会积累的方法
基于机会积累的方法主要用于计算一些非紧急设施。简而言之,就是我们通过对比在某个阈值之间公共服务设施的数量来测算公共服务设施的可达性。该方法最终获取的数据为周围XX公里接近XX阈值XX设施的总数量。(这里的阈值可以是时间也可以是距离)
基于空间交互作用的方法涉及到潜能模型上的一些内容,较为复杂,本文不做讨论(可能在后期专门写一篇文章来搞,看有没有空吧)
三、数据部分
本文数据来源主要选择高德地图API网页服务平台提供的路径规划工具
路径规划API是一套以HTTP形式提供的步行、公交、驾车查询及行驶距离计算接口,返回JSON 或 XML格式的查询数据,用于实现路径规划功能的开发。由于道路/数据/算法的变更,很可能存在间隔一段时间后请求相同起终点的经纬度返回不同结果。 [4]
本方法将主要获取该接口反馈的起点和终点的步行距离、步行时间预计两个参数。也就是我们常在使用高德地图查询路径时给我们反馈的下图两个参数:
该方法的好处其一在于得到的参数经过高德地图的测算,相对于依靠路网等矢量数据的Acrgis方法来说在准确度与数据的时效性上更为可靠;其二在于该方法在成本上相对于需要先获取POI数据与路网数据再进行OD建模的Acrgis来说更为快捷(当然,这是在你掌握了python语言的基础上)
本文选择了深圳市8个小区来分别测试其医疗设施中药店设施的可达性。
关于高德API密钥的申请,可以查阅其官方帮助文档[4]:
如何成为高德API开发者lbs.amap.com 如何申请Key?lbs.amap.com四、思路
思路上其实也很简单,主要分成三个部分。
第一个部分是获取8个小区周边的药店设施的信息,这一部分将通过高德地图API提供的圆形检索工具来实现,通过该接口,可以得到周边药店设施POI的名称、坐标、地址等属性。
第二个部分是将检索到的周边药店设施POI设置为D点,通过传入O-D,来计算两点之间的最短时间与最短距离。
第三个部分是将数据进行输出或者直接在Python中通过matpoltlib等优秀的图形库来进行可视化。
五、代码部分
精力原因,本教程将分成两部分,有空就更新第二部分。以下是整体数据采集的一个基础效果显示,不涉及到数据整理部分:
参考
- ^住房城乡建设部关于发布国家标准《城市居住区规划设计标准》的公告 http://www.mohurd.gov.cn/wjfb/201811/t20181130_238590.html
- ^宋正娜, 陈雯, 张桂香, 等. 公共服务设施空间可达性及其度量方法. 地理科学进展, 2010, 29(10): 1217-1224.
- ^ 孙斌栋, 阎宏, 张婷麟. 社区建成环境对健康的影响——基于居民个体超重的实证研究[J]. 地理学报, 2016, (10): 17-30
- ^ab高德地图API入门指南 https://lbs.amap.com/api/webservice/guide/api/direction