汪洪源，王秉文，毛曉楠,2

(1. 哈爾濱工業(yè)大學(xué) 空間光學(xué)工程研究中心·哈爾濱·150001；2.上海航天控制技術(shù)研究所·上海·201109)

0 引 言

全球定位系統(tǒng)(Global Positioning System，GPS)廣泛應(yīng)用于坦克、飛機(jī)以及艦船的導(dǎo)航，GPS不受戰(zhàn)場(chǎng)煙霧、云層以及沙塵的影響，可以全天候、全天時(shí)工作，定位精度高；但GPS易受電磁干擾，且在戰(zhàn)時(shí)一旦GPS的衛(wèi)星星座遭受敵方攻擊，GPS將會(huì)喪失部分功能，甚至有可能完全癱瘓[1]。慣性導(dǎo)航系統(tǒng)(Inertial Navigation System，INS)彌補(bǔ)了GPS的不足，INS完全自主、隱蔽性好、短時(shí)精度高且數(shù)據(jù)輸出率高，但其長(zhǎng)時(shí)間運(yùn)行會(huì)產(chǎn)生漂移誤差[2-3]。目前，坦克、飛機(jī)以及艦船主要使用GPS和INS組合進(jìn)行導(dǎo)航，GPS用于實(shí)現(xiàn)對(duì)INS漂移誤差的周期性校準(zhǔn)[4]。綜上，亟需尋找一個(gè)可以替代GPS的備用導(dǎo)航定位系統(tǒng)，以便在戰(zhàn)時(shí)當(dāng)GPS不可用時(shí)，該備用導(dǎo)航定位系統(tǒng)可以實(shí)現(xiàn)對(duì)INS漂移誤差的快速校準(zhǔn)，或者該備用導(dǎo)航定位系統(tǒng)可以不依賴(lài)INS而完成自主定姿定位。

作為光測(cè)設(shè)備的星敏感器受背景雜光的影響較大，所以以往星敏感器僅用于衛(wèi)星、宇宙飛船、火箭等空間航天器的導(dǎo)航[5-7]。近年來(lái)，隨著星慣組合導(dǎo)航技術(shù)的不斷發(fā)展，星敏感器導(dǎo)航技術(shù)正由單純的空間應(yīng)用逐漸向更廣泛的空中、近地面以及海平面領(lǐng)域發(fā)展。由于白天天空背景輻射強(qiáng)烈[8-9]，常見(jiàn)的可見(jiàn)光組件在近地空間的探測(cè)能力有限[10]。而相較于可見(jiàn)光相機(jī)，短波紅外相機(jī)在某些特定窗口處(J波段、H波段、Ks波段)的天空背景輻射強(qiáng)度更低、大氣透過(guò)率更高；短波紅外相機(jī)可探測(cè)的恒星數(shù)目更多；短波紅外相機(jī)的滿阱容量更大；短波紅外相機(jī)對(duì)霧、煙以及霾的敏感性更弱；短波紅外相機(jī)受白天海平面湍流的影響更小[11]。面向近地空間短波紅外星敏感器導(dǎo)航應(yīng)用的任務(wù)需求，針對(duì)不同大氣窗口，迫切需要綜合考慮恒星目標(biāo)以及天空背景特性、光學(xué)系統(tǒng)以及圖像傳感器參數(shù)，開(kāi)展恒星探測(cè)信噪比分析方法研究。而在星敏感器恒星探測(cè)信噪比定量計(jì)算方面，現(xiàn)有的研究工作大多集中在可見(jiàn)光波段[12-13]。因此，本文提出了一種近地空間短波紅外恒星探測(cè)信噪比分析方法。首先，分析了短波紅外恒星目標(biāo)與天空背景特性；其次，構(gòu)建了恒星目標(biāo)信號(hào)傳遞模型、天空背景信號(hào)轉(zhuǎn)換模型和恒星探測(cè)信噪比模型；最后，給定了仿真條件，設(shè)計(jì)了仿真流程，并完成了仿真試驗(yàn)。

1 恒星目標(biāo)與天空背景特性

1.1 恒星目標(biāo)特性

(1)波段及數(shù)量

2μm全天巡天(Two Micron All Sky Survey， 2MASS)點(diǎn)源目錄是近地空間短波紅外星敏感器的首選星表，其包括的恒星波段及星等不大于7的恒星數(shù)量如表1所示。

