全解析：波段表、ROLO模型與老化因子實(shí)現(xiàn))
簡(jiǎn)介面向海洋遙感、水色遙感研究人員及學(xué)生本手冊(cè)系統(tǒng)介紹GOCI2地球同步海洋水色成像儀的波段信息與性能指標(biāo)既適合新手入門(mén)也可作為課程教學(xué)或項(xiàng)目應(yīng)用的參考資料。文中對(duì)比GOCI與GOCI2列出各波段中心波長(zhǎng)、帶寬以及名義輻射亮度、最大海洋輻射亮度、閾值輻射亮度和最大云輻射亮度等關(guān)鍵參數(shù)并解釋NEdLSNR等定標(biāo)指標(biāo)同時(shí)梳理指向鏡、濾波輪、焦平面組件、遮光罩等設(shè)備結(jié)構(gòu)提及月球校準(zhǔn)計(jì)劃為理解衛(wèi)星觀測(cè)能力與數(shù)據(jù)質(zhì)量提供基礎(chǔ)。資源為1個(gè)PDF說(shuō)明手冊(cè)文件大小約2.93MB內(nèi)容緊湊、信息密度高。目前已有1158人學(xué)習(xí)下載特別適合從事海洋水色遙感數(shù)據(jù)處理、儀器定標(biāo)或GOCI系列數(shù)據(jù)應(yīng)用的讀者快速查閱與對(duì)照從而更好地開(kāi)展海洋環(huán)境監(jiān)測(cè)與遙感反演工作。1. GOCI-II 為什么要拿月球當(dāng)「定標(biāo)燈」波段表背后的在軌校準(zhǔn)邏輯一顆盯著海洋看的靜止軌道水色衛(wèi)星操作計(jì)劃里每個(gè)月卻要專門(mén)騰出十幾分鐘去拍月球——GOCI-II 的說(shuō)明手冊(cè)里甚至為「什么時(shí)候能拍月亮」列了七條優(yōu)先級(jí)判據(jù)。這不是天文學(xué)愛(ài)好而是遙感載荷最頭疼的問(wèn)題發(fā)射后儀器響應(yīng)會(huì)衰減太陽(yáng)漫射板自己也會(huì)老化而海洋水色的核心信號(hào)離水輻射只占傳感器入瞳輻射的百分之幾增益哪怕漂移 1%反演出的葉綠素濃度就跟著錯(cuò)。月球恰好是一面穩(wěn)定、無(wú)大氣、反射率長(zhǎng)期不變的「天然定標(biāo)板」。要搞懂 GOCI-II 的波段設(shè)計(jì)為什么是 13 個(gè)、帶寬為什么有 20nm 也有 483nm、SNR 為什么動(dòng)輒上千得先理解這套以月球?yàn)榛鶞?zhǔn)的在軌輻射定標(biāo)邏輯——而這份手冊(cè)正是把波段參數(shù)、ROLO 月球模型系數(shù)和定標(biāo)流程合在一份文檔里的完整底稿。2. 波段表不是參數(shù)堆砌13 個(gè)通道在補(bǔ)什么光譜信息2.1 從繼承到新增380865nm 的光譜覆蓋是怎么擴(kuò)展的對(duì)照手冊(cè)里的 GOCI 與 GOCI-II 波段表可以看到B1B4、B6、B8 等通道在兩代載荷中都有對(duì)應(yīng)關(guān)系而 510nm、620nm、680nm、709nm、745nm、865nm 這些行里 GOCI 波段欄是空的——它們是 GOCI-II 新增或重新分配的光譜通道。這份擴(kuò)展不是「多放幾個(gè)濾光片」那么簡(jiǎn)單每個(gè)新增波段都對(duì)應(yīng)一個(gè)具體的海洋光學(xué)問(wèn)題。海洋水色的難點(diǎn)在于傳感器看到的總輻射里大氣分子散射和氣溶膠散射貢獻(xiàn)了 80%90%真正攜帶水體信息的離水輻射只占很小一部分。要把這部分信號(hào)從大氣背景里「摳」出來(lái)就必須在近紅外波段同步測(cè)量氣溶膠的散射特性再外推回可見(jiàn)光波段做大氣校正。這就是為什么 745nm 和 865nm 這兩個(gè)近紅外通道是新增的核心865nm 處水體吸收極強(qiáng)、離水輻射近乎為零傳感器測(cè)到的基本全是大氣路徑輻射是氣溶膠校正的標(biāo)準(zhǔn)參考波段。709nm 則卡在浮游植物熒光的紅邊區(qū)域?qū)υ迦A暴發(fā)非常敏感。附帶提一句波段表末尾那個(gè) 643.5nm、帶寬 483nm 的通道在窄帶列表里顯得很突兀。常見(jiàn)的解讀是它承擔(dān)全色參考和云檢測(cè)功能——云在可見(jiàn)光寬波段有強(qiáng)且平坦的反射一個(gè)寬帶通道能把云掩膜做得更穩(wěn)同時(shí)為窄帶波段提供空間配準(zhǔn)的高信噪比基準(zhǔn)。2.2 SNR 與 NEdLSNR為什么海洋水色是「信噪比焦慮」的領(lǐng)域波段表里每行都有兩組性能指標(biāo)SNR標(biāo)稱輻射和 NEdLSNR標(biāo)稱輻射。海洋水色載荷的 SNR 普遍要求上千是因?yàn)殡x水輻射本身太弱而葉綠素濃度反演算法對(duì)輻射的微小變化極其敏感。GOCI-II 在 443nm 藍(lán)波段的 SNR 做到約 1050 量級(jí)——藍(lán)光在水體中穿透深、葉綠素吸收特征明顯但也容易被大氣散射污染所以既要有高 SNR又要有足夠的光譜純度。NEdLSNR 的全稱很長(zhǎng)本質(zhì)上就是「噪聲等效輻射差」在標(biāo)稱輻射水平下傳感器能分辨的最小輻射變化量?？磾?shù)值更有感覺(jué)表中大部分波段的 NEdLSNR 在 0.030.095 W/m2/μm/sr 之間和標(biāo)稱輻射亮度一比大約是千分之一的量級(jí)——這意味著儀器在標(biāo)準(zhǔn)場(chǎng)景下能分辨約 0.1% 的輻射差異。這個(gè)指標(biāo)對(duì)定標(biāo)特別關(guān)鍵如果每次月球成像的輻照度積分噪聲高于 NEdLSNR那老化因子的時(shí)間序列就會(huì)被噪聲淹沒(méi)根本看不出緩慢衰減趨勢(shì)。與 SNR 一起看這兩列構(gòu)成了一組互補(bǔ)信息SNR 是比值概念告訴你信號(hào)相對(duì)噪聲有多大NEdLSNR 是絕對(duì)分辨率概念告訴你光電探測(cè)鏈路到底能感知多小的輻射差。做定標(biāo)設(shè)計(jì)時(shí)前者決定單幀月球圖像的質(zhì)量夠不夠后者決定需要累積多少幀才能把老化因子的不確定度壓到可接受范圍。2.3 帶寬里的折中20nm、10nm、40nm 分別是怎么選的窄帶波段大多采用 20nm 帶寬這個(gè)數(shù)字本身是光譜純度與光子通量的折中。帶寬越窄光譜混疊越小對(duì)葉綠素?zé)晒獾忍卣鞣宓臏y(cè)量越干凈但到達(dá)焦平面的光子也越少SNR 會(huì)掉帶寬越寬信號(hào)強(qiáng)了卻可能把相鄰的大氣吸收特征或干擾譜段一起收進(jìn)來(lái)。20nm 是海洋水色領(lǐng)域多年實(shí)踐下來(lái)的平衡點(diǎn)。真正有意思的是那幾處刻意打破平衡的設(shè)計(jì)。709nm 通道收到 10nm是因?yàn)?720740nm 區(qū)間有水汽吸收帶如果帶寬太寬水汽吸收的時(shí)空變化會(huì)直接污染紅邊信號(hào)所以寧可犧牲光子通量也要把光譜窗口收窄。865nm 反過(guò)來(lái)放寬到 40nm因?yàn)樵撎幋髿馕杖?、氣溶膠信號(hào)本身平坦且較強(qiáng)放寬帶寬換 SNR 是劃算的買賣。從這兩個(gè)反例能看出來(lái)帶寬不是拍腦袋定的而是先確認(rèn)目標(biāo)特征的光譜位置再檢查周圍有沒(méi)有大氣吸收線最后才確定濾光片的通帶寬度。這套選擇和焦平面的機(jī)械結(jié)構(gòu)是配套的。手冊(cè)里列出了濾波輪組件——GOCI-II 是典型的濾波輪分時(shí)選波段方案不同波段不是同時(shí)成像而是靠旋轉(zhuǎn)濾波輪把目標(biāo)光譜依次放到焦平面上。這意味著月球定標(biāo)時(shí)13 個(gè)波段是逐通道采集的月球的相對(duì)位置和姿態(tài)在波段切換間隙仍在變化后文要講的過(guò)采樣校正就與此直接相關(guān)。焦平面組件和熱沉的設(shè)計(jì)同樣服務(wù)于定標(biāo)CMOS 暗電流隨溫度漂移熱沉把焦平面溫度穩(wěn)住NEdLSNR 才不會(huì)在長(zhǎng)時(shí)間曝光時(shí)被暗電流噪聲拖垮。3. 月球定標(biāo)鏈路ROLO 模型如何把「月亮亮度」變成已知數(shù)3.1 七條判據(jù)的優(yōu)先級(jí)月球不是想拍就能拍月球定標(biāo)的第一件事不是算模型而是確定「什么時(shí)候能拍」。手冊(cè)里列了七條優(yōu)先級(jí)判據(jù)完整還原了工程約束的層層嵌套。首先要滿足的是視場(chǎng)GOCI-II 的 FoR 是東西向 ±8.7°、南北向 ±11.0°而地球盤(pán)從靜止軌道看過(guò)去恰好也是 ±8.7°這意味著月球成像沒(méi)有任何額外的視場(chǎng)余量只能在月球進(jìn)入地球盤(pán)邊緣附近時(shí)捕獲。其次是操作時(shí)段約束GOCI-II 的標(biāo)稱運(yùn)行窗口是 05:1519:45KST月球定標(biāo)必須擠在這個(gè)時(shí)間段里不能為定標(biāo)單獨(dú)改變衛(wèi)星任務(wù)計(jì)劃。真正篩掉大部分日期的是月相判據(jù)月面照亮面積要大于 50%。這背后是 ROLO 模型的性質(zhì)——它是個(gè)經(jīng)驗(yàn)?zāi)Ｐ拖鄬?duì)不確定度隨相位角變化滿月前后模型外推最穩(wěn)而月牙階段的陰影效應(yīng)強(qiáng)、反射率擬合殘差大。再疊加「連續(xù)可觀測(cè)時(shí)間不少于 10 分鐘」「避開(kāi)地球雜散光 1° 區(qū)域」這些約束后實(shí)際情況是2019 年全年每個(gè)月的可觀測(cè)日只有 222 天不等推薦的單次成像窗口有的長(zhǎng)達(dá) 40 多分鐘有的只有 11 分鐘。定標(biāo)設(shè)計(jì)者必須提前數(shù)月把這份窗口表算好才能在不干擾常規(guī)海洋觀測(cè)的前提下把月球數(shù)據(jù)排進(jìn)任務(wù)時(shí)間線。3.2 ROLO 公式拆解13 個(gè)多項(xiàng)式系數(shù)在擬合什么物理過(guò)程ROLO 模型的核心是把月盤(pán)等效反射率寫(xiě)成幾個(gè)天文幾何變量的函數(shù)。手冊(cè)給出的形式是ln A_k Σ(i0..3) a_i·g^i Σ(j1..3) b_j·Φ^(2j-1) c1·θ c2·φ c3·Φθ c4·Φφ d1·e^(-g/p1) d2·e^(-g/p2) d3·(cos(g-p3))^p4逐項(xiàng)拆開(kāi)看就清楚它在擬合什么了。g 是絕對(duì)相位角即太陽(yáng)-月球-衛(wèi)星三者構(gòu)成的夾角月相變化本質(zhì)上就是 g 在變所以前四項(xiàng)多項(xiàng)式 a_i·g^i 描述了反射率隨月相的主趨勢(shì)。b_j·Φ^(2j-1) 里的 Φ 是太陽(yáng)的月面經(jīng)度描述的是同一月相下、陽(yáng)光從不同經(jīng)度方向照射月面時(shí)反射率的變化。θ 和 φ 是衛(wèi)星觀測(cè)的月面緯度、經(jīng)度對(duì)應(yīng)天平動(dòng)效應(yīng)——月球被潮汐鎖定但我們從不同角度看到的月面略有擺動(dòng)這部分描述了觀測(cè)幾何的微小擾動(dòng)。最有物理味道的是兩個(gè)指數(shù)項(xiàng)和一個(gè)余弦項(xiàng)。d1/d2 的指數(shù)項(xiàng)描述的是「對(duì)沖增強(qiáng)」相位角接近 0°滿月附近時(shí)月面顆粒的陰影幾乎消失反射率會(huì)突然升高這是個(gè)非線性很強(qiáng)的效應(yīng)純多項(xiàng)式擬合不好必須用指數(shù)項(xiàng)單獨(dú)捕捉。d3 的余弦項(xiàng)則是殘余的周期性起伏修正。整個(gè)公式取對(duì)數(shù)而不是直接擬合反射率是為了把動(dòng)態(tài)范圍壓縮讓低相位角的劇烈變化和高相位角的平緩變化在同一組系數(shù)下都能被線性描述。下面是把這個(gè)公式翻譯成可執(zhí)行邏輯的常見(jiàn)寫(xiě)法系數(shù)直接采用手冊(cè)公布的各波段 ROLO 版本數(shù)據(jù)# ROLO 模型計(jì)算某波段在給定觀測(cè)幾何下的月盤(pán)等效反射率 def rolo_reflectance(band_coeff, g, phi, theta, sel_lon, p1, p2, p3, p4): # band_coeff: 手冊(cè)系數(shù)表里該波段的 a0~a3, b1~b3, c1~c4, d1~d3 a, b, c, d band_coeff # 相位角主趨勢(shì)g 的四次多項(xiàng)式g 以弧度為單位 lnA a[0] a[1]*g a[2]*g**2 a[3]*g**3 # 太陽(yáng)月面經(jīng)度的奇次冪項(xiàng)描述陽(yáng)光照射方向的調(diào)制 for j in range(1, 4): lnA b[j-1] * phi**(2*j - 1) # 天平動(dòng)與交叉項(xiàng)觀測(cè)幾何的線性修正 lnA c[0]*theta c[1]*sel_lon c[2]*phi*theta c[3]*phi*sel_lon # 對(duì)沖增強(qiáng)與殘余周期項(xiàng)非線性最強(qiáng)、最需要單獨(dú)建模的部分 lnA d[0]*math.exp(-g/p1) d[1]*math.exp(-g/p2) d[2]*math.cos(g - p3)**p4 return math.exp(lnA)這段邏輯背后有個(gè)容易忽略的關(guān)鍵點(diǎn)系數(shù)表里的 a、b、d 等參數(shù)是逐波段擬合出來(lái)的B1380nm和 B12865nm的系數(shù)差異很大說(shuō)明月面反射率在不同波段對(duì)相位角的響應(yīng)不同——藍(lán)光受月壤多次散射影響更強(qiáng)近紅外則更接近單次反射主導(dǎo)。因此使用 ROLO 模型時(shí)不能把一組系數(shù)套到所有波段上必須按波段查表。手冊(cè)里 13 個(gè)波段各有獨(dú)立系數(shù)正是為了保留這種光譜依賴。3.3 從模型反射率到參考輻照度用儀器自己的光譜響應(yīng)去加權(quán)拿到月盤(pán)反射率 A(B,t) 之后下一步是把反射率轉(zhuǎn)換成傳感器入瞳處的參考輻照度。手冊(cè)里的積分公式是這個(gè)I_ROLO_Ref(B,t) (Ω_M / π) · ∫ E_SUN(λ) · R_λ(B) · A(B,t) dλ / ∫ R_λ(B) dλ這個(gè)公式最精髓的地方在于它沒(méi)有用單色光近似而是用 GOCI-II 各波段真實(shí)的光譜響應(yīng)函數(shù) R_λ(B) 做權(quán)重對(duì)太陽(yáng)光譜 E_SUN 和月面反射率的乘積做帶權(quán)積分。分母的 ∫R_λ(B)dλ 是歸一化項(xiàng)確保當(dāng) A 是常數(shù)時(shí)卷積結(jié)果等于 E_SUN·A 的平均值。Ω_M 是月盤(pán)張角對(duì)應(yīng)的立體角手冊(cè)給出的值是 6.4177×10?? sr——這個(gè)數(shù)由月球的角直徑?jīng)Q定是一個(gè)幾乎不變的幾何常量。積分結(jié)果 I_ROLO_Ref 是「參考距離下的月球光譜輻照度」但實(shí)際觀測(cè)時(shí)月球離衛(wèi)星的距離在變。月球軌道偏心導(dǎo)致月地距離有約 ±5% 的擺動(dòng)而輻照度按平方反比衰減對(duì)應(yīng)的修正量可達(dá) ±10% 量級(jí)。所以手冊(cè)緊接著給出了距離修正I_ROLO(B,t) I_ROLO_Ref(B,t) · C_DIS(t)其中 C_DIS(t) 必須用當(dāng)天的太陽(yáng)-月球、月球-衛(wèi)星實(shí)時(shí)距離計(jì)算。這一步看似只是乘個(gè)系數(shù)實(shí)際非常容易出錯(cuò)——后文避坑部分會(huì)專門(mén)展開(kāi)。4. 從月球圖像到老化因子過(guò)采樣、立體角與比率因子的實(shí)現(xiàn)細(xì)節(jié)4.1 過(guò)采樣校正為什么月球在 CMOS 上會(huì)被拉長(zhǎng)到 50 個(gè)像素GOCI-II 是步進(jìn)凝視成像機(jī)制衛(wèi)星不動(dòng)指向鏡快速轉(zhuǎn)向目標(biāo)區(qū)域后凝視積分。對(duì)靜止軌道衛(wèi)星來(lái)說(shuō)地球目標(biāo)基本不動(dòng)但月球是個(gè)例外——月球相對(duì)衛(wèi)星的角速度很明顯。手冊(cè)給出的數(shù)據(jù)很具體GOCI-II 的平均軌道角速度約 72.9 μrad/s單個(gè)像元的視場(chǎng)約 7 μrad月球成像的最長(zhǎng)積分時(shí)間 4.5 秒。簡(jiǎn)單乘一下就能算出1 秒積分內(nèi)月球在焦平面上移動(dòng)大約 10 個(gè)像元4.5 秒最壞情況是約 50 個(gè)像元。這意味著月球在圖像里不是靜止的圓盤(pán)而是一條拖長(zhǎng)的光痕。手冊(cè)明確說(shuō)由于采用 2D CMOS 探測(cè)器掃描速率校正不需要做但空間過(guò)采樣校正必須做。這背后的物理差別值得說(shuō)清楚CMOS 是全局曝光所有像素在同一時(shí)間窗口內(nèi)積分整幅月球圖像是被整體拉長(zhǎng)的不存在行間偽影但這個(gè)拉長(zhǎng)會(huì)直接影響月盤(pán)角大小的測(cè)量——圖像里月球的角尺寸 Y_Moon(t) 比理論值大逐像素輻照度積分時(shí)會(huì)把月球拖尾混入月盤(pán)范圍需要按比例收縮。手冊(cè)給出的修正系數(shù)是 C_SAM(t) arctan(3476.4 / D_MSat(t)) / Y_Moon(t)其中 3476.4km 是月球直徑D_MSat 是月衛(wèi)距。分子是星歷算出的理論角直徑分母是圖像實(shí)測(cè)角尺寸。下面的函數(shù)是落地這個(gè)公式的常見(jiàn)寫(xiě)法# 過(guò)采樣修正理論月球角直徑 / 圖像實(shí)測(cè)月球角直徑 def oversampling_correction(d_msat_km: float, y_moon_measured: float) - float: moon_diameter_km 3476.4 # 小角度近似下atan(直徑/距離) 即理論角直徑弧度 theoretical_angle math.atan(moon_diameter_km / d_msat_km) return theoretical_angle / y_moon_measuredC_SAM 小于 1 意味著圖像里的月球比理論大——典型的過(guò)采樣拖尾。處理時(shí)有兩種做法先對(duì)月球圖像做重采樣把拖尾壓縮回理論尺寸再做積分或者直接按 C_SAM 對(duì)積分輻照度做比例修正。手冊(cè)里特別強(qiáng)調(diào)詳細(xì)的過(guò)采樣模型還要考慮衛(wèi)星平臺(tái)指向穩(wěn)定度的貢獻(xiàn)不只是一個(gè)簡(jiǎn)單的標(biāo)量——平臺(tái)微振動(dòng)會(huì)把點(diǎn)擴(kuò)散函數(shù)再展寬一圈這部分在正式定標(biāo)處理時(shí)要用恒星像或月面邊緣的銳度來(lái)交叉驗(yàn)證。4.2 從 DN 到輻照度像素立體角與月盤(pán)積分掩膜完成過(guò)采樣修正后月球圖像里每個(gè)像素的 DN 值需要轉(zhuǎn)換成輻亮度再換算成月球總輻照度。手冊(cè)給出了最核心的積分式I_GC2(B,t) Σ(i1..N) L_MOON(B,t,i) · Ω(i)其中 Ω(i) 是第 i 個(gè)像素對(duì)應(yīng)的立體角由下式給出Ω(i) (GSD / 35786km)2GSD 是星下點(diǎn)地面采樣距離手冊(cè)說(shuō)明小于 250m35786km 是靜止軌道高度。這個(gè)做法的妙處在于它不依賴對(duì)月距離的絕對(duì)測(cè)量而是用儀器自身的幾何參數(shù)——每個(gè)像素在地球上對(duì)應(yīng)多大的地面腳印把這個(gè)腳印換算成像素立體角再乘上該像素測(cè)到的月球輻亮度逐像素累加得到月球總輻照度。量級(jí)感受一下250m 除以 35786km 大約是 7×10??平方后單個(gè)像素立體角約 5×10?11 sr月盤(pán)總立體角 6.4×10?? sr相當(dāng)于月盤(pán)覆蓋一千多像素。這里的工程細(xì)節(jié)在于月盤(pán)邊界的確定。月球圖像里需要一套掩膜算法把月盤(pán)像素和非月盤(pán)像素分開(kāi)直接用亮度閾值分割會(huì)受點(diǎn)擴(kuò)散函數(shù)翅膀影響邊界像素的能量既可能來(lái)自月球也可能來(lái)自背景更穩(wěn)的做法是用理論月盤(pán)半徑做掩膜再配合邊緣銳化。掩膜選大了把拖尾和雜散光算進(jìn)輻照度選小了月球邊緣的能量被切掉。這兩類錯(cuò)誤都會(huì)直接進(jìn)入最終的比率因子。4.3 比率因子與老化因子為什么月球定標(biāo)不碰逐像素增益把計(jì)算得到的儀器實(shí)測(cè)輻照度 I_GC2 與 ROLO 模型參考輻照度 I_ROLO 一比就得到月球比率因子R_Moon I_GC2 / I_ROLO這個(gè)比值是儀器在軌響應(yīng)相對(duì)模型基線的綜合漂移。如果儀器衰減I_GC2 變小R_Moon 小于 1。實(shí)際定標(biāo)里更常用的是它的倒數(shù)——老化因子 Δρ_Moon I_ROLO / I_GC2數(shù)值大于 1 表示需要把增益往上補(bǔ)償多少才能回到發(fā)射初期的響應(yīng)水平。手冊(cè)專門(mén)解釋了一個(gè)設(shè)計(jì)邊界為什么月球定標(biāo)不直接更新每個(gè)像素的輻射增益因?yàn)閱未卧虑虺上窭镌卤P(pán)只覆蓋探測(cè)器的一小部分視場(chǎng)無(wú)法覆蓋所有像素而 ROLO 模型給出的是月盤(pán)積分輻照度不是逐像素的月球輻亮度真值。所以在 GOCI-II 的定標(biāo)架構(gòu)里月球定標(biāo)只負(fù)責(zé)給出整體老化趨勢(shì)逐像素的非均勻性交給 PRNU像素響應(yīng)非均勻性數(shù)據(jù)來(lái)管。Δρ_Moon 描述的是整個(gè)探測(cè)器響應(yīng)的平均漂移PRNU 描述的是像素與像素之間的相對(duì)差異兩者相乘才是完整的逐像素輻射校正系數(shù)。這套分工是面陣探測(cè)器在軌定標(biāo)的典型設(shè)計(jì)——能拿到什么層次的信息就做什么層次的校正不越界。5. 月球定標(biāo)避坑指南相位角、距離與雜散光的五個(gè)真實(shí)陷阱月球定標(biāo)看起來(lái)就是「拍個(gè)月亮除一下模型」實(shí)際跑通全流程的人都知道坑幾乎全藏在邊界條件里。這里整理五條高頻翻車記錄每條都按「現(xiàn)象 → 原因 → 解決」梳理。5.1 拿殘?jiān)鲁鋈ザ?biāo)老化因子突然跳了 3%現(xiàn)象某次定標(biāo)時(shí)間窗太緊在月面照亮面積不足 50% 時(shí)硬拍了一組處理出來(lái)的 Δρ_Moon 比上一月的數(shù)據(jù)突兀跳變約 3%怎么看都像儀器出了故障。原因相位角偏離滿月區(qū)間后月面陰影效應(yīng)增強(qiáng)ROLO 模型的經(jīng)驗(yàn)擬合殘差被放大手冊(cè)給出的模型相對(duì)不確定度隨相位角的典型水平約 1%極端月相下誤差更大。解決嚴(yán)格執(zhí)行月相判據(jù)照亮面積低于 50% 的數(shù)據(jù)直接不納入老化趨勢(shì)序列如果實(shí)在缺數(shù)據(jù)至少把相位角記錄在元數(shù)據(jù)里后續(xù)用多相位角擬合曲線時(shí)給予降權(quán)。5.2 忘了乘距離修正同一個(gè)月兩次成像差出 8%現(xiàn)象間隔半個(gè)月的兩次月球成像算出的 I_GC2 差了將近 8%但儀器響應(yīng)不可能在這段時(shí)間里掉這么多。原因兩次成像時(shí)月球分別接近近地點(diǎn)和遠(yuǎn)地點(diǎn)月地距離相差約 5%輻照度按平方反比變化對(duì)應(yīng)約 10% 的擺動(dòng)處理鏈路里漏乘了 C_DIS(t) 距離修正項(xiàng)。解決把 C_DIS 的日期依賴參數(shù)太陽(yáng)-月球、月球-衛(wèi)星實(shí)時(shí)距離作為必填輔助數(shù)據(jù)與圖像文件一起歸檔處理腳本里對(duì)沒(méi)有星歷數(shù)據(jù)的圖像直接拒絕參與比率因子計(jì)算。5.3 月盤(pán)在圖像里大了一圈直接積分導(dǎo)致輻照度虛高現(xiàn)象按圖像里檢測(cè)到的月盤(pán)邊界做掩膜積分結(jié)果輻照度比 ROLO 模型參考值高 2%3%且每次成像偏高程度不穩(wěn)定。原因積分期間月球移動(dòng)和平臺(tái)微振動(dòng)造成過(guò)采樣月盤(pán)像被拉長(zhǎng)拖尾里的能量被劃進(jìn)月盤(pán)邊界點(diǎn)擴(kuò)散函數(shù)的翅膀也把月盤(pán)邊緣往外推了小半像素。解決先算 C_SAM若明顯偏離 1對(duì)圖像做重采樣收縮或用理論月盤(pán)半徑重新掩膜把超采樣修正強(qiáng)制放在輻照度積分之前而不是最后才乘因子的偷懶做法。5.4 把 R_Moon 當(dāng)逐像素增益直接乘圖像出了條紋現(xiàn)象某開(kāi)發(fā)者圖省事直接把 R_Moon 的倒數(shù)乘到每個(gè)像素的 DN 上做定標(biāo)結(jié)果海洋圖像出現(xiàn)規(guī)則條紋。原因R_Moon 是整體積分比值月盤(pán)只照亮探測(cè)器的一小部分像素它不包含像素間的響應(yīng)差異整體因子除到每一個(gè)像素上等于默認(rèn)所有像素衰減一致但 PRNU 的殘差立刻暴露成條紋。解決明確分工——R_Moon/Δρ_Moon 只進(jìn)趨勢(shì)監(jiān)測(cè)序列逐像素校正必須走 PRNU 路線兩者是乘法關(guān)系而不是等價(jià)關(guān)系。5.5 忽略地球雜散光月盤(pán)邊緣單側(cè)偏亮現(xiàn)象某次圖像的月盤(pán)邊緣有一個(gè)方向明顯偏亮輻照度積分偏高且不同月份偏亮方向不同。原因月球靠近地球盤(pán)邊緣時(shí)地球表面的云和太陽(yáng)反射光通過(guò)遮光罩衍射進(jìn)入視場(chǎng)形成約 1° 量級(jí)的雜散光區(qū)熱沉溫度漂移引起的暗電流背景也會(huì)在長(zhǎng)時(shí)間積分中疊加一個(gè)非均勻底。解決成像窗口選擇時(shí)執(zhí)行雜散光判據(jù)確保月球位置在地球盤(pán)邊緣 1° 以外圖像處理時(shí)用月盤(pán)周圍的暗區(qū)像素?cái)M合背景平面并扣除再對(duì)月盤(pán)做積分。這五條共同指向一個(gè)血淚經(jīng)驗(yàn)月球定標(biāo)的所有坑幾乎都發(fā)生在「把一個(gè)本應(yīng)做整體趨勢(shì)監(jiān)測(cè)的工具用在絕對(duì)測(cè)量或逐像素校正上」或者「省掉了某個(gè)幾何輔助項(xiàng)」的時(shí)候。把每條判據(jù)當(dāng)成物理約束而不是流程負(fù)擔(dān)才能真正把 R_Moon 的時(shí)間序列做出穩(wěn)定的趨勢(shì)價(jià)值。6. 為什么月球定標(biāo)只給「趨勢(shì)」不給「逐像素增益」一個(gè)值得守住的邊界回頭看 GOCI-II 的定標(biāo)架構(gòu)最值得琢磨的設(shè)計(jì)決策不是選了月球做定標(biāo)源而是明確劃定了月球定標(biāo)的能力邊界它只輸出整體老化因子不碰逐像素增益。這個(gè)邊界背后有三層理由。第一ROLO 模型本身的絕對(duì)精度有限——相對(duì)不確定度隨相位角在 1% 上下浮動(dòng)這個(gè)水平用來(lái)監(jiān)測(cè)一個(gè)緩慢的長(zhǎng)期衰減趨勢(shì)綽綽有余但用來(lái)標(biāo)定單個(gè)像素的絕對(duì)增益就不夠看了。第二月球在探測(cè)器上只覆蓋一小塊區(qū)域單次圖像連全幅探測(cè)器都鋪不滿想從中反演整幅焦平面的 PRNU 圖數(shù)據(jù)上就不成立。第三海洋水色產(chǎn)品的核心需求其實(shí)是穩(wěn)定性反演算法對(duì)儀器響應(yīng)的相對(duì)變化遠(yuǎn)比絕對(duì)響應(yīng)敏感只要保證逐月的 R_Moon 序列平滑、可預(yù)測(cè)產(chǎn)品的時(shí)間序列就不會(huì)出現(xiàn)虛假的年際波動(dòng)。理解了這層取舍再看整份手冊(cè)的落點(diǎn)就清楚了波段參數(shù)表定義了儀器能感知什么ROLO 系數(shù)定義了月球這個(gè)參考源怎么算過(guò)采樣和立體角積分定義了圖像怎么變成輻照度而比率因子和老化的分工定義了定標(biāo)的最終產(chǎn)出。四個(gè)環(huán)節(jié)每一處都在回答同一個(gè)問(wèn)題——如何讓一個(gè)不穩(wěn)定、會(huì)老化、有噪聲的光電系統(tǒng)持續(xù)輸出可比的科學(xué)數(shù)據(jù)。這套方案并不完美它放棄了逐像素絕對(duì)定標(biāo)的野心換來(lái)了一個(gè)簡(jiǎn)單、穩(wěn)健、可長(zhǎng)期運(yùn)行的穩(wěn)定性監(jiān)測(cè)鏈路。當(dāng)年某導(dǎo)師帶我做定標(biāo)處理時(shí)反復(fù)強(qiáng)調(diào)一句話定標(biāo)不是校正是先把誤差從哪兒來(lái)這件事搞明白。從那以后我每次拿到一組新的月球定標(biāo)數(shù)據(jù)都強(qiáng)制自己先走一遍「相位角檢查 → 距離修正 → 過(guò)采樣評(píng)估 → 雜散光背景減除」四步流程再談比率因子。這個(gè)習(xí)慣替我擋掉了不少看起來(lái)像儀器故障、其實(shí)是處理鏈路漏項(xiàng)的問(wèn)題。這套思路不限于 GOCI-II——任何用自然天體做在軌定標(biāo)的遙感載荷都逃不開(kāi)同樣的物理約束和同樣的邊界劃分。希望這些機(jī)制層面的拆解能幫你在看類似技術(shù)文檔時(shí)少踩幾個(gè)我們當(dāng)年踩過(guò)的坑。本文還有配套的精品資源點(diǎn)擊獲取 GOCI-II月球定標(biāo)全解析：波段表、ROLO模型與老化因子實(shí)現(xiàn)
簡(jiǎn)介面向海洋遙感、水色遙感研究人員及學(xué)生本手冊(cè)系統(tǒng)介紹GOCI2地球同步海洋水色成像儀的波段信息與性能指標(biāo)既適合新手入門(mén)也可作為課程教學(xué)或項(xiàng)目應(yīng)用的參考資料。文中對(duì)比GOCI與GOCI2列出各波段中心波長(zhǎng)、帶寬以及名義輻射亮度、最大海洋輻射亮度、閾值輻射亮度和最大云輻射亮度等關(guān)鍵參數(shù)并解釋NEdLSNR等定標(biāo)指標(biāo)同時(shí)梳理指向鏡、濾波輪、焦平面組件、遮光罩等設(shè)備結(jié)構(gòu)提及月球校準(zhǔn)計(jì)劃為理解衛(wèi)星觀測(cè)能力與數(shù)據(jù)質(zhì)量提供基礎(chǔ)。資源為1個(gè)PDF說(shuō)明手冊(cè)文件大小約2.93MB內(nèi)容緊湊、信息密度高。目前已有1158人學(xué)習(xí)下載特別適合從事海洋水色遙感數(shù)據(jù)處理、儀器定標(biāo)或GOCI系列數(shù)據(jù)應(yīng)用的讀者快速查閱與對(duì)照從而更好地開(kāi)展海洋環(huán)境監(jiān)測(cè)與遙感反演工作。1. GOCI-II 為什么要拿月球當(dāng)「定標(biāo)燈」波段表背后的在軌校準(zhǔn)邏輯一顆盯著海洋看的靜止軌道水色衛(wèi)星操作計(jì)劃里每個(gè)月卻要專門(mén)騰出十幾分鐘去拍月球——GOCI-II 的說(shuō)明手冊(cè)里甚至為「什么時(shí)候能拍月亮」列了七條優(yōu)先級(jí)判據(jù)。這不是天文學(xué)愛(ài)好而是遙感載荷最頭疼的問(wèn)題發(fā)射后儀器響應(yīng)會(huì)衰減太陽(yáng)漫射板自己也會(huì)老化而海洋水色的核心信號(hào)離水輻射只占傳感器入瞳輻射的百分之幾增益哪怕漂移 1%反演出的葉綠素濃度就跟著錯(cuò)。月球恰好是一面穩(wěn)定、無(wú)大氣、反射率長(zhǎng)期不變的「天然定標(biāo)板」。要搞懂 GOCI-II 的波段設(shè)計(jì)為什么是 13 個(gè)、帶寬為什么有 20nm 也有 483nm、SNR 為什么動(dòng)輒上千得先理解這套以月球?yàn)榛鶞?zhǔn)的在軌輻射定標(biāo)邏輯——而這份手冊(cè)正是把波段參數(shù)、ROLO 月球模型系數(shù)和定標(biāo)流程合在一份文檔里的完整底稿。2. 波段表不是參數(shù)堆砌13 個(gè)通道在補(bǔ)什么光譜信息2.1 從繼承到新增380865nm 的光譜覆蓋是怎么擴(kuò)展的對(duì)照手冊(cè)里的 GOCI 與 GOCI-II 波段表可以看到B1B4、B6、B8 等通道在兩代載荷中都有對(duì)應(yīng)關(guān)系而 510nm、620nm、680nm、709nm、745nm、865nm 這些行里 GOCI 波段欄是空的——它們是 GOCI-II 新增或重新分配的光譜通道。這份擴(kuò)展不是「多放幾個(gè)濾光片」那么簡(jiǎn)單每個(gè)新增波段都對(duì)應(yīng)一個(gè)具體的海洋光學(xué)問(wèn)題。海洋水色的難點(diǎn)在于傳感器看到的總輻射里大氣分子散射和氣溶膠散射貢獻(xiàn)了 80%90%真正攜帶水體信息的離水輻射只占很小一部分。要把這部分信號(hào)從大氣背景里「摳」出來(lái)就必須在近紅外波段同步測(cè)量氣溶膠的散射特性再外推回可見(jiàn)光波段做大氣校正。這就是為什么 745nm 和 865nm 這兩個(gè)近紅外通道是新增的核心865nm 處水體吸收極強(qiáng)、離水輻射近乎為零傳感器測(cè)到的基本全是大氣路徑輻射是氣溶膠校正的標(biāo)準(zhǔn)參考波段。709nm 則卡在浮游植物熒光的紅邊區(qū)域?qū)υ迦A暴發(fā)非常敏感。附帶提一句波段表末尾那個(gè) 643.5nm、帶寬 483nm 的通道在窄帶列表里顯得很突兀。常見(jiàn)的解讀是它承擔(dān)全色參考和云檢測(cè)功能——云在可見(jiàn)光寬波段有強(qiáng)且平坦的反射一個(gè)寬帶通道能把云掩膜做得更穩(wěn)同時(shí)為窄帶波段提供空間配準(zhǔn)的高信噪比基準(zhǔn)。2.2 SNR 與 NEdLSNR為什么海洋水色是「信噪比焦慮」的領(lǐng)域波段表里每行都有兩組性能指標(biāo)SNR標(biāo)稱輻射和 NEdLSNR標(biāo)稱輻射。海洋水色載荷的 SNR 普遍要求上千是因?yàn)殡x水輻射本身太弱而葉綠素濃度反演算法對(duì)輻射的微小變化極其敏感。GOCI-II 在 443nm 藍(lán)波段的 SNR 做到約 1050 量級(jí)——藍(lán)光在水體中穿透深、葉綠素吸收特征明顯但也容易被大氣散射污染所以既要有高 SNR又要有足夠的光譜純度。NEdLSNR 的全稱很長(zhǎng)本質(zhì)上就是「噪聲等效輻射差」在標(biāo)稱輻射水平下傳感器能分辨的最小輻射變化量?？磾?shù)值更有感覺(jué)表中大部分波段的 NEdLSNR 在 0.030.095 W/m2/μm/sr 之間和標(biāo)稱輻射亮度一比大約是千分之一的量級(jí)——這意味著儀器在標(biāo)準(zhǔn)場(chǎng)景下能分辨約 0.1% 的輻射差異。這個(gè)指標(biāo)對(duì)定標(biāo)特別關(guān)鍵如果每次月球成像的輻照度積分噪聲高于 NEdLSNR那老化因子的時(shí)間序列就會(huì)被噪聲淹沒(méi)根本看不出緩慢衰減趨勢(shì)。與 SNR 一起看這兩列構(gòu)成了一組互補(bǔ)信息SNR 是比值概念告訴你信號(hào)相對(duì)噪聲有多大NEdLSNR 是絕對(duì)分辨率概念告訴你光電探測(cè)鏈路到底能感知多小的輻射差。做定標(biāo)設(shè)計(jì)時(shí)前者決定單幀月球圖像的質(zhì)量夠不夠后者決定需要累積多少幀才能把老化因子的不確定度壓到可接受范圍。2.3 帶寬里的折中20nm、10nm、40nm 分別是怎么選的窄帶波段大多采用 20nm 帶寬這個(gè)數(shù)字本身是光譜純度與光子通量的折中。帶寬越窄光譜混疊越小對(duì)葉綠素?zé)晒獾忍卣鞣宓臏y(cè)量越干凈但到達(dá)焦平面的光子也越少SNR 會(huì)掉帶寬越寬信號(hào)強(qiáng)了卻可能把相鄰的大氣吸收特征或干擾譜段一起收進(jìn)來(lái)。20nm 是海洋水色領(lǐng)域多年實(shí)踐下來(lái)的平衡點(diǎn)。真正有意思的是那幾處刻意打破平衡的設(shè)計(jì)。709nm 通道收到 10nm是因?yàn)?720740nm 區(qū)間有水汽吸收帶如果帶寬太寬水汽吸收的時(shí)空變化會(huì)直接污染紅邊信號(hào)所以寧可犧牲光子通量也要把光譜窗口收窄。865nm 反過(guò)來(lái)放寬到 40nm因?yàn)樵撎幋髿馕杖?、氣溶膠信號(hào)本身平坦且較強(qiáng)放寬帶寬換 SNR 是劃算的買賣。從這兩個(gè)反例能看出來(lái)帶寬不是拍腦袋定的而是先確認(rèn)目標(biāo)特征的光譜位置再檢查周圍有沒(méi)有大氣吸收線最后才確定濾光片的通帶寬度。這套選擇和焦平面的機(jī)械結(jié)構(gòu)是配套的。手冊(cè)里列出了濾波輪組件——GOCI-II 是典型的濾波輪分時(shí)選波段方案不同波段不是同時(shí)成像而是靠旋轉(zhuǎn)濾波輪把目標(biāo)光譜依次放到焦平面上。這意味著月球定標(biāo)時(shí)13 個(gè)波段是逐通道采集的月球的相對(duì)位置和姿態(tài)在波段切換間隙仍在變化后文要講的過(guò)采樣校正就與此直接相關(guān)。焦平面組件和熱沉的設(shè)計(jì)同樣服務(wù)于定標(biāo)CMOS 暗電流隨溫度漂移熱沉把焦平面溫度穩(wěn)住NEdLSNR 才不會(huì)在長(zhǎng)時(shí)間曝光時(shí)被暗電流噪聲拖垮。3. 月球定標(biāo)鏈路ROLO 模型如何把「月亮亮度」變成已知數(shù)3.1 七條判據(jù)的優(yōu)先級(jí)月球不是想拍就能拍月球定標(biāo)的第一件事不是算模型而是確定「什么時(shí)候能拍」。手冊(cè)里列了七條優(yōu)先級(jí)判據(jù)完整還原了工程約束的層層嵌套。首先要滿足的是視場(chǎng)GOCI-II 的 FoR 是東西向 ±8.7°、南北向 ±11.0°而地球盤(pán)從靜止軌道看過(guò)去恰好也是 ±8.7°這意味著月球成像沒(méi)有任何額外的視場(chǎng)余量只能在月球進(jìn)入地球盤(pán)邊緣附近時(shí)捕獲。其次是操作時(shí)段約束GOCI-II 的標(biāo)稱運(yùn)行窗口是 05:1519:45KST月球定標(biāo)必須擠在這個(gè)時(shí)間段里不能為定標(biāo)單獨(dú)改變衛(wèi)星任務(wù)計(jì)劃。真正篩掉大部分日期的是月相判據(jù)月面照亮面積要大于 50%。這背后是 ROLO 模型的性質(zhì)——它是個(gè)經(jīng)驗(yàn)?zāi)Ｐ拖鄬?duì)不確定度隨相位角變化滿月前后模型外推最穩(wěn)而月牙階段的陰影效應(yīng)強(qiáng)、反射率擬合殘差大。再疊加「連續(xù)可觀測(cè)時(shí)間不少于 10 分鐘」「避開(kāi)地球雜散光 1° 區(qū)域」這些約束后實(shí)際情況是2019 年全年每個(gè)月的可觀測(cè)日只有 222 天不等推薦的單次成像窗口有的長(zhǎng)達(dá) 40 多分鐘有的只有 11 分鐘。定標(biāo)設(shè)計(jì)者必須提前數(shù)月把這份窗口表算好才能在不干擾常規(guī)海洋觀測(cè)的前提下把月球數(shù)據(jù)排進(jìn)任務(wù)時(shí)間線。3.2 ROLO 公式拆解13 個(gè)多項(xiàng)式系數(shù)在擬合什么物理過(guò)程ROLO 模型的核心是把月盤(pán)等效反射率寫(xiě)成幾個(gè)天文幾何變量的函數(shù)。手冊(cè)給出的形式是ln A_k Σ(i0..3) a_i·g^i Σ(j1..3) b_j·Φ^(2j-1) c1·θ c2·φ c3·Φθ c4·Φφ d1·e^(-g/p1) d2·e^(-g/p2) d3·(cos(g-p3))^p4逐項(xiàng)拆開(kāi)看就清楚它在擬合什么了。g 是絕對(duì)相位角即太陽(yáng)-月球-衛(wèi)星三者構(gòu)成的夾角月相變化本質(zhì)上就是 g 在變所以前四項(xiàng)多項(xiàng)式 a_i·g^i 描述了反射率隨月相的主趨勢(shì)。b_j·Φ^(2j-1) 里的 Φ 是太陽(yáng)的月面經(jīng)度描述的是同一月相下、陽(yáng)光從不同經(jīng)度方向照射月面時(shí)反射率的變化。θ 和 φ 是衛(wèi)星觀測(cè)的月面緯度、經(jīng)度對(duì)應(yīng)天平動(dòng)效應(yīng)——月球被潮汐鎖定但我們從不同角度看到的月面略有擺動(dòng)這部分描述了觀測(cè)幾何的微小擾動(dòng)。最有物理味道的是兩個(gè)指數(shù)項(xiàng)和一個(gè)余弦項(xiàng)。d1/d2 的指數(shù)項(xiàng)描述的是「對(duì)沖增強(qiáng)」相位角接近 0°滿月附近時(shí)月面顆粒的陰影幾乎消失反射率會(huì)突然升高這是個(gè)非線性很強(qiáng)的效應(yīng)純多項(xiàng)式擬合不好必須用指數(shù)項(xiàng)單獨(dú)捕捉。d3 的余弦項(xiàng)則是殘余的周期性起伏修正。整個(gè)公式取對(duì)數(shù)而不是直接擬合反射率是為了把動(dòng)態(tài)范圍壓縮讓低相位角的劇烈變化和高相位角的平緩變化在同一組系數(shù)下都能被線性描述。下面是把這個(gè)公式翻譯成可執(zhí)行邏輯的常見(jiàn)寫(xiě)法系數(shù)直接采用手冊(cè)公布的各波段 ROLO 版本數(shù)據(jù)# ROLO 模型計(jì)算某波段在給定觀測(cè)幾何下的月盤(pán)等效反射率 def rolo_reflectance(band_coeff, g, phi, theta, sel_lon, p1, p2, p3, p4): # band_coeff: 手冊(cè)系數(shù)表里該波段的 a0~a3, b1~b3, c1~c4, d1~d3 a, b, c, d band_coeff # 相位角主趨勢(shì)g 的四次多項(xiàng)式g 以弧度為單位 lnA a[0] a[1]*g a[2]*g**2 a[3]*g**3 # 太陽(yáng)月面經(jīng)度的奇次冪項(xiàng)描述陽(yáng)光照射方向的調(diào)制 for j in range(1, 4): lnA b[j-1] * phi**(2*j - 1) # 天平動(dòng)與交叉項(xiàng)觀測(cè)幾何的線性修正 lnA c[0]*theta c[1]*sel_lon c[2]*phi*theta c[3]*phi*sel_lon # 對(duì)沖增強(qiáng)與殘余周期項(xiàng)非線性最強(qiáng)、最需要單獨(dú)建模的部分 lnA d[0]*math.exp(-g/p1) d[1]*math.exp(-g/p2) d[2]*math.cos(g - p3)**p4 return math.exp(lnA)這段邏輯背后有個(gè)容易忽略的關(guān)鍵點(diǎn)系數(shù)表里的 a、b、d 等參數(shù)是逐波段擬合出來(lái)的B1380nm和 B12865nm的系數(shù)差異很大說(shuō)明月面反射率在不同波段對(duì)相位角的響應(yīng)不同——藍(lán)光受月壤多次散射影響更強(qiáng)近紅外則更接近單次反射主導(dǎo)。因此使用 ROLO 模型時(shí)不能把一組系數(shù)套到所有波段上必須按波段查表。手冊(cè)里 13 個(gè)波段各有獨(dú)立系數(shù)正是為了保留這種光譜依賴。3.3 從模型反射率到參考輻照度用儀器自己的光譜響應(yīng)去加權(quán)拿到月盤(pán)反射率 A(B,t) 之后下一步是把反射率轉(zhuǎn)換成傳感器入瞳處的參考輻照度。手冊(cè)里的積分公式是這個(gè)I_ROLO_Ref(B,t) (Ω_M / π) · ∫ E_SUN(λ) · R_λ(B) · A(B,t) dλ / ∫ R_λ(B) dλ這個(gè)公式最精髓的地方在于它沒(méi)有用單色光近似而是用 GOCI-II 各波段真實(shí)的光譜響應(yīng)函數(shù) R_λ(B) 做權(quán)重對(duì)太陽(yáng)光譜 E_SUN 和月面反射率的乘積做帶權(quán)積分。分母的 ∫R_λ(B)dλ 是歸一化項(xiàng)確保當(dāng) A 是常數(shù)時(shí)卷積結(jié)果等于 E_SUN·A 的平均值。Ω_M 是月盤(pán)張角對(duì)應(yīng)的立體角手冊(cè)給出的值是 6.4177×10?? sr——這個(gè)數(shù)由月球的角直徑?jīng)Q定是一個(gè)幾乎不變的幾何常量。積分結(jié)果 I_ROLO_Ref 是「參考距離下的月球光譜輻照度」但實(shí)際觀測(cè)時(shí)月球離衛(wèi)星的距離在變。月球軌道偏心導(dǎo)致月地距離有約 ±5% 的擺動(dòng)而輻照度按平方反比衰減對(duì)應(yīng)的修正量可達(dá) ±10% 量級(jí)。所以手冊(cè)緊接著給出了距離修正I_ROLO(B,t) I_ROLO_Ref(B,t) · C_DIS(t)其中 C_DIS(t) 必須用當(dāng)天的太陽(yáng)-月球、月球-衛(wèi)星實(shí)時(shí)距離計(jì)算。這一步看似只是乘個(gè)系數(shù)實(shí)際非常容易出錯(cuò)——后文避坑部分會(huì)專門(mén)展開(kāi)。4. 從月球圖像到老化因子過(guò)采樣、立體角與比率因子的實(shí)現(xiàn)細(xì)節(jié)4.1 過(guò)采樣校正為什么月球在 CMOS 上會(huì)被拉長(zhǎng)到 50 個(gè)像素GOCI-II 是步進(jìn)凝視成像機(jī)制衛(wèi)星不動(dòng)指向鏡快速轉(zhuǎn)向目標(biāo)區(qū)域后凝視積分。對(duì)靜止軌道衛(wèi)星來(lái)說(shuō)地球目標(biāo)基本不動(dòng)但月球是個(gè)例外——月球相對(duì)衛(wèi)星的角速度很明顯。手冊(cè)給出的數(shù)據(jù)很具體GOCI-II 的平均軌道角速度約 72.9 μrad/s單個(gè)像元的視場(chǎng)約 7 μrad月球成像的最長(zhǎng)積分時(shí)間 4.5 秒。簡(jiǎn)單乘一下就能算出1 秒積分內(nèi)月球在焦平面上移動(dòng)大約 10 個(gè)像元4.5 秒最壞情況是約 50 個(gè)像元。這意味著月球在圖像里不是靜止的圓盤(pán)而是一條拖長(zhǎng)的光痕。手冊(cè)明確說(shuō)由于采用 2D CMOS 探測(cè)器掃描速率校正不需要做但空間過(guò)采樣校正必須做。這背后的物理差別值得說(shuō)清楚CMOS 是全局曝光所有像素在同一時(shí)間窗口內(nèi)積分整幅月球圖像是被整體拉長(zhǎng)的不存在行間偽影但這個(gè)拉長(zhǎng)會(huì)直接影響月盤(pán)角大小的測(cè)量——圖像里月球的角尺寸 Y_Moon(t) 比理論值大逐像素輻照度積分時(shí)會(huì)把月球拖尾混入月盤(pán)范圍需要按比例收縮。手冊(cè)給出的修正系數(shù)是 C_SAM(t) arctan(3476.4 / D_MSat(t)) / Y_Moon(t)其中 3476.4km 是月球直徑D_MSat 是月衛(wèi)距。分子是星歷算出的理論角直徑分母是圖像實(shí)測(cè)角尺寸。下面的函數(shù)是落地這個(gè)公式的常見(jiàn)寫(xiě)法# 過(guò)采樣修正理論月球角直徑 / 圖像實(shí)測(cè)月球角直徑 def oversampling_correction(d_msat_km: float, y_moon_measured: float) - float: moon_diameter_km 3476.4 # 小角度近似下atan(直徑/距離) 即理論角直徑弧度 theoretical_angle math.atan(moon_diameter_km / d_msat_km) return theoretical_angle / y_moon_measuredC_SAM 小于 1 意味著圖像里的月球比理論大——典型的過(guò)采樣拖尾。處理時(shí)有兩種做法先對(duì)月球圖像做重采樣把拖尾壓縮回理論尺寸再做積分或者直接按 C_SAM 對(duì)積分輻照度做比例修正。手冊(cè)里特別強(qiáng)調(diào)詳細(xì)的過(guò)采樣模型還要考慮衛(wèi)星平臺(tái)指向穩(wěn)定度的貢獻(xiàn)不只是一個(gè)簡(jiǎn)單的標(biāo)量——平臺(tái)微振動(dòng)會(huì)把點(diǎn)擴(kuò)散函數(shù)再展寬一圈這部分在正式定標(biāo)處理時(shí)要用恒星像或月面邊緣的銳度來(lái)交叉驗(yàn)證。4.2 從 DN 到輻照度像素立體角與月盤(pán)積分掩膜完成過(guò)采樣修正后月球圖像里每個(gè)像素的 DN 值需要轉(zhuǎn)換成輻亮度再換算成月球總輻照度。手冊(cè)給出了最核心的積分式I_GC2(B,t) Σ(i1..N) L_MOON(B,t,i) · Ω(i)其中 Ω(i) 是第 i 個(gè)像素對(duì)應(yīng)的立體角由下式給出Ω(i) (GSD / 35786km)2GSD 是星下點(diǎn)地面采樣距離手冊(cè)說(shuō)明小于 250m35786km 是靜止軌道高度。這個(gè)做法的妙處在于它不依賴對(duì)月距離的絕對(duì)測(cè)量而是用儀器自身的幾何參數(shù)——每個(gè)像素在地球上對(duì)應(yīng)多大的地面腳印把這個(gè)腳印換算成像素立體角再乘上該像素測(cè)到的月球輻亮度逐像素累加得到月球總輻照度。量級(jí)感受一下250m 除以 35786km 大約是 7×10??平方后單個(gè)像素立體角約 5×10?11 sr月盤(pán)總立體角 6.4×10?? sr相當(dāng)于月盤(pán)覆蓋一千多像素。這里的工程細(xì)節(jié)在于月盤(pán)邊界的確定。月球圖像里需要一套掩膜算法把月盤(pán)像素和非月盤(pán)像素分開(kāi)直接用亮度閾值分割會(huì)受點(diǎn)擴(kuò)散函數(shù)翅膀影響邊界像素的能量既可能來(lái)自月球也可能來(lái)自背景更穩(wěn)的做法是用理論月盤(pán)半徑做掩膜再配合邊緣銳化。掩膜選大了把拖尾和雜散光算進(jìn)輻照度選小了月球邊緣的能量被切掉。這兩類錯(cuò)誤都會(huì)直接進(jìn)入最終的比率因子。4.3 比率因子與老化因子為什么月球定標(biāo)不碰逐像素增益把計(jì)算得到的儀器實(shí)測(cè)輻照度 I_GC2 與 ROLO 模型參考輻照度 I_ROLO 一比就得到月球比率因子R_Moon I_GC2 / I_ROLO這個(gè)比值是儀器在軌響應(yīng)相對(duì)模型基線的綜合漂移。如果儀器衰減I_GC2 變小R_Moon 小于 1。實(shí)際定標(biāo)里更常用的是它的倒數(shù)——老化因子 Δρ_Moon I_ROLO / I_GC2數(shù)值大于 1 表示需要把增益往上補(bǔ)償多少才能回到發(fā)射初期的響應(yīng)水平。手冊(cè)專門(mén)解釋了一個(gè)設(shè)計(jì)邊界為什么月球定標(biāo)不直接更新每個(gè)像素的輻射增益因?yàn)閱未卧虑虺上窭镌卤P(pán)只覆蓋探測(cè)器的一小部分視場(chǎng)無(wú)法覆蓋所有像素而 ROLO 模型給出的是月盤(pán)積分輻照度不是逐像素的月球輻亮度真值。所以在 GOCI-II 的定標(biāo)架構(gòu)里月球定標(biāo)只負(fù)責(zé)給出整體老化趨勢(shì)逐像素的非均勻性交給 PRNU像素響應(yīng)非均勻性數(shù)據(jù)來(lái)管。Δρ_Moon 描述的是整個(gè)探測(cè)器響應(yīng)的平均漂移PRNU 描述的是像素與像素之間的相對(duì)差異兩者相乘才是完整的逐像素輻射校正系數(shù)。這套分工是面陣探測(cè)器在軌定標(biāo)的典型設(shè)計(jì)——能拿到什么層次的信息就做什么層次的校正不越界。5. 月球定標(biāo)避坑指南相位角、距離與雜散光的五個(gè)真實(shí)陷阱月球定標(biāo)看起來(lái)就是「拍個(gè)月亮除一下模型」實(shí)際跑通全流程的人都知道坑幾乎全藏在邊界條件里。這里整理五條高頻翻車記錄每條都按「現(xiàn)象 → 原因 → 解決」梳理。5.1 拿殘?jiān)鲁鋈ザ?biāo)老化因子突然跳了 3%現(xiàn)象某次定標(biāo)時(shí)間窗太緊在月面照亮面積不足 50% 時(shí)硬拍了一組處理出來(lái)的 Δρ_Moon 比上一月的數(shù)據(jù)突兀跳變約 3%怎么看都像儀器出了故障。原因相位角偏離滿月區(qū)間后月面陰影效應(yīng)增強(qiáng)ROLO 模型的經(jīng)驗(yàn)擬合殘差被放大手冊(cè)給出的模型相對(duì)不確定度隨相位角的典型水平約 1%極端月相下誤差更大。解決嚴(yán)格執(zhí)行月相判據(jù)照亮面積低于 50% 的數(shù)據(jù)直接不納入老化趨勢(shì)序列如果實(shí)在缺數(shù)據(jù)至少把相位角記錄在元數(shù)據(jù)里后續(xù)用多相位角擬合曲線時(shí)給予降權(quán)。5.2 忘了乘距離修正同一個(gè)月兩次成像差出 8%現(xiàn)象間隔半個(gè)月的兩次月球成像算出的 I_GC2 差了將近 8%但儀器響應(yīng)不可能在這段時(shí)間里掉這么多。原因兩次成像時(shí)月球分別接近近地點(diǎn)和遠(yuǎn)地點(diǎn)月地距離相差約 5%輻照度按平方反比變化對(duì)應(yīng)約 10% 的擺動(dòng)處理鏈路里漏乘了 C_DIS(t) 距離修正項(xiàng)。解決把 C_DIS 的日期依賴參數(shù)太陽(yáng)-月球、月球-衛(wèi)星實(shí)時(shí)距離作為必填輔助數(shù)據(jù)與圖像文件一起歸檔處理腳本里對(duì)沒(méi)有星歷數(shù)據(jù)的圖像直接拒絕參與比率因子計(jì)算。5.3 月盤(pán)在圖像里大了一圈直接積分導(dǎo)致輻照度虛高現(xiàn)象按圖像里檢測(cè)到的月盤(pán)邊界做掩膜積分結(jié)果輻照度比 ROLO 模型參考值高 2%3%且每次成像偏高程度不穩(wěn)定。原因積分期間月球移動(dòng)和平臺(tái)微振動(dòng)造成過(guò)采樣月盤(pán)像被拉長(zhǎng)拖尾里的能量被劃進(jìn)月盤(pán)邊界點(diǎn)擴(kuò)散函數(shù)的翅膀也把月盤(pán)邊緣往外推了小半像素。解決先算 C_SAM若明顯偏離 1對(duì)圖像做重采樣收縮或用理論月盤(pán)半徑重新掩膜把超采樣修正強(qiáng)制放在輻照度積分之前而不是最后才乘因子的偷懶做法。5.4 把 R_Moon 當(dāng)逐像素增益直接乘圖像出了條紋現(xiàn)象某開(kāi)發(fā)者圖省事直接把 R_Moon 的倒數(shù)乘到每個(gè)像素的 DN 上做定標(biāo)結(jié)果海洋圖像出現(xiàn)規(guī)則條紋。原因R_Moon 是整體積分比值月盤(pán)只照亮探測(cè)器的一小部分像素它不包含像素間的響應(yīng)差異整體因子除到每一個(gè)像素上等于默認(rèn)所有像素衰減一致但 PRNU 的殘差立刻暴露成條紋。解決明確分工——R_Moon/Δρ_Moon 只進(jìn)趨勢(shì)監(jiān)測(cè)序列逐像素校正必須走 PRNU 路線兩者是乘法關(guān)系而不是等價(jià)關(guān)系。5.5 忽略地球雜散光月盤(pán)邊緣單側(cè)偏亮現(xiàn)象某次圖像的月盤(pán)邊緣有一個(gè)方向明顯偏亮輻照度積分偏高且不同月份偏亮方向不同。原因月球靠近地球盤(pán)邊緣時(shí)地球表面的云和太陽(yáng)反射光通過(guò)遮光罩衍射進(jìn)入視場(chǎng)形成約 1° 量級(jí)的雜散光區(qū)熱沉溫度漂移引起的暗電流背景也會(huì)在長(zhǎng)時(shí)間積分中疊加一個(gè)非均勻底。解決成像窗口選擇時(shí)執(zhí)行雜散光判據(jù)確保月球位置在地球盤(pán)邊緣 1° 以外圖像處理時(shí)用月盤(pán)周圍的暗區(qū)像素?cái)M合背景平面并扣除再對(duì)月盤(pán)做積分。這五條共同指向一個(gè)血淚經(jīng)驗(yàn)月球定標(biāo)的所有坑幾乎都發(fā)生在「把一個(gè)本應(yīng)做整體趨勢(shì)監(jiān)測(cè)的工具用在絕對(duì)測(cè)量或逐像素校正上」或者「省掉了某個(gè)幾何輔助項(xiàng)」的時(shí)候。把每條判據(jù)當(dāng)成物理約束而不是流程負(fù)擔(dān)才能真正把 R_Moon 的時(shí)間序列做出穩(wěn)定的趨勢(shì)價(jià)值。6. 為什么月球定標(biāo)只給「趨勢(shì)」不給「逐像素增益」一個(gè)值得守住的邊界回頭看 GOCI-II 的定標(biāo)架構(gòu)最值得琢磨的設(shè)計(jì)決策不是選了月球做定標(biāo)源而是明確劃定了月球定標(biāo)的能力邊界它只輸出整體老化因子不碰逐像素增益。這個(gè)邊界背后有三層理由。第一ROLO 模型本身的絕對(duì)精度有限——相對(duì)不確定度隨相位角在 1% 上下浮動(dòng)這個(gè)水平用來(lái)監(jiān)測(cè)一個(gè)緩慢的長(zhǎng)期衰減趨勢(shì)綽綽有余但用來(lái)標(biāo)定單個(gè)像素的絕對(duì)增益就不夠看了。第二月球在探測(cè)器上只覆蓋一小塊區(qū)域單次圖像連全幅探測(cè)器都鋪不滿想從中反演整幅焦平面的 PRNU 圖數(shù)據(jù)上就不成立。第三海洋水色產(chǎn)品的核心需求其實(shí)是穩(wěn)定性反演算法對(duì)儀器響應(yīng)的相對(duì)變化遠(yuǎn)比絕對(duì)響應(yīng)敏感只要保證逐月的 R_Moon 序列平滑、可預(yù)測(cè)產(chǎn)品的時(shí)間序列就不會(huì)出現(xiàn)虛假的年際波動(dòng)。理解了這層取舍再看整份手冊(cè)的落點(diǎn)就清楚了波段參數(shù)表定義了儀器能感知什么ROLO 系數(shù)定義了月球這個(gè)參考源怎么算過(guò)采樣和立體角積分定義了圖像怎么變成輻照度而比率因子和老化的分工定義了定標(biāo)的最終產(chǎn)出。四個(gè)環(huán)節(jié)每一處都在回答同一個(gè)問(wèn)題——如何讓一個(gè)不穩(wěn)定、會(huì)老化、有噪聲的光電系統(tǒng)持續(xù)輸出可比的科學(xué)數(shù)據(jù)。這套方案并不完美它放棄了逐像素絕對(duì)定標(biāo)的野心換來(lái)了一個(gè)簡(jiǎn)單、穩(wěn)健、可長(zhǎng)期運(yùn)行的穩(wěn)定性監(jiān)測(cè)鏈路。當(dāng)年某導(dǎo)師帶我做定標(biāo)處理時(shí)反復(fù)強(qiáng)調(diào)一句話定標(biāo)不是校正是先把誤差從哪兒來(lái)這件事搞明白。從那以后我每次拿到一組新的月球定標(biāo)數(shù)據(jù)都強(qiáng)制自己先走一遍「相位角檢查 → 距離修正 → 過(guò)采樣評(píng)估 → 雜散光背景減除」四步流程再談比率因子。這個(gè)習(xí)慣替我擋掉了不少看起來(lái)像儀器故障、其實(shí)是處理鏈路漏項(xiàng)的問(wèn)題。這套思路不限于 GOCI-II——任何用自然天體做在軌定標(biāo)的遙感載荷都逃不開(kāi)同樣的物理約束和同樣的邊界劃分。希望這些機(jī)制層面的拆解能幫你在看類似技術(shù)文檔時(shí)少踩幾個(gè)我們當(dāng)年踩過(guò)的坑。本文還有配套的精品資源點(diǎn)擊獲取