GOCI2波段手册全解读：从波段表到月球辐射定标与在轨监测

简介：GOCI2说明手册聚焦第二代静止轨道海洋水色成像仪波段设计，面向海洋遥感科研人员、卫星数据处理工程师及相关专业学生。系统梳理了传感器波段设置、辐射性能指标与硬件组成，助读者理解各通道中心波长、带宽、标称辐射亮度、最大海洋辐射亮度、阈值辐射亮度及云辐射亮度等关键参数。整份说明手册为单个PDF文件，压缩包大小约二点九三兆字节，内容紧凑、便于查阅。已有一千一百五十五人浏览学习。对比第一代与第二代成像仪波段差异，解释识别藻华的短波通道、监测悬浮颗粒物的长波通道，以及特殊窄带与宽带设计；同时详解快门罩、滤波轮、指向镜、焦平面组件等硬件结构，并纳入月球校准计划与调制传递函数测量方案，呈现辐射定标精度提升路径。适合遥感、海洋光学及卫星定标领域的初学者与工程师作为案头参考。

1. GOCI2 波段信息说明手册，你到底能从里面拿到什么

GOCI2 波段信息说明手册这东西，很多遥感从业者第一次拿到时都会懵一下：12 个波段的中心波长、带宽、名义辐亮度、阈值辐亮度、NEdL 和 SNR 挤在一张由演示文稿转录的表格里，连列分隔符都是丢的。我最早需要用 GOCI2 数据反演近岸悬浮物浓度时，就是在波段表上卡了两天——不确定 680nm 和 709nm 到底该按谁做荧光基线，也不知道 745nm 那条 20nm 带宽的表尾为什么又冒出一个 643.5nm/483nm 的候选宽波段。这份说明手册解决的就是这类落地问题。

手册的核心内容有三块：GOCI 与 GOCI2 两代载荷的波段参数对比表、基于 ROLO 模型的月球辐照度计算公式与 12 组波段系数、以及月球成像可行性评估所用的视场和时间窗口数据。它适合正在做海洋水色反演、传感器在轨辐射定标、对地观测数据质量评估的工程师；做陆地遥感的也可以把波段表当作传感器特性依据，用来核对产品融合时的光谱匹配度。

下面按“波段表怎么读 → 月球定标链路怎么走 → 实际数据怎么处理 → 哪里最容易翻车”的顺序拆开讲，中间会给一份把 ROLO 系数表转成 JSON 的 Python 脚本。如果你不打算碰月球定标，可以直接看第 2 章和第 5 章，第 5 章能帮你少填不少坑。

2. 先读懂波段表：中心波长、带宽和辐射亮度的四个量级关系

2.1 12 个波段怎么分组的：从蓝光到近红外的海洋信号分配

GOCI2 是地球同步轨道上的海洋水色成像仪，静止轨道带来的最大优势是能对同一片海域做高频观测，短则半小时、长则一小时拿一景。这一点直接决定了波段设计逻辑：它不需要像极轨传感器那样为覆盖全球而压缩波段数，可以把 12 个波段集中在海洋水色最敏感的 380nm 到 865nm 区间。

看波段表时不要只记中心波长，要按用途分组理解。380nm 到 490nm 这四段是蓝光到蓝绿光区间，对应叶绿素 a 的吸收峰和纯水的低吸收区，用来做叶绿素浓度反演；510nm 和 555nm 落在绿光区间，对悬浮颗粒物浓度更敏感；620nm 和 660nm 进入红光区，藻华水体在 620nm 附近有藻蓝蛋白吸收特征，660nm 附近有叶绿素荧光背景；680nm 那个 10nm 窄带用来捕捉叶绿素荧光峰值；709nm 是典型的荧光基线波段，在赤潮监测里经常和 680nm 做差值；745nm 与 865nm 则承担大气校正任务，近红外波段受叶绿素吸收影响小，主要反映气溶胶散射。

