为何Sentinel-2波段数值约为25000？如何正确可视化L1C级jp2遥感数据

一、Sentinel-2 L1C波段数值偏高的原因 #

你读取到的超过10000甚至20000的数值属于Sentinel-2 L1C产品的正常原始DN值，原因如下：

• Sentinel-2 L1C产品存储的是**大气表观反射率（TOA Reflectance）**的量化结果，欧空局官方默认设置的量化系数QUANTIFICATION_VALUE为10000，对应换算关系为：大气表观反射率 = 原始DN值 / 10000，该参数可以从你提到的MTD_MSIL1C.xml元数据文件中读取确认。
• 普通自然地表（植被、水体、土壤）的TOA反射率通常在0~1区间，但强反射目标（积雪、厚云、金属屋顶、盐滩、太阳耀斑区域）的反射率可以超过1，甚至达到2以上，对应原始DN值就会超过10000，你观测到的最大值23526对应反射率约2.35，属于高反射目标的正常数值范围。
• Sentinel-2 L1C采用16位无符号整型存储波段数据，可容纳的最大数值为65535，远高于常规地表的DN值范围，因此不会出现数值溢出问题。

二、使用3000作为缩放阈值的原因 #

你直接用全局最大值做归一化得到的图像偏暗，和使用3000阈值缩放效果正常的原因是一致的：

• 绝大多数无云、无积雪的常规地表目标，TOA反射率集中在0_{0.3区间，对应原始DN值范围就是0}3000。如果用全局最大值（比如你的23526）做归一化，常规地表的数值会被压缩到(3000/23526)*255 ≈ 32的低亮度区间，因此整体显示偏暗。
• 3000是遥感行业处理Sentinel-2真彩色可视化的通用经验阈值，将0_{3000的DN值线性映射到0}255的显示范围，超过3000的高反射目标直接截断为255（显示为白色），刚好可以最大化常规地表的显示对比度，符合人眼视觉习惯。
• 该阈值不是强制固定值，可以根据影像实际情况调整：雾霾较重的影像可以适当调低阈值，高反射地物（如裸露岩石、沙漠）占比高的影像可以适当调高到3500~4000区间。

补充参考操作 #

如果你需要先得到真实反射率再做可视化，可以执行如下转换：

import numpy as np
# 10000为官方量化系数
toa_reflectance = img / 10000
# 可视化时截断到0~0.3范围再映射到0~255
vis_img = (np.clip(toa_reflectance, 0, 0.3) * 255).astype(np.uint8)

上述代码和你使用的gdal_translate -scale 0 3000 0 255 -ot Byte命令效果完全一致。

内容的提问来源于stack exchange，提问作者Казыбек Аскарбек