gdal使用python翻译比例并保存为jpg

Question

我知道如何使用 gdal 翻译来缩放并通过命令行保存为 jpg：

gdal_translate image.bsq image.jpg -of JPEG -outsize 10% 10% -scale

这会产生（我称之为漂亮的图像）：

我想通过 python 生成类似的图像，例如：

from osgeo import gdal
img_bsq  = 'image.bsq'
img_jpg  = 'image.jpg'  
gdal.Translate(img_jpg, img_bsq, format='JPEG', width=1024, height=0, scaleParams=[[500,1000,10,20]])

我认为的问题是如何正确选择scaleParams。-scale命令行上似乎会自动计算值，如下所示man gdal_translate：

-scale [src_min src_max [dst_min dst_max]]:
           Rescale the input pixels values from the range src_min to src_max to the range dst_min to dst_max. If omitted the output range is 0
           to 255. If omitted the input range is automatically computed from the source data.

关于如何选择scaleParams（或其他相关选项）有什么建议吗？

Answer 1

您也应该可以将其留空，例如：

gdal.Translate(img_jpg, img_bsq, format='JPEG', width=1024, height=0, scaleParams=[[]])

这会导致 GDAL 自行猜测，如文档中所述：

-scale [src_min src_max [dst_min dst_max]]：将输入像素值从 src_min 到 src_max 范围重新缩放到 dst_min 到 dst_max 范围。如果省略，则输出范围为 0 到 255。如果省略，则根据源数据自动计算输入范围。

http://www.gdal.org/gdal_translate.html

或者，您也可以检索统计数据（每个频段）并自己编写一些内容。

获取统计数据：

ds = gdal.Open('img_bsq')
stats = [ds.GetRasterBand(i+1).GetStatistics(True, True) for i in range(ds.RasterCount)]
ds = None

vmin, vmax, vmean, vstd = zip(*stats)

有了这些统计数据，您应该能够想出一些所需的伸展运动。如果您想在每个频段的最小值和最大值之间进行缩放，您可以执行以下操作：

scaleParams = list(zip(*[vmin, vmax]))

或者，如果您想使用绝对最高和最低（所有频段）

scaleParams = [[min(vmin), max(vmax)]]

ETC。