(2)光譜輻照度

2MASS點(diǎn)源目錄中0等星的光譜輻照度如表2所示。

光譜輻照度和星等之間的關(guān)系如下

(1)

表1 2MASS點(diǎn)源目錄波段及數(shù)量

表2 2MASS點(diǎn)源目錄0等星光譜輻照度

式中，m表示星等；E0表示0等星的光譜輻照度；Em表示m星等的光譜輻照度。

1.2 天空背景特性

(1)大氣透過(guò)率

libRadtran是一款大氣輻射傳輸計(jì)算軟件，其可用于準(zhǔn)確分析天空背景特性。選用的求解器為MYSTIC，大氣模式為中緯度夏季，海拔高度為0km，天頂大氣透過(guò)率的仿真結(jié)果如圖1所示。J波段和H波段處的天頂大氣透過(guò)率均為98%，略高于Ks波段處的95%。

圖1 天頂大氣透過(guò)率Fig.1 Zenith atmospheric transmittance

(2)光譜輻亮度

設(shè)置觀測(cè)天頂角為30°，觀測(cè)方位角為180°，太陽(yáng)天頂角分別為0°、20°、40°、60°，太陽(yáng)方位角為0°，天空背景光譜輻亮度的仿真結(jié)果如圖2所示。同一太陽(yáng)天頂角下，J波段下的天空背景輻亮度最大，H波段次之，Ks波段最?。煌徊ǘ蜗?，太陽(yáng)天頂角越小，天空背景輻亮度越大。

圖2 不同太陽(yáng)天頂角下天空背景輻亮度Fig.2 Sky background radiance under different solar zenith angles

2 恒星探測(cè)信噪比模型構(gòu)建

2.1 恒星目標(biāo)信號(hào)傳遞

若已知恒星光譜輻照度、大氣透過(guò)率、光學(xué)系統(tǒng)參數(shù)以及圖像傳感器參數(shù)，則可定量計(jì)算出恒星目標(biāo)信號(hào)經(jīng)光電轉(zhuǎn)換后的等效電子數(shù)。

恒星星光信號(hào)受大氣衰減后，光譜照度水平為

Ea=Emτa

(2)

式中，τa表示大氣透過(guò)率。

再經(jīng)光學(xué)系統(tǒng)成像后，到達(dá)圖像傳感器平面的光譜輻射通量Φs為

(3)

式中，τo表示光學(xué)系統(tǒng)透過(guò)率；D表示光學(xué)系統(tǒng)有效孔徑。

圖像傳感器在積分時(shí)間t內(nèi)接收的恒星目標(biāo)信號(hào)光譜光子數(shù)Ns為

(4)

式中，t表示積分時(shí)間；Eph表示單光子能量；λ表示波長(zhǎng)；h表示普朗克常數(shù)；c表示光速。

考慮圖像傳感器的量子效率，則圖像傳感器生成的恒星星光信號(hào)光譜電子數(shù)Sλ為

(5)

式中，η表示圖像傳感器量子效率。

再考慮J、H、Ks波段的各自帶寬，則圖像傳感器生成的恒星星光信號(hào)電子數(shù)S為

(6)

式中，λ1、λ2分別表示帶寬下限和上限。

2.2 天空背景信號(hào)轉(zhuǎn)換

若已知天空背景輻亮度、光學(xué)系統(tǒng)參數(shù)以及圖像傳感器參數(shù)，則可定量計(jì)算出天空背景信號(hào)經(jīng)光電轉(zhuǎn)換后的等效電子數(shù)。

天空背景信號(hào)經(jīng)光學(xué)系統(tǒng)成像后，到達(dá)圖像傳感器平面的光譜輻射通量Φb為

(7)

式中，Lb表示天空背景信號(hào)光譜輻亮度；α表示圖像傳感器像元尺寸；f表示光學(xué)系統(tǒng)焦距。

圖像傳感器在積分時(shí)間t內(nèi)接收的天空背景信號(hào)光譜光子數(shù)Nb為

(8)

考慮圖像傳感器的量子效率，則圖像傳感器生成的天空背景信號(hào)光譜電子數(shù)Bλ為

(9)

再考慮J、H、Ks波段的各自帶寬，則圖像傳感器生成的天空背景信號(hào)電子數(shù)B為