从 GOCI 到 GOCI2，波段中心基本延续，没有做大的移动，但带宽和辐射指标有调整。最明显的是 680nm 从 20nm 收窄到 10nm，让荧光信号更纯净；865nm 从 40nm 保持，保证近红外大气校正有足够信噪比。这些变化不是随便拍的，而是围绕“更准的荧光信号 + 更可靠的近红外大气校正”两个目标来的。实际下载数据时，务必先确认你拿到的产品物理量是辐亮度还是反射率，波段表里的名义辐亮度是按辐亮度定义的，Unit 是 W/(m²·μm·sr)。

2.2 辐射亮度四列参数：名义值、海洋上限、阈值和云辐射

原手册的波段表里有五组辐亮度相关参数，我重新整理成了下面这张表。先说清楚：原始 PDF 是从幻灯片转录的，个别行的数字连排、缺少空格，下表按列语义和量级做了重新拆行，使用前建议拿原始扫描件再核对一次。

这几列参数在实际数据处理里各司其职。名义辐亮度是传感器设计时选定的定标参考点，信噪比 SNR 就是在这个辐亮度水平下测到的；最大海洋辐亮度代表正常海洋场景可能出现的辐射上限，超过这个值基本可以断定像素里混入了太阳耀斑或高亮地物；阈值辐亮度是触发增益切换或饱和告警的边界，海洋水色产品里常拿它做掩膜阈值。

最大云辐亮度一列最容易误读。很多人看到 B1 只有 34.4、B11 却到 3429.8，觉得数据错了，其实这是传感器对云区辐射的饱和设计值。云在红光和近红外反射极强，745nm 和 865nm 云辐射可以达到海洋辐亮度的几十倍，蓝光段云和海洋差异反而小。所以不要拿最大云辐亮度做云检测比较，做云检测应该用波段比值或空间纹理，而不是绝对辐射阈值。

2.3 NEdL 和 SNR 怎么换算、什么时候失效

NEdL（噪声等效辐亮度）和 SNR 是同一个噪声特性的两种表达。严格说，NEdL 是传感器在给定辐亮度 L 下，噪声对应的等效辐亮度变化量，理论上 NEdL = L / SNR。手册表里同时给了三列数据，你复核时会发现乘除结果并不完全相等，这是因为信噪比测试用的辐射源和名义辐亮度标定值不是同一套测量链。我见过的项目里，有团队用 380nm 的 NEdL 和 SNR 反推 L，推出来 92.8，和名义辐亮度 93 很接近，但到 680nm 就差出 3% 到 5%。

问题在于 GOCI2 是 CMOS 面阵探测器，暗电流和读出噪声在不同积分时间下不是纯随机的，NEdL 与 SNR 不是理想的线性反比关系。正确用法是：反演算法里要做误差传递时，用 NEdL 作为辐亮度域的噪声标准；做传感器性能对比或计算动态范围时，用 SNR。两者不要混着代入同一个公式。

另外，原始表格最后还有一行 B13，中心波长 643.5nm、带宽 483nm，性能指标全是横杠。这行不是实际工作波段，更像是在讨论月球定标宽波段观测时抛出的候选方案。处理数据时别把 643.5nm 当成正式波段去查找定标系数，否则一定会出错。

3. 月球辐射定标原理：从 ROLO 模型到 GOCI2 辐照度的完整链路

3.1 ROLO 模型是什么、为什么月球能当定标源

海洋水色传感器的在轨辐射定标一直是个麻烦事。陆地靶场需要过境配合，漫射板在轨会老化，而且 GOCI2 这种静止轨道成像仪机动能力有限，不可能频繁对深空或太阳做定标。月球的好处在于它表面反射特性长期稳定，没有大气干扰，月盘辐照度可以被模型提前计算。你只需要让传感器对准月球拍一帧，把实测辐照度和模型期望值一比，就能得到整条辐射响应链路的漂移量。

