如何将CF规范坐标系转换为Cartopy的CRS对象

Cartopy本身并没有内置直接从NetCDF的CF兼容元数据生成CRS对象的方法，不过你可以通过以下两种实用方式实现：

方法一：用xarray快速解析（推荐） #

xarray对NetCDF的CF元数据支持非常好，能自动识别grid_mapping相关信息：

1. 首先用xarray打开NetCDF文件：

import xarray as xr
import cartopy.crs as ccrs

# 替换为你的NetCDF文件路径
ds = xr.open_dataset('your_data.nc')

2. 提取CF规范中的grid_mapping变量（通常命名为crs、projection或对应投影名称，比如lambert_conformal）：

# 从数据集属性中获取grid_mapping变量名
grid_mapping_name = ds.attrs.get('grid_mapping', 'crs')
grid_mapping = ds[grid_mapping_name]

3. 根据元数据类型生成CRS对象：

• 如果元数据包含EPSG代码：

crs = ccrs.epsg(grid_mapping.epsg)

• 如果包含Proj4字符串：

crs = ccrs.Proj(grid_mapping.proj4)

• 如果包含WKT格式的坐标信息：

crs = ccrs.CRS.from_wkt(grid_mapping.wkt)

方法二：手动解析CF元数据（无依赖xarray时） #

如果不用xarray，你可以用netCDF4库读取元数据，然后根据CF规范的投影参数手动构建Cartopy CRS：
比如针对Lambert Conformal投影：

from netCDF4 import Dataset
import cartopy.crs as ccrs

nc = Dataset('your_data.nc', 'r')
# 获取grid_mapping变量
proj_var = nc.variables['lambert_conformal']

# 提取CF规范定义的投影参数
crs = ccrs.LambertConformal(
    central_longitude=proj_var.longitude_of_central_meridian,
    standard_parallels=(proj_var.standard_parallel, proj_var.second_standard_parallel),
    central_latitude=proj_var.latitude_of_projection_origin
)

不同投影类型对应不同的Cartopy CRS类，你需要根据CF元数据中的grid_mapping_name属性（比如lambert_conformal、mercator等）选择对应的类，并传入对应的参数。

注意点 #

• Cartopy的API确实没有专门在“输入输出能力”章节提供直接从CF元数据生成CRS的方法，所以需要借助xarray或手动解析参数来实现。
• 生成的CRS对象可以直接用作Cartopy绘图函数（比如ax.plot()、ax.pcolormesh()）的transform参数。

内容的提问来源于stack exchange，提问作者Stéphane