使用GDAL转换Uint16图像为Byte格式后图像偏暗的解决方法

Uint16格式像素值取值范围为0_{65535，Byte（8位无符号）像素值取值范围为0}255。默认调用gdal_translate执行位深转换时，若不额外指定缩放规则，GDAL会默认按0_65535到0255的全域线性规则做映射，也就是直接将像素值除以256得到结果。
但绝大多数Uint16影像（尤其是遥感影像、高动态范围拍摄的影像）的有效像素值并不会占满0_{65535的全量程，很多场景下有效像素区间仅在几十到几千的范围，直接做全域映射会把大部分有效像素压缩到0}255的极低值区间，最终输出图像就会出现明显偏暗的问题。

方法1：调用gdal_translate时直接指定拉伸规则（无需二次处理） #

转换时搭配-scale参数，按照影像实际有效像素区间做映射，避免无效的全域压缩：

• 先查询原始影像的实际像素值域，执行命令：
  gdalinfo -mm 输入影像路径
  从输出结果中获取影像的最小、最大像素值。如果需要更自然的视觉效果，可以取2%、98%分位值做截断，避免首尾极端异常值（比如噪点、过曝死白像素）拉伸缩放区间。
• 执行带缩放参数的转换命令：
  • 按实际最小最大值拉伸（亮度和原始影像完全匹配）：
    gdal_translate -ot Byte -scale 实际最小值 实际最大值 0 255 输入影像路径 输出影像路径
    例：若查询到影像实际有效像素范围为112~4200，命令为：
    gdal_translate -ot Byte -scale 112 4200 0 255 input.tif output.tif
  • 按2%百分比截断拉伸（GDAL 3.1及以上版本支持，视觉效果更均衡）：
    gdal_translate -ot Byte -scale 2% 98% 0 255 input.tif output.tif

  注意：如果是多波段彩色影像，建议每个波段单独统计值域、单独传入scale参数，避免全局拉伸导致的色偏。

方法2：代码调用GDAL实现可控转换 #

如果是自主编写代码调用GDAL做批量转换，可以自行实现像素映射逻辑，核心是先统计有效像素值域再做线性缩放，不要直接做值整除。Python实现示例：

from osgeo import gdal
import numpy as np

src_ds = gdal.Open("input_uint16.tif")
driver = gdal.GetDriverByName("GTiff")
out_ds = driver.CreateCopy(
    "output_byte.tif", src_ds, strict=0, options=["TILED=YES", "COMPRESS=LZW"]
)

for b in range(1, src_ds.RasterCount + 1):
    src_band = src_ds.GetRasterBand(b)
    src_arr = src_band.ReadAsArray()
    # 取2%、98%分位值作为拉伸区间，过滤极端值
    p2, p98 = np.percentile(src_arr, (2, 98))
    # 线性映射到0~255区间，超出范围的值截断到0/255
    scaled = np.clip((src_arr - p2) * 255 / (p98 - p2), 0, 255).astype(np.uint8)
    out_band = out_ds.GetRasterBand(b)
    out_band.WriteArray(scaled)
    out_band.SetColorInterpretation(src_band.GetColorInterpretation())

src_ds = None
out_ds = None

方法3：已生成的偏暗影像校正 #

如果已经完成转换得到偏暗的Byte影像，可以直接用GDAL自带的对比度拉伸工具校正，无需重新从Uint16源文件转换：
gdal_contrast_stretch -percentile 2 98 偏暗的输出影像.tif 校正后影像.tif
如果需要和原始Uint16影像视觉效果完全一致，可以提取原始影像各波段的直方图作为参考，对偏暗影像做直方图匹配即可。

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