查看原文
其他

geopandas,用python画地图原来这么简单!

朱卫军 Python大数据分析 2021-08-09

pandas应该是大家非常熟悉的Python第三方库,其主要用于数据整理和分析,这次来介绍pandas的一个近亲-geopandas


geopandas是用来处理地理空间数据的python第三方库,它是在pandas的基础上建立的,完美地融合了pandas的数据类型,并且提供了操作地理空间数据的高级接口,使得在python中进行GIS操作变成可能。


先看个示例,我们在python中显示世界地图

import pandas as pdimport geopandasimport matplotlib.pyplot as plt%matplotlib inline
world = geopandas.read_file(geopandas.datasets.get_path('naturalearth_lowres'))world.plot()plt.show()

当然除了上面的平面效果,你还可以做出这样的:

还可以给地图着色:



下面言归正传,正式讲解geopandas是如何工作的

前面说过,geopandas沿用了pandas的数据类型,所以geopandas中也有两种数据类型:

  • GeoSeries
  • GeoDataFrame
它们继承了pandas数据结构的大部分方法。
你可以把这两个数据结构当作地理空间数据的存储器,shapefile文件的pandas呈现。
Shapefile文件用于描述几何体对象:点,折线与多边形。例如,Shapefile文件可以存储井、河流、湖泊等空间对象的几何位置。除了几何位置,shp文件也可以存储这些空间对象的属性,例如一条河流的名字,一个城市的温度等等。
稍作解释:
  • GeoSeries对应Series,只有一列,里面的每个元素都是代表地理空间图形,有可能是点、线或者面。
  • GeoDataFrame是包含GeoSeries的数据结构,它是多列的,但其中一列必然是GeoSeries列,这个GeoSeries列被称作GeoDataFrame中的几何列。
  • GeoDataFrame的其他列,可以是几何图形的名字、属性等信息,比如国家的人口、面积、GDP等等。


解析一下前面的例子中完整代码:
import pandas as pdimport geopandasimport matplotlib.pyplot as plt%matplotlib inlineworld = geopandas.read_file(geopandas.datasets.get_path('naturalearth_lowres'))world.plot()plt.show()

对于下面一行:
world = geopandas.read_file(geopandas.datasets.get_path('naturalearth_lowres'))


geopandas的read_file方法可以读取shape文件,转化为GeoSeries和GeoDataFrame数据类型。


而geopandas.datasets.get_path('naturalearth_lowres')则是从geopandas自带的数据集中获取世界地图的shapefile文件。
你也可以用read_file方法读取自己的shapefile文件


所以说,world变量被赋予了一个GeoDataFrame数据列,它长这样:
这个数据列中,不仅有几何列geometry,还有其它属性列
比如:人口pop_est、大洲continent、国家名字name等


world.plot()则是将GeoDataFrame变成图形展示出来,得到世界地图:
假如说,我只想显示每个地区的中心点,像这样:
代码如下:
import pandas as pdimport geopandasimport matplotlib.pyplot as plt%matplotlib inline
world = geopandas.read_file(geopandas.datasets.get_path('naturalearth_lowres'))# 新增一列,每个国家的中心点world['centroid_column'] = world.centroid# 将新增列设置为几何列world = world.set_geometry('centroid_column') world.plot()plt.show()
world.centroid返回一个GeoSeries,里面包含每个地区的中心点空间数据。
.set_geometry函数则是将新增列设置为几何列,这样就会按照新的几何列显示地图。


GeoDataFrame文件的输出形式
前面讲到将shapefile读成GeoDataFrame格式,反过来你也可以将GeoDataFrame格式输出为shapefile文件。
比如前面的world:
world.to_file("countries.shp")
还可以输出为GeoJson、GeoPackage等地理空间文件类型。
world.to_file("countries.geojson", driver='GeoJSON')world.to_file("package.gpkg", layer='countries', driver="GPKG")


按照地区人口数量给世界地图着色
前面的世界地图都是一个颜色,如果我想按照每个地区人口数量给世界地图上色,该怎么办?
代码如下:
import pandas as pdimport geopandasimport matplotlib.pyplot as plt%matplotlib inline
world = geopandas.read_file(geopandas.datasets.get_path('naturalearth_lowres'))fig, ax = plt.subplots(1, 1)world.plot(column= 'pop_est', ax=ax, legend=True)plt.show()

其实简单,只需要在world.plot(column= 'pop_est', ax=ax, legend=True) 中加上分类参数column,分类标准为人口列:pop_est


你还可以改色系:
world.plot(column='gdp_per_cap', cmap='OrRd');
world.plot(column='gdp_per_cap', cmap='OrRd', scheme='quantiles');
总结:
geopandas是非常有用的地理空间处理工具,它很好的结合了pandas的功能,让你很轻松的探索空间数据,在python上快捷、批量进行gis操作。
geopandas还有许多的功能供爱好者探索,这里不做赘述,大家可以自行去官网学习。


最后展示一些官网案例:


福利时间:
推荐给大家一本比较不错的python地理分析教程:《Python地理空间分析指南》
该书内容详细讲解了如何利用python进行地理空间数据的分析。
想要PDF版本的小伙伴,在公众号后台回复:地理

END


往期精选

 


最全数据科学小抄,赶紧收藏吧!

教你如何使用Python制作酷炫二维码

numba,让python飞起来!
100行python代码爬取豆瓣《哪吒》短评

教你如何使用Python下载全网视频

小白如何入门Python爬虫

那些不为人知的优秀可视化库

Python机器学习·微教程

小白入门Python数据科学全教程

什么是matplotlib?

一文搞懂python匿名函数

一文读懂Python的map、reduce函数

一文搞懂python迭代器和生成器



Python大数据分析

data creat value

长按二维码关注

    您可能也对以下帖子感兴趣

    文章有问题?点此查看未经处理的缓存