靜軌紫外可見高光譜探測儀星上外定標方案

崔程光，李 凌，李云飛，郭永祥，李永強

（北京空間機電研究所，北京 100094）

1 引 言

星上定標直接反映了遙感器入軌以后儀器性能的實際情況，同時為了保證遙感器在軌工作的長期有效性，以及遙感信息的高精度獲取，遙感器需要配備星上定標系統，確定星上定標方案。星上定標方案是遙感器長期性能變化監測的重要手段，也是定量獲取遙感器在軌性能信息，校正在軌遙感數據的基礎［1］。

目前，國內外通常采用星上內定標系統和星上外定標系統。由于星上內定標系統只能完成部分光學系統監測，當未被標定的系統發生衰減或退化時會被認為是目標輻射量的變化。此外星上標準燈和太陽的光譜分布差異較大，需要將標準燈的定標結果進行光譜匹配處理，這會加大最終定標結果的合成不確定度［2-3］。對于長期在軌工作的遙感器，定標燈自身的衰減會引入其定標結果中，因此需要采用星上外定標系統。

太陽定標是目前比較成熟的星上定標方案［4］，但是由于幾何位置的約束，太陽定標不便于頻繁進行，只能在軌道的某幾個固定位置進行。月球定標同樣面臨這類問題，目前國際上通常采用交叉定標，但數據量有限［5］。美國通過Sea-WiFs 以 及 MODIS 中 的 Terra 和 Aqua 進 行 對 月交叉定標［6］。用“自動月球觀測”（Robotic Lunar Observatory，ROLO）模型標定 Terra MODIS 和Aqua MODIS 的響應度一致性在1%～3%，Sea-WiFs 與兩臺MODIS 的響應度一致性為3%～8%［7］。美國早在1991 年進行了第一代太陽恒星輻照度比照實驗（SOLSTICE I），以期以恒星監測太陽變化。他們采用31 顆恒星和太陽進行比對，其絕對精度達10%、相對精度達2%。美國又于2003 年發射了SOLSTICE II，其絕對精度為5%，相對精度為0.5%［8-9］。德國等也開展了太陽恒星比對研究工作［10］。星上外基準源逐漸豐富，但是由于遙感器多樣化和衛星平臺約束，單一空間外基準難以滿足遙感器長期在軌定標需求，因此需要對星上外定標方案進行可行性分析，制定合理的星上外定標方案。

恒星理論表明，B 型和 A 型恒星（B，A 源于以溫度分界的哈佛分類表示方法）在100～400 nm 波段內的穩定性優于1%。為了進一步提高在軌監測精度，本文提出一種采用恒星進行太陽輻射校正的方法。結合儀器特點和軌道幾何約束研究了紫外可見高光譜探測儀在軌高精度輻射性能監測方案，同時結合月球定標實現星上外定標核心部件——漫反射板的在軌監測。然后，利用衛星工具包（STK）驗證了在靜止軌道下太陽定標、月球定標和恒星定標方案的可行性，最終確定采用太陽、月球、恒星相結合的方式實現儀器星上定標。該定標方案可為靜軌紫外可見高光譜探測儀的在軌監測提供參考。

2 理論分析

紫外可見高光譜探測儀獲取的目標光譜輻射可表示為：

式中：I(λ)為入射光譜通量，S(λ)為響應信號，Sstray(λ)為雜散光響應輸出，Sdark(λ)為儀器暗信號。由于采用退偏器約束儀器的偏振響應，可忽略儀器的偏振影響［11］。

探測儀光譜響應度的在軌長期可靠性，或實時變化的在軌監測是獲取高精度目標光譜輻射信息的關鍵，因此需要對其在軌監測進行研究［12］。內定標方案不僅難以對全光路進行系統監測，而且成本較高，因此需要研究靜止軌道儀器的星上外定標方案，進而實現對探測儀關鍵輻射性能的監測。

2.1 儀器在軌運行t 時間后光譜輻亮度響應度

當進行太陽定標時，采用“太陽+俯仰漫反射板”的方式完成，此時儀器的光譜響應度為：

在同一星上外定標周期內完成恒星定標，假設探測儀在短時間內沒有衰減等影響，此時可得恒星的目標光譜輻照度為：

式（3）與式（2）聯立，可得太陽、恒星輻照度的關系，即：

通過分析可知，靜止軌道遙感器的特殊幾何位置不能實現每周一次的太陽定標，但對于部分恒星，遙感器每天均可見。因此，在遙感器對太陽和恒星均可見的時間內，太陽與恒星比對，并在太陽不可見周期內進行恒星定標。應用恒星對遙感器的光譜輻亮度每周進行一次定標。

