使用xarray绘制NDVI时序图遇Runtime Warning及空图问题求助

1. 处理NDVI计算中的无效值 #

Runtime Warning主要来自NDVI公式中nir+red=0导致的除以0操作，这类无效值（inf/NaN）会干扰均值计算和绘图：

• 计算NDVI时先过滤无效场景：

  import numpy as np
  import xarray as xr
  
  mask = (ds['nir'] + ds['red']) != 0
  ndvi = xr.where(mask, (ds['nir'] - ds['red']) / (ds['nir'] + ds['red']), np.nan)
  

• 计算均值时强制忽略NaN值：

  ndvi_mean = ndvi.mean(['x', 'y'], skipna=True)
  

2. 优化Dask分块配置 #

你的数据是Dask-backed的xarray对象，不合理的分块可能导致计算警告，即使load()/compute()后也无法消除：

• 查看当前分块结构：

  print(ndvi.data)
  

• 重新设置分块（优先按time维度分块，避免x/y维度分块过小）：

  ndvi = ndvi.chunk({'time': 1, 'x': -1, 'y': -1})  # 单时间步为一个块，x/y维度不分块
  ndvi_mean = ndvi.mean(['x', 'y'], skipna=True).compute()
  

3. 验证并修复数据有效性 #

如果均值结果全为NaN，说明原始影像数据的有效值占比极低：

• 检查均值数据是否有效：

  print(ndvi_mean.values)
  

• 转换数据类型避免数值溢出：

  ndvi = ndvi.astype('float64')
  

4. 强制指定Matplotlib绘图数据 #

若均值数据有效但仍无图形，绕过xarray的plot接口，直接用Matplotlib绘图：

import matplotlib.pyplot as plt

fig, ax = plt.subplots(figsize=(11, 4))
ax.plot(ndvi_mean['time'], ndvi_mean.values, 'b-^')
plt.show()

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