从数据获取到交通可达性分析【全流程】

嗨咯各位小伙伴们，小编又来跟大家分享GIS干货啦！

今天呢，我们一起探讨一下ArcGIS中基于最小阻抗的交通可达性的完整流程，当然这个教程在其他许多公众号里都比较常见，所以今天小编会从数据获取、数据处理、构建交通网络、计算交通可达性等几个方面出发，将交通可达性分析整个流程从无到有的跟大家一一落实。好啦，接下来就跟小编一起进入正题吧！

我们知道，交通可达性分析已经比较成熟应用于城市问题研究，

首先分享一个下载路网数据的网站www.openstreetmap.org。

1.路网数据获取。

路网的获取方式有很多种，比如全国基础地理信息库、各种下载数据的软件、osm网站等。那我们就以OSM路网为例，先看看怎么下载吧！首先，打开osm官网，搜索我们想要的地区，比如长沙，页面会自动跳转到长沙范围。

然后点击导出，选择手动框选区域，我们以五一广场周边作为研究区域（注意：这种下载方式不适合过大范围的选择，如需大范围的下载可选择下载镜像文件，在左侧边栏下面有几种方式可以选择），点击导出即可，数据会以osm的格式保存在本地。

2.ArcGIS加载osm文件。

从osm官网下载下来的文件是osm格式的，ArcGIS不支持怎么办？这时我们需要用到一个插件——ArcGIS editor for OSM（插件的下载链接我会附在文章后面，安装也很简单，一搜就有）。这里我们先用这个工具将osm数据加载进来，设置好目标文件名（放在默认地理数据库就行），系统会自动帮我们以shp格式保存下来，此过程比较慢，不是卡，需要耐心等待一会儿。

  可以看见，点线面都被加载进来啦，接下来我们需要把道路提取出来用于构建交通网络，直接按属性选择，然后根据我们分析目标选择需要的道路，右键图层导出数据即可。

3.构建交通路网

前几步主要介绍了怎样从osm网站下载数据并在ArcGIS中加载出来，进而提取出我们需要的道路网。那么这一步就是正式进入交通可达性分析的流程啦，首先我们先进行构建交通网络。

a.基础数据导入。

在工作文件夹中新建一个个人地理数据库“交通网络”，在“交通网络”下新建一个要素数据集“路网”，将之前整理好的道路导入到“路网”中（一定要记得先投影）。

b.基础数据编辑及检查

将数据库中的道路加载进来，打开编辑器，开始编辑，右键“道路”打开属性表，按属性选择道路类型为主要道路，下拉编辑器菜单选择合并，按照相同的步骤将其他同类型的道路一一合并。

全部合并后，将所有要素选中，打开高级编辑工具，点击打断相交线将要素在交点处打断，这是构建交通网络的需要，然后在编辑器中点击保存并停止编辑。

c.拓扑检查

在“路网”要素数据集中新建一个拓扑，按照提示一直往下走，直到设置拓扑规则的步骤，按照图中所示添加规则，设置完成后验证拓扑。

出现错误的地方，用“修剪”、“延伸”等工具进行修正，在验证拓扑直到没有错误即可构建交通路网。最终结果如图所示，路网中间已经没有拓扑错误，只有道路尽头存在悬挂点。

d.设置道路属性

右键道路属性，添加一个新字段“Drivetime”代表车行时间（分钟），打开道路属性表。利用字段计算器计算车行时间，公式为Drivetime=shape_length/1000。这里表示车行速度60公里每小时下的车行时间。当然每种类型的道路的车行时间不同，大家可以尝试。本文就假定是60公里每小时吧~

e.新建交通网络

右键“路网”数据集，新建网络数据集，按照提示往下点，依次设置路口转弯、连通性、高程建模、通行成本等属性，注意的是为网络制定通行成本的时候选择Drivetime为默认字段。即可完成一个简单的交通网络的构建。构建好的交通网络包括节点和边两种要素。

4.交通可达性计算

这里的交通可达性主要是区域各位置至其他任意位置的交通便捷程度，这里我们利用最小阻抗的可达性分析。我们继续吧~

a.新建OD成本矩阵

首先启动网络分析工具（Network Analyst）,下拉工具条菜单选择新建OD成本矩阵，加载起始点和目的地点，在这里我们默认将所有路口的交点设置为起始点和目的地点，设置好后点击工具条上的求解按钮。

b.计算可达性

打开od线的属性表可以看见，起始点、目的地点和车行时间是我们想看见的，接下来需要计算起始点的可达性。

右键点击“OriginID”,选择汇总，勾选“total_minutes”的总和选项，意思就是按照“OriginID”分类汇总车行时间，汇总方法是求和。最终生成可达性计算表。

将生成的可达性计算表添加到地图中显示，打开表格，添加一个新的字段“可达性”，利用字段计算器输入公式可达性=[Sum_Total_] /( [Cnt_Origin] -1)，即可求解成功。

c.可视化

将可达性表连接到起始点上，右键起始点，选择连接，按照图中提示进行设置即可。

右键起始点，选择属性，切换到符号系统中，对可达性进行分级显示，选择一个自己喜欢的色带，即可。

在工具中找到反距离权重插值工具，按照图中提示进行参数设置，其他的暂时先默认就行啦，点击确定生成可达性空间分布图。

   最后再通过一系列地图整饬进行制图表达就完成啦，最后给大家看看效果图~~~

    好啦，今天的内容到这里就结束啦，如果小伙伴们在实验过程中有什么地方遇到问题可以加小编提问，小编将竭力为您解答，我们下期见哦，后台回复【osm】领取插件链接点个赞呗！