如何将CSV转成的xarray数据集转为rioxarray以使用重格网功能

核心步骤：给xarray数据集添加空间参考信息 #

你通过CSV生成的xarray数据集ds缺少空间坐标系（CRS）信息，这是rioxarray正常工作的必要条件，按以下步骤转换：

1. 导入rioxarray库

import rioxarray

2. 为数据集指定坐标系
   假设你的经纬度数据采用WGS84（全球通用GPS坐标系，EPSG代码4326），直接调用rio.set_crs方法即可将xarray数据集转为rioxarray对象：

# 替换为你实际使用的EPSG代码
ds = ds.rio.set_crs("EPSG:4326")

针对散点数据的额外处理（若适用） #

如果你的CSV数据是离散经纬度点（非规则网格），需要先插值生成规则栅格才能使用rioxarray的重格网功能：

import numpy as np

# 创建规则网格的经纬度序列（可调整点数控制分辨率）
lon_grid = np.linspace(ds.longitude.min(), ds.longitude.max(), 100)
lat_grid = np.linspace(ds.latitude.min(), ds.latitude.max(), 100)

# 将散点数据插值到规则网格
ds_grid = ds.interp(longitude=lon_grid, latitude=lat_grid)

# 为规则栅格设置坐标系
ds_grid = ds_grid.rio.set_crs("EPSG:4326")

重格网操作示例 #

转换完成后，即可使用rioxarray的重格网相关方法，比如将数据投影到另一种坐标系：

# 重投影到Web墨卡托坐标系（EPSG:3857）
ds_regridded = ds_grid.rio.reproject("EPSG:3857")

内容的提问来源于stack exchange，提问作者thefrollickingnerd