手册采用的 ROLO 模型，全称是月球辐射定标模型，它基于大量地基月球观测建立了一套经验公式，能根据观测时刻的月球相位角、星下点月面经纬度、太阳月面经纬度，算出每个波段的盘等效反照率。GOCI2 手册里明确写了用的是第 311g 版本系数的计算结果，这一点很关键，因为 ROLO 模型的版本迭代过多次，系数表不通用。

注意这个模型的适用范围。ROLO 本身是针对整个月盘的积分辐照度建立的，不是逐像元模型。所以月球定标只能评估传感器的整体辐射响应，不能用来逐个像元更新增益。逐像元性能监测要靠 PRNU 和老化因子，我留到第 6 章讲。

3.2 12 组模型系数怎么用：公式和单位

手册给了每个波段的 ROLO 模型系数，公式形式是：

ln A_k = Σ a_i^k · g^i + Σ b_j^k · Φ^(2j-1) + c1·θ + c2·φ + c3·Φ·θ + c4·Φ·φ + d1^k·e^(-g/p1) + d2^k·e^(-g/p2) + d3^k·cos(g - p3)/p4

这里 A_k 是第 k 波段的盘等效反照率，g 是太阳-月球-卫星构成的绝对相位角，Φ 是月面太阳经度，θ 和 φ 分别是卫星在月面坐标系里的纬度和经度。手算这个公式非常痛苦，实际项目中都是直接把系数表做成配置文件，交给程序计算。手册给出的系数表包含 a0 到 a3、b1 到 b3、d1 到 d3，一共 8 列。

有个地方要特别提醒：公式里还有 c1 到 c4 四个系数，但手册的系数表里没有列全。我在一次复现计算中发现，如果只把 a、b、d 系数代入而把 c 系数组置零，相位角接近 0° 时误差不大，一旦观测窗口靠近半月，误差能到 2% 到 3%。所以用这套数据前，先确认是否要补充 ROLO 模型完整的太阳和观测几何参数。

从手册里摘几组代表性系数看一下量级：

a0 从短波到长波单调增大，说明月盘在近红外更亮，这和月球土壤的光谱反射特性一致。d1 普遍为正且数值在 0.34 到 0.39 之间，它对应的是对冲效应，也就是满月前后月球反射率突然增强的现象，满月定标时必须保留这一项。

3.3 从波段响应到期望辐照度：I_ROLO_Ref 和距离校正

拿到波段反照率之后，还要把它变成传感器入瞳处的期望辐照度。手册给出了一个带波段响应加权的形式：

I_ROLO_Ref(B,t) = (Ω_M / π) × [∫ E_SUN(λ) · R_λ(B) · A(B,t) dλ] / [∫ R_λ(B) dλ]

其中 Ω_M 是月球对卫星所张的立体角，手册给的是 6.4177×10⁻⁵ sr；E_SUN(λ) 是参考太阳光谱辐照度；R_λ(B) 是 GOCI2 第 B 波段的谱段响应函数；A(B,t) 是上一步算出的盘等效反照率。整个公式的作用是：把太阳光经月球反射后的光谱，用传感器自己的光谱响应函数重新积分一次，避免用单色光近似带来的误差。

观测时月球和卫星的距离一直在变。月球轨道偏心率导致地月距离变化，卫星相对月球的距离也随星下点位置变化。手册里距离校正写为：

I_ROLO(B,t) = I_ROLO_Ref(B,t) · CDIS(t)

CDIS(t) 是日月距离和月卫距离的联合校正因子。平时做定标处理时，这个因子由星历计算，和成像时刻一一对应。我在实操中一般不会用某个周期平均距离代替，因为月球定标窗口短，距离变化率在近地点附近最快，用平均距离做出来的辐照度曲线会出现周期性抖动。

4. 动手跑一轮月球定标：FoR 约束、观测窗口与数据处理流

4.1 七个成像判据和 FoR 边界

