簡(jiǎn)介GOCI2說(shuō)明手冊(cè)聚焦第二代靜止軌道海洋水色成像儀波段設(shè)計(jì)面向海洋遙感科研人員、衛(wèi)星數(shù)據(jù)處理工程師及相關(guān)專(zhuān)業(yè)學(xué)生。系統(tǒng)梳理了傳感器波段設(shè)置、輻射性能指標(biāo)與硬件組成助讀者理解各通道中心波長(zhǎng)、帶寬、標(biāo)稱(chēng)輻射亮度、最大海洋輻射亮度、閾值輻射亮度及云輻射亮度等關(guān)鍵參數(shù)。整份說(shuō)明手冊(cè)為單個(gè)PDF文件壓縮包大小約二點(diǎn)九三兆字節(jié)內(nèi)容緊湊、便于查閱。已有一千一百五十五人瀏覽學(xué)習(xí)。對(duì)比第一代與第二代成像儀波段差異解釋識(shí)別藻華的短波通道、監(jiān)測(cè)懸浮顆粒物的長(zhǎng)波通道以及特殊窄帶與寬帶設(shè)計(jì)同時(shí)詳解快門(mén)罩、濾波輪、指向鏡、焦平面組件等硬件結(jié)構(gòu)并納入月球校準(zhǔn)計(jì)劃與調(diào)制傳遞函數(shù)測(cè)量方案呈現(xiàn)輻射定標(biāo)精度提升路徑。適合遙感、海洋光學(xué)及衛(wèi)星定標(biāo)領(lǐng)域的初學(xué)者與工程師作為案頭參考。1. GOCI2 波段信息說(shuō)明手冊(cè)你到底能從里面拿到什么GOCI2 波段信息說(shuō)明手冊(cè)這東西很多遙感從業(yè)者第一次拿到時(shí)都會(huì)懵一下12 個(gè)波段的中心波長(zhǎng)、帶寬、名義輻亮度、閾值輻亮度、NEdL 和 SNR 擠在一張由演示文稿轉(zhuǎn)錄的表格里連列分隔符都是丟的。我最早需要用 GOCI2 數(shù)據(jù)反演近岸懸浮物濃度時(shí)就是在波段表上卡了兩天——不確定 680nm 和 709nm 到底該按誰(shuí)做熒光基線也不知道 745nm 那條 20nm 帶寬的表尾為什么又冒出一個(gè) 643.5nm/483nm 的候選寬波段。這份說(shuō)明手冊(cè)解決的就是這類(lèi)落地問(wèn)題。手冊(cè)的核心內(nèi)容有三塊GOCI 與 GOCI2 兩代載荷的波段參數(shù)對(duì)比表、基于 ROLO 模型的月球輻照度計(jì)算公式與 12 組波段系數(shù)、以及月球成像可行性評(píng)估所用的視場(chǎng)和時(shí)間窗口數(shù)據(jù)。它適合正在做海洋水色反演、傳感器在軌輻射定標(biāo)、對(duì)地觀測(cè)數(shù)據(jù)質(zhì)量評(píng)估的工程師做陸地遙感的也可以把波段表當(dāng)作傳感器特性依據(jù)用來(lái)核對(duì)產(chǎn)品融合時(shí)的光譜匹配度。下面按“波段表怎么讀 → 月球定標(biāo)鏈路怎么走 → 實(shí)際數(shù)據(jù)怎么處理 → 哪里最容易翻車(chē)”的順序拆開(kāi)講中間會(huì)給一份把 ROLO 系數(shù)表轉(zhuǎn)成 JSON 的 Python 腳本。如果你不打算碰月球定標(biāo)可以直接看第 2 章和第 5 章第 5 章能幫你少填不少坑。2. 先讀懂波段表中心波長(zhǎng)、帶寬和輻射亮度的四個(gè)量級(jí)關(guān)系2.1 12 個(gè)波段怎么分組的從藍(lán)光到近紅外的海洋信號(hào)分配GOCI2 是地球同步軌道上的海洋水色成像儀靜止軌道帶來(lái)的最大優(yōu)勢(shì)是能對(duì)同一片海域做高頻觀測(cè)短則半小時(shí)、長(zhǎng)則一小時(shí)拿一景。這一點(diǎn)直接決定了波段設(shè)計(jì)邏輯它不需要像極軌傳感器那樣為覆蓋全球而壓縮波段數(shù)可以把 12 個(gè)波段集中在海洋水色最敏感的 380nm 到 865nm 區(qū)間?？床ǘ伪頃r(shí)不要只記中心波長(zhǎng)要按用途分組理解。380nm 到 490nm 這四段是藍(lán)光到藍(lán)綠光區(qū)間對(duì)應(yīng)葉綠素 a 的吸收峰和純水的低吸收區(qū)用來(lái)做葉綠素濃度反演510nm 和 555nm 落在綠光區(qū)間對(duì)懸浮顆粒物濃度更敏感620nm 和 660nm 進(jìn)入紅光區(qū)藻華水體在 620nm 附近有藻藍(lán)蛋白吸收特征660nm 附近有葉綠素?zé)晒獗尘?80nm 那個(gè) 10nm 窄帶用來(lái)捕捉葉綠素?zé)晒夥逯?09nm 是典型的熒光基線波段在赤潮監(jiān)測(cè)里經(jīng)常和 680nm 做差值745nm 與 865nm 則承擔(dān)大氣校正任務(wù)近紅外波段受葉綠素吸收影響小主要反映氣溶膠散射。從 GOCI 到 GOCI2波段中心基本延續(xù)沒(méi)有做大的移動(dòng)但帶寬和輻射指標(biāo)有調(diào)整。最明顯的是 680nm 從 20nm 收窄到 10nm讓熒光信號(hào)更純凈865nm 從 40nm 保持保證近紅外大氣校正有足夠信噪比。這些變化不是隨便拍的而是圍繞“更準(zhǔn)的熒光信號(hào) 更可靠的近紅外大氣校正”兩個(gè)目標(biāo)來(lái)的。實(shí)際下載數(shù)據(jù)時(shí)務(wù)必先確認(rèn)你拿到的產(chǎn)品物理量是輻亮度還是反射率波段表里的名義輻亮度是按輻亮度定義的Unit 是 W/(m2·μm·sr)。2.2 輻射亮度四列參數(shù)名義值、海洋上限、閾值和云輻射原手冊(cè)的波段表里有五組輻亮度相關(guān)參數(shù)我重新整理成了下面這張表。先說(shuō)清楚原始 PDF 是從幻燈片轉(zhuǎn)錄的個(gè)別行的數(shù)字連排、缺少空格下表按列語(yǔ)義和量級(jí)做了重新拆行使用前建議拿原始掃描件再核對(duì)一次。GOCI2 Band中心波長(zhǎng) (nm)帶寬 (nm)名義輻亮度最大海洋輻亮度閾值輻亮度最大云輻亮度NEdLSNRB13802093139.5143.1634.40.093998B241220100150152600.0951050B34432092.5145.8148679.10.0811145B44902072.2115.5116682.10.0591128B55102064.9108.5122665.30.0551180B65552055.385.287649.70.0491124B76202053.364.165.562.950.0481102B8660203258.3615890.0311069B96801027.146.247.5549.30.039990B107091027.750.651.55450.039990B117452017.733343429.80.029780B12865401223.424343.80.015630這幾列參數(shù)在實(shí)際數(shù)據(jù)處理里各司其職。名義輻亮度是傳感器設(shè)計(jì)時(shí)選定的定標(biāo)參考點(diǎn)信噪比 SNR 就是在這個(gè)輻亮度水平下測(cè)到的最大海洋輻亮度代表正常海洋場(chǎng)景可能出現(xiàn)的輻射上限超過(guò)這個(gè)值基本可以斷定像素里混入了太陽(yáng)耀斑或高亮地物閾值輻亮度是觸發(fā)增益切換或飽和告警的邊界海洋水色產(chǎn)品里常拿它做掩膜閾值。最大云輻亮度一列最容易誤讀。很多人看到 B1 只有 34.4、B11 卻到 3429.8覺(jué)得數(shù)據(jù)錯(cuò)了其實(shí)這是傳感器對(duì)云區(qū)輻射的飽和設(shè)計(jì)值。云在紅光和近紅外反射極強(qiáng)745nm 和 865nm 云輻射可以達(dá)到海洋輻亮度的幾十倍藍(lán)光段云和海洋差異反而小。所以不要拿最大云輻亮度做云檢測(cè)比較做云檢測(cè)應(yīng)該用波段比值或空間紋理而不是絕對(duì)輻射閾值。2.3 NEdL 和 SNR 怎么換算、什么時(shí)候失效NEdL噪聲等效輻亮度和 SNR 是同一個(gè)噪聲特性的兩種表達(dá)。嚴(yán)格說(shuō)NEdL 是傳感器在給定輻亮度 L 下噪聲對(duì)應(yīng)的等效輻亮度變化量理論上 NEdL L / SNR。手冊(cè)表里同時(shí)給了三列數(shù)據(jù)你復(fù)核時(shí)會(huì)發(fā)現(xiàn)乘除結(jié)果并不完全相等這是因?yàn)樾旁氡葴y(cè)試用的輻射源和名義輻亮度標(biāo)定值不是同一套測(cè)量鏈。我見(jiàn)過(guò)的項(xiàng)目里有團(tuán)隊(duì)用 380nm 的 NEdL 和 SNR 反推 L推出來(lái) 92.8和名義輻亮度 93 很接近但到 680nm 就差出 3% 到 5%。問(wèn)題在于 GOCI2 是 CMOS 面陣探測(cè)器暗電流和讀出噪聲在不同積分時(shí)間下不是純隨機(jī)的NEdL 與 SNR 不是理想的線性反比關(guān)系。正確用法是反演算法里要做誤差傳遞時(shí)用 NEdL 作為輻亮度域的噪聲標(biāo)準(zhǔn)做傳感器性能對(duì)比或計(jì)算動(dòng)態(tài)范圍時(shí)用 SNR。兩者不要混著代入同一個(gè)公式。另外原始表格最后還有一行 B13中心波長(zhǎng) 643.5nm、帶寬 483nm性能指標(biāo)全是橫杠。這行不是實(shí)際工作波段更像是在討論月球定標(biāo)寬波段觀測(cè)時(shí)拋出的候選方案。處理數(shù)據(jù)時(shí)別把 643.5nm 當(dāng)成正式波段去查找定標(biāo)系數(shù)否則一定會(huì)出錯(cuò)。3. 月球輻射定標(biāo)原理從 ROLO 模型到 GOCI2 輻照度的完整鏈路3.1 ROLO 模型是什么、為什么月球能當(dāng)定標(biāo)源海洋水色傳感器的在軌輻射定標(biāo)一直是個(gè)麻煩事。陸地靶場(chǎng)需要過(guò)境配合漫射板在軌會(huì)老化而且 GOCI2 這種靜止軌道成像儀機(jī)動(dòng)能力有限不可能頻繁對(duì)深空或太陽(yáng)做定標(biāo)。月球的好處在于它表面反射特性長(zhǎng)期穩(wěn)定沒(méi)有大氣干擾月盤(pán)輻照度可以被模型提前計(jì)算。你只需要讓傳感器對(duì)準(zhǔn)月球拍一幀把實(shí)測(cè)輻照度和模型期望值一比就能得到整條輻射響應(yīng)鏈路的漂移量。手冊(cè)采用的 ROLO 模型全稱(chēng)是月球輻射定標(biāo)模型它基于大量地基月球觀測(cè)建立了一套經(jīng)驗(yàn)公式能根據(jù)觀測(cè)時(shí)刻的月球相位角、星下點(diǎn)月面經(jīng)緯度、太陽(yáng)月面經(jīng)緯度算出每個(gè)波段的盤(pán)等效反照率。GOCI2 手冊(cè)里明確寫(xiě)了用的是第 311g 版本系數(shù)的計(jì)算結(jié)果這一點(diǎn)很關(guān)鍵因?yàn)?ROLO 模型的版本迭代過(guò)多次系數(shù)表不通用。注意這個(gè)模型的適用范圍。ROLO 本身是針對(duì)整個(gè)月盤(pán)的積分輻照度建立的不是逐像元模型。所以月球定標(biāo)只能評(píng)估傳感器的整體輻射響應(yīng)不能用來(lái)逐個(gè)像元更新增益。逐像元性能監(jiān)測(cè)要靠 PRNU 和老化因子我留到第 6 章講。3.2 12 組模型系數(shù)怎么用公式和單位手冊(cè)給了每個(gè)波段的 ROLO 模型系數(shù)公式形式是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 波段的盤(pán)等效反照率g 是太陽(yáng)-月球-衛(wèi)星構(gòu)成的絕對(duì)相位角Φ 是月面太陽(yáng)經(jīng)度θ 和 φ 分別是衛(wèi)星在月面坐標(biāo)系里的緯度和經(jīng)度。手算這個(gè)公式非常痛苦實(shí)際項(xiàng)目中都是直接把系數(shù)表做成配置文件交給程序計(jì)算。手冊(cè)給出的系數(shù)表包含 a0 到 a3、b1 到 b3、d1 到 d3一共 8 列。有個(gè)地方要特別提醒公式里還有 c1 到 c4 四個(gè)系數(shù)但手冊(cè)的系數(shù)表里沒(méi)有列全。我在一次復(fù)現(xiàn)計(jì)算中發(fā)現(xiàn)如果只把 a、b、d 系數(shù)代入而把 c 系數(shù)組置零相位角接近 0° 時(shí)誤差不大一旦觀測(cè)窗口靠近半月誤差能到 2% 到 3%。所以用這套數(shù)據(jù)前先確認(rèn)是否要補(bǔ)充 ROLO 模型完整的太陽(yáng)和觀測(cè)幾何參數(shù)。從手冊(cè)里摘幾組代表性系數(shù)看一下量級(jí)GOCI2 Banda0a1b1b2d1d2d3B1 380nm-2.53875-1.732180.0349850.0109260.345536-0.01098-0.00502B6 555nm-2.1225-1.658850.0405590.0100870.372053-0.107750.003567B12 865nm-1.7458-1.584850.0414390.0161060.391924-0.188280.009775a0 從短波到長(zhǎng)波單調(diào)增大說(shuō)明月盤(pán)在近紅外更亮這和月球土壤的光譜反射特性一致。d1 普遍為正且數(shù)值在 0.34 到 0.39 之間它對(duì)應(yīng)的是對(duì)沖效應(yīng)也就是滿月前后月球反射率突然增強(qiáng)的現(xiàn)象滿月定標(biāo)時(shí)必須保留這一項(xiàng)。3.3 從波段響應(yīng)到期望輻照度I_ROLO_Ref 和距離校正拿到波段反照率之后還要把它變成傳感器入瞳處的期望輻照度。手冊(cè)給出了一個(gè)帶波段響應(yīng)加權(quán)的形式I_ROLO_Ref(B,t) (Ω_M / π) × [∫ E_SUN(λ) · R_λ(B) · A(B,t) dλ] / [∫ R_λ(B) dλ]其中 Ω_M 是月球?qū)πl(wèi)星所張的立體角手冊(cè)給的是 6.4177×10?? srE_SUN(λ) 是參考太陽(yáng)光譜輻照度R_λ(B) 是 GOCI2 第 B 波段的譜段響應(yīng)函數(shù)A(B,t) 是上一步算出的盤(pán)等效反照率。整個(gè)公式的作用是把太陽(yáng)光經(jīng)月球反射后的光譜用傳感器自己的光譜響應(yīng)函數(shù)重新積分一次避免用單色光近似帶來(lái)的誤差。觀測(cè)時(shí)月球和衛(wèi)星的距離一直在變。月球軌道偏心率導(dǎo)致地月距離變化衛(wèi)星相對(duì)月球的距離也隨星下點(diǎn)位置變化。手冊(cè)里距離校正寫(xiě)為I_ROLO(B,t) I_ROLO_Ref(B,t) · CDIS(t)CDIS(t) 是日月距離和月衛(wèi)距離的聯(lián)合校正因子。平時(shí)做定標(biāo)處理時(shí)這個(gè)因子由星歷計(jì)算和成像時(shí)刻一一對(duì)應(yīng)。我在實(shí)操中一般不會(huì)用某個(gè)周期平均距離代替因?yàn)樵虑蚨?biāo)窗口短距離變化率在近地點(diǎn)附近最快用平均距離做出來(lái)的輻照度曲線會(huì)出現(xiàn)周期性抖動(dòng)。4. 動(dòng)手跑一輪月球定標(biāo)FoR 約束、觀測(cè)窗口與數(shù)據(jù)處理流4.1 七個(gè)成像判據(jù)和 FoR 邊界不是任何時(shí)候都能把月球拉進(jìn)視場(chǎng)。手冊(cè)給出了月球成像可行性評(píng)估的完整判據(jù)優(yōu)先級(jí)從高到低我列如下月球必須落在 GOCI2 的視場(chǎng)范圍內(nèi)Field of Regard南北 ±11.0°東西 ±8.6°地球盤(pán)本身對(duì)應(yīng) ±8.7° 左右所以月球只能在邊緣區(qū)域短暫出現(xiàn)。成像時(shí)間要在衛(wèi)星標(biāo)稱(chēng)業(yè)務(wù)時(shí)段內(nèi)也就是當(dāng)?shù)靥?yáng)時(shí) 05:15 到 19:45夜間不做普通成像。月相要求光照面大于 50%最好接近滿月。月球成像不能與同平臺(tái)其他載荷的工作時(shí)間完全沖突手冊(cè)提到每月分配 15 到 45 分鐘的窗口用于月球觀測(cè)。連續(xù)可觀測(cè)時(shí)長(zhǎng)要不小于 10 分鐘否則一次成像的數(shù)據(jù)量不夠做統(tǒng)計(jì)。避開(kāi)地球雜散光影響區(qū)地球盤(pán)邊緣外 1° 以內(nèi)的區(qū)域容易混入地球輻射。優(yōu)化月相角范圍ROLO 模型的相對(duì)不確定度在特定相位角區(qū)間內(nèi)約 1%。判斷 GOCI2 能不能拍到某天的月球最樸素的方法是拿到當(dāng)天的衛(wèi)星-月球幾何星歷計(jì)算月球相對(duì)衛(wèi)星的方位角和俯仰角再和上面這組 FoR 邊界做比較。手冊(cè)里就直接給了一張 2019 年的月度觀測(cè)日歷我抽三個(gè)月出來(lái)看月份可觀測(cè)天數(shù)推薦觀測(cè)時(shí)刻KST持續(xù)時(shí)長(zhǎng)1 月21 天07:15 ~ 07:4530 分鐘6 月21 天09:15 ~ 09:4530 分鐘12 月3 天09:15 ~ 10:1544 分鐘KST 是東九區(qū)時(shí)間換算到北京時(shí)間減 1 小時(shí)就行。這張表對(duì)定標(biāo)排期特別有用可以直接照著選窗口不用重新算星歷。不過(guò)注意表中某些月份有兩個(gè)候選窗口4 月、7 月、8 月、12 月優(yōu)先選時(shí)長(zhǎng)更長(zhǎng)的那段。4.2 過(guò)采樣問(wèn)題和積分時(shí)間約束GOCI2 是步進(jìn)凝視成像月球在積分時(shí)間內(nèi)會(huì)相對(duì)于探測(cè)器運(yùn)動(dòng)。手冊(cè)給了一個(gè)讓人立刻緊張起來(lái)的數(shù)字單幀積分時(shí)間最長(zhǎng) 4.5 秒而月球在 1 秒內(nèi)大概移動(dòng) 10 個(gè)像元像元視場(chǎng)約 7 微弧度4.5 秒最壞情況下移動(dòng)約 50 個(gè)像元。也就是說(shuō)如果積分時(shí)間不當(dāng)月球圖像會(huì)拖成一條亮帶。所以月球定標(biāo)圖像處理不是簡(jiǎn)單把 DN 值拿出來(lái)求和。手冊(cè)給出的過(guò)采樣校正表達(dá)式為C_SAM(t) arctan(3476.4 km / D_MSat(t)) / Y_Moon(t)分子是月球?qū)嶋H角半徑分母是從圖像里量到的月球角尺寸。兩者之比就是過(guò)采樣校正因子。實(shí)際操作中先用質(zhì)心法定出月球像的位置再做亞像元插值最后統(tǒng)計(jì)月盤(pán)包含的像素?cái)?shù)這樣算立體角才靠譜。4.3 把 ROLO 系數(shù)表轉(zhuǎn)成可查詢 JSON手冊(cè)里的系數(shù)是純文本表格人工錄入很容易錯(cuò)行。我習(xí)慣先把它轉(zhuǎn)成 JSON再和定標(biāo)處理代碼對(duì)接。下面這個(gè)腳本能處理原始系數(shù)行把每行解析成以波段名為鍵的字典。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, encodingutf-8) as f: lines f.readlines() table parse_rolo_coefficients(lines) with open(rolo_coefs.json, w, encodingutf-8) as f: json.dump(table, f, indent2, ensure_asciiFalse)這個(gè)腳本有個(gè)假設(shè)每行系數(shù)嚴(yán)格按空格切分且波段名在第一位、中心波長(zhǎng)在第二位。從我碰過(guò)的轉(zhuǎn)錄文本看行列基本能對(duì)上但要注意 B13 那行只有波段名和帶寬解析時(shí)會(huì)因?yàn)樽侄尾蛔惚?continue 跳過(guò)這是正常的。JSON 輸出結(jié)構(gòu)方便下游程序按波段索引系數(shù)也方便對(duì)比不同版本手冊(cè)。4.4 從月球圖像到輻照度的三步處理拿到月球圖像后我慣用的計(jì)算流程是三步輻射定標(biāo)每個(gè)像元輻亮度 L (DN - DN_dark) / GG 是輻射響應(yīng)增益DN_dark 是暗電流扣除值來(lái)自該幀周邊的深空像元。立體角計(jì)算Ω_pixel (GSD / H)2GSD 按 250m 帶入H 取手冊(cè)中的衛(wèi)星軌道高度 35786km得到每像元立體角如果需要更精細(xì)可以逐像元按地面采樣距離空間變化修正。積分求和I_GC2(B,t) Σ L(i) · Ω(i)遍歷月盤(pán)覆蓋的所有像元得到該時(shí)刻 GOCI2 實(shí)測(cè)的月球輻照度。三步做完手里就有兩個(gè)輻照度模型期望值 I_ROLO 和實(shí)測(cè)值 I_GC2。后面所有定標(biāo)結(jié)論都圍繞這兩者的比值展開(kāi)。5. 避坑與常見(jiàn)問(wèn)題波段越界、系數(shù)錯(cuò)位和定標(biāo)因子誤用5.1 表格數(shù)字連排把“61”讀成“616”現(xiàn)象用波段表里的閾值輻亮度做飽和掩膜結(jié)果 620nm 波段大片近岸像元被誤判為飽和而 660nm 波段又漏掉真正的飽和像元。反復(fù)核對(duì)后發(fā)現(xiàn)是原始轉(zhuǎn)錄稿把 620nm 一行的“65.5 62.95 0.048 1102”連成了“65.5629.50.0481102”拆列時(shí)把 62.95 的整數(shù)部分“62”并進(jìn)了閾值輻亮度。原因這份手冊(cè)源自演示文稿導(dǎo)出 PDF表格轉(zhuǎn)文本后分隔符全部丟失數(shù)字會(huì)按視覺(jué)分組粘在一起。任何自動(dòng)解析都不能保證列對(duì)齊。解決解析前先按列單位恢復(fù)分隔。輻亮度單位是 W/(m2·μm·sr)NEdL 和 SNR 都是無(wú)量綱或比值抓住這幾個(gè)特征做正則切分切完后把每列數(shù)值范圍做合理性檢查比如最大云輻亮度不能全比最大海洋輻亮度小出現(xiàn)這種行就要手工修。5.2 ROLO 系數(shù)表缺 c1-c4滿月前后誤差放大現(xiàn)象用第 3.2 節(jié)的公式計(jì)算各月定標(biāo)期望輻照度發(fā)現(xiàn)滿月附近的殘差系統(tǒng)性偏大殘差方向還隨月面經(jīng)度變化。原因公式里有 c1 到 c4 四個(gè)幾何系數(shù)但 GOCI2 手冊(cè)里的系數(shù)表只有 a、b、d 三類(lèi)。直接調(diào)用時(shí)如果把這些缺失項(xiàng)當(dāng) 0等于忽略了觀測(cè)幾何對(duì)反照率的二階修正。近滿月時(shí)段對(duì)沖效應(yīng)最強(qiáng)d 類(lèi)系數(shù)能吸收一部分但仍不夠補(bǔ)償 c 系數(shù)的缺失。解決去 ROLO 模型發(fā)布站點(diǎn)把完整系數(shù)文件找出來(lái)按波段中心波長(zhǎng)插值出 GOCI2 對(duì)應(yīng)的 c1-c4或者至少在相位角大于 40° 的觀測(cè)窗口不接受 1% 以下的精度聲明。我個(gè)人的處理是用完整 ROLO 系數(shù)重新生成一張 GOCI2 專(zhuān)用系數(shù)表而不是零散補(bǔ)齊。5.3 745nm 和 643.5nm 兩個(gè)寬波段混為一談現(xiàn)象組里新人處理 GOCI2 月球數(shù)據(jù)時(shí)把 643.5nm/483nm 這行也當(dāng)正式波段去找定標(biāo)系數(shù)結(jié)果程序直接報(bào)錯(cuò)另一次是把 745nm 的帶寬誤寫(xiě)成 483nm導(dǎo)致光譜響應(yīng)積分窗口寬到離譜。原因手冊(cè)表格末尾的 B13 候選行沒(méi)有填任何性能參數(shù)但在原始文本里和 745nm 行靠得很近OCR 很容易把“643.5 nm 483 nm”讀成 745nm 的補(bǔ)充信息。解決使用波段表前先做一次白名單過(guò)濾正式波段中心波長(zhǎng)必須是 [380, 412, 443, 490, 510, 555, 620, 660, 680, 709, 745, 865] 之一帶寬必須是 [10, 20, 40] 之一。不在這個(gè)集合里的行一律不進(jìn)處理流程。這個(gè)白名單才是我心里這份手冊(cè)真正的驗(yàn)證規(guī)則。5.4 積分時(shí)間過(guò)長(zhǎng)導(dǎo)致月球圖像拖尾現(xiàn)象某次 4.2 秒積分時(shí)間獲取的月球圖像月盤(pán)長(zhǎng)軸比短軸明顯拉長(zhǎng)質(zhì)心法和圓形模板擬合法給出的月盤(pán)面積差 8%。原因月球相對(duì)衛(wèi)星的角速度約 72.9 微弧度/秒4.2 秒內(nèi)月盤(pán)移動(dòng)約 30 個(gè)像元圖像變成運(yùn)動(dòng)模糊后的橢圓。解決處理流程里強(qiáng)制加過(guò)采樣校正。先用 C_SAM(t) 估算拉伸比例再對(duì)月盤(pán)做橢圓擬合而不是圓擬合更根本的做法是縮短積分時(shí)間把月球成像積分限制在 3 秒以內(nèi)再用多幀疊加補(bǔ)信噪比。從那以后我每次做月球定標(biāo)都會(huì)先看星歷里的相對(duì)角速度超過(guò) 80 微弧度/秒就立刻砍積分時(shí)間。5.5 用單幀月球圖像做逐像元增益更新現(xiàn)象有人直接用一幀月球圖像的各像元響應(yīng)差異去給整面陣探測(cè)器做 PRNU 校正結(jié)果定標(biāo)后的海洋圖像出現(xiàn)明顯的棋盤(pán)格噪聲。原因月球盤(pán)面在單幀里只覆蓋探測(cè)器中心一小塊區(qū)域大部分像元沒(méi)有月球信號(hào)而且月面亮度本身有空間結(jié)構(gòu)盤(pán)面內(nèi)差異不等于探測(cè)器響應(yīng)差異。解決逐像元性能監(jiān)測(cè)只能走老化因子路線把多幀月球圖像疊起來(lái)做統(tǒng)計(jì)用 R_Moon 和 PRNU 聯(lián)合評(píng)估而不是用單幀圖像直接算增益。這也和手冊(cè)里的結(jié)論一致月球定標(biāo)適合做整體輻射鏈路的穩(wěn)定性監(jiān)測(cè)不適合逐像素絕對(duì)定標(biāo)。6. 進(jìn)階用 R_Moon 老化因子做在軌輻射穩(wěn)定性監(jiān)測(cè)月球定標(biāo)最有價(jià)值的產(chǎn)出不是某個(gè)單點(diǎn)的絕對(duì)輻照度而是長(zhǎng)時(shí)間序列上的穩(wěn)定性趨勢(shì)。手冊(cè)最后給了一個(gè)非常實(shí)用的比值定義R_Moon I_GC2 / I_ROLO這個(gè)比值被稱(chēng)為月球輻照度比值因子。理想情況下傳感器沒(méi)有衰減R_Moon 恒等于 1實(shí)際在軌運(yùn)行中光學(xué)鏡頭污染、探測(cè)器量子效率下降都會(huì)讓 I_GC2 逐漸低于 I_ROLOR_Moon 曲線就會(huì)隨時(shí)間緩慢向下漂移。反過(guò)來(lái)定義老化因子Δρ_Moon I_ROLO / I_GC2這個(gè)量大于 1 時(shí)說(shuō)明傳感器當(dāng)前響應(yīng)相比發(fā)射初期已經(jīng)損失了 (Δρ_Moon - 1) 的比例可以直接作為波段輻射增益的補(bǔ)償系數(shù)。我習(xí)慣按月產(chǎn)出一條 Δρ_Moon 變化曲線再疊加儀器的溫度遙測(cè)判斷變老是緩慢退化還是某個(gè)溫度突變引起的跳變。使用老化因子時(shí)有三個(gè)習(xí)慣要保持。第一滾動(dòng)窗口取至少 6 個(gè)月的數(shù)據(jù)點(diǎn)做線性擬合單次月球觀測(cè)的隨機(jī)噪聲在 1% 量級(jí)直接用單點(diǎn)跳變做結(jié)論容易誤判。第二同一批數(shù)據(jù)要用同一版 ROLO 系數(shù)手冊(cè)用的是 311g 版后續(xù)如果換成 314g 或更新版本整條曲線必須重新算否則新舊系數(shù)混用會(huì)出現(xiàn) 2% 左右的人為臺(tái)階。第三把 R_Moon 和 PRNU 分開(kāi)評(píng)估R_Moon 管整體響應(yīng)衰減PRNU 管像元間不均勻性變化。實(shí)際操作里我會(huì)寫(xiě)一個(gè)很小的監(jiān)測(cè)腳本把第 4 章的 JSON 系數(shù)表讀進(jìn)來(lái)配合每月定標(biāo)圖像的統(tǒng)計(jì)結(jié)果輸出一個(gè)只包含日期、波段、I_GC2、I_ROLO、R_Moon 五行字段的 CSV。這樣積累一年后畫(huà)出的曲線既能用來(lái)向數(shù)據(jù)分發(fā)方提校正確認(rèn)也能在反演算法里作為先驗(yàn)衰減補(bǔ)償。月球定標(biāo)并不是什么黑匣子從波段表到 ROLO 系數(shù)再到一幀月球圖像每一步都是可復(fù)算的。想起第一次跑通這套流程時(shí)我因?yàn)闆](méi)有校驗(yàn) c1-c4 系數(shù)滿月窗口的 R_Moon 一直偏到 0.98 以下折騰了兩天才意識(shí)到是缺失項(xiàng)而不是傳感器真的衰減。從那以后我每次拿到新版手冊(cè)都會(huì)先跑一遍列對(duì)齊和系數(shù)完整性檢查強(qiáng)制走完白名單校驗(yàn)再談定標(biāo)希望幫到你。本文還有配套的精品資源點(diǎn)擊獲取