由于恒星能量較弱，需要計算恒星定標的信噪比。恒星的具體位置、亮度和光譜特性及其表面溫度可以通過查閱星表得到，將恒星近似為黑體，根據黑體輻射定律，由Plank 公式表達恒星的光譜輻射亮度分布［13］，即：

式中：Bλ為恒星的光譜輻亮度分布，h為普朗克常數，c為光速，kB為玻爾茲曼常數，T為黑體溫度。

靜軌紫外可見高光譜探測儀單個像元接收到的光譜能量為：

式中：AD為探測器像元尺寸，τ(λ)為儀器的光譜透過率，D為探測儀入瞳大小，ES(λ)為探測儀入射窗口接收到的光譜輻照度。

探測器接收到的光子數為：

式中：η為探測器的量子效率，τ為積分時間。

探測儀的信噪比可以表示為［14］：

式中：Idark為暗信號，Dread為讀出噪聲。

2.2 儀器在軌運行t 時間后漫反射板BRDF

漫反射板衰減采用月球定標。假定漫反射板在衛星發射過程中其BRDF 不發生變化，即儀器在軌工作時第一次完成定標時漫反射板的BRDF 與實驗室標定值相同。完成太陽定標時有［15］：

完成對月定標時有：

儀器在軌工作一段時間t以后，漫反射板的衰減因子為：

式中，Lmoon(λ，t)Lmoon(λ，0)可由美國 ROLO 的月球模型確定。此時只需計算月球光譜輻亮度的相對變化。

3 可行性分析

紫外可見高光譜探測儀視場為±12°×6.4°，可以實現全球東西覆蓋和南北4 000 km覆蓋，為實現對我國的完全覆蓋，紫外高光譜探測儀的視場中心向北偏轉5°，遙感器的軌道參數如表1 所示。

考慮靜止軌道特點，需要分析太陽、恒星、月球各種定標方案的可行性，驗證各定標方案，為紫外可見可見高光譜探測儀的在軌應用提供支撐。STK 分析流程如圖1 所示。

表1 紫外可見高光譜探測儀軌道參數Tab.1 Orbit parameters of GUVI high-spectral instrument

圖1 衛星工具包（STK）分析流程Fig.1 Simulation analysis process of Satellite Tool Kit（STK）

3.1 太陽定標

在STK 中建立仿真場景，如圖2 所示，以遙感器為坐標系主體，x為遙感器飛行方向，z為遙感器中心視場方向，y由右手坐標系確定。

圖2 靜軌紫外可見高光譜探測儀太陽定標的STK 仿真場景Fig.2 Sun calibration simulation of GUVI in STK

