虹吸 NCSS returns nan 值

Question

我正在尝试使用来自 THREDDS 数据服务器的数据绘制最新的 HRRR 表面温度。

cat = TDSCatalog('https://thredds-test.unidata.ucar.edu/thredds/catalog/'
                 'grib/NCEP/HRRR/CONUS_2p5km/latest.xml')
dataset = cat.datasets[0]
ncss = dataset.subset()

sfctemp = ncss.query()
sfctemp.variables('Temperature_height_above_ground')
sfctemp.vertical_level(2.0)
sfctemp.add_lonlat().lonlat_box(north=55, south=20, east=281, west=230)
sfctemp.time(now)
sfctemp_data = ncss.get_data(sfctemp)

然而，当我从以下位置抓取并打印实际数据值时，效果很好：

sfctemp_vars = units.K * sfctemp_data.variables['Temperature_height_above_ground'][:].squeeze()

它returns:

[[nan nan nan ... nan nan nan] [nan nan nan ... nan nan nan] [nan nan nan ... nan nan nan] ... [nan nan nan ... nan nan nan] [nan nan nan ... nan nan nan] [nan nan nan ... nan nan nan]] kelvin

不完全确定我在这里做错了什么。帮助将不胜感激！

Answer 1

我认为实际上没有问题。如果我运行你的代码，我确实看到了你所看到的：

[[nan nan nan ... nan nan nan] [nan nan nan ... nan nan nan] [nan nan nan ... nan nan nan] ... [nan nan nan ... nan nan nan] [nan nan nan ... nan nan nan] [nan nan nan ... nan nan nan]] kelvin

这些只是缺失值，其中网格沿域的边缘具有 undefined/missing 值。如果我运行:

import numpy as np
np.nanmax(sfctemp_vars)

我得到：

308.1873

快速绘图：

import matplotlib.pyplot as plt
plt.imshow(sfctemp_vars)

给我一些看起来合理的东西：

注意顶部和底部的白色缺失数据块。

虹吸 NCSS returns nan 值

Siphon NCSS returns nan values

metpy

python-siphon