rioxarray栅格运算报格式错误、相加后网格改变，如何处理？

1 乘法运算报错"Wrong raster format when multiplying rasters using rioxarray" #

该报错由两个参与运算的栅格元数据不匹配触发，按以下步骤排查修复：

• 校验坐标系一致性：执行print(raster1.rio.crs, raster2.rio.crs)确认两者CRS完全一致，若不一致执行raster2 = raster2.rio.reproject(raster1.rio.crs)统一坐标系
• 校验nodata值一致性：执行print(raster1.rio.nodata, raster2.rio.nodata)确认两者nodata值相同，若不一致执行raster1 = raster1.rio.set_nodata(-9999)、raster2 = raster2.rio.set_nodata(-9999)统一设置
• 校验维度顺序：rioxarray要求栅格维度顺序为(波段, y, x)，若为单波段无波段维度的栅格，执行raster1 = raster1.expand_dims(dim="band", axis=0)补充维度后再运算

2 加法运算后网格维度改变的修复 #

该问题不是NaN值导致，本质是xarray默认的坐标对齐机制触发：两个栅格的网格分辨率、网格数量看似一致，实际x/y坐标存在浮点精度误差，运算时默认执行内连接，仅保留完全匹配的坐标，最终导致网格维度缩减。
按以下方法修复，完全保留原网格：

• 方案1：使用reproject_match强制对齐，直接将第二个栅格对齐到第一个栅格的坐标系统和网格，代码示例：

import rioxarray

# 对齐raster2到raster1的网格
raster2_aligned = raster2.rio.reproject_match(raster1)
# 执行加法运算，维度与raster1完全一致
result = raster1 + raster2_aligned

• 方案2：若确认两个栅格理论坐标完全一致，仅存在浮点精度误差，可直接替换坐标后运算：

# 把raster2的x/y坐标替换为raster1的坐标
raster2 = raster2.assign_coords(
    x = raster1.x,
    y = raster1.y
)
result = raster1 + raster2

• 额外优化：运算前可统一将nodata值替换为NaN，避免无效值干扰运算结果：

import numpy as np
raster1 = raster1.where(raster1 != raster1.rio.nodata, np.nan)
raster2 = raster2.where(raster2 != raster2.rio.nodata, np.nan)

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