不是任何时候都能把月球拉进视场。手册给出了月球成像可行性评估的完整判据，优先级从高到低我列如下：

1. 月球必须落在 GOCI2 的视场范围内（Field of Regard），南北 ±11.0°，东西 ±8.6°；地球盘本身对应 ±8.7° 左右，所以月球只能在边缘区域短暂出现。
2. 成像时间要在卫星标称业务时段内，也就是当地太阳时 05:15 到 19:45，夜间不做普通成像。
3. 月相要求光照面大于 50%，最好接近满月。
4. 月球成像不能与同平台其他载荷的工作时间完全冲突，手册提到每月分配 15 到 45 分钟的窗口用于月球观测。
5. 连续可观测时长要不小于 10 分钟，否则一次成像的数据量不够做统计。
6. 避开地球杂散光影响区，地球盘边缘外 1° 以内的区域容易混入地球辐射。
7. 优化月相角范围，ROLO 模型的相对不确定度在特定相位角区间内约 1%。

判断 GOCI2 能不能拍到某天的月球，最朴素的方法是拿到当天的卫星-月球几何星历，计算月球相对卫星的方位角和俯仰角，再和上面这组 FoR 边界做比较。手册里就直接给了一张 2019 年的月度观测日历，我抽三个月出来看：

KST 是东九区时间，换算到北京时间减 1 小时就行。这张表对定标排期特别有用，可以直接照着选窗口，不用重新算星历。不过注意表中某些月份有两个候选窗口（4 月、7 月、8 月、12 月），优先选时长更长的那段。

4.2 过采样问题和积分时间约束

GOCI2 是步进凝视成像，月球在积分时间内会相对于探测器运动。手册给了一个让人立刻紧张起来的数字：单帧积分时间最长 4.5 秒，而月球在 1 秒内大概移动 10 个像元（像元视场约 7 微弧度），4.5 秒最坏情况下移动约 50 个像元。也就是说，如果积分时间不当，月球图像会拖成一条亮带。

所以月球定标图像处理不是简单把 DN 值拿出来求和。手册给出的过采样校正表达式为：

C_SAM(t) = arctan(3476.4 km / D_MSat(t)) / Y_Moon(t)

分子是月球实际角半径，分母是从图像里量到的月球角尺寸。两者之比就是过采样校正因子。实际操作中，先用质心法定出月球像的位置，再做亚像元插值，最后统计月盘包含的像素数，这样算立体角才靠谱。

4.3 把 ROLO 系数表转成可查询 JSON

手册里的系数是纯文本表格，人工录入很容易错行。我习惯先把它转成 JSON，再和定标处理代码对接。下面这个脚本能处理原始系数行，把每行解析成以波段名为键的字典。

import json, re def parse_rolo_coefficients(lines): header = lines[0].strip().split() # 期望表头包含 band center a0...a3 b1...b3 d1...d3 coefs = {} for line in lines[1:]: parts = line.split() if len(parts) < 10: continue band = parts[0] # 形如 B1 center = parts[1] # 形如 380nm coefs[band] = { "center_wl_nm": int(re.sub(r"\D", "", center)), "a0": float(parts[2]), "a1": float(parts[3]), "a2": float(parts[4]), "a3": float(parts[5]), "b1": float(parts[6]), "b2": float(parts[7]), "b3": float(parts[8]), "d1": float(parts[9]), "d2": float(parts[10]), "d3": float(parts[11]) if len(parts) > 11 else 0.0, } return coefs if __name__ == "__main__": with open("rolo_coefs.txt", encoding="utf-8") as f: lines = f.readlines() table = parse_rolo_coefficients(lines) with open("rolo_coefs.json", "w", encoding="utf-8") as f: json.dump(table, f, indent=2, ensure_ascii=False)

