0


结合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里处理了

图片

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

图片

标签: arcgis

本文转载自: https://blog.csdn.net/weixin_45709844/article/details/120367970
版权归原作者 一点规划 所有, 如有侵权,请联系我们删除。

“结合Python的Mapbox绘制等时圈实践”的评论:

还没有评论