结合Python的Mapbox绘制等时圈实践

emm...对于我现在工作来讲，去学代码几乎没用，上一篇用matlab还是凭以前上课的记忆东问西凑给整的，但有时候会那么一点点能省不少事情，比如这个绘制等时圈，毕竟学习精力有限，贯彻拿来主义，虽然我不懂，但能看明白且会用就可以了，也去同事那偷师了，搜集整理了网络上的信息，放在这里备忘同时也供大家查看。

对于等时圈获取我这列举三种方式

第一种，是通过ArcGIS构建网络数据集，之前也做过GIS构建等时圈

第二种，老规划云了，朴实无华却又实用，最长60min，粘贴复制坐标提交就完事了，比如我输入成都市三圣乡坐标，提交后可以获取到其15min等时圈。

http://www.guihuayun.com/dengshiquan/index.php

但是json格式的，没法在Arcgis里处理（qgis可），推荐b站up主JackLIoydSmith（在b站直播）写了一个工具，适用于规划云json文件转shp文件，可以去关注他在网盘里下载，很方便。

第三种，Mapbox

Mapbox也可以在线获取等时圈，骑行、步行和驾车三种方式，可以拖动这个点。但没有下载矢量文件的，所以获取矢量数据需要通过代码来实现。

官方公众号：时空分析必备！一文了解 Mapbox 等时圈 API

mapbox的isochrone api：https://docs.mapbox.com/playground/isochrone/ mapbox的access token申请：https://docs.mapbox.com/help/glossary/access-token/

比如还是以成都三圣乡为例。

Step1：Anaconda和库安装：需要使用的库有GDAL（需要第一个安装），Fiona，Rtree，Shapely，descartes，pyproj，geopandas（让我同事装的）

Step2：获取你需要的坐标

Step3：打开Powershell。打开文件夹，按下Shift+鼠标右键，点击在“此处打开powershell窗口”，如下图所示。

输入以下命令，回车。

python .\Get_Isochrone_from_mapbox.py --C=104.13130703288624,30.580341209379554 --P walking --M=15,30,45 --F walk --Center sanshengxiang --Fig t

ps ：1. --C=后面坐标换成自己的坐标，建议用mapbox的坐标；

      2. --P walking 换成你的出行方式，分别是walking，driving，cycling。代表步行，开车，骑车；

      3. --M=10,20,30 换成你需要的时间，比如15,25,35。注意这里的逗号是英文的逗号。

      4. --F后面是矢量等时圈面层shp保存文件的名称；

      5. --Center后面是坐标点文件shp的保存名称；

      6. --Fig后面是等时圈图片的保存名称。

最后是导出点和面两个要素，最后可以直接输出图片，图片保存在预先设置的文件夹内。 

shp文件就可以放到ArcGIS里处理了

欢迎关注微信公众号 一点规划