如何使用rioxarray移除合并后栅格中的无数据值？

针对合并栅格后where方法无法移除-9999无数据值的问题，可以从以下几个方向排查和解决：

1. 打开单个栅格时显式指定无数据值 #

确保每个输入栅格在打开时就正确识别-9999为无数据值，避免合并过程中出现无数据标记不一致的情况：

elements = []
for item in files:
    # 打开文件时显式指定nodata参数
    ds = rxr.open_rasterio(item, nodata=-9999)
    elements.append(ds)

2. 处理浮点精度问题 #

如果栅格数据是浮点类型，直接用merged != -9999可能因精度误差导致匹配失败，改用np.isclose来判断：

import numpy as np
# 匹配接近-9999的值，设置合理的容差
merged = merged.where(~np.isclose(merged, -9999, atol=1e-3), drop=False)

3. 利用rioxarray原生无数据掩码方法 #

rioxarray提供了更直接的mask_nodata方法，可直接基于数据的nodata属性生成掩码数组，比手动where更可靠：

# 先为合并后的数据设置nodata属性
merged.rio.set_nodata(-9999, inplace=True)
# 生成掩码后的数组
merged_masked = merged.rio.mask_nodata()

4. 统一数据类型 #

检查合并后数据的类型，如果原本是整数型栅格合并后转为浮点，可能导致-9999存储为-9999.0，此时需统一比较值的类型或转换数据类型：

# 将数据转为整数类型（根据实际数据类型调整）
merged = merged.astype(np.int32)
merged = merged.where(merged != -9999, drop=False)

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