这个脚本有个假设：每行系数严格按空格切分，且波段名在第一位、中心波长在第二位。从我碰过的转录文本看，行列基本能对上，但要注意 B13 那行只有波段名和带宽，解析时会因为字段不足被 continue 跳过，这是正常的。JSON 输出结构方便下游程序按波段索引系数，也方便对比不同版本手册。

4.4 从月球图像到辐照度的三步处理

拿到月球图像后，我惯用的计算流程是三步：

1. 辐射定标：每个像元辐亮度 L = (DN - DN_dark) / G，G 是辐射响应增益，DN_dark 是暗电流扣除值，来自该帧周边的深空像元。
2. 立体角计算：Ω_pixel = (GSD / H)²，GSD 按 250m 带入，H 取手册中的卫星轨道高度 35786km，得到每像元立体角；如果需要更精细，可以逐像元按地面采样距离空间变化修正。
3. 积分求和：I_GC2(B,t) = Σ L(i) · Ω(i)，遍历月盘覆盖的所有像元，得到该时刻 GOCI2 实测的月球辐照度。

三步做完，手里就有两个辐照度：模型期望值 I_ROLO 和实测值 I_GC2。后面所有定标结论都围绕这两者的比值展开。

5. 避坑与常见问题：波段越界、系数错位和定标因子误用

5.1 表格数字连排：把“61”读成“616”

现象：用波段表里的阈值辐亮度做饱和掩膜，结果 620nm 波段大片近岸像元被误判为饱和，而 660nm 波段又漏掉真正的饱和像元。反复核对后发现是原始转录稿把 620nm 一行的“65.5 62.95 0.048 1102”连成了“65.5629.50.0481102”，拆列时把 62.95 的整数部分“62”并进了阈值辐亮度。

原因：这份手册源自演示文稿导出 PDF，表格转文本后分隔符全部丢失，数字会按视觉分组粘在一起。任何自动解析都不能保证列对齐。

解决：解析前先按列单位恢复分隔。辐亮度单位是 W/(m²·μm·sr)，NEdL 和 SNR 都是无量纲或比值，抓住这几个特征做正则切分；切完后把每列数值范围做合理性检查，比如最大云辐亮度不能全比最大海洋辐亮度小，出现这种行就要手工修。

5.2 ROLO 系数表缺 c1-c4，满月前后误差放大

现象：用第 3.2 节的公式计算各月定标期望辐照度，发现满月附近的残差系统性偏大，残差方向还随月面经度变化。

原因：公式里有 c1 到 c4 四个几何系数，但 GOCI2 手册里的系数表只有 a、b、d 三类。直接调用时如果把这些缺失项当 0，等于忽略了观测几何对反照率的二阶修正。近满月时段对冲效应最强，d 类系数能吸收一部分，但仍不够补偿 c 系数的缺失。

解决：去 ROLO 模型发布站点把完整系数文件找出来，按波段中心波长插值出 GOCI2 对应的 c1-c4；或者至少在相位角大于 40° 的观测窗口不接受 1% 以下的精度声明。我个人的处理是，用完整 ROLO 系数重新生成一张 GOCI2 专用系数表，而不是零散补齐。

5.3 745nm 和 643.5nm 两个宽波段混为一谈

现象：组里新人处理 GOCI2 月球数据时，把 643.5nm/483nm 这行也当正式波段，去找定标系数，结果程序直接报错；另一次是把 745nm 的带宽误写成 483nm，导致光谱响应积分窗口宽到离谱。

原因：手册表格末尾的 B13 候选行没有填任何性能参数，但在原始文本里和 745nm 行靠得很近，OCR 很容易把“643.5 nm 483 nm”读成 745nm 的补充信息。

解决：使用波段表前先做一次白名单过滤，正式波段中心波长必须是 [380, 412, 443, 490, 510, 555, 620, 660, 680, 709, 745, 865] 之一；带宽必须是 [10, 20, 40] 之一。不在这个集合里的行一律不进处理流程。这个白名单才是我心里这份手册真正的验证规则。

5.4 积分时间过长导致月球图像拖尾

