cartopy，一个非常好用的 Python 库！

更多资料获取

📚 个人网站：ipengtao.com

大家好，今天为大家分享一个非常好用的 Python 库 - cartopy。

Github地址：https://github.com/SciTools/cartopy

地图数据可视化在许多领域中都是至关重要的，无论是用于气象学、地理信息系统还是科学研究。Python Cartopy 库是一个功能强大的工具，用于创建地图可视化，它提供了丰富的功能和选项，可以满足各种地图可视化需求。本文将深入探讨 Python Cartopy 库，包括其基本概念、安装和使用方法，以及众多示例代码，以帮助大家更好地了解如何使用 Cartopy 创建令人印象深刻的地图可视化。

什么是 Python Cartopy？

Python Cartopy 是一个用于地图数据可视化的库，它建立在 Matplotlib 的基础之上，为用户提供了绘制各种地图投影的功能。Cartopy 的目标是使地图数据的可视化变得容易且高度可定制，同时提供了与地理数据交互的工具。

Python Cartopy 的关键特点

支持多种地图投影：Cartopy 支持多种常见的地图投影，如 Mercator、极坐标、等距圆柱等，以及用户自定义的投影。
地理特征数据：Cartopy 提供了许多地理特征数据集，包括国界、海岸线、河流、湖泊等，可用于增强地图可视化。
与 Matplotlib 集成：Cartopy 可与 Matplotlib 集成，可以在 Matplotlib 图形中添加地图元素。
地图互动性：Cartopy 允许用户轻松添加地图上的标签、注释和其他元素，以增加地图的信息含量。
地理坐标转换：Cartopy 提供了方便的方法来处理地理坐标的转换，例如将经纬度坐标转换为地图投影坐标。

安装 Python Cartopy

要开始使用 Python Cartopy，需要先安装它。

可以使用 pip 来进行安装：

pip install cartopy

安装完成后，就可以导入 Cartopy 并开始创建地图可视化了。

Python Cartopy使用

1. 基本地图绘制

首先，创建一个基本的地图可视化。

以下是一个示例代码，绘制了一个包含国界、海岸线和河流的地图：

import matplotlib.pyplot as plt
import cartopy.crs as ccrs# 创建地图投影
ax = plt.axes(projection=ccrs.PlateCarree())# 添加地理特征
ax.coastlines()
ax.add_feature(cartopy.feature.BORDERS, linestyle=':')
ax.add_feature(cartopy.feature.RIVERS)# 设置地图范围
ax.set_extent([-180, 180, -90, 90])# 显示地图
plt.show()

2. 自定义地图投影

Cartopy 可以使用不同的地图投影来满足不同需求。

以下是一个示例，绘制了一个极坐标地图：

import matplotlib.pyplot as plt
import cartopy.crs as ccrs# 创建极坐标地图投影
ax = plt.axes(projection=ccrs.NorthPolarStereo())# 添加地理特征
ax.coastlines()
ax.add_feature(cartopy.feature.BORDERS, linestyle=':')
ax.add_feature(cartopy.feature.LAND, edgecolor='black')# 设置地图范围
ax.set_extent([-180, 180, 60, 90], ccrs.PlateCarree())# 显示地图
plt.show()

3. 地图上的数据绘制

除了绘制基本地图外，Cartopy 还可以在地图上绘制数据。

以下是一个示例，绘制了一个温度数据的等温线图：

import numpy as np
import matplotlib.pyplot as plt
import cartopy.crs as ccrs# 创建地图投影
ax = plt.axes(projection=ccrs.PlateCarree())# 添加地理特征
ax.coastlines()
ax.add_feature(cartopy.feature.BORDERS, linestyle=':')# 创建示例数据
lons = np.linspace(-180, 180, 360)
lats = np.linspace(-90, 90, 180)
data = np.random.rand(180, 360)# 绘制等温线图
contour = ax.contourf(lons, lats, data, transform=ccrs.PlateCarree())# 添加颜色条
plt.colorbar(contour, ax=ax, shrink=0.6, label='Temperature (°C)')# 显示地图
plt.show()

更多功能和选项

当使用Python Cartopy库进行地图数据可视化时，有许多更高级的功能和选项可以探索，以满足各种复杂的地图可视化需求。

1. 自定义地图投影

Python Cartopy支持许多标准地图投影，但有时您可能需要使用自定义地图投影。可以使用cartopy.crs.Projection类创建自定义地图投影。

以下是一个示例，演示如何创建一个自定义地图投影并在其上绘制数据：

import matplotlib.pyplot as plt
import cartopy.crs as ccrs# 创建自定义地图投影
class MyProjection(ccrs.Projection):def __init__(self):super().__init__()@propertydef boundary(self):return (-180, 180, -90, 90)  # 自定义投影的边界坐标@propertydef threshold(self):return 1e3  # 阈值# 使用自定义投影创建地图
ax = plt.axes(projection=MyProjection())# 添加地理特征
ax.coastlines()
ax.add_feature(cartopy.feature.BORDERS, linestyle=':')# 绘制数据
# ...# 显示地图
plt.show()

2. 地图上的文本和标签

在地图上添加文本和标签可以增加地图的信息含量。Cartopy可以使用text()方法在地图上添加文本。

以下是一个示例，演示如何在地图上添加标签：

import matplotlib.pyplot as plt
import cartopy.crs as ccrs# 创建地图
ax = plt.axes(projection=ccrs.PlateCarree())# 添加地理特征
ax.coastlines()
ax.add_feature(cartopy.feature.BORDERS, linestyle=':')# 添加文本标签
ax.text(-100, 40, 'New York', transform=ccrs.PlateCarree(), fontsize=12, color='red')
ax.text(-118, 34, 'Los Angeles', transform=ccrs.PlateCarree(), fontsize=12, color='blue')# 显示地图
plt.show()

3. 与其他库的集成

Cartopy可以与其他科学计算库（如NumPy和Pandas）集成，以处理地理数据和执行更复杂的数据分析。

以下是一个示例，演示如何将Cartopy与NumPy和Matplotlib一起使用来创建等值线图：

import numpy as np
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature# 创建地图
ax = plt.axes(projection=ccrs.PlateCarree())# 添加地理特征
ax.coastlines()
ax.add_feature(cfeature.BORDERS, linestyle=':')# 创建示例数据
lons = np.linspace(-180, 180, 360)
lats = np.linspace(-90, 90, 180)
data = np.random.rand(180, 360)# 绘制等值线图
contour = ax.contourf(lons, lats, data, transform=ccrs.PlateCarree())# 添加颜色条
plt.colorbar(contour, ax=ax, shrink=0.6, label='Temperature (°C)')# 显示地图
plt.show()

4. 交互性

Cartopy还支持在地图上添加交互性元素，例如添加工具栏、缩放和平移等功能。可以使用Matplotlib的工具包来实现这些交互性。

以下是一个示例，演示如何添加缩放和平移工具：

import matplotlib.pyplot as plt
from mpl_toolkits.axes_grid1 import make_axes_locatable
import cartopy.crs as ccrs# 创建地图
ax = plt.axes(projection=ccrs.PlateCarree())# 添加地理特征
ax.coastlines()
ax.add_feature(cartopy.feature.BORDERS, linestyle=':')# 显示地图
divider = make_axes_locatable(ax)
cax = divider.append_axes("right", size="5%", pad=0.05)  # 创建颜色条轴# 添加缩放工具
plt.colorbar(contour, cax=cax, label='Temperature (°C)')
ax.set_title('Interactive Map')plt.show()