基于Python GDAL：16位RGB GeoTIFF转8位JPEG的缩放与验证方法

我来帮你搞定用Python GDAL把16位RGB GeoTIFF转成8位JPEG的这几个问题，这也是我日常处理遥感影像时经常用到的操作，咱们一个个说：

用gdal.Translate()的时候，你可以通过两个关键参数明确指定8位输出：

• outputType：直接设置为gdal.GDT_Byte（这是GDAL定义的8位无符号字节类型常量）
• format参数指定输出格式为'JPEG'

举个基础的代码片段：

from osgeo import gdal

# 输入输出文件路径
input_tif = "input_16bit_rgb.tif"
output_jpg = "output_8bit_rgb.jpg"

# 调用Translate完成格式指定
gdal.Translate(
    output_jpg,
    input_tif,
    format="JPEG",
    outputType=gdal.GDT_Byte
)

不过要注意：如果只指定输出类型不做数值缩放，16位中大于255的数值会直接被截断成255，大概率会导致影像过曝，所以接下来咱们说怎么处理数值缩放。

16位影像的数值范围通常是0-65535，而8位是0-255，需要把16位的线性范围映射到8位。用gdal.Translate()最方便的方式是通过options参数传入-scale选项，格式为-scale [原数据最小值] [原数据最大值] [目标最小值] [目标最大值]，GDAL会自动完成线性缩放计算。

比如把16位的0-65535完整映射到8位的0-255，代码可以写成：

from osgeo import gdal

input_tif = "input_16bit_rgb.tif"
output_jpg = "output_8bit_rgb.jpg"

# 定义转换选项：缩放范围+指定8位输出
translate_options = gdal.TranslateOptions(
    format="JPEG",
    outputType=gdal.GDT_Byte,
    options=["-scale", "0", "65535", "0", "255"]
)

gdal.Translate(output_jpg, input_tif, options=translate_options)

如果你的影像实际有效数值范围不是0-65535（比如有些遥感影像的有效范围是0-10000），可以先获取影像的统计值，再动态设置缩放范围，这样对比度会更自然：

from osgeo import gdal

ds = gdal.Open(input_tif)
# 获取第一个波段的统计值（RGB三波段范围一致的话可以复用，不一致则需分别处理）
band = ds.GetRasterBand(1)
min_val, max_val, _, _ = band.GetStatistics(True, True)

translate_options = gdal.TranslateOptions(
    format="JPEG",
    outputType=gdal.GDT_Byte,
    options=[f"-scale", str(min_val), str(max_val), "0", "255"]
)

gdal.Translate(output_jpg, input_tif, options=translate_options)
ds = None  # 关闭数据集释放资源

验证方法很简单，用GDAL打开输出文件后，查看每个波段的数据类型即可：

from osgeo import gdal

output_ds = gdal.Open("output_8bit_rgb.jpg")
if output_ds is None:
    print("无法打开输出文件，请检查路径是否正确")
else:
    for i in range(1, output_ds.RasterCount + 1):
        band = output_ds.GetRasterBand(i)
        data_type = band.DataType
        # GDAL常量对应：GDT_Byte=1（8位），GDT_UInt16=2（16位）
        if data_type == gdal.GDT_Byte:
            print(f"波段{i}是8位无符号字节类型")
        elif data_type == gdal.GDT_UInt16:
            print(f"波段{i}是16位无符号整数类型")
        else:
            print(f"波段{i}是其他类型：{data_type}")
    output_ds = None

你也可以用命令行工具gdalinfo快速验证，执行以下命令后，在输出结果里找DataType字段：

gdalinfo output_8bit_rgb.jpg

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