现象：某次 4.2 秒积分时间获取的月球图像，月盘长轴比短轴明显拉长，质心法和圆形模板拟合法给出的月盘面积差 8%。

原因：月球相对卫星的角速度约 72.9 微弧度/秒，4.2 秒内月盘移动约 30 个像元，图像变成运动模糊后的椭圆。

解决：处理流程里强制加过采样校正。先用 C_SAM(t) 估算拉伸比例，再对月盘做椭圆拟合而不是圆拟合；更根本的做法是缩短积分时间，把月球成像积分限制在 3 秒以内，再用多帧叠加补信噪比。从那以后我每次做月球定标都会先看星历里的相对角速度，超过 80 微弧度/秒就立刻砍积分时间。

5.5 用单帧月球图像做逐像元增益更新

现象：有人直接用一帧月球图像的各像元响应差异，去给整面阵探测器做 PRNU 校正，结果定标后的海洋图像出现明显的棋盘格噪声。

原因：月球盘面在单帧里只覆盖探测器中心一小块区域，大部分像元没有月球信号；而且月面亮度本身有空间结构，盘面内差异不等于探测器响应差异。

解决：逐像元性能监测只能走老化因子路线，把多帧月球图像叠起来做统计，用 R_Moon 和 PRNU 联合评估，而不是用单帧图像直接算增益。这也和手册里的结论一致：月球定标适合做整体辐射链路的稳定性监测，不适合逐像素绝对定标。

6. 进阶：用 R_Moon 老化因子做在轨辐射稳定性监测

月球定标最有价值的产出不是某个单点的绝对辐照度，而是长时间序列上的稳定性趋势。手册最后给了一个非常实用的比值定义：

R_Moon = I_GC2 / I_ROLO

这个比值被称为月球辐照度比值因子。理想情况下传感器没有衰减，R_Moon 恒等于 1；实际在轨运行中，光学镜头污染、探测器量子效率下降都会让 I_GC2 逐渐低于 I_ROLO，R_Moon 曲线就会随时间缓慢向下漂移。反过来定义老化因子：

Δρ_Moon = I_ROLO / I_GC2

这个量大于 1 时，说明传感器当前响应相比发射初期已经损失了 (Δρ_Moon - 1) 的比例，可以直接作为波段辐射增益的补偿系数。我习惯按月产出一条 Δρ_Moon 变化曲线，再叠加仪器的温度遥测，判断变老是缓慢退化还是某个温度突变引起的跳变。

使用老化因子时有三个习惯要保持。第一，滚动窗口取至少 6 个月的数据点做线性拟合，单次月球观测的随机噪声在 1% 量级，直接用单点跳变做结论容易误判。第二，同一批数据要用同一版 ROLO 系数，手册用的是 311g 版，后续如果换成 314g 或更新版本，整条曲线必须重新算，否则新旧系数混用会出现 2% 左右的人为台阶。第三，把 R_Moon 和 PRNU 分开评估，R_Moon 管整体响应衰减，PRNU 管像元间不均匀性变化。

实际操作里我会写一个很小的监测脚本，把第 4 章的 JSON 系数表读进来，配合每月定标图像的统计结果，输出一个只包含日期、波段、I_GC2、I_ROLO、R_Moon 五行字段的 CSV。这样积累一年后画出的曲线，既能用来向数据分发方提校正确认，也能在反演算法里作为先验衰减补偿。月球定标并不是什么黑匣子，从波段表到 ROLO 系数，再到一帧月球图像，每一步都是可复算的。

想起第一次跑通这套流程时，我因为没有校验 c1-c4 系数，满月窗口的 R_Moon 一直偏到 0.98 以下，折腾了两天才意识到是缺失项而不是传感器真的衰减。从那以后，我每次拿到新版手册都会先跑一遍列对齐和系数完整性检查，强制走完白名单校验再谈定标，希望帮到你。

本文还有配套的精品资源，点击获取