简介《GOCI2 说明手册 波段信息说明》面向海洋遥感、水色反演与卫星辐射定标方向的学习者和工程人员集中梳理 GOCI-II 的波段设置与关键性能指标。文档以单份 PDF 形式提供压缩包约 2.93MB内含 1 个文件便于快速检索与对照查阅。内容围绕 GOCI 与 GOCI-II 的波段对比展开列出各波段中心波长、带宽、名义辐射亮度、最大海洋辐射亮度、阈值辐射亮度、最大云辐射亮度及 NEdLSNR 等参数并涉及快门罩、指向镜、焦平面组件、热沉等设备结构。文中还结合 GSICS/IVOS 月球校准研讨会资料介绍 GOCI-II 月球校准计划、MTF 测量计划与月球成像可行性评估有助于理解仪器观测能力、辐射定标思路和海洋环境监测应用。目前已有 1161 人学习。1. 一句话说清 GOCI-II地球同步轨道上那台“看海水的分光计”手机对着海面拍一张出来就是一片蓝。可海洋研究者盯着同一片蓝想知道的是水里叶绿素有多少、悬浮泥沙有多密、有没有藻类正在暴发。把“一片蓝”拆成十几个能互相加减的数字靠的就是 GOCI 这类海洋色度成像仪——它固定在静止轨道上一天之内反复回看同一片海把颜色切成十几个很窄的光谱窗口。GOCI-II 是它的第二代波段从 380 nm 一直排到 865 nm还配了一整套月球定标方案用来盯住自己的“眼睛”有没有变花。搞懂这套东西等于搞懂遥感数据里“颜色”是怎么一步步变成“浓度”的。2. 波段表到底在说什么从 380 nm 到 865 nm 的十二个窗口翻开这份手册最先撞进眼里的往往是那张密密麻麻的波段表。十二行窄波段加一行宽通道每行后面跟着中心波长、带宽、四个辐亮度数值和一个噪声指标。第一次看很容易当成一堆天书数字划过去但这张表其实是整台仪器的“能力清单”哪些东西它看得见多亮的东西它会饱和多小的差异它分得出来全在表里写着。2.1 中心波长与带宽每一行都是一扇窄窗先说最基础的两个词。中心波长是这个通道最敏感的位置带宽是它愿意接收的有效光谱宽度。打个比方中心波长是窗框中间那条线带宽是窗户有多宽。窗户越窄看东西越“挑”越能分辨相邻两条光谱线的差别窗户越宽一次进来的光越多信号更足但会把相邻成分混在一起。下面这张表按手册给的数据整理中心波长和带宽两列是固定的波段中心波长带宽这个窗口主要关注什么B1380 nm20 nm短波端给大气校正托底B2412 nm20 nm黄色物质、悬浮物的强吸收区B3443 nm20 nm叶绿素吸收峰附近B4490 nm20 nm叶绿素与水体透明度B5510 nm20 nm低浓度叶绿素的细微变化B6555 nm20 nm悬浮颗粒物的反射峰B7620 nm20 nm藻华、色素信号B8660 nm20 nm叶绿素荧光附近B9680 nm10 nm荧光峰带宽收窄到 10 nmB10709 nm10 nm藻华在水面形成的“红边”B11745 nm20 nm大气校正的近红外参考B12865 nm40 nm气溶胶、近红外参考宽通道643.5 nm483 nm跨可见光的大窗口指标栏留空看这张表有两个容易忽略的细节。第一B9 和 B10 的带宽只有 10 nm比其它通道窄一半——这两个位置对应的是很窄的光谱特征窗子开大了就会被旁边的信号冲淡。第二B12 反过来带宽拉到 40 nm因为近红外这段本来就没有尖锐特征宽一点能多收集能量、换来更好的信噪比。一窄一宽不是随意定的都是拿“要分辨什么”和“能收到多少光”这两件事反复权衡的结果。至于最后那行 643.5 nm、带宽 483 nm 的宽通道它的辐亮度指标在表里是空的。这一点很值得注意指标留空通常意味着这个通道不参与常规的定量辐亮度流程它的角色更接近辅助或监视而不是出定量产品。2.2 四个辐亮度数值一台仪器的“量程刻度”波段表里最让人犯迷糊的是一行里出现四个名字都带“辐亮度”的数字。它们都是同一个物理量单位是 W·m⁻²·μm⁻¹·sr⁻¹差别只在于描述“在什么条件下”名义辐亮度正常海面情况下大概会遇到的亮度水平是设计时锚定的基准工作点。最大海洋辐射亮度海面在明亮条件下可能达到的上限超过它海洋目标的观测就开始失真。阈值辐亮度一条警戒线用来判断当前信号是否还在可靠范围内。最大云辐射亮度云顶反射的亮度能到多高。这个值远大于前三个因为白云反射的能量要比海水高一个量级。以 B1 为例名义辐亮度是 93最大海洋辐射亮度 139.5阈值辐亮度 143.16最大云辐射亮度 634.4。四个数字一摆量程的层次立刻清楚了正常观测在 93 附近海面再亮也就顶到 140 上下而云的亮度能冲到 630 以上。换句话说设计者心里同时装着“海”和“云”两种极端场景而不是只盯着海。这里有个很实用的读法名义辐亮度和最大海洋辐射亮度之间的差距反映的是这台仪器对“海面变亮”这件事留了多大余量最大云辐射亮度和它们之间的巨大落差则是留给云检测和大气校正的空间。如果只记住一个数就失去了这张表一半的信息。2.3 噪声指标能“看见”和能“分辨”是两件事表里最后一列是 NEdL中文常译成噪声等效辐亮度差。名字拗口意思其实很朴素仪器自身抖动产生的那点“底噪”折算成辐亮度是多少。这个数值越小说明它越能分辨两个靠得很近的亮度水平。B1 的 NEdL 是 0.09B12 是 0.015。单看数字没感觉跟同行的名义辐亮度放一起就有味道了B1 的名义辐亮度是 93底噪 0.09相当于在 93 这个量级上还能分辨千分之一左右的差别。这就是为什么波段表要同时给出名义辐亮度和噪声——光有量程只是说明“能看见多亮的东西”加上噪声才说明“能看清多细的差别”。一台量程很大但噪声也大的仪器拍出来的图可能很亮但一片糊。提示读波段表时先把中心波长和带宽扫一遍建立整体印象再回头只看你关心的那两三个波段对应的辐亮度和噪声最后才通读全表。反过来从第一个数字开始逐格读很容易读到第四行就失去方向。3. 为什么要用月球做定标一架不会自动变旧的“天上标尺”波段表解决的是“设计上能看多好”可仪器上了天实际表现会一点点漂。传感器会老化滤光片会变镜头会脏电子学增益也会飘。要判断数据还能不能信就必须有个稳定的参照物反复比对。问题是地面上那把“标准尺”没法搬到天上天上的东西自己也会变。3.1 星上定标器的难处参照物自己也会脏常规做法是在卫星上装一块漫反射板让它把太阳光均匀散射进镜头当作已知亮度的参照。听起来很完美但漫反射板暴露在太空环境里会被紫外辐照、污染、粒子轰击慢慢改变反射特性。于是出现一个尴尬的循环你想用板子校准仪器可板子自己也在退化你并不知道它退化了多少。所以工程师又给它配了专门盯板的装置——手册里提到的漫反射板老化监视装置就是干这个的。但这只是缓解不是根治。参照物自己会变旧这个根本矛盾始终存在。月球的好处恰好在这里。它没有大气表面在人类观测的时间尺度上几乎不变反射特性稳定得可以当作“宇宙标准灯”。而且它就在天上任何一颗能转动的卫星都有机会把它框进视场。用月球做定标相当于绕开了星上参照物会退化的死结。3.2 七条判据先看能不能拍到再看拍得好不好想法很美好但月球不是想拍就能拍。手册里列了一张分优先级的可行性判据读下来会发现它的排序逻辑非常“工程师”先保证能拍到再保证拍得好。优先级 1月球必须落在视场里。视场在东西方向是 ±8.7 度南北方向是 ±11.0 度地球圆盘正好对应 ±8.7 度。这条最硬进不了视场后面全都白搭。手册还特别注明不打算为此做整星姿态机动也就是说能不能拍到完全取决于自然几何不去为它额外调整姿态。优先级 2要在标称工作时段内。每天的运行窗口是 05:15 到 19:45出了这个区间仪器状态就不在标称条件下。优先级 3月球要有足够的亮面。被照亮的部分要超过一半满月当然最好。月相太小信号弱噪声占比就上来了。优先级 4不能和其它载荷抢时间。同一颗卫星上还有别的仪器在工作月球成像的时段不能完全和它们重叠。优先级 5单次可连续观测时间要够长。判据写的是不少于 10 分钟待定。时间太短积分不充分。优先级 6避开地球杂散光。地球本身很亮反射光会从视场边缘漏进来。手册把距地球圆盘 1 度以内的区域划成潜在污染区月球图像要尽量避开。优先级 7相位角要落在合适区间。这一条跟定标模型直接相关。月球反射率随相位角变化而模型本身在这个维度上有约 1% 的相对不确定度所以要挑一个模型算得最准的相位区间。七条判据摆在一起能看出一个很清楚的态度月球定标不是“有机会就拍”而是“在满足一堆约束的前提下挑一个最优窗口”。每一条都在把可选项往下砍砍到最后剩下的窗口其实相当窄。3.3 视场、指向与那堆机械件手册里那张仪器结构图上标了一串部件快门盖、支撑结构、卫星接口、四个双脚架、滤光轮、带热分流器的焦平面组件、指向镜、遮光罩、望远镜机构。还有一行小字注明某个辐射增强系统并不属于这台仪器——这种备注在手册里很常见作用是划清边界免得读者把别人的指标算到自己头上。跟月球成像关系最直接的是指向镜和遮光罩。指向镜决定镜头能不能转向月球所在的方向遮光罩负责挡住视场外的强光。四个双脚架把整机固定在卫星接口上热分流器负责把焦平面的热量导走——温度一飘暗电流就飘定标结果也跟着飘。看起来都是机械细节实际上每一条都在影响最后那个数字的可靠性。视场的设计也很有意思地球圆盘占 8.7 度仪器视场留 1.1 度给月球和恒星再留 0.5 度应对指向偏差。这个余量分配说明月球成像并不是“顺便蹭一下”而是从设计阶段就预留了位置。3.4 一张排期表透露的现实约束手册里附了一张 2019 年的逐月月球成像排期。数据很直白月份可观测天数推荐观测时长1 月2 天约 30 分钟2 月3 天约 30 分钟3 月1 天约 38 分钟4 月1 天约 46 分钟5 月1 天约 29 分钟6 月2 天约 30 分钟7 月2 天约 39 分钟8 月2 天约 43 分钟9 月1 天约 30 分钟10 月2 天约 11 分钟11 月4 天约 12 分钟12 月3 天约 44 分钟一个月能拍的日子只有一到四天可用的连续时间从 11 分钟到 46 分钟不等。10 月和 11 月特别短只有十一二分钟。这意味着整年的定标机会其实非常有限每一次拍摄都得当成正式任务来准备不能指望“这次没拍好下个月补”。这也解释了为什么判据要列得那么细——机会成本太高必须提前筛掉所有不合格的窗口。4. 从一张月球照片到一个数字ROLO 模型这条计算链拍到月球只是开始。真正的定标要回答一个更精确的问题仪器测到的月球有多亮和月球“应该有多亮”差了多少。这中间隔着一整条计算链每一环都需要单独交代清楚。4.1 第一步算出“月亮本来应该有多亮”月球不是一盏亮度恒定的灯。它的反射率会随相位角变化——满月时最亮月牙时最暗而且这个变化不是线性的。同时观测几何里月面的纬度、经度太阳照射的经度都会影响你看到的那部分月面反射多少光。国际上有一套长期维护的月球反射率模型手册里用的就是它。核心公式是圆盘等效反射率的对数值由几个多项式项和指数项叠加而成。系数按波段各自一套一共十二组。挑其中一组写成代码结构是这样的import numpy as np # 计算某个波段的圆盘等效反射率 Ak # g : 绝对相位角单位度 # theta : 观测点对应的月面纬度 # phi : 观测点对应的月面经度 # Phi : 太阳的月面经度 # coef : 该波段的一组系数含 a[4]、b[3]、c[4]、d[3]、p[4] def rolo_disk_reflectance(g, theta, phi, Phi, coef): a, b, c, d, p coef[a], coef[b], coef[c], coef[d], coef[p] g_rad np.deg2rad(g) # 相位角的多项式部分对应 a0 a1*g a2*g^2 a3*g^3 s sum(a[i] * g_rad ** i for i in range(4)) # 太阳经度的奇次幂部分Phi 的 1、3、5 次方 s sum(b[j] * Phi ** (2 * j 1) for j in range(3)) # 经纬度交叉项四个全局系数 s c[0] * theta c[1] * phi c[2] * Phi * theta c[3] * Phi * phi # 两个指数衰减项描述小相位角附近的陡峭变化 s d[0] * np.exp(-g_rad / p[0]) s d[1] * np.exp(-g_rad / p[1]) # 余弦项对应 cos(g - p3) / p4 s d[2] * np.cos(g_rad - p[2]) / p[3] return np.exp(s) # 反射率本身是对数结果的指数这段代码有两处值得单独说。第一np.exp(s)这一步说明原始公式给的是对数形式系数都是在对数域里拟合出来的算完必须取指数才能得到反射率。第二指数项里的p[0]和p[1]是衰减尺度它们决定了相位角从 0 度往上走时反射率下降得有多快。这两个参数在不同波段之间差异不小不能混用。系数本身按波段给全每个波段十四五个数字看起来庞杂但结构是统一的a管相位角的主趋势b管太阳经度c管交叉调制d和p管小角度附近的修正。只要结构对上了换波段只需要换系数。4.2 第二步把反射率积成辐照度有了反射率还不够仪器测的是辐照度。要把前者变成后者得走一遍积分积分月球参考辐照度 立体角 / π× 太阳光谱辐照度 × 仪器光谱响应 × 圆盘等效反射率的积分 ÷ 光谱响应的积分这条式子里太阳光谱辐照度是参考值仪器光谱响应是每个波段自己的响应曲线反射率就是上一步算出来的结果。除以光谱响应的积分本质是在做归一化——不同波段的响应宽度不一样不归一化就没法比较。手册里给了一个固定常数月球对应的立体角是 6.4177×10⁻⁵ 球面度。这个值在整个计算里是常量不需要每次重算。算完参考辐照度还要做距离校正。地球绕太阳的轨道是椭圆月球绕地球的轨道也不是正圆两者之间的距离每天都在变接收到的能量自然也跟着变。所以有一个距离校正因子把测量值折算到统一距离下距离校正后的测量辐照度 参考辐照度 × 距离校正因子这一步很容易被忽略因为它的影响不像噪声那么显眼但它是系统性的。少了这一步不同月份拍到的数据就没法直接对比。4.3 第三步把图像里的像素加回去前面两步算的都是“应该值”现在轮到“实测值”。月球在图像上占一片像素每个像素都有自己的立体角。把每个像素的辐亮度乘以它对应的立体角再全部加起来就得到了积分月球辐照度# 从月球图像反算积分月球辐照度 # 前提图像已经做过过采样校正 GSD 250.0 # 星下点地面采样距离单位 m手册给出优于 250 m H 35786.0e3 # 卫星轨道高度单位 m OMEGA_PIX (GSD / H) ** 2 # 单个像元张开的立体角单位球面度 def integrated_lunar_irradiance(radiance_img, moon_mask): # radiance_img : 二维数组每个像素的月球辐亮度 # moon_mask : 布尔数组标出落在月面上的像素 return float(np.sum(radiance_img[moon_mask] * OMEGA_PIX))这段代码的关键在OMEGA_PIX这个量。它是地面采样距离除以轨道高度再平方单位是球面度。地面采样距离越小单个像素张开的立体角越小同样的辐亮度贡献的能量也越小。手册里给出的地面采样距离优于 250 米所以这个值是个上限实际计算时如果手上有更精确的数字应该替换掉。另外要注意moon_mask。月球在图像上不是一个完整的圆盘边缘会因为曝光和指向偏差变得模糊所以哪些像素算“月面”需要单独判定。判得太松会把背景天光算进去判得太紧会漏掉边缘两种做法都会让结果偏低或偏高。4.4 第四步过采样校正、比值因子与老化因子到这一步还有两个坑。第一个坑是过采样。仪器在积分时间内不是静止的月球在动卫星也在动结果就是月面在图像上被“拖”出了一段。手册算过一笔账积分时间 1 秒月球大约移动 10 个像素最坏情况积分 4.5 秒能移动 50 个像素。卫星平均轨道角速度约 72.9 微弧度每秒单个像素视场约 7 微弧度两个数字一除正好对上。校正的办法是拿月球的真实角直径和图像上测到的角直径做比较过采样校正系数 arctan(3476.4 km / 月地距离) ÷ 图像上测得的月球角直径3476.4 公里是月球直径这个值是固定的。分母是图像上实测的角尺寸两者一比就知道图像被拉伸了多少倍再按比例修正回去。手册也坦白说如果判定确实需要校正更精细的模型要把月球运动、卫星运动和指向稳定性一起考虑那是下一步的工作。第二个坑是不能逐像素重算增益。这一点常被误解。月球在单张图像里只覆盖了一部分探测器像素而参考模型给的是整个月球的积分辐照度。用部分像素去反推每一个像素的绝对增益在数学上是不成立的。所以手册没有走这条路而是定义了比值因子比值因子 仪器测得的积分月球辐照度 ÷ 模型给出的参考辐照度这个比值反映的是整机层面的辐射性能变化不落到单个像素。把比值取倒数就得到月球老化因子月球老化因子 参考辐照度 ÷ 实测辐照度老化因子配合像素响应非均匀性就能评估每个像素的辐射性能变化。注意这里的逻辑顺序先用整体比值判断“仪器整体衰减了多少”再用非均匀性把整体衰减分摊到各个像素。这个两步走的做法正是被“月球只覆盖部分像素”这个技术限制逼出来的。5. 常见误区与易混概念这几个词最容易被混着用读到这儿概念大体串起来了但有几个地方特别容易在脑子里糊成一团。下面五条是初学者最常踩的坑每一条都值得停下来对一遍。5.1 误区一波段越多、带宽越窄就代表越先进易混点看到十二个窄波段加一个宽通道直觉上会认为波段数量多、带宽窄就是技术更强。为什么混带宽窄确实能提高光谱分辨能力波段多确实能覆盖更多现象这两件事单独看都对于是很容易直接画等号。怎么区分带宽是拿信噪比换来的。窗子开窄进来的人光少了噪声占比就上去了。B9 和 B10 敢用 10 nm是因为那两个位置有明确的光谱特征值得为之付出代价B12 用 40 nm是因为近红外没有尖锐特征宽一点反而更划算。所以判断一台仪器不能数波段个数要看每个波段的设计意图和它配套的噪声指标是否匹配。5.2 误区二把四个辐亮度数值当成同一个东西的不同叫法易混点名义辐亮度、最大海洋辐射亮度、阈值辐亮度、最大云辐射亮度四个名字都带“辐亮度”单位也一样看起来像是重复标注。为什么混它们确实是同一个物理量只是对应的观测条件不同名字又长得很像扫一眼就滑过去了。怎么区分把四个数按大小排一下层次立刻出来。以 B1 为例名义 93、最大海洋 139.5、阈值 143.16、最大云 634.4。名义是日常工作的锚点最大海洋是海面变亮的余量上限阈值是警戒线最大云是另一个量级的极端场景。四个数字描述的是一台仪器的完整量程剖面缺一个就少一块信息。5.3 误区三以为月球定标能给每个像素重新算增益易混点“定标”这个词天然带着“把每个像素都校准一遍”的联想。为什么混很多地面相机的平场校正确实是逐像素做的习惯成自然就会以为天上的定标也一样。怎么区分月球在单张图像里只覆盖部分探测器像素参考模型给的又是整月球的积分值用局部推全体在数学上不成立。所以手册走的是“先算整体比值因子再用非均匀性分摊”的路线。逐像素的绝对增益更新月球定标做不到。5.4 误区四把月球当成亮度恒定的标准灯易混点既然说月球表面稳定就容易进一步认为它每次看起来都一样亮。为什么混“稳定”指的是月面反射特性在长时间尺度上不变不是指观测到的亮度不变。怎么区分亮度随相位角剧烈变化满月和上下弦能差出好几倍。所以才有优先级 7 那条判据——要把相位角限制在模型算得最准的区间手册给的不确定度大约是 1%。除此之外还有距离校正、过采样校正。月球是稳定的参照物不是省事的参照物。5.5 误区五把 643.5 nm 那行宽通道和十二个窄波段当成一回事易混点表格里并排列着很容易顺口读成“十三个波段”。为什么混格式统一位置相邻视觉上就是同一类东西。怎么区分看指标栏。十二个窄波段每行都有完整的辐亮度和噪声指标643.5 nm 那行的辐亮度全是空的。指标留空意味着它不参与常规的定量辐亮度流程角色更接近辅助或监视。数波段的时候把它单独拎出来不要混进定量波段里算。6. 什么时候该翻这份手册什么时候别指望它手册的定位其实很清楚它是一份参数和流程的说明书不是一本遥感入门教材也不是一套可以直接运行的代码库。该翻它的时候是你已经知道自己在做什么只需要确认具体数字和流程。比如要核对某个波段的中心波长和带宽查 B10 到底是 709 nm 还是 710 nm比如要确认某个辐亮度指标的量级判断一组观测数据是不是已经饱和比如要复现月球定标链条需要拿到那套按波段分组的系数和立体角常数。这些场景下手册的价值在于准确和完整翻到对应那一页就能拿到想要的东西。它还把七条月球成像判据按优先级排好了序做观测计划时可以直接照着筛窗口省掉自己重新推一遍的功夫。别指望它的时候是当你连基本概念都还没建立。手册里的公式是给已经懂的人看的它不会解释为什么反射率要对数化、为什么立体角要平方、为什么过采样校正要用反正切。它也明确写了有些工作还在下一步——过采样校正的精细模型、逐像素辐照度的解析式都属于待完善的部分不能当成已有结论直接用。还有一点波段表里那行宽通道的指标是空的如果你需要它的定量参数这份材料给不了。下面这张表可以贴在显示器边上判断起来快一些你的需求手册能不能解决查某个波段的中心波长、带宽能直接查表判断观测是否接近饱和能对照四个辐亮度等级评估噪声水平够不够用能看噪声等效辐亮度差排一年的月球成像计划能照七条判据筛复现月球定标计算链能系数和常数都给全了学遥感定标的基本原理不能得另找入门材料拿宽通道的定量参数不能指标栏是空的直接跑出一套完整定标软件不能部分模型仍待完善我第一次系统读这类手册的时候犯过一个挺低级的错把波段表里四个辐亮度数值当成重复标注直接跳过去看噪声那一列结果后面算量程余量的时候怎么算都对不上来回折腾了半天才反应过来是自己漏了信息。从那以后我每次读参数表都强制自己先把同一行里所有同名的数字按大小排一遍再往下看宁可慢一点。希望这点经验能帮到你。这份 GOCI-II 波段与月球定标的完整说明材料已经整理好了需要对照原始表格和系数逐条核对的可以直接取用。本文还有配套的精品资源点击获取 
长河君 · 2026/10/31 13:06:22 · 阅读 1.2k GOCI-II波段参数与月球定标计算链解析
简介《GOCI2 说明手册 波段信息说明》面向海洋遥感、水色反演与卫星辐射定标方向的学习者和工程人员集中梳理 GOCI-II 的波段设置与关键性能指标。文档以单份 PDF 形式提供压缩包约 2.93MB内含 1 个文件便于快速检索与对照查阅。内容围绕 GOCI 与 GOCI-II 的波段对比展开列出各波段中心波长、带宽、名义辐射亮度、最大海洋辐射亮度、阈值辐射亮度、最大云辐射亮度及 NEdLSNR 等参数并涉及快门罩、指向镜、焦平面组件、热沉等设备结构。文中还结合 GSICS/IVOS 月球校准研讨会资料介绍 GOCI-II 月球校准计划、MTF 测量计划与月球成像可行性评估有助于理解仪器观测能力、辐射定标思路和海洋环境监测应用。目前已有 1161 人学习。1. 一句话说清 GOCI-II地球同步轨道上那台“看海水的分光计”手机对着海面拍一张出来就是一片蓝。可海洋研究者盯着同一片蓝想知道的是水里叶绿素有多少、悬浮泥沙有多密、有没有藻类正在暴发。把“一片蓝”拆成十几个能互相加减的数字靠的就是 GOCI 这类海洋色度成像仪——它固定在静止轨道上一天之内反复回看同一片海把颜色切成十几个很窄的光谱窗口。GOCI-II 是它的第二代波段从 380 nm 一直排到 865 nm还配了一整套月球定标方案用来盯住自己的“眼睛”有没有变花。搞懂这套东西等于搞懂遥感数据里“颜色”是怎么一步步变成“浓度”的。2. 波段表到底在说什么从 380 nm 到 865 nm 的十二个窗口翻开这份手册最先撞进眼里的往往是那张密密麻麻的波段表。十二行窄波段加一行宽通道每行后面跟着中心波长、带宽、四个辐亮度数值和一个噪声指标。第一次看很容易当成一堆天书数字划过去但这张表其实是整台仪器的“能力清单”哪些东西它看得见多亮的东西它会饱和多小的差异它分得出来全在表里写着。2.1 中心波长与带宽每一行都是一扇窄窗先说最基础的两个词。中心波长是这个通道最敏感的位置带宽是它愿意接收的有效光谱宽度。打个比方中心波长是窗框中间那条线带宽是窗户有多宽。窗户越窄看东西越“挑”越能分辨相邻两条光谱线的差别窗户越宽一次进来的光越多信号更足但会把相邻成分混在一起。下面这张表按手册给的数据整理中心波长和带宽两列是固定的波段中心波长带宽这个窗口主要关注什么B1380 nm20 nm短波端给大气校正托底B2412 nm20 nm黄色物质、悬浮物的强吸收区B3443 nm20 nm叶绿素吸收峰附近B4490 nm20 nm叶绿素与水体透明度B5510 nm20 nm低浓度叶绿素的细微变化B6555 nm20 nm悬浮颗粒物的反射峰B7620 nm20 nm藻华、色素信号B8660 nm20 nm叶绿素荧光附近B9680 nm10 nm荧光峰带宽收窄到 10 nmB10709 nm10 nm藻华在水面形成的“红边”B11745 nm20 nm大气校正的近红外参考B12865 nm40 nm气溶胶、近红外参考宽通道643.5 nm483 nm跨可见光的大窗口指标栏留空看这张表有两个容易忽略的细节。第一B9 和 B10 的带宽只有 10 nm比其它通道窄一半——这两个位置对应的是很窄的光谱特征窗子开大了就会被旁边的信号冲淡。第二B12 反过来带宽拉到 40 nm因为近红外这段本来就没有尖锐特征宽一点能多收集能量、换来更好的信噪比。一窄一宽不是随意定的都是拿“要分辨什么”和“能收到多少光”这两件事反复权衡的结果。至于最后那行 643.5 nm、带宽 483 nm 的宽通道它的辐亮度指标在表里是空的。这一点很值得注意指标留空通常意味着这个通道不参与常规的定量辐亮度流程它的角色更接近辅助或监视而不是出定量产品。2.2 四个辐亮度数值一台仪器的“量程刻度”波段表里最让人犯迷糊的是一行里出现四个名字都带“辐亮度”的数字。它们都是同一个物理量单位是 W·m⁻²·μm⁻¹·sr⁻¹差别只在于描述“在什么条件下”名义辐亮度正常海面情况下大概会遇到的亮度水平是设计时锚定的基准工作点。最大海洋辐射亮度海面在明亮条件下可能达到的上限超过它海洋目标的观测就开始失真。阈值辐亮度一条警戒线用来判断当前信号是否还在可靠范围内。最大云辐射亮度云顶反射的亮度能到多高。这个值远大于前三个因为白云反射的能量要比海水高一个量级。以 B1 为例名义辐亮度是 93最大海洋辐射亮度 139.5阈值辐亮度 143.16最大云辐射亮度 634.4。四个数字一摆量程的层次立刻清楚了正常观测在 93 附近海面再亮也就顶到 140 上下而云的亮度能冲到 630 以上。换句话说设计者心里同时装着“海”和“云”两种极端场景而不是只盯着海。这里有个很实用的读法名义辐亮度和最大海洋辐射亮度之间的差距反映的是这台仪器对“海面变亮”这件事留了多大余量最大云辐射亮度和它们之间的巨大落差则是留给云检测和大气校正的空间。如果只记住一个数就失去了这张表一半的信息。2.3 噪声指标能“看见”和能“分辨”是两件事表里最后一列是 NEdL中文常译成噪声等效辐亮度差。名字拗口意思其实很朴素仪器自身抖动产生的那点“底噪”折算成辐亮度是多少。这个数值越小说明它越能分辨两个靠得很近的亮度水平。B1 的 NEdL 是 0.09B12 是 0.015。单看数字没感觉跟同行的名义辐亮度放一起就有味道了B1 的名义辐亮度是 93底噪 0.09相当于在 93 这个量级上还能分辨千分之一左右的差别。这就是为什么波段表要同时给出名义辐亮度和噪声——光有量程只是说明“能看见多亮的东西”加上噪声才说明“能看清多细的差别”。一台量程很大但噪声也大的仪器拍出来的图可能很亮但一片糊。提示读波段表时先把中心波长和带宽扫一遍建立整体印象再回头只看你关心的那两三个波段对应的辐亮度和噪声最后才通读全表。反过来从第一个数字开始逐格读很容易读到第四行就失去方向。3. 为什么要用月球做定标一架不会自动变旧的“天上标尺”波段表解决的是“设计上能看多好”可仪器上了天实际表现会一点点漂。传感器会老化滤光片会变镜头会脏电子学增益也会飘。要判断数据还能不能信就必须有个稳定的参照物反复比对。问题是地面上那把“标准尺”没法搬到天上天上的东西自己也会变。3.1 星上定标器的难处参照物自己也会脏常规做法是在卫星上装一块漫反射板让它把太阳光均匀散射进镜头当作已知亮度的参照。听起来很完美但漫反射板暴露在太空环境里会被紫外辐照、污染、粒子轰击慢慢改变反射特性。于是出现一个尴尬的循环你想用板子校准仪器可板子自己也在退化你并不知道它退化了多少。所以工程师又给它配了专门盯板的装置——手册里提到的漫反射板老化监视装置就是干这个的。但这只是缓解不是根治。参照物自己会变旧这个根本矛盾始终存在。月球的好处恰好在这里。它没有大气表面在人类观测的时间尺度上几乎不变反射特性稳定得可以当作“宇宙标准灯”。而且它就在天上任何一颗能转动的卫星都有机会把它框进视场。用月球做定标相当于绕开了星上参照物会退化的死结。3.2 七条判据先看能不能拍到再看拍得好不好想法很美好但月球不是想拍就能拍。手册里列了一张分优先级的可行性判据读下来会发现它的排序逻辑非常“工程师”先保证能拍到再保证拍得好。优先级 1月球必须落在视场里。视场在东西方向是 ±8.7 度南北方向是 ±11.0 度地球圆盘正好对应 ±8.7 度。这条最硬进不了视场后面全都白搭。手册还特别注明不打算为此做整星姿态机动也就是说能不能拍到完全取决于自然几何不去为它额外调整姿态。优先级 2要在标称工作时段内。每天的运行窗口是 05:15 到 19:45出了这个区间仪器状态就不在标称条件下。优先级 3月球要有足够的亮面。被照亮的部分要超过一半满月当然最好。月相太小信号弱噪声占比就上来了。优先级 4不能和其它载荷抢时间。同一颗卫星上还有别的仪器在工作月球成像的时段不能完全和它们重叠。优先级 5单次可连续观测时间要够长。判据写的是不少于 10 分钟待定。时间太短积分不充分。优先级 6避开地球杂散光。地球本身很亮反射光会从视场边缘漏进来。手册把距地球圆盘 1 度以内的区域划成潜在污染区月球图像要尽量避开。优先级 7相位角要落在合适区间。这一条跟定标模型直接相关。月球反射率随相位角变化而模型本身在这个维度上有约 1% 的相对不确定度所以要挑一个模型算得最准的相位区间。七条判据摆在一起能看出一个很清楚的态度月球定标不是“有机会就拍”而是“在满足一堆约束的前提下挑一个最优窗口”。每一条都在把可选项往下砍砍到最后剩下的窗口其实相当窄。3.3 视场、指向与那堆机械件手册里那张仪器结构图上标了一串部件快门盖、支撑结构、卫星接口、四个双脚架、滤光轮、带热分流器的焦平面组件、指向镜、遮光罩、望远镜机构。还有一行小字注明某个辐射增强系统并不属于这台仪器——这种备注在手册里很常见作用是划清边界免得读者把别人的指标算到自己头上。跟月球成像关系最直接的是指向镜和遮光罩。指向镜决定镜头能不能转向月球所在的方向遮光罩负责挡住视场外的强光。四个双脚架把整机固定在卫星接口上热分流器负责把焦平面的热量导走——温度一飘暗电流就飘定标结果也跟着飘。看起来都是机械细节实际上每一条都在影响最后那个数字的可靠性。视场的设计也很有意思地球圆盘占 8.7 度仪器视场留 1.1 度给月球和恒星再留 0.5 度应对指向偏差。这个余量分配说明月球成像并不是“顺便蹭一下”而是从设计阶段就预留了位置。3.4 一张排期表透露的现实约束手册里附了一张 2019 年的逐月月球成像排期。数据很直白月份可观测天数推荐观测时长1 月2 天约 30 分钟2 月3 天约 30 分钟3 月1 天约 38 分钟4 月1 天约 46 分钟5 月1 天约 29 分钟6 月2 天约 30 分钟7 月2 天约 39 分钟8 月2 天约 43 分钟9 月1 天约 30 分钟10 月2 天约 11 分钟11 月4 天约 12 分钟12 月3 天约 44 分钟一个月能拍的日子只有一到四天可用的连续时间从 11 分钟到 46 分钟不等。10 月和 11 月特别短只有十一二分钟。这意味着整年的定标机会其实非常有限每一次拍摄都得当成正式任务来准备不能指望“这次没拍好下个月补”。这也解释了为什么判据要列得那么细——机会成本太高必须提前筛掉所有不合格的窗口。4. 从一张月球照片到一个数字ROLO 模型这条计算链拍到月球只是开始。真正的定标要回答一个更精确的问题仪器测到的月球有多亮和月球“应该有多亮”差了多少。这中间隔着一整条计算链每一环都需要单独交代清楚。4.1 第一步算出“月亮本来应该有多亮”月球不是一盏亮度恒定的灯。它的反射率会随相位角变化——满月时最亮月牙时最暗而且这个变化不是线性的。同时观测几何里月面的纬度、经度太阳照射的经度都会影响你看到的那部分月面反射多少光。国际上有一套长期维护的月球反射率模型手册里用的就是它。核心公式是圆盘等效反射率的对数值由几个多项式项和指数项叠加而成。系数按波段各自一套一共十二组。挑其中一组写成代码结构是这样的import numpy as np # 计算某个波段的圆盘等效反射率 Ak # g : 绝对相位角单位度 # theta : 观测点对应的月面纬度 # phi : 观测点对应的月面经度 # Phi : 太阳的月面经度 # coef : 该波段的一组系数含 a[4]、b[3]、c[4]、d[3]、p[4] def rolo_disk_reflectance(g, theta, phi, Phi, coef): a, b, c, d, p coef[a], coef[b], coef[c], coef[d], coef[p] g_rad np.deg2rad(g) # 相位角的多项式部分对应 a0 a1*g a2*g^2 a3*g^3 s sum(a[i] * g_rad ** i for i in range(4)) # 太阳经度的奇次幂部分Phi 的 1、3、5 次方 s sum(b[j] * Phi ** (2 * j 1) for j in range(3)) # 经纬度交叉项四个全局系数 s c[0] * theta c[1] * phi c[2] * Phi * theta c[3] * Phi * phi # 两个指数衰减项描述小相位角附近的陡峭变化 s d[0] * np.exp(-g_rad / p[0]) s d[1] * np.exp(-g_rad / p[1]) # 余弦项对应 cos(g - p3) / p4 s d[2] * np.cos(g_rad - p[2]) / p[3] return np.exp(s) # 反射率本身是对数结果的指数这段代码有两处值得单独说。第一np.exp(s)这一步说明原始公式给的是对数形式系数都是在对数域里拟合出来的算完必须取指数才能得到反射率。第二指数项里的p[0]和p[1]是衰减尺度它们决定了相位角从 0 度往上走时反射率下降得有多快。这两个参数在不同波段之间差异不小不能混用。系数本身按波段给全每个波段十四五个数字看起来庞杂但结构是统一的a管相位角的主趋势b管太阳经度c管交叉调制d和p管小角度附近的修正。只要结构对上了换波段只需要换系数。4.2 第二步把反射率积成辐照度有了反射率还不够仪器测的是辐照度。要把前者变成后者得走一遍积分积分月球参考辐照度 立体角 / π× 太阳光谱辐照度 × 仪器光谱响应 × 圆盘等效反射率的积分 ÷ 光谱响应的积分这条式子里太阳光谱辐照度是参考值仪器光谱响应是每个波段自己的响应曲线反射率就是上一步算出来的结果。除以光谱响应的积分本质是在做归一化——不同波段的响应宽度不一样不归一化就没法比较。手册里给了一个固定常数月球对应的立体角是 6.4177×10⁻⁵ 球面度。这个值在整个计算里是常量不需要每次重算。算完参考辐照度还要做距离校正。地球绕太阳的轨道是椭圆月球绕地球的轨道也不是正圆两者之间的距离每天都在变接收到的能量自然也跟着变。所以有一个距离校正因子把测量值折算到统一距离下距离校正后的测量辐照度 参考辐照度 × 距离校正因子这一步很容易被忽略因为它的影响不像噪声那么显眼但它是系统性的。少了这一步不同月份拍到的数据就没法直接对比。4.3 第三步把图像里的像素加回去前面两步算的都是“应该值”现在轮到“实测值”。月球在图像上占一片像素每个像素都有自己的立体角。把每个像素的辐亮度乘以它对应的立体角再全部加起来就得到了积分月球辐照度# 从月球图像反算积分月球辐照度 # 前提图像已经做过过采样校正 GSD 250.0 # 星下点地面采样距离单位 m手册给出优于 250 m H 35786.0e3 # 卫星轨道高度单位 m OMEGA_PIX (GSD / H) ** 2 # 单个像元张开的立体角单位球面度 def integrated_lunar_irradiance(radiance_img, moon_mask): # radiance_img : 二维数组每个像素的月球辐亮度 # moon_mask : 布尔数组标出落在月面上的像素 return float(np.sum(radiance_img[moon_mask] * OMEGA_PIX))这段代码的关键在OMEGA_PIX这个量。它是地面采样距离除以轨道高度再平方单位是球面度。地面采样距离越小单个像素张开的立体角越小同样的辐亮度贡献的能量也越小。手册里给出的地面采样距离优于 250 米所以这个值是个上限实际计算时如果手上有更精确的数字应该替换掉。另外要注意moon_mask。月球在图像上不是一个完整的圆盘边缘会因为曝光和指向偏差变得模糊所以哪些像素算“月面”需要单独判定。判得太松会把背景天光算进去判得太紧会漏掉边缘两种做法都会让结果偏低或偏高。4.4 第四步过采样校正、比值因子与老化因子到这一步还有两个坑。第一个坑是过采样。仪器在积分时间内不是静止的月球在动卫星也在动结果就是月面在图像上被“拖”出了一段。手册算过一笔账积分时间 1 秒月球大约移动 10 个像素最坏情况积分 4.5 秒能移动 50 个像素。卫星平均轨道角速度约 72.9 微弧度每秒单个像素视场约 7 微弧度两个数字一除正好对上。校正的办法是拿月球的真实角直径和图像上测到的角直径做比较过采样校正系数 arctan(3476.4 km / 月地距离) ÷ 图像上测得的月球角直径3476.4 公里是月球直径这个值是固定的。分母是图像上实测的角尺寸两者一比就知道图像被拉伸了多少倍再按比例修正回去。手册也坦白说如果判定确实需要校正更精细的模型要把月球运动、卫星运动和指向稳定性一起考虑那是下一步的工作。第二个坑是不能逐像素重算增益。这一点常被误解。月球在单张图像里只覆盖了一部分探测器像素而参考模型给的是整个月球的积分辐照度。用部分像素去反推每一个像素的绝对增益在数学上是不成立的。所以手册没有走这条路而是定义了比值因子比值因子 仪器测得的积分月球辐照度 ÷ 模型给出的参考辐照度这个比值反映的是整机层面的辐射性能变化不落到单个像素。把比值取倒数就得到月球老化因子月球老化因子 参考辐照度 ÷ 实测辐照度老化因子配合像素响应非均匀性就能评估每个像素的辐射性能变化。注意这里的逻辑顺序先用整体比值判断“仪器整体衰减了多少”再用非均匀性把整体衰减分摊到各个像素。这个两步走的做法正是被“月球只覆盖部分像素”这个技术限制逼出来的。5. 常见误区与易混概念这几个词最容易被混着用读到这儿概念大体串起来了但有几个地方特别容易在脑子里糊成一团。下面五条是初学者最常踩的坑每一条都值得停下来对一遍。5.1 误区一波段越多、带宽越窄就代表越先进易混点看到十二个窄波段加一个宽通道直觉上会认为波段数量多、带宽窄就是技术更强。为什么混带宽窄确实能提高光谱分辨能力波段多确实能覆盖更多现象这两件事单独看都对于是很容易直接画等号。怎么区分带宽是拿信噪比换来的。窗子开窄进来的人光少了噪声占比就上去了。B9 和 B10 敢用 10 nm是因为那两个位置有明确的光谱特征值得为之付出代价B12 用 40 nm是因为近红外没有尖锐特征宽一点反而更划算。所以判断一台仪器不能数波段个数要看每个波段的设计意图和它配套的噪声指标是否匹配。5.2 误区二把四个辐亮度数值当成同一个东西的不同叫法易混点名义辐亮度、最大海洋辐射亮度、阈值辐亮度、最大云辐射亮度四个名字都带“辐亮度”单位也一样看起来像是重复标注。为什么混它们确实是同一个物理量只是对应的观测条件不同名字又长得很像扫一眼就滑过去了。怎么区分把四个数按大小排一下层次立刻出来。以 B1 为例名义 93、最大海洋 139.5、阈值 143.16、最大云 634.4。名义是日常工作的锚点最大海洋是海面变亮的余量上限阈值是警戒线最大云是另一个量级的极端场景。四个数字描述的是一台仪器的完整量程剖面缺一个就少一块信息。5.3 误区三以为月球定标能给每个像素重新算增益易混点“定标”这个词天然带着“把每个像素都校准一遍”的联想。为什么混很多地面相机的平场校正确实是逐像素做的习惯成自然就会以为天上的定标也一样。怎么区分月球在单张图像里只覆盖部分探测器像素参考模型给的又是整月球的积分值用局部推全体在数学上不成立。所以手册走的是“先算整体比值因子再用非均匀性分摊”的路线。逐像素的绝对增益更新月球定标做不到。5.4 误区四把月球当成亮度恒定的标准灯易混点既然说月球表面稳定就容易进一步认为它每次看起来都一样亮。为什么混“稳定”指的是月面反射特性在长时间尺度上不变不是指观测到的亮度不变。怎么区分亮度随相位角剧烈变化满月和上下弦能差出好几倍。所以才有优先级 7 那条判据——要把相位角限制在模型算得最准的区间手册给的不确定度大约是 1%。除此之外还有距离校正、过采样校正。月球是稳定的参照物不是省事的参照物。5.5 误区五把 643.5 nm 那行宽通道和十二个窄波段当成一回事易混点表格里并排列着很容易顺口读成“十三个波段”。为什么混格式统一位置相邻视觉上就是同一类东西。怎么区分看指标栏。十二个窄波段每行都有完整的辐亮度和噪声指标643.5 nm 那行的辐亮度全是空的。指标留空意味着它不参与常规的定量辐亮度流程角色更接近辅助或监视。数波段的时候把它单独拎出来不要混进定量波段里算。6. 什么时候该翻这份手册什么时候别指望它手册的定位其实很清楚它是一份参数和流程的说明书不是一本遥感入门教材也不是一套可以直接运行的代码库。该翻它的时候是你已经知道自己在做什么只需要确认具体数字和流程。比如要核对某个波段的中心波长和带宽查 B10 到底是 709 nm 还是 710 nm比如要确认某个辐亮度指标的量级判断一组观测数据是不是已经饱和比如要复现月球定标链条需要拿到那套按波段分组的系数和立体角常数。这些场景下手册的价值在于准确和完整翻到对应那一页就能拿到想要的东西。它还把七条月球成像判据按优先级排好了序做观测计划时可以直接照着筛窗口省掉自己重新推一遍的功夫。别指望它的时候是当你连基本概念都还没建立。手册里的公式是给已经懂的人看的它不会解释为什么反射率要对数化、为什么立体角要平方、为什么过采样校正要用反正切。它也明确写了有些工作还在下一步——过采样校正的精细模型、逐像素辐照度的解析式都属于待完善的部分不能当成已有结论直接用。还有一点波段表里那行宽通道的指标是空的如果你需要它的定量参数这份材料给不了。下面这张表可以贴在显示器边上判断起来快一些你的需求手册能不能解决查某个波段的中心波长、带宽能直接查表判断观测是否接近饱和能对照四个辐亮度等级评估噪声水平够不够用能看噪声等效辐亮度差排一年的月球成像计划能照七条判据筛复现月球定标计算链能系数和常数都给全了学遥感定标的基本原理不能得另找入门材料拿宽通道的定量参数不能指标栏是空的直接跑出一套完整定标软件不能部分模型仍待完善我第一次系统读这类手册的时候犯过一个挺低级的错把波段表里四个辐亮度数值当成重复标注直接跳过去看噪声那一列结果后面算量程余量的时候怎么算都对不上来回折腾了半天才反应过来是自己漏了信息。从那以后我每次读参数表都强制自己先把同一行里所有同名的数字按大小排一遍再往下看宁可慢一点。希望这点经验能帮到你。这份 GOCI-II 波段与月球定标的完整说明材料已经整理好了需要对照原始表格和系数逐条核对的可以直接取用。本文还有配套的精品资源点击获取