利用高德API获取整个城市的公交路线并可视化（四）

副标题：公共汽电车站点覆盖率——以厦门市公交线路为例

书接上回，我们有了公交的线路、站点数据，并同时对数据质量进行了校验，但是不同城市情况不同，需要看当地对公交交通数据的开放程度，部分城市建设的有大数据平台，也可以检索到公共交通的一些标签数据，这篇文章我们来讨论一下公交覆盖率；

公交数据获取方式参考我上篇文章：利用高德API获取整个城市的公交路线并可视化（三）-CSDN博客

首先先根据行政区把在范围内的公交站点裁剪出来，最新2024年的行政区边界参考我发布的这个数据：【免费】2024全国行政区边界（省、市、县）_莫兰指数资源-CSDN文库，本文大部分操作由arcgispro进行演示，arcgis操作步骤大致一致；

因为我们拿到是数据是高德坐标系（GCJ-02）我们需要先转成WGS84地理坐标系，批量坐标转换这里就不再赘述，移步至往期文章，再投影成WGS 1984 Web Mercator (auxiliary sphere) 投影坐标系；

接下来就是创建缓冲区，参考国家发布的公共汽电车线网评价指标标准，一般选取300m或者500m，本文以500m为例，进行计算；

方法选择平面，平面用的是欧氏距离，融合选项选择全部，在arcgis中选择ALL，目的是使得叠加的站点缓冲区平面不会被反复计算，大部分公交站都是被多条线经停的，因为在点集成线的过程中需要保值每条线都要有这个站，那生成面是时候就自然这个站点的位置生成多个空间上叠加的面域，但是平面图是看不出来的 ，所以这一步骤通过融合将该站多个面域进行合并；

当然也可以用创建泰森多边形的方法，来切分每个站点不重叠的辐射范围，并且优化展示效果；

生成的结果图层与融合缓冲区图层进行【相交】即可，需要详细步骤的可以移步到我这篇文章：

ArcGIS应用指南：生成泰森多边形的缓冲区_交通站点以缓冲区范围结合泰森多边形形成的区域共同界定站域边界-CSDN博客

接下来就是计算缓冲区面积，单位选择平方千米，坐标系选取前面的投影坐标系，有需要也可以改成自己常用的坐标系，如国家2000等；

接着我们讨论一个比较复杂的问题，就是我们如何实现多对多裁剪，为什么讨论多对多裁剪呢，因为到这里我们的行政区属性表和融合面的属性表都是多个面层，我们需要计算多个行政区内的公交覆盖率，那就需要把那些跨几个行政区的面裁剪成几部分，比如下面这张图，这个面就比较调气横跨两个行政区；

那我们如何把每个行政区对应每个融合面呢，argis似乎并没有合适的工具解决这个问题，那需要我们把问题转化，把多对多的问题，再次转化为多对一，先看我们的聚集点，就是覆盖率，所以我们把所有的面合并成一个图层，再由行政区进行裁剪，这个问题就迎刃而解了；

这里把站点500m覆盖率不重叠部分融合成一张表，再由行政区进行分割，分割要素选择''行政区''，分割字段选择''县名''，这里需要建一个文件地理数据库，也就是.gdb文件，需要把分割的文件导入这个文件地理数据库；

然后我们根据数据保存的.gdb的文件路径，打开就可以找到他们；

我们可以看到每个面都被分割开了，PS：如果需要生成融合面的时候不产生多个面层，可以把厦门市行政区合成一个完整的面就不会出现上述问题了，也自然不需要这一步了；

接下来，就是我们的看图说话环节，乍看之下，我们看到湖里区和思明区的公交覆盖率很高，因为拿不到厦门市的用地规划数据，所以我们退而求其次，结合厦门经济特区年鉴-2023里面的湖里区和思明区国土面积来做计算，数据来源：厦门经济特区年鉴-2023-框架 - Computer Generated Web Page (xm.gov.cn)

截止2022年末，湖里区和思明区的国土面积为156.09平方千米，因为指标统计的是建成区面积，所以实际值是比这个小的，但即使如此公交500m覆盖率也有75%，可见这两个区的公交站点设置率普之高，表明这两个区的公共交通系统较为发达，居民出行便利度较高；

文章仅用于分享个人学习成果与个人存档之用，分享知识，如有侵权，请联系作者进行删除。所有信息均基于作者的个人理解和经验，不代表任何官方立场或权威解读。