Python地信专题 | 基于geopandas的空间数据分析-文件IO篇-linux运维-PHP中文网

Python地信专题 | 基于geopandas的空间数据分析-文件IO篇

蓮花仙者

发布： 2025-07-15 08:52:14

原创

328人浏览过

1 简介

在上一篇文章中我们对geopandas中的坐标参考系有了较为深入的学习，而在日常空间数据分析工作中矢量文件的读入和写出，是至关重要的环节。

作为基于geopandas的空间数据分析系列文章的第三篇，通过本文你将会学习到geopandas中的文件IO。

2 文件IO 2.1 矢量文件的读入

geopandas将fiona作为操纵矢量数据读写功能的后端。

使用geopandas.read_file()读取对应类型文件，而在后端实际上是使用fiona.open来读入数据，即两者参数是保持一致的，读入的数据自动转换为GeoDataFrame。

下面是geopandas.read_file()主要参数：

立即学习“Python免费学习笔记（深入）”；

下面结合上述参数，来介绍一下使用geopandas.read_file()在不同情况下读取常见格式矢量数据的方法。

使用到的示例数据为中国地图，CRS为EPSG:4326。

本文使用到的所有数据都可以在文章开头提及的Github仓库对应本文路径下找到：

图1

2.1.1 shapefile

作为非常常见的一种矢量文件格式，geopandas对shapefile提供了很好的读取和写出支持。

下面分为不同情况来介绍：

完整的shapefile

如图2，这是一个完整的shapefile：

图2

使用geopandas来读取这种形式的shapefile很简单：

代码语言：javascript代码运行次数：0运行复制

<code class="javascript">import geopandas as gpddata = gpd.read_file('geometry/china_provinces/china_provinces.shp')print(data.crs) # 查看数据对应的crsdata.head() # 查看前5行</code>

登录后复制

图3

缺少投影的shapefile

当shapefile中缺失.prj文件时，使用geopandas读入后形成的GeoDataFrame会缺失crs属性：

图4

如果已经知道数据对应的CRS，可以在读入数据后补充上crs信息以进行其他操作：

代码语言：javascript代码运行次数：0运行复制

<code class="javascript">import pyprojdata.crs = pyproj.CRS.from_user_input('EPSG:4326')data.crs</code>

登录后复制

图5

直接读取文件夹

当文件夹下只有单个shapefile时，可以直接读取该文件夹：

图6

读取zip压缩包中的文件

geopandas通过传入特定语法格式的文件路径信息，以支持直接读取.zip格式压缩包中的shapefile文件，主要分为两种情况。

当文件在压缩包内的根目录时，使用下面的语法规则来读取数据：

代码语言：javascript代码运行次数：0运行复制

<code class="javascript">zip://路径/xxx.zip</code>

登录后复制

譬如我们要读取图7所示的压缩包内文件：

图7

按照对应的语法规则，读取该类型数据方式如下：

图8

而当文件在压缩包内的文件夹中时，如图9：

图9

使用下面的语法规则来读取数据：

代码语言：javascript代码运行次数：0运行复制

<code class="javascript">zip://路径/xxx.zip!压缩包内指定文件路径</code>

登录后复制

将上述语法运用到上述文件：

图10

2.1.2 gdb与gpkg

对于Arcgis中的地理数据库gdb，以及QGIS中的GeoPackage，要读取其包含的矢量数据，就要涉及到图层的概念。

对应geopandas.read_file()的layer参数，只需要将gdb或gpkg文件路径作为filename参数，再将对应的图层名称作为layer参数传入：

gdb代码语言：javascript代码运行次数：0运行复制

<code class="javascript">data = gpd.read_file('geometry/china_provinces.gdb',                     layer='china_provinces')print(data.crs) # 查看数据对应的crsdata.head() # 查看前5行</code>

登录后复制

图11

gpkg

类似读入gdb文件：

硅基智能

基于Web3.0的元宇宙，去中心化的互联网，高质量、沉浸式元宇宙直播平台，用数字化重新定义直播

查看详情

代码语言：javascript代码运行次数：0运行复制

<code class="javascript">data = gpd.read_file('geometry/china_provinces.gpkg',                     layer='china_provinces',                     encoding='utf-8')print(data.crs) # 查看数据对应的crsdata.head() # 查看前5行</code>

