GDAL Python似乎以错误的方式读取栅格。
我试图将grib2文件读取为栅格并将其稍后转换为数组。用gdal.Open（file）读取文件后，我没有收到错误，一切似乎都很好。
我正在尝试读取以m / s为单位的包含风速u分量的文件，但是加载到数组的值在-280到-260之间。
我试图找到一些相关性，以为数据可能只是以另一种格式/单位保存了。我无法执行此操作，也无法在QGIS和/或ArcGIS PRO中打开grib2文件，从而生成的光栅图像的值范围为-16至16 m / s。
>>> import gdal
>>> import numpy
>>> a = gdal.Open(file)
>>> b = a.GetRasterBand(1)
>>> c = b.ReadAsArray()
>>> print (numpy.nanmax(c))
-257.3182262420654
>>> print (numpy.nanmin(c))
-287.0833629608154
gdalinfo命令行也说：
Minimum=-287.083, Maximum=-257.318, Mean=-272.451, StdDev=4.091
但是在ArcGIS Pro和/或QGIS 3.4中打开相同的文件将显示最小值-13.6516和最大值15.0198。
GDAL Python中的grib2读者似乎有问题吗？我欢迎您为解决这个奇怪的问题提供帮助。
谢谢您的时间