Panoply与rioxarray读取NetCDF文件数值不一致问题求助

针对你遇到的Panoply与rioxarray读取NetCDF数值不一致的问题，可从以下几个方向排查：

• 坐标索引方式错误
  你直接通过pr[0,:,:].values获取数值，本质是按行列索引取数，而非按经纬度定位。NetCDF的网格点经纬度与行列号并非直接对应，容易找错位置。
  解决方法：用rioxarray的sel方法按经纬度精准定位最近的网格点：

  import rioxarray
  
  # 读取文件时指定目标变量避免歧义
  pr = rioxarray.open_rasterio("pr_aus_CESM1-CAM5_rcp45_r1i1p1_CSIRO-DecCh-wrt-1986-2005-Scl_v1_day_2016-2045.nc", variable="pr")
  # 按经纬度选择最近网格点
  target_value = pr.sel(x=141.95, y=-20.65, method='nearest').values[0]
  print(target_value)
  

• 维度顺序混淆
  rioxarray.open_rasterio读取后默认维度为(band, y, x)，而原NetCDF的pr变量维度可能是(time, lat, lon)。若你误将band维度当成时间步，但实际映射关系不符，会导致取数错误。
  解决方法：先打印数据维度确认：

  print(pr.dims)
  print(pr.coords)
  

  确认time维度是否对应band，再调整索引方式。

• rioxarray自动缩放的重复操作
  rioxarray默认会自动读取并应用scale_factor和add_offset，你代码中再次执行*scale_factor + add_offset属于重复操作（虽当前参数为1和0不影响，但需确认是否有其他隐含处理）。可通过mask_and_scale=False关闭自动缩放，手动处理验证：

  pr = rioxarray.open_rasterio("pr_aus_CESM1-CAM5_rcp45_r1i1p1_CSIRO-DecCh-wrt-1986-2005-Scl_v1_day_2016-2045.nc", mask_and_scale=False)
  a = pr[0,:,:] * pr.attrs['scale_factor'] + pr.attrs['add_offset']
  target_value = a.sel(x=141.95, y=-20.65, method='nearest').values
  print(target_value)
  

• 数据类型精度问题
  检查原数据的 dtype：若数据存储为整数类型，可能存在精度截断（虽你说scale为1、offset为0，但需确认Panoply是否做了额外浮点转换）。执行print(pr.dtype)查看数据类型，若为整数，可尝试转换为浮点后再取数：

  pr_float = pr.astype('float32')
  target_value = pr_float.sel(x=141.95, y=-20.65, method='nearest').values[0]
  

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