Python地理数据处理：Geopandas实用教程

本篇文章向大家介绍《Python处理地理数据指南：Geopandas全攻略》，主要包括，具有一定的参考价值，需要的朋友可以参考一下。

Geopandas是地理数据处理首选，因它整合了Shapely、Fiona、Matplotlib和Pandas功能于一体。1. 它基于Pandas DataFrame扩展出GeoDataFrame和GeoSeries，支持空间数据操作；2. 提供统一API简化从加载、清洗到分析、可视化的流程；3. 通过.area、.intersects()等方法实现简便空间计算；4. 支持多种格式读取如Shapefile、GeoJSON；5. 允许CRS检查与转换避免操作错误；6. 内置buffer、dissolve等空间操作方法；7. 利用sjoin进行空间连接完成区域统计；8. plot方法可直接生成可视化地图，极大提升开发效率且降低学习门槛。

Geopandas是Python处理地理空间数据的核心库，它将地理空间数据结构与Pandas的数据操作能力结合起来，让空间数据处理、分析和可视化变得前所未有的简单和高效。在我看来，它就是Python地理信息处理领域的“瑞士军刀”。

Geopandas是一个基于Pandas DataFrame的扩展，它引入了GeoDataFrame和GeoSeries这两种核心数据结构。GeoDataFrame在普通的Pandas DataFrame基础上，增加了一个特殊的“几何图形”（geometry）列，这个列存储了地理要素（如点、线、面）的几何信息。你可以像操作普通DataFrame一样对GeoDataFrame进行数据筛选、分组和聚合，同时还能直接执行复杂的空间操作，比如计算面积、判断包含关系、进行空间连接等。

为什么Geopandas是地理数据处理的首选？

说实话，在我接触Geopandas之前，Python处理地理数据是个相当零散的任务。你需要用Shapely处理几何对象，用Fiona读写文件，用Matplotlib绘图，然后用Pandas管理属性数据。这就像是把一堆零件堆在一起，你需要自己组装。Geopandas的出现，彻底改变了这种局面。它把这些核心功能巧妙地整合到一个统一的API下，让数据科学家和GIS专业人士能够在一个框架内完成从数据加载、清洗、分析到可视化的全流程。

我个人认为，Geopandas最吸引人的地方在于它的“一体化”和“直观性”。它不仅继承了Pandas的强大数据处理能力，让你可以用熟悉的方式操作地理数据，还通过简洁的API暴露了底层的空间操作。比如，你想计算一个多边形的面积，直接调用.area属性就行；想判断两个几何图形是否相交，.intersects()方法就能搞定。这种设计极大地降低了学习曲线，也提升了开发效率。以前可能需要几行Shapely代码才能完成的操作，现在Geopandas一行就能搞定，这效率提升是实实在在的。

Geopandas核心数据结构与常见操作

Geopandas主要围绕两种核心数据结构展开：GeoSeries和GeoDataFrame。

GeoSeries可以看作是Pandas Series的地理空间版本，它专门用来存储几何对象（点、线、面）。每个几何对象都有自己的空间属性和方法。

GeoDataFrame则是GeoSeries的集合，它继承了Pandas DataFrame的所有功能，并额外包含一个特殊的几何列（通常命名为'geometry'），这个列存储了GeoSeries对象。它是你进行空间数据管理和分析的主要工作台。

加载地理数据： Geopandas能够轻松读取各种常见的地理空间数据格式，比如Shapefile、GeoJSON、KML、GPX等。最常用的函数是gpd.read_file()。

import geopandas as gpd

# 加载一个Shapefile
# 假设你有一个名为 'world_countries.shp' 的国家边界Shapefile
world = gpd.read_file("world_countries.shp")
print(world.head())

# 加载GeoJSON文件
# 假设你有一个名为 'cities.geojson' 的城市点数据
cities = gpd.read_file("cities.geojson")
print(cities.head())

坐标参考系统（CRS）管理： CRS是地理数据处理中一个非常关键但又容易出错的概念。它定义了地图上点的位置如何与地球表面的真实位置对应。Geopandas允许你检查、设置和转换数据的CRS。不一致的CRS是导致空间操作失败的常见原因。

# 查看GeoDataFrame的CRS
print(world.crs)

