在GDAL中制作的GeoTIFF在Arc / QGIS中加载时没有最小/最大范围

我正在使用GeoTIFF驱动程序在GDAL中创建一些浮点栅格。当我将生成的图像加载到QGIS或Arc中时，默认的符号系统将最小值设置为-3.40282e38，将最大值设置为3.40282e38，因此栅格看起来完全是灰色的。有没有一种方法可以直接将实际范围写入GeoTIFF，以便在将其加载到GIS程序中时自动很好地缩放直方图？我尝试过创建默认的直方图，如下所示：

rasterMin, rasterMax = raster.GetRasterBand(1).ComputeRasterMinMax()
raster.GetRasterBand(1).SetDefaultHistogram(rasterMin, rasterMax, 255)

但是我收到一个看起来像是绑定问题的错误：

  File "/usr/lib/python2.6/dist-packages/osgeo/gdal.py", line 846, in SetDefaultHistogram
    return _gdal.Band_SetDefaultHistogram(self, *args)
TypeError: not a sequence

关于我该怎么办，或者如果我做错了什么建议？

如果要生成数据，则可以随时跟踪最小值和最大值。如果您无法直接控制源数据，例如说要合成现有数据，则可以使用ComputeRasterMinMax(0)。

无论哪种方式，在获得最小值和最大值之后，都应调用SetStatistics()该带，并确保将数据集设置为None完成后才能确保将所有数据都写入磁盘。

另外，dataset.GetRasterBand(1).<some band function>由于GDAL Python处理波段引用的方式，我脑海中隐隐地告诉我使用将无法正常工作。我不记得详细信息，但是前段时间GDALDev邮件列表中包含一些内容。因此，最好显式创建一个band变量：band = dataset.GetRasterBand(1)，然后调用band.SetStatistics(min, max, mean, stddev)。

我不知道使用GDAL的方法，但是您可以使用python调用Arc中的“计算统计信息”工具。