3. 任务载荷图像数据获取与特性分析

星上图像预处理的第一步,不是滤波,也不是压缩,而是把“载荷到底给了我们什么数据”这件事弄清楚。同一颗卫星上,CCD、CMOS、红外焦平面和高光谱成像仪输出的原始数据在物理量纲、位宽、排列方式和噪声结构上差异很大。如果预处理算法直接套用可见光图像的假设,往往会在红外或高光谱通道上产生系统性偏差。本节围绕成像机理、原始数据格式与元数据、噪声来源与统计特性、典型场景退化因素四条线索展开,目标是让读者拿到一帧原始数据后,能判断它“从哪来、长什么样、脏在哪、为什么退化”。

一、四类载荷的成像机理差异

CCD 采用电荷转移读出,光生电荷在势阱中逐行转移到输出放大器,因此具有较高的填充因子和一致性,但存在电荷转移效率损失与拖尾(smear)风险,尤其在行转移期间仍有光照时。CMOS 每个像元自带放大与读出电路,支持窗口化读出和较高帧率,代价是像元间固定模式噪声(FPN)和列级放大器失配更明显。红外焦平面(如制冷型 InSb、MCT 或非制冷微测辐射热计)响应的是目标与背景的辐射差异,输出常为 14 bit 或 16 bit,且受探测器自身热辐射、1/f 噪声和坏元影响显著。高光谱成像仪在光谱维上把入瞳辐射色散到面阵探测器上,形成“空间×光谱”立方体,其数据本质是三维的,单波段信噪比通常低于全色通道。

这四类机理决定了原始数据的组织方式:CCD/CMOS 输出二维帧,红外输出二维帧但动态范围更大,高光谱输出按波段排列的二维切片序列或推扫式行序列。理解这一点,才能正确选择后续的暗电流扣除、非均匀性校正和光谱定标策略。

四类载荷成像机理与原始数据形态对照 CCD 电荷转移读出 势阱积累 → 行转移 → 输出放大器 二维帧 · 12~14 bit 风险:拖尾 / CTE 损失 填充因子高、一致性好 CMOS 像元内放大读出 光电转换 → 像元放大 → 列并行读出 二维帧 · 10~12 bit 风险:FPN / 列失配 支持窗口化、高帧率 红外焦平面 辐射差异响应 目标/背景辐射差 → 积分 → 读出 二维帧 · 14~16 bit 风险:1/f 噪声 / 坏元 动态范围大、需制冷 高光谱成像仪 色散分光 + 面阵探测 入瞳辐射 → 光栅色散 → 空间×光谱面阵 三维立方体 · 12~14 bit 风险:单波段 SNR 低 推扫/摆扫获取 原始数据层:像素值 + 元数据 像素值(DN) · 位宽 · 行/列排列 · 波段数 · 积分时间 · 增益 · 探测器温度 · 成像时刻 · 姿态/轨道参数 · 定标系数 缺少元数据 → 无法完成辐射定标与几何校正 噪声来源 光子散粒噪声 · 读出噪声 · 暗电流噪声 固定模式噪声 FPN · 1/f 噪声 · 量化噪声 坏元 / 列缺陷 / 光谱串扰 典型场景退化因素 低照度 → 信噪比下降 · 高动态场景 → 饱和/截断 运动模糊 · 大气散射/水汽吸收 · 云层遮挡 太阳耀斑 · 探测器温漂 · 光谱混叠

二、原始图像数据格式与元数据

星上原始数据通常以“帧头 + 像素数据 + 帧尾校验”的格式打包。帧头中必须包含成像时刻、积分时间、增益档位、探测器温度、波段编号和行号,这些字段直接决定后续定标能否进行。以某推扫式高光谱载荷为例,一帧原始数据可表示为:

帧头(64 Byte):
  frame_id      : uint32   帧序号
  timestamp     : uint64   成像时刻(UTC 微秒)
  integration   : uint16   积分时间(μs)
  gain          : uint8    增益档位(0~7)
  det_temp      : int16    探测器温度(0.01℃)
  band_id       : uint16   波段编号
  line_num      : uint32   推扫行号
像素数据:
  按行排列,每行 N 个像元,每像元 14 bit 左对齐存入 16 bit
帧尾:
  crc32         : uint32   校验

元数据缺失或错误是星上预处理最常见的“隐性故障”。例如积分时间字段若被误写为默认值,辐射定标后的辐亮度将整体偏移;探测器温度未记录,则非均匀性校正系数无法按温度插值。因此,预处理软件在读取原始数据时,应首先校验帧头字段的合法范围,再进入像素处理。

三、噪声来源与统计特性

原始图像的噪声并非单一类型,而是多种机制叠加。光子散粒噪声服从泊松分布,其方差等于均值,在低照度下主导;读出噪声近似高斯分布,与信号无关;暗电流噪声随积分时间和温度指数增长;固定模式噪声(FPN)在空间上固定、时间上缓变,表现为列条纹或像元偏置差异;1/f 噪声在红外探测器中尤为突出,功率谱密度与频率成反比。量化噪声由 ADC 位数决定,14 bit 量化时理论信噪比上限约为 86 dB。

实际分析时,可用“暗场帧”和“平场帧”分离噪声成分:暗场帧(无光照、相同积分时间)的均值反映暗电流偏置,标准差反映暗电流噪声与读出噪声之和;平场帧(均匀光照)的空间标准差反映 FPN 与散粒噪声的叠加。若暗场帧标准差随积分时间线性增长,说明暗电流噪声主导;若平场帧出现周期性列条纹,则 FPN 需要优先校正。

四、典型场景图像质量退化因素

退化因素可分为辐射域、几何域和光谱域三类。辐射域包括低照度信噪比下降、高动态场景饱和截断、太阳耀斑引起的局部过曝;几何域包括平台运动引起的运动模糊、姿态抖动导致的像移、大气折射引起的几何畸变;光谱域包括大气水汽吸收造成的波段间响应差异、光谱串扰导致的波段混叠。以低照度场景为例,若目标辐亮度低于探测器噪声等效辐亮度的 3 倍,单帧图像的信噪比将低于 10 dB,此时直接进行边缘检测会引入大量虚假边缘,应先做多帧累加或自适应降噪。

判断一帧原始数据是否“可用”,建议按以下顺序检查:帧头字段是否完整合法 → 暗场帧均值与标准差是否在预期范围 → 平场帧是否存在明显列条纹或坏元 → 目标区域是否饱和或过暗 → 波段间响应是否一致。任一环节异常,都应先定位原因,再决定是否进入预处理链。

综上,任务载荷图像数据获取与特性分析的核心,是把“载荷机理—数据格式—噪声统计—退化因素”串成一条可验证的链路。只有先看清原始数据的真实面貌,后续的辐射校正、去噪、几何校正和光谱复原才有可靠的输入基础。