ArcGIS Python API可视化及分析系列教程(一):入门与简介

0
分享 2018-04-04
虾神自己热烈祝贺一下虾神开新坑……





虾神2018年第一坑,开在了可视化上面,在2016年的时候,虾神完成了一套leaflet for R的系列教程,受到了很多同学的点赞,还有同学问,虾神你能不能写写Python的啊?我们毕竟用Python多啊……所以这次专门开了Python的web可视化教程,只不过可能有一点会让有些同学觉得不爽的是,用的是ArcGIS Python API。





其实以前也写过几篇有关ArcGIS Python API的文章,大家可以访问下面的地址:



如:



ArcGIS Python API获取以及对比landsat数据






作为Esri旗下的一款强大的API工具,ArcGIS Python API也是开源的,只不过里面有不少Esri公司的小花招,让大家都得依赖在ArcGIS产品的支持下而已。





整个系列的设计如下:






嗯,正如上面写的。。。整个课程不包括Python的语法,这个是另外的内容了,





虾神一直表示:Python最强大的地方就即学即用,随便在网络上看个三五十行代码和教程就学会了,当然有同学会问:






虾神的表态就是:小学生都要学Python了,你还有啥理由说这个东西太难?










第一节 主要做个简单的介绍,首先是ArcGIS Python API是个神马东东:










从这页可以看出一些端倪:





首先是这个API主要还是依赖于ArcGIS Web系统的东西(外网依托于ArcGIS Online,内网依赖于Portal for ArcGIS),如果没有ArcGIS产品,而且没有外网的话,就没法使用了。这个不像leaflet for R,不管有没有网络,不管有没有服务器,都能用,这算是一个劣势。





其次它利用的是ArcGIS Rest API提供的web服务能力,所以实际上是一种Python调用rest服务来实现的,而交互的时候,需要利用ArcGIS JavaScript API,所以配置有时候也挺麻烦。





好吧,以上两个问题就是ArcGIS Python API最大的劣势,如果接受不了的话,就没办法了,但是如果能够接受这两个劣势,你会发现,以后的全是优势,使用起来各种心潮澎湃,特别是让虾神这种JavaScript盲拥有了Web的能力,就冲这个,视为神器,名至实归!








啥叫神器:神器就像凹凸曼同学的灯一样,不管对手有多强大,也不管凹凸曼被对手揍成了啥怂样,只要这个灯一亮,所有对手的结局都一样……





说了这么多,这个神器主要拿来干嘛呢?它能干这些事情:






获取在线(以及企业内部)地图、数据、资源
• 地理编码(全球英文,中国国内中文)
• 管理Portal 、Server& online中的企业级账户
• 空间可视化与属性可视化
• 标准和扩展的空间分析
• 在web上进行资源发布&共享
• 更新数据
• GP调用和结果可视化






虾神作为一个天天要用地图的人来说,别说有这么多功能,就冲着一个可视化,就有足够的利用说服我将这个Python工具包作为必须打call的东东了。





下面给大家演示几个功能:

首先是获取在线数据,比如虾神想获取1997年-2007年十年中5.7级以上地震数据(当然,你可能还有其他的需求,可以继续关注课程),那么这个数据已经有人发布过了,我们仅需要获取就行:








那么在你的本地实现这个功能呢?Esri官方提供了一个在线云服务使用ArcGIS Python API,使用方式如下:





1、在浏览器输入地址:

https://notebooks.esri.com/



系统会自动给你建立一个服务,如下所示:






当然,Esri的服务器在亚马逊云上面……国内因为GFW,所以连接速度可能有些(你懂的)……





2、点击New——Python 3,新建一个自己的Notebook










输入如下代码,即可完成上面那幅图:
from arcgis.gis import GIS
gis = GIS("https://www.arcgis.com", "arcgis_python", "P@ssword123")
mymap = gis.map()
mymap.basemap = "osm"
items = gis.content.search("Earthquakes*", item_type="feature Layer",outside_org=True, max_items=1)
mymap.add_layer(items[0])
mymap





最后那句是把点数据渲染成热度图。注意,不要放在一起,否则会覆盖,在jupyter里面,最好能够每行脚本都分开编写,这里是为了节省,才写成了一行,以后不再对这写法进行强调。





是不是特别简单?





然后看看另外一个功能,比如地理编码:



代码如下:
from arcgis.gis import GIS
import pandas as pd
from arcgis.geocoding import Geocoder, get_geocoders,geocode
gis = GIS("https://www.arcgis.com", "arcgis_python", "P@ssword123")
mymap = gis.map()
mymap.basemap = "osm"
mymap


results = geocode('北京市海淀区莲花池西路国家测绘地理信息局')
pd.DataFrame(results)


for res in results:
popup = {
"title" : res["attributes"]["Region"],
"content" :"address:"+res["address"]
}
mymap.draw(res,popup=popup)

是不是特别嗨皮?很有点迫切想用的感觉?如何安装这个神器呢?我们下节继续,当然等不及的同学,可以自行查阅官方页面,安装步骤进行安装:

https://developers.arcgis.com/ ... t-up/

待续未完。


文章来源:http://blog.csdn.net/allenlu2008/article/details/79025823

0 个评论

要回复文章请先登录注册