登录后复制

图12

2.1.3 GeoJSON

作为web地图中最常使用的矢量数据格式，GeoJSON几乎被所有在线地图框架作为数据源格式，在geopandas中读取GeoJSON非常简单，只需要传入文件路径名称即可。

下面我们来读入图13所示的文件：

图13

图14

2.1.4 过滤

geopandas在0.1.0版本中新增了bbox过滤，在0.7.0版本中新增了蒙版过滤和行过滤功能，可以辅助我们根据自己的需要读入原始数据中的子集。

下面一一进行介绍：

bbox过滤

bbox过滤允许我们在read_file()中传入一个边界框作为参数bbox，格式为(左下角x, 左下角y, 右上角x, 右上角y)，这样在读入的过程中只会保留几何对象与bbox有相交的数据记录。

下面我们仍然以上文中使用过的中国地图数据为例，我们在读入的过程中，传入边界框：

代码语言：javascript代码运行次数：0运行复制

<code class="javascript">from shapely import geometrydata = gpd.read_file('geometry/china_provinces.json',                    bbox=(100, 20, 110, 30))%matplotlib widgetax = data.plot()# 绘制bbox框示意ax = gpd.GeoSeries([geometry.box(minx=100,                                 miny=20,                                 maxx=110,                                 maxy=30).boundary]).plot(ax=ax, color='red')</code>

登录后复制

图15

可以看到只有跟红色框有相交的几何对象被读入。

蒙版过滤

蒙版过滤和bbox过滤功能相似，都是筛选与指定区域相交的数据记录。

不同的是蒙版过滤通过mask参数可以传入任意形状的多边形，不再像bbox过滤那样只接受矩形：

代码语言：javascript代码运行次数：0运行复制

<code class="javascript">data = gpd.read_file('geometry/china_provinces.json',                    mask=geometry.Polygon([(100, 20), (110, 30), (120, 20)]))ax = data.plot()# 绘制bbox框示意ax = gpd.GeoSeries([geometry.Polygon([(100, 20),                                      (110, 30),                                      (120, 20)]).boundary]).plot(ax=ax, color='red')</code>

登录后复制

图16

可以看到只有跟红色多边形相交的几何对象被读入。

行过滤

行过滤的功能就比较简单，通过参数rows控制读入原数据的前若干行，可以用于在读取大型数据时先快速查看前几行以了解整个数据的格式：

图17

2.2 矢量文件的写出

在geopandas中使用to_file()来将GeoDataFrame或GeoSeries写出为矢量文件，主要支持shapefile、GeoJSON以及GeoPackage。

不像geopandas.read_file()可以根据传入的文件名称信息自动推断类型，我们在写出矢量数据时就需要使用driver参数来声明文件类型：

ESRI Shapefile

我们将上文最后一次读入的GeoDataFrame写出为ESRI Shapefile，设置driver参数为ESRI Shapefile。

如果你对文件编码有要求，这里可以使用encoding参数来指定，譬如这里我们指定为utf-8：

代码语言：javascript代码运行次数：0运行复制

<code class="javascript">'''在工程根目录下创建output文件夹'''import ostry:    os.mkdir('output')except FileExistsError:    pass    data.to_file('output/output.shp',             driver='ESRI Shapefile',             encoding='utf-8')</code>

登录后复制

可以看到在output文件夹下，成功导出了完整的shapefile：

图18

而如果导出的文件名不加后缀扩展名，则会生成包含在新目录下的shapefile：

代码语言：javascript代码运行次数：0运行复制

<code class="javascript">data.to_file('output/output_shapefile',             driver='ESRI Shapefile',             encoding='utf-8')</code>

登录后复制

图19

也可以向指定的文件夹下追加图层：

代码语言：javascript代码运行次数：0运行复制

<code class="javascript">data.to_file('output/output_shapefile_multi_layer',             driver='ESRI Shapefile',             layer='layer1',             encoding='utf-8')data.to_file('output/output_shapefile_multi_layer',             driver='ESRI Shapefile',             layer='layer2',             encoding='utf-8')data.to_file('output/output_shapefile_multi_layer',             driver='ESRI Shapefile',             layer='layer3',             encoding='utf-8')</code>

登录后复制