(10)

2.3 恒星探測(cè)信噪比

根據(jù)恒星目標(biāo)信號(hào)電子數(shù)、天空背景信號(hào)電子數(shù)、圖像傳感器讀出噪聲以及圖像傳感器暗電流，可知近地空間短波紅外星敏感器恒星探測(cè)信噪比rSN為

(11)

3 仿真試驗(yàn)及結(jié)果分析

3.1 仿真輸入條件

仿真輸入條件如表3所示。

表3 仿真輸入條件

3.2 仿真流程設(shè)計(jì)

仿真實(shí)現(xiàn)流程如圖3所示。

圖3 仿真試驗(yàn)流程Fig.3 Simulation test flow

3.3 仿真結(jié)果分析

基于仿真輸入條件及仿真實(shí)現(xiàn)流程，得到不同太陽(yáng)天頂角(0°、20°、40°、60°)、不同波段(J、H、Ks)、不同星等(4、5、6、7)下的恒星探測(cè)信噪比。其中，不同太陽(yáng)天頂角與不同波段下4等星恒星探測(cè)信噪比如圖4所示；不同太陽(yáng)天頂角與不同波段下5等星恒星探測(cè)信噪比如圖5所示；不同太陽(yáng)天頂角與不同波段下6等星恒星探測(cè)信噪比如圖6所示；不同太陽(yáng)天頂角與不同波段下7等星恒星探測(cè)信噪比如圖7所示。

圖4 不同太陽(yáng)天頂角下4等星恒星探測(cè)信噪比Fig.4 SNR for 4 magnitude at different solar zenith angles

圖5 不同太陽(yáng)天頂角下5等星恒星探測(cè)信噪比Fig.5 SNR for 5 magnitude at different solar zenith angles

圖6 不同太陽(yáng)天頂角下6等星恒星探測(cè)信噪比Fig.6 SNR for 6 magnitude at different solar zenith angles

圖7 不同太陽(yáng)天頂角下7等星恒星探測(cè)信噪比Fig.7 SNR for 7 magnitude at different solar zenith angles

同一星等及太陽(yáng)天頂角下，波長(zhǎng)越長(zhǎng)，恒星探測(cè)信噪比越大。以圖4為例，星等為4且太陽(yáng)天頂角為0°時(shí)，J、H、Ks波段的恒星探測(cè)信噪比分別為9.56、13.00、15.07；星等為4且太陽(yáng)天頂角為40°時(shí)，J、H、Ks波段的恒星探測(cè)信噪比分別為11.86、16.12、18.43。

同一太陽(yáng)天頂角及波段下，星等越小，恒星探測(cè)信噪比越大。以圖4～圖7為例，太陽(yáng)天頂角為20°且波段為J時(shí)，4、5、6、7等星的恒星探測(cè)信噪比分別為10.61、4.23、1.69、0.67；太陽(yáng)天頂角為60°且波段為H時(shí)，4、5、6、7等星的恒星探測(cè)信噪比分別為17.05、6.83、2.72、1.09。

同一星等及波長(zhǎng)下，太陽(yáng)天頂角越大，即恒星與太陽(yáng)之間的角距越大，恒星探測(cè)信噪比越大。以圖6為例，星等為6且波段為J時(shí)，0°、20°、40°、60°太陽(yáng)天頂角的恒星探測(cè)信噪比分別為1.52、1.69、1.89、1.99；星等為6且波段為Ks時(shí)，0°、20°、40°、60°太陽(yáng)天頂角的恒星探測(cè)信噪比分別為2.42、2.69、2.98、3.18。

4 結(jié) 論

面向近地空間短波紅外星敏感器導(dǎo)航任務(wù)需求，根據(jù)恒星目標(biāo)與天空背景特性，構(gòu)建了恒星探測(cè)信噪比模型。恒星探測(cè)信噪比的大小不僅與觀測(cè)波長(zhǎng)、太陽(yáng)位置、目標(biāo)星等、背景強(qiáng)度等相關(guān)，而且還受星敏感器光學(xué)系統(tǒng)以及短波紅外圖像傳感器參數(shù)的影響。本文尚未考慮恒星目標(biāo)與天空背景的偏振屬性，因此后續(xù)將研究矢量輻射傳輸效應(yīng)下的恒星探測(cè)問(wèn)題。