# 转换CRS到WGS84 (EPSG:4326)
# 很多Web地图服务都使用这个CRS
world_wgs84 = world.to_crs("EPSG:4326")
print(world_wgs84.crs)

# 转换到投影坐标系，例如Web Mercator (EPSG:3857)
world_mercator = world.to_crs("EPSG:3857")
print(world_mercator.crs)

基本空间操作： Geopandas提供了丰富的空间操作方法，它们直接作用于GeoDataFrame的几何列。

• buffer(): 为几何对象创建缓冲区。

  # 为城市点创建1度（地理坐标）或1000米（投影坐标）的缓冲区
  # 注意：在地理坐标系下，单位是度；在投影坐标系下，单位是投影单位（通常是米）
  cities_buffer = cities.to_crs("EPSG:3857").buffer(1000) # 1000米缓冲区
  print(cities_buffer.head())

• dissolve(): 根据某个属性将几何图形进行合并。

  # 假设world GeoDataFrame中有一个'continent'列，按大洲合并国家边界
  continents = world.dissolve(by='continent')
  print(continents.head())

• area / length: 计算几何图形的面积或长度。

  # 计算国家面积（注意CRS单位）
  world_area = world_mercator.area / 10**6 # 转换为平方公里
  print(world_area.head())

• 空间关系判断： intersects(), contains(), within(), touches(), crosses()等。 这些方法返回布尔值，用于判断几何图形之间的空间关系。

实践案例：用Geopandas进行空间分析与可视化

实际工作中，我们经常会遇到这样的需求：手头有一堆点位数据（比如商店、客户位置），想知道它们分别属于哪个区域（比如行政区划、销售区域），并对区域内的点进行统计。Geopandas的sjoin()（空间连接）功能简直是神来之笔，它能高效地完成这项任务。

案例：统计每个国家有多少个城市

假设我们有 world (国家边界) 和 cities (城市点) 两个GeoDataFrame，它们的CRS可能不一致，我们需要先统一它们。

import geopandas as gpd
import pandas as pd
import matplotlib.pyplot as plt

# 假设你已经加载了 'world_countries.shp' 和 'cities.geojson'
# world = gpd.read_file("world_countries.shp")
# cities = gpd.read_file("cities.geojson")

# 确保CRS一致，通常选择一个通用的地理坐标系（如WGS84）或合适的投影坐标系
# 这里我们都转到WGS84
world = world.to_crs("EPSG:4326")
cities = cities.to_crs("EPSG:4326")

# 执行空间连接：将城市点连接到它们所属的国家多边形上
# op='within' 表示只连接那些在国家边界内部的城市
# how='inner' 表示只保留有匹配的行
cities_with_countries = gpd.sjoin(cities, world, how="inner", op='within')

# 统计每个国家的城市数量
# 假设world GeoDataFrame有一个'name'列表示国家名称
city_counts_per_country = cities_with_countries.groupby('name').size().reset_index(name='city_count')

# 将统计结果合并回国家GeoDataFrame，以便可视化
world_with_counts = world.merge(city_counts_per_country, on='name', how='left')

# 填充没有城市的国家为0
world_with_counts['city_count'] = world_with_counts['city_count'].fillna(0)

# 可视化结果：绘制一个分级统计图（Choropleth Map）
fig, ax = plt.subplots(1, 1, figsize=(15, 10))
world_with_counts.plot(column='city_count', cmap='OrRd', linewidth=0.8, ax=ax, edgecolor='0.8', legend=True)

ax.set_title('全球各国城市数量分布', fontsize=15)
ax.set_axis_off() # 关闭坐标轴
plt.show()

这个例子展示了Geopandas在数据整合、空间分析和可视化方面的一体化能力。从加载不同类型的数据，到统一CRS，再到执行复杂的空间连接，最后通过简单的plot()方法生成直观的地图，整个流程都非常流畅。在处理这类任务时，Geopandas无疑是我的首选工具。它的设计理念和实用性，让地理数据处理不再是GIS专业人士的专属技能，而是更多数据科学工作者可以轻松掌握的利器。

今天关于《Python地理数据处理：Geopandas实用教程》的内容介绍就到此结束，如果有什么疑问或者建议，可以在golang学习网公众号下多多回复交流；文中若有不正之处，也希望回复留言以告知！