仿真場景的數學模型簡圖如圖3 所示，定義衛星本體豎直向下為Z軸，它與z軸的夾角為5°。定義α表示太陽和平面YZ的夾角，β表示太陽和平面XZ的夾角，θ為太陽和Z軸的夾角。當太陽處于探測儀視場中，同時太陽與地球位于同一象限內且不被地球所遮擋時可以進行太陽定標，即α∈[-12°，12°]，β∈[-3.2°，3.2°]，θ∈ (arcsin即θ∈[8.69°，90°)。 那 么 在太陽可見時間內存在可用角度，如圖4 所示。

圖3 太陽定標仿真數學簡圖Fig.3 Mathematical model in sun calibration simulation

圖4 太陽可見時刻內α，β，θ 角度變化Fig.4 Variations of α，β，θ for sun is visible

以 2016 至 2017 年為例，由圖 5 可知在一年內太陽僅在三月下旬、四月上旬，以及九月上旬和中旬可見，不能滿足長期太陽定標。

圖5 太陽的年可見性分析Fig.5 Annual access analysis of sun in solar calibration

3.2 月球定標

同理，建立月球定標的STK 仿真場景如圖6所示。定義α表示月球和平面YZ的夾角，β表示月球和平面XZ的夾角，θ為月球和Z軸的夾角。同樣要求月球定標時，有α∈[-12°，12°]，β∈[-3.2°，3.2°]，即θ∈[8.69°，90°)，可得在月球可見時間內相應角度的變化，如圖7 所示。

圖6 靜軌紫外可見高光譜探測儀的月球定標仿真場景Fig.6 Lunar calibration simulation of GUVI

圖7 月球可見時間內α，β，θ 角度變化Fig.7 Variations of α，β，θ for moon is visible

以2016～2017 年為例，由圖8 可知在一年每月均有6～8 天對月球可見，因此可以利用月球完成在軌漫反射板監測?？擅吭逻M行2 次同月相下月球定標，進而獲取對月定標數據［16］。

3.3 恒星定標

在STK 場景中建立星等在-10～1 之間的恒星模型，如圖9 所示。同太陽、月球定標分析相同，得到可見恒星如表2 所示。分別分析各恒星的可見性，如圖10～圖12 所示。經分析可得，參宿五、參宿四、南河三每天均有兩次對遙感器可見，可實現遙感器的長期星上定標。

圖8 月球的年可見性分析Fig.8 Annual access analysis of moon in lunar calibration

圖9 恒星仿真場景Fig.9 Stellar calibration simulation in STK

表2 可見恒星參數Tab.2 Parameters of visible stellars

為了研究各恒星的信噪比，本文采用黑體輻射定律估算它在儀器工作波段內，即250～500 nm 內的信噪比。查閱星表可得參宿五、參宿四、南河三的表面溫度分別為22 000，3 500，6 650 K。靜軌衛星遙感器的光譜分辨率為0.6 nm，占據3 個像元，單像元對應0.2 nm。此外，光學系統成像要求像面四個像元能量達80%，單像元達20%。光學系統中反射鏡共4 塊，單塊鏡反射率為0.9。光譜分光采用平面光柵，口徑為35 mm，衍射效率為0.7。探測器的量子效率取0.55。依相應公式分析可得，三顆恒星的信噪比如圖13 所示。由圖可知，南河三在250～500 nm 波段內信噪比優于25，滿足星上定標要求。

通過以上分析，結合探測儀及其軌道特點，可知太陽定標、恒星定標和月球定標的方案可行。因此在三月下旬、四月上旬，以及九月上旬和中旬，即在春分、秋分附近每周進行一次太陽定標。全年每月進行兩次月球定標，每周進行一次恒星定標，定標序列如圖14 所示。

圖10 恒星參宿四的年可見性分析Fig.10 Annual access analysis of Bellartix

圖11 參宿四的年可見性分析Fig.11 Annual access analysis of Betelgeuse

圖12 南河三的年可見性分析Fig.12 Annual access analysis of Procyon

圖13 可見恒星的信噪比分析Fig.13 Signal to noise analysis of visible stellars

4 不確定度分析

依據國際慣例對星上定標各模式進行了不確定度分析，結果如表3～表5 所示。根據國外太陽高精度測量儀器測得的在軌數據，及相關太陽模型得到250～500 nm 波段太陽標準光譜輻照度的不確定度為2%，日地距離變化引入的不確定為0.1%。為降低太陽活動周期的影響，這里采用與儀器定標時期相一致的太陽標準輻照度值，其不確定度可減小至1.50%，光譜定標引入的不確定度為0.50%。則δEsun為2.55%，儀器測量誤差為1.41%。漫反射板不確定度采用實驗室定標數據，即2.00%。太陽和恒星光闌系數實驗室的測量不確定度為0.50%。

表3 靜軌高光譜探測儀太陽定標不確定度Tab. 3 Uncertainty of solar calibration for GUVI（%）

表4 靜軌高光譜探測儀恒星定標不確定度Tab. 4 Uncertainty of stellar calibration for GUVI（%）

圖14 紫外可見高光譜探測儀星上外定標序列Fig.14 External calibration schedule of GUVI

表5 漫反射板監測不確定度Tab. 5 Uncertainty of diffuser （%）

靜軌紫外可見高光譜探測儀采用月球定標完成漫反射板衰減特性監測，每月在相同月相時進行定標，以降低月相變化對觀測結果的影響，約為0.80%。月球天平動影響為1%，月球觀測儀器的掃描精度誤差為0.5%，儀器探測光路的測量誤差為 1.41%［12］。

5 結 論

星上定標是實現遙感器在軌性能監測的重要手段，建立有效可行的星上定標方案是遙感器數據精細化應用的關鍵。本文通過分析太陽、月球、恒星等常用在軌全孔徑定標基準源，確定了靜軌紫外可見高光譜探測儀的星上外定標方案，其在軌不確定度均優于4%。其中，以月球完成漫反射板監測，是目前國內亟需研究的課題，特別是星上恒星定標研究，不僅為在軌輻射基準間比對校正奠定基礎，同時對在軌點擴散函數測量等以恒星作為基準源的進一步應用具有重要意義。